An Application of Upwind Difference Scheme with Preconditioned Numerical Fluxes to Gas-Liquid Two-Phase Flows

A time-consistent upwind difference scheme with a preconditioned numerical flux for unsteady gas-liquid multiphase flows is presented and applied to the analysis of cavitating flows. The fundamental equations were formulated in general curvilinear coordinates to apply to diverse flow fields. The preconditioning technique was applied specifically to the numerical dissipation terms in the upwinding process without changing the time derivative terms to maintain time consistency. This approach enhances numerical stability in unsteady multiphase flow computations, consistently delivering time-accurate solutions compared to conventional preconditioning methods. A homogeneous gas-liquid two-phase flow model, third-order Runge-Kutta method, and the flux difference splitting upwind scheme coupled with a third-order MUSCL TVD scheme were employed. Numerical tests of two-dimensional gas-liquid single- and two-phase flows over backward-facing step with different step height and flow conditions successfully demonstrated the capability of the present scheme. The calculations remained stable even for flows with a very low Mach number of 0.001, typically considered incompressible flows, and the results were in good agreement with the experimental data. In addition, we analyzed unsteady cavitating flows at high Reynolds numbers and confirmed the effectiveness and applicability of the present scheme for calculating unsteady gas-liquid two-phase flows.