AAAAbstract and Applied Analysis1687-04091085-3375Hindawi Publishing Corporation85761210.1155/2012/857612857612Research ArticleHomotopy Perturbation Method with an Auxiliary TermHeJi-Huan1,2NoorMuhammad Aslam1National Engineering Laboratory for Modern SilkCollege of Textile and EngineeringSoochow University199 Ren-ai Road, Suzhou 215123Chinasuda.edu.cn2Modern Textile InstituteDonghua UniversityShanghai 200051Chinadhu.edu.cn20121222012201202102011171120112012Copyright © 2012 Ji-Huan He.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.

The two most important steps in application of the homotopy perturbation method are to construct a suitable homotopy equation and to choose a suitable initial guess. The homotopy equation should be such constructed that when the homotopy parameter is zero, it can approximately describe the solution property, and the initial solution can be chosen with an unknown parameter, which is determined after one or two iterations. This paper suggests an alternative approach to construction of the homotopy equation with an auxiliary term; Dufing equation is used as an example to illustrate the solution procedure.

1. Introduction

The homotopy perturbation method  has been worked out over a number of years by numerous authors, and it has matured into a relatively fledged theory thanks to the efforts of many researchers, see Figure 1. For a relatively comprehensive survey on the concepts, theory, and applications of the homotopy perturbation method, the reader is referred to review articles .

Number of publications on homotopy perturbation according to web of science, August 20, 2011.

The homotopy perturbation method has been shown to solve a large class of nonlinear differential problems effectively, easily, and accurately; generally one iteration is enough for engineering applications with acceptable accuracy, making the method accessible to nonmathematicians.

In case of higher-order approximates needed, we can use parameter-expansion technology ; in this paper, we suggest an alternative approach by adding a suitable term in the homotopy equation.

2. Homotopy Perturbation Method

Consider a general nonlinear equationLu+Nu=0, where L and N are, respectively, the linear operator and nonlinear operator.

The first step for the method is to construct a homotopy equation in the form L̃u+p(Lu-L̃u+Nu)=0, where L̃ is a linear operator with a possible unknown constant and L̃u=0 can approximately describe the solution property. The embedding parameter p monotonically increases from zero to unit as the trivial problem (L̃u=0) is continuously deformed to the original one (Lu+Nu=0).

For example, consider a nonlinear oscillatoru+εu3=0,u(0)=A,u(0)=0. For an oscillator, we can use sine or cosine function. We assume that the approximate solution of (2.3) isu(t)=Acosωt, where ω is the frequency to be determined later. We, therefore, can chooseL̃u=ü+ω2u. Accordingly, we can construct a homotopy equation in the formü+ω2u+p(u3-ω2u)=0. When p=0, we haveü+ω2u=0,u(0)=A,u(0)=0, which describes the basic solution property of the original nonlinear equation, (2.3).

When p=1, (2.6) becomes the original one. So the solution procedure is to deform from the initial solution, (2.4), to the real one. Due to one unknown parameter in the initial solution, only one iteration is enough. For detailed solution procedure, refer to .

If a higher-order approximate solution is searched for, we can construct a homotopy equation in the formü+0u+pu3=0. We expand the solution and the coefficient, zero, of the linear term into a series in p:u=u0+pu1+p2u2+,0=ω2+pa1+p2a2+, where the unknown constant, ai, is determined in the (i+1)th iteration. The solution procedure is given in .

3. Homotopy Equation with an Auxiliary Term

In this paper, we suggest an alternative approach to construction of homotopy equation, which isL̃u+p(Lu-L̃u+Nu)+αp(1-p)u=0, where α is an auxiliary parameter. When α=0, (3.1) turns out to be that of the classical one, expressed in (2.2). The auxiliary term, αp(1-p)u, vanishes completely when p=0 or p=1; so the auxiliary term will affect neither the initial solution (p=0) nor the real solution (p=1). The homotopy perturbation method with an auxiliary term was first considered by Noor .

To illustrate the solution procedure, we consider a nonlinear oscillator in the formd2udt2+bu+cu3=0,u(0)=A,u(0)=0, where b and c are positive constants.

Equation (3.2) admits a periodic solution, and the linearized equation of (3.2) isu′′+ω2u=0,u(0)=A,u(0)=0, where ω is the frequency of (3.2).

We construct the following homotopy equation with an auxiliary term:u′′+ω2u+p[(b-ω2)u+cu3]+αp(1-p)u=0. Assume that the solution can be expressed in a power series in p as shown in (2.9).

Substituting (2.9) into (3.4), and processing as the standard perturbation method, we haveu0′′+ω2u0=0,u0(0)=A,u0(0)=0,u1′′+ω2u1+(b-ω2)u0+cu03+αu0=0,u1(0)=0,u1(0)=0,u2′′+ω2u2+(b-ω2)u1+3cu02u1+α(u1-u0)=0, with initial conditionsi=0ui(0)=A,i=0ui(0)=0. Solving (3.5), we haveu0=Acosωt.

Substituting u0 into (3.6) results intou1′′+ω2u1+A(α+b-ω2+34cA2)cosωt+14cA3cos3ωt=0.

Eliminating the secular term needsα+b-ω2+34cA2=0. A special solution of (3.10) isu1=-cA332ω2cos3ωt. If only a first-order approximate solution is enough, we just set α=0; this results inω=b+34cA2.

The accuracy reaches 7.6% even for the case cA2.

The solution procedure continues by submitting u1 into (3.7); after some simple calculation, we obtainu2′′+ω2u2-(αA+3c2A5128ω2)cosωt-(cA3(b-ω2)32ω2+3c2A564ω2+αcA332ω2)cos3ωt-3c2A5128ω2cos5ω=0. No secular term in u2 requiresαA+3c2A5128ω2=0. Solving (3.11) and (3.15) simultaneously, we obtainω=b+(3/4)cA2+(b+(3/4)cA2)2+(3/32)c2A42, and the approximate solution is u(t)=Acosωt, where ω is given in (3.16).

In order to compare with the perturbation solution and the exact solution, we set b=1. In case c1, (3.16) agrees with that obtained by classical perturbation method; when c, we havelimcω=(3/4)+(3/4)2+3/322cA2=0.8832cA2. The exact period readsTex=41+cA20π/2dx1-ksin2x, where k=cA2/2(1+cA2).

In case c, we havelimcTex=6.743cA2,ωex2π6.743cA2. Comparing between (3.17) and (3.20), we find that the accuracy reaches 5.5%, while accuracy of the first-order approximate frequency is 7.6%.

If a higher-order approximate solution is needed, we rewrite the homotopy equation in the formu′′+ω2u+p[(b-ω2)u+cu3]+1p(1-p)u=0. The coefficient, 1, in the auxiliary term, is also expanded in a series in p in the form1=α0+pα1+p2α2+, where αi is identified in the (i+2)th iteration. The solution procedure is similar to that illustrated above.

4. Discussions and Conclusions

Generally the homotopy equation can be constructed in the formL̃u+p(Lu-L̃u+Nu)+1f(p)g(p)h(u,u,u′′,)=0, where f and g are functions of p, satisfying f(0)=0 and g(1)=0, and h can be generally expressed in the formh=u+β1u+β2u′′+. For example, for the Blasius equationu′′′(η)+12u(η)u′′(η)=0,u(0)=u(0)=0,u(+)=1, where the superscript denotes derivative with respect to η, we can construct a homotopy equation in the form u′′′+au′′+p(12uu′′-au′′)+bp(1-p)u′′=0, where a and b are unknown constants to be determined.

The operator L̃() can be also a nonlinear one, for example, if we want to search for a solitonary solution, we can choose L̃u=ut+6uux+uxxx.

The homotopy equation can be easily constructed, and the solution procedure is simple. This paper can be considered a standard homotopy perturbation algorithm and can be used as a paradigm for many other applications.

Acknowledgment

The work is supported by PAPD (a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions).

HeJ.-H.Homotopy perturbation techniqueComputer Methods in Applied Mechanics and Engineering19991783-425726210.1016/S0045-7825(99)00018-31711041ZBL0956.70017HeJ.-H.A coupling method of a homotopy technique and a perturbation technique for non-linear problemsInternational Journal of Non-Linear Mechanics2000351374310.1016/S0020-7462(98)00085-71723761HeJ.-H.An elementary introduction to the homotopy perturbation methodComputers & Mathematics with Applications200957341041210.1016/j.camwa.2008.06.0032488613ZBL1165.65374HeJ.-H.Recent development of the homotopy perturbation methodTopological Methods in Nonlinear Analysis20083122052092432078ZBL1159.34333HeJ.-H.New interpretation of homotopy perturbation methodInternational Journal of Modern Physics B200620182561256810.1142/S02179792060348192251264RentoulM.ArielP. D.Extended homotopy perturbation method and the flow past a nonlinearly stretching sheetNonlinear Science Letters A201121730AbdouM. A.Zakharov-Kuznetsov equation by the homotopy analysis method and Hirota's Bilinear methodNonlinear Science Letters B2011199110HeJ.-H.Some asymptotic methods for strongly nonlinear equationsInternational Journal of Modern Physics B200620101141119910.1142/S02179792060337962220565ZBL1102.34039HeJ. H.An elementary introduction to recently developed asymptotic methods and nanomechanics in textile engineeringInternational Journal of Modern Physics B20082221348735782-s2.0-5094911226610.1142/S0217979208048668GanjiD. D.KachapiS. H. H.Analysis of nonlinear equations in fluidsProgress in Nonlinear Science201121293GanjiD. D.KachapiS. H. H.Analytical and numerical methods in engineering and applied sciencesProgress in Nonlinear Science201131579WangS. Q.HeJ. H.Nonlinear oscillator with discontinuity by parameter-expansion methodChaos, Solitons and Fractals20083546886912-s2.0-3474891756110.1016/j.chaos.2007.07.055XuL.Application of He's parameter-expansion method to an oscillation of a mass attached to a stretched elastic wirePhysics Letters A20073683-42592622-s2.0-3454766188710.1016/j.physleta.2007.04.004ShouD. H.HeJ. H.Application of parameter-expanding method to strongly nonlinear oscillatorsInternational Journal of Nonlinear Sciences and Numerical Simulation2007811211242-s2.0-34249996463NoorM. A.Some iterative methods for solving nonlinear equations using homotopy perturbation methodInternational Journal of Computer Mathematics2010871–314114910.1080/002071608019695132598731ZBL1182.65079