A discretization of phase mass balance in fractional step algorithms for the drift-flux model