F.A. Baidoo,I.M. Gamba,T.J.R. Hughes,M.R.A. Abdelmalik
Abstract:We derive a fourth order entropy stable extension of the Navier-Stokes-Fourier equations into the transition regime of rarefied gases. We do this through a novel reformulation of the closure of conservation equations derived from the Boltzmann equation that subsumes existing methods such as the Chapman-Enskog expansion. We apply the linearized version of this extension to the stationary heat problem and the Poiseuille channel and compare our analytical solutions to asymptotic and numerical solutions of the linearized Boltzmann equation. In both model problems, our solutions compare remarkably well in the transition regime. For some macroscopic variables, this agreement even extends far beyond the transition regime.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve the difficult problem of transition - zone modeling in rarefied gas transport problems. Specifically, when the gas is in the transition zone, its mean free path is relatively large compared to the macroscopic reference length scale, making traditional continuum models (such as the Navier - Stokes - Fourier equations) no longer applicable; but at the same time, the mean free path is not large enough to make statistical methods (such as the Direct Simulation Monte Carlo method) computationally efficient. Therefore, how to accurately describe the gas behavior in this transition zone becomes a challenge.
#### Main contributions
1. **Extension of the Navier - Stokes - Fourier equations**: By introducing the Variational Multiscale Moment Methods, the paper proposes a fourth - order entropy - stable extended form of the Navier - Stokes - Fourier equations, which is suitable for the transition zone of rarefied gases.
2. **Improvement of existing models**: In particular, in view of the instability and violation of the second law of thermodynamics in the Burnett equation, a new closure method is proposed to ensure the stability and physical rationality of the new model under different Knudsen numbers.
3. **Application verification**: Through the application to the steady - state heat conduction problem and the Poiseuille channel problem, it shows good consistency between the new model and the analytical solution of the linearized Boltzmann equation, especially excellent performance in the transition zone.
#### Representation of mathematical formulas
- **Boltzmann equation**:
\[
St\frac{\partial F(t, x, v)}{\partial t}+v_i\frac{\partial F(t, x, v)}{\partial x_i}=\frac{1}{\epsilon}C(F)(t, x, v)
\]
where \(St\) is the dimensionless Strouhal number, \(\epsilon\) is the Knudsen number, and \(C(F)\) is the collision operator.
- **Macroscopic conservation equations**:
\[
\begin{aligned}
&\int_{\mathbb{R}^D}F(t, x, v)dv = \rho\\
&\int_{\mathbb{R}^D}v_iF(t, x, v)dv=\rho u_i\\
&\int_{\mathbb{R}^D}v_i v_jF(t, x, v)dv = \rho u_i u_j+P_{ij}\\
&\int_{\mathbb{R}^D}v_i v_j v_kF(t, x, v)dv=\rho u_i u_j u_k+u_i P_{jk}+u_j P_{ik}+u_k P_{ij}+Q_{ijk}
\end{aligned}
\]
- **Stress tensor and heat flux tensor**:
\[
\begin{aligned}
&P_{ij}(t, x)=\int_{\mathbb{R}^D}(v_i - u_i)(v_j - u_j)F(t, x, v)dv\\
&Q_{ijk}(t, x)=\int_{\mathbb{R}^D}(v_i - u_i)(v_j - u_j)(v_k - u_k)F(t, x, v)dv
\end{aligned}
\]
- **Entropy - stable condition**:
\[
\text{For any}\ \bar{m}\in I,\ \langle\bar{m},\mu f\rangle = 0
\]
and the fine - scale equation:
\[
f = L^{-1}_\mu(I-\Pi_\mu)\left[v_i\frac{\partial}{\partial x}