One-dimensional multiple-temperature gas-kinetic Bhatnagar-Gross-Krook scheme for shock wave computation

Accurately computing the inner structure of normal shock waves or oblique shock waves is crucial for many hypersonic applications. As such, it will improve the prediction accuracy of aerodynamics properties and aerothermal effects on hypersonic vehicles and spacecraft during atmospheric entries. Because a shock wave usually has a thickness of a few mean free paths, it is very difficult to accurately compute the detailed nonequilibrium inner structure across a shock wave with a continuum method. This paper reports a gas-kinetic Bhatnagar-Gross-Krook scheme for computations of one-dimensional vibrationally nonequilibrium nitrogen flows through a planar shock wave. The present gas-kinetic Bhatnagar-Gross-Krook scheme solves for the shock structure with multiple temperatures, including two translational temperatures, one rotational temperature, and one vibrational temperature. The salient features of the present gas-kinetic Bhatnagar-Gross-Krook method are multifold. Its applicability covers a wide simulation regime, extending that of continuum flows to the transition flows; it is more computationally efficient in time than the traditional direct simulation Monte Carlo method for shock wave simulation. To provide proper downstream subsonic boundary conditions for very strong shock waves, it is required to determine a proper postshock equilibrium state in which all temperatures have accomplished relaxation processes to a common equilibrium temperature. Analytical expressions of a complete set of generalized Rankine-Hugoniot relations across a planar shock wave are obtained to account for the variant specific heat ratio due to inner energy excitations. Numerical simulation results by the present gas-kinetic Bhatnagar-Gross-Krook scheme and the direct simulation Monte Carlo method are found to be in good agreement.

[1]  Stefan Dietrich,et al.  Scalar and Parallel Optimized Implementation of the Direct Simulation Monte Carlo Method , 1996 .

[2]  V. Aristov.  Direct Methods for Solving the Boltzmann Equation and Study of Nonequilibrium Flows , 2001 .

[3]  Kun Xu,et al.  A multidimensional gas-kinetic BGK scheme for hypersonic viscous flow , 2005 .

[4]  S. P. Gill,et al.  Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena , 2002 .

[5]  R. Lindsay,et al.  Elements of gasdynamics , 1957 .

[6]  L. Talbot,et al.  Experimental Study of the Rotational Distribution Function of Nitrogen in a Shock Wave , 1966 .

[7]  Iain D. Boyd.  Analysis of vibrational‐translational energy transfer using the direct simulation Monte Carlo method , 1991 .

[8]  James N. Moss,et al.  Direct Simulation of Transitional Flow for Hypersonic Reentry Conditions , 1984 .

[9]  Mikhail Naumovich Kogan,et al.  Rarefied Gas Dynamics , 1969 .

[10]  Graham V. Candler,et al.  A multiple translational temperature gas dynamics model , 1994 .

[11]  Kun Xua.  Nonequilibrium Bhatnagar – Gross – Krook model for nitrogen shock structure , 2004 .

[12]  Kun Xu,et al.  Regularization of the Chapman-Enskog expansion and its description of shock structure , 2001 .

[13]  Li-Shi Luo,et al.  Gas-kinetic schemes for direct numerical simulations of compressible homogeneous turbulence. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  Michel A. Saad,et al.  Compressible Fluid Flow , 1985 .

[15]  Iain D. Boyd,et al.  Analysis of rotational nonequilibrium in standing shock waves of nitrogen , 1990 .

[16]  A. J. Eggers.  One-Dimensional Flows of an Imperfect Diatomic Gas , 1949 .

[17]  Kun Xu,et al.  A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method , 2001 .

[18]  G. A. Bird,et al.  Direct Molecular Simulation of a Dissociating Diatomic Gas , 1977 .

[19]  J. G. Parker.  Rotational and Vibrational Relaxation in Diatomic Gases , 1959 .

[20]  Forrest E. Lumpkin,et al.  Resolution of differences between collision number definitions in particle and continuum simulations , 1991 .

[21]  Ian D. Boyd,et al.  Rotational-Translational Energy Transfer in Rare? ed Nonequilibrium Flows , 1990 .

[22]  Kun Xu,et al.  Nonequilibrium Bhatnagar-Gross-Krook model for nitrogen shock structure , 2004 .

[23]  James N. Moss,et al.  Direct Simulation of Transitional Flow for Hypersonic Reentry Conditions , 2003 .

[24]  K. Xu.  Regularization of the Chapman-Enskog Expansion and Its Description of Shock Structure , 2000 .

[25]  H. Alsmeyer,et al.  Density profiles in argon and nitrogen shock waves measured by the absorption of an electron beam , 1976, Journal of Fluid Mechanics.