A second-order description of shock structure

Journal Article · Wed Mar 15 04:00:00 EST 1995 · Journal of Computational Physics

DOI:https://doi.org/10.1006/jcph.1995.1062· OSTI ID:105440

Reese, J M; Woods, L C ^[1]; Thivet, F J.P. ^[2]; Candel, S M ^[3]

---

1. Univ. Oxford (United Kingdom)
2. Centre d`Etudes et de Recherches de Toulouse (France)
3. Laboratoire d`Energetique Moleculaire et Marcroscopique, Chatenay-Malabry (France)

The structure of gas-dynamic shock waves is of interest in hypersonic flow studies and also constitutes a straightforward test for competing kinetic theories. The description of the shock profiles may be obtained from a second-order theory in the Knudsen number. The BGK approximation to the Boltzmann equation introduces additional terms in the transport of momentum and energy. These relations, known as the Burnett equations, improve the agreement between calculated shock profiles and experiment. However, for some formulations of these equations, the solution breaks down at a critical Mach number. In addition, certain terms in the Burnett equations allow unphysical effects in gas flow. A modified kinetic theory has been proposed by Woods (An Introduction to the Kinetic Theory of Gases and Magnetoplasmas, Oxford Univ. Press, Oxford, 1993) which eliminates the frame dependence of the standard kinetic theory and corrects some of the second-order terms. The method is used to solve the shock structure problem in one dimension. It is based on a finite difference global scheme (FDGS), in which a Newton procedure is applied to a discretized version of the governing equations and boundary conditions. The method is first applied to the Navier-Stokes formulation of the shock equations. It is then successfully used to integrate a modified version of the second-order equations derived by Woods for monatomic gases, up to a Mach number of 30. Results of the calculations are compared with experimental data for Argon gas flows characterized by upstream Mach numbers up to 10. The agreement is good, well within the data point spread. The FDGS method converges rapidly and it may be used to study other problems of the same general nature. 13 refs., 6 figs.

• OSTI ID: 105440
• Journal Information: Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: 2 Vol. 117; ISSN JCTPAH; ISSN 0021-9991
• Country of Publication: United States
• Language: English

Related Subjects

42 ENGINEERING
66 PHYSICS
99 GENERAL AND MISCELLANEOUS
FINITE DIFFERENCE METHOD
GAS FLOW
HYPERSONIC FLOW
NAVIER-STOKES EQUATIONS
NEWTON METHOD
SHOCK WAVES