A Robust Momentum Closure Approach For Multiphase Computational Fluid Dynamics Applications