Next Article in Journal
Interannual to Decadal Variations of Submesoscale Motions around the North Pacific Subtropical Countercurrent
Next Article in Special Issue
Scale-Resolving Simulation of a Propane-Fuelled Industrial Gas Turbine Combustor Using Finite-Rate Tabulated Chemistry
Previous Article in Journal
Uncertainty Quantification Methodologies Applied to the Rotor Tip Clearance Effect in a Twin Scroll Radial Turbine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Gaseous and Gaseous-Dusty, Premixed Flame Propagation in Obstructed Passages with Tightly Placed Obstacles

by
Furkan Kodakoglu
1,
Sinan Demir
1,
Damir Valiev
2,3 and
V’yacheslav Akkerman
1,*
1
Center for Innovation in Gas Research and Utilization (CIGRU), Department of Mechanical and Aerospace Engineering, West Virginia University, 1306 Evansdale Drive, Morgantown, WV 26506, USA
2
Center for Combustion Energy, Key Laboratory for Thermal Science and Power Engineering of the Ministry of Education of China, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, China
3
Department of Applied Physics and Electronics, Umeå University, 901 87 Umeå, Sweden
*
Author to whom correspondence should be addressed.
Fluids 2020, 5(3), 115; https://doi.org/10.3390/fluids5030115
Submission received: 9 June 2020 / Revised: 9 July 2020 / Accepted: 11 July 2020 / Published: 17 July 2020

Abstract

:
A recent predictive scenario of premixed flame propagation in unobstructed passages is extended to account for obstructions that can be encountered in facilities dealing with explosive materials such as in coalmines. Specifically, the theory of globally-spherical, self-accelerating premixed expanding flames and that of flame acceleration in obstructed conduits are combined to form a new analytical formulation. The coalmining configuration is imitated by two-dimensional and cylindrical passages of high aspect ratio, with a comb-shaped array of tightly placed obstacles attached to the walls. It is assumed that the spacing between the obstacles is much less or, at least, does not exceed the obstacle height. The passage has one extreme open end such that a flame is ignited at a closed end and propagates to an exit. The key stages of the flame evolution such as the velocity of the flame front and the run-up distance are scrutinized for variety of the flame and mining parameters. Starting with gaseous methane-air and propane-air flames, the analysis is subsequently extended to gaseous-dusty environments. Specifically, the coal (combustible, i.e., facilitating the fire) and inert (such as sand, moderating the process) dust and their combinations are considered, and the impact of the size and concentration of the dust particles on flame acceleration is quantified. Overall, the influence of both the obstacles and the combustion instability on the fire scenario is substantial, and it gets stronger with the blockage ratio.

1. Introduction

Historically, coal mining exhibits one of the highest fatality and injury rates among the industries dealing with flammable gases and combustible dust, with occupational hazards occurring on a regular basis [1]. For instance, not long ago an unfortunate mining accident in Soma, Turkey took hundreds of lives [2]. Among the causes of disasters, a spontaneous accidental explosion can lead to a catastrophic consequence. This happens, in particular, because methane is often present in the coal seams and can be accumulated inside a coal mining passage during the coal extraction process. If sufficient ventilation is not provided, a risk of violent methane/air/coal-dust explosion is viable. Consequently, in order to foresee potential risks and to install any safety measures against a fire accident, an understanding of the processes in a burning accident is essential. A predictive scenario of a coal mining fire [3] was a step in this direction. Specifically, the key characteristics of a fire such as the evolution of the shape and velocity of a premixed flame front, as well as the flame run-up distance were predicted by combining the mechanism of flame acceleration due to a “finger” flame shape [4] with another flame acceleration mechanism [5] devoted to self-similar accelerative outward propagation of a global spherical flame front, which is wrinkled because of the Darrieus–Landau (DL) instability (see [5] for more details on the dynamic and scalar fields associated with flame propagation and the fluid ahead of the flame front). However, [3] assumed smooth walls of a passage, while it was subsequently recognized by [6] that obstructions, which are inevitably present in a coal mine (such as mining equipment, belt conveyor systems, and piles of rubble), can block a noticeable part of the passage, thereby providing a significant impact on a fire scenario. Therefore, there is, a critical need to account for obstacles in a predictive coal mining fire scenario as addressed in the present work.
Previously, Silvestrini et al. [7] developed a simplified formulation for evaluating the flame run-up distance for smooth tubes, and then extended it to account for obstacles, and obtained that the obstacles promote flame acceleration, thereby leading to shorter run-up distances. Ciccarelli and Dorofeev [8] reviewed the knowledge on flame acceleration in smooth and obstructed channels, emphasizing the formation of a cellular flame structure due to the DL instability, as observed in the small-scale experiments for the flames in the obstructed channels [8]. Houim and Oran [9] considered a channel with smooth walls, filled with a coal dust or ash particles layered at the bottom of the channel. They demonstrated, computationally, that the flame interaction with the coal dust formed hot spots ahead of the flame front; and autoignition of these hot spots produced the jumps in the flame position and sharp spikes in the flame velocity. Seshadri et al. [10] considered the structure of a premixed flame propagating in environments where combustible particles were uniformly distributed and identified the effective equivalence ratio and laminar burning velocity of the mixture in terms of particle size and concentration. Xie et al. [11] showed, experimentally, that a small size of the combustible coal particles promoted the overall equivalence ratio of the fuel mixture, as well as the laminar flame speed in this mixture. Zheng et al. [12] investigated the effect of a distribution of a methane-air fuel mixture composition inside a two-dimensional (2D) obstructed channel on the deflagration-to-detonation transition (DDT), with the finding that the inhomogeneous distributions either did not produce a detonation or resulted in a detonation that was decoupled into a flame and a shock earlier as compared with the homogeneous distribution. The wall restrictions in the confined enclosures, as well as the shape and layout of the obstacles can significantly influence flame acceleration in the obstructed channels. In unconfined spaces, Ogawa et al. [13] showed that an array of square obstacles led to continuous detonation propagation. Bychkov et al. [14] provided an explanation of a physical mechanism of flame acceleration in a 2D obstructed channel, equipped with a comb-shaped array of tightly placed obstacles. Specifically, delayed burning in the pockets between the obstacles caused a jet flow along the centerline of the channel, thereby providing extremely powerful flame acceleration [14]. Subsequently, Valiev et al. [15] extended a 2D study [14] into cylindrical-axisymmetric geometry, with much stronger flame acceleration occurring in the cylindrical tubes [15] as compared with the 2D channels [14]. Unlike finger-flame acceleration [4], acceleration due to obstacles [14,15] is unlimited in time and it can eventually trigger a DDT constituting an extra, shock-based disaster for personnel and equipment in underground enclosures such as coal mines or subway tunnels.
In the present work, the coal mining passage is imitated by an obstructed channel/tube [14,15]. Specifically, we combine the formulation of a globally-spherical, self-accelerating expanding flame front, distorted by the DL instability [5] and the formulation of fast flame acceleration due to the obstacles in a channel/tube [14,15], thus, revisiting a predictive scenario of a burning accident in a coal mining passage [3]. It is shown that the obstacles play a key role in the coal mining fire scenario, and that the DL instability can significantly promote obstacle-based acceleration.
Starting with a homogeneous methane-air fuel mixture, the analysis is subsequently extended to a gaseous-dusty environment using a modified Seshadri formulation [3,10]. The parametric study involves variations of the blockage and equivalence ratios, α and ϕ , as well as the size of the dust particles and their type (inert such as sand, combustible such as coal, and their combination).

2. Formulation

Here, we present our formulation, starting with the 2D geometry in Section 2.1, and then its counterpart, the approach for cylindrical-axisymmetric geometry in Section 2.2. First, we recall the basics of flame acceleration mechanism for a globally-spherical, expanding flame front corrugated due to the DL instability. Indeed, any large-scale premixed flame front is prone to this mode of instability such that the radius of a globally-spherical, expanding flame obeys the time-dependent power law [5]:
R f ( t ) = R 0 + C t n C t n ,
where a small quantity R 0 can be considered to be a fitting parameter that can account for the ignition uncertainty in the experiments or a critical flame radius at which self-similar acceleration visually starts [5]. In the present work, we neglect this parameter as negligible as compared with a passage height or an obstacle length, R 0 α H < H . The factor C in Equation (1) can be evaluated as [3,5]:
C = k D L n 1 ( Θ U f / n ) n ,
with the thermal expansion ratio Θ ρ f u e l / ρ b u r n t and the DL cutoff wavenumber k D L 2 π / λ D L ( λ D L is the DL critical wavelength), which is assumed to obey the relation (see [16] for details):
λ D L = 2 π L f ( 1 + ( Θ + 1 ) ( Θ 1 ) 2 Θ ln Θ ) ; k D L = L f 1 ( 1 + ( Θ + 1 ) ( Θ 1 ) 2 Θ ln Θ ) 1 ,
where L f D t h / U f is the conventional definition of the flame thickness, and D t h is the thermal diffusivity coefficient of the fuel mixture. In the present work, Equation (3) is chosen because it was used in our previous works, such as in [3], and therefore it is natural to use the same formula here, for better comparison. It is nevertheless noted that there exist alternative formulas for λ D L as discussed in [17] and references therein. In this light, the present formulation is flexible with respect to the choice of the formula for λ D L , such that Equation (3) can be readily replaced by any alternative. The DL cutoff wavelength λ D L is a quantity, proportional to the flame thickness, λ D L L f , with a proportionality factor of the order of 10 1 ~ 10 2 depending on thermal-chemical flame properties such as the expansion ratio Θ [16,17]. Equation (3) is an example of such a dependence. Obviously, if L f 0 , then λ D L 0 ,
In [5], the authors developed an analytical theory describing self-similar accelerative outward propagation of a globally-spherical flame front. It is wrinkled because of the DL instability and obeys Equation (1). In particular, the dynamic and scalar fields associated with flame propagation, and the flow ahead of the flame, were analyzed, with a particular focus on (although not limited to) the flame-generated compression waves and the time evolution of the radial flow velocity. The exponent n in Equation (1) has been reported to be in the range 4 / 3 n 3 / 2 [18,19,20,21,22]. Similar to [3], in the present work, we employed n = 1.4 in the majority of cases, although the impact of varying n was also analyzed. According to Equation (1), instead of the unstretched laminar flame velocity U f , the instantaneous radial flame velocity with respect to the fuel mixture is
U D L ( t ) = n C Θ t n 1 .

2.1. Two-Dimensional (2D) Geometry

We consider a 2D passage (channel) of width 2 H = 2.1   m as illustrated in Figure 1, which is closed at one end, with a premixed flame front propagating towards the open end. The passage is blocked by the obstacles of length, α H , such that the central segment of the channel of width, 2 ( 1 α ) H , remains unobstructed. Theoretically, we adopt a limit of tightly placed obstacles, Δ z α H , which allows treating the flow between the obstacles as laminar. Nevertheless, it is noted that the Bychkov model [14] is applicable even far beyond the restriction Δ z α H ; namely, it is valid as long as Δ z α H or even Δ z H as validated in the computational study [23]. Thus, it should be enough to keep the obstacle spacing less than the obstacle height, thereby having their ratio not exceeding unity, Δ z / α H < 1 (although the tighter the better, of course). In this respect, the assumption of a laminar flow between the obstacles can be waved, to some extent. For instance, the Eulerian assumption of a non-viscous flow might lead to a similar x-direction velocity profile without imposing tightly placed obstacles. However, in that case, there would be additional effects of a vertical motion between the obstacles, presumably, leading to pulsations, which would need to be studied. Therefore, as of now, we assume a laminar flow in the pockets between the obstacles.
With such a laminar flow approach, the flame front inside the pockets can be taken as locally planar all the time, thus, spreading in the x-direction with the instantaneous global flame velocity U D L ( t ) . As the burnt matter expands with a thermal expansion ratio Θ , the flow is pushed out of the “pockets”. Coming into a central part of the passage, the flow changes its direction and pushes the flame forward in the z-direction towards the exit. This creates a positive feedback between the flame and the pockets as the flame is pushed forward, thereby creating new pockets behind it.
It is emphasized that the flame dynamics and morphology are controlled by the processes occurring in the obstructed segment of the passage and not in its unobstructed, central part. Indeed, although the flame tip is driven by a powerful jet flow; such a jet flow is generated due to delayed burning in the pockets between the obstacles. Moreover, the total burning rate (which grows extremely fast) correlates well with the total surface area of the flame front, which is mainly provided by the flame segments in the pockets between the obstacles. Consequently, the shape (surface area) of the flame tip does not seem to play a role in the flame surface area, and therefore the total burning rate (even if the tip would be very corrugated, this would not change the total flame surface area and the total burning rate). For this reason, the shape of the flame tip can be assumed to be planar all the time [14,15]. With this assumption, we can potentially neglect the effect of the flame growth at the initial stages of burning. Nevertheless, we accounted for the evolution of flame “skirt” R f described by Demir et al. [3], with a point ignition in our analysis,
R f ( t ) = Θ H ( Θ 1 ) { 1 exp [ ( Θ 1 ) Θ H ( k D L n 1 ( Θ U f n ) n ) t n ] } .
With methane-air and propane-air laminar burning velocities not exceeding U f = 0.3 ~ 0.4   m / s [24] and the respective speeds of sound circa c 0 ~ 350   m / s , the initial Mach number associated with flame propagation appears to be very small, M a U f / c 0 ~ 10 3 1 . This guides us to employ a low Mach number model of an incompressible flow, M a 0 , similar to [3]. This model has been widely employed analytically and verified by fully compressible computational simulations in our previous works [4,15]. The model has been justified to work well at the initial, essentially subsonic stage of flame acceleration, and to provide reasonable evaluations even in the case of compressible flows. Nevertheless, it is recognized that the model would break when the burning velocity approaches near-sonic values. Such a limitation of the model and the way to overcome is commented on in Section 3.1. For a 2D geometry, the incompressible continuity equation reads as:
u x x + u z z = 0
Similar to the Bychkov model [14,15], here, we employed the assumption of a potential flow in the burnt gas, which has been validated by the computational simulations in [15]. It is nevertheless noted that, in a real passage, obstacles generate vorticity at least in the central region close to the obstacles. At the same time, according to [23], although the flow in the burnt gas is generally rotational in practice, because of the curved flame shape, it can be treated as potentially close to the wall when the flame front is locally planar, which is enough to certify the Bychkov approach [23].
The question about vorticity and rotational region has been addressed, in part, in [23,25]. According to the simulations [25], vorticity provides a notable impact on the Bychkov mechanism only if the obstacle spacing exceeds the passage radius, Δ Z > R , with the Bychkov model being fully broken if Δ Z = 2 R . Obviously, such Δ Z is much larger than what usually means tightly placed obstacles, Δ Z α R , say Δ Z / R = 0.2 ~ 0.6 [23], thereby justifying the potential flow approach. As a result, considering the potential flow in the free part of the passage, with the boundary condition | u x | = ( Θ 1 ) U D L ( t ) at x = ± ( 1 α ) H , we find the velocity profiles in the burnt gas,
( u x ; u z ) = ( Θ 1 ) U D L ( t ) ( 1 α ) H ( x ; z ) .
The flame tip Z t i p obeys the evolution equation
d Z t i p d t = u z ( Z t i p ) + Θ U D L ( t ) .
which mathematically expresses the fact that the flame tip velocity in the laboratory reference frame equals the velocity of the burnt matter flow taken at the flame tip location, u z ( Z t i p ) , plus the flame tip velocity with respect to this flow, Θ U D L ( t ) . Accounting for Equation (7), Equation (8) can be rewritten as:
d Z t i p d t = ( Θ 1 ) ( 1 α ) U D L ( t ) H Z t i p ( t ) + Θ U D L ( t ) .
While both terms in the right-hand-side (RHS) of Equation (9) are, initially, of the same order, the first (“flow”) term soon dominates over the second (“flame”) term, because the value Z t i p ( t ) grows promptly with time. With the initial condition Z t i p ( 0 ) = 0 , the solution to Equation (9) reads as:
Z t i p = Θ ( 1 α ) H ( Θ 1 ) { exp [ ( Θ 1 ) ( 1 α ) H C Θ t n ] 1 } = Θ ( 1 α ) H ( Θ 1 ) { exp [ ( Θ 1 ) ( 1 α ) H ( k D L Θ ) n 1 ( U f n ) n t n ] 1 }
Then, the velocity of the flame tip is given by:
d Z t i p d t = U t i p = n C t n 1 exp [ ( Θ 1 ) ( 1 α ) H C Θ t n ] .
In the case of n = 1 , the DL instability disappears and Equations (7)–(11) reproduce their counterparts of [14]. We next recognize that the DL instability should not provide an increase in the flame surface/velocity infinitely because the characteristic cross-sectional size of the flame eventually stops growing. Specifically, the characteristic length scale of the flame radius does not grow further after a flame skirt touches an obstacle in the case of point ignition. This instant t o b s can be calculated from a condition R f ( t o b s ) = ( 1 α ) H in Equation (5), namely:
t o b s = n Θ U f { Θ H k D L 1 n Θ 1 ln [ Θ 1 + α ( Θ 1 ) ] } 1 / n
after which the instantaneous global flame velocity U D L remains constant. The respective flame tip position and the instantaneous global flame velocity at this instant are
Z t i p ( t o b s ) = Θ ( 1 α ) H ( Θ 1 ) [ ( Θ 1 + α ( Θ 1 ) ) 1 / ( 1 α ) 1 ] ,
U D L ( t o b s ) = n C Θ t o b s n 1 .
Then, instead of the time-dependent quantity U D L ( t ) in Equation (9), we have a constant value U D L ( t o b s )
d Z t i p d t = ( Θ 1 ) ( 1 α ) U D L ( t o b s ) H Z t i p + Θ U D L ( t o b s ) .
Integrating Equation (15) for t > t o b s , with the matching condition (13), we find
Z t i p = Θ ( 1 α ) H ( Θ 1 ) { [ Θ 1 + α ( Θ 1 ) ] 1 / ( 1 α ) exp [ ( Θ 1 ) ( 1 α ) H U D L ( t o b s ) ( t t o b s ) ] 1 } ,
d Z t i p d t = U t i p = Θ U D L ( t o b s ) [ Θ 1 + α ( Θ 1 ) ] 1 / ( 1 α ) exp [ ( Θ 1 ) ( 1 α ) H U D L ( t o b s ) ( t t o b s ) ] .
We also determine the flame run-up time, t r u d , which we evaluate as the time at which the flame velocity reaches the speed of sound of the reactants, c 0 . Namely, Equation (17) gives the run-up time as:
t r u d = t o b s + ( 1 α ) H ( Θ 1 ) U D L ( t o b s ) ln { c 0 Θ U D L ( t o b s ) [ 1 + α ( Θ 1 ) Θ ] 1 / ( 1 α ) } ,
with the corresponding run-up distance Z r u d being
Z r u d = Θ ( 1 α ) H ( Θ 1 ) { [ Θ 1 + α ( Θ 1 ) ] 1 / ( 1 α ) exp [ ( Θ 1 ) ( 1 α ) H U D L ( t o b s ) ( t r u d t o b s ) ] 1 }

2.2. Cylindrical Geometry

We next derive a similar formulation for the cylindrical geometry. Similarly, the passage radius is taken as R = 1.05   m . Evolution of the flame skirt in an unobstructed section is described as [3]:
R f ( t ) = Θ R β tanh ( β Θ C t n R ) ,
where β = Θ ( Θ 1 ) . The incompressible continuity equation in a cylindrical-axisymmetric geometry reads:
1 r r ( r u r ) + u z z = 0 ,
with a symmetric boundary condition | u r | = ( Θ 1 ) U D L ( t ) at r = ± ( 1 α ) R . Assuming that the flow in the unobstructed part of the passage is potential (and the shape of the flame front is not of significant importance), the velocity profile reads:
( u r ; u z ) = ( Θ 1 ) U D L ( t ) ( 1 α ) R ( r ; 2 z ) .
The propagation velocity of the flame tip with respect to burnt gas is given by:
d Z t i p d t = u z ( Z t i p ) + Θ U D L ( t ) = 2 ( Θ 1 ) U D L ( t ) ( 1 α ) R Z t i p + Θ U D L ( t ) .
Furthermore, an additional increase in the flame shape occurs in a cylindrical geometry, because the flame in an axisymmetric pocket expands with the radius as [15]:
R f , o = ( 1 α ) R + U f [ t t f ( z ) ] ,
where t f ( z ) is the instant at which the fresh gas in the pocket between the obstacles at the position z starts burning and R f , o represents the radial coordinate of flame skirt in an obstructed passage. Averaging the last term in Equation (24) as t t f ( z ) = ( 1 α ) R / 2 U f ( Θ 1 ) (see [15] for details), we arrive at the following modified version of the evolution Equation (23) for Z t i p :
d Z t i p d t = 2 ( Θ 1 ) U D L ( t ) ( 1 α ) R ( 1 + 1 2 ( Θ 1 ) ) Z t i p + Θ U D L ( t ) = ( 2 Θ 1 ) U D L ( t ) ( 1 α ) R Z t i p + Θ U D L ( t ) .
Integrating Equation (25) with the initial condition Z t i p ( 0 ) = 0 , we find the solution in the form:
Z t i p = Θ ( 1 α ) R ( 2 Θ 1 ) { exp [ ( 2 Θ 1 ) ( 1 α ) R C Θ t n ] 1 } ,
with the velocity of the flame tip being
d Z t i p d t = U t i p = n C t n 1 exp [ ( 2 Θ 1 ) ( 1 α ) R C Θ t n ] .
Again, as in the 2D case, if n = 1 , then the DL instability disappears, and the present formulation reproduces that of [15]. Furthermore, similar to the 2D formulation of the previous subsection, the average-in-mean radial flame size does not grow after the instant t o b s , which can be determined from the condition R f ( t o b s ) = ( 1 α ) R in Equation (20), namely:
t o b s = n Θ U f { Θ R k D L 1 n 2 β ln [ Θ + ( 1 α ) β Θ ( 1 α ) β ] } 1 / n .
The respective flame tip position and the instantaneous global flame velocity at this instant are
Z t i p ( t o b s ) = Θ ( 1 α ) R ( 2 Θ 1 ) { [ Θ + ( 1 α ) β Θ ( 1 α ) β ] ( 2 Θ 1 ) 2 ( 1 α ) β 1 }
U D L ( t o b s ) = n C Θ t o b s n 1 .
Then, instead of the time-dependent quantity U D L ( t ) in Equation (25), we have a constant value U D L ( t o b s )
d Z t i p d t = ( 2 Θ 1 ) U D L ( t o b s ) ( 1 α ) H + Θ U D L ( t o b s ) .
Integrating the Equation (31) for t > t o b s , with the matching condition (29), we find
Z t i p = Θ ( 1 α ) R ( 2 Θ 1 ) { [ Θ + ( 1 α ) β Θ ( 1 α ) β ] ( 2 Θ 1 ) 2 ( 1 α ) β exp [ ( 2 Θ 1 ) ( 1 α ) R U D L ( t o b s ) ( t t o b s ) ] 1 } ,
d Z t i p d t = U t i p = Θ U D L ( t o b s ) [ Θ + ( 1 α ) β Θ ( 1 α ) β ] ( 2 Θ 1 ) 2 ( 1 α ) β exp [ ( 2 Θ 1 ) ( 1 α ) R U D L ( t o b s ) ( t t o b s ) ] .
Similar to the 2D analysis, from the condition U t i p ~ c 0 the run-up time in the cylindrical geometry can be evaluated from Equation (33) as:
t r u d = t o b s + ( 1 α ) R ( 2 Θ 1 ) U D L ( t o b s ) ln { c 0 Θ U D L ( t o b s ) [ Θ ( 1 α ) β Θ + ( 1 α ) β ] ( 2 Θ 1 ) 2 ( 1 α ) β } ,
with the corresponding run-up distance
Z r u d = Θ ( 1 α ) R ( 2 Θ 1 ) { [ Θ + ( 1 α ) β Θ ( 1 α ) β ] ( 2 Θ 1 ) 2 ( 1 α ) β exp [ ( 2 Θ 1 ) ( 1 α ) R U D L ( t o b s ) ( t r u d t o b s ) ] 1 } .
Finally, it is recognized that a number of effects remained beyond consideration in the present work. For instance, only the DL flame instability is accounted for as the agent triggering unsteadiness, without considering turbulence or other combustion instability modes such as the diffusional-thermal (DT), Rayleigh–Taylor (RT), or Kevin–Helmholtz (KH) instabilities. In the future, these features could be potentially incorporated into the analysis by replacing the quantities λ D L (Equation (3)) and U D L ( t ) (Equation (4)) by the respective DT/RT/KH/turbulence-induced parameters. In addition, analysis of partly open obstructed duct, as in the pioneering Taylor–Bimson (TB) model [26], would be of interest. However, all such extensions of our work will be performed elsewhere.

3. Results and Discussion

In this section, we present and discuss the results of an intensive parametric study performed. Specifically, the first subsection is devoted to homogeneously-gaseous, methane-air and propane-air burning with the thermal-chemical parameters (such as Θ and U f ) being the functions of the equivalence ratio ϕ , as tabulated in Table 1; see also [3,24]. In the second subsection, the analysis is extended to a gaseous-dusty environment.

3.1. Gaseous Combustion

We start with Figure 2 presenting the time evolutions of the flame tip position, Z t i p , Figure 2a, and its velocity, U t i p , Figure 2b, in a 2D geometry, for stoichiometric methane (CH4)-air burning, as well as that of propane (C3H8)-air, for comparison, with various blockages α = 0 , 1 / 3 , 1 / 2 , 2 / 3 employed. The horizontal dotted lines in Figure 2b show the speeds of sound for the methane-air, 354   m / s , and propane-air, 340   m / s , mixtures (Table 1). The case of no obstacles, α = 0 , 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 t ~ 0.089   s and 0.072   s for methane-air and propane-air burning, respectively. It is also noted that by these times, the propane-air flame overcomes the sound threshold 340   m / s , whereas the methane-air flame stops at U t i p   ~ 288   m / s . In contrast, in an obstructed channel, α > 0 , 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 M a 0 , and then extending it to account for the finite M a 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 U f 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 M a , 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 α = 0 and methane-air burning in a 2D geometry, here, at the same conditions, the flame would overcome the sound speed at t ~ 0.058   s , i.e., slightly prior the deceleration stage starting at t ~ 0.0735   s . 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 n of Equation (1) to be the same as that in [3], i.e., n = 1.4 . 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 n = 1 by default). It is noted, in this respect, that an appropriate choice of the factor n 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 n 1.5 , the Princeton experiments in a dual chamber at elevated pressures suggested n = 1.3 ~ 1.33 ; 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 n . Moreover, the quantity of n 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 n in the range 4 / 3 n 3 / 2 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 n = 1.33 ,   1.4   and   1.5 , and both 2D (Figure 5a,b) and cylindrical (Figure 5c,d), geometries are studied. It is seen that the variations of n impact the flame position, velocity and acceleration substantially, although the effect is quantitative but not qualitative. As expected, flame acceleration in the case of n = 1.5 proceeds noticeably faster than for n = 1.4 , whereas a flame with n = 1.33 accelerates noticeably slower as compared with that with n = 1.4 (of course, provided that other combustion characteristics are kept the same). Nevertheless, without the final answer to the question about an appropriate choice for n , in the rest of this work we keep using its median value, n = 1.4 , 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 ( ϕ = 1 ) gaseous methane-air combustion considered in Figure 2, Figure 3 and Figure 4 to the equivalence ratios in the range 0.6 ϕ 1.4 , see Figure 6, Figure 7, Figure 8 and Figure 9. In particular, Figure 6 presents the 2D configuration with various α = 0 , 1 / 3 , 1 / 2 , 2 / 3 for ϕ = 0.8 ,   1 ,   1.2 . It is seen that a slightly fuel-lean flame with ϕ = 0.8 accelerates much slower than the ϕ 1 flames, especially in the cases of α = 1 / 3 and 1 / 2 . This is because of a much lower U f (and thereby higher L f and lower k D L ) 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 ϕ = 0.8 flame are evaluated as t r u d ~ 0.18   s for α = 1 / 3 and t r u d ~ 0.1   s for α = 1 / 2 . In the case of no obstacles, α = 0 , no detonation is predicted for a ϕ = 0.8 flame in a 2D passage. Overall, among all equivalence ratios considered, the fastest flame acceleration is observed for a slightly rich flame of ϕ = 1.1 .
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 ϕ = 0.6 and ϕ = 1.4 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 ϕ = 0.6 flame stops accelerating after advancing 5.8   m and attaining the maximal velocity of 24.6   m / s , 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 ϕ = 0.6 flames than for the ϕ = 1.4 flames: the lines for ϕ = 0.6 with various α go more widely than those for ϕ = 1.4 . It is also seen in Figure 7 that the rich flames, ϕ = 1.4 , accelerate faster than the lean flames, ϕ = 0.6 . 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 ϕ = 0.8 flames appear t r u d ~ 0.1   s for α = 1 / 3 and t r u d ~ 0.08   s for α = 1 / 2 . A key difference between the two geometries is that whereas no detonation is predicted for the ϕ = 0.8 flame in a 2D case without obstacles, α = 0 , in the cylindrical configuration with α = 0 , the ϕ = 0.8 methane-air flame was able to reach the speed of sound, 353   m / s , at t ~ 0.12   s , thereby making a detonation possible. Again, further away from stoichiometry, see Figure 9, the ϕ = 0.6 and ϕ = 1.4 flames accelerate slower and the acceleration time lasts longer as compared to the ϕ = 0.8 ~ 1.2 flames in Figure 8. For example, in the case of α = 0 , the ϕ = 0.6 flame stops accelerating when advancing 11.2   m and reaches the maximal velocity of 84.7   m / 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 ( ϕ = 1.4 ) accelerate faster than the highly lean flames ( ϕ = 0.6 ) in the cylindrical case, and the impact of obstacles seems to be more important for the ϕ = 0.6 flames (the lines corresponding to various blockage ratios go more widely).
Using Table 1, we next analyze the run-up times t r u d , Equations (18) and (34), and the respective run-up distances Z r u d , Equations (19) and (35), for the methane (CH4)-air and propane (C3H8)-air flames of various equivalence ratios. Specifically, Figure 10 presents Z r u d versus ϕ for various blockage ratios, including the case of no obstacles, α = 0 , 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 ϕ ~ 1.1 . In the 2D case, Figure 10a, we have Z r u d ~ 8.19   m , 7   m , and 5.46   m for α = 1 / 3 , 1 / 2 , and 2 / 3 , respectively. The case of α = 0 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 80   m for ϕ = 0.6 and up to 35   m for ϕ = 1.4 (still in the 2D geometry). In the cylindrical-axisymmetric configuration, for the fastest methane-air flames with ϕ = 1.1 we found the run-up distances as small as Z r u d ~ 5.31   m , 4.11   m , 3.45   m , and 2.64   m for α = 1 / 3 , 1 / 2 , and 2 / 3 , respectively. For lean or rich methane-air burning, the run-up distances are much higher in the cylindrical geometry, i.e., up to 40   m for ϕ = 0.6 and up to 18   m for ϕ = 1.4 . 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 α = 0 , when the equivalence ratio is in the range 0.8 ϕ 1.3 , which generally agrees with [3]. Overall, for the same geometry, α and ϕ , the run-up distances are dramatically shorter for the C3H8-air flames as compared with the CH4-air flames; this would overcome the sound threshold for additional equivalence ratios in the case of α = 0 , being within the ranges 1 ϕ 1.2 and 0.7 ϕ 1.4 for the 2D and cylindrical-axisymmetric geometries, respectively.

3.2. Validation of Gaseous Formulation

We have compared our analytical predictions with the data available in the literature such as the experiments by [31] and computational simulations (FAST and ALLA) in [6]. The results are shown in Figure 11 and Figure 12, where the flame tip velocity is shown versus its position, for various blockage ratios, such that the linear trend represents exponential acceleration. The simulation results [6] are shown by circles while the square markers present the experimental data [31]. The solid/dashed lines are devoted to the theoretical predictions. We observe reasonably good quantitative agreement between the 2D theory and the literature data. As for the cylindrical theory, our qualitative linear theoretical trend fits that of the markers very well.

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 U f U d , f [3]:
U d , f = U f ϕ s ϕ C P C T ( T f T b ) 2 ( T b T u T f T u ) exp ( E ( T f T b ) T f T b R u )
where ϕ s is the modified equivalence ratio of the gaseous-dusty-air mixture in the presence of combustible dust particles:
ϕ s = [ ( ( m f u e l m ) / ( M C H 4 ) ) / ( ( m a i r m ) / ( M a i r ) ) ] a c t [ ( ( m C H 4 m ) / ( M C H 4 ) ) / ( ( m a i r m ) / ( M a i r ) ) ] s t .
M C H 4 and M a i r are the respective molar masses; m C H 4 m , m a i r m and m f u e l m are the original masses per unit volume for a given equivalence ratio; C T = C P + C s n s V s ρ s / ρ is the specific heat of the whole mixture, containing the components for the gas, C P , and dust particles, C s ; ρ s is the density of a single dust particle, while ρ = ρ u + c s is that for the gaseous-dusty fuel-air mixture, with the density of the gas ρ u and the concentration of the dust particles c s ; n s = ( c s / ρ s ) / V s is the number of particles per unit volume, with V s = 4 π r s 3 / 3 being the volume of a single particle of radius r s ; R u = 8.314   J / ( mol · K ) is the universal gas constant; T u = 300   K is the reactants temperature and T b is the adiabatic flame temperature based on the purely methane-air equivalence ratio. Similar to [3], here, T b is calculated as a fifth-order polynomial function of the equivalence ratio ϕ [32]:
T b = ( 2.21 × 10 4 ) ϕ 5 + ( 8.042 × 10 4 ) ϕ 4 + ( 1.171 × 10 5 ) ϕ 3 + ( 8.471 × 10 4 ) ϕ 2 + ( 2.854 × 10 4 ) ϕ + 4.89 × 10 3
valid in the range 0.6 ϕ 1.6 [32]. Knowing ϕ s from Equation (37) and calculating a new flame temperature for the dusty-gaseous mixture, T f * , one can find U d , f , 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 ϕ 1 , methane-air combustion, the global chemical reaction is described by the equation:
ϕ C H 4 + 2 ( O 2 + 3.76 N 2 ) ϕ C O 2 + 2 ϕ H 2 O + 7.52 N 2 + 4 ( 1 ϕ ) O 2
such that the heat release in the process of burning of ϕ moles of C H 4 and 2 × ( 1 + 3.76 ) = 9.52 moles of air is given by [3,11]:
Q ϕ = ( T b T u ) C p   n p r o d u c t ,
where n p r o d u c t 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 C H 4 -air combustion of a given ϕ as follows [3]:
Q = [ ( T b T u ) C P . n p r o d u c t ] n a i r 9.52 ( V C H 4 + V a i r ) .
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]:
Q = [ ( T f * * T u ) C P . n p r o d u c t ] n a i r 9.52 ( V C H 4 + V a i r ) + c s C s ( T f * * + T u ) + L v ,
where L ν = 0.01 w ν Δ h C H 4 [11] is the heat of gasification per unit volume. From the last equation, the secondly revised flame temperature, T f * * , is calculated as:
T f * * = Q L v n a i r 9.52 ( V C H 4 + V a i r ) C p n p r o d u c t + c s C s + T u .
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, T f * * * = ( T f * + T f * * ) / 2 . Similar to a combustible dust case, T f * * and T f * * * are used in Equation (36) to find a new laminar flame speed U d , f 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 λ D L 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, λ D L is considered to be the quantity devoted to such an effective fluid. Consequently, λ D L is given by the same formula, Equation (3), although the variables in this formula ( Θ , L f ) are adjusted due to the presence of dust particles.
We can also provide another justification for this approach, namely, because λ D L L f , if the flame thickness does not change much when the dust particles are laden, then the value λ D L , 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 c s = 50 g / m 3 and of dust particles radius r s = 75   μ m in the 2D (a, b) and cylindrical-axisymmetric (c, d) geometries. We employed the lean ( ϕ = 0.7 ) methane-air fuel mixture and various blockage ratios α , including the case of no obstacles, α = 0 . 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 c s = 120   g / m 3 keeping the same particle radius, r s = 75   μ m , 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 r s = 75   μ m . The impact of the blockage ratio is noticeable, especially in the cylindrical geometry. The flame velocity in the case of inert particles and α = 2 / 3 is equivalent to the event of combustible dust but with no obstacles, α = 0 , for t ~ 0.136   s ; 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 c s = 250   g / m 3 , with the particle size r s = 75   μ m , 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 c s = 50 and 120   g / m 3 , the combustible particles promote flame acceleration, whereas the combined (combustible + inert) and inert particles also suppress it in the situation of c s = 250   g / m 3 .
Figure 16a–d is the counterpart of Figure 14a–d, respectively, for a smaller particle radius, r s = 10   μ m . 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 38   m / s for the particles of size r s = 75   μ m , in the case of r s = 10   μ m , the sound threshold of 352   m / s for ϕ = 0.7 methane-air burning was reached in the combustible coal gaseous-dusty environment with α = 2 / 3 at the time instant t ~ 0.114   s (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, 1 / 3 α 2 / 3 , we obtained faster acceleration than in both respective cases of other dusts and no dust (Figure 16d). In contrast to the r s = 75   μ m case in Figure 14, the combined combustible-inert particles promoted flame acceleration in the case of r s = 10   μ m , 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 r s = 75   μ m in Figure 14, they also suppress flame acceleration for r s = 10   μ m . The aforementioned effects of all particles, i.e., combustible, inert, and combined, and particle sizes are found to grow with the blockage ratio α .

4. Conclusions

Therefore, we have undertaken a step towards a predictive scenario of a burning accident in obstructed coal mining passage by developing a theoretical formulation, which combines the mechanism of flame acceleration in obstructed pipes [14,15] with that due to the DL instability [5]. Specifically, the 2D planar and cylindrical-axisymmetric geometries were considered, along with an assumption of incompressible flow for various passage configurations and the compositions of the combustible premixture. The newly identified flame propagation scenario was studied in terms of the evolution of the flame tip position Z t i p and its velocity (in the laboratory reference frame) U t i p . The role of the obstacles, as well as the DL instability, on a fire scenario are found to be significant, with a stronger effect observed in the cylindrical geometry. Closer to the stoichiometric condition, a flame propagates and accelerates faster; and acceleration is also promoted with an increase in the blockage ratio α . Starting with homogeneously gaseous combustion, the analysis was subsequently extended to incorporate dust particles in the passage. Namely, inert and combustible dust, as well as their combination, were considered. We found that the combustible dust particles of radii 10 ~ 75   μ m facilitate flame acceleration, while acceleration is moderated by the inert particles. The effect of particle size is significant in a manner that smaller particles lead to faster flame acceleration. While a combination of combustible and inert dust, of dust particle radius r s = 75   μ m , moderates flame acceleration as compared with purely gaseous combustion, a combustible-inert combination of the particles of smaller radius, r s = 10   μ m , promotes flame acceleration. The effect of a type and size of the particles increases with the blockage ratio.

Author Contributions

Conceptualization, F.K. and V.A.; funding acquisition, V.A.; investigation, F.K. and S.D.; methodology, F.K. and V.A.; project administration, V.A.; resources, V.A.; supervision, V.A. and D.V.; writing, F.K., V.A., and D.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by the Alpha Foundation for the Improvement of Mine Safety and Health, Inc., through Award #ASTI14FO-82 (V.A.), prior to which it was funded by the U.S. National Science Foundation (NSF) through CAREER Award No. 1554254 (V.A.), as well as by the West Virginia Higher Education Policy Commission through Grant #HEPC.dsr.18.7 (V.A.). The views, opinions, and recommendations expressed in this work are solely those of the authors and do not imply any endorsement by the Alpha Foundation, its directors and staff.

Acknowledgments

The authors dedicate this work to deceased Vitaly Bychkov (1968–2015) of Umea University, Sweden, whose role in analytical studies of premixed flame acceleration is hard to overestimate.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

R f radius of a global spherical expanding flame front
C constant defined in Equation (1)
t time
n Darrieus–Landau instability exponent
k D L Darrieus–Landau cutoff wavenumber
U f laminar flame velocity
ρ f u e l density of fuel mixture
ρ b u r n t density of burnt gas
L f flame thickness
D t h thermal diffusivity coefficient
U D L instantaneous radial flame velocity
H half-width of a two-dimensional (2D) passage
R radius of a cylindrical passage
Δ z obstacle spacing
x , rradial direction
z axial direction
R f ( t ) flame “skirt”
Δ z obstacle spacing
u x , u r radial velocity
u z axial velocity
Z t i p flame tip position
U t i p flame tip velocity
t o b s the time flame skirt touches an obstacle
R f ( t o b s ) flame “skirt” radius at t o b s
Z t i p ( t o b s ) flame tip position at t o b s
U D L ( t o b s ) global flame velocity at t o b s
t r u d flame run-up time
c 0 speed of sound
Z r u d flame run-up distance
R f , o flame “skirt” in [15]
t f ( z ) the time instant at which the fresh gas between obstacles at the position z starts burning
U d , f laminar flame velocity in an “effective” gaseous-dusty environment
C p specific heat of gaseous air-fuel mixture
C T entire specific heat
T f flame temperature with particles
T b adiabatic flame temperature based on purely methane-air equivalence ratio
T u unburnt gas temperature
E activation energy
R u gas constant
M molar masses
m original masses
C s specific heat of dust particles
n s number of particles per unit volume
V s volume of a single particle
ρ s density of dust
ρ density of gaseous-dusty fuel-air mixture
ρ u density of gas
c s concentration of the particles
r s radius of a particle
Q volumetric heat release
L v heat of gasification per unit volume
n a i r number of moles of air per unit volume
V C H 4 volume of methane
V a i r volume of air
n p r o d u c t number of moles of the burning products
α blockage ratio
ϕ equivalence ratio
λ D L Darrieus–Landau critical wavelength
Θ ρ u / ρ b thermal expansion ratio
β defined as Θ ( Θ 1 )
ϕ s modified equivalence ratio in the gaseous-dusty air mixture
f flame
D L Darrieus–Landau
f u e l fuel mixture
b u r n t burnt gas
t h thermal
t i p tip
o b s obstacle
r u d run-up distance
0 initial
f , o flame in obstructed passage
d dust
s particle
T total
b burnt
u unburnt
a c t actual
s t stoichiometric
p r o d u c t product

References

  1. Ramani, R.V. Mining disasters caused and controlled by mankind: The case for coal mining and other mineral—Part 1: Causes of mining disasters. Nat. Resour. Forum 1995, 19, 233–242. [Google Scholar] [CrossRef]
  2. Lowen, M. Turkey Mine Disaster: Raw Anger in Soma a Year on. BBC World News. 13 May 2015. Available online: https://www.bbc.com/news/world-europe-32709431 (accessed on 1 June 2020).
  3. Demir, S.; Bychkov, V.; Chalagalla, S.H.R.; Akkerman, V. Towards a predictive scenario of a burning accident in a mining passage. Combust. Theory Model. 2017, 21, 997–1022. [Google Scholar] [CrossRef]
  4. Bychkov, V.; Akkerman, V.; Fru, G.; Petchenko, A.; Eriksson, L.-E. Flame acceleration at the early stages of burning in tubes. Combust. Flame 2007, 150, 263–276. [Google Scholar] [CrossRef]
  5. Akkerman, V.; Law, C.K.; Bychkov, V. Self-similar accelerative propagation of expanding wrinkled flames and explosion triggering. Phys. Rev. E 2011, 83, 026305. [Google Scholar] [CrossRef]
  6. Oran, E.S. Numerical Tools for Mitigation of Methane Explosions in Coalmines; Grant No. AFC215-20; The Alpha Foundation for the Improvement of Mine Safety and Health, Inc.: Philadelphia, PA, USA, 2017. [Google Scholar]
  7. Silvestrini, M.; Genova, B.; Parisi, G.; Trujillo, F.J.L. Flame acceleration and DDT run-up distance for smooth and obstacles filled tubes. J. Loss Prev. Process Ind. 2008, 21, 555–562. [Google Scholar] [CrossRef]
  8. Ciccarelli, G.; Dorofeev, S. Flame acceleration and transition to detonation in ducts. Prog. Energy Combust. Sci. 2008, 34, 499–550. [Google Scholar] [CrossRef]
  9. Houim, R.W.; Oran, E.S. Structure and flame speed of dilute and dense layered coal-dust explosion. J. Loss Prev. Process Ind. 2015, 36, 214–222. [Google Scholar] [CrossRef]
  10. Seshadri, K.; Berlad, A.L.; Tangirala, V. The structure of premixed particle-cloud flames. Combust. Flame 1992, 89, 333–342. [Google Scholar] [CrossRef]
  11. Xie, Y.; Raghavan, V.; Rangwala, A.S. Study of interaction of entrained coal dust particles in lean methane-air premixed flames. Combust. Flame 2012, 159, 2449–2456. [Google Scholar] [CrossRef] [Green Version]
  12. Zheng, W.; Kaplan, C.R.; Houim, R.W.; Oran, E.S. Flame acceleration and transition to detonation: Effects of a composition gradient in a mixture of methane and air. Proc. Combust. Inst. 2019, 37, 3521–3528. [Google Scholar] [CrossRef]
  13. Ogawa, T.; Gamezo, V.; Oran, E.S. Flame acceleration and transition to detonation in an array of square obstacles. J. Loss Prev. Process Ind. 2013, 26, 355–362. [Google Scholar] [CrossRef]
  14. Bychkov, V.; Valiev, D.; Eriksson, L.-E. Physical mechanism of ultrafast flame acceleration. Phys. Rev. Lett. 2008, 101, 164501. [Google Scholar] [CrossRef] [PubMed]
  15. Valiev, D.; Bychkov, V.; Akkerman, V.; Law, C.K.; Eriksson, L.-E. Flame acceleration in channels with obstacles in the deflagration-to-detonation transition. Combust. Flame 2010, 157, 1012–1021. [Google Scholar] [CrossRef] [Green Version]
  16. Pelce, P.; Clavin, P. Influence of hydrodynamics and diffusion upon the stability limit of laminar premixed flames. J. Fluid Mech. 1982, 124, 219–237. [Google Scholar] [CrossRef]
  17. Akkerman, V.; Bychkov, V. Velocity of weakly turbulent flames of finite thickness. Combust. Theory Model. 2005, 9, 323–351. [Google Scholar] [CrossRef]
  18. Gostintsev, Y.; Istratov, A.; Shulenin, Y. Self-similar propagation of a free turbulent flame in mixed gas mixtures. Combust. Explos. Shock Waves 1988, 24, 563–569. [Google Scholar] [CrossRef]
  19. Bradley, D.; Cresswell, T.; Puttock, J.S. Flame acceleration due to flame-induced instabilities in large-scale explosions. Combust. Flame 2001, 124, 551–559. [Google Scholar] [CrossRef]
  20. Molkov, V.V.; Makarov, D.V.; Schneider, H. Hydrogen-air deflagrations in open atmosphere: Large eddy simulation analysis of experimental data. Int. J. Hydrog. Energy 2007, 32, 2198–2205. [Google Scholar] [CrossRef]
  21. Kim, W.K.; Mogi, T.; Kuwana, K.; Dobashi, R. Self-similar propagation of expanding spherical flames in large-scale gas explosions. Proc. Combust. Inst. 2015, 35, 2051–2058. [Google Scholar] [CrossRef]
  22. Jomaas, G.; Law, C.K.; Bechtold, J. On transition to cellularity in expanding spherical flames. J. Fluid Mech. 2007, 583, 1–26. [Google Scholar] [CrossRef]
  23. Ugarte, O.; Bychkov, V.; Sadek, J.; Valiev, D.; Akkerman, V. Critical role of blockage ratio for flame acceleration in channels with tightly-spaced obstacles. Phys. Fluids 2016, 28, 093602. [Google Scholar] [CrossRef]
  24. Davis, S.G.; Quinard, J.; Searby, G. Markstein numbers in counterflow, methane- and propane-air flames: A computational study. Combust. Flame 2002, 130, 123–136. [Google Scholar] [CrossRef]
  25. Adebiyi, A.; Alkandari, R.; Valiev, D.; Akkerman, V. Effect of surface friction on ultrafast flame acceleration in obstructed cylindrical pipes. Aip Adv. 2019, 9, 035249. [Google Scholar] [CrossRef] [Green Version]
  26. Taylor, P.; Bimson, S.J. Flame propagation along a vented duct containing grids. Proc. Combust. Inst. 1988, 22, 1355–1362. [Google Scholar] [CrossRef]
  27. Demir, S.; Calavay, A.R.; Akkerman, V. Influence of gas compressibility on a burning accident in a mining passage. Combust. Theory Model. 2018, 22, 338–358. [Google Scholar] [CrossRef]
  28. Akkerman, V.; Valiev, D. Moderation of flame acceleration in obstructed cylindrical pipes due to Gas Compression. Phys. Fluids 2018, 30, 106101. [Google Scholar] [CrossRef]
  29. Kodakoglu, F.; Akkerman, V. Analytical study of an effect of gas compressibility on a burning accident in an obstructed passage. Phys. Fluids 2020, 32, 074602. [Google Scholar] [CrossRef]
  30. Akkerman, V.; Law, C.K. Effect of acoustic coupling on power-law flame acceleration in spherical confinement. Phys. Fluids 2013, 25, 013602. [Google Scholar] [CrossRef]
  31. Kuznetsov, M.; Ciccarelli, G.; Dorofeev, S.; Alekseev, V.; Yankin, Y.; Kim, T. DDT in methane-air mixtures. Shock Waves 2002, 12, 215–220. [Google Scholar] [CrossRef]
  32. Morley, C. GASEQ—A Chemical Equilibrium Program. 2005. Available online: http://www.gaseq.co.uk/ (accessed on 1 May 2020).
  33. Demir, S.; Sezer, H.; Bush, T.; Akkerman, V. Promotion and mitigation of premixed flame propagation in a gaseous-dusty environment with various dust distributions. Fire Saf. J. 2019, 105, 270–276. [Google Scholar] [CrossRef]
Figure 1. Illustration of flame acceleration inside an obstructed passage (a) and delayed burning between obstacles (b).
Figure 1. Illustration of flame acceleration inside an obstructed passage (a) and delayed burning between obstacles (b).
Fluids 05 00115 g001
Figure 2. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a 2D geometry for stoichiometric ( ϕ = 1 ) methane (CH4)-air and propane (C3H8)-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Figure 2. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a 2D geometry for stoichiometric ( ϕ = 1 ) methane (CH4)-air and propane (C3H8)-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Fluids 05 00115 g002
Figure 3. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a cylindrical-axisymmetric geometry for stoichiometric CH4-air and C3H8-air burning with various blockage ratios, α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Figure 3. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a cylindrical-axisymmetric geometry for stoichiometric CH4-air and C3H8-air burning with various blockage ratios, α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Fluids 05 00115 g003
Figure 4. Comparison of the 2D and cylindrical-axisymmetric geometries. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) for stoichiometric CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Figure 4. Comparison of the 2D and cylindrical-axisymmetric geometries. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) for stoichiometric CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Fluids 05 00115 g004
Figure 5. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometry for stoichiometric CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 and various power factors n = 1.33 ,   1.4 ,   1.5 .
Figure 5. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometry for stoichiometric CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 and various power factors n = 1.33 ,   1.4 ,   1.5 .
Fluids 05 00115 g005
Figure 6. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a 2D geometry for lean ( ϕ = 0.8 ), stoichiometric ( ϕ = 1 ), and rich ( ϕ = 1.2 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Figure 6. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a 2D geometry for lean ( ϕ = 0.8 ), stoichiometric ( ϕ = 1 ), and rich ( ϕ = 1.2 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Fluids 05 00115 g006
Figure 7. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a 2D geometry for highly-lean ( ϕ = 0.6 ) and highly-rich ( ϕ = 1.4 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Figure 7. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a 2D geometry for highly-lean ( ϕ = 0.6 ) and highly-rich ( ϕ = 1.4 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Fluids 05 00115 g007
Figure 8. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a cylindrical-axisymmetric geometry for lean ( ϕ = 0.8 ), stoichiometric ( ϕ = 1 ), and rich ( ϕ = 1.2 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Figure 8. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a cylindrical-axisymmetric geometry for lean ( ϕ = 0.8 ), stoichiometric ( ϕ = 1 ), and rich ( ϕ = 1.2 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Fluids 05 00115 g008
Figure 9. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a cylindrical-axisymmetric geometry for highly-lean ( ϕ = 0.6 ) and highly-rich ( ϕ = 1.4 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Figure 9. Time evolution of the flame tip position Z t i p (a) and velocity U t i p (b) in a cylindrical-axisymmetric geometry for highly-lean ( ϕ = 0.6 ) and highly-rich ( ϕ = 1.4 ) CH4-air burning with various blockage ratios α = 0 , 1 / 3 , 1 / 2 , 2 / 3 .
Fluids 05 00115 g009
Figure 10. The flame run-up distance versus the equivalence ratio ϕ for CH4-air and C3H8-air burning at various blockage ratios: α = 0 , 1 / 3 , 1 / 2 , 2 / 3 in a 2D (a) and cylindrical-axisymmetric (b) geometries.
Figure 10. The flame run-up distance versus the equivalence ratio ϕ for CH4-air and C3H8-air burning at various blockage ratios: α = 0 , 1 / 3 , 1 / 2 , 2 / 3 in a 2D (a) and cylindrical-axisymmetric (b) geometries.
Fluids 05 00115 g010
Figure 11. Comparison of our theoretical formulation with the experiments [31] and computational simulations [6] (FAST, yellow, and ALLA, blue). The flame tip velocity vs. its position for the blockage ratios α = 0.3 (a) and 0.6 (b).
Figure 11. Comparison of our theoretical formulation with the experiments [31] and computational simulations [6] (FAST, yellow, and ALLA, blue). The flame tip velocity vs. its position for the blockage ratios α = 0.3 (a) and 0.6 (b).
Fluids 05 00115 g011
Figure 12. Comparison of our theoretical formulation with the experiments [31] and computational simulations [6] (FAST, yellow, and ALLA, blue). The flame tip velocity vs. its position for the blockage ratios α = 0.3 (a) and 0.6 (b).
Figure 12. Comparison of our theoretical formulation with the experiments [31] and computational simulations [6] (FAST, yellow, and ALLA, blue). The flame tip velocity vs. its position for the blockage ratios α = 0.3 (a) and 0.6 (b).
Fluids 05 00115 g012
Figure 13. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 , without and with dust particles (combustible, inert, and combined) of particle radius r s = 75   μ m and concentration c s = 50   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Figure 13. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 , without and with dust particles (combustible, inert, and combined) of particle radius r s = 75   μ m and concentration c s = 50   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Fluids 05 00115 g013
Figure 14. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 , without and with dust particles (combustible, inert, and combined) of particle radius r s = 75   μ m and concentration c s = 120   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Figure 14. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 , without and with dust particles (combustible, inert, and combined) of particle radius r s = 75   μ m and concentration c s = 120   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Fluids 05 00115 g014
Figure 15. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 , without and with dust particles (combustible, inert, and combined) of particle radius r s = 75   μ m and concentration c s = 250   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Figure 15. Time evolution of the flame tip position Z t i p (a,c) and velocity U t i p (b,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 , without and with dust particles (combustible, inert, and combined) of particle radius r s = 75   μ m and concentration c s = 250   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Fluids 05 00115 g015
Figure 16. Time evolution of the flame tip position Z t i p (a,b) and velocity U t i p (c,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 with and without dust particles (combustible, inert, and combined) of particle radius r s = 10   μ m and concentration c s = 120   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Figure 16. Time evolution of the flame tip position Z t i p (a,b) and velocity U t i p (c,d) in a 2D (a,b) and cylindrical-axisymmetric (c,d) geometries for lean CH4-air burning of ϕ = 0.7 with and without dust particles (combustible, inert, and combined) of particle radius r s = 10   μ m and concentration c s = 120   g / m 3 , for various blockage ratios, α = 0 , 1 / 3 , 2 / 3 .
Fluids 05 00115 g016
Table 1. The parameters for methane-air and propane-air combustion: the thermal expansion ratio, Θ , the laminar burning velocity, U f , and the sound speed in the fuel mixture, c 0 , versus the equivalence ratio, ϕ .
Table 1. The parameters for methane-air and propane-air combustion: the thermal expansion ratio, Θ , the laminar burning velocity, U f , and the sound speed in the fuel mixture, c 0 , versus the equivalence ratio, ϕ .
Methane-Air Fuel Mixtures
ϕ 0.60.70.80.911.11.21.31.4
Θ ρ u / ρ b 5.546.116.657.127.487.557.437.287.09
U f   ( m / s )0.0890.1690.2540.3250.3710.3830.3450.2500.137
c 0   ( m / s )351.5352.1352.7353.3353.9354.5355.1355.6356.2
Propane-Air Fuel Mixtures
ϕ 0.630.70.80.911.11.21.31.4
Θ ρ u / ρ b 6.046.567.157.668.028.0887.887.74
U f   ( m / s )0.1470.2170.3030.3740.4180.4290.3990.3220.226
c 0   ( m / s )343342.3341.6340.8340.1339.5338.8338.1337.5

Share and Cite

MDPI and ACS Style

Kodakoglu, F.; Demir, S.; Valiev, D.; Akkerman, V. Analysis of Gaseous and Gaseous-Dusty, Premixed Flame Propagation in Obstructed Passages with Tightly Placed Obstacles. Fluids 2020, 5, 115. https://doi.org/10.3390/fluids5030115

AMA Style

Kodakoglu F, Demir S, Valiev D, Akkerman V. Analysis of Gaseous and Gaseous-Dusty, Premixed Flame Propagation in Obstructed Passages with Tightly Placed Obstacles. Fluids. 2020; 5(3):115. https://doi.org/10.3390/fluids5030115

Chicago/Turabian Style

Kodakoglu, Furkan, Sinan Demir, Damir Valiev, and V’yacheslav Akkerman. 2020. "Analysis of Gaseous and Gaseous-Dusty, Premixed Flame Propagation in Obstructed Passages with Tightly Placed Obstacles" Fluids 5, no. 3: 115. https://doi.org/10.3390/fluids5030115

APA Style

Kodakoglu, F., Demir, S., Valiev, D., & Akkerman, V. (2020). Analysis of Gaseous and Gaseous-Dusty, Premixed Flame Propagation in Obstructed Passages with Tightly Placed Obstacles. Fluids, 5(3), 115. https://doi.org/10.3390/fluids5030115

Article Metrics

Back to TopTop