3.1. Gaseous Combustion
We start with
Figure 2 presenting the time evolutions of the flame tip position,
,
Figure 2a, and its velocity,
,
Figure 2b, in a 2D geometry, for stoichiometric methane (CH
4)-air burning, as well as that of propane (C
3H
8)-air, for comparison, with various blockages
employed. The horizontal dotted lines in
Figure 2b show the speeds of sound for the methane-air,
, and propane-air,
, mixtures (
Table 1). The case of no obstacles,
, reproduces, completely, the situation of “finger + DL” flame acceleration [
3]. It is noted that this acceleration is limited in time such that the flame would start decelerating when its skirt contacts a sidewall at
and
for methane-air and propane-air burning, respectively. It is also noted that by these times, the propane-air flame overcomes the sound threshold
, whereas the methane-air flame stops at
. In contrast, in an obstructed channel,
, acceleration is unlimited in time until the flame approaches the speed of sound and can eventually trigger a detonation. We should recall, at this point, that approaching the near-sonic values by the flame front will eventually break the incompressible approach, adopted in Equation (6), and the entire present formulation. Indeed, to describe the DDT stage accurately, we have to incorporate the impacts of gas compressibility into the present analysis. This can be done by considering the formulation of
Section 2 as the zeroth-order approach in
, and then extending it to account for the finite
according to the methodology employed earlier for unobstructed [
27] and obstructed [
28] geometries. However, such a rigorous extension of the present formulation to account for the compressibility effects requires a separate work and is presented elsewhere [
29]. It is also noticed that the analytical incompressible formulation of
Section 2 does not involve pressure as a parameter (except for the fact that pressure comes indirectly through the thermal-chemical parameters such as
or
, which are taken for ambient, atmospheric pressure indeed in the present section). However, if the present formulation is extended to account for a finite
, as discussed above, then pressure, its variations and gradient, are directly involved in a revised formulation, with the ambient, atmospheric pressure imposed as the initial conditions and a boundary condition in the open side of the passage.
Figure 3 is a counterpart of
Figure 2 for the cylindrical geometry. It is seen that the flame accelerates faster in this case. Opposite to the situation of
and methane-air burning in a 2D geometry, here, at the same conditions, the flame would overcome the sound speed at
, i.e., slightly prior the deceleration stage starting at
. To demonstrate this more clearly, in
Figure 4 we have compared the results obtained in the 2D and cylindrical-axisymmetric geometries. Overall,
Figure 2,
Figure 3 and
Figure 4 show that the obstacles influence a coal mine fire scenario significantly, making acceleration potentially unlimited in contrast to the case of no obstacles, which was considered in [
3]. In addition, flame acceleration in
Figure 2,
Figure 3 and
Figure 4 exceeds that of the original theories [
14,
15] by orders of magnitude, thereby certifying that the DL instability facilitates obstacle-based acceleration.
It should be mentioned that everywhere except for
Figure 5 we employed the exponent
of Equation (1) to be the same as that in [
3], i.e.,
. This is in order to compare the present work with the “unobstructed burning accident” theories [
3] (as well as with the original Bychkov theories [
14,
15], which did not consider the DL instability, thereby having
by default). It is noted, in this respect, that an appropriate choice of the factor
is not yet finalized. Indeed, while various experimental studies such as Gostintsev et al. [
18], Bradley et al. [
19], Molkov et al. [
20], Kim et al. [
21] have reported
, the Princeton experiments in a dual chamber at elevated pressures suggested
; see Jomaas et al. [
22] and numerous references therein. The theoretical work [
30] has provided an attempt to explain such a discrepancy between various experiments by the flame-acoustic coupling, which modifies
. Moreover, the quantity of
can also be potentially modified by turbulence and other combustion instability modes such as Rayleigh–Taylor instability (assuming that the general trend of Equation (1) remains).
We have considered
in the range
and scrutinized its impact on the present formulation. Specifically,
Figure 5 compares the time evolutions of the flame tip positions (a, c) and velocities (b, d) for stoichiometric methane-air burning considering
, and both 2D (
Figure 5a,b) and cylindrical (
Figure 5c,d), geometries are studied. It is seen that the variations of
impact the flame position, velocity and acceleration substantially, although the effect is quantitative but not qualitative. As expected, flame acceleration in the case of
proceeds noticeably faster than for
, whereas a flame with
accelerates noticeably slower as compared with that with
(of course, provided that other combustion characteristics are kept the same). Nevertheless, without the final answer to the question about an appropriate choice for
, in the rest of this work we keep using its median value,
, the same as that in the unobstructed theories [
3]. This allow us to make a comparison with [
3], thereby separating the impact of obstacles.
Next, we extend the stoichiometric (
) gaseous methane-air combustion considered in
Figure 2,
Figure 3 and
Figure 4 to the equivalence ratios in the range
, see
Figure 6,
Figure 7,
Figure 8 and
Figure 9. In particular,
Figure 6 presents the 2D configuration with various
for
. It is seen that a slightly fuel-lean flame with
accelerates much slower than the
flames, especially in the cases of
and
. This is because of a much lower
(and thereby higher
and lower
) inherent to such a slightly lean condition. However, it is recalled that flame acceleration in an obstructed passage is unlimited in time, and therefore it can eventually trigger a detonation in the case of sufficiently long passage and time. In particular, in this geometry, the approximate run-up times until the detonation initiation for the
flame are evaluated as
for
and
for
. In the case of no obstacles,
, no detonation is predicted for a
flame in a 2D passage. Overall, among all equivalence ratios considered, the fastest flame acceleration is observed for a slightly rich flame of
.
The flame will propagate slower if we deviate further from stoichiometry, as depicted in
Figure 7, showing the evolutions of the position and velocity of the
and
methane-air flames (to avoid messy, we split the plots for various equivalence ratios between
Figure 6 and
Figure 7). It is also seen that the duration of acceleration of the highly lean/rich flames in
Figure 7 exceeds that of the stoichiometric or slightly lean/rich flames in
Figure 6. In particular, without obstructions, the
flame stops accelerating after advancing
and attaining the maximal velocity of
, before the flame skirt contacts the wall. Obviously, this acceleration scenario does not end thereafter, if the obstacles are added to the passage wall. In contrast, the flame front will keep accelerating until its speed reaches the speed of sound and, eventually, the detonation is triggered. It is noted that the obstacles facilitate flame acceleration as compared with that in an unobstructed passage, and the impact of obstacles seems to be more important for the
flames than for the
flames: the lines for
with various
go more widely than those for
. It is also seen in
Figure 7 that the rich flames,
, accelerate faster than the lean flames,
. With respect to the latter, the situation of highly non-stoichiometric combustion in
Figure 7 qualitatively resembles slightly non-stoichiometric burning in
Figure 6.
Figure 8 and
Figure 9 are the cylindrical-axisymmetric counterparts of
Figure 6 and
Figure 7, respectively. Here, most of the 2D results discussed above remain qualitatively the same, although, quantitatively, the flames accelerate faster in the cylindrical passages. Here, the approximate run-up times until the detonation initiation for the
flames appear
for
and
for
. A key difference between the two geometries is that whereas no detonation is predicted for the
flame in a 2D case without obstacles,
, in the cylindrical configuration with
, the
methane-air flame was able to reach the speed of sound,
, at
, thereby making a detonation possible. Again, further away from stoichiometry, see
Figure 9, the
and
flames accelerate slower and the acceleration time lasts longer as compared to the
flames in
Figure 8. For example, in the case of
, the
flame stops accelerating when advancing
and reaches the maximal velocity of
s at the instant when the flame skirt contacts the side wall of the cylindrical passage. Similar to a 2D geometry, the highly rich flames (
) accelerate faster than the highly lean flames (
) in the cylindrical case, and the impact of obstacles seems to be more important for the
flames (the lines corresponding to various blockage ratios go more widely).
Using
Table 1, we next analyze the run-up times
, Equations (18) and (34), and the respective run-up distances
, Equations (19) and (35), for the methane (CH
4)-air and propane (C
3H
8)-air flames of various equivalence ratios. Specifically,
Figure 10 presents
versus
for various blockage ratios, including the case of no obstacles,
, in the 2D (
Figure 10a) and cylindrical (
Figure 10b) geometries. Overall,
Figure 10 agrees with our analysis above in that the shortest run-up distances are observed for a slightly fuel-rich methane-air flame of
. In the 2D case,
Figure 10a, we have
,
, and
for
, respectively. The case of
in a 2D geometry is not relevant because a flame skirt contacts a sidewall and stops accelerating before the DDT event for all
considered, which is in line with the findings of [
3]. For the lean or rich methane-air flames, the run-up distances are much higher, namely, up to
for
and up to
for
(still in the 2D geometry). In the cylindrical-axisymmetric configuration, for the fastest methane-air flames with
we found the run-up distances as small as
,
,
, and
for
, respectively. For lean or rich methane-air burning, the run-up distances are much higher in the cylindrical geometry, i.e., up to
for
and up to
for
. It is noted that, unlike the 2D configuration, in the cylindrical-axisymmetric case, the detonation is predicted for methane-air burning even in the case of
, when the equivalence ratio is in the range
, which generally agrees with [
3]. Overall, for the same geometry,
and
, the run-up distances are dramatically shorter for the C
3H
8-air flames as compared with the CH
4-air flames; this would overcome the sound threshold for additional equivalence ratios in the case of
, being within the ranges
and
for the 2D and cylindrical-axisymmetric geometries, respectively.
3.3. Extension to Gaseous-Dusty Environment
Starting with homogeneously gaseous combustion in the previous subsection, we next extend our analysis to a gaseous-dusty environment by using a modified version of the Seshadri formulation [
10], which expresses the laminar flame velocity as a function of local thermal-chemical properties of the gas and dust particles (inert, such as sand; combustible, i.e., coal; and combined) in the form
[
3]:
where
is the modified equivalence ratio of the gaseous-dusty-air mixture in the presence of combustible dust particles:
and
are the respective molar masses;
,
and
are the original masses per unit volume for a given equivalence ratio;
is the specific heat of the whole mixture, containing the components for the gas,
, and dust particles,
;
is the density of a single dust particle, while
is that for the gaseous-dusty fuel-air mixture, with the density of the gas
and the concentration of the dust particles
;
is the number of particles per unit volume, with
being the volume of a single particle of radius
;
is the universal gas constant;
is the reactants temperature and
is the adiabatic flame temperature based on the purely methane-air equivalence ratio. Similar to [
3], here,
is calculated as a fifth-order polynomial function of the equivalence ratio
[
32]:
valid in the range
[
32]. Knowing
from Equation (37) and calculating a new flame temperature for the dusty-gaseous mixture,
, one can find
, Equation (36), to be employed in the formulation in the previous section.
We next employ the methodology of [
3,
11]. Unlike a combustible dust, an inert particle acts as a heat sink because it absorbs some heat from the flame and reduces the flame temperature. For
, methane-air combustion, the global chemical reaction is described by the equation:
such that the heat release in the process of burning of
moles of
and
moles of air is given by [
3,
11]:
where
is the number of moles of the burning products, which depends on the equivalence ratio
. Assuming that the entire heat released from the reaction is used to raise the temperature of the mixture, Equation (40) can be modified into an expression for the volumetric heat release from
-air combustion of a given
as follows [
3]:
Next, it is assumed that a flame with particles releases the same amount of heat while it is also influenced by the temperature rise of particles [
11]. Then Equation (41) can be extended as [
3]:
where
[
11] is the heat of gasification per unit volume. From the last equation, the secondly revised flame temperature,
, is calculated as:
Finally, the effect of the combination of the inert and combustible dust particles are accounted for by averaging the adiabatic flame temperature over those values associated with both effects separately,
. Similar to a combustible dust case,
and
are used in Equation (36) to find a new laminar flame speed
for the inert and combined dust particle incorporations, respectively. As of now, all the particles are assumed to be distributed uniformly inside a coal mining passage; non-uniform dust distributions can be considered elsewhere, for instance, following the method of [
33].
It is even a more interesting question with respect to what happens with the DL cutoff wavelength when we go from the gaseous to the gaseous-dusty environment. To be self-consistent with the Seshadri formulation (which actually imitates a multi-phase system by an “effective” fluid with modified properties, adjusted due to the presence of solid particles), here, is considered to be the quantity devoted to such an effective fluid. Consequently, is given by the same formula, Equation (3), although the variables in this formula (, ) are adjusted due to the presence of dust particles.
We can also provide another justification for this approach, namely, because , if the flame thickness does not change much when the dust particles are laden, then the value , as well as the onset and emergence of the DL instability would also not change much. Nevertheless, we recognize that a rigorous analysis of the DL instability of gaseous-dusty flames would require a separate study, with extra features to be incorporated. In particular, the transport processes such as heat transfer could be modified in an intriguing way, in particular, due to radiation.
Figure 13 depicts the situation of gaseous-dusty combustion, with the dust of concentration
and of dust particles radius
in the 2D (a, b) and cylindrical-axisymmetric (c, d) geometries. We employed the lean (
) methane-air fuel mixture and various blockage ratios
, including the case of no obstacles,
. It is seen that combined (combustible + inert) and inert dust moderate flame acceleration, whereas combustible particles slightly facilitate flame propagation. The effect of obstacles is also noticeable (as compare with an unobstructed passage).
Next, we increase the dust concentration. Namely, in
Figure 14, it is promoted to
keeping the same particle radius,
, and other characteristics and geometry, as in
Figure 13. It is observed that combustible dust promotes flame acceleration, whereas inert dust and its combination with combustible dust moderate the acceleration process for the particles of radius
. The impact of the blockage ratio is noticeable, especially in the cylindrical geometry. The flame velocity in the case of inert particles and
is equivalent to the event of combustible dust but with no obstacles,
, for
; and thereafter the flame accelerates faster in the presence of inert particles (
Figure 14d). A relatively high concentration was also considered. Specifically,
Figure 15 presents the case of
, with the particle size
, the same as in
Figure 13 and
Figure 14. This investigation reveals that if the dust particle size is kept constant while increasing the number of particles (i.e., the concentration), the effect of particles becomes more important. This is observed in
Figure 13,
Figure 14 and
Figure 15 that the flames in gaseous-dusty environments departs from the case of no dust particles more widely as the concentration increases. Similar to the cases of
and
, the combustible particles promote flame acceleration, whereas the combined (combustible + inert) and inert particles also suppress it in the situation of
.
Figure 16a–d is the counterpart of
Figure 14a–d, respectively, for a smaller particle radius,
. It is observed, here, that smaller particles have a stronger impact on flame propagation. In the 2D geometry, while the flame velocities did not exceed
for the particles of size
, in the case of
, the sound threshold of
for
methane-air burning was reached in the combustible coal gaseous-dusty environment with
at the time instant
(
Figure 16b). In the cylindrical-axisymmetric geometry, the coal particles deviate noticeably larger from the case of no particles (
Figure 16c,d). In fact, a particle type appears to be the most influential factor for flame acceleration in the obstructed passages. Namely, in any case of combustible dust present in any obstructed passages considered,
, we obtained faster acceleration than in both respective cases of other dusts and no dust (
Figure 16d). In contrast to the
case in
Figure 14, the combined combustible-inert particles promoted flame acceleration in the case of
,
Figure 16. These findings show that the impact of heat release on flame acceleration in a coal mining passage is significant, i.e., it facilitates the fire process, and, furthermore, it dominates over the impact of a heat sink when the particles are smaller. As for the inert particles, similar to the case of
in
Figure 14, they also suppress flame acceleration for
. The aforementioned effects of all particles, i.e., combustible, inert, and combined, and particle sizes are found to grow with the blockage ratio
.