XiaoCheng Mi

Geometric scaling for a detonation wave governed by a pressure-dependent reaction rate and yielding confinement

The propagation of detonation waves in reactive media bounded by an inert, compressible layer is examined via computational simulations in two different geometries, axisymmetric cylinders, and two dimensional, planar slabs. For simplicity, an ideal gas equation of state is used with a pressure-dependent reaction rate that results in a detonation wave structure that does not exhibit cellular instability. The detonation is initiated as an ideal Chapman-Jouguet (CJ) detonation with a one-dimensional structure, and then allowed to propagate into a finite diameter or thickness layer of explosive surrounded by an inert layer. The yielding confinement of the inert layer results in the detonation wave decaying to a sub-CJ steady state velocity or failing entirely. Simulations are performed with different values of the reaction rate pressure exponent (n = 2 and 3) and different impedance confinement (greater than, less than, and equal to that of the explosive). The velocity decrement and critical dimension (critic...

Influence of discrete sources on detonation propagation in a Burgers equation analog system.

An analog to the equations of compressible flow that is based on the inviscid Burgers equation is utilized to investigate the effect of spatial discreteness of energy release on the propagation of a detonation wave. While the traditional Chapman-Jouguet (CJ) treatment of a detonation wave assumes that the energy release of the medium is homogeneous through space, the system examined here consists of sources represented by δ functions embedded in an otherwise inert medium. The sources are triggered by the passage of the leading shock wave following a delay that is either of fixed period or randomly generated. The solution for wave propagation through a large array (10^{3}-10^{4}) of sources in one dimension can be constructed without the use of a finite difference approximation by tracking the interaction of sawtooth-profiled waves for which an analytic solution is available. A detonation-like wave results from the interaction of the shock and rarefaction waves generated by the sources. The measurement of the average velocity of the leading shock front for systems of both regular, fixed-period and randomized sources is found to be in close agreement with the velocity of the equivalent CJ detonation in a uniform medium, wherein the sources have been spatially homogenized. This result may have implications for the applicability of the CJ criterion to detonations in highly heterogeneous media (e.g., polycrystalline, solid explosives) and unstable detonations with a transient and multidimensional structure (e.g., gaseous detonation waves).

Front roughening of flames in discrete media

The morphology of flame fronts propagating in reactive systems composed of randomly positioned, pointlike sources is studied. The solution of the temperature field and the initiation of new sources is implemented using the superposition of the Greens function for the diffusion equation, eliminating the need to use finite-difference approximations. The heat released from triggered sources diffuses outward from each source, activating new sources and enabling a mechanism of flame propagation. Systems of 40000 sources in a 200×200 two-dimensional domain were tracked using computer simulations, and statistical ensembles of 120 realizations of each system were averaged to determine the statistical properties of the flame fronts. The reactive system of sources is parameterized by two nondimensional values: the heat release time (normalized by interparticle diffusion time) and the ignition temperature (normalized by adiabatic flame temperature). These two parameters were systematically varied for different simulations to investigate their influence on front propagation. For sufficiently fast heat release and low ignition temperature, the front roughness [defined as the root mean square deviation of the ignition temperature contour from the average flame position] grew following a power-law dependence that was in excellent agreement with the Kardar-Parisi-Zhang (KPZ) universality class (β=1/3). As the reaction time was increased, lower values of the roughening exponent were observed, and at a sufficiently great value of reaction time, reversion to a steady, constant-width thermal flame was observed that matched the solution from classical combustion theory. Deviation away from KPZ scaling was also observed as the ignition temperature was increased. The features of this system that permit it to exhibit both KPZ and non-KPZ scaling are discussed.

Effect of spatial discretization of energy on detonation wave propagation

Detonation propagation in the limit of highly spatially discretized energy sources is investigated. The model of this problem begins with a medium consisting of a calorically perfect gas with a prescribed energy release per unit mass. The energy release is collected into sheet-like sources that are now embedded in an inert gas that fills the spaces between them. The release of energy in the first sheet results in a planar blast wave that propagates to the next source, which is triggered after a prescribed delay, generating a new blast, and so forth. The resulting wave dynamics as the front passes through hundreds of such sources is computationally simulated by numerically solving the governing one-dimensional Euler equations in the lab-fixed reference frame. The average wave speed for each simulation is measured once the wave propagation has reached a quasi-periodic solution. Velocities in excess of the CJ speed are found as the sources are made increasingly discrete, with the deviation above CJ being as great as

Propagation of gaseous detonation waves in a spatially inhomogeneous reactive medium

15\%

Propagation of a Pressure-Dependent Detonation with Different Acoustic Impedance Confinements

. The total energy release, delay time, and whether the sources remain lab-fixed or are convected with the flow do not have a significant influence on the deviation of the average wave speed away from CJ. Such continuous waves can also be shown to have a time-averaged structure consistent with the classical ZND structure of a detonation. In the limit of highly discrete sources, temporal averaging of the wave structure shows that the effective sonic surface does not correspond to an equilibrium state. The average state of the flow leaving the wave in this case does eventually reach the equilibrium Hugoniot, but only after the effective sonic surface has been crossed. Thus, the super-CJ waves observed in the limit of highly discretized sources can be understood as weak detonations due to the non-equilibrium state at the effective sonic surface.

Investigation of detonation velocity in heterogeneous explosive system using the reactive Burgers’ analog

Detonation propagation in a compressible medium wherein the energy release has been made spatially inhomogeneous is examined via numerical simulation. The inhomogeneity is introduced via step functions in the reaction progress variable, with the local value of energy release correspondingly increased so as to maintain the same average energy density in the medium, and thus a constant Chapman Jouguet (CJ) detonation velocity. A one-step Arrhenius rate governs the rate of energy release in the reactive zones. The resulting dynamics of a detonation propagating in such systems with one-dimensional layers and two-dimensional squares are simulated using a Godunov-type finite-volume scheme. The resulting wave dynamics are analyzed by computing the average wave velocity and one-dimensional averaged wave structure. In the case of sufficiently inhomogeneous media wherein the spacing between reactive zones is greater than the inherent reaction zone length, average wave speeds significantly greater than the corresponding CJ speed of the homogenized medium are obtained. If the shock transit time between reactive zones is less than the reaction time scale, then the classical CJ detonation velocity is recovered. The spatio-temporal averaged structure of the waves in these systems is analyzed via a Favre averaging technique, with terms associated with the thermal and mechanical fluctuations being explicitly computed. The analysis of the averaged wave structure identifies the super-CJ detonations as weak detonations owing to the existence of mechanical non-equilibrium at the effective sonic point embedded within the wave structure. The correspondence of the super-CJ behavior identified in this study with real detonation phenomena that may be observed in experiments is discussed. ∗ Corresponding author: andrew.higgins@mcgill.ca

Dual-stage ignition of boron particle agglomerates

A series of numerical simulations of a pressure-dependent detonation propagating with losses in charges was conducted to study the effect of the acoustic impedance of confinements (greater than, less than, and equal to that of the explosive). The acoustic impedance of the confinement was varied by inversely changing the temperature and density of the confinement while maintaining the pressure as constant. The results showed that the confinement has a strong influence on the structures of steady detonation waves. By quantitatively comparison, the results showed the ability of higher impedance confinement to enable the detonation wave to propagate at greater velocity for the same charge thickness. With different impedance confinement and pressure exponent values, all the measured shock front curvature and detonation velocity relation agrees with the classic model of Wood and Kirkwood. Additionally, the results also showed that the feasibility of different theoretical models should be carefully validated for different values of acoustic impedance confinement.

Effect of spatial heterogeneity on near-limit propagation of a pressure-dependent detonation

In this study, the applicability of the Chapman-Jouguet (CJ) criterion is tested numerically for heterogeneous explosive media using a simple detonation analog. The analog system consists of a reactive Burgers’ equation coupled with an Arrhenius type reaction wave, and the heterogeneity of the explosive media is mimicked using a discrete energy source approach. The governing equation is solved using a second order, finite-volume approach and the average propagation velocity of the discrete detonation is determined by tracking the leading shock front. Consistent with previous studies, the averaged velocity of the leading shock front from the unsteady numerical simulations is also found to be in good agreement with the velocity of a CJ detonation in a uniform medium wherein the energy source is spatially homogenized. These simulations have thus implications for whether the CJ criterion is valid to predict the detonation velocity in heterogeneous explosive media.