Unsteady features on one-dimensional hydrogen-air detonations

Yu Daimon, Akiko Matsuo

Research output: Contribution to journalArticle

14 Citations (Scopus)

Abstract

The features of one-dimensional unsteady detonations are studied numerically using a hydrogen-air detailed chemical reaction model. A series of simulations are carried out while degree of overdrive, initial pressure, initial temperature, and equivalence ratio are varied. The oscillation modes and mechanisms of the one-dimensional detonations are discussed with reference to shock pressure histories and x-t diagrams of density distributions. As the degree of overdrive is reduced with a stoichiometric mixture of hydrogen-air at P0=0.421 atm and T0=293 K, a steady state appears, along with a high-frequency mode and a low-frequency mode. The oscillation mechanism of the high-frequency mode is the same as that of the regular regime of unsteady shock-induced combustion observed around a spherical projectile flying at hypersonic velocity in detonable gases. The degree of overdrive threshold between the steady and unsteady region increases monotonically with initial pressure and decreases monotonically with initial temperature. When the equivalence ratio is changed, the threshold has a minimum value around φ=1. We focus attention on a nondimensional effective activation energy, which is generally used for linear stability analysis. The oscillation mode depends highly on the nondimensional effective activation energy. The oscillation of the detonation front appears as the nondimensional effective activation energy goes past a threshold value of 5.2. Furthermore, we investigate the failed regime and possible reignition in this regime. In the failed regime, a detonation wave breaks up into a leading shock, a contact discontinuity, and a rarefaction wave. When the shock is weak, reignition time becomes very long. Therefore, the reignition after the failed regime is difficult to reproduce in the restricted computational domain. The reignition process in the failed regime is investigated by means of analysis consisting of integration along the point of intersection between a Rayleigh line for weak leading shock and a partially burnt Hugoniot curve. The reignition time increases dramatically with decreasing temperature behind the shock wave, when the gas condition goes past the second explosion limit. The second explosion limit is one of the characteristics of the hydrogen-air detailed chemical reaction model and does not exist in the one-step chemical reaction model. Lastly, the reignition time obtained by the analysis is compared with that obtained by the simulation results. The simulation results agree well with the analytical results.

Original languageEnglish
Article number116101
JournalPhysics of Fluids
Volume19
Issue number11
DOIs
Publication statusPublished - 2007 Nov

Fingerprint

Detonation
detonation
ignition
Hydrogen
Chemical reactions
shock
Activation energy
air
hydrogen
Air
Explosions
chemical reactions
oscillations
Gases
activation energy
Linear stability analysis
thresholds
equivalence
explosions
Hypersonic aerodynamics

ASJC Scopus subject areas

  • Mechanics of Materials
  • Computational Mechanics
  • Physics and Astronomy(all)
  • Fluid Flow and Transfer Processes
  • Condensed Matter Physics

Cite this

Unsteady features on one-dimensional hydrogen-air detonations. / Daimon, Yu; Matsuo, Akiko.

In: Physics of Fluids, Vol. 19, No. 11, 116101, 11.2007.

Research output: Contribution to journalArticle

@article{f9c51ed12f72426e81d8d72348e70b1b,
title = "Unsteady features on one-dimensional hydrogen-air detonations",
abstract = "The features of one-dimensional unsteady detonations are studied numerically using a hydrogen-air detailed chemical reaction model. A series of simulations are carried out while degree of overdrive, initial pressure, initial temperature, and equivalence ratio are varied. The oscillation modes and mechanisms of the one-dimensional detonations are discussed with reference to shock pressure histories and x-t diagrams of density distributions. As the degree of overdrive is reduced with a stoichiometric mixture of hydrogen-air at P0=0.421 atm and T0=293 K, a steady state appears, along with a high-frequency mode and a low-frequency mode. The oscillation mechanism of the high-frequency mode is the same as that of the regular regime of unsteady shock-induced combustion observed around a spherical projectile flying at hypersonic velocity in detonable gases. The degree of overdrive threshold between the steady and unsteady region increases monotonically with initial pressure and decreases monotonically with initial temperature. When the equivalence ratio is changed, the threshold has a minimum value around φ=1. We focus attention on a nondimensional effective activation energy, which is generally used for linear stability analysis. The oscillation mode depends highly on the nondimensional effective activation energy. The oscillation of the detonation front appears as the nondimensional effective activation energy goes past a threshold value of 5.2. Furthermore, we investigate the failed regime and possible reignition in this regime. In the failed regime, a detonation wave breaks up into a leading shock, a contact discontinuity, and a rarefaction wave. When the shock is weak, reignition time becomes very long. Therefore, the reignition after the failed regime is difficult to reproduce in the restricted computational domain. The reignition process in the failed regime is investigated by means of analysis consisting of integration along the point of intersection between a Rayleigh line for weak leading shock and a partially burnt Hugoniot curve. The reignition time increases dramatically with decreasing temperature behind the shock wave, when the gas condition goes past the second explosion limit. The second explosion limit is one of the characteristics of the hydrogen-air detailed chemical reaction model and does not exist in the one-step chemical reaction model. Lastly, the reignition time obtained by the analysis is compared with that obtained by the simulation results. The simulation results agree well with the analytical results.",
author = "Yu Daimon and Akiko Matsuo",
year = "2007",
month = "11",
doi = "10.1063/1.2801478",
language = "English",
volume = "19",
journal = "Physics of Fluids",
issn = "1070-6631",
publisher = "American Institute of Physics Publising LLC",
number = "11",

}

TY - JOUR

T1 - Unsteady features on one-dimensional hydrogen-air detonations

AU - Daimon, Yu

AU - Matsuo, Akiko

PY - 2007/11

Y1 - 2007/11

N2 - The features of one-dimensional unsteady detonations are studied numerically using a hydrogen-air detailed chemical reaction model. A series of simulations are carried out while degree of overdrive, initial pressure, initial temperature, and equivalence ratio are varied. The oscillation modes and mechanisms of the one-dimensional detonations are discussed with reference to shock pressure histories and x-t diagrams of density distributions. As the degree of overdrive is reduced with a stoichiometric mixture of hydrogen-air at P0=0.421 atm and T0=293 K, a steady state appears, along with a high-frequency mode and a low-frequency mode. The oscillation mechanism of the high-frequency mode is the same as that of the regular regime of unsteady shock-induced combustion observed around a spherical projectile flying at hypersonic velocity in detonable gases. The degree of overdrive threshold between the steady and unsteady region increases monotonically with initial pressure and decreases monotonically with initial temperature. When the equivalence ratio is changed, the threshold has a minimum value around φ=1. We focus attention on a nondimensional effective activation energy, which is generally used for linear stability analysis. The oscillation mode depends highly on the nondimensional effective activation energy. The oscillation of the detonation front appears as the nondimensional effective activation energy goes past a threshold value of 5.2. Furthermore, we investigate the failed regime and possible reignition in this regime. In the failed regime, a detonation wave breaks up into a leading shock, a contact discontinuity, and a rarefaction wave. When the shock is weak, reignition time becomes very long. Therefore, the reignition after the failed regime is difficult to reproduce in the restricted computational domain. The reignition process in the failed regime is investigated by means of analysis consisting of integration along the point of intersection between a Rayleigh line for weak leading shock and a partially burnt Hugoniot curve. The reignition time increases dramatically with decreasing temperature behind the shock wave, when the gas condition goes past the second explosion limit. The second explosion limit is one of the characteristics of the hydrogen-air detailed chemical reaction model and does not exist in the one-step chemical reaction model. Lastly, the reignition time obtained by the analysis is compared with that obtained by the simulation results. The simulation results agree well with the analytical results.

AB - The features of one-dimensional unsteady detonations are studied numerically using a hydrogen-air detailed chemical reaction model. A series of simulations are carried out while degree of overdrive, initial pressure, initial temperature, and equivalence ratio are varied. The oscillation modes and mechanisms of the one-dimensional detonations are discussed with reference to shock pressure histories and x-t diagrams of density distributions. As the degree of overdrive is reduced with a stoichiometric mixture of hydrogen-air at P0=0.421 atm and T0=293 K, a steady state appears, along with a high-frequency mode and a low-frequency mode. The oscillation mechanism of the high-frequency mode is the same as that of the regular regime of unsteady shock-induced combustion observed around a spherical projectile flying at hypersonic velocity in detonable gases. The degree of overdrive threshold between the steady and unsteady region increases monotonically with initial pressure and decreases monotonically with initial temperature. When the equivalence ratio is changed, the threshold has a minimum value around φ=1. We focus attention on a nondimensional effective activation energy, which is generally used for linear stability analysis. The oscillation mode depends highly on the nondimensional effective activation energy. The oscillation of the detonation front appears as the nondimensional effective activation energy goes past a threshold value of 5.2. Furthermore, we investigate the failed regime and possible reignition in this regime. In the failed regime, a detonation wave breaks up into a leading shock, a contact discontinuity, and a rarefaction wave. When the shock is weak, reignition time becomes very long. Therefore, the reignition after the failed regime is difficult to reproduce in the restricted computational domain. The reignition process in the failed regime is investigated by means of analysis consisting of integration along the point of intersection between a Rayleigh line for weak leading shock and a partially burnt Hugoniot curve. The reignition time increases dramatically with decreasing temperature behind the shock wave, when the gas condition goes past the second explosion limit. The second explosion limit is one of the characteristics of the hydrogen-air detailed chemical reaction model and does not exist in the one-step chemical reaction model. Lastly, the reignition time obtained by the analysis is compared with that obtained by the simulation results. The simulation results agree well with the analytical results.

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

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

U2 - 10.1063/1.2801478

DO - 10.1063/1.2801478

M3 - Article

AN - SCOPUS:36849081221

VL - 19

JO - Physics of Fluids

JF - Physics of Fluids

SN - 1070-6631

IS - 11

M1 - 116101

ER -