Analytical Study of Some Important Generalized Nonlinear Partial Differential Equations

Marwan Alquran, Mahmoud Mohammad, and Ahmad Ababneh Department of Mathematics and Statistics, Jordan University of Science and Technology, Irbid 22110, Jordan Correspondence should be addressed to Marwan Alquran; marwan04@just.edu.jo Received 4 January 2013; Accepted 1 March 2013 Academic Editor: Anjan Biswas Copyright © 2013 Marwan Alquran et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 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.


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 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 (, ) : Ω × [0, 1] → , which satisfies which is equivalent to where  ∈ Ω,  ∈ [0, 1], that is, called embedding parameter,  0 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 () − ( 0 ) to () − ().In topology, this is called deformation, and () − ( 0 ) 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  = 1 gives the solution of ( 7) To study the convergence of the method, rewrite (11) in the following form: applying the inverse operator  −1 , 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,  → 1.
(2) The norm of  −1 (/) must be smaller than one so that the series converges.

Fisher Equation with Convection Term
Consider the Fisher equation with convection term as follows: By means of homotopy perturbation technique, the homotopy construction of ( 19) is or where the initial approximation  0 = (, 0).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  0 ,  1 ,  2 , and so forth: )) , )) , ))) . . .
Considering the first 11 terms of (23), then the approximate solution of (19) by setting  = 1 is )) The exact solution of (19) is The behaviors of (, ) and  app (, ) 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.

Sharma-Tasso-Olver Equation
Consider the Sharma-Tasso-Olver equation as follows: By means of homotopy perturbation technique, the homotopy construction of ( 28) is or where the initial approximation  0 = (, 0).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  0 ,  1 ,  2 , and so forth: Considering the first 8 terms of (32), then the approximate solution of (28) by setting  = 1 is

− 𝛼
The exact solution of (28) is 4.1.=1.The behavior of (, ) and  app (, ) is shown in Figure 4 for the regions −10 ≤  ≤ 10 and 0 ≤  ≤ 0.3.Absolute errors for differences between the exact solution and the 8th order approximate solution given by HPM for  = 1 are given in Table 2.

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

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  0 = (, 0).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  0 ,  1 ,  2 , and so forth: Journal of Mathematics Considering the first 11 terms of (41), then the approximate solution of (37) by setting  = 1 is  The exact solution of (37) is Figure 6 shows that the best choice of  is −5 ≤  ≤ 5, which gives solution that converges rapidly than other choices for .When  = −1, 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.
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.

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.