Effect of Spatial Heterogeneity on Near-Limit Propagation of a Stable Detonation
The effect of introducing a spatial heterogeneity into an explosive medium is studied computationally by examining the detonation velocity near the limit to propagation in a thin explosive layer. The explosive system studied is an ideal gas with a single exothermic reaction governed by a pressure-dependent reaction rate ($p^n$) with a pressure exponent of $n = 3$. A pressure-dependent reaction rate, rather than the exponential dependence of reaction on temperature of Arrhenius kinetics, is used so that the detonation wave is stable in the homogeneous case and can be modelled with simple, analytical techniques, and thus the effect of introducing heterogeneity can be clearly identified. The two-dimensional slab of explosive is bounded by a layer of inert gas with the same thermodynamic properties as the explosive. The heterogeneity is introduced into the explosive via a large-amplitude, two-dimensional sinusoidal ripple in density in the initialization of the simulation, while maintaining a constant pressure. The computational simulations are initialized with a ZND solution for the ideal CJ detonation, and the detonation is allowed to propagate into the explosive layer. The simulations show that the detonation in the heterogeneous media exhibits a cellular-like structure of complex shock interactions. The detonation is able to propagate into a significantly thinner layer of explosive and can exhibit a greater velocity than the corresponding homogeneous case. A parametric study of varying the wavelength of the sinusoid shows the existence of an optimal size of heterogeneity at which the favorable effect is the greatest corresponding to a wavelength that is approximately 10 to 50 times the half reaction zone length of the ideal CJ detonation. As the wavelength is decreased to the size of the reaction zone length, the behavior reverts back to the homogeneous case.