`Mathematical Problems in EngineeringVolume 2008 (2008), Article ID 795838, 5 pageshttp://dx.doi.org/10.1155/2008/795838`
Research Article

## Homotopy Perturbation Method for Solving Reaction-Diffusion Equations

School of Information Engineering, Zhejiang Forestry College, Lin'an 311300, Zhejiang, China

Received 16 November 2007; Revised 11 February 2008; Accepted 27 February 2008

Copyright © 2008 Yu-Xi Wang 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.

#### Abstract

The homotopy perturbation method is applied to solve reaction-diffusion equations. In this method, the trial function (initial solution) is chosen with some unknown parameters, which are identified using the method of weighted residuals. Some examples are given. The obtained results are compared with the exact solutions, revealing that this method is very efficient and the obtained solutions are of high accuracy.

#### 1. Introduction

In this paper, we consider a reaction-diffusion process governed by the nonlinear ordinary differential equation [1]: with boundary conditions where represents the steady-state temperature for the corresponding reaction-diffusion equation with the reaction term is the power of the reaction term (heat source), generally it follows is the length of the sample (heat conductor). The physical interpretation of (1.1) was given in [1].

Recently, various different analytical methods were applied to nonlinear equations arising in engineering applications, such as the homotopy perturbation method [210], and exp-function method [11, 12], a complete review is available in [13]. This problem was studied by Lesnic using Adomian method [1], and by Mo [14] using variational method. In this paper, the homotopy perturbation method [2, 3, 13] is applied to the discussed problem, and the obtained results show that the method is very effective and simple.

#### 2. Homotopy Perturbation Method

In order to use the homotopy perturbation, we construct a homotopy in the form [2, 3, 13] with initial approximation where is an unkown constant to be further determined. It is obvious that (2.2) satisfies the boundary conditions.

We rewrite (2.1) in the form of We suppose the solution of (2.3) has the form Substituting (2.4) into (1.1) and equating the terms with the identical powers of , we can solve sequentially with ease. Setting , we obtain the approximate solution of (1.1) in the form of To illustrate its solution procedure, we consider some special cases.

Case 1 (). Under such case, we can easily obtain sequentially We, therefore, obtain the approximate solution in the form of In order to identify the unknown constant in (2.7), we apply the method of weighted residuals. Subsituting (2.7) into (1.1) results in the following residual: It is obvious that and . We locate at , and set , yielding the result

Case 2 (). The solution procedure is the same as that for Case 1. We can easily obtain the following linear equations: We obtain the following second-order approximate solution: Similarly, we locate at , and set to identify the unknown constant, which reads .

Case 3 (). By the same manuplation as illustrated in above cases, we obtain Using the method of weighted residuals, we set , resulting in .

Figure 1 shows the remarkable accuracy of the obtained results.

Figure 1: Comparison of approximate solutions with exact ones. Continued line: approximate solution; discontinued line: exact solution.

#### 3. Conclusion

The homotopy perturbation method deforms a complex problem under study to a simple problem routinely. If initial guess is suitably chosen, one iteration is enough, making the method a most attractive one. The method is of remarkable simplicity, while the obtained results are of utter accuracy on the whole solution domain. The method can be applied to various other nonlinear problems without any difficulty.

#### References

1. D. Lesnic, “A nonlinear reaction-diffusion process using the Adomian decomposition method,” International Communications in Heat and Mass Transfer, vol. 34, no. 2, pp. 129–135, 2007.
2. J.-H. He, “New interpretation of homotopy perturbation method,” International Journal of Modern Physics B, vol. 20, no. 18, pp. 2561–2568, 2006.
3. J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.
4. T. Ozis and A. Yildirim, “A comparative study of He's homotopy perturbation method for determining frequency-amplitude relation of a nonlinear oscillator with discontinuities,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 243–248, 2007.
5. A. Belendez, A. Hernandez, T. Belendez et al., “Application of He's homotopy perturbation method to the Duffing-harmonic oscillator,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 1, pp. 79–88, 2007.
6. X.-C. Cai, W.-Y. Wu, and M.-S. Li, “Approximate period solution for a kind of nonlinear oscillator by He's perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 109–112, 2006.
7. M. Dehghan and F. Shakeri, “Solution of an integro-differential equation arising in oscillating magnetic fields using He's homotopy perturbation method,” Progress in Electromagnetics Research, vol. 78, pp. 361–376, 2008.
8. F. Shakeri and M. Dehghan, “Inverse problem of diffusion equation by He's homotopy perturbation method,” Physica Scripta, vol. 75, no. 4, pp. 551–556, 2007.
9. M. Dehghan and F. Shakeri, “Solution of a partial differential equation subject to temperature overspecification by He's homotopy perturbation method,” Physica Scripta, vol. 75, no. 6, pp. 778–787, 2007.
10. D. D. Ganji and A. Sadighi, “Application of He's homotopy-perturbation method to nonlinear coupled systems of reaction-diffusion equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 4, pp. 411–418, 2006.
11. S.-D. Zhu, “Exp-function method for the Hybrid-Lattice system,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 3, pp. 461–464, 2007.
12. S.-D. Zhu, “Exp-function method for the discrete mKdV lattice,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 3, pp. 465–468, 2007.
13. J.-H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.
14. L.-F. Mo, “Variational approach to reaction-diffusion process,” Physics Letters A, vol. 368, no. 3-4, pp. 263–265, 2007.