Research Article | Open Access
Afgan Aslanov, "Homotopy Perturbation Method for Solving Wave-Like Nonlinear Equations with Initial-Boundary Conditions", Discrete Dynamics in Nature and Society, vol. 2011, Article ID 534165, 10 pages, 2011. https://doi.org/10.1155/2011/534165
Homotopy Perturbation Method for Solving Wave-Like Nonlinear Equations with Initial-Boundary Conditions
The homotopy perturbation method is employed to obtain approximate analytical solutions of the wave-like nonlinear equations with initial-boundary conditions. An efficient way of choosing the auxiliary operator is presented. The results demonstrate reliability and efficiency of the method.
In this paper, we consider the equation with initial conditions and boundary condition where , , , , and are known functions.
Problems like (1.1)-(1.2)-(1.3) model many problems in classical and quantum mechanics, solitons, and matter physics [1, 2]. If is a function of only, we obtain a Klein-Gordon or sine-Gordon-type equations.
In the last decade, some various approximate methods have been developed, such as the homotopy perturbation method (HPM) [3–13] and Adomian’s decomposition method (ADM) [14–20] to solve linear and nonlinear differential equations.
Unlike the various approximation techniques for solving nonlinear wave type problems, which are usually valid for initial value problems (without boundary conditions) or some special type of problems (homogenous, etc.), our technique is applicable for all initial-boundary problems of type (1.1)-(1.2)-(1.3). Chowdhury and Hashim  applied the HPM for solving Klein-Gordon and sine-Gordon equations, with initial conditions (1.2). El-Sayed  and Wazwaz and Gorguis  used ADM for solving wave-like and heat-like problems. Their approaches cannot be applied for all wave-like equations with initial-boundary conditions since the operator cannot control the boundary condition (1.3)—see Example 2.1 below.
The central idea here is that the problem has a unique solution (see, e.g., ) and therefore there exists an inverse of the operator .
The main idea of HPM is to introduce a homotopy parameter, say , which takes values from 0 to 1. When the equation usually reduces to a sufficiently simplified form (linear or very easy nonlinear). As increases to 1, the equation goes through a sequence of “deformations” (homotopics) and at takes the original form of the equation.
We rewrite (1.1) as and construct the following homotopy: Usually we take as a solution of the problem (1.1)-(1.2)-(1.3) with or simply . Assume that the solution of (1.6) is in the form and substituting (1.7) into (1.6) and equating terms of the same powers of we obtain a system of equations for . Solving these system of equations we obtain a solution in the form
The HPM and ADM offer excellent choices for obtaining the closed-form analytical solutions of wave-like equations. Chowdhury and Hashim , El-Sayed , and Wazwaz and Gorguis  recently showed how the HPM and ADM can be applied to find an analytic approximate solution of wave-like equations with initial conditions (1.2). They mainly used the operator for solving the wave-like problems. But in case of inhomogeneous nonlinear or even linear equations with initial-boundary conditions, these approaches have some difficulties. If we construct the standard homotopy with , for solving wave-like equations with initial-boundary conditions, usually in the second or even in the first stage of HPM, we obtain an “overdetermined” or very difficult problem. To explain these difficulties we consider the following example.
Example 2.1. Consider the linear problem The exact solution is Using our HPM technique we can easily find this solution. Indeed, let us take , and construct the homotopy Now substituting (1.7) into and equating the coefficients of like powers of , we get a system of linear equations Solving correspondingly we obtain . Thus we obtain an exact solution .
Now let us show that this problem can not be solved when is taken as .
Indeed if we choose or (which seems most natural and appropriate) and and construct the homotopy
we obtain an overdetermined problem for in the form
which has no solution.
If in the first stage we consider the boundary conditions , , , we obtain and in the second stage we need to solve the problem
which is overdetermined again and has no solution.
Example 2.2. Consider the problem The exact solution is . We take again , , and construct the homotopy Substituting (1.7) into (2.10) we obtain and equating terms of the coefficients of like powers of gives
Solving these equations we obtain (see [21, Chapter 3.4])for and for . In a similar manner we have and therefore
for . In a similar manner we obtain for , where and so The absolute errors between the exact and three term approximation of the series solution for some values of are shown in Table 1.
Example 2.3. Now we consider the problem
Solving equations yields In a similar manner we get
for . Taking only two-term approximation for we have for . The absolute errors between the exact and two-term approximation of the series solution for some values of are shown in Table 2. A higher accuracy level can be attained by evaluating some more terms.
Example 2.4. Now we consider the problem It is easy to show that the HPM can not be applied if is taken as . The exact solution is We construct the homotopy with . Substituting (1.7) into and equating the coefficients of like powers of , we get a system of linear equations Solving we obtain Continuing we obtain Thus the three-term approximation is and the error is less than for .
Homotopy perturbation method has been successful for solving many linear and nonlinear wave-type problems. However, it has difficulties in dealing with initial boundary problems, namely, in including all initial and boundary conditions together into the process of homotopy perturbation and computation. Our main goal is to construct the homotopy perturbation scheme containing all initial and boundary conditions together. The goal is achieved by involving an auxiliary operator which includes both variables and .
- P. J. Caudray, I. C. Eilbeck, and J. D. Gibbon, “The sine-Gordon equation as a model classical field theory,” Nuovo Cimento, vol. 25, pp. 497–511, 1975.
- R. K. Dodd, I. C. Eilbeck, and J. D. Gibbon, Solitons and Nonlinear Wave Equations, Academic, London, UK, 1982.
- S. J. Liao, On the proposed homotopy analysis technique for nonlinear problems and its applications, Ph.D. dissertation, Shangai Jio Tong University, Shangai, China, 1992.
- J. H. He, “Variational iteration method-a kind of non-linear technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, pp. 699–708, 1999.
- 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.
- J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Caos, Solitons & Fractals, vol. 26, pp. 695–700, 2005.
- J. H. He, Non-Perturbative Methods for Strongly Nonlinear Problems, Die Deutsche Bibliothek, Leipzig, Germany, 2006.
- J. H. He, “New interpretation of homotopy perturbation method,” International Journal of Modern Physics B, vol. 20, no. 18, pp. 2561–2568, 2006.
- M. S. H. Chowdhury and I. Hashim, “Application of homotopy-perturbation method to Klein-Gordon and sine-Gordon equations,” Chaos, Solitons & Fractals, vol. 39, no. 4, pp. 1928–1935, 2009.
- M. J. Ablowitz, B. M. Herbst, and C. Schober, “Homotopy perturbation method and axisymmetric flow over a stretching sheet,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 4, pp. 399–406, 2006.
- M. E. Berberler and A. Yildirim, “He's homotopy perturbation method for solving the shock wave equation,” Applicable Analysis, vol. 88, no. 7, pp. 997–1004, 2009.
- L. Cveticanin, “The homotopy-perturbation method applied for solving complex-valued differential equations with strong cubic nonlinearity,” Journal of Sound and Vibration, vol. 285, no. 4-5, pp. 1171–1179, 2005.
- D. D. Ganji and A. Rajabi, “Assessment of homotopy-perturbation and perturbation methods in heat radiation equation,” International Communications in Heat and Mass Transfer, vol. 33, pp. 391–400, 2006.
- D. D. Ganji and A. Sadighi, “Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 24–34, 2007.
- G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.
- G. Adomian, Nonlinear Stochastic Systems and Applications to Physics, vol. 46 of Mathematics and its Applications, Kluwer Academic Publishers, Dordrecht, Germany, 1989.
- A. M. Wazwaz, “A comparison between Adomian decomposition method and Taylor series method in the series solutions,” Applied Mathematics and Computation, vol. 97, no. 1, pp. 37–44, 1998.
- D. Kaya, “A numerical solution of the sine-Gordon equation using the modified decomposition method,” Applied Mathematics and Computation, vol. 143, no. 2-3, pp. 309–317, 2003.
- S. M. El-Sayed, “The decomposition method for studying the Klein-Gordon equation,” Chaos, Solitons & Fractals, vol. 18, no. 5, pp. 1025–1030, 2003.
- 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.
- W. A. Strauss, Partial Differential Equations, John Wiley & Sons, New York, NY, USA, 1992.
Copyright © 2011 Afgan Aslanov. 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.