A Physical-Constraint-Preserving Discontinuous Galerkin Method for Weakly Compressible Two-Phase Flows

This work focuses on the robust and high-order numerical simulations of weakly compressible two-phase flows by using the discontinuous Galerkin (DG) method combined with an explicit strong-stability-preserving Runge–Kutta scheme. In order to improve the computational robustness under large density ratios, a nonlinear weighted essentially non-oscillatory (WENO) limiter and a positivity-preserving limiter are specially designed and applied with the aim of dampening the nonphysical oscillations around the phase interface and preventing the occurrence of negative density, respectively. More importantly, we theoretically prove that the present method is able to satisfy the uniform-pressure–velocity criterion which states that uniform pressure and velocity profiles around an isolated phase interface should be preserved during the simulation. The performance of the present method is validated by a range of benchmark test cases with density ratios up to 1000:1. The results demonstrate that the present method possesses a good capability of simulating weakly compressible two-phase flows with large density ratios.