Implicit Simplified Unified Gas-Kinetic Scheme for Diatomic Rarefied Flows Based on OpenFOAM
Preprint
- 19 June 2026
- preprint
- Published by Elsevier in SSRN Electronic Journal
Abstract
An implicit simplified unified gas-kinetic scheme (ISUGKS) for diatomic gas flows is developed and implemented within the OpenFOAM framework. The solver integrates the Rykov kinetic model to characterize translational–rotational energy exchange and nonequilibrium effects, adopting an implicit macro–micro coupling strategy to eliminate the time-step limitation of explicit schemes. The macroscopic conservation equations are discretized using vector flux splitting, while the microscopic distribution function is solved with a first-order upwind scheme. Both systems are solved by a matrix-free point-relaxation symmetric Gauss–Seidel (PR-SGS) method for low cost and fast convergence. A dual-time-stepping approach is employed to preserve multiscale flow characteristics and accelerate steady-state convergence, and an error compensation technique is introduced to suppress conservation errors caused by velocity-space discretization. The interface flux is reconstructed via the right rectangular integration rule to simplify computation and improve efficiency. Validations are conducted through benchmark cases including a supersonic rarefied flow over a flat plate, a 3D supersonic/hypersonic flow over a sphere, an Apollo re-entry capsule, and a lifting-body vehicle, covering supersonic to hypersonic rarefied flow regimes. Numerical results show good agreement with experimental data, implicit UGKS, and DSMC solutions. Compared with the explicit counterpart, the proposed solver achieves more than one order of magnitude speedup under the same parallel configuration, demonstrating significantly improved efficiency for large-scale engineering simulations. The developed ISUGKS solver fills the gap of implicit multiscale kinetic solvers for diatomic gases in OpenFOAM and provides a reliable tool for hypersonic rarefied aerothermodynamic predictions.Keywords
This publication has 42 references indexed in Scilit:
- A simplified unified wave-particle method for diatomic gases with rotational and vibrational non-equilibriumAerospace Science and Technology, 2025
- Simplified unified wave-particle method for diatomic gases based on Rykov modelActa Mechanica Sinica, 2025
- Performance evaluation of flux-reconstruction schemes in the discrete unified gas-kinetic scheme for low-speed continuum flowsComputers & Mathematics with Applications, 2025
- Implicit unified gas-kinetic scheme for steady state solution of hypersonic thermodynamic non-equilibrium flowsCommunications in Nonlinear Science and Numerical Simulation, 2024
- A multi-prediction implicit scheme for steady state solutions of gas flow in all flow regimesCommunications in Nonlinear Science and Numerical Simulation, 2021
- An open source, parallel DSMC code for rarefied gas flows in arbitrary geometriesComputers & Fluids, 2010
- Direct Methods for Solving the Boltzmann Equation and Study of Nonequilibrium Flows2001
- Discrete-Velocity Models and Numerical Schemes for the Boltzmann-BGK Equation in Plane and Axisymmetric GeometriesJournal of Computational Physics, 2000
- A model kinetic equation for a gas with rotational degrees of freedomFluid Dynamics, 1975
- Generalization of the Krook kinetic relaxation equationFluid Dynamics, 1968