ISRN Condensed Matter Physics

Volume 2013 (2013), Article ID 129571, 8 pages

http://dx.doi.org/10.1155/2013/129571

Research Article

Temperature Variations Analysis for Condensed Matter Micro- and Nanoparticles Combustion Burning in Gaseous Oxidizing Media by DTM and BPES

M. Hatami,¹ D. D. Ganji,¹ and K. Boubaker²

¹Department of Mechanical Engineering, Babol University of Technology, P.O. Box 484, 47148-71167 Babol, Iran

²Equipe de Physique des Dispositifs à Semiconducteurs, Faculté des Sciences de Tunis, Campus Universitaire, 2092 Tunis, Tunisia

Received 29 May 2013; Accepted 23 July 2013

Academic Editors: C. Andreani and J. Tempere

Copyright © 2013 M. Hatami et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Combustion process for iron particles burning in the gaseous oxidizing medium is investigated using the Boubaker polynomial expansion scheme (BPES) and the differential transformation method (DTM). Effects of thermal radiation from the external surface of burning particle and alterations of density of iron particle with temperature are considered. The solutions obtained using BPES technique and DTM are compared with those of the fourth-order Runge-Kutta numerical method. Results reveal that BPES is more accurate and reliable method than DTM. Also the effects of some physical parameters that appeared in mathematical section on temperature variations of particles as a function of time are studied.

1. Introduction

Combustion of metallic particles is one of the most challenging issues in industries that manufacture, process, generate, or use combustible dusts, and an accurate knowledge of their explosion hazards is essential. Many studies have been done for estimating and modeling the particle and dust combustion. Haghiri and Bidbadi [1] investigated the dynamic behavior of particles across flame propagation through a two-phase mixture consisting of micro-iron particles and air. They assumed three zones for flame structure: preheat, reaction, and postflame (burned). Liu et al. [2] analyzed the flame propagation through hybrid mixture of coal dust and methane in a combustion chamber. A one-dimensional, steady-state theoretical analysis of flame propagation mechanism through microiron dust particles based on dust particles’ behavior with special remark on the thermophoretic force in small Knudsen numbers is presented by Bidabadi et al. [3]. Haghiri and Bidabadi [4] performed a mathematical model to analyze the structure of flame propagating through a two-phase mixture consisting of organic fuel particles and air. In contrast to previous analytical studies, they take thermal radiation effect into consideration, which has not been attempted before. Recently, Bidbadi and Mafi [5] solved the nonlinear energy equation that resulted from particle combustion modeling by using homotopy perturbation method (HPM), and they presented equations for calculating the convective heat transfer coefficient and burning time for iron particles. Because HPM needs perturbation and a small parameter, it can be solved by other high-accuracy analytical methods where in the present study two of them are presented and compared with fourth-order Runge-Kutta numerical method.

Polynomial expansion methods are extensively used in many mathematical and engineering fields to yield meaningful results for both numerical and analytical analyses. Among the most frequently used polynomials, the Boubaker polynomial expansion scheme which is firstly introduced by Boubacker [6] is one of the interesting tools which were associated with several applied physical and mathematical problems [7–26].

The differential transformation method is an alternative procedure for obtaining an analytic Taylor series solution of differential equations. The main advantage of this method is that it can be applied directly to nonlinear differential equations without requiring linearization, discretization, and therefore, it is not affected by errors associated with discretization. The concept of DTM was first introduced by Zhou [27], who solved linear and nonlinear problems in electrical circuits. Chen and Ho [28] developed this method for partial differential equations, and Ayaz [29] applied it to the system of differential equations; this method is very powerful [30]. Jang et al. [31] applied the two-dimensional DTM to partial differential equations. This method was successfully applied to various application problems [32, 33].

Motivated by previously mentioned works, this paper aims to introduce two analytical methods for obtaining the temperature of iron particle during combustion. So BPES and DTM are presented. These methods have an excellent agreement with numerical Runge-Kutta method; they also have very low errors without any need for perturbation or discretization compared to previous analytical methods in the literature.

2. Problem Description

Consider a spherical particle which due to high reaction with oxygen will be combusted. Since the thermal diffusivity of substance is large and Biot number is small (), it is assumed that the particle is isothermal. In this state, a lumped system analysis is applicable. When this criterion is satisfied, the variation of temperature with location within the particle will be slight and can be approximated as being uniform, so particle has a spatially uniform temperature, and therefore, the temperature of particle is a function of time only, , and is not a function of radial coordinate, . The assumptions used in this modeling are [5] as follows.

(1)The spherical particle burns in a quiescent, infinite ambient medium, and there are no interactions with other particles, and also the effects of forced convection are ignored.(2)Thermophysical properties for the particle and ambient gaseous oxidizer are assumed to be constant.(3)The particle radiates as a gray body to the surroundings without the contribution of the intervening medium.

By these assumptions and considering the particle as a thermodynamic system and by the use of principle of conservation of energy (first law of thermodynamics), the energy balance equation for this particle can be written as

where is the rate of energy entering the system which is due to absorption of total radiation incident on the particle surface from the surrounding, is the rate of energy leaving the system by mechanisms of convection on the particle surface and thermal radiation that emits from the outer surface of particle, is the rate of generation of energy inside the particle due to the combustion process and is equal to the heat released from the chemical reaction, and is the rate of change in total energy of particle. These energy terms can be calculated by [5]

By substituting (2)–(5) in (1),

Three reasonable assumptions are used for improving (6).

(I)Both absorptivity and emissivity of the surface depend on the temperature and the wavelength of radiation. Kirchhoff’s law of radiation states that the absorptivity and the emissivity of a surface at a given temperature and wavelength are equal ().(II)The initial temperature of the particle at the beginning of combustion can be regarded as the initial condition. This temperature is known as ignition temperature .(III)The density of particle is a function of particle temperature, so it can be considered as a linear function ().

By applying these assumptions, (6) will be converted to the following:

For solving this nonlinear differential equation, It is more suitable that all the terms be converted to the dimensionless form. The following set of dimensionless variables are defined:

Consequently, the nonlinear differential equation and its initial condition can be expressed in the dimensionless form

3. Applied Analytical Methods

In this section two analytical methods called differential transformation method (DTM) and Boubaker polynomial expansion scheme (BPES) with their application in the problem are presented.

3.1. Differential Transformation Method (DTM)

For understanding the method’s concept, suppose that is an analytic function in domain and represents any point in the domain. The function is then represented by one power series whose center is located at . The Taylor series expansion function of is in form of

The Maclaurin series of can be obtained by taking in (11) expressed as

As explained in [27], the differential transformation of the function is defined as follows:

where represents the transformed function and is the original function. The differential spectrum of is confined within the interval, where is a constant value. The differential inverse transformation of is defined as follows:

It is clear that the concept of differential transformation is based the Taylor series expansion. The values of function at values of argument are referred to as discrete; that is, is known as the zero discrete, as the first discrete, and so forth. The more discreteness available, the more precise it is possible to restore the unknown function. The function consists of the -function , and its value is given by the sum of the -functions with as its coefficient. In real applications, at the right choice of constant , the larger values of argument the discrete of spectrum reduce rapidly. Some important mathematical operations performed by differential transformation method are listed in Table 1. By applying DTM from Table 1, transformed form of (9) will be

where is transformed form of and

Table 1: Some fundamental operations of the differential transformation method.

Origin function
Transformed function

Transformed form of initial condition (see (10)) will be

For example, for an iron particle with 20 μm diameter (see Table 2) solving (15) makes

Table 2: Properties and conditions for combustion of iron particles.

Particle diameter (μm)

20
0.051595
0.002630
0.98579
1.17647
0.35294

60
0.051595
0.007567
2.83636
1.17647
0.35294

100
0.051595
0.011493
4.30774
1.17647
0.35294

By substituting DTM transformed terms of (see (18)) into (13), can be determined as

3.2. Boubaker Polynomials Expansion Scheme (BPES)

The resolution of system (9) and (10) along with boundary conditions has been achieved using the Boubaker polynomials expansion scheme BPES [7–26]. This scheme is a resolution protocol which has been successfully applied to several applied physics and mathematics problems. The BPES protocol ensures the validity of the related boundary conditions regardless of main equation features. The Boubaker polynomials expansion scheme BPES is based on the Boubaker polynomials first derivatives properties:

Several solutions have been proposed through the BPES in many fields such as numerical analysis [7–10], theoretical physics [10–13], mathematical algorithms [14], heat transfer [15], homodynamic [16, 17], material characterization [18], fuzzy systems modeling [19–24], and biology [25, 26].

The resolution protocol is based on settingas an estimator to the-dependent variable:

where are the -order Boubaker polynomials [7, 8], are minimal positive roots, is a prefixed integer, andare unknown pondering real coefficients.

The main advantage of this formulation is the verification of boundary conditions, expressed in (9), in advance to resolution process. In fact, thanks to the properties expressed in (21)-(22), these conditions are reduced to the inherently verified linear equations:

The BPES solution for (9) is obtained, according to the principles of the BPES, by determining the nonnull set of coefficients that minimize the absolute difference between left and right sides of the following equations, which follow a majoring the quadrature term within the sum:

The final solution is obtained by substituting the obtained values of the coefficients in (22).

4. Results and Discussion

As described in problem description section, the equation of combustion for a particle was introduced as (9) with boundary condition equation (10). Some constant parameters for this equation are introduced in Table 2 for iron particle. For these three particle diameters, (9) is solved by DTM and BPES whose results are presented through Figures 1(a) and 1(b), respectively. As seen in these figures, both of the two analytical methods have very excellent agreement with numerical fourth-order Runge-Kutta method, and also these figures reveal that when particle diameter increases, an increase in combustion temperature due to higher energy realized in combustion is observed. The figures indicate that the temperature of iron particle (for all sizes) reaches the maximum value at the end of the burning time, and the chemical reaction, is finished, and larger particles have higher maximum temperature than that of smaller particles. This result is similar to the experimental result presented by Tang et al. [34] and analytical method by Bidabadi and Mafi [5]. Obtained data from applied methods for a special case ( μm) are compared in Table 3. This table confirms that BPES has lower errors compared to DTM, but both methods are convenient and accurate. Effects of and on nondimensional temperature profile for nanoparticles are shown in Figures 2(a) and 2(b), respectively. As seen, the increase in makes an increase in temperature profile, but makes a decrease due to the increase in radiation heat transfer term in the particle (see (8)).

Table 3: Comparison of DTM and BPES results with fourth-order Runge-Kutta numerical method for diameter 20 μm.

Num
DTM
BPES
% Error DTM
% Error BPES

0.0
1.00000
1.00000
1.00000000000
—
—

0.05
1.05696555223261
1.056965552
1.05696555214

0.10
1.11094655557423
1.110946451
1.11094651374

0.15
1.16210940722201
1.162109321
1.16210937273

0.20
1.21060998272779
1.210609952
1.21060997044

0.25
1.25659429223923
1.256594206
1.25659425774

0.30
1.30019895535529
1.300198823
1.30019890241

0.35
1.34155222687316
1.341552153
1.34155219732

0.40
1.38077485817965
1.380774801
1.38077483531

0.45
1.41798035901816
1.417980211
1.41798029981

0.50
1.45327532764821
1.453275203
1.45327527779

(a)

(b)

Figure 1: (a) Comparison between DTM and numerical method in different particle diameters (see Table 2), (b) BPES results for particle temperature for  μm.

(a)

(b)

Figure 2: (a) Effect of on nondimensional temperature profile for micro- and nanoparticles, (b) effect of on nondimensional temperature profile for micro- and nanoparticles.

An important point for particle combustion is surrounding temperature in which particle is combusted. As (2) reveals, energy inputs from the surrounding occur through radiation heat transfer, and when surrounding has higher temperature, input energy increases and consequently temperature in combustion of particle will increase. This effect is depicted in Figure 3(a). Figure 3(b) demonstrates the effect of heat realized parameter (Ψ) on temperature profile versus time. It is completely evident that by increasing the generated heat in the combustion, temperature will increase significantly.

(a)

(b)

Figure 3: (a) Effect of surrounding temperature on nondimensional temperature profile, (b) effect of heat realized parameter (Ψ) on nondimensional temperature profile.

5. Conclusion

In this study, the equation of temperature variation in combustion process for particles is presented. The effects of thermal radiation from the external surface of burning particle and alterations of density of iron particle with temperature are considered. Due to nonlinearity of described equation, Boubaker polynomial expansion scheme (BPES) and differential transformation method (DTM) have been presented in order to obtain analytical solutions. Results show that both methods have good agreement with numerical method, but BPES had lower errors than DTM. Also the effects of surrounding temperature, realized heat parameter, particle diameter, and convection and radiation heat transfer parameter on combustion temperature have been investigated.

References

A. Haghiri and M. Bidabadi, “Dynamic behavior of particles across flame propagation through micro-iron dust cloud with thermal radiation effect,” Fuel, vol. 90, no. 7, pp. 2413–2421, 2011.

Y. Liu, J. Sun, and D. Chen, “Flame propagation in hybrid mixture of coal dust and methane,” Journal of Loss Prevention in the Process Industries, vol. 20, no. 4–6, pp. 691–697, 2007.

M. Bidabadi, A. Haghiri, and A. Rahbari, “Mathematical modeling of velocity and number density profiles of particles across the flame propagation through a micro-iron dust cloud,” Journal of Hazardous Materials, vol. 176, no. 1–3, pp. 146–153, 2010.

A. Haghiri and M. Bidabadi, “Modeling of laminar flame propagation through organic dust cloud with thermal radiation effect,” International Journal of Thermal Sciences, vol. 49, no. 8, pp. 1446–1456, 2010.

M. Bidabadi and M. Mafi, “Time variation of combustion temperature and burning time of a single iron particle,” International Journal of Thermal Sciences, vol. 65, pp. 136–147, 2013.

K. Boubaker, “The Boubaker polynomials, a new function class for solving bi-varied second order differential equations,” Far East Journal of Applied Mathematics, vol. 31, pp. 299–320, 2008.

M. Agida and A. S. Kumar, “A Boubaker Polynomials Expansion Scheme solution to random Love's equation in the case of a rational kernel,” Electronic Journal of Theoretical Physics, vol. 7, no. 24, pp. 319–326, 2010.

A. Yildirim, S. T. Mohyud-Din, and D. H. Zhang, “Analytical solutions to the pulsed Klein-Gordon equation using Modified Variational Iteration Method (MVIM) and Boubaker Polynomials Expansion Scheme (BPES),” Computers and Mathematics with Applications, vol. 59, no. 8, pp. 2473–2477, 2010.

J. Ghanouchi, H. Labiadh, and K. Boubaker, “An attempt to solve the heat transfer equation in a model of pyrolysis spray using 4q-order m-Boubaker polynomials,” International Journal of Heat and Technology, vol. 26, no. 1, pp. 49–53, 2008.

S. Slama, J. Bessrour, K. Boubaker, and M. Bouhafs, “A dynamical model for investigation of A₃ point maximal spatial evolution during resistance spot welding using Boubaker polynomials,” The European Physical Journal of Applied Physics, vol. 44, no. 3, pp. 317–322, 2008.

S. Slama, M. Bouhafs, and K. B. Ben Mahmoud, “A Boubaker polynomials solution to heat equation for monitoring A₃ point evolution during resistance spot welding,” International Journal of Heat and Technology, vol. 26, no. 2, pp. 141–145, 2008.

S. Lazzez, K. B. Ben Mahmoud, S. Abroug, F. Saadallah, and M. Amlouk, “A Boubaker Polynomials Expansion Scheme (BPES)-related protocol for measuring sprayed thin films thermal characteristics,” Current Applied Physics, vol. 9, no. 5, pp. 1129–1133, 2009.

T. Ghrib, K. Boubaker, and M. Bouhafs, “Investigation of thermal diffusivity-microhardness correlation extended to surface-nitrured steel using Boubaker polynomials expansion,” Modern Physics Letters B, vol. 22, no. 29, pp. 2893–2907, 2008.

S. Fridjine, K. B. Ben Mahmoud, M. Amlouk, and M. Bouhafs, “A study of sulfur/selenium substitution effects on physical and mechanical properties of vacuum-grown ZnS1-xSex compounds using Boubaker Polynomials Expansion Scheme (BPES),” Journal of Alloys and Compounds, vol. 479, no. 1-2, pp. 457–461, 2009.

C. Khélia, K. Boubaker, T. Ben Nasrallah, M. Amlouk, and S. Belgacem, “Morphological and thermal properties of β-SnS2 sprayed thin films using Boubaker polynomials expansion,” Journal of Alloys and Compounds, vol. 477, no. 1-2, pp. 461–467, 2009.

K. B. Ben Mahmoud and M. Amlouk, “The 3D Amlouk-Boubaker expansivity-energy gap-Vickers hardness abacus: a new tool for optimizing semiconductor thin film materials,” Materials Letters, vol. 63, no. 12, pp. 991–994, 2009.

M. Dada, O. B. Awojoyogbe, and K. Boubaker, “Heat transfer spray model: an improved theoretical thermal time-response to uniform layers deposit using Bessel and Boubaker polynomials,” Current Applied Physics, vol. 9, no. 3, pp. 622–624, 2009.

S. Tabatabaei, T. Zhao, O. Awojoyogbe, and F. Moses, “Cut-off cooling velocity profiling inside a keyhole model using the Boubaker Polynomials Expansion Scheme,” International Journal of Heat and Mass Transfer, vol. 45, no. 10, pp. 1247–1251, 2009.

A. Belhadj, J. Bessrour, M. Bouhafs, and L. Barrallier, “Experimental and theoretical cooling velocity profile inside laser welded metals using keyhole approximation and Boubaker polynomials expansion,” Journal of Thermal Analysis and Calorimetry, vol. 97, no. 3, pp. 911–915, 2009.

A. Belhadj, O. F. Onyango, and N. Rozibaeva, “Boubaker polynomials expansion scheme-related heat transfer investigation inside keyhole model,” Journal of Thermophysics and Heat Transfer, vol. 23, no. 3, pp. 639–640, 2009.

P. Barry and A. Hennessy, “Meixner-type results for Riordan arrays and associated integer sequences,” Journal of Integer Sequences, vol. 13, no. 9, pp. 1–34, 2010.

A. S. Kumar, “An analytical solution to applied mathematics-related Love's equation using the Boubaker Polynomials Expansion Scheme,” Journal of the Franklin Institute, vol. 347, no. 9, pp. 1755–1761, 2010.

S. Fridjine and M. Amlouk, “A new parameter: an abacus for optimizing PVT hybrid solar device functional materials using the Boubaker Polynomials Expansion Scheme,” Modern Physics Letters B, vol. 23, no. 17, pp. 2179–2191, 2009.

M. Benhaliliba, C. E. Benouis, K. Boubaker, M. Amlouk, and A. Amlouk, “A new guide to thermally optimized doped oxides monolayer spray-grown solar cells: the Amlouk-Boubaker optothermal expansivity ψab,” in Solar Cells—New Aspects and Solutions, L. A. Kosyachenko, Ed., pp. 27–41, InTech, 2011.

A. Milgram, “The stability of the Boubaker Polynomials Expansion Scheme (BPES)-based solution to Lotka-Volterra problem,” Journal of Theoretical Biology, vol. 271, no. 1, pp. 157–158, 2011.

H. Rahmanov, “A solution to the non linear Korteweg-De-Vries equation in the particular case dispersion-adsorption problem in porous media using the spectral Boubaker Polynomials Expansion Scheme (BPES),” Studies in Nonlinear Sciences, vol. 2, no. 1, pp. 46–49, 2011.

K. Zhou, Differential Transformation and Its Applications for Electrical CircuIts, Huazhong University Press, Wuhan, China, 1986.

C. K. Chen and S. H. Ho, “Solving partial differential equations by two-dimensional differential transform method,” Applied Mathematics and Computation, vol. 106, no. 2-3, pp. 171–179, 1999.

F. Ayaz, “Solutions of the system of differential equations by differential transform method,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 547–567, 2004.

I. H. A.-H. Hassan, “Comparison differential transformation technique with adomian decomposition method for linear and nonlinear initial value problems,” Chaos, Solitons and Fractals, vol. 36, no. 1, pp. 53–65, 2008.

M.-J. Jang, C.-L. Chen, and Y.-C. Liu, “Two-dimensional differential transform for partial differential equations,” Applied Mathematics and Computation, vol. 121, no. 2-3, pp. 261–270, 2001.

A. A. Joneidi, D. D. Ganji, and M. Babaelahi, “Differential transformation method to determine fin efficiency of convective straight fins with temperature dependent thermal conductivity,” International Communications in Heat and Mass Transfer, vol. 36, no. 7, pp. 757–762, 2009.

M. Momeni, N. Jamshidi, A. Barari, and G. Domairry, “Numerical analysis of flow and heat transfer of a viscoelastic fluid over a stretching sheet,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 21, no. 2, pp. 206–218, 2011.

F. D. Tang, S. Goroshin, and A. J. Higgins, “Modes of particle combustion in iron dust flames,” Proceedings of the Combustion Institute, vol. 33, no. 2, pp. 1975–1982, 2011.