S. Ansumali,I.V. Karlin,S. Arcidiacono,A. Abbas,N.I. Prasianakis
Abstract:Exact solution to the hierarchy of nonlinear lattice Boltzmann (LB) kinetic equations in the stationary planar Couette flow is found at non-vanishing Knudsen numbers. A new method of solving LB kinetic equations which combines the method of moments with boundary conditions for populations enables to derive closed-form solutions for all higher-order moments. Convergence of results suggests that the LB hierarchy with larger velocity sets is the novel way to approximate kinetic theory.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to verify and improve the physical validity of the lattice Boltzmann method (LBM) in simulating micro - scale fluid flow, especially its application under non - equilibrium and high Knudsen number conditions. Specifically:
1. **Verify the physical validity of LBM**: Many numerical studies rely on simulations with limited precision, and the question of whether the LBM method can be used as an effective physical model to describe micro - flow has not been clearly answered. Therefore, the author verifies the validity of the LBM method by finding exact solutions.
2. **Improve the LBM model**: The traditional D2Q9 model performs poorly in some cases (such as the transition zone, Kn ≳ 0.1). For this reason, the author introduces the higher - order D2Q16 model to improve the ability to describe non - equilibrium effects, especially the prediction of the boundary Knudsen layer.
3. **Provide exact solutions**: The author provides a closed - form analytical solution for two - dimensional planar Couette flow, including not only the slip velocity, but also higher - order moments such as shear stress and normal stress difference. These solutions can be compared with the results of molecular dynamics (MD) and direct simulation Monte Carlo (DSMC) methods, thereby verifying the accuracy of the LBM method.
4. **Emphasize the importance of boundary conditions**: The author points out that using the standard "bounce - back" boundary condition will result in a zero slip velocity at all Knudsen numbers, indicating that choosing appropriate boundary conditions is crucial for the application of the LBM method in micro - flow.
Through these works, the author shows that the LBM hierarchy has the potential as a new approximation method in describing micro - scale flow, and its accuracy can be improved by increasing the number of discrete velocity sets.
### Formula summary
- **Velocity distribution function of D2Q9 model**:
\[
f_{\text{eq},i}(\rho, j_x, j_y) = w_i \left( \rho + \frac{j_\alpha c_{i\alpha}}{c_s^2} + \frac{j_\alpha j_\beta}{2\rho c_s^4} (c_{i\alpha} c_{i\beta} - c_s^2 \delta_{\alpha\beta}) \right)
\]
where \( c_s=\sqrt{\frac{k_B T_0}{m}} \) is the speed of sound, and \( T_0 \) is the reference temperature.
- **x - direction velocity solution of D2Q9 model**:
\[
u_x=\frac{1}{\Theta_9} \left( \frac{\Delta U}{L} y + U \right)
\]
where \( \Theta_9 = 1 + 2\text{Kn} \), \( \Delta U=U_2 - U_1 \) is the relative velocity between the plates, and \( U=\frac{U_1 + U_2}{2} \) is the centerline velocity.
- **x - direction velocity solution of D2Q16 model**:
\[
u_x=\frac{1}{Z_{16}} \sinh\left( \frac{y}{\text{Kn}L} \right) \Delta U+\frac{1}{\Theta_{16}} \left( \frac{\Delta U}{L} y + U \right)
\]
where
\[
\Theta_{16}=1 + 2\text{Kn} \left( 2\cosh\left( \frac{1}{2\text{Kn}} \right)+\mu \sinh\left( \frac{1}{2\text{Kn}} \right) \right)/\left( \mu \cosh\left( \frac{1}{2\text{Kn}} \right)+2\sqrt{3} \sinh\left( \frac{