Accuracy analysis of high-order lattice Boltzmann models for rarefied gas flows

In this work, we have theoretically analyzed and numerically evaluated the accuracy of high-order lattice Boltzmann (LB) models for capturing non-equilibrium effects in rarefied gas flows. In the incompressible limit, the LB equation is shown to be able to reduce to the linearized Bhatnagar–Gross–Krook (BGK) equation. Therefore, when the same Gauss–Hermite quadrature is used, LB method closely resembles the discrete velocity method (DVM). In addition, the order of Hermite expansion for the equilibrium distribution function is found not to be directly correlated with the approximation order in terms of the Knudsen number to the BGK equation for incompressible flows. Meanwhile, we have numerically evaluated the LB models for a standing-shear-wave problem, which is designed specifically for assessing model accuracy by excluding the influence of gas molecule/surface interactions at wall boundaries. The numerical simulation results confirm that the high-order terms in the discrete equilibrium distribution function play a negligible role in capturing non-equilibrium effect for low-speed flows. By contrast, appropriate Gauss–Hermite quadrature has the most significant effect on whether LB models can describe the essential flow physics of rarefied gas accurately. Our simulation results, where the effect of wall/gas interactions is excluded, can lead to conclusion on the LB modeling capability that the models with higher-order quadratures provide more accurate results. For the same order Gauss–Hermite quadrature, the exact abscissae will also modestly influence numerical accuracy. Using the same Gauss–Hermite quadrature, the numerical results of both LB and DVM methods are in excellent agreement for flows across a broad range of the Knudsen numbers, which confirms that the LB simulation is similar to the DVM process. Therefore, LB method can offer flexible models suitable for simulating continuum flows at the Navier–Stokes level and rarefied gas flows at the linearized Boltzmann model equation level

Lattice ellipsoidal statistical BGK model for thermal non-equilibrium flows

A thermal lattice Boltzmann model is constructed on the basis of the ellipsoidal statistical Bhatnagar-Gross-Krook (ES-BGK) collision operator via the Hermite moment representation. The resulting lattice ES-BGK model uses a single distribution function and features an adjustable Prandtl number. Numerical simulations show that using a moderate discrete velocity set, this model can accurately recover steady and transient solutions of the ES-BGK equation in the slip-flow and early transition regimes in the small Mach number limit that is typical of microscale problems of practical interest. In the transition regime in particular, comparisons with numerical solutions of the ES-BGK model, direct Monte Carlo and low-variance deviational Monte Carlo simulations show good accuracy for values of the Knudsen number up to approximately 0:5. On the other hand, highly non-equilibrium phenomena characterized by high Mach numbers, such as viscous heating and force-driven Poiseuille flow for large values of the driving force, are more difficult to capture quantitatively in the transition regime using discretizations that have been chosen with computational efficiency in mind such as the one used here, although improved accuracy is observed as the number of discrete velocities is increased

Analytical solution for the lattice Boltzmann model beyond Naviers-Stokes

To understand lattice Boltzmann model capability for capturing nonequilibrium effects, the model with first-order expansion of the equilibrium distribution function is analytically investigated. In particular, the velocity profile of Couette flows is exactly obtained for the D2Q9 model, which shows retaining the first order expansion can capture rarefaction effects in the incompressible limit. Meanwhile, it clearly demonstrates that the D2Q9 model is not able to reflect flow characteristics in the Knudsen layer