An improved fast iterative moment method for the simulation of near-continuum gas flows with Shakhov model
Abstract
The fast iterative moment (FIM) method developed in [G. Li, C. Wang, and Z. Hu, Comput. Fluids, 311 (2026), 107031] achieves fast convergence for near-continuum gas flows modeled by the Boltzmann-BGK equation through alternating iterations between the high-order moment system and the consistent hydrodynamic equations. In this paper, we extend the FIM solver to the Shakhov model and improve it in the following aspects. A finite volume method with linear reconstruction is employed for high-resolution spatial discretization of both systems. Following the original FIM method, an implicit-explicit scheme treating the Shakhov collision term implicitly is adopted for the high-order moment system. A local Fourier analysis of a linearized one-dimensional BGK model reveals the need for more hydrodynamic iterations at smaller Knudsen numbers and their slower convergence on refined grids. Motivated by these findings, a Newton-multigrid method with a new regularization strategy is introduced for the hydrodynamic equations, generally providing an effective hydrodynamic correction with one step per alternating iteration, thereby reducing the need for empirical tuning. Numerical experiments for planar Couette flow, lid-driven cavity flow, and heat transfer in a bottom-heated cavity are carried out to investigate the performance of the improved FIM solver. The results demonstrate improved efficiency and robustness over the original FIM method, particularly at small Knudsen numbers, and favorable performance with both first- and second-order spatial discretizations.
Related
- Moment obstructions and continuum-to-discrete bounds for checkerboard no-three-in-line sets
- A Novel Solving Method for Multi-Constrained Optimization of Ceramic Composite Armor with Failure Considerations Based on an Improved Response Surface Methodology–Adaptive Moment Estimation (RSM-Adam) Framework
- Grad-div stabilized finite element Picard-Newton iterative methods for the Navier-Stokes equations
- Solving the Yamabe problem by an iterative method on a small Riemannian domain
- Iterative Methods for Solving New Generalized Absolute Value Equations
BibTeX
@misc{cai2026improved,
title = {{An improved fast iterative moment method for the simulation of near-continuum gas flows with Shakhov model}},
author = {Zhenning Cai and Zhicheng Hu and Guanghan Li and Chunwu Wang and Xiaowen Wang},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.33537},
eprint = {2610.33537},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.33537},
}
RIS
TY - UNPB
TI - An improved fast iterative moment method for the simulation of near-continuum gas flows with Shakhov model
AU - Zhenning Cai
AU - Zhicheng Hu
AU - Guanghan Li
AU - Chunwu Wang
AU - Xiaowen Wang
PY - 2026
DA - 2026/10/24
PB - arXiv
DO - 10.48550/arxiv.2610.33537
UR - https://arxiv.org/abs/2610.33537
AB - The fast iterative moment (FIM) method developed in [G. Li, C. Wang, and Z. Hu, Comput. Fluids, 311 (2026), 107031] achieves fast convergence for near-continuum gas flows modeled by the Boltzmann-BGK equation through alternating iterations between the high-order moment system and the consistent hydrodynamic equations. In this paper, we extend the FIM solver to the Shakhov model and improve it in the following aspects. A finite volume method with linear reconstruction is employed for high-resolution spatial discretization of both systems. Following the original FIM method, an implicit-explicit scheme treating the Shakhov collision term implicitly is adopted for the high-order moment system. A local Fourier analysis of a linearized one-dimensional BGK model reveals the need for more hydrodynamic iterations at smaller Knudsen numbers and their slower convergence on refined grids. Motivated by these findings, a Newton-multigrid method with a new regularization strategy is introduced for the hydrodynamic equations, generally providing an effective hydrodynamic correction with one step per alternating iteration, thereby reducing the need for empirical tuning. Numerical experiments for planar Couette flow, lid-driven cavity flow, and heat transfer in a bottom-heated cavity are carried out to investigate the performance of the improved FIM solver. The results demonstrate improved efficiency and robustness over the original FIM method, particularly at small Knudsen numbers, and favorable performance with both first- and second-order spatial discretizations.
ER -