Abstract

The aim of this paper is to extend the use of homotopy perturbation method (HPM) to study the solutions for some important generalized nonlinear partial differential equations (PDEs) such as Fisher equation with convection term, Sharma-Tasso-Olver (STO) equation, and Fitzhugh-Nagumo (FN) equation.

1. Introduction

The main purpose of this paper is to apply the homotopy perturbation method with the help of symbolic computation to obtain approximate solution of the following nonlinear equations.

Fisher equation with convection term serves as a basic model for population dynamics or chemical kinetics in mathematical biology. The solution shows exponential growth counteracted by nonlinear damping (where the damping rate is proportional to ). Fisher equation has the form subject to

Fisher equation was studied by many authors using different methods such as the tanh method [1] and variational iteration method [2].

Sharma-Tasso-Olver (STO) equation is a good example to show the fusion of soliton solutions. It has the form subject to

Hirota’s direct method [3] and extended tanh method [4] were used to solve the STO equation.

Fitzhugh-Nagumo (FN) equation is an important nonlinear reaction-diffusion equation and is usually used to model the transmission of nerve impulses; also it is used in circuit theory, biology, and the area of population genetics. FN equation has the form subject to

FN equation was studied by means of improved sine-cosine method [5] and Exp-function method [6].

2. Survey of Homotopy Perturbation Method

To illustrate the basic ideas of this method, we consider the following nonlinear differential equation: with the boundary conditions where is a general differential operator, is a boundary operator, is a known analytical function, and is the boundary of the domain . The operator can be divided into two parts, which are and , where is the linear, and is the nonlinear operator. Therefore, (7) can be written as follows:

Liao in [7] constructed a homotopy , which satisfies which is equivalent to where , , that is, called embedding parameter, is an initial approximation solution of (7), which satisfies the boundary conditions. Obviously, from (10) and (11) we obtain and the changing process of from 0 to 1, is just that of from to . In topology, this is called deformation, and and are called homotopic. Applying the perturbation technique, we can first use the embedding parameter as a small parameter and assume that the solutions of (10) or (11) can be expressed as a power series in as follows:

Setting gives the solution of (7)

To study the convergence of the method, rewrite (11) in the following form: applying the inverse operator , to both sides of (15), we obtain substituting (13) into the right-hand side of (16), we get the following form:

The exact solution may be obtained by using (14)

The series (14) is convergent for most cases. However the following suggestions have been made by He [8], to find the convergence rate on nonlinear operator.(1)The second derivative of with respect to must be small because the parameter may be relatively large, that is, .(2)The norm of must be smaller than one so that the series converges.

More about this method and illustrative examples, the reader may refer to the following articles [917].

3. Fisher Equation with Convection Term

Consider the Fisher equation with convection term as follows: subject to

By means of homotopy perturbation technique, the homotopy construction of (19) is or where the initial approximation . Suppose that the solution of (19) has the form

Substituting (23) in (22) and comparing the coefficients of identical degrees of , we obtain the following linear equations:

Solving this system, we obtain the following solutions for , , , and so forth:

Considering the first 11 terms of (23), then the approximate solution of (19) by setting is

The exact solution of (19) is

The behaviors of and are shown in Figure 1 and the absolute error is shown in Table 1.

Figures 2 and 3 show that the best choice of is 1.

4. Sharma-Tasso-Olver Equation

Consider the Sharma-Tasso-Olver equation as follows: subject to

By means of homotopy perturbation technique, the homotopy construction of (28) is or where the initial approximation . Suppose that the solution of (28) has the form

Substituting (32) in (31) and comparing the coefficients of identical degrees of , we obtain the following linear equations:

Solving this system, we obtain the following solutions for , , , and so forth:

Considering the first 8 terms of (32), then the approximate solution of (28) by setting is

The exact solution of (28) is

4.1.

The behavior of and is shown in Figure 4 for the regions and . Absolute errors for differences between the exact solution and the 8th order approximate solution given by HPM for are given in Table 2.

4.2.

In this part, we consider and compare it with the above case . Figure 5 shows this comparison. One can see that for the case the solution converges rapidly than .

5. Fitzhugh-Nagumo Equation

Consider the Fitzhugh-Nagumo (FN) equation as follows: subject to

By means of homotopy perturbation technique, the homotopy construction of (37) is or where the initial approximation . Suppose that the solution of (37) has the form

Substituting (41) in (40) and comparing the coefficients of identical degrees of , we obtain the following linear equations:

Solving this system, we obtain the following solutions for , , , and so forth:

Considering the first 11 terms of (41), then the approximate solution of (37) by setting is

The exact solution of (37) is

Figure 6 shows that the best choice of is , which gives solution that converges rapidly than other choices for . When , the FN equation reduces to the Newell-Whitehead (NW) equation, which is an important nonlinear reaction-diffusion equation and is usually used to model the transmission of nerve impulses, also is used in circuit theory, biology, and the area of population genetics as mathematical models.

The behavior of and is shown in Figure 7 for the regions and , when .

Absolute errors for differences between the exact solution and the 11th order approximate solution given by HPM for NW equation are given in Table 3.

6. Conclusion

In this paper, the homotopy perturbation method has been used for finding the approximate solutions of Fisher equation, Sharma-Tasso-Olver equation, and Fitzhugh-Nagumo equation. HPM provides highly accurate numerical solutions for our problems and they also do not require large computer memory and discretization of variable . The approximations are not valid only for small parameters but also for larger ones and the initial approximation can be arbitrarily chosen with unknown constants.

The small size of computations in comparison with the computational size required in characteristics method and the rapid convergence show that the homotopy perturbation method is more reliable and introduces a significant improvement in solving nonlinear partial differential equations.