Mass-conserved volumetric lattice Boltzmann method for complex flows with willfully moving boundaries

Huidan Yu, Xi Chen, Zhiqiang Wang, Debanjan Deep, Everton Lima, Ye Zhao, Shawn D. Teague

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

In this paper, we develop a mass-conserved volumetric lattice Boltzmann method (MCVLBM) for numerically solving fluid dynamics with willfully moving arbitrary boundaries. In MCVLBM, fluid particles are uniformly distributed in lattice cells and the lattice Boltzmann equations deal with the time evolution of the particle distribution function. By introducing a volumetric parameter P(x,y,z,t) defined as the occupation of solid volume in the cell, we distinguish three types of lattice cells in the simulation domain: solid cell (pure solid occupation, P=1), fluid cell (pure fluid occupation, P=0), and boundary cell (partial solid and partial fluid, 0<P<1). The formulation of volumetric lattice Boltzmann equations are self-regularized through P and consist of three parts: (1) collision taking into account the momentum exchange between the willfully moving boundary and the flow; (2) streaming accompanying a volumetric bounce-back procedure in boundary cells; and (3) boundary-induced volumetric fluid migration moving the residual fluid particles into the flow domain when the boundary swipes over a boundary cell toward a solid cell. The MCVLBM strictly satisfies mass conservation and can handle irregular boundary orientation and motion with respect to the mesh. Validation studies are carried out in four cases. The first is to simulate fluid dynamics in syringes focusing on how MCVLBM captures the underlying physics of flow driven by a willfully moving piston. The second and third cases are two-dimensional (2D) peristaltic flow and three-dimensional (3D) pipe flow, respectively. In each case, we compare the MCVLBM simulation result with the analytical solution and achieve quantitatively good agreements. The fourth case is to simulate blood flow in human aortic arteries with a very complicated irregular boundary. We study steady flow in two dimensions and unsteady flow via the pulsation of the cardiac cycle in three dimensions. In the 2D case, both vector (velocity) and scalar (pressure) fields are compared to computation results from a well-established Navier-Stokes solver and reasonably good agreements are obtained. In the 3D case, the unsteady flow pattern and wall shear stress are well captured at the representative time instants during the pulsation. The validations demonstrate that the MCVLBM is a relatively simple but reliable computation scheme to deal with static or moving irregular boundaries.

Original languageEnglish
Article number063304
JournalPhysical Review E - Statistical, Nonlinear, and Soft Matter Physics
Volume89
Issue number6
DOIs
StatePublished - Jun 11 2014

Fingerprint

Moving Boundary
Lattice Boltzmann Method
Cell
cells
Fluid
fluids
Lattice Boltzmann Equation
Irregular
occupation
Unsteady Flow
Fluid Dynamics
unsteady flow
fluid dynamics
Peristaltic Flow
Partial
Pipe Flow
Wall Shear Stress
Bounce
syringes
Mass Conservation

ASJC Scopus subject areas

  • Condensed Matter Physics
  • Statistical and Nonlinear Physics
  • Statistics and Probability

Cite this

Mass-conserved volumetric lattice Boltzmann method for complex flows with willfully moving boundaries. / Yu, Huidan; Chen, Xi; Wang, Zhiqiang; Deep, Debanjan; Lima, Everton; Zhao, Ye; Teague, Shawn D.

In: Physical Review E - Statistical, Nonlinear, and Soft Matter Physics, Vol. 89, No. 6, 063304, 11.06.2014.

Research output: Contribution to journalArticle

Yu, Huidan ; Chen, Xi ; Wang, Zhiqiang ; Deep, Debanjan ; Lima, Everton ; Zhao, Ye ; Teague, Shawn D. / Mass-conserved volumetric lattice Boltzmann method for complex flows with willfully moving boundaries. In: Physical Review E - Statistical, Nonlinear, and Soft Matter Physics. 2014 ; Vol. 89, No. 6.
@article{92f74aa9209c4c7c84135c9eae554cc2,
title = "Mass-conserved volumetric lattice Boltzmann method for complex flows with willfully moving boundaries",
abstract = "In this paper, we develop a mass-conserved volumetric lattice Boltzmann method (MCVLBM) for numerically solving fluid dynamics with willfully moving arbitrary boundaries. In MCVLBM, fluid particles are uniformly distributed in lattice cells and the lattice Boltzmann equations deal with the time evolution of the particle distribution function. By introducing a volumetric parameter P(x,y,z,t) defined as the occupation of solid volume in the cell, we distinguish three types of lattice cells in the simulation domain: solid cell (pure solid occupation, P=1), fluid cell (pure fluid occupation, P=0), and boundary cell (partial solid and partial fluid, 0<P<1). The formulation of volumetric lattice Boltzmann equations are self-regularized through P and consist of three parts: (1) collision taking into account the momentum exchange between the willfully moving boundary and the flow; (2) streaming accompanying a volumetric bounce-back procedure in boundary cells; and (3) boundary-induced volumetric fluid migration moving the residual fluid particles into the flow domain when the boundary swipes over a boundary cell toward a solid cell. The MCVLBM strictly satisfies mass conservation and can handle irregular boundary orientation and motion with respect to the mesh. Validation studies are carried out in four cases. The first is to simulate fluid dynamics in syringes focusing on how MCVLBM captures the underlying physics of flow driven by a willfully moving piston. The second and third cases are two-dimensional (2D) peristaltic flow and three-dimensional (3D) pipe flow, respectively. In each case, we compare the MCVLBM simulation result with the analytical solution and achieve quantitatively good agreements. The fourth case is to simulate blood flow in human aortic arteries with a very complicated irregular boundary. We study steady flow in two dimensions and unsteady flow via the pulsation of the cardiac cycle in three dimensions. In the 2D case, both vector (velocity) and scalar (pressure) fields are compared to computation results from a well-established Navier-Stokes solver and reasonably good agreements are obtained. In the 3D case, the unsteady flow pattern and wall shear stress are well captured at the representative time instants during the pulsation. The validations demonstrate that the MCVLBM is a relatively simple but reliable computation scheme to deal with static or moving irregular boundaries.",
author = "Huidan Yu and Xi Chen and Zhiqiang Wang and Debanjan Deep and Everton Lima and Ye Zhao and Teague, {Shawn D.}",
year = "2014",
month = "6",
day = "11",
doi = "10.1103/PhysRevE.89.063304",
language = "English",
volume = "89",
journal = "Physical Review E",
issn = "2470-0045",
publisher = "American Physical Society",
number = "6",

}

TY - JOUR

T1 - Mass-conserved volumetric lattice Boltzmann method for complex flows with willfully moving boundaries

AU - Yu, Huidan

AU - Chen, Xi

AU - Wang, Zhiqiang

AU - Deep, Debanjan

AU - Lima, Everton

AU - Zhao, Ye

AU - Teague, Shawn D.

PY - 2014/6/11

Y1 - 2014/6/11

N2 - In this paper, we develop a mass-conserved volumetric lattice Boltzmann method (MCVLBM) for numerically solving fluid dynamics with willfully moving arbitrary boundaries. In MCVLBM, fluid particles are uniformly distributed in lattice cells and the lattice Boltzmann equations deal with the time evolution of the particle distribution function. By introducing a volumetric parameter P(x,y,z,t) defined as the occupation of solid volume in the cell, we distinguish three types of lattice cells in the simulation domain: solid cell (pure solid occupation, P=1), fluid cell (pure fluid occupation, P=0), and boundary cell (partial solid and partial fluid, 0<P<1). The formulation of volumetric lattice Boltzmann equations are self-regularized through P and consist of three parts: (1) collision taking into account the momentum exchange between the willfully moving boundary and the flow; (2) streaming accompanying a volumetric bounce-back procedure in boundary cells; and (3) boundary-induced volumetric fluid migration moving the residual fluid particles into the flow domain when the boundary swipes over a boundary cell toward a solid cell. The MCVLBM strictly satisfies mass conservation and can handle irregular boundary orientation and motion with respect to the mesh. Validation studies are carried out in four cases. The first is to simulate fluid dynamics in syringes focusing on how MCVLBM captures the underlying physics of flow driven by a willfully moving piston. The second and third cases are two-dimensional (2D) peristaltic flow and three-dimensional (3D) pipe flow, respectively. In each case, we compare the MCVLBM simulation result with the analytical solution and achieve quantitatively good agreements. The fourth case is to simulate blood flow in human aortic arteries with a very complicated irregular boundary. We study steady flow in two dimensions and unsteady flow via the pulsation of the cardiac cycle in three dimensions. In the 2D case, both vector (velocity) and scalar (pressure) fields are compared to computation results from a well-established Navier-Stokes solver and reasonably good agreements are obtained. In the 3D case, the unsteady flow pattern and wall shear stress are well captured at the representative time instants during the pulsation. The validations demonstrate that the MCVLBM is a relatively simple but reliable computation scheme to deal with static or moving irregular boundaries.

AB - In this paper, we develop a mass-conserved volumetric lattice Boltzmann method (MCVLBM) for numerically solving fluid dynamics with willfully moving arbitrary boundaries. In MCVLBM, fluid particles are uniformly distributed in lattice cells and the lattice Boltzmann equations deal with the time evolution of the particle distribution function. By introducing a volumetric parameter P(x,y,z,t) defined as the occupation of solid volume in the cell, we distinguish three types of lattice cells in the simulation domain: solid cell (pure solid occupation, P=1), fluid cell (pure fluid occupation, P=0), and boundary cell (partial solid and partial fluid, 0<P<1). The formulation of volumetric lattice Boltzmann equations are self-regularized through P and consist of three parts: (1) collision taking into account the momentum exchange between the willfully moving boundary and the flow; (2) streaming accompanying a volumetric bounce-back procedure in boundary cells; and (3) boundary-induced volumetric fluid migration moving the residual fluid particles into the flow domain when the boundary swipes over a boundary cell toward a solid cell. The MCVLBM strictly satisfies mass conservation and can handle irregular boundary orientation and motion with respect to the mesh. Validation studies are carried out in four cases. The first is to simulate fluid dynamics in syringes focusing on how MCVLBM captures the underlying physics of flow driven by a willfully moving piston. The second and third cases are two-dimensional (2D) peristaltic flow and three-dimensional (3D) pipe flow, respectively. In each case, we compare the MCVLBM simulation result with the analytical solution and achieve quantitatively good agreements. The fourth case is to simulate blood flow in human aortic arteries with a very complicated irregular boundary. We study steady flow in two dimensions and unsteady flow via the pulsation of the cardiac cycle in three dimensions. In the 2D case, both vector (velocity) and scalar (pressure) fields are compared to computation results from a well-established Navier-Stokes solver and reasonably good agreements are obtained. In the 3D case, the unsteady flow pattern and wall shear stress are well captured at the representative time instants during the pulsation. The validations demonstrate that the MCVLBM is a relatively simple but reliable computation scheme to deal with static or moving irregular boundaries.

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

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

U2 - 10.1103/PhysRevE.89.063304

DO - 10.1103/PhysRevE.89.063304

M3 - Article

C2 - 25019909

AN - SCOPUS:84902449523

VL - 89

JO - Physical Review E

JF - Physical Review E

SN - 2470-0045

IS - 6

M1 - 063304

ER -