A stable SPH model with large CFL numbers for multi-phase flows with large density ratios

The discontinuity across interface of multi-phase flows with large density ratios usually poses great challenges for numerical simulations. The smoothed particle hydrodynamics (SPH) is a meshless method with inherent advantages in dealing with multi-phase flows without the necessity of tracking the moving interfaces. In this paper, we develop a new weakly-compressible SPH model for multi-phase flows with large density ratios while allowing large CFL numbers. In the present SPH model, the continuity equation is first modified by eliminating the influence from particles of different phases based on the simple fact that different phases will not contribute when calculating the density for immiscible multi-phase flows; thus, the modified continuity equation will only consider the influence from neighboring particles of the same phase. The pressure and density of the particles of other phases are then re-initialized by using the Shepard interpolation function. The present multi-phase SPH model has been tested by four numerical examples, including the two-phase hydrostatic water, standing waves, liquid sloshing, and dam breaking. It has been demonstrated that the present multi-phase SPH model can obtain satisfactory results stably, even at large CFL numbers, and this means that large time steps can be employed. Therefore, the present multi-phase SPH model can significantly save computational cost through using large time steps, especially for large-scale problems with a large number of particles.