An Efficient Decomposition Method for Solving Bratu’s Boundary Value Problem

Abstract

The purpose of this paper is to employ the Adomian Decomposition Method (ADM) and Restarted Adomian Decomposition Method (RADM) with new useful techniques to resolve Bratu’s boundary value problem by using a new integral operator. The solutions obtained in this way require the use of the boundary conditions directly. The obtained results indicate that the new techniques give more suitable and accurate solutions for the Bratu-type problem, compared with those for the ADM and its modification.

Share and Cite:

Al-Mazmumy, M. , Al-Mutairi, A. and Al-Zahrani, K. (2017) An Efficient Decomposition Method for Solving Bratu’s Boundary Value Problem. American Journal of Computational Mathematics, 7, 84-93. doi: 10.4236/ajcm.2017.71007.

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] .

$\begin{array}{l}{u}^{″}+\lambda {\text{e}}^{u}=0,\text{}0 (1)

The exact solution to Equation (1) is given in [3] [4] [5] by

$u\left(x\right)=-2\mathrm{ln}\left[\frac{\mathrm{cosh}\left(\left(x-\frac{1}{2}\right)\frac{\theta }{2}\right)}{\mathrm{cosh}\left(\frac{\theta }{4}\right)}\right].$ (2)

where $\theta$ satisfies

$\theta =\sqrt{2\lambda }\mathrm{cosh}\left(\frac{\theta }{4}\right).$

There are three cases for $\lambda$ :

1) If $\lambda >{\lambda }_{c}$ , then the Bratu problem has zero solution.

2) If $\lambda ={\lambda }_{c}$ , then the Bratu problem has one solution.

3) If $\lambda <{\lambda }_{c}$ , then the Bratu problem has two solutions.

where the critical value ${\lambda }_{c}$ satisfies the equation

$1=\frac{1}{4}\sqrt{2{\lambda }_{c}}\mathrm{sinh}\left(\frac{{\theta }_{c}}{4}\right),$

where

${\lambda }_{c}=3.513830719.$

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.

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 ${\text{e}}^{u}$ , and we use Taylor series of ${\text{e}}^{u}$ in second case.

Case 1. We consider the boundary value problem of Bratu-type with nonlinear term as ${\text{e}}^{u}$

$\begin{array}{l}{u}^{″}+\lambda {\text{e}}^{u}=0,\text{}0 (3)

The operator form of Equation (3) can be written as

$\begin{array}{l}Lu=-\lambda {\text{e}}^{u},\hfill \\ u\left(0\right)=u\left(1\right)=0.\hfill \end{array}$ (4)

where the differential operator $L$ is $L=\frac{{\text{d}}^{2}}{\text{d}{x}^{2}}.$

The inverse ${L}^{-1}$ is assumed a new two-fold integral operator [29] given by

${L}^{-1}\left(\cdot \right)={\int }_{0}^{x}{\int }_{0}^{x}\left(\cdot \right)\text{d}x\text{d}x-x{\int }_{0}^{1}{\int }_{0}^{x}\left(\cdot \right)\text{d}x\text{d}x.$

By applying ${L}^{-1}$ on both sides of Equation (4) and using the boundary condition $u\left(0\right)=u\left(1\right)=0$ we get

$u\left(x\right)=-\lambda {L}^{-1}\left({\text{e}}^{u}\right).$ (5)

The decomposition method assumes $u\left(x\right)$ and the nonlinear term

$N\left(u\right)={\text{e}}^{u}$ as infinite series given by

$u\left(x\right)=\underset{n=0}{\overset{\infty }{\sum }}{u}_{n},N\left(u\right)=\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}$ (6)

where ${A}_{n}$ are the so-called Adomian polynomials of ${u}_{0},\cdots ,{u}_{n}$ given by

${A}_{n}=\frac{1}{n!}\frac{{\text{d}}^{n}}{\text{d}{\lambda }^{n}}{\left[N\left(\underset{i=0}{\overset{n}{\sum }}{\lambda }^{i}{u}_{i}\right)\right]}_{\lambda =0},n=0,1,2,\cdots$

substituting Equation (6) into Equation (5) yields

$\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}=-\lambda {L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}\right).$

Identifying the zeroth component ${u}_{0}$ by 0, the remaining components

${u}_{n},n\ge 1$ can be obtained recurrently by using the relation

$\begin{array}{l}{u}_{0}=0,\\ {u}_{n+1}=-\lambda {L}^{-1}\left({A}_{n}\right),n\ge 0,\end{array}$ (7)

where ${A}_{n}$ of the nonlinear term ${\text{e}}^{u}$ , we can obtain

$\begin{array}{l}{A}_{0}={\text{e}}^{{u}_{0}},\\ {A}_{1}={u}_{1}{\text{e}}^{{u}_{0}},\\ {A}_{2}=\left({u}_{2}+\frac{1}{2}{u}_{1}^{2}\right){\text{e}}^{{u}_{0}},\\ {A}_{3}=\left({u}_{3}+{u}_{1}{u}_{2}+\frac{1}{6}{u}_{1}^{3}\right){\text{e}}^{{u}_{0}},\\ ⋮\end{array}$

The solution components ${u}_{n+1}$ from Equation (7) can be calculated as

$\begin{array}{l}{u}_{1}=-\frac{1}{2}\lambda {x}^{2}+\frac{1}{2}\lambda x,\\ {u}_{2}=\frac{1}{24}{\lambda }^{2}{x}^{4}-\frac{1}{12}{\lambda }^{2}{x}^{3}+\cdots ,\\ {u}_{3}=-\frac{1}{180}{\lambda }^{3}{x}^{6}+\frac{1}{60}{\lambda }^{3}{x}^{5}-\cdots ,\\ {u}_{4}=\frac{17}{20160}{\lambda }^{4}{x}^{8}-\frac{17}{5040}{\lambda }^{4}{x}^{7}+\cdots ,\\ ⋮\end{array}$

The approximate solution $u\left(x\right)$ is obtained in a series form

$u\left(x\right)=\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}={u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+{u}_{4}+\cdots .$

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 ${\text{e}}^{u}$ as

${\text{e}}^{u}\simeq 1+u+\frac{{u}^{2}}{2!}+\frac{{u}^{3}}{3!}.$

The Equation (3) can be written as

$\begin{array}{l}{u}^{″}+\lambda \left(1+u+\frac{{u}^{2}}{2!}+\frac{{u}^{3}}{3!}\right)=0,0 (8)

Applying the standard ADM with new integral operator in Equation (8), we get

$\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}=-\lambda \left({L}^{-1}\left(1\right)+{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}\right)+\frac{1}{2!}{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}\right)+\frac{1}{3!}{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{B}_{n}\right)\right),$

where

${A}_{n}=\frac{1}{n!}\frac{{\text{d}}^{n}}{\text{d}{\lambda }^{n}}{\left[{\left(\underset{i=0}{\overset{n}{\sum }}{\lambda }^{i}{u}_{i}\right)}^{2}\right]}_{\lambda =0},{B}_{n}=\frac{1}{n!}\frac{{\text{d}}^{n}}{\text{d}{\lambda }^{n}}{\left[{\left(\underset{i=0}{\overset{n}{\sum }}{\lambda }^{i}{u}_{i}\right)}^{3}\right]}_{\lambda =0},n=0,1,2,\cdots$

Identifying the zeroth component ${u}_{0}$ by $–\lambda {L}^{-1}\left(1\right)$ , the remaining components ${u}_{n},n\ge 1$ can be obtained recurrently by using the relation

$\begin{array}{l}{u}_{0}=-\frac{1}{2}\lambda {x}^{2}+\lambda \frac{1}{2}x,\hfill \\ {u}_{n+1}=-\lambda \left({L}^{-1}\left({u}_{n}\right)+\frac{1}{2}{L}^{-1}\left({A}_{n}\right)+\frac{1}{6}{L}^{-1}\left({B}_{n}\right)\right),n\ge 0.\hfill \end{array}$ (9)

The solution components ${u}_{n+1}$ from Equation (9) can be calculated as

$\begin{array}{l}{u}_{1}=\frac{1}{2688}{\lambda }^{4}{x}^{8}-\frac{1}{672}{\lambda }^{4}{x}^{7}+\cdots ,\\ {u}_{2}=-\frac{1}{3913728}{\lambda }^{7}{x}^{14}+\frac{1}{559104}{\lambda }^{7}{x}^{13}-\cdots ,\\ {u}_{3}=\frac{5}{28554559488}{\lambda }^{10}{x}^{20}-\frac{25}{14277279744}{\lambda }^{10}{x}^{19}+\cdots ,\\ {u}_{4}=-\frac{449}{3741789475307520}{\lambda }^{13}{x}^{26}+\frac{449}{287829959639040}{\lambda }^{13}{x}^{25}-\cdots ,\\ ⋮\end{array}$

The approximate solution $u\left(x\right)$ is obtained in a series form

$u\left(x\right)=\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}={u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+{u}_{4}+\cdots .$

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

$\begin{array}{l}{u}^{″}+\lambda \left(1+u+\frac{{u}^{2}}{2!}+\frac{{u}^{3}}{3!}\right)=0,\text{}0 (10)

Applying the standard ADM with new integral operator in Equation (10), we get

$\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}=-\lambda \left({L}^{-1}\left(1\right)+{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}\right)+\frac{1}{2!}{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}\right)+\frac{1}{3!}{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{B}_{n}\right)\right).$ (11)

This gives

$\begin{array}{l}{u}_{0}=-\frac{1}{2}\lambda {x}^{2}+\lambda \frac{1}{2}x,\\ {u}_{n+1}=-\lambda \left({L}^{-1}\left({u}_{n}\right)+\frac{1}{2}{L}^{-1}\left({A}_{n}\right)+\frac{1}{6}{L}^{-1}\left({B}_{n}\right)\right),n\ge 0.\end{array}$ (12)

Adding and subtracting $g\left(x\right)$ to right side of Equation (11) to obtain

$\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}=g\left(x\right)-\lambda \left({L}^{-1}\left(1\right)+{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{u}_{n}\right)+\frac{1}{2!}{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}\right)$ $+\frac{1}{3!}{L}^{-1}\left(\underset{n=0}{\overset{\infty }{\sum }}{B}_{n}\right)\right)-g\left(x\right).$

By equating the terms we can get

$\begin{array}{l}{u}_{0}=g\left(x\right),\\ {u}_{1}=-\lambda \left({L}^{-1}\left(1\right)+{L}^{-1}\left({u}_{0}\right)+\frac{1}{2}{L}^{-1}\left({A}_{0}\right)+\frac{1}{6}{L}^{-1}\left({B}_{0}\right)\right)-g\left(x\right),\\ {u}_{n+2}=-\lambda \left({L}^{-1}\left({u}_{n+1}\right)+\frac{1}{2}{L}^{-1}\left({A}_{n+1}\right)+\frac{1}{6}{L}^{-1}\left({B}_{n+1}\right)\right),n\ge 0.\end{array}$ (13)

Step 1: In this step, $g\left(x\right)$ is calculated from Equation (12) as follows:

$\begin{array}{l}{u}_{0}=-\frac{1}{2}\lambda {x}^{2}+\frac{1}{2}\lambda x,\\ {u}_{1}=\frac{1}{2688}{\lambda }^{4}{x}^{8}-\frac{1}{672}{\lambda }^{4}{x}^{7}+\cdots ,\end{array}$

So

$g\left(x\right)={\varphi }^{1}\left(x\right)={u}_{0}+{u}_{1}.$

Step 2: Now, components of the RADM is computed from Equation (13) as follows:

$\begin{array}{l}{u}_{0}=g\left(x\right),\\ {u}_{1}=-\frac{1}{75744726220800}{\lambda }^{13}{x}^{26}+\frac{1}{5826517401600}{\lambda }^{13}{x}^{25}-\cdots ,\\ {u}_{2}=\frac{1}{2070913964648082230476800}{\lambda }^{22}{x}^{44}-\frac{1}{94132452938549192294400}{\lambda }^{22}{x}^{43}+\cdots ,\\ {u}_{3}=-\frac{1921}{110351113887194570160875411425198080000}{\lambda }^{31}{x}^{62}+\cdots ,\\ {u}_{4}=\frac{18470161}{29478690975630536257743382096194416238070333440000000}{\lambda }^{40}{x}^{80}-\cdots ,\end{array}$

So

${\varphi }^{2}\left(x\right)={u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+{u}_{4}.$

The approximate solution $u\left(x\right)$ is obtained in a series form

$u\left(x\right)={\varphi }^{2}\left(x\right)={u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+{u}_{4}.$

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 $\lambda =1$

We can consider the Bratu’s problem

$\begin{array}{l}{u}^{″}+{\text{e}}^{u}=0,0

In Table 1, shows the exact solution for case $\lambda =1$ 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 $\lambda =2$

We can consider the Bratu’s problem

$\begin{array}{l}{u}^{″}+{\text{e}}^{u}=0,\text{}0

Table 2 shows the exact solution for case $\lambda =2$ 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.

Conflicts of Interest

The authors declare no conflicts of interest.