Fokker-Planck Central Moment Lattice Boltzmann Method for Effective Simulations of Fluid Dynamics

William Schupbach,Kannan Premnath

2024-10-13

Abstract:We present a new formulation of the central moment lattice Boltzmann (LB) method based on a continuous Fokker-Planck (FP) kinetic model, originally proposed for stochastic diffusive-drift processes (e.g., Brownian dynamics), by adapting it as a collision model for the continuous Boltzmann equation (CBE) for fluid dynamics. The FP collision model has several desirable properties, including its ability to preserve the quadratic nonlinearity of the CBE, unlike that based on the common Bhatnagar-Gross-Krook model. Rather than using an equivalent Langevin equation as a proxy, we construct our approach by directly matching the changes in different discrete central moments independently supported by the lattice under collision to those given by the CBE under the FP-guided collision model. This can be interpreted as a new path for the collision process in terms of the relaxation of the various central moments to 'equilibria', which we term as the Markovian central moment attractors that depend on a diffusion coefficient tensor. The construction of the method using central moments rather than via distribution functions facilitates its numerical implementation and analysis. We show its consistency to the Navier-Stokes equations via a Chapman-Enskog analysis and elucidate the choice of the diffusion coefficient based on the second order moments in accurately representing flows at relatively low viscosities. We will demonstrate the accuracy and robustness of our new central moment FP-LB formulation, termed as the FPC-LBM, using the D3Q27 lattice for simulations of a variety of flows, including wall-bounded turbulent flows. We show that the FPC-LBM is more stable than other existing LB schemes based on central moments, while avoiding numerical hyperviscosity effects in flow simulations at relatively very low physical fluid viscosities through a refinement to a model founded on kinetic theory.

Fluid Dynamics,Computational Engineering, Finance, and Science,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the numerical instability and numerical hyper - viscosity problems existing in the existing Lattice Boltzmann Method (LBM) when simulating fluid flows with low viscosity or high Reynolds numbers. Specifically, traditional LBM collision models such as the Single Relaxation Time (SRT) model and the Multiple Relaxation Time (MRT) model are prone to generate numerical errors and instabilities when dealing with fluid flows under these extreme conditions. These problems limit the effectiveness and accuracy of LBM in practical applications.

To overcome these problems, the author proposes a Central - Moment Lattice Boltzmann Method based on the Fokker - Planck (FP) collision model (FPC - LBM). This method constructs a new collision process by directly matching the changes of discrete central moments with the changes of the continuous Boltzmann equation under the FP collision model. The FP collision model has the advantage of retaining the quadratic non - linearity of the Boltzmann collision integral, which makes FPC - LBM more stable and accurate when dealing with fluid flows with low viscosity or high Reynolds numbers.

### Main contributions of the paper:
1. **Propose a new FPC - LBM**: Based on the FP collision model, a new LBM method is constructed through the matching principle of central moments.
2. **Improve numerical stability**: By adjusting the diffusion coefficient tensor, the numerical hyper - viscosity problem is avoided, and the stability of the method in simulating fluid flows with low viscosity or high Reynolds numbers is improved.
3. **Theoretical analysis**: Through Chapman - Enskog analysis, the consistency of FPC - LBM is proved, and the change rate of central moments is derived in detail.
4. **Numerical verification**: Through a series of standard benchmark tests, the accuracy and stability of FPC - LBM are verified, and it is compared with other existing LBM methods.

### Mathematical formulas:
- **FP collision model**:
\[
\left(\frac{\delta f}{\delta t}\right)_{\text{FP coll}}=\omega_{\text{FP}}\left[\frac{\partial}{\partial \xi_i}\left((\xi_i - u_i) f\right)+D'_{ij}\frac{\partial^2 f}{\partial \xi_i\partial \xi_j}\right]
\]
where \(\omega_{\text{FP}}\) is the collision frequency and \(D'_{ij}\) is the diffusion tensor parameter.

- **Definition of central moments**:
\[
\Pi_{mn}=\langle f, W_{mn}\rangle
\]
where \(W_{mn}=(\xi_x - u_x)^m(\xi_y - u_y)^n\).

- **Change rate of central moments**:
\[
\left(\frac{\delta \Pi_{mn}}{\delta t}\right)_{\text{FP coll}}=\omega_{\text{FP}}\left[\left(\frac{\delta \Pi_{mn}^{\text{FP1}}}{\delta t}\right)+\left(\frac{\delta \Pi_{mn}^{\text{FP1}}}{\delta t}\right)\right]
\]

Through these improvements, FPC - LBM can provide more stable and accurate simulation results under a wider range of fluid flow conditions.