Abstract
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 . 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 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 . Therefore, we assume. 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 and , respectively. Following the realistic assumptions above, the system of governing equations has the form: 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:
2.1. Nondimensional Analysis
In this study, we introduce the following dimensionless variables: We assume that at initial stage the temperature of gas and solid particles are the same, that is, . Therefore, using (2.6) in (2.1)β(2.5), the dimensionless system of governing equations has the following form: with the initial conditions: In (2.7)β(2.11), the following dimensionless parameters have been introduced: 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 and (see, e.g., [22, 23, 25] and the references therein). In the events that and parameters , , and 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 , we get Adding (3.1) and (2.10) together gives Rearranging, we get From (2.8) and (2.9), we obtain Substituting (3.4) into (3.3) yields Integrating (3.5) with respect to , we obtain where is constant of integration.
Applying the initial condition (2.11) in (3.6), we get Using (3.7) in (3.6) gives exact integral of the system of equations (2.7)β(2.10) as Now, by substituting in (3.8) into (2.7), it is possible to obtain the following equation: 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 , solid particle temperature , gas concentration and liquid-fuel droplet radius for typical practical rate of reactions with various values of heat exchange parameter when , , , , , , ,. The qualitative impacts of interphase heat exchange between gas and solid particles in this study are ascertained in the figures.
(a)
(b)
(c)
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 , 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 and , 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 , , , and with numeric exponent for various values of at ignition point. Important facts to note from these results are the following. (i)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. (ii)It is also obvious from Figuresββ5(a) and 5(b) that an increase in numeric exponent 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.
Table 1 shows the variation of dimensionless ignition times with interphase heat exchange parameter 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.
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: (1)occurrence of convectional thermal explosions and quasi-steady solution after explosion. (2)significant delay in ignition time. (3)reduction in inert solid particles temperature. (4)remarkable reduction in the size of the liquid-fuel droplet leading to complete evaporation. (5)Depletion of gas concentration during the combustion process.
The results highlighted above do guarantee safety in practical explosion and combustion devices.
Nomenclature
English Symbols: | Preexponential factor (1/s) |
: | Dimensionless parameter expressing ratio of latent heat of evaporation and specific combustion energy |
: | Molar concentration of combustible gaseous mixture (kmol/m3) |
: | Specific heat capacity of the gas phase at constant pressure (JKβ1βkgββ1) |
: | Specific heat capacity of the solid phase at constant pressure (JKβ1βkgββ1) |
: | Activation energy (J/kmol) |
: | Latent heat of evaporation (J/kg) |
: | Molar mass (kg/kmol) |
: | Numerical exponent |
: | Number of droplets per unit volume (mβ3) |
: | Number of solid particles per unit volume (mβ3) |
: | Specific combustion energy (J/kg) |
: | Droplet radius (m) |
: | Solid particle radius (m) |
: | Universal gas constant Jkmolβ1βKβ1 |
: | Time (s) |
: | Combustible gas initial temperature (K) |
: | Dusty particle initial temperature (K) |
: | Combustible gas temperature (K) |
: | Dusty particle temperature (K). |
: | 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 |
: | Dimensionless parameter expressing heat loss via convection from gas phase |
: | Dimensionless parameter expressing heat loss via radiation from gas phase |
: | 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. |
0: | Initial |
: | Gas mixture |
: | Combustible gas component of the mixture (fuel) |
: | Liquid droplets |
: | Constant pressure |
: | Ignition. |