
Research Article


Perturbation Method as a Powerful Tool to Solve Highly Nonlinear Problems: The Case of Gelfand’s Equation 

U. FilobelloNino,
H. VazquezLeal,
K. Boubaker,
Y. Khan,
A. PerezSesma,
A. SarmientoReyes,
V.M. JimenezFernandez,
A. DiazSanchez,
A. HerreraMay,
J. SanchezOrea
and
K. PereyraCastro



ABSTRACT

Solving nonlinear ordinary differential equations is relevant because phenomena on the frontiers of modern sciences are often nonlinear in nature; therefore this article proposes Perturbation Method (PM) to solve nonlinear problems. As case study PM is employed to obtain a handy approximate solution for Gelfand’s differential equation which governing combustible gas dynamics. Comparing figures between approximate and exact solutions, it is shown that PM method result extremely efficient.




How
to cite this article:
U. FilobelloNino, H. VazquezLeal, K. Boubaker, Y. Khan, A. PerezSesma, A. SarmientoReyes, V.M. JimenezFernandez, A. DiazSanchez, A. HerreraMay, J. SanchezOrea and K. PereyraCastro, 2013. Perturbation Method as a Powerful Tool to Solve Highly Nonlinear Problems: The Case of Gelfand’s Equation. Asian Journal of Mathematics & Statistics, 6: 7682. DOI: 10.3923/ajms.2013.76.82 URL: https://scialert.net/abstract/?doi=ajms.2013.76.82



Received: February 28, 2013;
Accepted: April 27, 2013;
Published: July 18, 2013


INTRODUCTION As it is known, Gelfand’s equation (also known as Bratu’s problem in 1D) models the chaotic dynamics in combustible gas thermal ignition. Therefore, it is important to search for accurate solutions for this equation. Unfortunately, it is difficult to solve, like many others nonlinear differential equations that appear in engineering.
The Perturbation Method (PM) is a well established method; it is among the
pioneer techniques to approach various kinds of nonlinear problems. This procedure
was originated by S.D. Poisson and extended by J.H. Poincare. Although the method
appeared in the early 19th century, the application of a perturbation procedure
to solve nonlinear differential equations was performed later on that century.
The most significant efforts were focused on celestial mechanics, fluid mechanics
and aerodynamics (Chow, 1995; Holmes,
1995).
In a broad sense, it is possible to express a nonlinear differential equation
in terms of one linear part and other nonlinear. The nonlinear part is considered
as a small perturbation through a small parameter (the perturbation parameter).
The assumption that the nonlinear part is small compared to the linear is considered
as a disadvantage of the method. There are other modern alternatives to find
approximate solutions to the differential equations that describe some nonlinear
problems such as those based on: variational approaches (Kazemnia
et al., 2008; Noorzad et al., 2008),
tanh method (Evans and Raslan, 2005), expfunction (Mahmoudi
et al., 2008), Adomian’s decomposition method (Kooch
and Abadyan, 2011, 2012; Vanani
et al., 2011; Chowdhury, 2011), homotopy
perturbation method (Ganji et al., 2008, 2009;
Sharma and Methi, 2011; VazquezLeal
et al., 2012a, b; FilobelloNino
et al., 2012a, c; Khan
and Wu, 2011; Mirgolbabaei and Ganji, 2009; Tolou
et al., 2008), homotopy analysis method (Patel
et al., 2012) and Boubaker polynomials expansion scheme (Agida
and Kumar, 2010; Ghanouchi et al., 2008),
variational iteration method (Saravi et al., 2013),
among many others.
Although the PM method provides in general, better results for small perturbation
parameters ε<<1 , it will be seen that the approximation obtained,
besides to be handy has a good accuracy, even for relatively large values of
the perturbation parameter (FilobelloNino et al.,
2012b).
BASIC IDEA OF PERTURBATION METHOD Let the differential equation of one dimensional nonlinear system be in the form: where, it is assumed that x is a function of one variable x = x(t), L(x) is a linear operator which, in general, contains derivatives in terms of t, N(x) is a nonlinear operator and ε is a small parameter. Considering the nonlinear term in (1) to be a small perturbation and assuming that the solution for (1) can be written as a power series in the small parameter ε: After substituting (2) into (1) and equating terms having identical powers of ε, it is possible to obtain a number of differential equations that can be integrated, recursively, to find the values for the functions:
x_{0}(t), x_{1}(t), x_{2}(t)...
APPROXIMATE SOLUTION OF GELFAND’S EQUATION The equation to solve is: where, ε is a positive parameter. It is possible to find a handy solution for (3) by applying the PM method. Identifying terms: where, prime denotes differentiation respect to x. To solve (3), first it is required expanding the exponential term of Gelfand^{,}s problem, resulting: identifying ε with the PM parameter, then it is assumed a solution for (6) in the form: Equating the terms with identical powers of ε it can be solved for y_{0}(x), y_{1}(x), y_{2}(x).. and so on. Later it will be seen that, a very good handy result is obtained, by keeping up to third order approximation:
Thus, the solution for the lowest order approximation is:
On the other hand the solution of (9) is given by: The solving process for (10), leads to result: In the same way, from (11): and so on. _{}By substituting (1215) into (7) it is obtained a handy third order approximation for the solution of (3), as it is shown: Considering as a case study, the value of the Gelfand’s parameter ε = 1, (16) adopts the form: DISCUSSION The fact that the PM depends on a parameter which is assumed small, suggests that the method is limited. In this work, the PM method has been applied to the problem of finding an approximate solution for the highly nonlinear Gelfand’s differential equation. This equation is relevant because it describes the dynamics in combustible gas thermal ignition. Figure 1 shows the comparison between approximation (17) with the four order Runge Kutta (RK4) numerical solution. It can be noticed that (17) sketches successfully the qualitative paraboliclike behavior of Bratu’s problem. This characteristics can be deduced from (3), because it implies that y”(x)<0 for 0≤x≤1 and hence the solution of Gelfand’s equation has to be concave downward in the same interval. As a matter of fact Fig. 2, shows that the maximum absolute error is about 5E4 at x = 0.5, which proves the efficiency of PM method, especially because it was only considered the thirdorder approximation; by using four terms in series expansion of nonlinear part εe^{y(x)}.
The PM method provides in general, better results for small perturbation parameters
ε<<1, (Eq. 1) and when are included the most number
of terms from (Eq. 2).

Fig. 1: 
Approximation of Gelfand’s equation and RK4 comparison
for e = 1 

Fig. 2: 
Absolute error for approximation (17) (e = 1) 

Fig. 3: 
Approximation of Gelfand’s equation and RK4 comparison
for e = 1.5 
To be precise, ε is a parameter of smallness, that measure how greater
is the contribution of linear term L(x) than the one of N(x) in (1). Figure
1 and 2 show a noticeable fact, that (Eq.
17) provides a good approximation as a solution of (3), despite of the fact
that perturbation parameter ε = 1 cannot be considered small. Moreover,
the Fig. 3 shows that (Eq. 16) is in good
agreement with RK4, for ε = 1.5.
Saravi et al. (2013) employed VIM method to
provide an approximate solution of (Eq. 3). The maximum absolute
error reported for the VIM solution is 1.59355E3 at x = 0.5, which shows the
high accuracy of PM in comparison with VIM. In addition, VIM solution generates
a more complicated solution for similar number of iterations.
Finally, our approximate solution (Eq. 16) does not depend
of any adjustment parameter, for which, it is in principle, a general expression
for Gelfand^{,}s problem.
CONCLUSION In this study, PM was presented to construct analytical approximate solutions of Gelfand’s equation in the form of rapidly convergent series. The success of the method for this case it has to be considered as a possibility to apply it in other non linear problems, instead of using other sophisticated and difficult methods. From Fig. 13, it is deduced that the proposed solutions are highly accurate. ACKNOWLEDGMENTS The authors gratefully acknowledge the financial support provided by the National Council for science and Technology of Mexico (CONACyT) through Grant CB201001 no. 157024. The authors would like to express their gratitude to RogelioAlejandro CallejasMolina and Roberto RuizGomez for their contribution to this study.

REFERENCES 
1: Agida, M. and A.S. Kumar, 2010. A boubaker polynomials expansion scheme solution to random Love equation in the case of a rational kernel. Electronic J. Theoretical Phys., 7: 319326. Direct Link 
2: Chow, T.L., 1995. Classical Mechanics. John Wiley and Sons Inc., USA
3: Chowdhury, S.H., 2011. A comparison between the modified homotopy perturbation method and adomian decomposition method for solving nonlinear heat transfer equations. J. Applied Sci., 11: 14161420. CrossRef 
4: Evans, D.J. and K.R. Raslan, 2005. The Tanh function method for solving some important nonlinear partial dfferential equation. Int. J. Computat. Math., 82: 897905. CrossRef  Direct Link 
5: FilobelloNino, U., H. VazquezLeal, R. CastanedaSheissa, A. Yildirim and L. HernandezMartinez et al., 2012. An approximate solution of blasius equation by using HPM method. Asian J. Math. Stat., 5: 5059. CrossRef 
6: FilobelloNino, U., H. VazquezLeal, Y. Khan, A. Yyldyrym and D. PereyraDiaz et al., 2012. HPM applied to solve nonlinear circuits: A study case. Applied Math. Sci., 6: 43314344. Direct Link 
7: FilobelloNino, U., H. VazquezLeal, Y. Khan, A. Yildirim and V.M. JimenezFernandez et al., 2012. Perturbation method and LaplacePade approximation to solve nonlinear problems. Miskolc Mathematical Notes. (In Press). Direct Link 
8: Ganji, D.D., H. Babazadeh, F. Noori, M.M. Pirouz and M. Janipour, 2009. An application of homotopy perturbation method for non linear blasius equation to boundary layer flow over a flat plate. Int. J. Nonlinear Sci., 7: 309404. Direct Link 
9: Ganji, D.D., H. Mirgolbabaei, M. Miansari and M. Miansari, 2008. Application of homotopy perturbation method to solve linear and nonlinear systems of ordinary differential equations and differential equation of order three. J. Applied Sci., 8: 12561261. CrossRef  Direct Link 
10: Ghanouchi, J., H. Labiadh and K. Boubaker, 2008. An Attempt to solve the heat transfer equation in a model of pyrolysis spray using 4qOrder Boubaker polynomials. Int. J. Heat Technol., 26: 4953.
11: Holmes, M.H., 1995. Introduction to Perturbation Methods. SpringerVerlag, New York
12: Kazemnia, M., S.A. Zahedi, M. Vaezi and N. Tolou, 2008. Assessment of modified variational iteration method in BVPs highorder differential equations. J. Applied Sci., 8: 41924197. CrossRef  Direct Link 
13: Khan, Y. and Q. Wu, 2011. Homotopy perturbation transform method for nonlinear equations using He's polynomials. Comput. Math. App., 61: 19631967. CrossRef 
14: Kooch, A. and M. Abadyan, 2011. Evaluating the ability of modified Adomian decomposition method to simulate the instability of freestanding carbon nanotube: comparison with conventional decomposition method. J. Applied Sci., 11: 34213428. CrossRef 
15: Kooch, A. and M. Abadyan, 2012. Efficiency of modified adomian decomposition for simulating the instability of nanoelectromechanical switches: Comparison with the conventional decomposition method. Trends Applied Sci. Res., 7: 5767. CrossRef  Direct Link 
16: Mahmoudi, J., N. Tolou, I. Khatami, A. Barari and D.D. Ganji, 2008. Explicit solution of nonlinear ZKBBM wave equation using expfunction method. J. Applied Sci., 8: 258363. CrossRef  Direct Link 
17: Mirgolbabaei, H. and D.D. Ganji, 2009. Application of homotopy perturbation method to solve combined korteweg de vriesmodified korteweg de vries equation. J. Applied Sci., 9: 35873592. CrossRef  Direct Link 
18: Noorzad, R., A.T. Poor and M. Omidvar, 2008. Variational iteration method and homotopyperturbation method for solving burgers equation in fluid dynamics. J. Applied Sci., 8: 369373. CrossRef  Direct Link 
19: Patel, T., M.N. Mehta and V.H. Pradhan, 2012. The numerical solution of burger's equation arising into the irradiation of tumour tissue in biological diffusing system by homotopy analysis method. Asian J. Appl. Sci., 5 : 6066. CrossRef 
20: Saravi, M., M. Hermann and D. Kaiser, 2013. Solution of Bratu's equation by He's variational iteration method. Am. J. Comput. Applied Mathe., 3: 4648.
21: Sharma, P.R. and G. Methi, 2011. Applications of homotopy perturbation method to partial differential equations. Asian J. Math. Stat., 4: 140150. CrossRef 
22: Tolou, N., J. Mahmoudi, M. Ghasemi, I. Khatami, A. Barari and D.D. Ganji, 2008. On the nonlinear deformation of elastic beams in an analytic solution. Asian J. Scient. Res., 1: 437443. CrossRef  Direct Link 
23: Vanani, S.K., S. Heidari and M. Avaji, 2011. A lowcost numerical algorithm for the solution of nonlinear delay boundary integral equations. J. Applied Sci., 11: 35043509. CrossRef  Direct Link 
24: VazquezLeal, H., U. FilobelloNino, R. CastanedaSheissa, L. Hernandez Martinez and A. SarmientoReyes, 2012. Modified HPMs inspired by homotopy continuation methods. Math. Problems Eng.,, Vol. 2012. CrossRef  Direct Link 
25: VazquezLeal H., R. CastanedaSheissa, U. FilobelloNino, A. SarmientoReyes and J. SanchezOrea, 2012. High accurate simple approximation of normal distribution related integrals. Math. Problems Eng., Vol. 2012. CrossRef  Direct Link 



