A computational study of the interaction of gaseous detonations with a compressible layer

Abstract

The propagation of two-dimensional cellular gaseous detonation bounded by an inert layer is examined via computational simulations. The analysis is based on the high-order integration of the reactive Euler equations with a one-step irreversible reaction. To assess whether the cellular instabilities have a significant influence on a detonation yielding confinement, we achieved numerical simulations for several mixtures from very stable to mildly unstable. The cell regularity was controlled through the value of the activation energy, while keeping constant the ideal Zel’dovich - von Neumann - Döring (ZND) half-reaction length. For stable detonations, the detonation velocity deficit and structure are in accordance with the generalized ZND model, which incorporates the losses due to the front curvature. The deviation with this laminar solution is clear as the activation energy is more significant, increasing the flow field complexity, the variations of the detonation velocity, and the transverse wave strength. The chemical length scale gets thicker, as well as the hydrodynamic thickness. The sonic location is delayed due to the presence of hydrodynamic fluctuations, for which the intensity is increased with the activation energy as well as with the losses to a lesser extent. The flow field has been studied through numerical soot foils, detonation velocities, and 2D detonation front profiles, which are consistent with experimental findings. The velocity deficit increases with the cell irregularity. Moreover, the relation between the detonation limits obtained numerically and in detonation experiments with losses is discussed.