Ahmed, Y., Gbeminiyi, S. (2019). On The Transient Combustion Analysis of Iron Micro-particles in a Gaseous Oxidizing Medium Using Variation of Parameter Method. Reliability Engineering and Resilience, 1(1), 1-11. doi: 10.22115/rer.2019.85314

Yinusa Amoo Ahmed; Sobamowo M Gbeminiyi. "On The Transient Combustion Analysis of Iron Micro-particles in a Gaseous Oxidizing Medium Using Variation of Parameter Method". Reliability Engineering and Resilience, 1, 1, 2019, 1-11. doi: 10.22115/rer.2019.85314

Ahmed, Y., Gbeminiyi, S. (2019). 'On The Transient Combustion Analysis of Iron Micro-particles in a Gaseous Oxidizing Medium Using Variation of Parameter Method', Reliability Engineering and Resilience, 1(1), pp. 1-11. doi: 10.22115/rer.2019.85314

Ahmed, Y., Gbeminiyi, S. On The Transient Combustion Analysis of Iron Micro-particles in a Gaseous Oxidizing Medium Using Variation of Parameter Method. Reliability Engineering and Resilience, 2019; 1(1): 1-11. doi: 10.22115/rer.2019.85314

On The Transient Combustion Analysis of Iron Micro-particles in a Gaseous Oxidizing Medium Using Variation of Parameter Method

^{}Department of Mechanical Engineering, University of Lagos, Akoka, Lagos, Nigeria

Receive Date: 08 March 2019,
Revise Date: 19 March 2019,
Accept Date: 19 March 2019

Abstract

Iron micro-particles when eject to the atmosphere from industries, factory and production enterprises pollute the atmosphere and generally affect the respiratory system. This paper presents a solution to such transient problem analytically by employing Variation of Parameter Method (VPM). For proper understanding of the problem chemistry and associated time of burning, the parameters involved are meticulously studied. The radiation property of the iron particle at high temperature as well as the impact of this temperature level on density are considered by incorporating radiation and linear micro particles density variation in the governing model. The obtained exact solution using VPM is verified with Runge-Kutta and also compared with the results of other works. Furthermore, error analysis is performed and discussed. The obtained result shows how the surrounding temperature and the heat realization term continue to influence combustible temperature history until an asymptotic behaviour is attained. It is envisaged that the present study will gain application in trying to annul some of the challenges the industries and firms have to overcome on combustion of solid combustibles like iron particles and also in the optimization of different particles burning time.

Factories and industries have been faced with the demerits associated with metallic combustibles like iron dust and also the challenges involve in predicting the burning time of the combustibles in question. The explosion hazards associated with micro and nano combustibles need serious attention so as to make different industrial operations applaudable. Consequently, researchers have attempted to better understand and model the particle combustion phenomenon. One of the attempt include the recent study of Sun et al. [1]. In their work, they carried out an investigation on the possible output when iron micro-particles are heated and suspended in air simultaneously. They made a conclusion on the possible range of combustion of the combustibles without phase flame by employing photomicrographs. The information gathered from the photomicrographs were then used to generate a relationship between the combustion period of iron micro-particles and diameter. Haghiri and Bidbadi [2] studied the dynamics of two phase mixture using the law of propagation of flame and accounted for thermal radiation. The mixture used includes iron micro dust and air. The result of their work shows that once the thermal radiation is augmented, particle burning speed and vaporization are exponentially enhanced. The work further shows the consideration of thermal radiation parameters when adequately considered compared with the case where this effect is correspondingly neglected. An examination was carried out by Hertzberg et al. [3] on the explosion limits of combustible micro-metallic when temperature and pressure are critically varied from the ambient condition. Argument were also presented on how the research can benefit industries that are directly affected by the research scope. Meanwhile, after the combustible dusts have been properly modelled, there is always need to solve the obtained nonlinear transient model in order to visualize and predict the thermal behaviour of the combustibles in question. As a result of the aforementioned, different analytical schemes have been adopted to obtain exact solution to the nonlinear model. An interesting attempt to solve the nonlinear model was presented by Hatami et al. [4]. They applied three different method of weighted residual and verified the results numerically. It was concluded that least square method displays better predictions as it agrees perfectly with numerical. However, in their work, no close form solution was obtained with all the parameters preserved. He [5–9] used perturbation method (HPM) to analyze and solve different nonlinear equations as applied to mathematics and engineering. He also shows ways of augmenting the scheme to handle strongly nonlinear engineering and mathematical problems. Other interesting approximate analytical methods such as HPM [10], DTM [11–13] have also been employed to investigate heat transfer in porous longitudinal fin as well as in convective heat transfer analysis considering internal heat generation. Sobamowo et al. [14] did a further study by considering fin with similar characteristics but with internal heat generation that is temperature-dependent and with the inclusion of a magnetic field. They presented interesting results and an elaborate parametric studies. The efficiency of the fin analyzed was also depicted and discussed. The obtained results using HPM were finally verified and a good agreement was established. In an attempt to obtain exact solution to partial differential equations, Yaseen et al. [15] solved the Laplace equation analytically by using DJ method. They considered different boundary conditions and also compared the obtained results with those of previous studies. Sobamowo and Yinusa [16,17] adopted the Adomian decomposition method and DJ method to study the thermal and environmental behaviour associated with combustible metallic micro-particles in an oxidizing medium. Furthermore, Yinusa and Sobamowo [18] did a recent study on the analysis of dynamic behaviour of a tensioned Carbon nanotube in thermal and pressurized environments. In their work, integral transform was employed to obtain the dynamic behaviour and stability criteria of the carbon nanotube. Nur et al. carried out an elementary investigation on the influence of tea dust geometry and sizeon its explosion by considering some food particles [19].

Motivated by the previous works, this work introduces an approximate analytical scheme for generating and monitoring combustibles temperature history during and after burning, so VPM is applied. The method agrees excellently with numerical Forth order Runge–Kutta and other methods with minimal error as compared to those of previous works; hence the verification of the schemes.

2. Description of problem and governing the equation

Consider a spherical shaped iron micro-particle as presented in Fig. 1. The iron dust has the capability to combust when in high reaction with an oxidizer.

By incorporating the necessary assumptions [16,17] and performing energy balance:

(1)

Fig. 1. Schematic of combusted iron particle in gaseous media [4].

Representing the energy term [4], we have as shown below:

(2)

(3)

(4)

(5)

put Eqs. (2), (3), (4) and (5) into Eq. (1), we obtain

(6)

In order to improve Eq. (6), the remaining assumptions are used:

Applying the above assumptions (j - l), Eq. (6) becomes,

(7)

Introducing the following dimensionless parameters,

(8)

Using Eq. (8) in Eq. (7), the dimensionless form of Eq. (7) becomes as presented by [4,16,17],

(9)

with initial condition,

(10)

3. Methods of solution

Due to the nonlinear terms in Eq. (9), it is very difficult to develop a closed form or an exact analytical solution to the nonlinear equation. Therefore, the common practice is to make recourse to numerical method. However, in recent time, several semi- or approximate analytical methods have been developed to solve nonlinear equations. In this present study, the nonlinear equation in Eq. (9) is be solved analytically using Variation of Parameter Method (VPM)

3.1. Basic principle: Variation of Parameter Method (VPM)

The basic concept and the procedure of variation of parameter method for solving differential equations are outlined as follows:

Generally, nonlinear equation may be operated to possess the form

(11)

The linear terms are decomposed into L + R, where

L is the highest order derivative which is easily invertible

R is the remainder of the linear operator of order less than L.

g is the system input or the source term

f is the system output,

Nf represents the nonlinear terms.

The variation of parameter method provides the general iterative scheme for Eq. (17) as:

(12)

where the initial approximation is given by

(13)

m is the order of the given differential equation,

k_{i }s are the unknown constants that can be determined by initial/boundary conditions and

is the multiplier that reduces the order of the integration and can be determined using Wronskian technique.

(14)

3.2. Method of solution: Variation of Parameter Method (VPM)

Recall that the nonlinear governing equation as shown in Eq. (9) may be expressed as,

(15)

(16)

Grouping the coefficients, the above equation becomes:

(17)

re-representing the coefficients, we have

(18)

where

(19)

Since the combustion model is:

Applying the VPM principle described above, we have for the governing equation,

(20)

Using the initial condition and performing the VPM iteration, the term by term solution becomes,

(21)

(22)

(23)

4. Results and Discussion

The verification of the analytical method used (VPM) with a numerical Runge-Kutta is depicted in Fig. 2. The scheme, VPM ascertain a very close agreement with numerical method. In order to visualize the accuracy of the scheme, a super-imposed plot which shows the temperature profile of a 20μm combusting iron particle is inspected as displayed in Fig. 2 and also Table 1. From the figure, it is evident that VPM agrees excellently with the Numerical scheme and as such is efficient for the problem in concern.

Fig. 2. Verification of VPM with Numerical.

4.1. Effect of particle diameter on the temperature history

Fig. 3 and Fig. 4 depict the influence of combusting particle diameter on temperature profile and burning rate using VPM. From the graphs, it is shown that particle diameters have evident influence on the temperature profile. A particle with 60μm diameter was observed to possess a higher temperature profile which means that increasing the combusting particle diameter causes a corresponding increase in the temperature profile and the burning time. As a result of this evident impact, the particle diameter may be used as a controlling agent in reducing the hazardous effects that normally propagate from iron particle combustion.

Fig. 3. Effect of particle diameter on the temperature Profile with DJM.

4.2. Effect of and on the temperature history

Fig. 5 and Fig. 6 depict the influence of and on the temperature profile. From the figures, it can be seen that increasing and decreases the combustion temperature with this effect more pronounced with . The decrease in combustion temperature with a corresponding increase in and is as a result of an increase in the radiation heat transfer term in the combustion particle.

Fig. 5. Effect of on the temperature profile.

Fig. 6. Effect of on the temperature profile.

4.3. Effect of heat realization and surrounding temperature on the temperature history

Fig. 7 and Fig. 8 display the impact of the heat realization term and the temperature of the surrounding on the combustion temperature. From the plots, it can be concluded that that increasing the heat realization term and the temperature of the surrounding increases the combustion temperature. This increase is significant for the heat realization term variation than that of the surrounding temperatures except for high values of surrounding temperature. The physical implication of the combined increase in both components result in micro dusts of very high temperature and latent energy which when inhale may be hazardous. As a result, it is very important to release such combustibles only into properly controlled or monitored domain.

Fig. 7. Effect of heat realized term on the temperature profile.

Fig. 8. Effect of surrounding temperature on the temperature profile.

Table. 1. Comparism of the analytical scheme with a numerical method for a 20μm iron particle.

τ

θ (τ) for a Particle diameter of (μm)

Numerical

ADM [15]

DJM [16]

VPM (Present study)

0.0

1.000000000000000

1.000000000000000

1.000000000000000

1.000000000000000

0.1

1.113333969181095

1.113333981177264

1.113333969181095

1.113333969181095

0.2

1.215117348008129

1.215117711150980

1.215117348008133

1.215117348008129

0.3

1.306590963764452

1.306593345739770

1.306590963764455

1.306590963764452

0.4

1.388844932815013

1.388852820269113

1.388844932815017

1.388844932815014

0.5

1.462843666959360

1.462858827578177

1.462843666959367

1.462843666959360

0.6

1.529449065177065

1.529458070391791

1.529449065177069

1.529449065177064

0.7

1.589444597054137

1.589387100510296

1.589444597054139

1.589444597054137

0.8

1.643562984179010

1.643271057270020

1.643562984179017

1.643562984179010

0.9

1.692520185796611

1.691615617727115

1.692520185796617

1.692520185796612

1.0

1.737058395009000

1.734792471017498

1.737058395009000

1.737058395009000

In addition,

The residual error associated with the method of solution for the 20μm case is plotted as shown below:

Fig. 8. Residual Error.

5. Conclusion

In this work, a study has been carried out for the determination of the temperature history of iron particle during combustion process. The results of the VPM solution were verified numerically. It was established that VPM gives a good result and is efficient for the problem investigated. Also, parametric studies were performed to fully understand how the combusting particle diameter, density, radiative term, heat realized term and other parameters affect the burning time as well as the combustion temperature. The results revealed that by increasing the heat realized parameter, combustion temperature increased until a steady state was reached. It was also achieved that particle diameter have evident influence on the temperature profile. A particle with 60μm diameter was observed to possess a higher temperature profile which means that an increase in the combusting particle diameter causes a corresponding increase in the temperature profile as well as the burning time. As a result of this evident impact, the particle diameter may be used as a controlling agent in reducing the hazardous effects that normally propagate from iron particle combustion. It is hoped that the present study will enhance the understanding of the combustion of the particle and also obviate the challenges facing industries on combustion of metallic particles such as iron particles as well as in the determination of different particles burning time.

References

[1] SUN J-H, DOBASHI R, HIRANO T. Combustion Behavior of Iron Particles Suspended in Air. Combust Sci Technol 2000;150:99–114. doi:10.1080/00102200008952119.

[2] Haghiri A, Bidabadi M. Dynamic behavior of particles across flame propagation through micro-iron dust cloud with thermal radiation effect. Fuel 2011;90:2413–21. doi:10.1016/j.fuel.2011.03.032.

[3] Hertzberg M, Zlochower IA, Cashdollar KL. Metal dust combustion: Explosion limits, pressures, and temperatures. Symp Combust 1992;24:1827–35. doi:10.1016/S0082-0784(06)80214-6.

[4] Hatami M, Ganji DD, Jafaryar M, Farkhadnia F. Transient combustion analysis for iron micro-particles in a gaseous media by weighted residual methods (WRMs). Case Stud Therm Eng 2014;4:24–31. doi:10.1016/j.csite.2014.06.003.

[5] He J-H. Homotopy perturbation technique. Comput Methods Appl Mech Eng 1999;178:257–62. doi:10.1016/S0045-7825(99)00018-3.

[6] HE J-H. ADDENDUM: NEW INTERPRETATION OF HOMOTOPY PERTURBATION METHOD. Int J Mod Phys B 2006;20:2561–8. doi:10.1142/S0217979206034819.

[7] He J-H. A coupling method of a homotopy technique and a perturbation technique for non-linear problems. Int J Non Linear Mech 2000;35:37–43. doi:10.1016/S0020-7462(98)00085-7.

[8] HE J-H. SOME ASYMPTOTIC METHODS FOR STRONGLY NONLINEAR EQUATIONS. Int J Mod Phys B 2006;20:1141–99. doi:10.1142/S0217979206033796.

[9] HE J-H. A NEW PERTURBATION TECHNIQUE WHICH IS ALSO VALID FOR LARGE PARAMETERS. J Sound Vib 2000;229:1257–63. doi:10.1006/jsvi.1999.2509.

[10] Saedodin S, Shahbabaei M. Thermal Analysis of Natural Convection in Porous Fins with Homotopy Perturbation Method (HPM). Arab J Sci Eng 2013;38:2227–31. doi:10.1007/s13369-013-0581-6.

[11] Darvishi MT, Gorla R, Khani F, Aziz A. Thermal performance of a porus radial fin with natural convection and radiative heat losses. Therm Sci 2015;19:669–78. doi:10.2298/TSCI120619149D.

[12] Moradi A, Hayat T, Alsaedi A. Convection-radiation thermal analysis of triangular porous fins with temperature-dependent thermal conductivity by DTM. Energy Convers Manag 2014;77:70–7. doi:10.1016/j.enconman.2013.09.016.

[13] Hoshyar H, Ganji DD, Abbasi M. Determination of Temperature Distribution for Porous Fin with Temperature-Dependent Heat Generation by Homotopy Analysis Method. J Appl Mech Eng 2015;04. doi:10.4172/2168-9873.1000153.

[14] Sobamowo G, Adeleye O, Yinusa A. Analysis of convective-radiative porous fin With temperature-dependent internal heat Generation and magnetic field using Homotopy Perturbation method. J Comput Appl Mech 2017;12:127–45.

[15] Yaseen M, Samraiz M, Naheed S. Exact solutions of Laplace equation by DJ method. Results Phys 2013;3:38–40. doi:10.1016/j.rinp.2013.01.001.

[16] Sobamowo GM, Yinusa A. Transient Combustion Analysis for Iron Micro-Particles in a Gaseous Oxidizing Medium Using Adomian Decomposition Method. J Comput Eng Phys Model 2018;1:1–15. doi:10.22115/cepm.2018.122052.1012.

[17] Sobamowo GM, Yinusa AA. Transient Combustion Analysis of Iron Micro-Particles in a Gaseous Oxidizing Medium Using a New Iterative Method. J Comput Eng Math 2018;5:3–16. doi:10.14529/jcem180301.

[18] Yinusa A, Sobamowo G. Analysis of Dynamic Behaviour of a Tensioned Carbon Nanotube in Thermal and Pressurized Environments. Karbala Int J Mod Sci 2019;5. doi:10.33640/2405-609X.1015.

[19] Nur Hikmah S, Siti Zubaidah S, Noor Sabrina AM, Siti Kholijah AM, Rohaida CM, Shalyda MS, et al. Preliminary study on the tea dust explosion: the effect of tea dust size. MATEC Web Conf 2019;255:02014. doi:10.1051/matecconf/201925502014.