Revised Chapman-Enskog analysis for a class of forcing schemes in the lattice Boltzmann method

Kinetic Simulations of Compressible Non-Ideal Fluids: From Supercritical Flows to Phase-Change and Exotic Behavior

We investigate a kinetic model for compressible non-ideal fluids. The model imposes the local thermodynamic pressure through a rescaling of the particle’s velocities, which accounts for both long- and short-range effects and hence full thermodynamic consistency. The model is fully Galilean invariant and treats mass, momentum, and energy as local conservation laws. The analysis and derivation of the hydrodynamic limit is followed by the assessment of accuracy and robustness through benchmark simulations ranging from the Joule–Thompson effect to a phase-change and high-speed flows. In particular, we show the direct simulation of the inversion line of a van der Waals gas followed by simulations of phase-change such as the one-dimensional evaporation of a saturated liquid, nucleate, and film boiling and eventually, we investigate the stability of a perturbed strong shock front in two different fluid mediums. In all of the cases, we find excellent agreement with the corresponding theoretical analysis and experimental correlations. We show that our model can operate in the entire phase diagram, including super- as well as sub-critical regimes and inherently captures phase-change phenomena.

Density gradient calculation in a class of multiphase lattice Boltzmann models

The multiphase lattice Boltzmann (LB) models based on pairwise interactions show great potential for simulating multiphase flows due to the conceptual and computational simplicity. Although the dynamics of multiphase flows are reproduced by the pairwise interaction force, the gradient of density (or effective density, i.e., pseudopotential) is implicitly involved in these models via the specialized forcing scheme or the consistent scheme for ɛ^{3}-order term. This work focuses on the calculation of density gradient in this class of multiphase LB models. Theoretical analyses are first carried out to reveal the involvement and calculation of density gradient. On the basis of a low Mach number approximation, an improved scheme is then proposed to calculate the density gradient for the recent LB model with self-tuning equation of state. Analytical and numerical calculations show that the improved scheme is more accurate and can help to reduce the numerical error when the reduced temperature is relatively low.

Consistent forcing scheme in the cascaded lattice Boltzmann method

In this paper, we give an alternative derivation for the cascaded lattice Boltzmann method (CLBM) within a general multiple-relaxation-time (MRT) framework by introducing a shift matrix. When the shift matrix is a unit matrix, the CLBM degrades into an MRT LBM. Based on this, a consistent forcing scheme is developed for the CLBM. The consistency of the nonslip rule, the second-order convergence rate in space, and the property of isotropy for the consistent forcing scheme is demonstrated through numerical simulations of several canonical problems. Several existing forcing schemes previously used in the CLBM are also examined. The study clarifies the relation between MRT LBM and CLBM under a general framework.