Stable multiphase moving particle semi-implicit method for incompressible interfacial flow

• Published in: Computer methods in applied mechanics and engineering 2017-05, Vol.318, p.636-666
• Main Authors: Duan, Guangtao, Chen, Bin, Koshizuka, Seiichi, Xiang, Hao
• Format: Article
• Language: English
• Subjects: Acceleration; Computational fluid dynamics; Computer simulation; Continuity (mathematics); Density ratio; Flow stability; Fluid flow; Harmonic analysis; High acceleration; Implicit methods; Incompressible flow; Interface stability; Laminar flow; Multi-density model; Multi-viscosity model; Multiphase flow; Multiphase moving particle semi-implicit (MMPS) method; Poisson distribution; Poisson equation; Pressure gradient model; Stability; Studies; Taylor instability; Viscosity; Viscosity ratio
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2017.01.002

The moving particle semi-implicit (MPS) method is extended into a multiphase MPS (MMPS) method, where multiphase fluids are modeled as a multi-viscosity and multi-density fluid. Interparticle viscosity and density are adopted to model the interaction between particles of different phases. However, such a straightforward extension is prone to instability because the light particles at the interface suffer from exceptionally high acceleration. Therefore, two approaches, MMPS-HD (harmonic density) and MMPS-CA (continuous acceleration), are proposed to suppress the instability. In the first approach, harmonic mean interparticle density is applied to discretize the multiphase pressure Poisson equation to avoid the exceptionally high acceleration at the interface. In the second approach, new MMPS formulations are derived from the locally weighted average of interaction acceleration between particles to guarantee the continuity of acceleration and velocity. The particle stabilizing term (PST) is then decoupled from the original pressure gradient model and adopts the single-phase formulation to guarantee stability. The developed multi-viscosity and -density models are verified using the multi-fluid Poiseuille flow and Rayleigh–Taylor instability, respectively. Furthermore, two benchmark cases of rising bubbles in 2D and 3D with a wide range of density and viscosity ratios are simulated to demonstrate the capability and robustness of the proposed methods in complex multiphase flows. The proposed method can produce stable and reliable results up to a high density ratio of approximately 1000 and viscosity ratio of approximately 100. •Multi-viscosity and multi-density models are developed for Moving Particle Semi-implicit (MPS) method.•Two approaches, MMPS-HD and MMPS-CA, are proposed to stabilize simulations with high density ratio.•Harmonic mean density is utilized to avoid extremely large acceleration at interface in MMPS-HD.•New stable formulations are developed to keep acceleration continuous at interface in MMPS-CA.•Simulations of multiphase flows with high density and viscosity ratios verify the proposed methods.