A mathematical model for thermal explosion in a combustible dusty gas containing fuel droplets with general Arrhenius reaction-rate laws, convective and radiative heat losses, and interphase heat exchange between gas and inert solid particles is investigated. The objective of the study is to examine the effects of interphase heat exchange between the gas and solid particles on (i) ignition of reacting gas, (ii) accumulation of heat by the solid particles during combustion process (iii) evaporation of the liquid fuel droplets, and (iv) consumption of reacting gas concentration. The equations governing the physical model with realistic assumptions are stated and nondimensionalised leading to an intractable system of first-order coupled nonlinear differential equations, which is not amenable to exact methods of solution. Therefore, we present numerical solutions as well as different qualitative effects of varying interphase heat exchange parameter. Graphs and Table feature prominently to explain the results obtained.
1. Introduction
The concept of thermal explosion in combustible dusty gas containing fuel droplets is of great importance in safety aspect of nuclear facilities, furnaces, gas turbines and internal combustion engines, coal mine, and so forth. The essence of the inert solid particles in the combustible gas is to delay the ignition or explosion which could cause catastrophe. The phase of research on combustion and explosions that began at the end of the nineteenth century and continues to this day is associated with the invention of the internal combustion engine, with the development of explosive technology and of internal ballistics for artillery, and, in recent decade, with the extensive introduction of jet and diesel engines. In many respects, these have stimulated the rapid development of combustion science [1]. The procedure for thermal explosion in gases, which contain fuel droplets, has been of much interest. After Semenov [2] developed the basic theory of phenomenon of thermal explosion, models that are more complicated have been suggested in [3, 4].
The focus of this research work is therefore on the long-standing problem of thermal explosion and ignition in a combustible gas containing fuel droplets and its numerous applications to furnaces, gas turbines, and internal combustion engines [3–7]. Over recent years, the theoretical analysis of this problem has been performed mainly by the use of the geometrical asymptotic method of integral manifolds [8–17]. An alternative approach to the problem is to analyse the equations in some limiting cases. Klammer et al. [18] investigated ignition, combustion, and detonation processes in dusty gases with combustion reaction. The dusty gas was considered as a two-continuum medium taking into account transport effects in the phases and nonequilibrium chemical reactions. Two-dimensional problems of ignition and detonation were developed in a plane gallery caused by a supersonic inflow stream and heating of the closed end of the gallery are studied with the analytical method of catastrophe and two finite-difference numerical methods. Krainov and Shaurman [19] studied the limits of flame propagation in a gas with suspended inert particles in the presence of external heat removal. The mathematical model used was based on an unsteady heat-diffusion two-temperature model of gas combustion in the presence of inert particles. The problem was solved by a numerical method. A parametric analysis was performed, and critical values of the parameter that characterizes external heat removal were obtained. Dispersed-phase parameters were determined for which the two-temperature nature of the medium was insignificant. Thus, the critical conditions of flame propagation in the gas with inert particles were obtained. Comparison of the numerical and analytic solutions showed that they agreed satisfactorily with small particles of the dispersed phase. Ben-Dor and Igra [20] considered the relaxation zone behind normal shock waves in a reacting dusty gas. It was assumed that the gas is monatomic. The conservation equations for a suspension composed of an ionized gas and small solid dust particles were formulated and then solved numerically. The solution revealed that the presence of the dust has a significant effect on the postshock flow field. Because of the dust, the relaxation zone was longer than in the pure plasma case; the equilibrium values for the suspension pressure and density were higher than in the dust-free case, whereas the values obtained for the temperature, degree of ionization, and velocity was lower. The numerical solution was executed for shock Mach numbers ranging from 10 to 17. It was found that the thermal relaxation length for the plasma decreases rapidly with increasing shock Mach number, whereas the thermal relaxation length for the suspension increases slightly with increasing M. Gol’dshtein et al. [8] studied criteria for thermal explosion with reactant consumption in a dusty gas. The dynamical regimes of the system were classified as slow regimes, thermal explosion with delay, and thermal explosion (without delay). The critical transition conditions for the different dynamical regimes were analysed. They emphasized that the critical conditions for transition between slow regimes and explosion with delay was a thermal explosion limit. The Thermal explosion limit was described in the phase space by a so-called duck-trajectory. El-Sayed [21] investigated the critical conditions of the adiabatic explosion problem of a gas-solid (dusty gas) mixture. The definitions used for the homogeneous gas to determine criticality were used for the gas-solid mixture. The analysis revealed that the classical definition of the critical point can be adopted and modified to determine the critical condition in γ-η and γ-θ domains. It was also found that using the definition of criticality as an inflection point in the critical trajectory in the θ-η plane gives the same results as given by the classical definition of criticality. It was interesting to see that the critical and maximum points in the gas temperature-concentration domain could coincide. It was found that the presence of a solid produces more than one critical temperature. The limiting cases of the problem were also offered. The numerical solution showed that the supercritical trajectory shows a thermal runaway for the gas over the solid at the end of reaction.
In this paper, an attempt has been made to extend the problem of thermal explosion in a combustible gas mixture with addition of inert solid particles and to generalise the problem based on temperature dependence of the reaction rate (i.e., Arrhenius Power-law model equation [22–24]) given as K(Tg)=A0(Tg/Tg0)nexp(-E/RuTg) while taking into account convective and radiative heat losses, temperature dependence of density, and thermal conductivity of the gas. Therefore, a generalised physical model for thermal explosion in combustible dusty gas mixture containing fuel droplets is developed in the present paper. The main interest is focused on numerical solutions of the system of coupled nonlinear ordinary differential equations governing the physical model. Parametric analysis is performed, and the results show that the delay before ignition of the reacting gas, accumulation of heat by dusty particles, evaporation of liquid fuel droplets, and consumption of reacting gas concentration depend significantly on interphase heat exchange between the gas and solid particles.
2. Mathematical Model
In this study, we consider combustible dusty gas with addition of fuel droplets. Dusty gas is a combustible gas with addition of solid particles. It is assumed that the solid phase is inert, monosized, and uniformly heated, and the dusty gas is optically thick. The dependence of the convective heat transfer coefficient on both gas temperature and droplet radii is taken into account. Thermal conductivity of both monoatomic and polyatomic gases is proportional to Tg. Therefore, we assumeλg(Tg)=λg0Tg/Tg0. The rate of reaction is based on generalised temperature-dependent Arrhenius’s equation. Temperature dependence of the gas density and preexponential factor is given as ρg(Tg)=ρg0Tg0/Tg and A(Tg)=A0(Tg/Tg0)n, respectively. Following the realistic assumptions above, the system of governing equations has the form:Cpgϕgρg0Tg0Tg-1dTgdt=QfmfϕgcfA0Tg0-nTgnexp(-ERuTg)-4πRdndλg0TgTg0(Tg-Tg0)-4πRd2ndσ(Tg4-Tg04)-4πRsnsλs(Tg-Ts),ρsCpsdTsdt=4πRSnsλs(Tg-Ts),dRddt=-1Lρdλg0RdTgTg0(Tg-Tg0)-1Lρdσ(Tg4-Tg04),φgdcfdt=-φgcfA0Tg0-nTgnexp(-ERuTg)+4πRdndλg0LmfTgTg0(Tg-Tg0)+4πRd2ndσLmf(Tg4-Tg04).
The system (2.1)–(2.4) includes energy equations for the reacting gas (2.1) and solid particles (2.2), a mass equation for the liquid droplet (2.3) and a concentration equation for the reacting gas mixture (2.4). Initial conditions for the system:
Tg(0)=Tg0,Ts(0)=Ts0,Rd(0)=Rd0,cf(0)=cf0.
2.1. Nondimensional Analysis
In this study, we introduce the following dimensionless variables:θg=(Tg-Tg0)ERuTg02,θs=(Tg-Ts)ERuTs02,r=RdRd0,η=cfcf0,τ=tt*.
We assume that at initial stage the temperature of gas and solid particles are the same, that is, Tg0=Ts0. Therefore, using (2.6) in (2.1)–(2.5), the dimensionless system of governing equations has the following form:(1+βθg)-1dθgdτ=δη(1+βθg)nexp(θg1+βθg)-r{α1θg(1+βθg)+α2r[(1+βθd)4-1]}-α3(θg-θs),dθsdτ=α3μ(θg-θs),drdτ=-μ*r{α1θg(1+βθg)+α2r[(1+βθg)4-1]},dηdτ=-α4δη(1+βθg)nexp(θg1+βθg)+ψr{α1θg(1+βθg)+α2r[(1+βθg)4-1]},
with the initial conditions:θg(0)=0,θs(0)=0,r(0)=1,η(0)=1.
In (2.7)–(2.11), the following dimensionless parameters have been introduced:β=RuTg0E,B=LQf,δ=EQfmfcf0A0t*ρg0CpgRuTg02exp(-ERuTg0),α1=4πRd0ndλg0t*ρg0Cpgϕg,α2=4πRd02ndσTg02Et*ρg0CpgϕgRu,α3=4πRsnsλst*ρg0Cpgϕg,α4=Bψ,μ*=ρg0CpgϕgRuTg024πRd03ndρdEL,μ=ρg0CpgϕgρsCps,ψ=ρg0CpgϕgRuTg0Emfcf0L.
Some special cases of the system of coupled nonlinear ordinary differential equations (2.7)–(2.10) with initial conditions (2.11) and related problems have been studied for n≠0 and α3=0 (see, e.g., [22, 23, 25] and the references therein). In the events that n=0 and parameters α1, α2, and α3 are varied, analyses have been performed both analytically and numerically in ([2, 8–17, 26, 27], and references cited therein). In this study, attention is focused on the generalization of the problem and effect of the interphase heat exchange between the reacting gas and solid particles on the ignition of reacting gas, accumulation of heat by dusty particles, rate of evaporation of liquid fuel droplets, and depletion of reacting gas concentration.
3. Numerical Computation
In this section, we present numerical solutions of the system of governing equations (2.7)–(2.10) with initial conditions (2.11). Meanwhile, before actual computation, we begin by deriving the exact integral for the system of ordinary differential equations as follows:
Multiplying (2.7) by α4, we getα4(1+βθg)-1dθgdτ=α4δη(1+βθg)nexp(θg1+βθg)-α4r{α1θg(1+βθg)+α2r[(1+βθd)4-1]}-α3α4(θg-θs).
Adding (3.1) and (2.10) together givesα4(1+βθg)-1dθgdτ+dηdτ=-α4r{α1θg(1+βθg)+α2r[(1+βθd)4-1]}+ψr{α1θg(1+βθg)+α2r[(1+βθd)4-1]}-α3α4(θg-θs).
Rearranging, we getα4(1+βθg)-1dθgdτ+dηdτ=(ψ-α4)r{α1θg(1+βθg)+α2r[(1+βθd)4-1]}-α3α4(θg-θs).
From (2.8) and (2.9), we obtainα3(θg-θs)=1μdθsdτ{α1θg(1+βθg)+α2r[(1+βθg)4-1]}=-rμ*drdτ.
Substituting (3.4) into (3.3) yieldsα4(1+βθg)-1dθgdτ+dηdτ=(α4-ψμ*)r2drdτ-α4μdθsdτ.
Integrating (3.5) with respect to τ, we obtainα4βln(1+βθg)+η=(α4-ψ3μ*)r3-α4μθs+N,
where N is constant of integration.
Applying the initial condition (2.11) in (3.6), we getN=1-ψ3μ*(B-1),sinceα4=Bψ.
Using (3.7) in (3.6) gives exact integral of the system of equations (2.7)–(2.10) asη=η(θg,θs,r)=1-ψ{Bβln(1+βθg)-(1-B3μ*)(1-r3)+Bμθs}.
Now, by substituting η in (3.8) into (2.7), it is possible to obtain the following equation:dθgdτ=δ{1-ψ[Bβln(1+βθg)-(1-B3μ*)(1-r3)+Bμθs]}(1+βθg)n+1exp(θg1+βθg)-r(1+βθg){α1θg(1+βθg)+α2r[(1+βθd)4-1]}-α3(1+βθg)(θg-θs).
The resulting systems of (2.8), (2.9), and (3.9) are coupled nonlinear differential equations and do not possess closed form solutions. Therefore, a numerical technique based on finite difference approximation is used to solve the problems. The effect of the interphase heat exchange parameter is adequately examined on the ignition time, combustible gas temperature, solid particle temperature, droplet radius, and combustible gas concentration as shown in the figures.
4. Results
The general characteristic features of Figures 1–4 show time histories of dimensionless gas temperature θg, solid particle temperature θs, gas concentration η and liquid-fuel droplet radius r for typical practical rate of reactions with various values of heat exchange parameter α3 when δ=0.3, β=0.05, B=0.008, μ*=5, ψ=0.3, α1=3E-04, α2=2E-05,μ=1. The qualitative impacts of interphase heat exchange between gas and solid particles α3 in this study are ascertained in the figures.
Variation of dimensionless combustible gas temperature with interphase heat exchange parameters when n={-2,0,0.5}.
Variation of dimensionless solid particle temperature with interphase heat exchange parameters when n={-2,0,0.5}.
Variation of dimensionless droplet radius with interphase heat exchange parameters when n={-2,0,0.5}.
Variation of dimensionless combustible gas concentration with interphase heat exchange parameters when n={-2,0,0.5}.
We notice in Figure 1 the occurrence of convectional thermal explosions in sensitized, Arrhenius, and bimolecular rate of reactions, respectively. Meanwhile in bimolecular rate of reaction, when α3=1.0, it is quite interesting to note that the problem has a quasi-steady solution after explosion. In this case, the medium is no longer homogeneous and diffusion term has to be taken into consideration. The medium could not be regarded as well stirred. Of a particular interest is the notable delay in the ignition time as the interphase heat exchange between gas and solid particles increases.
In Figure 2, it is observed that the heat exchange between the gas and solid particles raises the temperature of the solid particles at different time for sensitized, Arrhenius, and bimolecular rate of reactions. This implies that the solid particles accumulate heat over a period of time from reacting gas, since the solid particles are inert.
Figure 3 shows that the size of the liquid-fuel droplet radius significantly reduces as heat exchange between the gas and solid particles increases for sensitized, Arrhenius, and bimolecular rate of reactions. This reveals that the rate of evaporation of the liquid-fuel droplet substantially depends on heat exchange between the reacting gas and solid particles. We note that when α3=0.5 and α3=1.0, the liquid-fuel droplet completely evaporated just before ignition.
Figure 4 reveals a remarkable reduction in gas concentration as heat exchange between the gas and solid particles increases for sensitized, Arrhenius, and bimolecular rate of reactions. This depicts depletion of gas concentration during the combustion process.
In Figure 5, we present variation of θg, θs, r, and η with numeric exponent for various values of α3 at ignition point. Important facts to note from these results are the following.
As can be seen from Figures 5(a), 5(c), and 5(d), an increase in interphase heat exchange between the gas and solid particles causes significant reduction in the values of gas temperature, droplet radius, and gas concentration whereas in Figure 5(b) it raises the value of solid particle temperature at ignition point for sensitized, Arrhenius, and bimolecular reaction rates, respectively.
It is also obvious from Figures 5(a) and 5(b) that an increase in numeric exponent n ascertains reduction in the values of gas and solid particles temperature at ignition point, whereas in Figures 5(d) and 5(c) we observe increment in the values of gas concentration and insignificant reduction in the values of droplet radius at ignition point, respectively.
Variation of θg, θs, r, and η with numeric exponent for various values of α3 at ignition point.
Table 1 shows the variation of dimensionless ignition times τig with interphase heat exchange parameter α3 for typical sensitized, Arrhenius, and bimolecular rate of reactions at ignition point. It is significant to note from the table that increase in heat exchange between the gas and solid particles causes a notable delay in the ignition time of the combustible gas mixture for sensitized, Arrhenius, and bimolecular rate of reactions. The results also reveal that responsiveness to delay in ignition time is characterized by numeric exponent, which means that delay is more enunciate in sensitized reaction compared to Arrhenius and bimolecular reactions, respectively.
Variation of ignition times with interphase heat exchange parameter at ignition point whenn={-2,0,0.5}.
n
τig
α3=0
α3=0.1
α3=0.5
α3=1.0
−2.0
4.567872
4.938240
5.296256
7.160448
0.0
4.074048
4.444416
5.555520
6.296256
0.5
3.950592
4.320960
5.432064
6.049344
5. Conclusion
The problem of thermal explosion in combustible dusty gas mixtures containing fuel droplets have been extended to permit a more general temperature dependent rate of reaction for most typical practical reactions under physically reasonable assumptions. The mathematical formulation involves a system of four highly nonlinear ordinary differential equations, which was solved numerically by finite difference approximation scheme. In multiple phase processes, interphase heat exchange plays the role of heat losses in homogeneous combustible gas mixtures. The model with heat losses provides more detailed information about the effects of interphase heat exchange and numeric exponent characterising sensitized, Arrhenius, and bimolecular rate of reactions. The study revealed that interphase heat exchange between gas and solid particles, and numeric exponent play important roles in determining the following:
occurrence of convectional thermal explosions and quasi-steady solution after explosion.
significant delay in ignition time.
reduction in inert solid particles temperature.
remarkable reduction in the size of the liquid-fuel droplet leading to complete evaporation.
Depletion of gas concentration during the combustion process.
The results highlighted above do guarantee safety in practical explosion and combustion devices.
NomenclatureEnglish SymbolsA:
Preexponential factor (1/s)
B:
Dimensionless parameter expressing ratio of latent heat of evaporation and specific combustion energy
cf:
Molar concentration of combustible gaseous mixture (kmol/m3)
Cpg:
Specific heat capacity of the gas phase at constant pressure (JK−1 kg −1)
Cps:
Specific heat capacity of the solid phase at constant pressure (JK−1 kg −1)
E:
Activation energy (J/kmol)
L:
Latent heat of evaporation (J/kg)
mf:
Molar mass (kg/kmol)
n:
Numerical exponent
nd:
Number of droplets per unit volume (m−3)
ns:
Number of solid particles per unit volume (m−3)
Q:
Specific combustion energy (J/kg)
Rd:
Droplet radius (m)
Rs:
Solid particle radius (m)
Ru:
Universal gas constant Jkmol−1 K−1
t:
Time (s)
Tg0:
Combustible gas initial temperature (K)
Ts0:
Dusty particle initial temperature (K)
Tg:
Combustible gas temperature (K)
Ts:
Dusty particle temperature (K).
Greek Symbolsρ:
Density of the combustible gaseous mixture (kg/m3)
λ:
Thermal conductivity (Wm−1 K−1)
φ:
Volumetric phase content (dimensionless)
σ:
Stefan-Boltzmann’s constant
β:
Dimensionless activation energy
τ:
Dimensionless time
δ:
Dimensionless parameter expressing reciprocal of the characteristic time for adiabatic temperature rise
α1:
Dimensionless parameter expressing heat loss via convection from gas phase
α2:
Dimensionless parameter expressing heat loss via radiation from gas phase
α3:
Dimensionless parameter expressing interphase heat exchange between gas and solid particle
μ:
Dimensionless parameter expressing energy needed to transfer heat from gas phase to solid phase
μ*:
Dimensionless parameter expressing energy needed to evaporate all fuel droplets
ψ:
Dimensionless parameter expressing energy needed to consume all gas concentration.
Subscripts0:
Initial
g:
Gas mixture
f:
Combustible gas component of the mixture (fuel)
d:
Liquid droplets
p:
Constant pressure
ig:
Ignition.
ZeldovichYa. B.BarenblattG. I.LibrovichV. B.MakhviladzeG. M.1985New York, NY, USAConsultants Bureauxxi+597781350SemenovN. N.19581&2New York, NY, USAPergamon PressKuoK. K.1986New York, NY, USAJohn Wiley & SonsGrayJ. J.2000SecondBoston, Mass, USABirkhäuserxx+3381751835HeywoodJ. B.1988InternationalSingaporeMcGraw-HillStoneR.1992New York, NY, USAMacmillanWarnatzJ.MassU.DibbleR. W.20013rdBerlin, GermanySpringerGol'dshteinV.vladimir@bgumail.bgu.ac.ilZlnovievA.annzi@ivory.bgu.ac.ilSobolevV.ShchepakinaE.Criterion for thermal explosion with reactant consumption in a dusty gas199645219522103211910.1098/rspa.1996.0111ZBL0874.34052GoldfarbI.SazhinS.ZinovievA.Delayed thermal explosion in flammable gas containing fuel droplets: asymptotic analysis2004504399414211676110.1007/s10665-004-9388-5ZBL1178.76377GoldfarbI.Gol’dshteinV.ZinovievA.Delayed thermal explosion in porous media: method of invariant manifolds2002673263280190150510.1093/imamat/67.3.263ZBL1008.80004GoldfarbI.goldfarb@cs.bgu.ac.ilZinovievA.annaz@nace.ac.ilA study of delayed spontaneous insulation fires2003311649150010.1016/S0375-9601(03)00506-1GoldfarbI.Gol'dshteinV.GreenbergJ. B.KuzmenkoG.Thermal explosion in a droplet-gas cloud20004328931610.1088/1364-7830/4/3/304ZBL0995.76096GoldfardI.Gol’dshteinV.KuzmenkoG.SazhinS .S.Thermal radiation effect on thermal explosion in a gas containing fuel droplets1999376978710.1088/1364-7830/3/4/309Gol’dshteinV.SobolevV.Integral manifolds in chemical kinetics and combustion1992Providence, RI, USAAmerican Mathematical Society7392Gol’dshteinV.GoldfardI.ShreiberI.ZinovievA.Oscillations in a combustible gas bubble1998211710.1080/713665366ZBL1040.80501GoldfarbI.Gol’dshteinV.KatzD.SazhinS.Radiation effect on thermal explosion in gas containing evaporated fuel droplets20074635837010.1016/j.ijthermalsci.2006.06.014Gol’dshteinV.ZinovievA.Thermal explosion in multiple phase media199730847714780172534010.1016/S0362-546X(96)00356-2ZBL0896.34062KlammerJ.KlemensR.KorobeinikovV. P.MarkovV. V.Men'shovI. S.WolanskiP.On ignition and unsteady flows of dusty gases with combustion reactions19991421819010.1080/00102209908924189KrainovA. Yu.ShaurmanV. A.On the limits of flame propagation in a dusty gas199733440340810.1007/BF02671832Ben-DorG.IgraO.The relaxation zone behind normal shock waves in a reacting dusty gas. I Monatomic gases198227337739510.1017/S0022377800010965El-SayedS. A.shamad53@hotmail.comAdiabatic thermal explosion of a gas-solid mixture2004176223725610.1080/00102200490256144AdegbieK. S.2008Ogbomoso, NigeriaLautechAdegbieK. S.On delay-type behaviour of thermal explosion in a combustible gas mixture containing fuel droplets with arrhenius power-law model2008686982International Union of Pure and Applied Chemistry (IUPAC) Goldbook Definition of Modified Arrhenius Equation, http://en.wikipedia.org/wiki/Arrhenius/equationAdegbieK. S.Some remarks on Arrhenius Power-law Model for thermal explosion in a combustible gas mixture containing fuel droplets2008101358370AyeniR. O.OtolorinO.AyodeleO. T.AdegbieK. S.Effect of frank-kamenetskii parameter on thermal explosion of flammable gas containing fuel droplets2005915562Frank-KamenetskiiD. A.1969New York, NY, USAPlenum Press