Research Article  Open Access
A HomotopyAnalysis Approach for Nonlinear WaveLike Equations with Variable Coefficients
Abstract
We are interested in the approximate analytical solutions of the wavelike nonlinear equations with variable coefficients. We use a wave operator, which provides a convenient way of controlling all initial and boundary conditions. The proposed choice of the auxiliary operator helps to find the approximate series solution without any discretization, linearization, or restrictive assumptions. Several examples are given to verify the reliability and efficiency of the method.
1. Introduction
We consider the equationwith initial conditionsand boundary conditionwhere , , , , , and are known functions.
Note that the proposed method can be applied for equations like with the same type of initialboundary conditions.
Problems like (1)–(3) model many problems in classical and quantum mechanics, solitons, and matter physics [1, 2]. If is a function of only and we obtain KleinGordon or sineGordon type equations. These models can describe some nonlinear phenomena; for example, wavelike equation can describe earthquake stresses [3], coupling currents in a flat multistrand twolayer superconducting cable [4], and nonhomogeneous elastic waves in soils [5]. Typical examples of the wavelike equations with variable coefficients are EulerTricomi equation [6] or Chaplygin equation [7] given bywhich are useful in the study of transonic flow, where is the stream function of a planeparallel steadystate gas flow, is positive at subsonic and negative at supersonic speed, and is the angle of inclination of the velocity vector. Some Chaplygin type of equation of the special formwhere is the speed of sound, has also applications in the study of transonic flow [8].
Note that we use the term wavelike equation to describe the partial differential equations with the terms and ; that is, the term “wavelike” may not correspond to the real physical waves, in general.
Recently, there has been a growing interest for obtaining the explicit solutions to wavelike and heatlike models by analytic techniques. Wazwaz [9] used the tanh method to obtain the exact solution of the sineGordon equation. Kaya [10] applied the modified decomposition method to solve the sineGordon equation. Aslanov [11] used homotopy perturbation method to solve KleinGordon type of equations with unbounded righthand side. ElSayed [12] and Wazwaz and Gorguis [13] used Adomian decomposition method for solving wavelike and heatlike problems.
The homotopy analysis method [14–18] is developed to search the accurate asymptotic solutions of nonlinear problems. Liao [17] proved that the homotopy analysis method (HAM) contains some other nonperturbation techniques, such as Adomian’s decomposition method and Lyapunov’s artificial small parameter. Hayat and Sajid [19] and Abbasbandy [20] pointed out that the homotopy perturbation method is only a special case of the HAM. Aslanov [21] used homotopy perturbation method to solve wavelike equations with initialboundary conditions. Rajaraman [22] and Alomari at al. [23] applied HAM for solving nonlinear equations with initial conditions. Öziş and Ağırseven [24] used homotopy perturbation method for solving heatlike and wavelike equations with variable coefficients.
Various methods for obtaining exact and approximate solutions to nonlinear partial differential equations have been proposed. Among these methods are the homotopy perturbation and Adomian decomposition methods [25–29], the variational iteration method [30], homotopy analysis method [31], and others.
Here we will further extend the applications of HAM to obtain an approximate series solution for the nonlinear wavelike equations with variable coefficients and with initialboundary conditions. The difficulty in the use of standard HAM is that the choice of the linear operator in standard form (like or ) cannot control the boundary conditions (2)(3).
Unlike the various approximation techniques, which are usually valid for problems with (only) initial conditions, our technique is applicable for a wide range of initialboundary problems of types (1)–(3). The central idea here is that the problem,has a unique solution (see, e.g., [32]) and therefore there exists an inverse of the operator . This operator can control all initialboundary conditions in each step of HAM.
Therefore we rewrite (1) asfor some appropriate constant , and construct the socalled zeroorder deformation equation [14, 15]:where is an embedding parameter, is a nonzero auxiliary function, is an initial guess, and is an unknown function. The conditions and correspond torespectively. Thus as increases from 0 to 1, the solution varies from the initial guess to the solution [14, 15].
Expanding in Taylor series with respect to , one haswhereIf the auxiliary linear operator, the initial guess, the auxiliary parameter , and the auxiliary function are so properly chosen, then series (11) converges at andwhich must be one of solutions of the original equation, as proved by Liao [17].
According to definition (12), the governing equation can be deduced from the zeroorder deformation equation (9). Define the vectorDifferentiating equation (9) times with respect to embedding parameter and then setting , we have the socalled thorder deformation equationwhere
2. Applications
To demonstrate the advantages of our approach first we consider the wavelike linear equation with constant coefficients.
Example 1. We consider the initialboundary value problemThe exact solution is .
The traditional methods (with ) do not work for this kind of problems. For example, the operator can not control the condition in every iteration step, the same for the operator ; that is, we need an operator that can control all initial/boundary conditions. Clearly the most appropriate operator should be the wave operator. We take . Equation (1) suggests that we define the nonlinear operator asUsing the above definition, we construct the zerothorder deformation equationand the th order deformation equationwith the initial/boundary conditions , , and .
Now it follows from the theory of wave equations that [32] the solution of the equation with the same initialboundary conditions isAccording to (21) we now successively obtain and it follows fromthat [32].
Example 2. We consider the initialboundary value problemThe exact solution isFirst we rewrite the equation asThe solution of the equation with the same initialboundary conditions will be taken as an initial approximation:For we haveand therefore and continuing in this way, we obtain and .
The above example demonstrates the importance of the proposed technique with the use of wave operator. In fact, all traditional approaches with some auxiliary operator do not work, since the (exact) solution is not analytical function on the whole region. For example, the approaches used in [23, 24] can not be applied to solve this problem (note that the approaches used in [23, 24] are effective and simple for the problems with analytical solutions and/or for the problems with initial conditions).
Example 3. We consider the nonlinear initialboundary value problemThe exact solution is . We take and as a solution of the problem with the initial/boundary conditions , , and [32]:And now we obtain from (16)where . Hence we obtainfor . The absolute errors between the exact and the twoterm approximation of the series solution for some values of with are shown in Table 1. The exact solution for Example 3 is shown in Figure 1 and the approximate solution is shown in Figure 2. A higher accuracy level can be attained by evaluating some more terms.

Example 4. Consider the nonlinear initialboundary value problemThe exact solution is . We take , and consequentlyThe absolute errors between the exact and the twoterm approximation of the series solution for some values of for are shown in Table 2. The exact solution for Example 4 is shown in the Figure 3 and the approximate solution is shown in Figure 4.

Example 5. Consider the nonlinear initialboundary value problemwhose exact solution isFirst we rewrite the equation asThe solution of the equation with the same initialboundary conditions isFor we havefor , and for the exact solution we obtain .
3. Discussion
The main goal of this work was to propose a reliable method for solving wavelike equations with variable coefficients. The proposed equations may not be solved by the method of separation of variables or by HAM or ADM in standard form. The main difficulty in the use of previous methods is related to the choice of an auxiliary linear operator . Traditional methods work effectively in case of analytical solutions in the whole region, that is, in fact, when some initial/boundary condition is supposed by another one or in case of some special type of equations (homogeneous, etc.).
The proposed method was applied directly without any need for restrictive assumptions, and this gives it a wider applicability. This method is capable of greatly reducing the volume of computational work compared to standard approaches while still maintaining high accuracy of the approximate solution. A higher accuracy level can be attained even by evaluating some twothree terms in the series solution.
The approach was tested by employing the method to obtain solutions for several problems. The results obtained in all cases demonstrate the efficiency of this approach.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
References
 P. J. Caudrey, J. C. Eilbeck, and J. D. Gibbon, “The sineGordon equation as a model classical field theory,” Il Nuovo Cimento B, vol. 25, no. 2, pp. 497–512, 1975. View at: Publisher Site  Google Scholar
 R. K. Dodd, I. C. Eilbeck, and J. D. Gibbon, Solitons and Nonlinear Wave Equations, Academic Press, London, UK, 1982.
 J. R. Holliday, J. B. Rundle, K. F. Tiampo, W. Klein, and A. Donnellan, “Modification of the pattern informatics method for forecasting large earthquake events using complex eigenfactors,” Tectonophysics, vol. 413, no. 12, pp. 87–91, 2006. View at: Publisher Site  Google Scholar
 A. A. Akhmetov, “Long current loops as regular solutions of the equation for coupling currents in a flat twolayer superconducting cable,” Cryogenics, vol. 43, no. 3–5, pp. 317–322, 2003. View at: Publisher Site  Google Scholar
 G. D. Manolis and T. V. Rangelov, “Nonhomogeneous elastic waves in soils: notes on the vector decomposition technique,” Soil Dynamics and Earthquake Engineering, vol. 26, no. 10, pp. 952–959, 2006. View at: Publisher Site  Google Scholar
 A. R. Manwell, The Tricomi Equation with Applications to the Theory of Plane Transonic Flow, Pitman, Marshfield, Mass, USA, 1979.
 A. V. Bitsadze, Some Classes of Partial Differential Equations, Gordon & Breach, 1988, (Russian).
 L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Pergamon Press, 1982.
 A.M. Wazwaz, “The tanh method: exact solutions of the sineGORdon and the sinhGORdon equations,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1196–1210, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 D. Kaya, “A numerical solution of the sineGordon equation using the modified decomposition method,” Applied Mathematics and Computation, vol. 143, no. 23, pp. 309–317, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 A. Aslanov, “The homotopyperturbation method for solving kleingordontype equations with unbounded righthand side,” Zeitschrift für Naturforschung Section A, vol. 64, no. 12, pp. 149–152, 2009. View at: Google Scholar
 S. M. ElSayed, “The decomposition method for studying the KleinGordon equation,” Chaos, Solitons & Fractals, vol. 18, no. 5, pp. 1025–1030, 2003. View at: Publisher Site  Google Scholar
 A.M. Wazwaz and A. Gorguis, “Exact solutions for heatlike and wavelike equations with variable coefficients,” Applied Mathematics and Computation, vol. 149, no. 1, pp. 15–29, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, Shanghai, China, 1992.
 S.J. Liao, “An approximate solution technique not depending on small parameters: a special example,” International Journal of NonLinear Mechanics, vol. 30, no. 3, pp. 371–380, 1995. View at: Publisher Site  Google Scholar
 S.J. Liao, “A kind of approximate solution technique which does not depend upon small parameters—II. An application in fluid mechanics,” International Journal of NonLinear Mechanics, vol. 32, no. 5, pp. 815–822, 1997. View at: Publisher Site  Google Scholar
 S. J. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, Chapman & Hall, CRC, Boca Raton, Fla, USA, 2003.
 S. J. Liao, “An optimal homotopyanalysis approach for strongly nonlinear differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 8, pp. 2003–2016, 2010. View at: Publisher Site  Google Scholar
 T. M. Hayat and M. Sajid, “On analytic solution for thin film flow of a fourth grade fluid down a vertical cylinder,” Physics Letters A: General, Atomic and Solid State Physics, vol. 361, no. 45, pp. 316–322, 2007. View at: Publisher Site  Google Scholar
 S. Abbasbandy, “The application of homotopy analysis method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 360, no. 1, pp. 109–113, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 A. Aslanov, “Homotopy perturbation method for solving wavelike nonlinear equations with initialboundary conditions,” Discrete Dynamics in Nature and Society, vol. 2011, Article ID 534165, 10 pages, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 R. Rajaraman, “Analytical solutions for some of the nonlinear hyperboliclike equations with variable conditions,” Global Journal of Science Frontier Research Mathematics and Decision Sciences, vol. 12, no. 5, pp. 54–60, 2012. View at: Google Scholar
 A. K. Alomari, M. S. M. Noorani, and R. Nazar, “Solutions of heatlike and wavelike equations with variable coefficients by means of the homotopy analysis method,” Chinese Physics Letters, vol. 25, no. 2, pp. 589–592, 2008. View at: Publisher Site  Google Scholar
 T. Öziş and D. Ağırseven, “He's homotopy perturbation method for solving heatlike and wavelike equations with variable coefficients,” Physics Letters A, vol. 372, no. 38, pp. 5944–5950, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 E. M. E. Zayed, T. A. Nofal, and K. A. Gepreel, “The homotopy perturbation method for solving nonlinear Burgers and new coupled MKdV equations,” Zeitschrift für Naturforschung A, vol. 63, pp. 627–633, 2008. View at: Google Scholar
 E. M. E. Zayed, T. A. Nofal, and K. A. Gepreel, “Homotopy perturbation and Adomain decomposition methods for solving nonlinear Boussinesq equations,” Communications on Applied Nonlinear Analysis, vol. 15, no. 3, pp. 57–70, 2008. View at: Google Scholar  MathSciNet
 E. M. E. Zayed, T. A. Nofal, and K. A. Gepreel, “On using the homotopy perturbation method for finding the traveling wave solutions of generalized nonlinear HirotaSatsuma coupled KdV equations,” International Journal of Nonlinear Science, vol. 7, no. 2, pp. 159–166, 2009. View at: Google Scholar
 E. M. Zayed, T. A. Nofal, and K. A. Gepreel, “The travelling wave solutions for nonlinear initialvalue problems using the homotopy perturbation method,” Applicable Analysis, vol. 88, no. 4, pp. 617–634, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 E. M. E. Zayed and H. M. AbdelRahman, “The homotopy perturbation method and the Adomian decomposition method for the nonlinear coupled equations,” Journal of Partial Differential Equations, vol. 22, no. 4, pp. 334–351, 2009. View at: Publisher Site  Google Scholar
 E. M. E. Zayed and H. M. Abdel Rahman, “The variational iteration method and the variational homotopy perturbation method for solving the KdVBurgers equation and the SharmaTassoOlver equation,” Zeitschrift für Naturforschung A, vol. 65, no. 1, pp. 25–33, 2010. View at: Google Scholar
 E. M. E. Zayed and H. M. Abdel Rahman, “The homotopy analysis method for solving the nonlinear evolution equations in mathematical physics,” Communications on Applied Nonlinear Analysis, vol. 18, no. 3, pp. 53–70, 2011. View at: Google Scholar  Zentralblatt MATH
 W. A. Strauss, Partial Differential Equations, John Wiley & Sons, New York, NY, USA, 1992. View at: MathSciNet
Copyright
Copyright © 2015 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.