On the Exact Solution of Burgers-Huxley Equation Using the Homotopy Perturbation Method

The Homotopy Perturbation Method (HPM) is used to solve the Burgers-Huxley non-linear differential equations. Three case study problems of Burgers-Huxley are solved using the HPM and the exact solutions are obtained. The rapid convergence towards the exact solutions of HPM is numerically shown. Results show that the HPM is efficient method with acceptable accuracy to solve the Burgers-Huxley equation. Also, the results prove that the method is an efficient and powerful algorithm to construct the exact solution of non-linear differential equations.

• Burgers-Huxley Equation
• Homotopy Perturbation Method
• Nonlinear Differential Equations

1. Introduction

Most of the nonlinear differential equations do not have an analytical solution. Recently, semi-analytical solutions of real-life mathematical modeling are considered as a key tool to solve nonlinear differential equations.

The idea of the Homotopy Perturbation Method (HPM) which is a semi-analytical method was first pioneered by He [1] . Later, the method is applied by He [2] to solve the non-linear non-homogeneous partial differential equations. Nourazar et al. [3] used the homotopy perturbation method to find exact solution of Newell-Whitehead- Segel equation. Krisnangkura et al. [4] obtained exact traveling wave solutions of the generalized Burgers- Huxley equation by using hyperbolic tangent method. New exact solutions of the generalized Burgers-Huxley equation are also obtained by Gao and Zhao [5] . Hashim et al. [6] have solved the generalized Burgers-Huxley equation by using adomian decomposition method. Wang et al. [7] studied the solitary wave solutions of the generalized Burgers-Huxley equation. Also, Darvishi et al. [8] have used spectral collocation method and Darvishi’s preconditionings to solve the generalized Burgers-Huxley equation.

The generalized Burgers-Huxley equation is investigated by Satsuma [9] as:

where are real constants and is a positive integer and. Equation (1.1) models the interaction between reaction mechanisms, convection effects and diffusion transports [9] .

When and, Equation (1.1) is reduced to the Huxley equation which describes nerve pulse propagationin nerve fibers and wall motion in liquid crystals [10] . When and, Equation (1.1) is reduced to the Burgers equation describing the far field of wave propagation in nonlinear dissipative systems [11] . When, and, Equation (1.1) is turned into the Burgers-Huxley equation. The Burgers-Hux- ley equation is used to model the interaction between reaction mechanisms, convection effects and diffusion transport, nerve pulse propagation in nerve fibers as well as wall motion in liquid crystals.

In the present research work, the Homotopy Perturbation Method (HPM) is applied to obtain the closed form solution of the non-linear Burgers-Huxley equation. Three case study problems of non-linear Burgers-Huxley equations are solved by using the HPM. The trend of the rapid convergence towards the exact solution is shown when compared to the exact solution.

The idea of homotopy perturbation method is presented in Section 2. Application of the homotopy perturbation method to the exact solution of Burgers-Huxley equation is presented in Section 3.

2. The Idea of Homotopy Perturbation Method

The Homotopy Perturbation Method (HPM) is originally initiated by He [1] . This is a combination of the classical perturbation technique and homotopy technique. The basic idea of the HPM for solving nonlinear differential equations is as follow; consider the following differential equation:

where is any differential operator. We construct a homotopy as follow:

where is a functional operator with the known solution. It is clear that when is equal to zero then, and when is equal to 1, then. It is worth noting that as the embedding parameter increases monotonically from zero to unity the zero order solution continuously deforms into the original problem. The embedding parameter, , is considered as an expanding parameter [2] . In the homotopy perturbation method the embedding parameter is used to get series expansion for solution as:

when, then Equation (2.2) becomes the approximate solution to Equation (2.1) as:

The series Equation (2.4) is a convergent series and the rate of convergence depends on the nature of Equation (2.1) [1] [2] . It is also assumed that Equation (2.2) has a unique solution and by comparing the like powers of the solution of various orders is obtained. These solutions are obtained using the Maple package.

3. The Burgers-Huxley Equation

To illustrate the capability and reliability of the method, three cases of nonlinear diffusion equations are presented.

Case І: in this case we will examine the Burgers-Huxley equation for, so, the equation is written as:

Subject to initial condition:

We construct a homotopy for Equation (3.1) in the following form:

The solution of Equation (3.1) can be written as a power series in p as:

Substituting Equation (3.4) and Equation (3.2) into Equation (3.3) and equating the terms with identical powers of p:

Using the Maple package to solve recursive sequences, Equation (3.5), we obtain the followings:

By setting p = 1 in Equation (3.4) the solution of Equation (3.1) can be obtained as

Therefore the solution of Equation (3.1) is written as:

The Taylor series expansion for is written as:

Combining Equation (3.8) with Equation (3.7), we get as follow:

This is the exact solution of the problem, Equation (3.1). Table 1 shows the trend of rapid convergence of the results of to using the HPM. The rapid convergence of the solution

toward the exact solution, the maximum relative error of less than 0.0000058% is achieved as shown in Table 1.

Case ІІ: In Equation (1.1) for, the Burgers-Huxley equation is written as:

Table 1

The percentage of relative errors of the results of to of the HPM solution of Equation (3.1)

		Percentage of relative error (%RE)
		0.01693168743	0.01002710463	0.005488150424
		0.000002337346256	2.025644856e−7	9.527700575e−7
		2.153484215e−10	3.464089104e−10	3.155246333e−10
		9.744801271e−12	3.920421861e−11	9.937961084e−11
		0.05344388963	0.03164997413	0.01732302882
		0.00006806597676	0.000003967422423	0.00002580410708
		6.184344787e−8	9.467317545e−8	5.516785201e−8
		1.008793018e−8	1.166422821e−9	2.026398250e−9
		0.07311570399	0.04329980772	0.02369935005
		0.0001675410515	0.000007498069748	0.00006123279256
		2.799248652e−7	4.009848425e−7	2.364525175e−7
		5.775483086e−8	6.758498122e−9	1.111128458e−8

Subject to initial condition:

To solve Equation (3.10) we construct a homotopy in the following form:

The solution of Equation (3.10) can be written as a power series in p as:

Substituting Equation (3.13) and Equation (3.11) in to Equation (3.12) and equating the term with identical powers of p, leads to:

Using the Maple package to solve recursive sequences, Equation (3.14), we obtain the followings:

By setting p = 1 in Equation (3.13), the solution of Equation (3.10) can be obtained as

Therefore the solution of Equation (3.10) is written as:

The Taylor series expansion for is written as:

By substituting Equation (3.17) into Equation (3.16), Equation (3.16) can be reduced to:

This is the exact solution of the problem, Equation (3.10). Table 2 shows the trend of rapid convergence of the results of to using the HPM solution toward the exact solution.

The maximum relative error of less than 0.00014% is achieved in comparison to the exact solution as shown in Table 2.

Case ІІІ: In Equation (1.1) for, the Burgers-Huxley equation becomes:

Subject to initial condition:

Table 2

The percentage of relative errors of the results of to of the HPM solution of Equation (3.10)

		Percentage of relative error (%RE)
		0.0484797171	0.056937877	0.063676094
		0.0000184239461	0.000009125432	0.000006989146
		7.8094040e−9	3.9049212e−9	1.37134980e−8
		3.61301281e−10	6.1915442e−11	1.3095951e−10
		0.1570606291	0.184462686	0.206292613
		0.000524561340	0.00023856284	0.00024584133
		0.00000177771365	0.00000129074423	0.00000380612758
		2.19608021e−7	2.10498615e−7	9.144252e−9
		0.2177728801	0.255767283	0.2860356311
		0.001277016710	0.00055367038	0.00065831599
		0.0000075558572	0.0000060480463	0.0000170599007
		0.000001302473287	0.000001221861195	8.044206e−8

We construct a homotopy for Equation (3.19) in the following form:

The solution of Equation (3.19) can be written as a power series in p as:

Substituting Equation (3.22) and Equation (3.20) into Equation (3.21) and equating the terms with identical powers of p:

Using the Maple package to solve recursive sequences, Equation (3.23), we obtain the followings:

By setting p = 1 in Equation (3.22) the solution of Equation (3.19) can be obtained as

Thus the solution of Equation (3.19) can be written as:

The Taylor series expansion for is written as:

Comparing Equation (3.26) with Equation (3.25), thus Equation (3.25) can be reduced to:

This is the exact solution of the problem, Equation (3.19). Table 3 shows the trend of rapid convergence of the

results of to using the HPM solution toward the exact solution. The

maximum relative error of less than 0.038% is achieved in comparison to the exact solution as shown in Table 3.

Deng [12] obtained some travelling solitary wave solutions of Equation (1.1) by applying the first-integral method as follows:

Table 3

The percentage of relative errors of the results of to of the HPM solution of Equation (3.19)

		Percentage of relative error (%RE)
		0.1473751972	0.1768549738	0.1894968996
		0.00008115001396	0.0004703355304	0.0008489173163
		5.853207295e−7	0.000001136972415	2.175815927e−7
		4.468762836e−8	5.878826669e−8	2.556253934e−8
		0.5347070619	0.6416656548	0.6875331484
		0.001445157793	0.01763206139	0.03053744832
		0.0002129909216	0.0003484139717	0.00009299132649
		0.00003726778757	0.00005691935052	0.00002679413014
		0.7871664200	0.9446250035	1.012148624
		0.002172311991	0.04914756514	0.08362907281
		0.001086228788	0.001651183957	0.0005006633300
		0.0002239213593	0.0003721014272	0.0001680489489

where

and is arbitrary constant.

This is in full agreement of the closed form solutions of Equation (3.1), Equation (3.10) and Equation (3.19) for differences value of parameters and in the three cases. So, it can be concluded that the HPM is a powerful and efficient technique to solve the non-linear Burgers-Huxley equation.

4. Conclusion

In the present research work, the exact solution of the Burgers-Huxley nonlinear diffusion equation is obtained using the HPM. The validity and effectiveness of the HPM is shown by solving three non-homogenous non-linear differential equations and the very rapid convergence to the exact solutions is also numerically demonstrated. The trend of rapid and monotonic convergence of the solution toward the exact solution is clearly shown by obtaining the relative error in comparison to the exact solution. The rapid convergence towards the exact solutions of HPM indicates that, using the HPM to solve the non-linear differential equations, a reasonable less amount of computational work with acceptable accuracy may be sufficient. Moreover, it can be concluded that the HPM is a very powerful and efficient technique which can construct the exact solution of nonlinear differential equations.

References

1. Deng, X.J. (2008) Travelling Wave Solutions for the Generalized Burgers-Huxley Equation. Applied Mathematics and Computation, 204, 733-737. http://dx.doi.org/10.1016/j.amc.2008.07.020

2. Whitham, G.B. (1974) Linear and Nonlinear Waves. Wiley, New York.

3. Wang, X.Y. (1985) Nerve Propagation and Wall in Liquid Crystals. Physics Letters A, 112, 402-406.http://dx.doi.org/10.1016/0375-9601(85)90411-6

4. Satsuma, J. (1987) Topics in Soliton Theory and Exactly Solvable Nonlinear Equations. World Scientific, The Singapore City.

5. Darvishi, M.T., Kheybari, S. and Khani, F. (2008) Spectral Collocation Method and Darvishi’s Preconditionings to Solve the Generalized Burgers-Huxley Equation. Communications in Nonlinear Science and Numerical Simulation, 13, 2091-2103. http://dx.doi.org/10.1016/j.cnsns.2007.05.023

6. Wang, X.Y., Zhu, Z.S. and Lu, Y.K. (1990) Solitary Wave Solutions of the Generalised Burgers-Huxley Equation. Journal of Physics A: Mathematical and General, 23, 271. http://dx.doi.org/10.1088/0305-4470/23/3/011

7. Hashim, I., Noorani, M.S.M. and Said Al-Hadidi, M.R. (2006) Solving the Generalized Burgers-Huxley Equation Using the Adomian Decomposition Method. Mathematical and Computer Modelling, 43, 1404-1411.http://dx.doi.org/10.1016/j.mcm.2005.08.017

8. Gao, H. and Zhao, R.X. (2010) New Exact Solutions to the Generalized Burgers-Huxley Equation. Applied Mathematics and Computation, 217, 1598-1603. http://dx.doi.org/10.1016/j.amc.2009.07.020

9. Krisnangkura, M., Chinviriyasit, S. and Chinviriyasit, W. (2012) Analytic Study of the Generalized Burger’s-Huxley Equation by Hyperbolic Tangent Method. Applied Mathematics and Computation, 218, 10843-10847.http://dx.doi.org/10.1016/j.amc.2012.04.044

10. Nourazar, S.S., Soori, M. and Nazari-Golshan, A. (2011) On the Exact Solution of Newell-Whitehead-Segel Equation Using the Homotopy Perturbation Method. Australian Journal of Basic and Applied Sciences, 5, 1400-1411.

11. He, J.H. (2005) Application of Homotopy Perturbation Method to Nonlinear Wave Equations. Chaos Solitons & Fractals, 26, 695-700. http://dx.doi.org/10.1016/j.chaos.2005.03.006

12. He, J.H. (1999) Homotopy Perturbation Technique. Computer Methods in Applied Mechanics and Engineering, 178, 257-262. http://dx.doi.org/10.1016/S0045-7825(99)00018-3