Adomian decomposition method is presented as a method for the solution of the Burger’s equation, a popular PDE model in the fluid mechanics. The method is computationally simple in application. The approximate solution is obtained by considering only the first two terms of the decomposition in this paper. Numerical experimentation shows accuracy of a minimum error of order five for various space steps and coefficient of kinematic viscosity. The method is considered high in accuracy.
Burger’s equation is a fundamental partial differential equation in fluid mechanics. It is also a very important model encountered in several areas of applied mathematics such as heat conduction, acoustic waves, gas dynamics and traffic flow [1] . Analytical solutions of the partial differential equations modeling physical phenomena exist only in few of the cases. Therefore the need for the construction of efficient numerical methods for the approximate solution of these models always exists. Many of the analytical solutions to the Burger’s equation involve Fourier series. According to [2] , the convergence of such Fourier-series based solutions is very slow. Several researchers have proposed various numerical methods for the solution of the Burger’s equation. [3] and [4] used the finite element method for the solution of the Burger’s equation. [5] and [6] used the finite difference method. [7] used the direct variational method while [8] used the projection method by B-spline.
A decomposition method which provides convergent solutions to nonlinear stochastic operator equations was developed in [9] . [10] gave proof of convergence of Adomian decomposition method when applied to differential equations. [11] constructed an algorithm for solving nonlinear equations based on Newton-Raphson method using Adomian decomposition approach. In an effort to extend the usefulness of Adomian’s decomposition method into other areas of mathematics, [12] developed a technique for calculating Adomian’s polynomials for nonlinear operators. [13] presented a method for the solution of homogeneous differential equations using Adomian’s decomposition without noisy terms. [14] extended the application of the Adomian decomposition method for the calculations of non-differential functions.
In this paper, presentation of a numerical method for the solution of the nonlinear reaction-diffusion Burger’s equation using Adomian decomposition is made. The organization of this paper is as follows: in Section 2 the theoretical approach is presented. In Section 3 the Adomian’s polynomials for the Burger’s equation are determined. In Section 4 computational results for the Burger’s equation using Adomian’s decomposition are presented while conclusion is presented in Section 5.
2. Theoretical Approach
Consider the Burger’s equation
Subject to initial condition
And boundary conditions:
Defining the operators:
Equation (1) can be written as:
Solving Equation (6) for and respectively gives
and
where and are the inverse operators of and given respectively as:
Operating both sides of Equations (7) and (8) with the inverse operators (9) we obtain:
where
and
Adding Equations (10) and (11) and dividing by 2, gives the canonical form:
The parameterized form of (14) is
where
According to [15] the parameterized decomposition forms of and Nu are:
An are the Adomian’s polynomials which can be generated for all types of nonlinearities [16] .
Substituting Equations (17) and (18) into Equation (15) gives
Expanding both sides of Equation (19) gives
By comparing coefficients of both hand sides of Equation (20), it is obtained that:
From which we establish the recursive relation
3. Determination of Adomian’s Polynomials
The Adomian’s polynomials, are defined in such a way that
[17]
Substituting Equation (17) into Equation (18) and expanding the RHS gives
From Equation (22), we establish that the Adomian’s polynomials have the following forms:
The Adomian’s polynomials for Equation (1) are obtained from the recurrent relation
4. Numerical Experiment
We present numerical results to illustrate the effectiveness of the proposed method. Consider Burger’s Equation (1) with the following initial and boundary conditions.
The exact solution of Equation (1) with the above conditions was given in [1] as
Computations of Absolute and Relative Errors
Absolute errors of the method were computed by use of the formula:
where the numerical solution at the grid point is, and the exact solution at the same grid point is.
Relative errors were computed by use of the formula:
where the numerical solution at the grid point is, and the exact solution at the same grid point is.
Figure 1 is a 3 dimensional plot of the surface while Figure 2 is a 3 dimensional plot of the surface
5. Conclusion
In this paper, Adomian decomposition method is used to solve the burger’s equation numerically. From Table 1 and Table 2, this method is considered to give accurate results for specified values of the parameter, ε. We considered the first two terms of the decomposition to make the approximate solution. Accuracy of the method
A 3 dimensional plot of the surface<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/11-1100460x57.png" xlink:type="simple"/></inline-formula>
A 3 dimensional plot of the surface<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/11-1100460x59.png" xlink:type="simple"/></inline-formula>
Approximate solutions by Adomian’s decomposition, errors and relative errors for t = 0.2, e = 0.001, a = 2
x
Adomian
Theoretical solution
Error
Relative Errors
0.1
0.0006610924400
0.0003891849100
2.719075297 × 10−4
2.718017485 × 10−4
0.2
0.0013184569550
0.0008129882596
5.054686954 × 10−4
5.050580892 × 10−4
0.3
0.001968724990
0.0013210998120
6.47625178 × 10−4
6.467707292 ×10−4
0.4
0.002596743099
0.0020106848070
5.86058292 × 10−4
5.84882278 × 10−4
0.5
0.0031486447293
0.0031366545250
1.9922768 × 10−5
1.986047255 × 10−5
0.6
0.003546100003
0.005799752849
2.253652846 × 10−3
2.240657585 × 10−3
0.7
0.003616277199
0.03291862860
2.930235148 × 10−2
2.836849937 × 10−2
0.8
0.003425426033
−0.006782588856
1.020801489 × 10−2
1.027772447 × 10−2
Approximate solutions by Adomian’s decomposition, errors and relative errors for t = 0.05, e = 0.01, a = 2
x
Adomian
Theoretical Solution
Error
Relative Errors
0.2
0.01342133395
0.008119261273
5.302072677 × 10−3
5.259370475 × 10−3
0.4
0.02618512105
0.02006671613
6.11840492 × 10−3
5.998043877 × 10−3
0.6
0.03577004231
0.05766474009
2.189469778 × 10−2
2.070098109 × 10−2
0.8
0.03132012403
−0.06901554553
1.003356696 × 10−1
1.077737325 × 10−1
can be improved by considering more terms in the solution approximation. All computations were carried out using Maple 15.
Cite this paper
Iyakino P.Akpan, (2015) Adomian Decomposition Approach to the Solution of the Burger’s Equation. American Journal of Computational Mathematics,05,329-335. doi: 10.4236/ajcm.2015.53030
ReferencesMaleknejad, K., Babolian, E., Shaerlar, A.J. and Jahnagir, M. (2011) Operational Matrices for Solving Burgers Equation by Using Block-Pulse Functions with Error Analysis. Australian Journal of Basic and Applied Sciences, 5, 602-609.Aksan, E.N., Ozdes, A. and Ozia, T. (2006) A Numerical Solution of Burger’s Based Squares Approximation. Applied Mathematics and Computation, 176, 270-279. http://dx.doi.org/10.1016/j.amc.2005.09.045Varoglu, E. and Finn, W.D.L. (1980) Space-Time Finite Elements Incorporating Characteristics for Burger’s Equation. International Journal for Numerical Methods in Engineering, 16, 171-184. http://dx.doi.org/10.1002/nme.1620160112Caldwell, J., Wanless, P. and Cook, A.E. (1981) A Finite Element Approach to Burger’s Equation. Applied Mathematical Modelling, 5, 189-193. http://dx.doi.org/10.1016/0307-904X(81)90043-3Smith, G.D. (1989) Numerical Solution of Partial Differential Equations. Oxford Press, Oxford.Ames, W.F. (1992) Numerical Methods for Partial Differential Equations. Academic Press, New York.Ozas, T. and Ozdes, A. (1996) A Direct Variational Method Applied to Burgers Equation. Journal of Computational and Applied Mathematics, 71, 163-175. http://dx.doi.org/10.1016/0377-0427(95)00221-9Dag, J., Saka, B. and Boz, A. (2005) B-Splin Galekin Method for Numerical Solutions of the Burgers Equation. Applied Mathematics and Computation, 166, 506-522. http://dx.doi.org/10.1016/j.amc.2004.06.078Adomian, G. (1986) Nonlinear Stochastic Operator Equations. Academic Press, New York.Abbaoui, K. and Cherruault, Y. (1994) Convergence of Adomian’s Method Applied to Differential Equations. Computers & Mathematics with Applications, 28, 103-109. http://dx.doi.org/10.1016/0898-1221(94)00144-8Abbasbandy, S. (2003) Improving Newton-Raphson Method for Nonlinear Equations by Modified Adomian Decomposition Method. Applied Mathematics and Computation, 145, 887-893.
http://dx.doi.org/10.1016/S0096-3003(03)00282-0Chen, W. and Lu, Z. (2004) An Algorithm for Adomian Decomposition Method. Applied Mathematics and Computation, 170, 570-583. http://dx.doi.org/10.1016/j.amc.2003.10.037Gharsseldien, Z.M. and Hemida, K. (2009) New Technique to Avoid Noise Terms on the Solution of Inhomogeneous Differential Equations by Using Adomian Decomposition Method. Communication in Nonlinear Science and Numerical Simulation, 14, 685-696. http://dx.doi.org/10.1016/j.cnsns.2007.11.018Wu, G.C. (2011) Adomian Decomposition Method for Non-Smooth Initial Value Problems. Mathematical and Computer Modeling, 54, 2104-2108. http://dx.doi.org/10.1016/j.mcm.2011.05.018Mamaloukas, C. (2000) An Approximate Solution of Burgers Equation Using Adomian’s Decomposition Method. Proceedings of the Conference of Applied Differential Geometry—General Relativity and the Workshop on Global Analysis, Differential Geometry and Lie Algebra, Balkan Press, 88-98.Adomian, G. (1989) Nonlinear Stochastic Systems Theory and Application to Physics. Kluwer Academic Publishers, Netherlands. http://dx.doi.org/10.1007/978-94-009-2569-4Bhadauria, R., Singh, A.K. and Singh, D.P. (2012) Solution of Reaction-Diffusion Equation by Adomian Decomposition Method. International Journal of Engineering Science and Technology, 4, 2626-2630.