Numerical Solutions of Three-Dimensional Coupled Burgers’ Equations by Using Some Numerical Methods

Abstract

In this paper, we found the numerical solution of three-dimensional coupled Burgers’ Equations by using more efficient methods: Laplace Adomian decomposition method, Laplace transform homotopy perturbation method, variational iteration method, variational iteration decomposition method and variational iteration homotopy perturbation method. Example is examined to validate the efficiency and accuracy of these methods and they reduce the size of computation without the restrictive assumption to handle nonlinear terms and it gives the solutions rapidly.

Share and Cite:

Ahmad Alhendi, F. and Alderremy, A. (2016) Numerical Solutions of Three-Dimensional Coupled Burgers’ Equations by Using Some Numerical Methods. Journal of Applied Mathematics and Physics, 4, 2011-2030. doi: 10.4236/jamp.2016.411201.

1. Introduction

The Burgers Equation was first presented by Bateman  and treated later by J. M. Burgers (1895-1981) then it is widely named as Burgers’ Equation  . Burgers’ Equation is nonlinear partial differential equation of second order which is used in various fields of physical phenomena such as boundary layer behaviour, shock weave formation, turbulence, the weather problem, mass transport, traffic flow and acoustic transmission   . In addition, coupled Burgers’ Equations has played an important role in many physical applications such as hydrodynamic turblence, vorticity transport, skock wave, dispersion in porous media and wave processes. In particular, the three-dimentional coupled Burgers’ Equations are important in large scale structure formation in the universe   . In order to make a great application for burgers’ Equations, many researchers have been interested in solving it by various techniques. Analytic solution of one-dimensional Burgers’ Equation is got by many standard methods such as Backland transformation method, differential transformation method and tanh-coth method  , while an analytical solution of two-dimensional coupled Burgers’ Equations is first presented by Fletcher using the Hopf-Cole transformation  and an analytical solution of three-dimensional coupled Burgers’ Equations is derived by Srivastava et al using the variable of separable and Hopf-Cole transformation  . The finite difference, finite element, spectral methods, Adomian decomposition method, the variational iteration method, homotopy perturbation method HPM and Eulerian-Lagrangian method, etc. gave an numerical solution of one- and two-dimentional Burgers’ Equations   -  .

Shukla et al are proposed a numerical solutions of three-dimensional coupled viscous Burgers’ Equations by using a modified cubic B-spline differential quadrature method  .

The motive of this paper is to find the numerical solution of three-dimensional coupled Burgers’ Equations by using more efficiently methods: Laplace Adomian decomposition method, Laplace transform homotopy perturbation method, the variational iteration method, variational iteration decomposition method and variational iteration homotopy perturbation method. We consider three-dimensional couple Burgers’ Equations as the following  : (1)

with the initial conditions: (2)

and the boundary conditions: (3)

where and is its boundary, and are the velocity components to be determined, g1, g2, g3, f1, f2 and f3 are known functions and R is the Reynolds number.

This paper is organized into five sections. Each method is in one section. We showed an overview of these methods then explained methodology and finally illustrated the methods by using examples. It is clear to see that numerical methods are reasonably in good covenant with the exact solution.

The Laplace transform LT is an integral transform discovered by Pierre-Simon Laplace. LT is a very powerful technique for solving ordinary and partial differential Equations, which transforms the original differential equation into an elementary algebraic expression  .

Definition: the Laplace transform of where , denoted by , is given by: Adomian decomposition method ADM is proposed in 1980 by Grorge Adomian. ADM has been encountered much attention in recent years in applied mathematics in general and in solving Burgers’ Equations in particular. A wide class of linear and nonlinear, ordinary and partial differential Equations solved easily and more accurately via ADM. It has successfully used to handle most type of physical models of partial differential Equations without dependence on linearization or any restrictive assumptions that may change physical behavior of the models under study     .

ADM consists of decomposing the unknown functions of any equations into a sums of an infinite number of components defined by the decomposition series: (4)

where are to be determined in recursive manner. The nonlinear term can be expressed by an infinite series of Adomian polynomial which is given as: (5)

where Adomian polynomial using the form (6)

There is a growing interest of researchers has been to study the Adomian decomposition method ADM      . In this work, we will use Laplace transform- Adomian decomposition method (LT-ADM) introduced by Khuri khuri  . Some time it is known as Laplace Adomian decomposition method (LADM). This numerical technique explains how the Laplace transform may be used to approximate the solutions of the nonlinear partial differential equations (PDEs) including Burgers’ Equations with the decomposition method     .

2.1. Methodology of LT-ADM for Three-Dimensional Couple Burgers’ Equations

We consider the system (1) and apply LT on both side of it:

(7)

We can write (7) as:

(8)

By using (2) we get:

(9)

Using inverse Laplace transform on both sides of (9), we have

(10)

from (4) we can write the solutions as:

(11)

Now, we assume the nonlinear terms as:

(12)

where and are the Adomian polynomial given as (6). From (11) and (12) into (10), we have:

(13)

where. Then, using (2), (6), (12) into (13), we have

(14)

2.2. Example

Consider (1) if the exact solution is  :

(15)

To solve this example by LT-ADM, we follow the methodology which discussed in subsection (2.1). The accuracy of LT-ADM for the three-dimensional coupled Burgers’ Equations agrees very well with the exact solution and absolute errors are very small for the current choice of x, y, z and t. The results are shown in Tables 1-3 for R = 100, x = 0.1, y = 0.02 and z = 0.03.

3. Laplace Transform Homotopy Perturbation Method

Homotopy perturbation method HPM was first proposed by He. HPM is combination of traditional perturbation method and homotopy method. The important advantage of HPM is that the nonlinear term can be easily handled. It is easy to calculate the solution

Table 1. The absolute error (AEs) of by LT-ADM for example 2.2.

Table 2. The AEs of by LT-ADM for example 2.2.

Table 3. The AEs of by LT-ADM for example 2.2.

with this method. Linear or nonlinear ODEs and PDEs are studied successfully by using LT-HPM  ,      . To figure out how HPM works  , consider n-dimensional Burgers’ equation

(16)

We construct the following homotopy:

(17)

or

(18)

where, , , is an embedding

parameter while are initial approximations of (16). Assume the solution of (16) has the form

(19)

Now, substituting from Equation (19) in Equation (18) and comparing coefficients of terms with identical powers of p, we get:

(20)

The solution is

(21)

In  , Aminikhah presented LT-HPM to solve nonlinear Blasius' viscous flow equation. In  , some application examples of LT-HPM for nonlinear ODEs with Dirichlet, mixed, and Neumann boundary conditions were presented. It has used to solve PDEs then to solve Burgers’ Equations      .

Now, we are going to study LT-HPM for (1) as following.

3.1. Methodology of LT-HPM for Three-Dimensional Couple Burgers’ Equations

We consider the system (1) and apply HPM on it. We construct the following homo- topy:

(22)

where and are the initial approximation values of (1), and and have the following forms:

(23)

we can write (22) as:

(24)

Applying LT on both side of (24)

(25)

By applying inverse of LT on (25), we get

(26)

By substituting U, V and W from (23) in (26) and comparing coefficients of terms with identical powers of p:

(27)

Finally, the approximate solutions are:

(28)

3.2. Example

To solve (1) by LT-HPM, we follow the methodology which discussed in subsection (3.1). The results in Tables 4-6 show that LT-HPM is more effective and high accuracy when compared with the exact solutions for R = 100, x = 0.1, y = 0.02 and z = 0.03.

Table 4. The AEs of by LT-HPM for example 3.2.

Table 5. The AEs of by LT-HPM for example 3.2.

Table 6. The AEs of by LT-HPM for example 3.2.

4. The Variational Iteration Method

The variational iteration method (VIM) was proposed by Ji-Huan He in 1997     . The VIM has been applied successfully for the most important problems in physically important phenomena including Burgers’ Equation      The VIM solve linear or nonlinear ODEs and PDEs without needing small parameter or linearization and by few iterations lead to high accurate solutions.

To explain the basic concepts of the VIM, we consider n-dim of burgers’ Equation (16). We can write the correction functional for it as  :

(29)

where, , are the general Lagrangian multipliers which can be find via variational theory, are restricted variation which means. The solution is given by

(30)

4.1. Methodology of VIM for Three-Dimensional Couple Burgers’ Equations

Consider the system (1), we can construct a the following correction functional:

(31)

. Then, the iteration formulae are given as:

(32)

4.2. Example

To solve (1) by VIM, we follow the methodology which discussed in subsection (4.1). The results in Tables 7-9 are shown that the efficiency and accuracy of the VIM. It reduces the size of computation without the restrictive assumption to handle nonlinear terms and it gives the solutions rapidly.

5. Variational Iteration Decomposition Method

The variational iteration decomposition method (VIDM) is technique combination of two the most powerful mathematical methods for solving a large class of differential Equations, namely variational iteration method and Adomian decomposition method. In 2007 VIDM has been used to solve quadratic Riccati differential Equation problems

Table 7. The AEs of by VIM for example 4.2.

Table 8. The AEs of by VIM for example 4.2.

Table 9. The AEs of by VIM for example 4.2.

 . Noor et al.   used this method for solving eighth-order boundary value problems, sixth-order boundary value problems and higher dimensional initial boundary value problems. Grover and Tomer solved twelfth order boundary value problems by using VIDM. In 2013, the fractional Riccati differential Equation is solved by VIM by using Adomian polynomials for nonlinear terms  .

To illustrate the general concept of VIDM by using (5), (6) in (29), hence, we have the correction functional for (16) as:

(33)

We are solved three-dimensional coupled Burgers’ Equations by using VIDM as following.

5.1. Methodology of VIDM for Three-Dimensional Couple Burgers’ Equations

We consider the system (1). Next, by using (11) and (12) in (32), we obtain the iterative scheme to find the approximate solutions by VIDM as following:

(34)

5.2. Example

To solve (1) by VIDM by using (34). The numerical results in Tables 10-12 show that VIDM is an effective and powerful method to find better results.

Table 10. The AEs of by VIDM for example 5.2.

Table 11. The AEs of by VIDM for example 5.2.

Table 12. The AEs of by VIDM for example 5.2.

6. Variational Iteration Homotopy Perturbation Method

The variational iteration homotopy perturbation method (VIHPM) is combination of two well-known methods, namely variational iteration method and homotopy perturbation method. VIHPM has been applied in   -  for solving a large class of differential Equations.

To illustrate the concept of VIHPM  , we consider (16) and assume the solution of (16) has the form

(35)

from (35), (16) can be written as:

(36)

we can from the correction functional for(36) we can write

(37)

where, , from (35) in (37) and by comparing the coefficients of like powers of p, we get

(38)

The approximate solutions are give by

(39)

we used this method to solve three-dimensional coupled Burgers’ Equations as fol- lowed.

6.1. Methodology of VIHPM for Three-Dimensional Couple Burgers’ Equations

we consider the correction functional (31) with by assuming that

(40)

by VIHPM, we have

(41)

(42)

by comparing the coefficients of like power of p, we get

(43)

The approximate solutions are given by

(44)

6.2. Example

To solve (1) by VIHPM, we follow the methodology discussed in subsection (6.1). The accuracy of VIHPM for the three-dimensional coupled Burgers’ Equations agrees very well with the exact solution and absolute errors are very small for the current choice of x, y, z and t. The result are shown in Tables 13-15 for R = 100, x = 0.1, y = 0.02 and z = 0.03.

7. Conclusion

In this work, these previous methods mentioned above have been successfully used for finding the solution of three-dimensional coupled Burgers’ Equations. The numerical

Table 13. The AEs of by VIHPM for example 6.2.

Table 14. The AEs of by VIHPM for example 6.2.

Table 15. The AEs of by VIHPM for example 6.2.

results are obtained for approximation and compared with the exact solutions and the results show that we achieve an excellent approximation to the actual solution of the equations by using only two iterations. The results show that these methods are powerful mathematical tools to solving a three-dimensional coupled Burgers’ Equation. In our work, we use the Maple to calculate approximate solutions in our systems by using those very efficient methods.

Acknowledgements

This paper was funded by King Abdulaziz City for Science and Technology (KACST) in Saudi Arabia. The authors therefore thank them for their full collaboration.

Conflicts of Interest

The authors declare no conflicts of interest.     customer@scirp.org +86 18163351462(WhatsApp) 1655362766  Paper Publishing WeChat 