Reach Us
+44-1474-556909

**Ophir Nave ^{1*}, Yaron Lehavi^{2}, Vladimir Gol’dshtein^{1} and Suraj Ajadi^{3}**

^{1}Department of Mathematics, Ben-Gurion University of the Negev, PO Box 653, Beer-Sheva, 84105, Israel

^{2}Jerusalem College of Technology (JCT), Department of Physics David Yellin, Jerusalem, Israel

^{3}Faculty of Science, Department of Mathematics, O.A.U, Ile-Ife, Nigeria

- *Corresponding Author:
- Ophir Nave

Department of Mathematics

Ben-Gurion University of the Negev

PO Box 653, Beer-Sheva, 84105, Israel

**Fax:**972-8-647-7648

**E-mail:**[email protected]

**Received Date:** April 20, 2012; **Accepted Date:** June 20, 2012; **Published Date:** June 23, 2012

**Citation:** Nave O, Lehavi Y, Gol’dshtein V, Ajadi S (2012) Application of the Homotopy Perturbation Method (HPM) and the Homotopy Analysis Method (HAM) to the Problem of the Thermal Explosion in a Radiation Gas with Polydisperse Fuel Spray. J Appl Computat Math 1:115. doi: 10.4172/2168-9679.1000115

**Copyright:** © 2012 Nave O, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Applied & Computational Mathematics

The aim of this work is to apply the homotopy perturbation method and homotopy analysis method to the problem of thermal explosion in a flammable gas mixture with the addition of volatile fuel droplets. The system of equations that describes the effects of heating, evaporation, and combustion of fuel in a polydisperse spray is simplified. Both convective and radiative heating of droplets is taken into account in the model. The model for the radiative heating of droplets takes into account the semitransparency of the droplets. The results of the analysis have been applied to the modeling of the thermal explosion in diesel engines. We applied the Homotopy Perturbation Method and the Homotopy Analysis Method to the new model and we found the region of the convergence of the considered solutions of the relevant physical parameters. The results demonstrate that these methods are very effective for solving nonlinear problems in science and engineering.

Homotopy perturbation method (HPM); Homotopy analysis method (HAM); Nonlinear differential equations; Thermal explosion; Polydisperse fuel spray

• *A* pre-exponential rate factor ( *s*^{−1} )

• *B *universal gas constant ( *Jkmol*^{−1}K^{−1} )

•* C* molar concentration ( *kmolm*^{−3} )

• *c* specific heat capacity ( *Jkg*^{−1}K^{−1})

• *E* activation energy ( *Jkmol*^{−1} )

• *F* −*M* refer to the model (21)-(24) with b = 0

• *k* number of droplets

•* L* liquid evaporation energy (i.e., latent heat of evaporation, Enthalpy of evaporation) ( *Jkg*^{−1} )

• *m* droplet mass

• *n* number of droplets per unit volume (*m*^{−3} )

• *Q* combustion energy ( *Jkg*^{−1} )

• *q* heat flux (*Wm*^{−2} )

• *R* radius of droplet (*m* )

• *r *dimensionless radius

• *T* temperature ( *K* )

• *t* time ( s )

• t_{react} characteristic reaction time ( *s* ) defined in Equation (25)

**Greek symbols**

• α dimensionless volumetric phase content

• β dimensionless reduced initial temperature (with respect to the so-called activation temperature E / B )

• *γ* dimensionless parameter that represents the reciprocal of the final dimensionless adiabatic temperature of the thermally insulated system after the explosion has been completed

• ε_{1i} i =1,..., k dimensionless parameters defined in Equation (25) and describes the competition between the combustion and the evaporation processes

• ε_{2i} i =1,..., k dimensionless parameters defined in Equation (25). This parameter relates the heat released during combustion and energy that is needed to evaporate all the fuel droplets

• ε_{3i} i =1,..., k dimensionless parameters defined in Equation (25). This parameter describes the ratio of treact and the characteristic droplet heating time

• ε_{4i} i =1,..., k dimensionless parameters defined in Equation (25). This parameter is proportional to the ratio of radiative and convective fluxes

• *η *dimensionless fuel concentration

• *θ* dimensionless temperature

• *λ* thermal conductivity ( *Wm*^{−1}*K*^{−1} )

• *μ* molar mass ( *kgkmol*^{−1} )

• *ν* dimensionless parameter defined in Equation (25)

• *ρ* density ( *kgm*^{−3} )

• *σ* Stefan-Boltzmann constant (*W* *m*^{−2}*K*^{−4} )

• *τ* dimensionless time

• *ψ* represents the internal characteristics of the fuel (the ratio of the specific combustion energy and the latent heat of evaporation) defined in Equation (25) and for diesel fuel *ψ* >>1

**Subscripts**

• *con* convection

• *d* liquid fuel droplets

•* f* combustible gas component of the mixture

• *g* gas mixture

• *l* liquid phase

• *p* under constant pressure

•* rad *adiation

• 0 initial state

In this paper we investigated the problem of thermal explosion in a fuel mixture and gas. This problem, in most cases, has been studied based on the application of computational fluid dynamics (CFD) packages [1]. This method could take into account the complicated geometry of the enclosure and the chemistry of the processes. Hence, this makes it particularly attractive for engineering applications including the modeling of combustion processes in diesel engines. Other approaches to this problem are based on a asymptotic analysis of equations describing the limiting cases of the processes. One of these approaches is based on the application of the zero-order approximation of the geometric version of the asymptotic method of integral manifolds (MIM), developed for combustion applications in [2,3]. For example, in [4] the method of integral manifold was applied to a specific problem of modeling of ignition process in a diesel engines by using the P-1 model. The chemical term was presented in the Arrhenius form with the pre-exponential factor calculated from the enthalpy equation, using the well known Shell autoignition model. The results predicted by the analytical solution were compared with those that predicted by the computational fluid dynamics package VECTIS. The effects of the thermal radiation were shown to be very significant especially at high temperatures.

Other asymptotic methods were proposed by [5] in 1992, known as the homotopy analysis method (HAM) and by [6] in 1998, known as the homotopy perturbation method (HPM), which is a special case of HAM. The HPM and the HAM methods are mathematical tools that are based on homotopy, a fundamental concept in topology and differential geometry. They are analytical approaches to formulate the series solution of linear and nonlinear partial differential equations. We refer the reader to [7] for an enlightening comparison between HAM and HPM.

HPM couples the homotopy technology and perturbation method including the modified Lindstedt-Poincare method [8]. The authors of [9] modified the Multiple Scales method by incorporating the time transformation of Lindstedt Poincare method. In [10] the authors contrasted two different approaches of Lindstedt-Poincare methods using the duffing equation. The main deficiencies in applying perturbation methods is that a small parameter is needed in the equations.

The HPM was further developed and improved and applied to nonlinear oscillators with discontinuities [11], nonlinear wave equations [12], boundary value problem [13], limit cycle and bifurcation of nonlinear problems [14] and many other subjects. In recent years, the application of the homotopy perturbation method (HPM) in nonlinear problems has been developed by scientists and engineers, because this method deforms the difficult problem under study into a simple problem which is easy to solve [15-17]. Most perturbation methods assume a small parameter exists, but most nonlinear problems have no small parameter at all. Unlike analytical perturbation methods, the HPM and HAM do not depend on a small parameter which is difficult to find.

These two methods also provide a simple way to ensure the convergence of the series solution. Moreover, these methods provide a large degree of freedom to choose an appropriate base functions to approximate the linear and nonlinear problems [18]. Another important advantage of this method is that one can construct a continuous mapping of an initial guess approximation to the exact solution of the given problem through an auxiliary linear operator. To ensure the convergence of the series solution an auxiliary parameter is used. In [19] Liao has substantiated that the HAM differs from the other analytical methods in that it ensures the convergence of the series solution by choosing a proper value for the convergence-control parameter.

In this paper we have rewritten the model that was proposed by [20] for polydisperse fuel spray and applied the HPM and HAM to the problem of thermal explosion in a fuel mixture and gas. Based on these two methods, we present an analytical solutions for various values of the relevant physical parameter and we discuss the convergence of these solutions. We also compare our results to numerical solutions.

To explain this method, let us consider the following equation:

(1)

with the boundary conditions of:

(2)

where *A* , *B* , *f (r)* are a general differential operator, a boundary operator, a known analytical function respectively. Ω is the domain. Generally, the operator *A* can be decomposed into a linear part *L* and a nonlinear part *N(u)* . Hence, Equation (1) can be written as:

(3)

By the homotopy technique, we construct a homotopy which satisfies:

(4)

where is an embedding parameter and *α *(*r, p*) is a function of r and p , and u_{0}(r) denote the initial approximation of α (r) .

When *p* = 0 we have

(5)

and when *p* =1 we have

(6)

As we mention before, L denote an auxiliary linear operator. In addition L have the property:

(7)

Using (7), it is clear that for *p* = 0

(8)

is the solution of the equation:

(9)

And for p =1

(10)

is the solution of the equation:

(11)

When the embedding parameter* p* increase from 0 to1, the solution *α* (*r*, *p*) of the equation:

(12)

depends upon the embedding parameter p and the varies from the initial approximation u_{0} (r) to the solution u(r) of equation (1). In topology, such a kind of continuous variation is called deformation.

According to the HPM, we can first use the embedding parameter *p* as a “small parameter”, and assume that the solutions of Equation (4) can be written as a power series in *p* :

(13)

Setting *p* =1 yields in the approximate solution of (1) to

(14)

The combination of the perturbation method and the homotopy method is called the HPM, which eliminates the drawbacks of the traditional perturbation methods while keeping all their advantages. The rate of convergent of series (13) depends on the nonlinear operator *A*(*u*) .

The physical assumptions are as follows: The combustible gas mixture contains evaporating ideal spherical droplets of fuel. The liquid droplets form a polydisperse spray. The medium is assumed to be spatially homogeneous. The variations in pressure in the enclosure, and their influence on the combustion processes are ignored. The heat flux from the burning gas to the droplets is assumed to consist of two components: convection and radiation, and the form of these two components is as follows:

(15)

(16)

where for *μm* units the value of and are: and . The energy that is needed for heating fuel vapor from the droplet temperature to gas temperature is ignored. The thermal conductivity of the liquid phase is much greater than that of the gas phase. The volume fraction of the liquid phase is much less than that of the gas phase. The heat transfer coefficient in the liquidgas mixture is assumed to be controlled by the thermal properties of the gas phase. External heat losses are ignored. Fuel droplets are semitransparent. Combustion takes place in the gas phase only. Combustion is modeled as a one-step first-order exothermic reaction with gaseous fuel as a deficient reactant. The droplets are assumed to be stationary. Under these assumptions, we rewrite the model as in [21], which is in the form of monodisperse fuel spray, to a polydisperse fuel spray as follows:

(17)

(18)

(19)

(20)

In non-dimensional parameters and by applying the Frank- Kamenetskii approximation [22], the model has the form of:

(21)

(22)

(23)

(24)

where the following dimensionless parameters have been introduced:

(25)

The non-dimensional initial conditions are:

(26)

For simplicity, in our numerical simulations and when applying the HPM and the HAM we assume *b* = 0 .

**Application of HPM to the Problem of the Thermal Explosion in Polydisperse Fuel Spray**

By applying the HPM method to the system of Equations (21)-(24) we obtain the following HPM-system:

(27)

(28)

(29)

(30)

with the initial conditions:

(31)

According to HPM method the terms α_{1} , α_{2}, α_{3}^{i} and α_{4}^{i} for i =1,..., k has the form of:

(32)

Suppose that the solution of (27)-(30) takes the form of

(33)

(34)

(35)

(36)

Substituting Equations (32) with the initial conditions (31) into Equations (27)-(30) for three different size of droplets i.e., *k* = 3 , using the Taylor expansion for the exponent [23,24] and finally collecting the terms in power of *p* up to order 3 we obtain:

equations for gas temperature θ_{g}

(37)

(38)

(39)

equations for concentration *η*

(40)

(41)

(42)

(43)

(44)

(45)

(46)

(47)

(48)

We derived a system of (24) ordinary differential equation with 24 with unknown functions: α_{1,1} , α_{1,2} , α_{1,3} , α_{2,1} , α_{2,2} , α_{2,3} , α^{1}_{3,1,} α^{1}_{3,2, }α^{1}_{3,3, }α^{2}_{3,1, }α^{2}_{3,2, }α^{2}_{3,3. }α^{3}_{3,1, }α^{3}_{3,2, }α^{3}_{3,3, }α^{1}_{4,1}, α^{1}_{4,2, }α^{1}_{4,3, }α^{2}_{4,1,} α^{2}_{4,2,} α^{2}_{4,3,} α^{3}_{4,1,} α^{3}_{4,2,} α^{3}_{4,3 . } The initial conditions are:

(49)

In this section, we discuss the convergence of the HAM solutions. The convergence depends on the so-called convergence - control - parameter , and so, we plot the -curve for θ_{g}(0) , η (0) , θ_{d}(0), and r_{d} (0). The interval of convergence is determined by the flat portion of the -curve. In order to plot the -curve we applied the HAM as given in [24] to our new model (21)-(24).

**An introduction to Homotopy Analysis Method (HAM)**

Consider the following differential equation:

(50)

where *N* is a nonlinear operator, is a vector of spatial variables, *t* denotes time and *u* is an unknown function.

**Zero order deformation of HAM:** By means of generalizing the traditional concept of homotopy, Liao [24] constructs the so-called zero-order deformation equation:

(51)

where is a non-zero auxiliary parameter, *H* is an auxiliary function, is an auxiliary linear operator, is an initial guess of ; is a unknown function. The degree of freedom is to choose the initial guess, the auxiliary linear operator, the auxiliary parameter, and the auxiliary function *H* . Expanding in Taylor series with respect to the embedding parameter p , one has

(52)

where

(53)

If the auxiliary linear operator, the initial guess, the auxiliary parameter, and the auxiliary function are so properly chosen that the above series converges at *p* =1, one has

(54)

which must be one of the solutions of the original nonlinear equation, as proved in [24]. If the same initial guess and the same auxiliary linear operator are chosen, the approximations given by the homotopy perturbation method are exactly a special case of those given by the homotopy analysis method when = −1 and *H* =1 . The series (54) itself is in principle a kind of Taylor series (at p = 1). Hence, mathematically, homotopy perturbation method itself is also a kind of generalized Taylor technique.

**mth-order deformation:** Define the vector:

(55)

Differentiating Equation (51) *m* -times with respect to the embedding parameter *p* and then setting *p* = 0 and finally dividing the terms by *m!* , we obtain the *m* th-order deformation equation in the form of:

(56)

where,

(57)

and m χ_{m} is the unit step function. Applying the inverse operator on both side of Equation (56), we get

(58)

In this way, it is easy to obtain um for *m* ≥1 , at *m* th-order and finally get the solution as:

(59)

In our model we choose the initial guess to be , , , and which satisfied the initial conditions. The linear operator will be:

(60)

with the property , where c_{1} and c_{2} are constants of integration. According to the system (21)-(24) and the terms as in Equation (32) the nonlinear operators will be defined as follows:

(61)

(62)

(63)

(64)

By substituting the series (32) into Equations (61)-(65) correspondingly we get the terms for R_{m} accordingly to Equation (57) as follows:

(65)

(66)

(67)

(68)

Now, the solution of the *m* th-order deformation can be expressed according to Equation (58) which can be solved by symbolic software such as Mathematica 8.0, Maple, Matlab and so on. We obtain a family of solutions that depends on the auxiliary parameter . So, regarding as independent variable, it is easy to plot the curves. For example, we can plot the curves:

(69)

The curves (i =1, 2,3) are a function of and thus can be plotted by a curve . According to [20] there exists a horizontal line segment (flat portion of the -curve) in the figure of and called *the valid region* of which corresponds to a region of convergent of the solutions. Thus, if we choose any value in the valid region of we are sure that the corresponding solutions series are convergent. For given initial approximation , the auxiliary linear operator , and the auxiliary function , the valid region of for different special quantities are often nearly the same for a given problem. Hence, the so-called -curve provides us a convenient way to show the influence of on the convergence region of the solutions series.

We compared the system dynamics of the models (21)-(24) and (37)-(48) with and without the impact of the thermal radiation. The results are based on the following diesel engines parameter values:

*Diesel − engines*

(70)

We studied the problem of the the effect of fuel spray polydispersity on the ignition process in a fuel cloud by applying numerical simulation, the homotopy perturbation method and the homotopy analysis method for 30th order deformation. We compared between the homotopy perturbation method and by solving the full system of the model i.e., the system of Equations: (21)-(24) for the solution profiles of the of the gas temperature, droplet temperature, radius and concentration numerically. Although we take into account only three different size of droplets, our results show that the homotopy perturbation method provides an excellent approximation of the solutions of the system with high accuracy (**Figures** **1**-**4**).

**Figure 1:** Solution profiles of the gas temperature θ_{g} −τ ; 1: Full model solved
numerically with the impact of the thermal radiation, 2: HPM model with the impact
of the thermal radiation for the = −1, (2): HAM model with the impact of
the thermal radiation for the = 0.04 , 3: Full model solved numerically without
the impact of the thermal radiation, 4: HPM model without the impact of the
thermal radiation, (4): HAM model without the impact of the thermal radiation
for = 0.04 .

**Figure 2:** Solution profiles of the droplet temperature
θ_{d} −τ ; 1: Full model
solved numerically with the impact of the thermal radiation, 2: HPM model
with the impact of the thermal radiation for the = −1, (2): HAM model with
the impact of the thermal radiation for the = 0.04 , 3: Full model solved
numerically without the impact of the thermal radiation, 4: HPM model without
the impact of the thermal radiation, (4): HAM model without the impact of the
thermal radiation for = 0.04 .

**Figure 3:** Solution profiles of the radius r −τ ; 1: Full model solved numerically
with the impact of the thermal radiation, 2: HPM model with the impact
of the thermal radiation for the = −1 , (2): HAM model with the impact
of the thermal radiation for the = 0.04 , 3: Full model solved numerically
without the impact of the thermal radiation, 4: HPM model without the impact
of the thermal radiation, (4): HAM model without the impact of the thermal
radiation for = 0.04 .

**Figure 4:** Solution profiles of the concentration η −τ ; 1: Full model solved
numerically with the impact of the thermal radiation, 2: HPM model with the
impact of the thermal radiation for the = −1, (2): HAM model with the impact
of the thermal radiation for the = 0.04 , 3: Full model solved numerically
without the impact of the thermal radiation, 4: HPM model without the impact
of the thermal radiation, (4): HAM model without the impact of the thermal
radiation for = 0.04 .

The gas temperature trajectory, **Figure 1**, for all models with and without the impact of the thermal radiation, starts with a fast increase in the temperature until then the temperature decreases from , which means cooling before ignition, until . This continuous process of cooling before ignition is summarized in **Table 1** and corresponding in **figures** **1**-**4**. This dimensionless time, *τ* , refers to the ignition time and it is compatible for all the solution profiles for the gas and droplets temperature, radius and concentration. According to these results, the F −M has the smallest ignition time with and without the impact of the thermal radiation when comparing to HPM and HAM. The HPM results are closer to the F −M than the HAM results with and without the impact of the thermal radiation.

Model | radiation | no-radiation |

1: Full-model | =0.446306 | =0.562209 |

2: HPM | =0.482483 | =0.585683 |

3: HAM | =0.625222 | =0.726507 |

**Table 1:** The process of cooling before the ignition time for the different models. *F* −*M* refer to the full-model, HPM with = −1 and HAM with = 0.04 .

The presence of the small parameter *γ* in the gas temperature equations, such that Equations (37)-(39) form a singularly perturbed system, enables one to exploit the geometrical version of the method of the integral manifold and hence to separate the model into fast and slow subsystems. According to this method, the above results are also sustained with MIM [25].

In order to clarify the influence of the thermal radiation on the dimensionless time before the final explosion of the system we introduce the term *impact of the thermal radiation* as given in [26] which is measured in percent and defined as:

(71)

(72)

(73)

here the subscript *F* −*M* refers to the full model solved numerically, *τ* is the ignition time, and the superscript *rad* and *no − rad* refer to the model with and without the impact of the thermal radiation respectively. We also defined the parameters , , and which show the difference between the different models in percent. The results are as follows: , , and , , and . As we can see from these results, the HPM model is closer to the model solved with numerical simulations than the HAM model (the comparison between the *F* −*M* and the HPM model based on the value of as equal to −1, and the comparison between the *F* −*M* and the HAM model based on the value of as equal to 0.04 ).

As we mentioned in the previous section, the convergence depends on the convergence-control-parameter , and so we plot the -curve for and as shown in **Figure 5** for m = 30 in Equation (59) i.e. 30 th order approximation. According to **Figure 5** , the interval of convergence that agrees for all of the corresponding solutions is . In order to emphasize the impact of the convergence-control parameter on the solutions profiles we defined the terms:

(74)

(75)

(76)

(77)

which point out the difference (in dimensionless-time) between the profiles solutions of the HPM and HAM methods for different with the impact of the thermal radiation and without the impact of the thermal radiation respectively from the numerical results, and τ refer to the ignition time. The results are summarized in **Table 2**. According to these results, the HPM and the HAM methods are closed to the numerical results for both with and without the impact of the thermal radiation.

Parameter | |

1: _rad^=0.04=|_HAM^=0.04-_F-M| | 0.17891 |

2: _rad^=-1=|_HPM^=-1-_F-M| | 0.03617 |

3: _no-rad^=0.04=|_HAM^=0.04-_F-M| | 0.16429 |

4: _no-rad^=-1=|_HPM^=-1-_F-M| | 0.02347 |

**Table 2:** The impact of the convergence-control parameter on the solutions
profiles with and without the impact of the thermal radiation comparing to the numerical
simulations.

We have shown the the solutions obtained by HPM and HAM are convergent and that they extremely well with numerical simulations. It has also been shown that the homotopy perturbation method, which is a special case of the homotopy analysis method when = −1 , yields convergent solutions for all of the cases considered. These results demonstrate that HPM and HAM are very effective analytical methods for solving nonlinear problems in science and engineering.

Our next step in this direction is to apply the HPM to the continuous model as in our previous work [27,28].

- Aggrawal SK (1998) A review of spray ignition phenomena: present status and future research. Prog Energy Combust Sci 24: 565-600.
- Babushok VI, Gol'dshtein V (1988) Structure of the thermal explosion limit.Combust and Flame 72: 221-226.
- Gol'dshtein VM, Sobolev VA (1992) Singularity theory and some problems of functional analysis.American Math Soci 153 :73-92.
- Sazhin SS, Feng G, Heikal MR (1992) A model for fuel spray penetration. Fuel 80: 2171-2180.
- Liao SJ (1992) The proposed homotopy analysis technique for the solution of nonlinear problems. Ph.D Thesis Shanghai Jiao Tong University.
- He JH (1998) Approximate analytical solution for seepage flow with fractional derivatives in porous media. Comput Meth Appl Mech Eng 16: 57-68.
- Liao SJ (2005) Comparison between the homotopy analysis method and homotopy perturbation method. Appl Math Comp169: 1186-1194.
- Liao SJ(2007) A short communication on Dr. He's modified Lindstedt-Poincare method NONLINEAR DYNAMICS49: 317-318.
- Pakdemirli M, Karahan MMF, Boyac H (2009) A new perturbation algorithm with better convergence properities: multiple scales Lindstedt-Poincare method Math Comp App 14: 31-44.
- Hu H, Xiong ZG (2004) Comparison of two Lindstedt-Poincare type perturbation methods. J Sound and Vibration 278: 437-444.
- He JH (2004) The homotopy perturbation method for non-linear oscillators with discontinuities. Appl Math Comput 151: 287-292.
- He JH (2005) Application of homotopy perturbation method to nonlinear wave equations. Chaos Solitons Fractals 26: 695-700.
- He JH (2006) Homotopy perturbation method for solving boundary problems. Phys Lett A 350: 87-88.
- He JH (2005) Limit cycle and bifurcation of nonlinear problems. Chaos Solitons Fractals 26: 827-833.
- He JH (2004) Asymptotology by homotopy perturbation method. Appl Math Comput 156: 591-596.
- He JH (2005) Limit cycle and bifurcation of nonlinear problems. Chaos Solitons Fractals 26: 827-833.
- He JH (2005) Homotopy-perturbation method for bifurcation of nonlinear problems. Internat J Nonlinear Sci Numer Simul 6: 207-208.
- Liao SJ (2004) On the homotopy analysis method for nonlinear problems. App Math Comp 147: 499-513.
- Liao SJ (2009) Notes on the homotopy analysis method: Some definition and theorems. Comm in Nonlinear Sci Numerical Simu 14: 983-997.
- Goldfarb I, Gol'dshtein V, Zinoviev A (2002) Delayed thermal explosion in porous media: Method of invariant manifolds. IMA J App Math 67: 263-280.
- Frank DA-Kamenetskii (1969) Diffusion and heat exchange in chemical kinetics. Plenum Press New York
- Ajadi SO, Zuilino M (2011) Approximate analytical solutions of reaction–diffusion equations with exponential source term: Homotopy perturbation method (HPM). App Math Let 24: 1634-1639
- Turgut O, Ahmet Y (2008) Comparison between Adomian's method and He's homotopy perturbation method. Comp Math App 56: 1216-1224.
- Liao SJ (2003) Beyond Perturbation: Introduction to the Homotopy Analysis Method. Chapman and Hall/CRC Press. Boca Raton.
- Goldfarb I, Gol'dshtein V, Greenberg BJ, Kuzmenko G (2002) Analysis of Radiative Heat-Loss Effects in Thermal Explosion of a Gas using the Integral Manifold Methods. J Engng Math 44: 229-243.
- Goldfarb ,Gol'dshtein VM, Katz D, Sazhin S (2007) Radiation effect on thermal explosion in a gas containing evaporation fuel droplets. Int Journal of Thermal Sci 46: 358-370.
- Nave O, Gol'dshtein VM , Vitcheslav B (2010) A probabilistic model of thermal explosion in polydisperse fuel spray. Applied Math Comp 217: 2698-2709.
- Nave O,Gol'dshtein VM (2011) The flammable spray effect on thermal explosion of combustible gas-fuel mixture. Combus Sci Tech 183: 519-539.

Select your language of interest to view the total content in your interested language

- Adomian Decomposition Method
- Algebraic Geometry
- Analytical Geometry
- Applied Mathematics
- Axioms
- Balance Law
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Complex Analysis
- Computational Model
- Convection Diffusion Equations
- Cross-Covariance and Cross-Correlation
- Differential Equations
- Differential Transform Method
- Fourier Analysis
- Fuzzy Boundary Value
- Fuzzy Environments
- Fuzzy Quasi-Metric Space
- Genetic Linkage
- Hamilton Mechanics
- Hypothesis Testing
- Integrated Analysis
- Integration
- Large-scale Survey Data
- Matrix
- Microarray Studies
- Mixed Initial-boundary Value
- Molecular Modelling
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Nonlinear Differential Equations
- Number Theory
- Numerical Solutions
- Physical Mathematics
- Quantum Mechanics
- Quantum electrodynamics
- Quasilinear Hyperbolic Systems
- Regressions
- Relativity
- Riemannian Geometry
- Robust Method
- Semi Analytical-Solution
- Sensitivity Analysis
- Smooth Complexities
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Theoretical Physics
- Theory of Mathematical Modeling
- Three Dimensional Steady State
- Topology
- mirror symmetry
- vector bundle

- Total views:
**13174** - [From(publication date):

September-2012 - Jul 21, 2019] - Breakdown by view type
- HTML page views :
**9257** - PDF downloads :
**3917**

**Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals**

International Conferences 2019-20