Mathematical Problems in Engineering

Volume 2009 (2009), Article ID 427516, 12 pages

http://dx.doi.org/10.1155/2009/427516

## Solving Heat and Wave-Like Equations Using He's Polynomials

Basic Sciences, HITEC University, 44000 Taxila Cantt, Pakistan

Received 21 May 2009; Revised 30 June 2009; Accepted 6 July 2009

Academic Editor: Ji Huan He

Copyright © 2009 Syed Tauseef Mohyud-Din. 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

We use He's polynomials which are calculated form homotopy perturbation method (HPM) for solving heat and wave-like equations. The proposed iterative scheme finds the solution without any discretization, linearization, or restrictive assumptions. Several examples are given to verify the reliability and efficiency of the method. The fact that suggested technique solves nonlinear problems without using Adomian's polynomials is a clear advantage of this algorithm over the decomposition method.

#### 1. Introduction

The heat and wave-like models are the integral part of applied sciences and arise in various physical phenomena. Several techniques including spectral, characteristic, modified variational iteration and Adomian’s decomposition have been used for solving these problems; see [1–3] and the references therein. Most of these techniques encounter a considerable size of difficulty. He [4–11] developed and formulated homotopy perturbation method (HPM) by merging the standard homotopy and perturbation. The homotopy perturbation method (HPM) proved to be compatible with the versatile nature of the physical problems and has been applied to a wide class of functional equations; see [1, 4–19] and the references therein. In this technique, the solution is given in an infinite series usually converging to an accurate solution; see [1, 4–19] and the references therein. It is worth mentioning that HPM is applied without any discretization, restrictive assumption or transformation and is free from round off errors. The HPM is applied for all the nonlinear terms in the problem without discretizing either by finite difference or by spline techniques at the nodes and involves laborious calculations coupled with a strong possibility of the ill-conditioned resultant equations which are a complicated problem to solve. Moreover, unlike the method of separation of variables that requires initial and boundary conditions, the homotopy perturbation method (HPM) provides an analytical solution by using the initial conditions only. The fact that HPM solves nonlinear problems without using Adomian’s polynomials is a clear advantage of this technique over the decomposition method. It is worth mentioning that [12, 13] introduced He’s polynomials by splitting the nonlinear term and proved their compatibility with Adomian’s polynomials coupled with the conclusion that He’s polynomials are easier to calculate, are more user friendly, and are independent of the complexities arising in calculating the so-called Adomian’s polynomials. It is to be highlighted that He’s polynomials are calculated from homotopy perturbation method (HPM). Inspired and motivated by the ongoing research in this area, we use He’s polynomials for solving heat and wave-like equations. It is worth mentioning that Noor and Mohyud-Din [20] introduced a homotopy approach which involves an additional term and consequently leads towards laborious and redundant calculations, whereas the approach used in Section 2 is more precise and easier to implement. Moreover, it reduces the huge unnecessary calculation arising in [20]. Several examples are given to verify the reliability and efficiency of the algorithm.

#### 2. Homotopy Perturbation Method and He’s Polynomials

To explain the homotopy perturbation method, we consider a general equation of the type where is any integral or differential operator. We define a convex homotopy by where is a functional operator with known solutions , which can be obtained easily. It is clear that, for we have This shows that continuously traces an implicitly defined curve from a starting point to a solution function . The embedding parameter monotonically increases from zero to unit as the trivial problem is continuously deforms the original problem . The embedding parameter can be considered as an expanding parameter [1, 2, 4–19]. The homotopy perturbation method uses the homotopy parameter as an expanding parameter [4–11] to obtain if , then (2.5) corresponds to (2.2) and becomes the approximate solution of the form It is well known that series (2.6) is convergent for most of the cases and also the rate of convergence is dependent on ; see [4–11]. We assume that (2.6) has a unique solution. The comparisons of like powers of give solutions of various orders. In sum, according to [12, 13], He’s HPM considers the solution, , of the homotopy equation in a series of as follows: and the method considers the nonlinear term as where ’s are the so-called He’s polynomials [12, 13], which can be calculated by using the formula

#### 3. Numerical Applications

In this section, we use He’s polynomials which are calculated from homotopy perturbation method (HPM) for solving heat and wave-like equations.

*Example 3.1 ([1, 2]). *Consider the one-dimensional initial boundary value problem which describes the heat-like models
with boundary conditions
and initial conditions
Apply the convex homotopy
Compare the coefficient of like powers of
where are He’s polynomials. The series solution is given by
and in a closed form by

*Example 3.2 ([1, 2]). *Consider the two-dimensional initial boundary value problem which describes the heat-like models
with boundary conditions
and initial conditions
Apply the convex homotopy method
Compare the coefficient of like powers of
where are He’s polynomials. The series solution is given by
and in a closed form by

*Example 3.3 ([1, 2]). *Consider the three-dimensional inhomogeneous initial boundary value problem which describes the heat-like models
subject to the boundary conditions
and the initial conditions
Apply the convex homotopy method
Compare the coefficient of like powers of
where are He’s polynomials. The series solution is given by
and in a closed form by

*Example 3.4. *Consider the following nonlinear heat-like model:
subject to the initial conditions
Apply the convex homotopy method
Compare the coefficient of like powers of
where are He’s polynomials. The series solution is given by
and the closed form solution is given as

*Example 3.5 ([1, 2]). *Consider the one-dimensional initial boundary value problem which describes the wave-like models
subject to the boundary conditions
and initial conditions
Apply the convex homotopy method
Compare the coefficient of like powers of
where are He’s polynomials. The series solution is given by
and in a closed form by

*Example 3.6 ([1, 2]). *Consider the two-dimensional initial boundary value problem which describes the wave-like models
subject to the Neumann boundary conditions
and initial conditions
Apply the convex homotopy method
Compare the coefficient of like powers of
where are He’s polynomials. The series solution is given by
and in a closed form by

*Example 3.7 ([1, 2]). *Consider the three-dimensional inhomogeneous initial boundary value problem which describes the wave-like models
subject to the boundary conditions
and the initial conditions
Apply the convex homotopy method
Compare the coefficient of like powers of and proceede as before, the series solution is given by
and in a closed form by

#### 4. Conclusion

In this paper, we use He’s polynomials which are calculated from homotopy perturbation method (HPM) for solving heat and wave-like equations. The method is applied in a direct way without using linearization, transformation, discretization, or restrictive assumptions. It may be concluded that the proposed scheme is very powerful and efficient in finding the analytical solutions for a wide class of boundary value problems. The method gives more realistic series solutions that converge very rapidly in physical problems. It is also observed that He’s polynomials are compatible with Adomian’s Polynomials but are easier to calculate and are more user friendly.

#### Acknowledgments

The author is highly grateful to the referees and Professor Dr. Ji-Huan He for their very constructive comments. The author would like to thank Brig (R) Qamar Zaman, Vice Chancellor, HITEC University, Taxila Cantt, Pakistan for the provision of conducive environs for research.

#### References

- M. A. Noor and S. T. Mohyud-Din, “Modified variational iteration method for heat and wave-like equations,”
*Acta Applicandae Mathematicae*, vol. 104, no. 3, pp. 257–269, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A.-M. Wazwaz and A. Gorguis, “Exact solutions for heat-like and wave-like equations with variable coefficients,”
*Applied Mathematics and Computation*, vol. 149, no. 1, pp. 15–29, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Wilcox, “Closed-form solution of the differential equation $({\partial}^{2}/\partial x\partial y+ax(\partial /\partial x)+by(\partial /\partial y)+cxy+\partial /\partial t)P(x,y,t)=0$ by normal-ordering exponential operators,”
*Journal of Mathematical Physics*, vol. 11, no. 4, pp. 1235–1237, 1970. View at Publisher · View at Google Scholar · View at MathSciNet - J.-H. He, “An elementary introduction to recently developed asymptotic methods and nanomechanics in textile engineering,”
*International Journal of Modern Physics B*, vol. 22, no. 21, pp. 3487–3578, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-H. He, “Recent development of the homotopy perturbation method,”
*Topological Methods in Nonlinear Analysis*, vol. 31, no. 2, pp. 205–209, 2008. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “Some asymptotic methods for strongly nonlinear equations,”
*International Journal of Modern Physics B*, vol. 20, no. 10, pp. 1141–1199, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “Homotopy perturbation method for solving boundary value problems,”
*Physics Letters A*, vol. 350, no. 1-2, pp. 87–88, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - J.-H. He, “Comparison of homotopy perturbation method and homotopy analysis method,”
*Applied Mathematics and Computation*, vol. 156, no. 2, pp. 527–539, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 6, no. 2, pp. 207–208, 2005. View at Google Scholar - J.-H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,”
*Applied Mathematics and Computation*, vol. 151, no. 1, pp. 287–292, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - A. Ghorbani and J. Saberi-Nadjafi, “He's homotopy perturbation method for calculating adomian polynomials,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 8, no. 2, pp. 229–232, 2007. View at Google Scholar - A. Ghorbani, “Beyond Adomian polynomials: He polynomials,”
*Chaos, Solitons & Fractals*, vol. 39, no. 3, pp. 1486–1492, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - S. T. Mohyud-Din, M. A. Noor, and K. I. Noor, “Traveling wave solutions of seventh-order generalized KdV equations using He's polynomials,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 10, no. 2, pp. 227–233, 2009. View at Google Scholar - S. T. Mohyud-Din and M. A. Noor, “Homotopy perturbation method for solving fourth-order boundary value problems,”
*Mathematical Problems in Engineering*, vol. 2007, Article ID 98602, 15 pages, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. T. Mohyud-Din, M. A. Noor, and K. I. Noor, “Homotopy perturbation method for unsteady flow of gas through a porous medium,”
*International Journal of Modern Physics B*. In press. - S. T. Mohyud-Din, M. A. Noor, and K. I. Noor, “On the coupling of polynomials with correction functional,”
*International Journal of Modern Physics B*. In press. - M. A. Noor and S. T. Mohyud-Din, “Homotopy perturbation method for solving nonlinear higher-order boundary value problems,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 9, no. 4, pp. 395–408, 2008. View at Google Scholar - L. Xu, “He's homotopy perturbation method for a boundary layer equation in unbounded domain,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1067–1070, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - M. A. Noor and S. T. Mohyud-Din, “Modified variational iteration method for heat and wave-like equations,”
*Acta Applicandae Mathematicae*, vol. 104, no. 3, pp. 257–269, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet