An Analytical Approximate Solution for the Bratu Problem by using Nonlinearities Distribution Homotopy Perturbation Method

DOI : 10.17577/IJERTV10IS050432

Download Full-Text PDF Cite this Publication

  • Open Access
  • Authors : Mario A. Sandoval-Hernandez, U. Filobello-Nino, H. Vazquez-Leal , S. E. Torreblanca Bouchan, O. Alvarez-Gasca, A. D. Contreras-Hernandez , B. E. Palma-Grayeb, Elisa De Leo Baquero, Alexis C. Bielma Perez , Julio C. Vichi Mendoza, J. Sanchez-Orea, R. A. Callejas-Molina, C. Hoyos-Reyes
  • Paper ID : IJERTV10IS050432
  • Volume & Issue : Volume 10, Issue 05 (May 2021)
  • Published (First Online): 03-06-2021
  • ISSN (Online) : 2278-0181
  • Publisher Name : IJERT
  • License: Creative Commons License This work is licensed under a Creative Commons Attribution 4.0 International License

Text Only Version

An Analytical Approximate Solution for the Bratu Problem by using Nonlinearities Distribution Homotopy Perturbation Method

Mario A. Sandoval-Hernandez1, U. Filobello-Nino2, H. Vazquez-Leal2*, S. E. Torreblanca Bouchan1,

  1. Alvarez-Gasca2, A. D. Contreras-Hernandez2, B. E. Palma-Grayeb2, Elisa De Leo Baquero1, Alexis C. Bielma Pérez1, Julio C. Vichi Mendoza1, J. Sanchez-Orea2, R. A. Callejas-Molina3,

    C. Hoyos-Reyes2

    1 CBtis 190, Av. 15 Col. Venustiano Carranza 2da Secc, 94297, Boca del Rio Ver., México.

    2 Facultad de Instrumentación Electrónica, Universidad Veracruzana, Circuito Gonzalo Aguirre Beltrán S/N, Xalapa, Veracruz, 91000, México.

    3 Doctorado en Ciencias de la Ingeniería, Tecnológico Nacional de México en Celaya, Antonio García Cubas #600, Celaya, Guanajuato, 38010, México.

    Abstract This work introduces the nonlinearities distribution homotopy perturbation method in order to get an accurate analytical approximate solution for the relevant Bratu problem which has many applications in science and engineering. We will show that the distribution of the nonlinearities in the different iterations of the proposed method is a convenient scheme which eases the obtaining of the approximation based on elementary integrals. We will conclude that the proposed method is very convenient for practical applications.

    Keywords:-Bratu problem, nonlineardifferential equation, approximate solution.


      The Bratu problem in one dimension is called Liouville- Bratu-Gelfand and has many applications in science and engineering where physical and chemical phenomena are modeled. One of the applications of this equation is in solid fuel ignition models of the thermal combustion theory [1-2]. Gelfand's differential equation with boundary value in one dimension in planar coordinates is given by [3-8].

      u eu(x) 0. 0<x<1. (1)

      with boundary conditions (0) = (1) = 0.

      The analytical solution of (1)is given by [3, 4]


      Different approximate methods have been applied to solve the Bratu problem. For example in [9] Chebyshev polynomials was employed in order to obtain an approximate solution for the proposed problem. In [10] the Adomian decomposition method was used to solve the equation. In [11] the B-spline method was proposed to solve the equation. In [12] the Leal- polynomials were proposed to find an accurate analytical approximate solution for (1). In [13] Perturbation Method was employed to get a handy approximate solution and so on.

      This research paper is organized as follows In Section 2 we introduce the idea of Homotopy Perturbation Method (HPM). For Section 3 we briefly introduce Nonlinearities Distribution Homotopy Perturbation Method (NDHPM). Additionally Section 4 presents the application of the proposed method in the search for an analytical approximate solution for Bratu problem. Besides a discussion on the results is presented in Section 5. Finally, a brief conclusion is given in Section 6.


      To show how homotopy perturbation method (HPM) works, first, consider a general nonlinear differential equation as follows

      () () = 0, (5)

      () = 2ln

      cosh (( 2) 2)


      cosh ()


      With the boundary conditions

      (, /) = 0, ,



      where is the solution of



      where is a general differential operator, is a boundary

      operator, () a known analytical function and is the domain boundary for . In general terms, can be divided

      = 2 cosh ( ). (3) into two operators and , where is linear and is


      There are three cases for the solution; if > then the

      nonlinear; so that (5) can be rewritten as

      problem has zero solution, if = there is one solution, and < there are two solutions, where = 3.513830719 is the critical value and satisfies the equation.

      1 = 1 2 sinh () . (4)

      4 4

      () + () () = 0.

      A homotopy can be constructed as [14, 15]

      (, ) = (1 )[() (0)] +

      [() + () ()] = 0, ,



      where is a homotopy parameter, which value varies continuously from 0 to 1, 0 is the first approximation for the

      epy( x) 1 py

      p2 y2


      • p3 y3



      solution of (5) that satisfies the boundary conditions. The

      solution of (8) can be written as power series of


      = 0 + 1 + 22 + .

      After substituting (16) into (15) we obtain

      2 2 3 3

      2 2 3 3

      y p p y p y

      (9) 1

      py 2 6 …

      0, (18)

      Substituting (9) into (8) and equating identical powers of p terms, there can found values for 0, 1, 2, . when tends to (5), we obtain the approximate solution for (1) in the form

      let it be

      = 0 + 1 + 2 + .


      y pn y (x),







      In [16] was introduced nonlinearities distribution homotopy perturbation method (NDHPM) with the purpose to reduce the complexity of solving differential equations in terms of power series.

      To carry out the procedure a homotopy is written in the form (Vazquez-Leal et al. 2012b) [16]

      (, ) = (1 )[() ( )] +

      so that after substituting (19) into (18) and equating identical powers of p terms we obtain



      p0 ) y 0,



      p1 ) y 0,

      2 0

      2 0

      p2 ) y y 0,

      p3 ) y y 2 0,


      [() + (, ) (, )] = 0, .


      and so on.

      3 2 0


      Homotopy function (11) is essentially the same as (8), except because the non-linear operator N and the non homogeneous function , contain embedded the homotopy parameter . The rest of the procedure coincides with that explained in the last Section.

      From [14,16] we have

      The solutions of equations (20) which satisfy the initial conditions of the problem in accordance with the proposed problem are:



      = 0 + 1 + 22

      + .


      p0 ) y




      When tends to 1, it is expected to get an approximate solution for (7) in the form

      = 0 + 1 + 2 + .

      p1 ) y








      p2 ) y






      A2 x4

      Next, we will employ NDHPM in order to find an analytical approximate solution for the Bratu problem.

      p3 ) y




      In accordance with the proposed problem we proposed the following homotopy:

      1 p y p y epy( x) 0, (14)


      and so on.

      Substituting (21) into (13)

      2 3 2 4

      2 3 2 4

      y(x) A K x x Ax A x ,

      2 6 24

      where the constant K is given by


      y p epy( x) 0.


      We note exponential function contain embedded the homotopy parameter .

      K c D E.


      After using Taylor approximation for the exponential function we obtain the following approximation as follows:

      With the end to get a better precision, we will determine the value of K, from the requiring that the approximate soltion satisfy the real value for x=0.5. As a consequence

      (22) adopts the following form:




      A2 x4

      y(x) 384 0.0208A 0.256817 x 2


      6 24


      With the purpose to verify the accuracy of the proposed approximation we propose the value =0.5 in such a way that

      (24) assumes the form




      A2 x4

      y(x) 384 0.0208A 0.256817 x 4 12 48 .


      With the end to determinate A we apply the condition y(1)=0 to (25) to generate the following algebraic equation

      7 A2 3A

      Figure 1. Comparison between numerical solution and approximate solution

      (27) for = 0.5

      On the other hand, Table 1 shows a more detailed comparison between numerical solution for (1) and

      384 48

      0.006817 0.


      approximation (27) of this work for =0.5. The mentioned table shows that the relative error is less than one percent for

      The solution to this equation is substituted into (25) to obtain

      x2 3 4

      most of the interval. It is clear that (27), is competitive taking into account that the proposed method required to solve essentially elementary integrals and an algebraic

      y(x) 0.2590x 0.008817x


      0.0002332x . (27)

      equation. From above it is very important to emphasize that, it is possible to improve the accuracy of our approximation

      Figure 1 compares (27) with the numerical solution for (1) with = 0.5.

      Although the precision of (27) is clear at sight, from Table 1 is clear the precision of the proposed solution.


In this work NDHPM was employed in order to find an approximate solution, for the nonlinear ordinary differential equation that describes the relevant Bratu problem that models the solid fuel ignition which rises of the thermal combustion theory. One of the great advantages of the proposed method is that it systematically distribute the nonlinearity term in the different iterations which ease the obtaining of an analytical approximate solution. The procedure consisted in ensure that the approximate solution satisfy the real value for x=0.5 and pose an algebraic equation for determining an unknown parameter A. The solution of this equation was obtained from the boundary condition (1) = 0 This solution provides the sought analytical approximation solution.

Figure 1 shows the comparison between numerical solution and approximate solution (27) for =0.5. It can be noticed that curves are in good agreement, whereby it is clear that the proposed method is potentially useful in the search for approximate solutions of nonlinear problems definite with boundary conditions.

performing more iterations from the homotopy (18).


This work presented method NDHPM with the purpose to calculate an analytical approximate solution for the Bratu problem. The method basically works calculating elementary integrals after distributing the nonlinear term in the several iterations of the method in order to ease the procedure. As a matter of fact the method expressed the solution of a differential equation in terms of the solution of one or more algebraic equations. The application of the proposed method in this work showed the convenience of employ NDHPM as a practical tool with the purpose to obtain accuracy solutions for boundary value problems instead of using other more sophisticated and cumbersome procedures.




Error relative












0.59 %

























Table 1: Comparison between (27), and exact solution using =0.5


Authors would like to thank Roberto Ruiz Gomez for his contribution to this project.


  1. A.S. Mounim and BM de Dormale, From the fitting techniques to accurate schemes for the liouville-bratu-gelfand problem, Numerical Methods for Partial Differential Equations, 2006, 22(4), 761775.

  2. S. Li and S.J. Liao, An analytic approach to solve multiple solutions of a strongly nonlinear problem, Applied mathematics and computation, 2005, 169(2), 854865.

  3. Jacobsen, J. and Schmitt, K. (2002) The Liouville-Bratu-Gelf and Problem for Radial Operators. Journal of Differential Equations, 184, 283-298.

  4. Ascher, U.M., Matheij, R. and Russell, R.D. (1995) Numerical Solution of Boundary Value Problems for Ordinary Differential Equations. Society for Industrial and Applied Mathematics, Philadelphia.

  5. Boyd, J.P. (2003) Chebyshev Polynomial Expansions for Simultaneous Approximation of Two Branches of a Function with Application to the One-Dimensional Bratu Equation. Applied Mathematics and Computation, 142, 189-200.

  6. Boyd, J.P. (1986) An Analytical and Numerical Study of the Two- Dimensional Bratu Equation. Journal of Scientific Computing, 2, 183- 206.

  7. Buckmire, R. (2003) Investigations of Nonstandard, Mickens-Type, Finite-Difference Schemes for Singular Boundary Value Problems in Cylindrical or Spherical Coordinate. Numerical Methods for Partial Differential Equations, 19, 380-398.

  8. Aregbesola, Y. (2003) Numerical Solution of Bratu Problem Using the Method of Weighted Residual. Electronic Journal of Southern African Mathematical Sciences Association, 3, 1-7.

  9. J.P. Boyd, Chebyshev polynomial expansions for simultaneous approximation of two branches of a function with application to the one-dimensional bratu equation, Applied mathematics and computation, 2003, 143(2-3), 189200.

  10. A.M. Wazwaz, Adomian decomposition method for a reliable treatment of the bratutype equations. Applied Mathematics and Computation, 2005, 166(3), 652663.

  11. H. Caglar, N. Caglar, M. Ä Ozer, A. Valar3stos, and A.N. Anagnostopoulos. B-spline method for solving Bratus problem. International Journal of Computer Mathematics, 87(8): 18851891, 2010.

  12. Vazquez-Leal, Hector, Mario Alberto Sandoval-Hernandez, Uriel Filobello-Nino, and Jesus Huerta-Chua. "The novel Leal-polynomials for the multi-expansive approximation of nonlinear differential equations." Heliyon 6, no. 4 (2020): e03695.

  13. Filobello-Nino, Uriel, Hector Vázquez-Leal, K. Boubaker, Y. Khan, A. Perez-Sesma, A. Sarmiento-Reyes, V. M. Jimenez-Fernandez et al. "Perturbation method as a powerful tool to solve highly nonlinear problems: the case of Gelfand's equation." Asian Journal of Mathematics & Statistics 6, no. 2 (2013): 76

  14. He JH (1998) A coupling method of a homotopy technique and a perturbationtechnique for nonlinear problems. Int J Non-Linear Mech 351:3743,doi:10.1016/S0020-7462(98)00085-7

  15. He JH (1999) Homotopy perturbation technique. Comput Methods Applied Mech Eng 178:257262, doi:10.1016/S0045-7825(99)00018-3

  16. Vazquez-Leal H, Sarmiento-Reyes A, Khan Y, Filobello-Nino U,Diaz- Sanchez A (2012b) Rational biparameter homotopy perturbation method and Laplace-Padé coupled version. J Appl Math 2012(923975):21, doi:10.1155/ 2012/923975

Leave a Reply