Steady non-ideal detonations in cylindrical sticks of expolsives

Numerical simulations of detonations in cylindrical rate-sticks of highly non-ideal explosives are performed, using a simple model with a weakly pressure dependent rate law and a pseudo-polytropic equation of state. Some numerical issues with such simulations are investigated, and it is shown that very high resolution (hundreds of points in the reaction zone) are required for highly accurate (converged) solutions. High resolution simulations are then used to investigate the qualitative dependences of the detonation driving zone structure on the diameter and degree of confinement of the explosive charge. The simulation results are used to show that, given the radius of curvature of the shock at the charge axis, the steady detonation speed and the axial solution are accurately predicted by a quasi-one-dimensional theory, even for cases where the detonation propagates at speeds significantly below the Chapman-Jouguet speed. Given reaction rate and equation of state models, this quasi-one-dimensional theory offers a significant improvement to Wood-Kirkwood theories currently used in industry

A numerical study of detonation diffraction

An investigation of detonation diffraction through an abrupt area change has been carried out via a set of two-dimensional numerical simulations parameterized by the activation energy of the reactant. Our analysis is specialized to a reactive mixture with a perfect gas equation of state and a single-step reaction in the Arrhenius form. Lagrangian particles are injected into the flow as a diagnostic tool for identifying the dominant terms in the equation that describes the temperature rate of change of a fluid element, expressed in the shock-based reference system. When simplified, this equation provides insight into the competition between the energy release rate and the expansion rate behind the diffracting front. The mechanism of spontaneous generation of transverse waves along the diffracting front is carefully analysed and related to the sensitivity of the reaction rate to temperature. We study in detail three highly resolved cases of detonation diffraction that illustrate different types of behaviour, super-, sub- and near-critical diffraction

The role of unsteadiness in direct initiation of gaseous detonations

An analytical model is presented for the direct initiation of gaseous detonations by a blast wave. For stable or weakly unstable mixtures, numerical simulations of the spherical direct initiation event and local analysis of the one-dimensional unsteady reaction zone structure identify a competition between heat release, wave front curvature and unsteadiness. The primary failure mechanism is found to be unsteadiness in the induction zone arising from the deceleration of the wave front. The quasi-steady assumption is thus shown to be incorrect for direct initiation. The numerical simulations also suggest a non-uniqueness of critical energy in some cases, and the model developed here is an attempt to explain the lower critical energy only. A critical shock decay rate is determined in terms of the other fundamental dynamic parameters of the detonation wave, and hence this model is referred to as the critical decay rate (CDR) model. The local analysis is validated by integration of reaction-zone structure equations with real gas kinetics and prescribed unsteadiness. The CDR model is then applied to the global initiation problem to produce an analytical equation for the critical energy. Unlike previous phenomenological models of the critical energy, this equation is not dependent on other experimentally determined parameters and for evaluation requires only an appropriate reaction mechanism for the given gas mixture. For different fuel–oxidizer mixtures, it is found to give agreement with experimental data to within an order of magnitude