The Physics, Chemistry and Dynamics of Explosions: Doi:10.1098/rsta.2011.0385
The Physics, Chemistry and Dynamics of Explosions: Doi:10.1098/rsta.2011.0385
I N T RO D U C T I O N
1. Introduction
Explosions are ubiquitous in the Universe. Beginning with the Big Bang, the
history of the Universe has been determined largely by explosions, which disrupt
orderly progression and initiate evolution in new directions. While intergalactic
and galactic explosions are of many different types, explosions on the Earth
involve a narrower range of phenomena that are somewhat more amenable to
human understanding. Effects of explosions on humankind, both beneficial and
detrimental, have been investigated from the beginning of science, resulting
today in a degree of knowledge of their properties—their physics, chemistry
and dynamics. During the past century, this knowledge has been used both to
design explosions for specific purposes and to teach us how to avoid those that
are harmful. Any knowledge, however, can be used or misused, for actual or
misconceived benefit. Possible future advances include mitigation of explosion
hazards, improved and new combustion engines and use of inertial confinement
for the development of nuclear fusion as an energy source.
One contribution of 12 to a Theme Issue ‘The physics, chemistry and dynamics of explosions’.
2. Fundamentals
The basis of our present understanding began with the experimental observations
of Berthelot & Vieille [1] and Mallard & LeChatelier [2], and the theoretical
explanations of Chapman [3] and Jouguet [4], more than a century ago. Notable
advances were made by Zel’dovich [5], von Neumann [6] and Döring [7], among
many others, who identified elements of the structure (the ZND structure) of
detonations, which are the stronger of the two classical types of explosions.
The names of these nine scientists have been thoroughly associated with the
fundamental concepts that they developed. These concepts identify two distinct
steady-flow regimes, to which time-dependent explosions may evolve, or between
which time-dependent explosive transitions may occur.
It is instructive to consider the pressure–volume (p–v) diagram of a fluid as
shown in figure 1. The pressure p0 and specific volume v0 of the initial mixture
are indicated as a point in this diagram. For steady, planar, one-dimensional flow,
principles of conservation of mass, momentum and energy can be used to find the
locus of possible final states of the fluid after the heat release has occurred in
an explosion. That locus, called the Hugoniot curve, as illustrated in the figure,
consists of two disconnected parts, an upper branch and a lower branch. The
two branches must be separated for exothermic processes because application
of steady mass and momentum conservation alone requires that a straight line
p
upper branch
Hugoniot curve
v0 v
Figure 1. Traditional pressure–volume (p–v) diagram showing the Hugoniot curve with possible
states of deflagrations and detonations.
connecting the initial and final states in this diagram, called the Rayleigh line
(known in gas dynamics as a trajectory for constant-area compressible flow with
heat addition), cannot have a positive slope. When the final state lies along the
lower branch of the Hugoniot curve, the process is called a deflagration, and when
it lies along the upper branch, it is a detonation.
Thus, there are two very different and distinct types of one-dimensional,
steady reaction waves propagating through homogeneous reactive material,
consistent with conservation conditions. The figure illustrates that detonations
exhibit large pressure increases and small volume decreases, while deflagrations
have large volume increases and small pressure decreases. Propagation velocities
of such fronts are proportional to the square root of the negative of the
slope of the Rayleigh line, two representative examples of which are indicated
on the figure. Propagation velocities of detonations are thus seen to exceed
those of deflagrations; detonations are found to be supersonic and deflagrations
subsonic. The figure also shows that there is a minimum velocity for detonations
and a maximum velocity for deflagrations, limiting conditions that correspond to
tangency of the Rayleigh line and the Hugoniot curve, and the associated final
conditions are called the (upper and lower) Chapman–Jouguet (CJ ) points of
the curve. Further analysis demonstrates that final fluid velocities are sonic at
the CJ points.
Except for CJ conditions, when there is any intersection at all between a
given Rayleigh line and the Hugoniot curve, there are two, as can be seen in
the figure. Therefore, there are two different types of each process. From the
slopes, it is seen that both types of deflagrations travel slower than CJ, and
both types of detonations faster. The first intersection is called weak and the
second strong. Final fluid velocities are supersonic for strong deflagrations and
for weak detonations, while they are subsonic in the other two cases. Although
strong deflagrations have been postulated for certain galactic winds, all known
deflagrations on the Earth are weak, approaching CJ only for highly turbulent,
energetic flames. On the other hand, most detonations are CJ, although both
weak and strong detonations occur. Strong detonations often arise for periods of
time in transitions from deflagrations to detonations (DDTs). Laboratory studies
of chemical kinetics in shock tubes, when exothermic overall, in principle also
are strong detonations for the entire shock-tube test time. On the other hand,
condensation ‘shocks’ in wind tunnels are weak detonations.
Explosions, however, are dynamic processes that often are neither the classical
deflagrations nor the classical detonations described by steady-state flows
in homogeneous backgrounds with perfectly unobtrusive geometries. As an
example, consider the preceding conclusion that detonations are supersonic and
deflagrations subsonic; this will be seen in a later article to not always be true
for porous media. Although pressure changes associated with deflagrations are
smaller than those associated with detonations, the expansion and heating that
occur in deflagrations cause pressure build-up, especially in confined spaces.
This can be quite damaging and thereby causes the processes to be classified
as explosions, even though detonations do not develop. In many instances, these
pressure waves and their associated heating produce DDT. The most damaging
explosions, however, are generally detonations because they have the highest
pressures. Pressures even higher than those in CJ detonations can occur in
transient processes, such as shock–shock interactions and strong detonations,
which precede the completion of DDT. Depending on the specific configuration
and both the reactivity and the energetics of the combustible mixtures, periods
may occur in which there is supersonic front propagation slower than CJ, and
such slow (sometimes called ‘weak’) detonations are types of explosions that may
not transform fully to CJ and that cannot be classical steadily propagating fronts.
Nevertheless, they may be very damaging and therefore need to be understood.
There thus exist many non-classical types of explosions requiring further study.
The first article in this issue presents an examination of one of the largest
recent accidental chemical explosions. The explosion at Buncefield, UK, in 2005
is unique in that it is unusual to have such a thorough investigation of a large-
scale, accidental, statistically improbable event. Questions raised by the forensics
committee included the following. What caused the fuel spillage? How did it
ignite? Do we know enough to understand the dynamics of the turbulent flame
development and propagation? What was the role of the obstacles in the path of
the explosion front? And finally, did a detonation develop at any stage? This first
article, by Bradley et al. [8], summarizes the work of the technical investigation
committee, presents the current state of the analysis of this event, puts it in
the context of other historic large-scale accidental explosions and then describes
where future research is required. This needed research, such as the research
reported in the subsequent articles in this issue, goes well beyond the idealizations
identified earlier.
The ZND detonation structure envisions steady, planar flow with heat release
initiated by the temperature rise across a strong planar shock wave (as in
many shock-tube experiments). Evidence that this picture is much too simplified
for most real-world explosions can be at least traced to the observations
of spinning detonations by Campbell & Woodhead [9], and clear recognition
of the prevalent existence of multi-dimensional structures was provided by
Denisov & Troshin [10]. We now know that the ZND structure is generally
unstable, evolving into cellular detonation structures that range from regular
diamond patterns to very complex dynamic interactions of shock waves, reaction
regions and shear layers behind leading shock fronts. In the second article,
Kessler et al. [11] use multi-dimensional numerical simulations to study these
phenomena. By extending the computations to detonation propagation normal
to a gradient of stoichiometry, they identify upstream (behind the detonation
front) reaction regions into which reactants diffuse from opposite sides. This
region is bounded by rich and lean reaction regions, and the result is a triple-
flame detonation wave. This type of gradient of fuel concentration could develop
in the more complex spatial and temporal fuel distributions that occur in
explosion-prone situations, such as in tunnels in coal mines or during industrial
fuel leakages.
Complementing this computational work are the analytical studies in the third
article; Clavin & Williams [12] present a unified formulation that covers a range
of detonation topics, including initiation, propagation, quenching, pulsation and
cellular structure. Asymptotic analyses are discussed here for the limit of strongly
temperature-sensitive rates of heat release and for the Newtonian limit in which
the difference between constant-pressure and constant-volume heat capacities
vanishes (assumed by Isaac Newton and found to be surprisingly accurate for
high-temperature real gases in hypersonics). Planar detonations are shown to be
unstable to shear-wave disturbances, even when the heat-release rate is insensitive
to temperature, with a weakly nonlinear analysis indicating that this bifurcation
leads to distinct diamond patterns. This article also clarifies conflicting literature
concerning the stability of shock waves in general non-ideal fluids, deriving for
the first time the boundary between decaying and non-radiating, non-decaying
stable shock discontinuities.
The next two articles are theoretical and experimental investigations,
respectively, of explosions in porous media. The first, by Brailovsky et al. [13],
summarizes implications of a simplified mathematical model based on the concept
of hydraulic resistance. Introduction of the hydraulic-resistance approximation
reduces the order of the differential equations that describe the detonation and
deflagration processes, thereby facilitating their solution and making it easier
to understand their predictions. For these reasons, it is important to determine
experimentally the extent to which the model can correspond to reality, and
a new favourable comparison is identified (called Shchelkin’s effect there). The
experimental observations and measurements reported by Ciccarelli [14], in the
second of these articles, expose experimental similarities and differences between
explosion propagation in porous media and in chambers filled with obstacles.
These results are quite revealing; they show, in particular, structures of explosions
in the low-velocity slow-detonation regime, where the propagation velocities of
the fronts are between cs and the CJ detonation velocity. It is made clear
here that acquisition of detailed, local data on explosions in porous media is
of the rate of expansion and size of the Universe. Understanding why SNIa are so
uniform involves understanding the evolution of turbulent flames and detonations
in the thermonuclear environment in the star, which certainly is a challenging
topic in explosions.
25
1/d (m–1) 20
15
10
0
4 6 8 10 12 14 16
% CH4
Figure 2. Measured lean and rich detonation limits for methane in air as a function of the per cent
of methane by volume. Limits are shown for four values of 1/d, the inverse of the diameter of the
test device. Names and symbols refer to earlier studies (squares, Kogarko [20]; triangles, Matsui
[21]; circles, Gamezo et al. [22]).
Relief valves and blow-out caps are common explosion-control devices. There have
been studies of detonation arrestors, for example, employing narrow channels,
porous materials, dust dispersions, water sprays, water-filled passages, etc. There
have been methods for extinguishing explosions that attempt to disperse inert
rock dust or placing buckets of water to be overturned during an explosion. But
the task is confusing and difficult. Methods that appear to damp explosions under
some conditions enhance them under other conditions. We would like to be able to
control transients, to change the explosion state, so that we can ignite, enhance or
quench at will, turn detonations into deflagrations or taylor the process of DDT.
Much more understanding is necessary to achieve these objectives.
Finally, more research is needed on the role of turbulence in explosions.
A half-century ago, when the ZND detonation structure was widely accepted
as universal, fine-scale shadowgraph observations surprisingly suggested that
detonations were turbulent. Improvements in resolution soon showed that this
turbulence, in reality, was the now-familiar regular cellular diamond pattern.
Today, with greatly improved experimental techniques, especially for detonations
having highly irregular, somewhat chaotic cell structures, reacting shear layers
behind leading shocks are often described as turbulent. The true nature of
this kind of turbulence, however, remains to be determined. As high-speed
deflagrations interact with shock waves and boundary layers, they are subject
to shock-related and shear-flow instabilities that may generate turbulence of
character different from that which is commonly understood. More thought is
thus warranted about turbulence in all types of explosions.
There is now a great deal of knowledge about turbulence in fluids, which has
been employed to construct diagrams for turbulent combustion, with different
types of combustion processes reasoned to occur in different parts of the diagrams.
One selection of coordinates for such a diagram is a turbulence Reynolds number
and a turbulence Damköhler number, the ratio of a turbulence time scale to
a chemical time scale [23]. Turbulent combustion is considered to occur in a
distributed manner at low Damköhler numbers and in sheets at high Damköhler
numbers, for example. Questions arise, however, as to whether such a classification
is sufficient for explosions.
In a recent computational study showing DDT in unconfined systems [24],
compressibility effects suggest that a Mach number should be considered as
a coordinate for classifying regimes of turbulent combustion. This possibility
remains to be explored. Turbulent combustion computations for SNIa employ
classical turbulence concepts, but in that application, the associated Lewis
number, the ratio of the heat diffusivity to the reactant diffusivity, is very large, of
the order of 1000 or more. Even laminar deflagrations, however, are known to be
unstable to pulsating instabilities at large Lewis numbers, approaching a pulsating
type of chaos, a turbulence quite different from that known for fluids, for Lewis
numbers greater than about 10, thereby calling into question the basis of the
standard candle. Consideration therefore perhaps should be given to including
the Lewis number as a coordinate for diagrams of turbulent combustion. Thus,
much remains to be learned about turbulence in explosions.
5. Concluding comment
References
1 Berthelot, M. & Vieille, P. 1881 On the propagation velocity of explosion phenomena. Acad.
Sci. Comptes Rendus Math. 93, 18–22.
2 Mallard, E. & LeChatelier, H. L. 1881 On the flame propagation velocity in explosive gas
mixtures. Acad. Sci. Comptes Rendus Math. 93, 145–148.
3 Chapman, D. L. 1899 On the rate of explosion in gases. Philos. Mag. 47, 90–104.
4 Jouguet, E. 1905 On the propagation of chemical reactions in gases. J. Math. Pures Appl. 6,
347–425.
5 Zel’dovich, Y. B. 1940 On the theory of the propagation of detonation in gaseous systems.
J. Exp. Theor. Phys. 10, 542–568.
6 von Neumann, J. 1942 Theory of detonation waves. Report no. 549, Progress Report no. 238
(April), Office of Scientific Research and Development (OSRD), Washington, DC.
7 Döring, W. 1943 On the detonation process in gases. Ann. Phys. 43, 421–436.
8 Bradley, D., Chamberlain, G. A. & Drysdale, D. D. 2012 Large vapour cloud explosions, with
particular reference to that at Buncefield. Phil. Trans. R. Soc. A 370, 544–566. (doi:10.1098/
rsta.2011.0419)
9 Campbell, C. & Woodhead, D. W. 1926 The ignition of gases an explosive wave. J. Chem. Soc.
2, 3010–3021. (doi:10.1039/jr9262903010)
10 Denisov, Y. N. & Troshin, Y. K. 1959 Pulsating and spinning detonation of gaseous mixtures
in tubes. Dokl. Akad. Nauk. 125, 110–113.
11 Kessler, D. A., Gamezo, V. N. & Oran, E. S. 2012 Gas-phase detonation propagation in mixture
composition gradients. Phil. Trans. R. Soc. A 370, 567–596. (doi:10.1098/rsta.2011.0342)
12 Clavin, P. & Williams, F. A. 2012 Analytical studies of the dynamics of gaseous detonations.
Phil. Trans. R. Soc. A 370, 597–624. (doi:10.1098/rsta.2011.0345)
13 Brailovsky, I., Kagan, L. & Sivashinsky, G. 2012 Combustion waves in hydraulically resisted
systems. Phil. Trans. R. Soc. A 370, 625–646. (doi:10.1098/rsta.2011.0341)
14 Ciccarelli, G. 2012 Explosion propagation in inert porous media. Phil. Trans. R. Soc. A 370,
647–667. (doi:10.1098/rsta.2011.0346)
15 Chiquete, C. & Tumin, A. 2012 Stability of detonation in a circular pipe with porous walls.
Phil. Trans. R. Soc. A 370, 668–688. (doi:10.1098/rsta.2011.0343)
16 Bradley, D. 2012 Autoignitions and detonations in engines and ducts. Phil. Trans. R. Soc. A
370, 689–714. (doi:10.1098/rsta.2011.0367)
17 Thomas, G. 2012 Some observations on the initiation and onset of detonation. Phil. Trans. R.
Soc. A 370, 715–739. (doi:10.1098/rsta.2011.0368)
18 Starikovskiy, A., Aleksandrov, N. & Rakitin, A. 2012 Plasma-assisted ignition and deflagration-
to-detonation transition. Phil. Trans. R. Soc. A 370, 740–773. (doi:10.1098/rsta.2011.0344)
19 Wheeler, J. C. 2012 Astrophysical explosions: from solar flares to cosmic gamma-ray bursts.
Phil. Trans. R. Soc. A 370, 774–799. (doi:10.1098/rsta.2011.0351)
20 Kogarko, S. M. 1959 Detonation of methane–air mixtures and the detonation limits of
hydrocarbon–air mixtures in a large-diameter pipe. Sov. Phys. Tech. Phys. 3, 1904–1916.
21 Matsui, H. 2002 Detonation propagation limits in homogeneous and heterogeneous systems.
J. Phys. IV 12, 11–17. (doi:10.1051/jp4:20020262)
22 Gamezo, V. N. et al. In press. Detonability of natural gas–air mixtures. Combust. Flame.
(doi:10.1016/j.combustflame.2011.08.009)
23 Williams, F. A. 2001 Some recent studies in turbulent combustion. In Smart control of turbulent
combustion (ed. A. Yoshida), pp. 1–11. Berlin, Germany: Springer.
24 Poludnenko, A. Y., Gardner, T. A. & Oran, E. S. 2011 Spontaneous transition of turbulent
flames to detonations in unconfined media. Phys. Rev. Lett. 107, 054 501–054 514. (doi:10.1103/
PhysRevLett.107.054501)