AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 10.1155/2015/628310 628310 Research Article A Homotopy-Analysis Approach for Nonlinear Wave-Like Equations with Variable Coefficients http://orcid.org/0000-0001-5858-9848 Aslanov Afgan Naz R. Mathematics and Computing Department, Beykent University, 34396 Istanbul Turkey beykent.edu.tr 2015 1102015 2015 20 02 2015 23 04 2015 1102015 2015 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.

We are interested in the approximate analytical solutions of the wave-like 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 equation(1)utt-kx,tuxx=fu,ux,uy,uxt+gx,t,x>0,t>0with initial conditions(2)ux,0=φx,utx,0=ψxand boundary condition(3)u0,t=ht,where f, g, k, φ, ψ, and h are known functions.

Note that the proposed method can be applied for equations like (4)kx,tutt-uxx=fu,ux,uy,uxt+gx,t,with the same type of initial-boundary conditions.

Problems like (1)–(3) model many problems in classical and quantum mechanics, solitons, and matter physics [1, 2]. If f is a function of u only and k(x,t)=const we obtain Klein-Gordon or sine-Gordon type equations. These models can describe some nonlinear phenomena; for example, wave-like equation can describe earthquake stresses , coupling currents in a flat multistrand two-layer superconducting cable , and nonhomogeneous elastic waves in soils . Typical examples of the wave-like equations with variable coefficients are Euler-Tricomi equation  or Chaplygin equation  given by(5)utt-tuxx=0,utt-ktuxx=0,which are useful in the study of transonic flow, where u=u(t,x) is the stream function of a plane-parallel steady-state gas flow, k(t) is positive at subsonic and negative at supersonic speed, and x is the angle of inclination of the velocity vector. Some Chaplygin type of equation of the special form(6)utt+x21-x2/c2uxx+xux=0,where c=c(x) is the speed of sound, has also applications in the study of transonic flow .

Note that we use the term wave-like equation to describe the partial differential equations with the terms utt and uxx; that is, the term “wave-like” may not correspond to the real physical waves, in general.

Recently, there has been a growing interest for obtaining the explicit solutions to wave-like and heat-like models by analytic techniques. Wazwaz  used the tanh method to obtain the exact solution of the sine-Gordon equation. Kaya  applied the modified decomposition method to solve the sine-Gordon equation. Aslanov  used homotopy perturbation method to solve Klein-Gordon type of equations with unbounded right-hand side. El-Sayed  and Wazwaz and Gorguis  used Adomian decomposition method for solving wave-like and heat-like problems.

The homotopy analysis method  is developed to search the accurate asymptotic solutions of nonlinear problems. Liao  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  and Abbasbandy  pointed out that the homotopy perturbation method is only a special case of the HAM. Aslanov  used homotopy perturbation method to solve wave-like equations with initial-boundary conditions. Rajaraman  and Alomari at al.  applied HAM for solving nonlinear equations with initial conditions. Öziş and Ağırseven  used homotopy perturbation method for solving heat-like and wave-like 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 , the variational iteration method , homotopy analysis method , and others.

Here we will further extend the applications of HAM to obtain an approximate series solution for the nonlinear wave-like equations with variable coefficients and with initial-boundary conditions. The difficulty in the use of standard HAM is that the choice of the linear operator L in standard form (like Lu=utt or Lu=uxx) 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 initial-boundary problems of types (1)–(3). The central idea here is that the problem,(7)utt-c2uxx=gx,t,ux,0=φx,utx,0=ψx,u0,t=ht,has a unique solution (see, e.g., ) and therefore there exists an inverse of the operator L:utt-c2uxx. This operator can control all initial-boundary conditions in each step of HAM.

Therefore we rewrite (1) as(8)Lu=-c2uxx+kx,tuxx+fu,ux,uy,uxy+gx,t,for some appropriate constant c, and construct the so-called zero-order deformation equation [14, 15]:(9)1-qLϕx,t;q-u0x,t=qNϕx,t;q,where q[0,1] is an embedding parameter, is a nonzero auxiliary function, L(u)=utt-c2uxx,N(ϕ)=ϕtt-k(x,t)ϕxx-f(ϕ,ϕx,ϕy,ϕxt)-g(x,t),u0(x,t) is an initial guess, and ϕ(x,t;q) is an unknown function. The conditions q=0 and q=1 correspond to(10)ϕx,t;0=u0x,t,ϕx,t;1=ux,t,respectively. Thus as q increases from 0 to 1, the solution ϕ(x,t;q) varies from the initial guess u0(x,t) to the solution u(x,t) [14, 15].

Expanding ϕ(x,t;q) in Taylor series with respect to q, one has(11)ϕx,t;q=u0x,t+m=1umx,tqm,where(12)umx,t=1m!mϕx,t,qqmq=0.If the auxiliary linear operator, the initial guess, the auxiliary parameter , and the auxiliary function are so properly chosen, then series (11) converges at q=1 and(13)ϕx,t;1=u0x,t+m=1umx,t,which must be one of solutions of the original equation, as proved by Liao .

According to definition (12), the governing equation can be deduced from the zero-order deformation equation (9). Define the vector(14)un=u0,u1,,un.Differentiating equation (9) m times with respect to embedding parameter q and then setting q=0, we have the so-called mth-order deformation equation(15)Lum-χmum-1=Rmum-1,where(16)Rmum-1=1m-1!m-1Nϕx,t,qqm-1q=0,(17)χm=0,m11,m>1.

2. Applications

To demonstrate the advantages of our approach first we consider the wave-like linear equation with constant coefficients.

Example 1.

We consider the initial-boundary value problem(18)utt+uxx=-2uxt,x>0,t>0,ux,0=x3,utx,0=-3x2,u0,t=t3.The exact solution is u=(x-t)3.

The traditional methods (with Lu=utt) do not work for this kind of problems. For example, the operator Lu=utt can not control the condition u(0,t)=0 in every iteration step, the same for the operator Lu=uxx; 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 L(ϕ)=2ϕ(x,t,q)/t2-2ϕ(x,t,q)/x2. Equation (1) suggests that we define the nonlinear operator as(19)Nϕ=22ϕx,t,qx2+22ϕx,t,qxt.Using the above definition, we construct the zeroth-order deformation equation(20)1-qLϕx,t;q-u0x,t=qNϕx,t;qand the mth order deformation equation(21)Lum-χmum-1=Rmum-1,with the initial/boundary conditions um(x,0)=0,  (um)t(x,0)=0, and um(0,t)=0.

Now it follows from the theory of wave equations that  the solution of the equation utt-uxx=0 with the same initial-boundary conditions is(22)u0x,t=12x+t3+x-t3-12x-tx+t3y2dy=x-t3forx>t,u0x,t=12x+t3+x-t3-12t-xx+t3y2dy+t-x3=t-x3,forx<t.According to (21) we now successively obtain (23)Lu1=2u0xx+2u0xtgx,t=0and it follows from(24)u1=20tx-t+sx+t-sgy,sdyds,forx>t,u1=20t-xt-x-sx+t-sgy,sdyds+2t-xtx-t+sx+t-sgy,sdyds,forx<tthat u1=u2==0 .

Example 2.

We consider the initial-boundary value problem(25)utt=6-x-t26uxx+u,x>0,t>0,ux,0=x3,utx,0=-3x2,u0,t=0.The exact solution is(26)u=x-t3ifx>t0ifx<t.First we rewrite the equation as(27)utt-uxx=-x-t26uxx+u.The solution of the equation utt-uxx=0 with the same initial-boundary conditions will be taken as an initial approximation:(28)u0=12x+t3+x-t3-12x-tx+t3y2dy=x-t3forx>t,u0=12x+t3-t-x3-12t-xx+t3y2dy=0forx<t.For u1 we have(29)Lu1=x-t26u0xx-u0and therefore (30)u1=20tx-t+sx+t-sy-s266y-s-y-s3dyds=0,forx>t,u1=20t-xt-x-sx+t-s0dyds+2t-xtx-t+sx+t-s0dyds=0,forx<t,and continuing in this way, we obtain u2=u3==0 and uexact=u0+u1.

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 Lu=utt 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 initial-boundary value problem(31)utt-x2uxx=xet-ux2+e2t,ux,0=x,utx,0=x,u0,t=0.The exact solution is u=xet. We take L(u)=utt-uxx and u0 as a solution of the problem L(u)=0 with the initial/boundary conditions u(x,0)=x, ut(x,0)=x, and u(0,t)=0 :(32)u0=12x+t+x-t+12x-tx+tydy=xt+1forx>t,u0=12x+t-t-x+12t-xx+tydy=xt+1for  x<t.And now we obtain from (16)(33)Lu1=--1-x2u0xx+xet-u0x2+e2t=-xet-t+12+e2t,u1=-20tx-t+sx+t-sgy,sdyds,forx>t,u1=-20t-xt-x-sx+t-sgy,sdyds-2t-xtx-t+sx+t-sgy,sdydsforx<t,where g(y,s)=yes-(s+1)2+e2s. Hence we obtain(34)u1=-xet-x-tx-12t-12t2-13t3-112t4+14e2t-14,forx>t,u1=-tx2-3x2-t2x-tx33-t3x3+t2x22-e2t4e2x-2tx+xet+x22-x33+x412+e2t4for x<t. The absolute errors between the exact and the two-term approximation of the series solution for some values of (x,t)[0,1]×[0,1] with =-1 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.

Maximum errors for Example 3.

t x 0.1 0.2 0.3 0.4 0.5 0.6 0.8 0.9
0.1 9.0 · 10−6 9.0 · 10−6 9.0 · 10−6 9.0 · 10−6 9.0 · 10−6 9.0 · 10−6 9.0 · 10−6 9.0 · 10−6
0.2 1.4 · 10−4 1.5 · 10−4 1.5 · 10−4 1.5 · 10−4 1.5 · 10−4 1.5 · 10−4 1.5 · 10−4 1.5 · 10−4
0.4 2.0 · 10−3 2.7 · 10−3 2.9 · 10−3 2.9 · 10−3 2.9 · 10−3 2.9 · 10−3 2.9 · 10−3 2.9 · 10−3
0.5 4.7 · 10−3 6.8 · 10−3 7.5 · 10−3 7.6 · 10−3 7.6 · 10−3 7.6 · 10−3 7.6 · 10−3 7.6 · 10−3
0.7 1.7 · 10−2 2.6 · 10−2 3.1 · 10−2 3.3 · 10−2 3.4 · 10−2 3.4 · 10−2 3.4 · 10−2 3.4 · 10−2
0.8 2.9 · 10−2 4.6 · 10−2 5.5 · 10−2 6.0 · 10−2 6.2 · 10−2 6.3 · 10−2 6.3 · 10−2 6.3 · 10−2
0.9 4.6 · 10−2 7.5 · 10−2 9.2 · 10−2 0.1020 0.10682 0.10888 0.10973 0.10974

The exact solution for Example 3.

The approximate solution for Example 3.

Example 4.

Consider the nonlinear initial-boundary value problem(35)utt-tuxx=u2-t2x4-2t2x3-t2x2-2t2,ux,0=0,utx,0=x+x2,u0,t=0.The exact solution is u=xt+x2t. We take L=utt-uxx, and consequently(36)u0=12x-tx+ty+y2dy=13tt2+3x2+3x,forx>t,u0=12t-xx+ty+y2dy=13x3t2+3t+x2,forx<t,u1=-1180t8+245t6x2+245t6x-23t3,forx>t,u1=--7945t7x+649t6x2+145t6x-47345t5x3+25336t4x4-89t4x3-7345t3x5+236t3x4-13t3x-145t2x6-356t2x5+t2x2-t2x+145tx7+7118tx6+103tx3-1360x8-125126x7-236x4,forx<t.The absolute errors between the exact and the two-term approximation of the series solution for some values of (x,t)[0,1]×[0,1] for =-1 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.

Maximum errors for Example 4.

t x 0.1 0.2 0.3 0.4 0.5 0.6 0.8 0.9
0.1 3.2 · 10−4 3.3 · 10−4 3.3 · 10−4 3.3 · 10−4 3.3 · 10−4 3.3 · 10−4 3.3 · 10−4 3.3 · 10−4
0.2 1.2 · 10−3 2.6 · 10−3 2.6 · 10−3 2.6 · 10−3 2.6 · 10−3 2.6 · 10−3 2.6 · 10−3 2.6 · 10−3
0.4 3.3 · 10−3 6.6 · 10−3 1.4 · 10−2 2.1 · 10−2 2.1 · 10−2 2.1 · 10−2 2.1 · 10−2 2.1 · 10−2
0.5 5.5 · 10−3 8.7 · 10−3 1.2 · 10−2 2.6 · 10−2 4.0 · 10−2 4.1 · 10−2 4.0 · 10−2 4.0 · 10−2
0.7 1.8 · 10−2 2.2 · 10−2 1.5 · 10−2 4.8 · 10−3 8.9 · 10−3 5.2 · 10−2 0.10648 0.10507
0.8 3.7 · 10−2 4.4 · 10−2 2.9 · 10−2 1.9 · 10−3 2.0 · 10−2 1.4 · 10−2 0.19175 0.14981
0.9 0.07335 9.0 · 10−2 6.5 · 10−2 1.5 · 10−2 4.1 · 10−2 8.0 · 10−2 1.6 · 10−2 0.1997

Exact solution for Example 4.

Approximate two-term approximation for Example 4.

Example 5.

Consider the nonlinear initial-boundary value problem(37)utt=6+x-t56uxx-u2,x>0,t>0,ux,0=x3,utx,0=-3x2,u0,t=t3,whose exact solution is(38)u=x-t3ifx>t,t-x3ifx<t.First we rewrite the equation as(39)utt-uxx=x-t56uxx-u2.The solution of the equation utt-uxx=0 with the same initial-boundary conditions is(40)u0=12x+t3+x-t3-12x-tx+t3y2dy=x-t3,forx>t,(41)u0=12x+t3-t-x3-12t-xx+t3y2dy+t-x3=t-x3,forx<t.For u1 we have(42)Lu1=-x-t56u0xx-u02,u1=-20tx-t+sx+t-sy-s566y-s-y-s6dyds=0,forx>t,u1=-20t-xt-x-sx+t-ss-y566s-y-s-y6dyds-2t-xtx-t+sx+t-ss-y566s-y-s-y6dyds=0for x<t, and for the exact solution we obtain u=u0+u1.

3. Discussion

The main goal of this work was to propose a reliable method for solving wave-like 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 L. 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 two-three 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.

Caudrey P. J. Eilbeck J. C. Gibbon J. D. The sine-Gordon equation as a model classical field theory Il Nuovo Cimento B 1975 25 2 497 512 10.1007/bf02724733 2-s2.0-51249191762 Dodd R. K. Eilbeck I. C. Gibbon J. D. Solitons and Nonlinear Wave Equations 1982 London, UK Academic Press Holliday J. R. Rundle J. B. Tiampo K. F. Klein W. Donnellan A. Modification of the pattern informatics method for forecasting large earthquake events using complex eigenfactors Tectonophysics 2006 413 1-2 87 91 10.1016/j.tecto.2005.10.008 2-s2.0-30944453549 Akhmetov A. A. Long current loops as regular solutions of the equation for coupling currents in a flat two-layer superconducting cable Cryogenics 2003 43 3–5 317 322 10.1016/s0011-2275(03)00028-6 2-s2.0-0038298713 Manolis G. D. Rangelov T. V. Non-homogeneous elastic waves in soils: notes on the vector decomposition technique Soil Dynamics and Earthquake Engineering 2006 26 10 952 959 10.1016/j.soildyn.2006.01.024 2-s2.0-33746710898 Manwell A. R. The Tricomi Equation with Applications to the Theory of Plane Transonic Flow 1979 Marshfield, Mass, USA Pitman Bitsadze A. V. Some Classes of Partial Differential Equations 1988 Gordon & Breach (Russian) Landau L. D. Lifshitz E. M. Fluid Mechanics 1982 Pergamon Press Wazwaz A.-M. The tanh method: exact solutions of the sine-GORdon and the sinh-GORdon equations Applied Mathematics and Computation 2005 167 2 1196 1210 10.1016/j.amc.2004.08.005 MR2169761 2-s2.0-25144450080 Kaya D. A numerical solution of the sine-Gordon equation using the modified decomposition method Applied Mathematics and Computation 2003 143 2-3 309 317 10.1016/s0096-3003(02)00363-6 MR1981698 2-s2.0-0037507444 Aslanov A. The homotopy-perturbation method for solving klein-gordon-type equations with unbounded right-hand side Zeitschrift für Naturforschung Section A 2009 64 1-2 149 152 2-s2.0-61449170106 El-Sayed S. M. The decomposition method for studying the Klein-Gordon equation Chaos, Solitons & Fractals 2003 18 5 1025 1030 10.1016/s0960-0779(02)00647-1 2-s2.0-0038577058 Wazwaz A.-M. Gorguis A. Exact solutions for heat-like and wave-like equations with variable coefficients Applied Mathematics and Computation 2004 149 1 15 29 10.1016/S0096-3003(02)00946-3 MR2030981 ZBL1038.65103 2-s2.0-0346119929 Liao S. J. The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis] 1992 Shanghai, China Shanghai Jiao Tong University Liao S.-J. An approximate solution technique not depending on small parameters: a special example International Journal of Non-Linear Mechanics 1995 30 3 371 380 10.1016/0020-7462(94)00054-e 2-s2.0-0000615571 Liao S.-J. A kind of approximate solution technique which does not depend upon small parameters—II. An application in fluid mechanics International Journal of Non-Linear Mechanics 1997 32 5 815 822 10.1016/s0020-7462(96)00101-1 2-s2.0-0031232480 Liao S. J. Beyond Perturbation: Introduction to the Homotopy Analysis Method 2003 Boca Raton, Fla, USA Chapman & Hall, CRC Liao S. J. An optimal homotopy-analysis approach for strongly nonlinear differential equations Communications in Nonlinear Science and Numerical Simulation 2010 15 8 2003 2016 10.1016/j.cnsns.2009.09.002 2-s2.0-76049117931 Hayat T. M. Sajid M. 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 2007 361 4-5 316 322 10.1016/j.physleta.2006.09.060 2-s2.0-33845432502 Abbasbandy S. The application of homotopy analysis method to nonlinear equations arising in heat transfer Physics Letters A 2006 360 1 109 113 10.1016/j.physleta.2006.07.065 MR2288118 2-s2.0-33750339270 Aslanov A. Homotopy perturbation method for solving wave-like nonlinear equations with initial-boundary conditions Discrete Dynamics in Nature and Society 2011 2011 10 534165 10.1155/2011/534165 MR2838432 2-s2.0-84855406707 Rajaraman R. Analytical solutions for some of the nonlinear hyperbolic-like equations with variable conditions Global Journal of Science Frontier Research Mathematics and Decision Sciences 2012 12 5 54 60 Alomari A. K. Noorani M. S. M. Nazar R. Solutions of heat-like and wave-like equations with variable coefficients by means of the homotopy analysis method Chinese Physics Letters 2008 25 2 589 592 10.1088/0256-307x/25/2/064 2-s2.0-42649113047 Öziş T. Ağırseven D. He's homotopy perturbation method for solving heat-like and wave-like equations with variable coefficients Physics Letters A 2008 372 38 5944 5950 10.1016/j.physleta.2008.07.060 MR2458895 2-s2.0-50049115246 Zayed E. M. E. Nofal T. A. Gepreel K. A. The homotopy perturbation method for solving nonlinear Burgers and new coupled MKdV equations Zeitschrift für Naturforschung A 2008 63 627 633 Zayed E. M. E. Nofal T. A. Gepreel K. A. Homotopy perturbation and Adomain decomposition methods for solving nonlinear Boussinesq equations Communications on Applied Nonlinear Analysis 2008 15 3 57 70 MR2441806 2-s2.0-51349121599 Zayed E. M. E. Nofal T. A. Gepreel K. A. On using the homotopy perturbation method for finding the traveling wave solutions of generalized nonlinear Hirota-Satsuma coupled KdV equations International Journal of Nonlinear Science 2009 7 2 159 166 Zayed E. M. Nofal T. A. Gepreel K. A. The travelling wave solutions for non-linear initial-value problems using the homotopy perturbation method Applicable Analysis 2009 88 4 617 634 10.1080/00036810902943604 MR2541144 2-s2.0-67651229170 Zayed E. M. E. Abdel-Rahman H. M. The homotopy perturbation method and the Adomian decomposition method for the nonlinear coupled equations Journal of Partial Differential Equations 2009 22 4 334 351 10.4208/jpde.v22.n4.3 Zayed E. M. E. Abdel Rahman H. M. The variational iteration method and the variational homotopy perturbation method for solving the KdV-Burgers equation and the Sharma-Tasso-Olver equation Zeitschrift für Naturforschung A 2010 65 1 25 33 2-s2.0-77955964657 Zayed E. M. E. Abdel Rahman H. M. The homotopy analysis method for solving the nonlinear evolution equations in mathematical physics Communications on Applied Nonlinear Analysis 2011 18 3 53 70 ZBL1236.35161 2-s2.0-80052008466 Strauss W. A. Partial Differential Equations 1992 New York, NY, USA John Wiley & Sons MR1159712