A weakly compressible SPH method for violent multi-phase flows with high density ratio
The weakly compressible SPH (WCSPH) method is known suffering from low computational efficiency, or unnatural voids and unrealistic phase separation when it is applied to simulate highly violent multi-phase flows with high density ratio, such as that between water and air. In this paper, to remedy these issues, we propose a multi-phase WCSPH method based on a low-dissipation Riemann solver and the transport-velocity formulation. The two-phase Riemann problem is first constructed to handle the pairwise interaction between fluid particles, then modified for the fluid-wall interaction to impose the solid wall boundary condition. Since the method uses the same artificial speed of sound for both heavy and light phases, the computational efficiency increases greatly. Furthermore, due to the transport-velocity formulation employed for the light phase and application of the two-phase Riemann problem, the unnatural voids and unrealistic phase separation are effectively eliminated. The method is validated with several 2- and 3D cases involving violent water-air flows. The results have been compared with existing experimental data, previous numerical and analytical solutions, where the proposed method demonstrates good robustness, improved or comparable accuracy, respectively, comparing to previous methods with same choice of sound speed or those with much less computational efficiency.