Thermal Structure and Burning Velocity of Flames in Non-volatile
Fuel Suspensions
http://dx.doi.org/10.1016/j.proci.2016.06.043
Michael J. Soo*, Keishi Kumashiro, Samuel Goroshin, David L. Frost, Jeffrey M. Bergthorson
arXiv:1609.00347v1 [physics.flu-dyn] 1 Sep 2016
Department of Mechanical Engineering, McGill University
817 Sherbrooke Street West, Montreal, QC, H3A 0C3, Canada
*Email: michael.soo@mail.mcgill.ca
Abstract
Flame propagation through a non-volatile solid-fuel suspension is studied using a simplified, time-dependent numerical model
that considers the influence of both diffusional and kinetic rates on the particle combustion process. It is assumed that particles react via a single-step, first-order Arrhenius surface reaction with an oxidizer delivered to the particle surface through gas
diffusion. Unlike the majority of models previously developed for flames in suspensions, no external parameters are imposed,
such as particle ignition temperature, combustion time, or the assumption of either kinetic- or diffusion-limited particle combustion regimes. Instead, it is demonstrated that these parameters are characteristic values of the flame propagation problem
that must be solved together with the burning velocity, and that the a priori imposition of these parameters from single-particle
combustion data may result in erroneous predictions. It is also shown that both diffusive and kinetic reaction regimes can
alternate within the same flame and that their interaction may result in non-trivial flame behavior. In fuel-lean mixtures, it is
demonstrated that this interaction leads to certain particle size ranges where burning velocity increases with increasing particle
size, opposite to the expected trend. For even leaner mixtures, the interplay between kinetic and diffusive reaction rates leads to
the appearance of a new type of flame instability where kinetic and diffusive regimes alternate in time, resulting in a pulsating
regime of flame propagation.
Keywords
Two-phase flames, particle combustion regime, ignition, extinction, instability
1 Introduction
Understanding the mechanisms controlling flame propagation in suspensions of non-volatile particles is crucial to obtaining
efficient combustion of metallized propellants, slurry fuels, pulverized coal, and powdered metals as carbon-free chemical
energy carriers [1]. It is also necessary for the mitigation of catastrophic explosions in coal mines or in process industries that
involve handling metallic dusts and other combustible solid powders.
Like gas flames, flames in particulate suspensions at the laboratory scale are primarily driven by molecular heat diffusion
and have comparable burning velocities [1]. Nevertheless, they exhibit several significant differences in their structure and
behavior from homogeneous flames due to their multiphase nature.
The main distinctive feature of a flame in a solid suspension is the ability of particles to ignite–that is, to transition from
a combustion regime limited by reaction kinetics to a regime limited by diffusion of the oxidizing gas towards the particle
surface, or in the case of evaporating particles, towards the micro-flame enveloping each individual particle.
After ignition, the temperature of the particle or micro-flame can exceed the gas temperature by several hundred degrees,
often surpassing the adiabatic flame temperature for fuel-lean mixtures. The particle combustion rate in the diffusion combustion regime is a weak, non-Arrhenius, function of gas temperature. Unlike gas flames, the width of the flame reaction zone in
particle suspensions can span a large temperature range and can be comparable to, or even exceed, that of the preheat zone [2].
The existence of diffusion micro-flames within a global flame-front (in effect, flames within the flame), which are insensitive
to the bulk gas temperature, makes dust flames resistant to heat loss [1, 3, 4, 5] and also serves to maintain a constant burning
velocity with increasing fuel concentration in fuel-rich mixtures [6]. The ability of particles to ignite, together with low ignition
temperatures, may result in much wider flame propagation limits for particle suspensions than for gaseous fuels.
Despite the overall qualitative understanding of the crucial role of the particle combustion regime on burning velocity and
thermal structure, the theoretical description of flames in particulate suspensions has been limited to simple semi-empirical
1
models that postulate either purely diffusion or purely kinetic modes of particle combustion [6, 7]. The diffusive combustion models presume that particles within the suspension ignite and transit to the diffusion regime instantaneously when they
reach the ignition temperature of a single, isolated particle. The common assumption is that, after ignition, the particle within
the suspension will have a combustion time equal to that of an isolated single particle. Using this approach, particle ignition
temperature and combustion time are considered to be external parameters that are independent of the flame-propagation mechanism. As such, they often are taken from experiments with individual particles or calculated using theoretical models for single
particle combustion. These assumptions are useful for estimation, but are, in general, not justified and may lead to erroneous
predictions. In reality, the particle ignition temperature, particle reaction time, and the actual regime of particle combustion
are all characteristic values of the flame propagation problem directly linked to the burning velocity eigenvalue. Moreover, the
particle combustion regime may alternate throughout the flame. For example, the particle may start to react in a kinetically
limited regime, then transit to a diffusion-dominated combustion mode, before returning to a kinetics-dominated mode [3]. As
a result, a non-negligible fraction of the particle mass may be consumed during both diffusive and kinetic combustion, leading
to a complex dependence of the burning velocity on particle size and concentration, as demonstrated in this study.
In this paper, the thermal structure of a flame in a particulate suspension is investigated using a simplified, transparent
model that assumes that non-volatile solid fuel particles react via a single-step Arrhenius surface reaction with gaseous oxidizer
delivered to the particle surface by diffusion. Besides incorporating heterogeneous reaction kinetics, the model does not impose
the particle combustion mode or any other external combustion parameters. The flame propagation problem is solved numerically in a non-stationary formulation developed by Spalding [8]. This approach avoids the difficulties inherent in a steady-state
formulation and permits the investigation of flame stability, which has led to the discovery of a new type of oscillating flame in
heterogeneous mixtures.
2 Model Formulation
2.1 Combined Kinetic-Diffusive Reaction Rate
Following the quasi-stationary approach of Frank-Kamenetskii [4], the overall reaction rate per unit surface area of a particle
in an oxidizing gas, accounting for both kinetic and diffusion “resistances” can be written as
ω̇ = γ
κβ
C0
κ+β
(1)
where ω̇ is the particle mass consumption per unit surface area, γ is the stoichiometric coefficient, and C0 is the concentration
of oxidizer in the bulk gas far from the particle surface [9]. The kinetic term, κ, is the overall Arrhenius surface reaction
rate, κ = k0 exp (−Ea /Ru Ts ), where k0 is the pre-exponential factor, Ea is the activation energy, and Ru is the universal gas
constant. β is the mass transfer coefficient between a particle and the gas.
For simplicity, the contribution of Stefan flow to the heat and mass transfer between the particle and the surrounding gas
is assumed to be negligible. Stefan flow is small when the molecular weight of the reaction products is close to that of the
consumed oxidizer or when the initial oxidizer concentration is relatively low [4]. In the absence of Stefan flow, the mass
transfer coefficient takes the simple form, β = ShDi /2r. Here, Di denotes the oxidizer diffusivity at the particle-gas interface,
and r, the particle radius. For a spherical particle that is stationary relative to the gas, the exact solution to the steady-state
diffusion equation yields a Sherwood number, Sh, equal to 2 [4].
The reaction rate in Eq. (1) can then be written in terms of the normalized Damköhler number, Da∗ , which is related to the
traditional Damköhler number, Da = κ/β, by Da∗ = Da/(1 + Da) [9]. The resulting expression for the heat release rate per
unit particle surface area is:
Qr = q ω̇ = qγβDa∗ C0
(2)
where q is the heat of reaction.
This formulation naturally incorporates the two limiting kinetic and diffusive regimes. In a kinetically controlled regime
(β ≫ κ), Da∗ approaches zero. In the diffusion-limited regime (κ ≫ β), the reaction rate is primarily limited by diffusive
transport of oxidizer, and Da∗ approaches unity.
2.2 Ignition and Combustion of a Single Particle
The reaction of a suspension of particles within the flame is fundamentally different from the reaction of an isolated single particle [9]. However, mapping the reaction regimes for a single particle is crucial to understanding this difference and interpreting
the flame structure in a suspension. The interplay between kinetic and diffusion reaction rates of heterogeneous chemical reactions leading to the processes of ignition and extinction was first investigated by Frank-Kamenetskii [4] and then analyzed
2
further by Vulis [5]. Their analyses are adapted here in a simplified form to interpret the reaction behavior of a single particle
injected into a hot oxidizing gas.
The modified Semenov diagram shown in Fig. 1 plots the reaction heat release rate that accounts for both kinetic and
diffusive rates as a function of the particle temperature, Ts . At low temperatures, the reaction rate increases exponentially with
Ts due to Arrhenius kinetics. At high temperatures, Da∗ approaches unity and the reaction rate in Eq. (2) becomes practically
independent of temperature. If the heat loss from radiation is negligible, as is the case for moderate temperatures and small
particles, then the heat loss rate is proportional to the temperature difference between the particle and gas and is plotted in Fig.
1 as straight lines: Ql = h(Ts − Tg ). Here h is the heat transfer coefficient, given by h = Nuλi /2r, where λi is the thermal
conductivity of the gas at the gas-solid interface. The Nusselt Number, Nu, is equal to 2 for a spherical particle that is stationary
relative to the gas [4].
r
r
1
2
diffusion limit
Heat production/loss
D
E
Q
l
T
Q
r
K
I
i
r
crit
T
g0
r
0
T
s
Fig. 1: Modified Semenov diagram illustrating particle ignition (I) at some critical temperature (Tg0 for r2 ) and transition to near the diffusion limit (D).
Extinction (E) occurs as the particle shrinks and undergoes the reverse transition to the kinetic limit (K). Stable states are shown by (•) while unstable states are
denoted by (◦). The heat production curve is solid, and the heat loss curve is dashed. Inset: ignition temperature as a function of particle size showing critical
radius below which ignition is impossible [9].
At some critical gas temperature, known as the ignition temperature, the heat loss and the heat release curves become tangent
(point I on Fig. 1). The particle starts to accumulate heat and promptly transitions from a predominantly kinetic combustion
regime to a predominantly diffusive regime with temperatures close to the adiabatic temperature of the stoichiometric mixture
(point D). Once ignited, the burning particle inevitably extinguishes when its radius reduces to the point (r1 ) where the heat
loss and heat release curves become tangent again, albeit at a higher particle temperature (point E). After extinction, the particle
will transit to a kinetic burning regime (point K). For large particles, the extinction radius is small in comparison to the initial
particle size, and, thus, the residual particle mass at extinction is negligible since m ∼ r3 . For smaller particles, the mass of the
extinguished particle may not be negligible relative to the initial particle mass. If the initial particle size is reduced even further,
the slope of the heat loss curve becomes so steep that tangency of the heat release and heat loss curves becomes impossible.
Below this critical radius, the particle cannot ignite at any gas temperature (see inset in Fig. 1) [9].
The particle temperature history for three different initial sizes injected into a hot oxidizing gas is illustrated in Fig. 2. The
first case (r0 ≫ rcrit ) corresponds to a large particle that undergoes heating, then ignites and reacts to completion almost entirely
in the diffusive regime. The second case (r0 ∼ rcrit ) corresponds to an intermediate-sized particle that ignites but extinguishes
soon afterwards without achieving full-fledged diffusive combustion. After extinction, the particle continues to react in the
kinetic regime at a temperature close to the gas temperature. The third case (r0 < rcrit ) corresponds to a particle whose size is
below the critical value, and, as such, reacts entirely in the kinetic regime without ignition.
This analysis shows that the particle reaction regime is a strong function of the particle radius as well as the surrounding
temperature and oxidizer concentration. In a flame, a reacting particle with constantly changing radius is exposed to rising
temperature and falling oxidizer concentration fields that are functions of the burning velocity. Therefore, postulation of the
particle reaction regime and “ignition temperature” and “combustion time” parameters is, generally, unjustified in flame modelling. As demonstrated in this paper, only by solving the flame propagation problem as a whole can the particle combustion
regime and other fundamental flame parameters be determined.
2.3 Governing Equations and Numerical Method
The flame propagation problem is cast in a time-dependent formulation, based on the method developed by Spalding [8], and
solved numerically. This method has been used to analyze gas-phase flames [10]. Smoot et al. also used this method to
model flames in coal dust-air mixtures [11], but no attempt was made to draw general conclusions beyond the specifics of coal
combustion.
3
T
s
r
0
>> r
crit
r
T
g
r
r
0
0
~ r
< r
crit
crit
t
Fig. 2: Particle temperature and radius history of a particle injected into a hot oxdizing gas at temperature Tg burning in a predominantly diffusive combustion
regime (r0 ≫ rcrit ), an intermediate regime (r0 ∼ rcrit ), and a kinetic combustion regime (r0 < rcrit ) [9].
The model presented here employs several simplifying assumptions: (i) the solid fuel can be modelled using the continuum
approach described in [12, 13], (ii) the velocity slip between the particles and gas is negligible, (iii) the molecular weight and
heat capacity of the products and oxidizer are similar, and (iv) the solid particles do not undergo any phase transitions [14]. All
of these are second-order effects.
The governing equations for the flame are simplified by introducing a density-weighted coordinate, x, which is related to
R x′
the physical coordinate, x′ , by x = 0 [ρg (x̄′ )/ρg0 ]dx̄′ [8]. In this coordinate system, the continuity equation is automatically
satisfied, and advection effects induced by thermal expansion of the gas are absent. Once found, the solution in x-space can be
transformed into x′ -space, as described by Margolis [15].
For simplicity, the temperature dependences of transport and kinetic properties are chosen to eliminate thermal expansion effects. With these assumptions, the governing equations can be reduced to the following expressions for the gas-phase
temperature (Eq. (3)), particle temperature (Eq. (4)), particle mass, ms (Eq. (5)), and normalized oxidizer mass fraction,
Y = Yox /Yox,0 (Eq. (6)). The combined kinetic-diffusion reaction rate, ω̇, is given by Eq. (1).
ρg0 cg
cs
∂ 2 Tg
∂Tg
+ N0 Ah(Ts − Tg )
= λ0
∂t
∂x2
∂
(ms Ts ) = Aq ω̇Y − Ah(Ts − Tg )
∂t
∂ms
= −Aω̇Y
∂t
∂2Y
N0 ω̇A
∂Y
Y
= D0 2 −
∂t
∂x
γρg0
(3)
(4)
(5)
(6)
Here, A is the instantaneous particle surface area, cs and cg are the solid and gas-phase specific heats, D0 and λ0 are the
bulk mass diffusivity and thermal conductivity of the gas, respectively, ρg0 is the initial gas density, and N0 is the initial particle
number density. The numerical values for the thermodynamic, transport and reaction properties used are given in Table 1.
These values are similar to those for an aluminum particle in air, but the model is not meant to describe aluminum combustion
in particular. The reaction parameters are chosen such that the switch from kinetics- to diffusion-dominated particle combustion
can be observed in the micrometer particle size range.
The equations were numerically integrated in MATLAB by the method of lines using a second-order finite difference
formulation and a multi-step, variable order, implicit time integration scheme [16]. The integration was performed over an
effectively semi-infinite domain with an adiabatic and impermeable hot wall condition [15]. Because long-term solutions were
of interest, the exact initial conditions used in this study turned out to be immaterial [17].
3 Results and Discussion
3.1 Flame Thermal Structure
After an initial transient period, a steady-state flame, propagating at a constant speed, typically develops. As discussed below,
for some values of φ and r0 , no steady-state solution is observed. Steady-state flame profiles of gas and particle temperatures,
oxidizer mass fraction and particle radius are shown for the fuel-lean case (φ = 0.5) in Fig. 3 for three different initial particle
4
Table 1: Numerical values of gas and solid-fuel parameters.
λ0 , λi
cg
D0 , Di
ρg0
γ
C0
0.02 W1 m−1 K−1
1010 J1 kg−1 K−1
2.0 × 10−5 m2 s−1
1.17 kg1 m−3
1.12
0.27 kg1 m−3
Ea /Ru
cs
k0
ρs
q
15 000 K
979 J1 kg−1 K−1
70 m1 s−1
2700 kg1 m−3
3.10×107 J1 kg−1
sizes corresponding to overall kinetic (r0 = 1 µm), diffusive (r0 = 10 µm), and intermediate (r0 = 4 µm) particle combustion
regimes.
0.6
r
0.5
= 1
m
a
r = 4
0.5
m
0
0.4
0.3
0.3
0.2
0.2
0.1
0.1
0.0
0.0
0.0
1.0
1.0
1.0
0.5
0.5
0.3
T
T
0.2
0
m
s
g
r
Y
0.5
r/ r
0
and
Y
0.1
r = 10
0.5
0.4
0.4
T/T
0.6
0.6
0
Da* and
q
0.0
1.0
0.0
0.0
1.0
1.0
0.5
0.5
Da*
q
0.5
0.0
0.0
0.0
-4
-2
0
x (mm)
2
4
-4
-2
0
x (mm)
2
4
-4
-2
0
2
4
x (mm)
Fig. 3: Flame structure profiles for three different particle sizes in a fuel-lean flame (φ = 0.5). Top: gas and particle temperature non-dimensionalized by
Ta = Ea /Ru ; Middle: normalized oxidizer mass fraction, Y , and particle radius, r/r0 ; Bottom: reaction heat release rate, q̇, and Da∗ . The arrow indicates
direction of propagation.
For small particle sizes, the flame structure is similar to that of a homogeneous gaseous flame. The gas and particle
temperatures remain close to each other throughout the flame, and the reaction zone is thin and located near the point where the
temperature reaches the adiabatic flame temperature of the mixture.
In contrast, the flame structure for large particles is significantly different. After ignition, the particles attain a temperature
close to the adiabatic flame temperature of the stoichiometric mixture, which greatly exceeds that of the fuel-lean mixture.
In addition, the reaction zone is wide, encompassing a large gas temperature range. There is also a slight lag of the particle
temperature profile behind the gas temperature in the flame preheat zone.
In the case of intermediate particle sizes, particles ignite, but then extinguish before achieving full-fledged diffusion-limited
combustion. Thus, the particle temperature surpasses the adiabatic flame temperature of the mixture, but does not reach that
of the stoichiometric mixture. In contrast to the reaction of a single particle of the same size, the kinetic reaction time after
the transition back to the kinetic combustion mode is relatively short due to the high flame temperatures and, thus, does not
significantly increase the overall combustion time.
Flames in fuel-rich suspensions, as shown in Fig. 4, have a different structure to that of fuel-lean mixtures, since the
limiting reactant changes from the solid fuel, which has a negligible mass diffusivity, to the oxidizing gas, which readily
diffuses across the flame. The primary difference is that the temperature of the particles after ignition does not exceed the
adiabatic flame temperature of the mixture. The reason for this is that the temperature separation between particles and gas in
the diffusive combustion regime is proportional to the concentration of the oxidizer in a particular location, which falls to zero
as the temperature rises. Due to complete oxidizer consumption at the end of the combustion zone and, therefore, a negligibly
small reaction rate, the value of Da∗ remains constant after the oxidizer is depleted.
3.2 Burning Velocity
The burning velocity is calculated based on the flame-front displacement as a function of time. The burning velocity is plotted
in Fig. 5 as a function of initial particle size for different equivalence ratios. The general trend for stoichiometric, rich, and
somewhat lean mixtures is that the burning velocity increases monotonically with decreasing particle size. This is simply due to
the increasing specific reaction surface area with decreasing particle size. In contrast, this dependence becomes non-monotonic
for very lean mixtures. For equivalence ratios below φ = 0.4, there is also a region where the burning velocity curves split into
two separate branches between which no steady solution exists.
5
0.6
r
0.5
m
= 1
a
r = 4
0.5
m
0
0.4
0.3
0.3
0.2
0.2
0.1
0.1
0.0
0.0
0.0
1.0
1.0
1.0
0.5
0.5
0.3
T
T
0.2
g
0
m
r
Y
0.5
r/ r
0
and
Y
0.1
s
r = 10
0.5
0.4
0.4
T/T
0.6
0.6
0
Da* and
q
0.0
1.0
Da*
q
0.5
0.0
0.0
0.0
1.0
1.0
0.5
0.5
0.0
0.0
-4
-2
0
2
4
-4
-2
0
x (mm)
2
-4
4
-2
0
2
4
x (mm)
x (mm)
Fig. 4: Flame structure profiles for three different particle sizes in a fuel-rich mixture (φ = 1.5). The same variables are shown as in Fig. 3.
10
8
4
u
b
(cm/s)
6
2
0
2
4
6
r
0
8
10
( m)
Fig. 5: Burning velocity, ub , as a function of particle size for different equivalence ratios. Hatched area corresponds to regions where unstable flame propagation
is observed.
The nontrivial dependence of the burning velocity on particle size in lean mixtures is the result of the intricate interplay
between the kinetic and diffusive particle combustion modes. Because of the relatively low flame temperatures in lean mixtures,
the combustion time of a larger particle reacting diffusively can be shorter than that of a smaller particle reacting kinetically.
This behavior is counter-intuitive, but is a direct consequence of the ability of large particles to ignite and burn with temperatures
much higher than the adiabatic flame temperature of the lean mixture. The appearance of “flames within the flame”, results in
higher reaction rates for suspensions of larger particles that offset, and potentially even overcome, the reduced specific reaction
surface area as compared to finer suspensions. This behavior can be seen first as a plateau in the dependence of burning velocity
on particle radius (Fig. 5, φ = 0.4) and then a decline as particle size decreases (Fig. 5, φ < 0.4). It can also be seen in Fig.
6 as intersections of the burning velocity versus equivalence ratio curves for different particle sizes. As the particle size is
increased, eventually the effect of decreasing specific reaction surface area again dominates, and the burning velocity decreases
with increasing particle size. These trends have been observed experimentally in the combustion of coal suspensions [11] and
liquid sprays [18], although no clear physical explanation was given at the time.
3.3 Kinetic-Diffusive Instabilities
For equivalence ratios below φ = 0.4, there is a region where no steady-state flame propagation occurs (see hatched area of
Fig. 5). Instead, these flames oscillate both in terms of burning velocity and flame structure. An example of instantaneous
profiles of particle and gas temperature, and Da∗ at different instances over a single oscillation period are shown in Fig. 7 for
an equivalence ratio of φ = 0.4 and initial particle size of r0 = 4 µm.
The gas temperature is observed to oscillate about the adiabatic flame temperature, while the particle temperature oscillates
with a much greater amplitude than that of the gas. At t3 , the combustion regime is more diffusive, as indicated by the large
temperature separation and Da∗ , while at t1 and t5 , the regime is more kinetic, as indicated by the low values of temperature
separation and Da∗ .
At first glance, these oscillations resemble the well-known thermo-diffusive flame instability often observed in homogeneous gas [19] and quasi-homogeneous solid mixtures [20], and also recently observed in aluminum dust flames [21]. The
6
4
r = 1
(cm/s)
0
(cm/s)
0
r = 3
m
2
m
r = 10
m
0
1
0.30
0.32
0.34
0.36
0.38
m
r = 6
0
u
r = 1
m
0
b
10
3
0.40
u
b
5
r = 10
0
m
0
2
4
6
8
10
Fig. 6: Burning velocity as a function of equivalence ratio for different initial particle sizes. The expanded view shows the intersections of burning velocities
for lean equivalence ratios for various particle sizes.
= 0.325 r = 4
0.4
m
0
T
s
T
g
T/T
a
0.3
T
0.2
ad
t
3
t
0.1
2
t
1
= 0.70
t
= 0.63
4
= 0.80
t
= 0.45
5
= 0.97
0.0
t
3
Da
*
0.4
0.2
t
t
4
t
2
1
t
5
4
8
12
16
20
x (mm)
Fig. 7: Instantenous profiles of gas and particle temperatures and Da∗ at different instances in time over a single oscillation period of an oscillating flame. The
units of t are seconds. The arrow indicates the direction of propagation.
physical mechanisms underlying the thermo-diffusive instability and the heterogeneous pulsating flame shown here are, however, fundamentally different. In homogeneous mixtures, pulsations are induced by changes in the gas molecular transport
properties (Lewis number), activation energy, and overall heat production, as required by the thermo-diffusive instability criteria formulated by Matkowsky and coworkers [22, 23]. In the particle flame, the fuel has zero diffusivity, and, for a fixed
equivalence ratio, the gas-phase transport, kinetic parameters, and overall heat production are fixed. In particular, they do not
change with particle size. Instead, the cause of the observed instability is evidently rooted in the processes of particle ignition
and extinction that are, indeed, very sensitive to particle size. Thus, it may be called a kinetic-diffusive instability, reflecting
the fact that the particle combustion regime in the pulsating flame oscillates between diffusive and kinetic combustion modes.
This interesting new phenomenon will be studied in detail in subsequent work.
3.4 Comparison to Semi-Empirical Flame Models
Most existing models for flames in suspensions are semi-empirical and postulate either a purely kinetic or purely diffusive
regime of particle combustion. Plots of the normalized burning velocity versus particle size predicted by such models are
compared with the result of the current numerical simulations in Fig. 8 for the lean mixture with an equivalence ratio of
φ = 0.3.
For the case of the purely kinetic regime, the burning velocity is calculated using an analytical expression formulated by
Zeldovich and Frank-Kamenetskii, modified for a combustible mixture with a single-step, first order heterogeneous surface
reaction [17]. The result is normalized to match that of numerical simulations for the smallest particles, which burn in the
kinetic regime.
The burning velocity in the case of the purely diffusive regime is calculated using an expression given in [2, 6] that requires
two externally defined combustion parameters: the ignition temperature and combustion time of the single particle. Ignition
temperatures at different particle sizes are calculated using the Semenov ignition criterion as shown in Fig. 1 and have a
dependence with particle size similar to that shown in the inset of Fig. 1. The particle combustion time is calculated using a
simple diffusive reaction model that predicts a tc ∼ r02 dependence [24]. The result of the purely diffusive model is normalized
7
2.5
current model
(cm/s)
diffusive regime
1.5
1.0
u
b
kinetic regime
0.30
2.0
0.5
0.0
5
10
15
20
25
r ( m)
0
Fig. 8: Comparison of computed burning velocity to semi-empirical asymptotic Zeldovich and Frank-Kamenetskii (kinetic) and modified Mallard-Le Chatelier
(diffusive) flame models. The semi-empirical models are normalized to the numerical model at r0 = 1 µm for the kinetic model and r0 = 25 µm for the
diffusive model.
to match that of the numerical simulations for the largest particles, which burn in the diffusive regime.
As can be seen in Fig. 8, these semi-empirical models only qualitatively predict the asymptotic behavior of burning velocity
on particle size where the combustion regimes are either purely kinetic or diffusive. Over a wide particle size range, where
diffusion and kinetic rates are comparable, these models cannot capture the physics of the flame propagation and yield erroneous
results.
4 Conclusion
The comprehensive analysis of flame propagation in non-volatile particle suspensions presented in this paper has a historical
parallel. It can be compared with the replacement of the semi-empirical Mallard-Le Chatelier flame model, based on the notion
of ignition temperature and reaction time as external parameters, to flame models developed in the 20th century based on the
Arrhenius reaction law. Here, it is similarly shown that the particle ignition temperature, particle reaction time, and regime
of particle combustion are, in general, characteristic values of the flame problem that must be solved for together with the
burning velocity. Moreover, this analysis demonstrates that the thermal structure of flames in suspensions is defined primarily
by the interplay between kinetic and diffusive reaction modes that also yields non-trivial dependence of the burning velocity on
particle size. This interplay leads to the emergence of a new type of pulsating instability that is related to the process of ignition
and extinction of particles within the flame.
Acknowledgments
Support for this work was provided by the Defense Threat Reduction Agency under contract HDTRA1-11- 1-0014, the NSERC
CREATE on Clean Combustion Engines, and a McGill Engineering Doctoral Award. Additional funding was provided by the
Panda Faculty Scholarship in Sustainable Engineering & Design and a William Dawson Scholarship.
References
[1] J. Bergthorson, S. Goroshin, M. Soo, P. Julien, J. Palecka, D. Frost, D. Jarvis, Appl. Energy 160 (2015) 368–382.
[2] S. Goroshin, M. Bidabadi, J. H. S. Lee, Combust. Flame 105 (1) (1996) 147–160.
[3] F.-D. Tang, S. Goroshin, A. J. Higgins, Proc. Combust. Inst. 33 (2) (2011) 1975–1982.
[4] D. A. Frank-Kamenetskii, Diffusion and Heat Transfer in Chemical Kinetics, Plenum Press, New York, 1969.
[5] L. A. Vulis, Thermal Regimes of Combustion, McGraw-Hill, New York, 1961.
[6] S. Goroshin, I. Fomenko, J. H. S. Lee, Proc. Combust. Inst. 26 (2) (1996) 1961–1967.
[7] K. Seshadri, A. Berlad, V. Tangirala, Combust. Flame 89 (3-4) (1992) 333–342.
[8] D. B. Spalding, Philos. Trans. R. Soc. London, Ser. A 249 (957) (1956) 1–25.
[9] M. Soo, S. Goroshin, J. M. Bergthorson, D. L. Frost, Propellants Explos. Pyrotech. 40 (4) (2015) 604–612.
[10] Y. B. Zeldovich, G. I. Barenblatt, Combust. Flame 3 (1959) 61–74.
8
[11] L. Douglas Smoot, M. Duane Horton, Prog. Energy Combust. Sci. 3 (4) (1977) 235–258.
[12] S. Goroshin, F.-D. Tang, A. J. Higgins, Phys. Rev. E 84 (2011) 027301.
[13] S. Goroshin, J. H. S. Lee, Y. Shoshin, Proc. Combust. Inst. 27 (1) (1998) 743–749.
[14] S. V. Goroshin, Y. L. Shoshin, A. N. Zolotko, Akademiia Nauk SSSR, Doklady 321 (3) (1991) 548–551.
[15] S. B. Margolis, J. Comput. Phys. 27 (3) (1978) 410–427.
[16] W. E. Schiesser, G. W. Griffiths, A Compendium of Partial Differential Equation Models: Method of Lines Analysis with
Matlab, Cambridge University Press, 2009.
[17] Y. B. Zeldovich, G. I. Barenblatt, V. B. Librovich, G. M. Makhviladze, Mathematical Theory of Combustion and Explosions, Consultants Bureau, New York, 1985.
[18] K.-K. Chan, C.-S. Jou, Fuel 67 (9) (1988) 1223–1227.
[19] E. Christiansen, C. Sung, C. Law, Proc. Combust. Inst. 27 (1) (1998) 555–562.
[20] H. Yi, J. Moore, J. Mater. Sci. 25 (2) (1990) 1159–1168.
[21] P. Julien, J. Vickery, S. Goroshin, D. L. Frost, J. M. Bergthorson, Combust. Flame 162 (11) (2015) 4241–4253.
[22] B. J. Matkowsky, G. I. Sivashinsky, SIAM J. Appl. Math. 35 (3) (1978) 465–478.
[23] B. J. Matkowsky, D. O. Olagunju, SIAM J. Appl. Math. 39 (2) (1980) 290–300.
[24] F.-D. Tang, S. Goroshin, A. Higgins, J. Lee, Proc. Combust. Inst. 32 (2) (2009) 1905–1912.
9