Stability of detonations for an idealized condensed-phase model

• Published in: Journal of fluid mechanics Vol. 595; pp. 45 - 82
• Main Authors: SHORT, M., ANGUELOVA, I. I., ASLAM, T. D., BDZIL, J. B., HENRICK, A. K., SHARPE, G. J.
• Format: Journal Article
• Language: English
• Published: Cambridge, UK Cambridge University Press 25.01.2008
• Subjects: Applied sciences; Chemical industry and chemicals; Exact sciences and technology; Fluid dynamics; Industrial chemicals; Powders, propellants, explosives; Stability analysis; Velocity; Detonation; Wave propagation; Stability; Liquid explosive; Solid explosive; Travelling wave; Algorithm; Modeling; Shock wave
• ISSN: 0022-1120; 1469-7645
• DOI: 10.1017/S0022112007008750

The stability of travelling wave Chapman–Jouguet and moderately overdriven detonations of Zeldovich–von Neumann–Döring type is formulated for a general system that incorporates the idealized gas and condensed-phase (liquid or solid) detonation models. The general model consists of a two-component mixture with a one-step irreversible reaction between reactant and product. The reaction rate has both temperature and pressure sensitivities and has a variable reaction order. The idealized condensed-phase model assumes a pressure-sensitive reaction rate, a constant-γ caloric equation of state for an ideal fluid, with the isentropic derivative γ=3, and invokes the strong shock limit. A linear stability analysis of the steady, planar, ZND detonation wave for the general model is conducted using a normal-mode approach. An asymptotic analysis of the eigenmode structure at the end of the reaction zone is conducted, and spatial boundedness (closure) conditions formally derived, whose precise form depends on the magnitude of the detonation overdrive and reaction order. A scaling analysis of the transonic flow region for Chapman–Jouguet detonations is also studied to illustrate the validity of the linearization for Chapman–Jouguet detonations. Neutral stability boundaries are calculated for the idealized condensed-phase model for one- and two-dimensional perturbations. Comparisons of the growth rates and frequencies predicted by the normal-mode analysis for an unstable detonation are made with a numerical solution of the reactive Euler equations. The numerical calculations are conducted using a new, high-order algorithm that employs a shock-fitting strategy, an approach that has significant advantages over standard shock-capturing methods for calculating unstable detonations. For the idealized condensed-phase model, nonlinear numerical solutions are also obtained to study the long-time behaviour of one- and two-dimensional unstable Chapman–Jouguet ZND waves.