An Efficient Decomposition Method for Solving Bratu’s Boundary Value Problem ()
1. Introduction
In 1914, the Bratu problem was first set up by Bratu, then the problem was named after him [1] . Gelfand and the French mathematician Liouville [2] gave the problem’s generalization so-called Liouville-Bratu-Gelfand problem. The summary of the history and significance of the problem is in [3] .
The classical nonlinear Bratu’s boundary value problem in one-dimensional planar coordinates is in the form [3] - [8] .
(1)
The exact solution to Equation (1) is given in [3] [4] [5] by
(2)
where
satisfies
There are three cases for
:
1) If
, then the Bratu problem has zero solution.
2) If
, then the Bratu problem has one solution.
3) If
, then the Bratu problem has two solutions.
where the critical value
satisfies the equation
where
The Bratu problem was gained the attention of researchers for two important reasons:
First: it appears in a wide variety of application areas such as physical, chemical and engineering. Through these applications, we can mention the reaction of thermal, the theory of chemical reactor and nanotechnology. Specific models include Chandrasekhar model of the expansion of the universe [9] , and solid fuel ignition model for the thermal reaction process in a combustible, nondeformable material of constant density during the ignition period [3] . In recent years, the problem has found applications in Engineering such as electro-spinning process for the manufacturing of the Nano-fibers introduced in [10] and [11] .
Second: because of its simplicity, the equation is widely used as a benchmarking tool for various numerical techniques. Several numerical methods, such as the finite difference method [12] , finite element approximation [7] , weighted residual method [8] , and the shooting method [13] , a variational iteration scheme [14] [15] , differential transformation method [16] , homotopy analysis [17] and some additional papers about numerical and analytical methods of the Bratu-type problem [18] [19] [20] .
In 1980, the Adomian decomposition method (ADM) was first introduced and developed by Adomian [21] [22] . The main advantage of this method is that the solution is considered as an infinite series which can be rapidly converging. In addition, many kinds of applications, such as the algebraic, integral, differential, partial and integro-differential equations or system of equations can be solved by using the ADM effectively and efficiently. This is independent of the fact that the equation was linear/non- linear, homogeneous/non-homogeneous, or with constant/variable coefficient.
As far as our information, few of researches have used ADM and its modifications to obtain the approximate solution of the Bratu’s boundary value problem, such as ADM used in [23] [24] , the decomposition method with Green function used in [25] and the Laplace transform decomposition method used in [26] [27] .
In present study, we apply the ADM and restarted Adomian decomposition method (RADM) [28] by modify this approach and propose a new integral operator in [29] for solving boundary value problem of the Bratu-type. Whereas the boundary conditions are used directly to obtain the solution and demonstrate its efficiency on two examples compared to the exact solution and another methods in ADM and its application.
This research consists of an introduction and three four sections as follows: In Section 2, the ADM and the RADM with new integral operator will be presented as it applies to Bratu problem. In Section 3, the Comparison of this method to two test problems is presented, together with the Adomian decomposition method and its modifications. Finally, the conclusion is as follows in Section 4.
2. Description of the Method
In this section, the standard ADM and RADM by using new integral operator are introduced, and the scheme is implemented for the solution of Bratu’s boundary value problem.
2.1. Adomian Decomposition Method
We have two cases to apply the standard ADM with new integral operator for solving boundary value problem of the Bratu-type. For first case, we consider the nonlinear term as
, and we use Taylor series of
in second case.
Case 1. We consider the boundary value problem of Bratu-type with nonlinear term as
(3)
The operator form of Equation (3) can be written as
(4)
where the differential operator
is
The inverse
is assumed a new two-fold integral operator [29] given by
By applying
on both sides of Equation (4) and using the boundary condition
we get
(5)
The decomposition method assumes
and the nonlinear term
as infinite series given by
(6)
where
are the so-called Adomian polynomials of
given by
substituting Equation (6) into Equation (5) yields
Identifying the zeroth component
by 0, the remaining components
can be obtained recurrently by using the relation
(7)
where
of the nonlinear term
, we can obtain
The solution components
from Equation (7) can be calculated as
The approximate solution
is obtained in a series form
Case 2. For accelerating the convergence of the ADM when used for nonlinear differential equations, we replaced nonlinear terms by their Taylor expansion.
Toward this end, we can consider
as
The Equation (3) can be written as
(8)
Applying the standard ADM with new integral operator in Equation (8), we get
where
Identifying the zeroth component
by
, the remaining components
can be obtained recurrently by using the relation
(9)
The solution components
from Equation (9) can be calculated as
The approximate solution
is obtained in a series form
2.2. Restarted Adomian Decomposition Method
Generally, the RADM has the same structure as the ADM, but the ADM is used more than once. In this section, we extend the RADM for solving the Bratu’s boundary value problem. Noting that we can apply the RADM with new integral operator for the translated equation
(10)
Applying the standard ADM with new integral operator in Equation (10), we get
(11)
This gives
(12)
Adding and subtracting
to right side of Equation (11) to obtain
By equating the terms we can get
(13)
Step 1: In this step,
is calculated from Equation (12) as follows:
So
Step 2: Now, components of the RADM is computed from Equation (13) as follows:
So
The approximate solution
is obtained in a series form
3. Numerical Results
In this section, we present our numerical results for Bratu’s problem. In order to demonstrate the robustness of the schemes, we will consider of two parameters for the eigenvalue λ. The two methods will be illustrated through the following examples.
3.1. Bratu’s Problem for
We can consider the Bratu’s problem
In Table 1, shows the exact solution for case
derived from Equation (2) and the absolute errors, i.e. the difference between the exact and the numerical solutions obtained by the ADM, the ADM with Taylor and the RADM, using only five terms.
The approximation’s accuracy is reflected in Figure 1, while the solid curve represents the exact solution, the dotted curve is the ADM’s approximation solution, the dashed curve is the ADM’s approximation solution with Taylor, but the dash dotted curve is the approximation solution of RADM. We can realize the almost perfect match of these solutions.
3.2. Bratu’s Problem for
We can consider the Bratu’s problem
Table 2 shows the exact solution for case
derived from Equation (2) and the absolute errors use only five terms.
The approximation’s accuracy is reflected in Figure 2. We can see the almost
Table 1. Error values of the approximate solutions to the Bratu’s boundary value problem when λ = 1.
Figure 1. Comparison between the exact and approximate solutions for λ = 1.
Table 2. Error values of the approximate solutions to the Bratu’s boundary value problem when λ = 2.
Figure 2. Comparison between exact and approximate solutions for λ = 2.
perfect match of those solutions, bating the approximate solutions of ADM.
4. Conclusion
The ADM and RADM with new integral operator as presented in this paper have been shown to be more efficient for solving Bratu’s boundary value problem. The main advantage of the method is that it can use the boundary conditions directly. In addition, it is capable of greatly reducing the size of computational work while still maintaining high accuracy of the numerical solution. Furthermore, this paper shows the validity of the ADM and RADM with new operator for nonlinear problems in science and engineering.