Abstract

The homotopy analysis method (HAM) is applied to obtain the approximate analytic solution of the Korteweg-de Vries (KdV) and Burgers equations. The homotopy analysis method (HAM) is an analytic technique which provides us with a new way to obtain series solutions of such nonlinear problems. HAM contains the auxiliary parameter , which provides us with a straightforward way to adjust and control the convergence region of the series solution. The resulted HAM solution at 8th-order and 14th-order approximation is then compared with that of the exact soliton solutions of KdV and Burgers equations, respectively, and shown to be in excellent agreement.

1. Introduction

It is difficult to solve nonlinear problems, especially by analytic technique. The homotopy analysis method (HAM) [1, 2] is an analytic technique for nonlinear problems, which was first introduced by Liao in 1992. This method has been successfully applied to many nonlinear problems in engineering and science, such as the magnetohydrodynamics flows of non-Newtonian fluids over a stretching sheet [3], boundary layer flows over an impermeable stretched plate [4], nonlinear model of combined convective and radiative cooling of a spherical body [5], exponentially decaying boundary layers [6], and unsteady boundary layer flows over a stretching flat plate [7]. Thus the validity, effectiveness, and flexibility of the HAM are verified via all of these successful applications. Also, many types of nonlinear problems were solved with HAM by others [822].

The Korteweg-de Vries equation (KdV equation) describes the theory of water wave in shallow channels, such as canal. It is an important mathematical model in nonlinear wave's theory and nonlinear optics. The same examples are widely used in solid-state physics, fluid physics, plasma physics, and quantum field theory.

The Burgers equation is a fundamental partial differential equation from fluid mechanics. It occurs in various areas of applied mathematics, such as modeling of gas dynamics and traffic flow. The first steady-state solution of Burgers equation was given by Bateman [23] in 1915. Although, the equation gets its name from the immense research of Burgers [24] beginning in 1939. The study of the general properties of the Burgers equation can be used as a model for any nonlinear wave diffusion problem subject to destruction [25]. Depending on the problem being modeled, this destruction may result from elasticity, gas dynamics, heat conduction, chemical reaction, or other resource.

In this paper, we employ the homotopy analysis method to obtain the solutions of the Korteweg-de Vries (KdV) and Burgers equations so as to provide us a new analytic approach for nonlinear problems.

2. Basic Ideas of Homotopy Analysis Method (HAM)

Consider a nonlinear equation in a general form: where is a nonlinear operator, is unknown function. Let denote an initial guess of the exact solution , an auxiliary parameter an auxiliary function, and an auxiliary linear operator, as an embedding parameter by means of homotopy analysis method, we construct the so-called zeroth-order deformation equation It is very significant that one has great freedom to choose auxiliary objects in HAM. Clearly, when it holds that respectively. Then as long as increase from 0 to 1, the solution varies from initial guess to the exact solution .

Liao [2] by Taylor theorem expanded in a power series of as follow: where

The convergence of the series (2.4) depends upon the auxiliary parameter , auxiliary function , initial guess , and auxiliary linear operator . If they were chosen properly, the series (2.4) is convergence at one has According to definition (2.5), the governing equation can be inferred from the zeroth-order deformation equation (2.2). Define the vector

Differentiating the zero-order deformation equation (2.2) -times with respect to and dividing them by and finally setting we obtain the so-called -order deformation equation where

Theorem 2.1 (Liao [2]). As long as the series (2.6) is convergent, it is convergent to exact solution of (2.1).
Note that homotopy analysis method contains the auxiliary parameter , which provide us with that control and adjustment of the convergence of the series solution (2.6).

3. Exact Solution

The Korteweg-de Vries equation (KdV equation) describes the theory of water wave in shallow channels, such as canal. It is a nonlinear equation which governed by subject to We will suppose that the solution with its derivative, tends to zero [26, 27] when .

In 2001, Wazwaz [28] provided an exact solution or equivalently The Burgers equation is describe by subject to The exact solution of this equation is [29]

4. HAM Solution

4.1. The KdV Equation

For HAM solution of KdV equation we choose as the initial guess and as the auxiliary linear operator satisfying where is a constant.

We consider auxiliary function zeroth-order deformation problemth-order deformation problem

We can use software Mathematica for solving the set of linear equation (4.6) with condition (4.7). It is found that the solution in a series form is given by

The analytical solution given by (4.8) contains the auxiliary parameter , which influences the convergence region and rate of approximation for the HAM solution. In Figure 1, the -curves are plotted for when at 8th-order approximation.

As pointed out by Liao [2], the valid region of is a horizontal line segment. It is clear that the valid region for this case is . According to Theorem 2.1, the solution series (4.8) must be exact solution, as long as it is convergent. In this case, for and , the exact solution and HAM solution are the same, as shown in Figure 2. The obtained numerical results are summarized in Table 1.

In Figure 3, we study the diagrams of the results obtained by HAM for , and in comparison with the exact solution (3.1); we can see the best value for in this case is .

4.2. The Burgers Equation

In this section, for HAM solution of the Burgers equation we choose as the initial guess and as the auxiliary linear operator satisfying where is a constant.

We consider auxiliary function zeroth-order deformation problemth-order deformation problem

We can use software Mathematica for solving the set of linear equation (4.14) with condition (4.15). It is found that the solution in a series form is given by

The analytical solution given by (4.16) contains the auxiliary parameter , which influences the convergence region and the rate of approximation for the HAM solution. In Figure 4, the -curve is plotted for , when at 14th-order approximation.

It is clear that the valid region for this case is . According to Theorem 2.1, the solution series (4.16) must be exact solution, as long as it is convergent. In this case, for and , the exact solution and HAM solution are the same, as shown in Figure 5. The obtained numerical results are summarized in Table 2.

In Figure 6, we study the diagrams of the results obtained by HAM for and in comparison with the exact solution (3.5); we can see the best value for in this case is .

5. Conclusion

In this paper, the homotopy analysis method (HAM) [2] is applied to obtain the solitary solution of the KdV and Burger equations. HAM provides us with a convenient way to control the convergence of approximation series, which is a fundamental qualitative difference in analysis between HAM and other methods. So, these examples show the flexibility and potential of the homotopy analysis method for complicated nonlinear problems in engineering.

Acknowledgments

This research is partially funded by MOHE FRGS Vot no. 78675 and UTM RUG Vot no. PY/2011/02418. The authors wish to thank the anonymous reviewers whose comments led to some improvement to the presentation of the results.