SciELO - Scientific Electronic Library Online

vol.40 número4Esterification and transesterification assisted by microwaves of crude palm oil: Heterogeneous catalysisKinetic modeling of a photocatalytic reactor designed for removal of gas-phase benzene: a study on limiting resistances using design of experiments índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados



  • No hay articulos citadosCitado por SciELO

Links relacionados

  • No hay articulos similaresSimilares en SciELO


Latin American applied research

versión impresa ISSN 0327-0793

Lat. Am. appl. res. vol.40 no.4 Bahía Blanca oct. 2010


An analytical study of radiation effect on the ignition of magnesium particles using perturbation theory

M. Bidabadi, S. Montazerinejad and S.A. Fanaee

Iran University of Science and Technology
Department of Mechanical Engineering, Combustion Research Laboratory, Narmak, Tehran 16844, Iran,,

Abstract — An analytical method is proposed to solve the heat Energy equation of combustion of magnesium particles with allowance for the heterogeneous chemical reactions and the region of the thermal influence of the particle on the gas. In this model the solution of the problem in a steady formulation is found and the flame propagation mechanism is considered to be radiation and conduction.Radiative heat transfer plays a major role in single particle combustion. Flame equations of the single particle combustion are solved via the new nonlinear differential equation by using perturbation theorem. In this paper we can apportion the magnesium-particle ignition to regular regimes and also regime of particle extinction and the ignition. Under the steady approximation, the values of particle and gas phase temperatures are calculated and by adding radiation term in conservative equations, both of these parameters will be increased.

Keywords — Heterogeneous. Radiative Heat Transfer. Nonlinear Equation. Perturbation Theorem.


The problem of physicomathematical modeling of ignition and combustion of metal samples is of considerable interest for various branches of industry (Fedorovet al., 2003).

The main objects described in Fedorov et al. (2003) are the pointwise and partly distributed models of ignition of small metal particles with low-temperature oxidation proceeding on the particle surface. The heat dissipated in to the gas phase was ignored. This means that the thickness of the so called surrounding film is negligibly small. It seems of interest to study the effect of this parameter on the thermal history of the reacting particle.

Particle diameter plays a significant role in determining the relevant combustion mechanisms by affecting the characteristic transport diffusion time relative to the chemical kinetics time. A large diameter particle at high pressure may burn under diffusion controlled conditions, whereas a small particle at low pressure may burn under kinetically controlled conditions (Yetter and Dryer, 2001; Huang et al., 2006). The oxidizer type also has strong effects on magnesium particle combustion, since the flame and surface temperatures can be affected by transport in the surrounding gas. In the last two decades with the rapid development of nonlinear science, an ever increasing interest of scientists and engineers in the analytical techniques for nonlinear problems appeared. The widely applied techniques are perturbation methods. Perturbation method is one of the well-known methods to solve the nonlinear equations which was studied by a large number of researchers such as Bellman (1964), Cole (1968) and O'Malley (1974). A valuable model of ignition of single magnesium particle, proposed by Fedorov and Shul'gin (2006), is expanded in this study. In previous theoretical models (for example, Fedorov et al., 2003; Fedorov and Shul'gin, 2006; Fedorov, 1994; Fedorov, 1996) it is assumed that radiative heat transfer was neglected. The radiation term is added to the previous model (Fedorov and Shul'gin, 2006), and its effect on the combustion of magnesium particle is investigated. The initial investigation of the combustion of bimodal aluminum and iron particles is also studied with help of the presented model.


Combustion of fine metal particles, which takes place in the gas phase (where evaporated sample particle reacts with oxidizer) is formulated here.

We consider a metal particle of radius rp surrounded by a gas layer of thickness L-rp. We assume an exothermic chemical reaction of oxidation to proceed on the sample surface.

The main focus of this research is made on the effect of radiation on the combustion of metal dust particles. This paper develops the previous model published by Fedorov and Shul'gin (2006).

Then the mathematical model that describes the temperature fields in general case Ti has the form


The temperature of the ambient gas and sample particle are indicated by Ti.

In the previous relationship, υ is the factor of symmetry equal to 0, 1, and 2 for the planar, cylindrical, and spherical cases, respectively, ρi, λi and c i are the density, thermal conductivity , and specific heat of the phase; the subscript i=1 and 2 refers to parameters of the gas and the particle respectively. Qrad is the radiation heat transfer, which was not taken into account in previous models but is considered in the present article. In fact, it distinguishes the present study from the previous works in this field.

Radiation heat transfer is given by the Stefan-Boltzmann law,


where σ and ei are the Stefan-Boltzmann constant and emissivity of the flame, respectively. In addition is the temperature at the boundary of the gas region.

In the present study, the steady state problem is investigated and therefore the governing equations for gas and particle are represented as follows:




The three needed boundary conditions for Eqs. (3) and (4) are in the core of particle, the interface between particle and gas oxidizer and the interface between the oxidizer and the surrounding environment.

The boundary conditions for the Eq. 3 and Eq. 4 are as follows:


where q0 is the heat release per unit mass of the oxide, ρ3 is the oxide density, h is the oxide-film thickness, dh/dt is the rate of variation of the oxide-film thickness.

Moreover, K represents the frequency factor, E the activation energy of low-temperature oxidation, R is the universal gas constant, the heat-transfer coefficient is calculated as follows (Fedorov et al., 2003),

After some simplifications, this coefficient is calculated by α=Nu/2, where the Nusselt number is taken as a constant equal to 2 (Fedorov and Shul'gin, 2006). In this research, the radiation effect is considered from the gas and particle to the surrounding environment with the temperature of .

Since the ambient temperature is lower than the gas and particle temperatures, the radiative heat transfer between the oxidizer and particle is neglected compared to the radiative heat transfer from the gas and particle to the surrounding environment.


We can write particle and gas phase energy conservation Eq. 3 and Eq. 4 in a non-dimensional form by using the following non-dimensional parameters.


Since the energy equation in each system (i.e., gas and particle) is conserved, the total energy of this system (gas+ particles) in the preheated zone of the gas-particles mixture is conserved.

As stated in Fedorov and Shul'gin (2006) the parameters have the following values:

The value of these parameters for the other particles (Al, Fe) are obtained from references.

By substituting dimensionless parameters in the governing equations, we can write the equations obtained for each phase.

Particle phase:


Gas phase:


The values of the non-dimensional emissivity (ε) for particle and gas are defined as:



Exact solutions are rare in many branches of fluid mechanics, solid mechanics and other physical phenomena because of nonlinearities, non-homogeneities, etc. Hence in engineering and physics we are forced to determine approximate solutions of the problem they are facing. These approximations may be purely numerical, purely analytical, or combination of numerical techniques.

In this work, we use the perturbation theory to determine approximate solutions of differential equations that govern the phenomenon. Perturbation techniques (Nayfeh, 1973; Naives, 1979) are widely applied for obtaining approximate solutions to these equations involving a small parameter ε. According to these techniques, the solution is represented by the first few terms of an asymptotic expansion, usually not more than two terms. The expansions may be carried out in terms of a parameter (small or large) which appears naturally in the equations, or which may be artificially introduced for convenience. In this article, the system of Eq. 3 and Eq. 4 for the particle and the gas is nonlinear and heterogeneous. By using perturbation series we need to expand our differential equation in power series of a small parameter ε.


Considering the accuracy that is needed for this work, the final solution of Eq. 3 and Eq. 4 is considered up to first order

Zero order

Particle phase:


Gas phase:


The boundary conditions for Eqs. (12) and (13) are obtained from the boundary condition (5). The boundary condition (14) for zero-order equation is determined by the perturbation theory. The leading term of two phase temperatures , satisfies the boundary conditions


The concepts of zero-order and first order are exactly related to perturbation derivation which has been used for solution of nonlinear equations.

The physical meaning of these boundary conditions is described as below:

  • The symmetric condition is assumed at = 0.
  • Two boundary conditions are needed at = 1. The first one is the matching condition for the temperature profile and the next one is the equality of heat flow from the particle and heat flow from the ignition of the particles.
  • At = L,the external heat flow from the oxidizer is equal to the convective heat transfer from the surrounding gas. According to Fedorov et al. (1996) and based on the present radiative heat transfer, the heat release due to the oxidation can be neglected with a good approximation.

First - order


The boundary conditions for Eq. 15 and Eq.16 are as follows:


Since the final solution is continuous for both particle and oxidizer and the ignition process is assumed homogeneous, the temperature equilibrium of gas and particles seems to be a logical assumption.

Since it is considered that the surrounding gas is heated by the output source (such as electrical power), the surrounding gas temperature remains constant. It is a reasonable assumption which is supposed by lots of researchers in the experimental apparatus, so = cte.

Under the steady approximation, the heat transfer equation in each phase is solved by matching the heat fluxes on the boundary of solid particle phase and gas phase; and the temperature distribution is consequently obtained.

The equations that describe the zero order of temperatures of the solid particle and the gas have solutions of the form

where υ= 0, 2 and ai and bi (i = 1, 2) are constants. Note that the case of cylindrical symmetry is considered in a similar manner; a logarithmic function is used as the basic solution. It follows from the condition at the point of symmetry that =const. Using the conjugation of temperatures and heat inflow at the interface between the particle and the gas (Eq. 5), and also the condition of heat transfer at the outer boundary and after some transformations, we find a function that describes the catastrophe/ignition manifold in the form


This expression determines the leading term of particle temperature; and after substituting the defined temperature into the related boundary conditions we find the coefficients a1 and b1, i.e., the gas temperature in the region surrounding the particle, in the form


Using the methods of the elementary catastrophe theory, we can obtain the turning points of the ignition manifold from the equation which has two roots:

Then, in accordance with (Fedorov, 1994; Fedorov, 1996), such a manifold can be constructed in the plane which is schematically shown in (Fedorov and Shul'gin, 2006), where are the values calculated at the turning points, i.e., for

Rigorously speaking, it is seen from (Fedorov and Shul'gin, 2006) that the steady solution is not unique. On the lower branch I with ββ_, the solution describes the regime of regular heating. For β+ββ_, we have ; it should be noticed that branch II is unstable (Fedorov and Shul'gin, 2006). Finally, for 0 ≥ ββ_ we have and this solution belongs to branch III, which corresponds to the phenomenological criterion of ignition. Based on these facts, we can formulate the following statement. The quantity β, which defines the dynamics of particle heating (regular or with ignition) is a decreasing function tending to a certain constant as

This is in agreement with an actually observed fact that injection of an additional mass of the gas into the system increases the ignition limit of the sample.

This assumption is illustrated in Fig. 1 and Fig. 2, shows the dependency of βq0 on for Nu =2 for the spherical and planar cases respectively. It can also be noted that the ignition limit is higher in the spherical case, because the criterion of regular development of the process 0≥ββ_ can be more readily satisfied. Indeed,

Figure 1. Dependency of on L for spherical case

Figure 2. Dependency of on L for planar case

The next important step is the computation of the first order temperature of particle and gas phase. By superposing the leading term of two phase temperatures into Eq. 15 and Eq. 16 the system of differential equations in each phase can be solved.

By solving the system of differential equations in each phase, the first order of particle and gas temperature in the cases of planar and spherical coordinates can be found.

Particle Temperature



Gas Temperature

Planar coordinate


Spherical coordinate


From the procedure of matching boundary conditions obtained from solutions on both sides of the boundary ( = 1) and also from system of boundary conditions (17) at = 0 and = , the parameters C2, C3, C4, C5, C6 for Eq. 20 and Eq. 21 are found.


In order to find the influence of the oxide-film thickness on the ignition of magnesium particle, some computations with different sizes of were performed.

The results plotted in Fig. 3 show that the main changes in temperature of the system occur mainly near the contact boundary, with r < 4. In this region, the temperature gradients acquire the highest values. If the region occupied by the gas is large, the temperature changes only weakly. This allows us to argue that the boundary of the region occupied by the gas can be assumed to be approximately equal to four particle radii.

Increasing the values of gas thickness () gives rise to the gas temperature gradient; This phenomenon is understood considering a layer of insulation which encloses the spherical solid particle.The inner temperature of the insulation is the particle phase temperature T2, and the outer surface is exposed to a convection environment (boundary of the gas region) at .

It can be proven that the "critical radius of insulation" for the spherical particle is found as follows:

If the outer radius is less than the value given by this equation, then the heat transfer will be increased by adding more insulation. For outer radii greater than the critical value an increase in insulation thickness will cause a decrease in heat transfer.

The results shown in Fig .4 are related to the previous work in which the effect of radiation term was ignored (Fedorov and Shul'gin, 2006) and the results in Fig. 3 demonstrate the impact of this term in ignition of magnesium particle temperature.As seen, the highest temperature is achieved when the radiation term is considered in energy conservation equation and the lowest temperature is achieved without radiation term. Differences between values of gas temperature in Fig. 4 and Fig. 3 are also considerable.

In the previous model, heat transfer mechanism was just by conduction, however in this model, the conduction heat transfer and the heat transfers from surrounding gas radiation to particle is investigated.

Since the heat radiation is transferred from the surrounding gas to the gas and particle fuel, the gas and particle temperature increase due to the energy absorption. So most amount of energy transfers to the particle surface, and it is obvious that the temperature of the particle surface increases, due to gaining more heat, and physically this phenomenon is acceptable.

The investigation of variation of particle temperature versus the solid particle radius (), with various surrounding film thicknesses in the planar case shows that the influence of the radius of the metal partical on the variation of particle temperature appears to be weak for different values of .On the other hand we can assume that the particle temperature is approximately constant.

The non-dimensional gas phase temperature as a function of non-dimensional radius for values of non-dimensional gas layer thickness in planar case is shown in Fig.5; it can be noted that increasing the value of the radius results in a decrease of the gas phase temperature.

In Fig. 6, the value of the particle temperature in branch I (regular heating regime), is plotted as a function of for given values of in case of spherical coordinate.

Figure 3. Distribution of gas temperature for branch I of the ignition manifold for different values of in spherical case.

Figure 4. Distribution of gas temperature for branch I for different values of in spherical case without radiation in [6]

Figure 5. Distribution of gas temperature for branch I for in planar case

Figure 6. Distribution of particle temperature for branch I for different values of in spherical case.

Figure 7. Calculated values of the gas temperature as a function of particle radius for various values of in spherical case

Figure 8. Calculated values of the gas temperature as a function of particle radius for various values of in spherical case

In Fig. 7 the calculated gas temperature for Aluminium particle combustion is plotted as a function of the particle radius, ; it shows that for given value of the gas phase temperature gradient decreases with decreasing values of .

The gas phase temperature profile of heating of iron particle is plotted against the sample radius in Fig.8. In accordance with the condition in case of aluminium and magnesium particles, the gas temperature gradient decreases with the decrease of surrounding film thickness.


In this paper, an analytical model was presented to investigate the effect of radiation in the combustion of fine metal particles, in both planar and spherical coordinates. Also the mathematical model with allowance for the heterogeneous chemical reaction and the region of the thermal influence of particle on gas is developed to describe the thermal behaviour of a magnesium particle in a gas medium.

Considering the complications in the modeling of radiation term in the combustion of fine particles, some simplifications were made in order to find the influence of radiation on temperature profile of gas and particle phase.

By solving the energy equation in each phase, and equating the heat fluxes and temperatures at one side of the boundaries with their equivalents at the other side, we can obtain an algebraic relation for temperature and draw temperature profiles of gas and particle phases.

Due to the addition of the term of radiative heat transfer to the particle and gas phase energy conservation equations, the influence of radiation on temperature profiles was found to be an increase in temperature of the gas and particle phase as strong and weak respectively.

1. Bellman, R., Perturbation Techniques in Mathematics, Physics and Engineering, Holt, Rinehart & Winston, New York (1964).         [ Links ]
2. Cole, J.D., Perturbation Methods in Applied Mathematics, Blaisedell, Waltham, MA (1968).         [ Links ]
3. Fedorov, A.V., Physicomathematical modeling of ignition of fine magnesium particles, Preprint No. 12-94, Inst. Theor. Appl. Mech., Sib. Div., Acad. of Sci. of the USSR, Novosibirsk (1994).         [ Links ]
4. Fedorov, A.V., "Numerical and analytical study of magnesium particle ignition," Combust., Expl., Shock Waves, 32, 64-72 (1996).         [ Links ]
5. Fedorov, A.V., V.M. Fomin and Y.A. Gosteev, Theoretical Description of Ignition of Reactive Gas Suspensions, Izd. Nov. Gos. Arkh. Str. Univ., Novosibirsk (2003).         [ Links ]
6. Fedorov, A.V. and A.V.Shul'gin, "Conjugate Mathematical Model of Ignition of Magnesium Samples," Combustion, Explosion, and shock Waves, 42, 295-301 (2006).         [ Links ]
7. Naives, A.H., Introduction to Perturbation Techniques. Wiley, New York (1979).         [ Links ]
8. Nayfeh, A.H., Perturbation Methods. Wiley, New York (1973).         [ Links ]
9. O'Malley Jr., R.E., Introduction to Singular Perturbation, Academic Press, New York (1974).         [ Links ]
10. Yetter, R.A., and F.L. Dryer,Metal Particle Combustion and Classification, Micro-gravity Combustion: Fire in Free Fall, Academic Press (2001).         [ Links ]
11. Huang, Y., G.A. Risha, V. Yang and R.A. Yetter, Flame Propagation in Bimodal Nano/Micro-sized Aluminum Particles/Air Mixtures, 44th AIAA Aerospace Sciences Meeting and Exhibit (2006).         [ Links ]

Received: June 9, 2009.
Accepted: November 28, 2009.
Recommended by Subject Editor: Walter Ambrosini.

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License