Abstract:With discretized particle velocity space, a unified gas-kinetic scheme for entire Knudsen number flows has been constructed based on the Bhatnagar-Gross-Krook (BGK) model (2010. J. Comput. Phys., 229, 7747-7764). In comparison with many existing kinetic schemes for the Boltzmann equation, besides accurate capturing of non-equilibrium flows, the unified method has no difficulty to get accurate solution in the continuum flow regime as well, such as the solution of the Navier-Stokes equations. More importantly, in the continuum flow regime, the time step used by the unified scheme is determined by the Courant-Friedrichs-Lewy condition, which can be many orders larger than the particle collision time. In some sense, the unified method overcomes the time step barrier for many kinetic methods, such as direct simulation Monte Carlo (DSMC), direct Boltzmann solver and many other kinetic solvers. The unified scheme is a multiscale method, where the macroscopic flow variables and microscopic gas distribution function are updated simultaneously. In an early approach in the unified BGK (U-BGK) scheme, the heat flux in the BGK model is modified through the update of macroscopic flow variables, then this modification feeds back into the update of non-equilibrium gas distribution function. In this paper, we are going to develop a unified scheme for the Shakhov model, the so-called unified Shakhov (U-Shk), where the heat flux is corrected directly through the modification of gas distribution function. Theoretically, it will be shown that current U-Shk is more consistent than the U-BGK for the highly non-equilibrium flow computations. The study of shock structures from low to high Mach numbers will be presented and the simulation results will be compared with DSMC solutions as well as possible experimental measurements. The result improvement of U-Shk over U-BGK is clearly achieved. Based on the simulation results, now we fully believe that the unified scheme is an accurate and efficient flow solver in all Knudsen number flow regime.

Regularization of the Chapman-Enskog Expansion and Its Description of Shock Structure

Study of Shock Structures Using the Unified Gas-Kinetic Wave-Particle Method with Various BGK Models

One-Dimensional Non-Equilibrium Shock Wave Structure Analysis by Simplified Conventional Burnett Equations

Computation of 1-D Shock Structure in A Gas in Rotational Non-Equilibrium Using A New Set of Simplified Burnett Equations

Non-equilibrium Shock Structure Calculation with High-order Modified Navier–Stokes Equations

Computation of Shock Wave Structure Using a Simpler Set of Generalized Hydrodynamic Equations Based on Nonlinear Coupled Constitutive Relations

Computation of 1-D Shock Structure Using Nonlinear Coupled Constitutive Relations and Generalized Hydrodynamic Equations

Shock Structure Singularity Smoothed by Adiabatic Approximation of Generalized Hydrodynamics

Continuum Formulation For Non-Equilibrium Shock Structure Calculation

Three-dimensional linearized stability analysis of Burnett equations for a monatomic gas

Modified Gas-Kinetic Scheme for Shock Structures in Argon

Hydrodynamic and Thermodynamic Nonequilibrium Effects around Shock Waves: Based on a Discrete Boltzmann Method

Chapman-Enskog theory for nearly integrable quantum gases

On the stability problem in relativistic thermodynamics: implications of the Chapman-Enskog formalism

Multiple Translational Temperature Model and Its Shock Structure Solution

An Improved Unified Gas-Kinetic Scheme and the Study of Shock Structures

A multi-temperature model of shock structure for aerothermodynamic calculation

The kinetic scheme for the full-Burnett equations

Appendix: Chapman-Enskog Expansion in the Lattice Boltzmann Method

Nonequilibrium Bhatnagar–Gross–Krook model for nitrogen shock structure

The NRxx Method for Polyatomic Gases