Modeling flows with curved detonation waves

Ralph Menikoff, Klaus Lackner, Bruce G. Bukiet

Research output: Contribution to journalArticle

21 Citations (Scopus)

Abstract

The nonzero width of a curved detonation front has a significant effect on its propagation. Physically, the curvature effect is a small but important correction. In numerical simulations, the curvature effect tends to be greatly exaggerated due to an artificially large width of an underresolved wave. In the context of reactive fluid flow, a detonation wave consists of a lead shock followed by a thin reaction zone. The curvature effect is determined by the dynamics within the reaction zone; in particular, the competition between a source term for the rate of chemical energy release and a geometric source term due to front curvature. When the width of the reaction zone is small compared with the radius of curvature of the front, the reaction zone can be approximated as quasi-steady and be modeled locally by a system of ordinary differential equations (ODEs) in which the front curvature enters as a parameter. Front curvature breaks the Galilean invariance and the quasi-steady approximation is only valid in a distinguished frame determined by both the front curvature and the tangential velocity divergence ahead of the front. The quasi-steady ODEs determine the reaction zone profile to leading order and can be viewed as an extension of the ZND model. However, the source terms in the ODEs lead to modified Hugoniot jump conditions that take into account front curvature and reaction zone width. When the reaction zone is underresolved, a calculation in effect reduces to a 'capturing algorithm' in which the burn model plays an analogous role for detonation waves as artificial viscosity does for shock waves. In particular, the reaction zone has an unphysically large width that is proportional to the cell size. As a consequence of the modified jump conditions, the numerical reaction zone width gives rise to an artificial curvature effect. This causes numerical solutions to depend on the cell size and orientation of a detonation front relative to the grid. Two algorithms that eliminate numerical curvature effects are discussed, detonation shock dynamics and front tracking.

Original languageEnglish (US)
Pages (from-to)219-240
Number of pages22
JournalCombustion and Flame
Volume104
Issue number3
DOIs
StatePublished - Feb 1996
Externally publishedYes

Fingerprint

detonation waves
Detonation
curvature
Ordinary differential equations
detonation
differential equations
Invariance
Shock waves
Flow of fluids
Lead
Viscosity
shock
chemical energy
Computer simulation
cells
fluid flow
shock waves
invariance
divergence
grids

ASJC Scopus subject areas

  • Energy Engineering and Power Technology
  • Fuel Technology
  • Mechanical Engineering

Cite this

Modeling flows with curved detonation waves. / Menikoff, Ralph; Lackner, Klaus; Bukiet, Bruce G.

In: Combustion and Flame, Vol. 104, No. 3, 02.1996, p. 219-240.

Research output: Contribution to journalArticle

Menikoff, Ralph ; Lackner, Klaus ; Bukiet, Bruce G. / Modeling flows with curved detonation waves. In: Combustion and Flame. 1996 ; Vol. 104, No. 3. pp. 219-240.
@article{5caa31ff81d84929b57cf96b81cfca7d,
title = "Modeling flows with curved detonation waves",
abstract = "The nonzero width of a curved detonation front has a significant effect on its propagation. Physically, the curvature effect is a small but important correction. In numerical simulations, the curvature effect tends to be greatly exaggerated due to an artificially large width of an underresolved wave. In the context of reactive fluid flow, a detonation wave consists of a lead shock followed by a thin reaction zone. The curvature effect is determined by the dynamics within the reaction zone; in particular, the competition between a source term for the rate of chemical energy release and a geometric source term due to front curvature. When the width of the reaction zone is small compared with the radius of curvature of the front, the reaction zone can be approximated as quasi-steady and be modeled locally by a system of ordinary differential equations (ODEs) in which the front curvature enters as a parameter. Front curvature breaks the Galilean invariance and the quasi-steady approximation is only valid in a distinguished frame determined by both the front curvature and the tangential velocity divergence ahead of the front. The quasi-steady ODEs determine the reaction zone profile to leading order and can be viewed as an extension of the ZND model. However, the source terms in the ODEs lead to modified Hugoniot jump conditions that take into account front curvature and reaction zone width. When the reaction zone is underresolved, a calculation in effect reduces to a 'capturing algorithm' in which the burn model plays an analogous role for detonation waves as artificial viscosity does for shock waves. In particular, the reaction zone has an unphysically large width that is proportional to the cell size. As a consequence of the modified jump conditions, the numerical reaction zone width gives rise to an artificial curvature effect. This causes numerical solutions to depend on the cell size and orientation of a detonation front relative to the grid. Two algorithms that eliminate numerical curvature effects are discussed, detonation shock dynamics and front tracking.",
author = "Ralph Menikoff and Klaus Lackner and Bukiet, {Bruce G.}",
year = "1996",
month = "2",
doi = "10.1016/0010-2180(95)00106-9",
language = "English (US)",
volume = "104",
pages = "219--240",
journal = "Combustion and Flame",
issn = "0010-2180",
publisher = "Elsevier Inc.",
number = "3",

}

TY - JOUR

T1 - Modeling flows with curved detonation waves

AU - Menikoff, Ralph

AU - Lackner, Klaus

AU - Bukiet, Bruce G.

PY - 1996/2

Y1 - 1996/2

N2 - The nonzero width of a curved detonation front has a significant effect on its propagation. Physically, the curvature effect is a small but important correction. In numerical simulations, the curvature effect tends to be greatly exaggerated due to an artificially large width of an underresolved wave. In the context of reactive fluid flow, a detonation wave consists of a lead shock followed by a thin reaction zone. The curvature effect is determined by the dynamics within the reaction zone; in particular, the competition between a source term for the rate of chemical energy release and a geometric source term due to front curvature. When the width of the reaction zone is small compared with the radius of curvature of the front, the reaction zone can be approximated as quasi-steady and be modeled locally by a system of ordinary differential equations (ODEs) in which the front curvature enters as a parameter. Front curvature breaks the Galilean invariance and the quasi-steady approximation is only valid in a distinguished frame determined by both the front curvature and the tangential velocity divergence ahead of the front. The quasi-steady ODEs determine the reaction zone profile to leading order and can be viewed as an extension of the ZND model. However, the source terms in the ODEs lead to modified Hugoniot jump conditions that take into account front curvature and reaction zone width. When the reaction zone is underresolved, a calculation in effect reduces to a 'capturing algorithm' in which the burn model plays an analogous role for detonation waves as artificial viscosity does for shock waves. In particular, the reaction zone has an unphysically large width that is proportional to the cell size. As a consequence of the modified jump conditions, the numerical reaction zone width gives rise to an artificial curvature effect. This causes numerical solutions to depend on the cell size and orientation of a detonation front relative to the grid. Two algorithms that eliminate numerical curvature effects are discussed, detonation shock dynamics and front tracking.

AB - The nonzero width of a curved detonation front has a significant effect on its propagation. Physically, the curvature effect is a small but important correction. In numerical simulations, the curvature effect tends to be greatly exaggerated due to an artificially large width of an underresolved wave. In the context of reactive fluid flow, a detonation wave consists of a lead shock followed by a thin reaction zone. The curvature effect is determined by the dynamics within the reaction zone; in particular, the competition between a source term for the rate of chemical energy release and a geometric source term due to front curvature. When the width of the reaction zone is small compared with the radius of curvature of the front, the reaction zone can be approximated as quasi-steady and be modeled locally by a system of ordinary differential equations (ODEs) in which the front curvature enters as a parameter. Front curvature breaks the Galilean invariance and the quasi-steady approximation is only valid in a distinguished frame determined by both the front curvature and the tangential velocity divergence ahead of the front. The quasi-steady ODEs determine the reaction zone profile to leading order and can be viewed as an extension of the ZND model. However, the source terms in the ODEs lead to modified Hugoniot jump conditions that take into account front curvature and reaction zone width. When the reaction zone is underresolved, a calculation in effect reduces to a 'capturing algorithm' in which the burn model plays an analogous role for detonation waves as artificial viscosity does for shock waves. In particular, the reaction zone has an unphysically large width that is proportional to the cell size. As a consequence of the modified jump conditions, the numerical reaction zone width gives rise to an artificial curvature effect. This causes numerical solutions to depend on the cell size and orientation of a detonation front relative to the grid. Two algorithms that eliminate numerical curvature effects are discussed, detonation shock dynamics and front tracking.

UR - http://www.scopus.com/inward/record.url?scp=0029914236&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0029914236&partnerID=8YFLogxK

U2 - 10.1016/0010-2180(95)00106-9

DO - 10.1016/0010-2180(95)00106-9

M3 - Article

VL - 104

SP - 219

EP - 240

JO - Combustion and Flame

JF - Combustion and Flame

SN - 0010-2180

IS - 3

ER -