Central Moment Lattice Boltzmann Method using a Pressure-based Formulation for Multiphase Flows at High Density Ratios and including Effects of Surface Tension and Marangoni Stresses

Farzaneh Hajabdollahi,Kannan Premnath,Samuel W.J. Welch

DOI: https://doi.org/10.48550/arXiv.1909.04166

2019-09-10

Abstract:Simulation of multiphase flows require coupled capturing or tracking of the interfaces in conjunction with the solution of fluid motion often occurring at multiple scales. We will present unified cascaded LB methods based on central moments for the solution of the incompressible two-phase flows at high density ratios and for capturing of the interfacial dynamics. Based on a modified continuous Boltzmann equation (MCBE) for two-phase flows, where a kinetic transformation to the distribution function involving the pressure field is introduced to reduce the associated numerical stiffness at high density gradients, a central moment cascaded LB formulation for computing the fluid motion will be constructed. In this LB scheme, the collision step is prescribed by the relaxation of various central moments to their equilibria that are reformulated in terms of the pressure field obtained via matching to the continuous equilibria based on the transformed Maxwell distribution. Furthermore, the differential treatments for the effects of the source term representing the change due to the pressure field and of the source term due to the interfacial tension force and body forces appearing in the MCBE on different moments are consistently accounted for in this cascaded LB solver that computes the pressure and velocity fields. In addition, another cascaded LB scheme via modified equilibria will be developed to solve for the interfacial dynamics represented by a phase field model based on the conservative Allen-Cahn equation. Based on numerical simulations of a variety of two-phase flow benchmark problems at high density ratios and involving the effects of surface tension and its tangential gradients (Marangoni stresses), we will validate our unified cascaded LB approach and also demonstrate improvements in numerical stability.

Computational Physics,Fluid Dynamics

What problem does this paper attempt to address?

This paper aims to solve several key challenges in multiphase fluid simulation at high density ratios, especially in cases involving surface tension and Marangoni stress effects. Specifically, the problems that the paper attempts to solve include:

1. **Numerical stiffness problem**: Under high - density gradients, the traditional Lattice Boltzmann method (LBM) will encounter numerical stiffness problems, leading to computational instability. The paper reduces numerical stiffness by introducing the modified continuous Boltzmann equation (MCBE) and using the pressure field to transform the distribution function.

2. **Accurate capture of interface dynamics**: The interface dynamics in multiphase fluids are very complex and require accurate capture of the interface motion. The paper proposes a phase - field model based on the conservative Allen - Cahn equation (ACE), which describes the evolution of the interface through the competing effects of the diffusion term and the phase - separation flux term.

3. **Combination of multi - relaxation - time (MRT) and central - moment methods**: To improve numerical stability, the paper develops a cascaded LBM method based on central moments, using multiple relaxation times (MRT) to handle the relaxation processes of different central moments. This method can better maintain Galilean invariance and improve numerical stability.

4. **Treatment of source and force terms**: There are multiple source and force terms in multiphase fluids, such as pressure changes, surface tension forces, and body forces. The paper proposes a consistent source/force treatment scheme, which separately considers the effects of these source and force terms on different moments to ensure the accuracy of numerical simulations.

5. **Multiphase fluid simulation at high density ratios**: The paper verifies the effectiveness and stability of the proposed method through a series of benchmark test problems, especially under high - density - ratio and low - viscosity conditions. These test problems include the effects of surface tension and Marangoni stress.

### Formula summary

- **Conservative Allen - Cahn equation (ACE)**:
\[
\frac{\partial\phi}{\partial t}+\nabla\cdot(\phi\mathbf{u})=\nabla\cdot\left[M_{\phi}(\nabla\phi - \theta\mathbf{n})\right]
\]
where $\mathbf{u}$ is the fluid velocity, $\mathbf{n}$ is the unit normal vector, $M_{\phi}$ is the mobility, and $\theta$ is expressed as:
\[
\theta = -\frac{4(\phi - \phi_{A})(\phi - \phi_{B})}{W(\phi_{A}-\phi_{B})}
\]

- **Surface tension force $F_{s}$**:
\[
F_{s}=-\tilde{\kappa}|\nabla\phi|^{2}(\nabla\cdot\mathbf{n})\mathbf{n}
\]
where $\tilde{\kappa}$ is related to the surface tension $\sigma$ as $\tilde{\kappa}=\gamma\sigma W$.

- **Continuous central moment**:
\[
\hat{\Pi}_{M}^{(m,n)}=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}f_{M}(\rho,\mathbf{u})(\xi_{x}-u_{x})^{m}(\xi_{y}-u_{y})^{n}d\xi_{x}d\xi_{y}
\]
\[
\hat{\Pi}_{M}^{(0,m,n)}=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}f_{M}(\rho,0)(\xi_{x}-u_{x})^{m}(\xi_{y}-u_{y})^{n}d\xi_{x}d\xi_{y}
\]

- **Pressure distribution function $g$**:
\[
g = f c_{s}^{2}