A second order description of shock structure

Jason Reese, L.C. Woods, F.J.P. Thivet, S.M. Candel

Research output: Contribution to journalArticle

26 Citations (Scopus)

Abstract

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 Math 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. This article describes a novel method devised to solve the time-independent conservation equations, including 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 derired 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 up-stream 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.
LanguageEnglish
Pages240-250
JournalJournal of Computational Physics
Volume117
Issue number2
DOIs
Publication statusPublished - 15 Mar 1995

Fingerprint

Kinetic theory
kinetic theory
shock
Burnett equations
Mach number
Flow of gases
Wood
Magnetoplasma
Kinetic theory of gases
Hypersonic flow
Boltzmann equation
Gas dynamics
gas flow
Shock waves
Argon
Conservation
Momentum
formulations
hypersonic flow
Boundary conditions

Cite this

Reese, J., Woods, L. C., Thivet, F. J. P., & Candel, S. M. (1995). A second order description of shock structure. Journal of Computational Physics, 117(2), 240-250. https://doi.org/10.1006/jcph.1995.1062
Reese, Jason ; Woods, L.C. ; Thivet, F.J.P. ; Candel, S.M. . / A second order description of shock structure. In: Journal of Computational Physics. 1995 ; Vol. 117, No. 2. pp. 240-250.
@article{6fedef38719046a5ad1969ec176094cd,
title = "A second order description of shock structure",
abstract = "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 Math 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. This article describes a novel method devised to solve the time-independent conservation equations, including 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 derired 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 up-stream 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.",
author = "Jason Reese and L.C. Woods and F.J.P. Thivet and S.M. Candel",
year = "1995",
month = "3",
day = "15",
doi = "10.1006/jcph.1995.1062",
language = "English",
volume = "117",
pages = "240--250",
journal = "Journal of Computational Physics",
issn = "0021-9991",
number = "2",

}

Reese, J, Woods, LC, Thivet, FJP & Candel, SM 1995, 'A second order description of shock structure' Journal of Computational Physics, vol. 117, no. 2, pp. 240-250. https://doi.org/10.1006/jcph.1995.1062

A second order description of shock structure. / Reese, Jason; Woods, L.C.; Thivet, F.J.P.; Candel, S.M. .

In: Journal of Computational Physics, Vol. 117, No. 2, 15.03.1995, p. 240-250.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A second order description of shock structure

AU - Reese, Jason

AU - Woods, L.C.

AU - Thivet, F.J.P.

AU - Candel, S.M.

PY - 1995/3/15

Y1 - 1995/3/15

N2 - 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 Math 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. This article describes a novel method devised to solve the time-independent conservation equations, including 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 derired 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 up-stream 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.

AB - 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 Math 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. This article describes a novel method devised to solve the time-independent conservation equations, including 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 derired 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 up-stream 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.

UR - http://www.sciencedirect.com/science/article/pii/S0021999185710625

U2 - 10.1006/jcph.1995.1062

DO - 10.1006/jcph.1995.1062

M3 - Article

VL - 117

SP - 240

EP - 250

JO - Journal of Computational Physics

T2 - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -