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

Guangtao Duan, Bin Chen, Seiichi Koshizuka, Hao Xiang

Research output: Contribution to journalArticle

25 Citations (Scopus)

Abstract

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.

Original languageEnglish
Pages (from-to)636-666
Number of pages31
JournalComputer Methods in Applied Mechanics and Engineering
Volume318
DOIs
Publication statusPublished - 2017 May 1
Externally publishedYes

Fingerprint

Viscosity
viscosity
high acceleration
fluids
Fluids
Particle interactions
Multiphase flow
Poisson equation
Pressure gradient
harmonics
formulations
Flow of fluids
multiphase flow
laminar flow
pressure gradients
continuity
bubbles
interactions

Keywords

  • Multi-density model
  • Multi-viscosity model
  • Multiphase flow
  • Multiphase moving particle semi-implicit (MMPS) method
  • Pressure gradient model
  • Stability

ASJC Scopus subject areas

  • Computational Mechanics
  • Mechanics of Materials
  • Mechanical Engineering
  • Physics and Astronomy(all)
  • Computer Science Applications

Cite this

Stable multiphase moving particle semi-implicit method for incompressible interfacial flow. / Duan, Guangtao; Chen, Bin; Koshizuka, Seiichi; Xiang, Hao.

In: Computer Methods in Applied Mechanics and Engineering, Vol. 318, 01.05.2017, p. 636-666.

Research output: Contribution to journalArticle

Duan, Guangtao ; Chen, Bin ; Koshizuka, Seiichi ; Xiang, Hao. / Stable multiphase moving particle semi-implicit method for incompressible interfacial flow. In: Computer Methods in Applied Mechanics and Engineering. 2017 ; Vol. 318. pp. 636-666.
@article{a70ddf3a056741828c99940f05e073f6,
title = "Stable multiphase moving particle semi-implicit method for incompressible interfacial flow",
abstract = "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.",
keywords = "Multi-density model, Multi-viscosity model, Multiphase flow, Multiphase moving particle semi-implicit (MMPS) method, Pressure gradient model, Stability",
author = "Guangtao Duan and Bin Chen and Seiichi Koshizuka and Hao Xiang",
year = "2017",
month = "5",
day = "1",
doi = "10.1016/j.cma.2017.01.002",
language = "English",
volume = "318",
pages = "636--666",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0374-2830",
publisher = "Elsevier",

}

TY - JOUR

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

AU - Duan, Guangtao

AU - Chen, Bin

AU - Koshizuka, Seiichi

AU - Xiang, Hao

PY - 2017/5/1

Y1 - 2017/5/1

N2 - 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.

AB - 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.

KW - Multi-density model

KW - Multi-viscosity model

KW - Multiphase flow

KW - Multiphase moving particle semi-implicit (MMPS) method

KW - Pressure gradient model

KW - Stability

UR - http://www.scopus.com/inward/record.url?scp=85013833936&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85013833936&partnerID=8YFLogxK

U2 - 10.1016/j.cma.2017.01.002

DO - 10.1016/j.cma.2017.01.002

M3 - Article

AN - SCOPUS:85013833936

VL - 318

SP - 636

EP - 666

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0374-2830

ER -