《Advances In The Homotopy Analysis Method》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1ADVANCESINTHEHOMOTOPYANALYSISMETHOD
2May2,201314:6BC:8831-ProbabilityandStatisticalTheoryPST˙wsThispageintentionallyleftblank
3ADVANCESINTHEHOMOTOPYANALYSISMETHODEditorShijunLiaoShanghaiJiaoTongUniversity,ChinaWorldScientificNEWJERSEY•LONDON•SINGAPORE•BEIJING•SHANGHAI•HONGKONG•TAIPEI•CHENNAI
4PublishedbyWorldScientificPublishingCo.Pte.Ltd.5TohTuckLink,Singapore596224USAoffice:27WarrenStreet,Suite401-402,Hackensack,NJ07601UKoffice:57SheltonStreet,CoventGarden,LondonWC2H9HELibraryofCongressCataloging-in-PublicationDataAdvancesinthehomotopyanalysismethod/editedbyShijunLiao,professor,deputydirectoroftheStateKeyLabofOceanEngineering,ShanghaiJiaoTongUniversity,China.pagescmIncludesbibliographicalreferences.ISBN978-9814551243(hardcover:alk.paper)1.Homotopytheory.I.Liao,Shijun,1963–editorofcompilation.QA612.7.A3752014514'.24--dc232013028624BritishLibraryCataloguing-in-PublicationDataAcataloguerecordforthisbookisavailablefromtheBritishLibrary.Copyright©2014byWorldScientificPublishingCo.Pte.Ltd.Allrightsreserved.Thisbook,orpartsthereof,maynotbereproducedinanyformorbyanymeans,electronicormechanical,includingphotocopying,recordingoranyinformationstorageandretrievalsystemnowknownortobeinvented,withoutwrittenpermissionfromthePublisher.Forphotocopyingofmaterialinthisvolume,pleasepayacopyingfeethroughtheCopyrightClearanceCenter,Inc.,222RosewoodDrive,Danvers,MA01923,USA.Inthiscasepermissiontophotocopyisnotrequiredfromthepublisher.In-houseEditor:AngelineFongPrintedinSingapore
5October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.0PrefaceThehomotopyisafundamentalconceptintopology,whichcanbetracedbacktoJulesHenriPoincar´e(1854–1912),aFrenchmathematician.Basedonthehomotopy,twomethodshavebeendeveloped.Oneisthehomotopycontinuationmethoddatingbackto1930s,whichisaglobalconvergentnumericalmethodmainlyfornonlinearalgebraicequations.Theotheristhehomotopyanalysismethod(HAM)proposedin1990sbyShijunLiao,theeditorofthisbook,whichisananalyticapproximationmethodwithguaranteeofconvergence,mainlyfornonlineardifferentialequations.Differentfromperturbationtechniqueswhicharestronglydependentuponsmall/largephysicalparameters(i.e.perturbationquantities),theHAMhasnothingtodowithanysmall/largephysicalparametersatall.Besides,manyanalyticapproximationmethods,suchas“Lyapunovartifi-cialsmallparametermethod”,“Adomiandecompositionmethod”andsoon,areonlyspecialcasesoftheHAM.Unlikeotheranalyticapproximationtechniques,theHAMprovidesusgreatfreedomandflexibilitytochooseequation-typeandsolutionexpressionofhigh-orderapproximationequa-tions.Noticethat“theessenceofmathematicsliesentirelyinitsfreedom”,aspointedoutbyGeorgCantor(1845–1918).Mostimportantly,differentfromallofotheranalyticapproximationmethods,theHAMprovidesusaconvenientwaytoguaranteetheconvergenceofapproximationseriesbymeansofintroducingtheso-called“convergence-controlparameter”.Infact,itistheconvergence-controlparameterthatdifferstheHAMfromallotheranalyticapproximationmethods.Asaresult,theHAMisgenerallyvalidforvarioustypesofequationswithhighnonlinearity,especiallyforthosewithoutsmall/largephysicalparameters.Since1992whentheearlyHAMwasfirstproposedbyLiao,theHAMhasbeendevelopinggreatlyintheoryandappliedsuccessfullytonumer-oustypesofnonlinearequationsinlotsofdifferentfieldsbyscientists,researchers,engineersandgraduatedstudentsindozensofcountries.Alloftheseverifytheoriginality,novelty,validityandgeneralityoftheHAM.v
6October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.0viPrefaceSo,itisnecessarytodescribe,althoughbriefly,thecurrentadvancesoftheHAMinboththeoryandapplications.Thisisthefirstmotivationofthebook,whosechaptersarecontributedbytheleadingresearchersintheHAMcomingfromsevencountries.Anytrulynewmethodshouldgivesomethingnoveland/orbetter.Inthepast20years,hundredsofarticlesrelatedtotheHAMwerepublishedinvariousfields,andsomenewsolutionswereindeedfoundbymeansoftheHAM.Thus,itisnowthetimetosuggestsomevaluablebutchallengingnonlinearproblemstotheHAMcommunity.Thisisthesecondmotivationofthebook.Someoftheseproblemsareveryfamous,withalonghistory.Hopefully,theabove-mentionedfreedomandflexibilityoftheHAMmightcreatesomenovelideasandinspirebrave,enterprising,youngresearcherswithstimulatedimaginationtoattackthemwithsatisfactoryresults.IpersonallybelievethattheapplicationsoftheHAMonthesefamous,chal-lengingproblemsmightnotonlyindicatethegreatpotentialoftheHAM,butalsoleadtogreatmodificationsoftheHAMintheory.AbriefreviewoftheHAMisgiveninChapter1,withsomesuggestedchallengingproblems.Thefascinating“PredictorHAM”and“SpectralHAM”aredescribedinChapters2and3,respectively.Someinterestingtheoreticalworksontheauxiliarylinearoperator,convergence-controlpa-rameterandconvergenceofapproximationseriesaredescribedinChapters4and5.AnattractiveapplicationoftheHAMaboutflowsofnanofluidisgiveninChapter6.AcharmingapplicationoftheHAMfortime-fractionalboundary-valueproblemisillustratedinChapter7.TheHAM-basedMaplepackageNOPH1.0.2(http://numericaltank.sjtu.edu.cn/NOPH.htm)forperiodicoscillationsandlimitcyclesofnonlineardynamicsystemswithvariousapplicationsisdescribedinChapter8.TheHAM-basedMathemat-icapackageBVPh2.0(http://numericaltank.sjtu.edu.cn/BVPh.htm)forcouplednonlinearordinarydifferentialequationsanditsapplicationsaregiveninChapter9.Bothofthemareeasy-to-use,user-friendly,andfreeavailableonlinewithuser’sguide.TheycangreatlysimplifysomeapplicationsoftheHAM.Itisagreatpitythatitisimpossibletodescribe,evenbriefly,thewholeadvancesoftheHAMintheoryandapplicationsinsuchabook.Here,IwouldliketoexpressmysincereandtruthfulacknowledgementstoalloftheHAMcommunityfortheirgreatcontributionstotheHAM.ShijunLiaoJune2013,Shanghai
7October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.0ContentsPrefacev1.ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod1S.-J.Liao2.PredictorHomotopyAnalysisMethod(PHAM)35S.AbbasbandyandE.Shivanian3.SpectralHomotopyAnalysisMethodforNonlinearBoundaryValueProblems85S.MotsaandP.Sibanda4.StabilityofAuxiliaryLinearOperatorandConvergence-ControlParameter123R.A.VanGorder5.AConvergenceConditionoftheHomotopyAnalysisMethod181M.Turkyilmazoglu6.HomotopyAnalysisMethodforSomeBoundaryLayerFlowsofNanofluids259T.HayatandM.Mustafavii
8October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.0viiiContents7.HomotopyAnalysisMethodforFractionalSwift–HohenbergEquation291S.DasandK.Vishal8.HAM-BasedPackageNOPHforPeriodicOscillationsofNonlinearDynamicSystems309Y.-P.Liu9.HAM-BasedMathematicaPackageBVPh2.0forNonlinearBoundaryValueProblems361Y.-L.ZhaoandS.-J.Liao
9October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1Chapter1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethodShijunLiaoShanghaiJiaoTongUniversity,Shanghai200240,Chinasjliao@sjtu.edu.cnAbriefreviewofthehomotopyanalysismethod(HAM)andsomeofitscurrentadvancesaredescribed.Weemphasizethattheintroductionofthehomotopy,abasicconceptintopology,isamilestoneofthean-alyticapproximationmethods,sinceitisthehomotopywhichprovidesusgreatfreedomandflexibilitytochooseequationtypeandsolutionex-pressionofhigh-orderapproximationequations.Besides,theso-called“convergence-controlparameter”isamilestoneoftheHAM,too,sinceitistheconvergence-controlparameterthatprovidesusaconvenientwaytoguaranteetheconvergenceofsolutionseriesandthatdifferstheHAMfromallotheranalyticapproximationmethods.RelationsoftheHAMtothehomotopycontinuationmethodandotheranalyticapproximationtechniquesarebrieflydescribed.Someinterestingbutchallengingnon-linearproblemsaresuggestedtotheHAMcommunity.AspointedoutbyGeorgCantor(1845–1918),“theessenceofmathematicsliesentirelyinitsfreedom”.Hopefully,theabove-mentionedfreedomandgreatflex-ibilityoftheHAMmightcreatesomenovelideasandinspirebrave,enterprising,youngresearcherswithstimulatedimaginationtoattackthemwithsatisfactory,betterresults.1
10October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.12S.-J.LiaoContents1.1.Background.....................................21.2.AbriefhistoryoftheHAM............................51.3.SomeadvancesoftheHAM............................101.3.1.Generalizedzeroth-orderdeformationequation..............101.3.2.SpectralHAMandcomplicatedauxiliaryoperator............131.3.3.PredictorHAMandmultiplesolutions..................151.3.4.ConvergenceconditionandHAM-basedsoftware.............161.4.Relationshipstoothermethods..........................181.5.Chanceandchallenge:somesuggestedproblems................201.5.1.Periodicsolutionsofchaoticdynamicsystems..............211.5.2.PeriodicorbitsofNewtonianthree-bodyproblem............221.5.3.Viscousflowpastasphere.........................241.5.4.Viscousflowpastacylinder........................251.5.5.Nonlinearwaterwaves...........................26References.........................................291.1.BackgroundPhysicalexperiment,numericalsimulationandanalytic(approximation)methodarethreemainstreamtoolstoinvestigatenonlinearproblems.Withoutdoubt,physicalexperimentisalwaysthebasicapproach.How-ever,physicalexperimentsareoftenexpensiveandtime-consuming.Be-sides,modelsforphysicalexperimentsareoftenmuchsmallerthantheoriginalones,butmostlyitisveryhardtosatisfyallsimilaritycriteri-ons.Bymeansofnumericalmethods,nonlinearequationsdefinedinrathercomplicateddomaincanbesolved.However,itisdifficulttogainnumeri-calsolutionsofnonlinearproblemswithsingularityandmultiplesolutionsordefinedinaninfinitydomain.Bymeansofanalytic(approximation)methods,onecaninvestigatenonlinearproblemswithsingularityandmul-tiplesolutionsinaninfinityinterval,butequationsshouldbedefinedinasimpleenoughdomain.So,physicalexperiments,numericalsimulationsandanalytic(approximation)methodshavetheirinherentadvantagesanddisadvantages.Therefore,eachofthemisimportantandusefulforustobetterunderstandnonlinearproblemsinscienceandengineering.Ingeneral,exact,closed-formsolutionsofnonlinearequationsarehardlyobtained.Perturbationtechniques[1–4]arewidelyusedtogainanalyticapproximationsofnonlinearequations.Usingperturbationmethods,manynonlinearequationsaresuccessfullysolved,andlotsofnonlinearphenom-enaareunderstoodbetter.Withoutdoubt,perturbationmethodsmakegreatcontributiontothedevelopmentofnonlinearscience.Perturbation
11October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod3methodsaremostlybasedonsmall(orlarge)physicalparameters,calledperturbationquantity.Usingsmall/largephysicalparameters,perturba-tionmethodstransferanonlinearequationintoaninfinitenumberofsub-problemsthataremostlylinear.Unfortunately,manynonlinearequationsdonotcontainsuchkindofperturbationquantitiesatall.Moreimpor-tantly,perturbationapproximationsoftenquicklybecomeinvalidwhentheso-calledperturbationquantitiesenlarge.Inaddition,perturbationtech-niquesaresostronglydependentuponphysicalsmallparametersthatwehavenearlynofreedomtochooseequationtypeandsolutionexpressionofhigh-orderapproximationequations,whichareoftencomplicatedandthusdifficulttosolve.Duetotheserestrictions,perturbationmethodsarevalidmostlyforweaklynonlinearproblemsingeneral.Ontheotherside,somenon-perturbationmethodswereproposedlongago.Theso-called“Lyapunov’sartificialsmall-parametermethod”[5]cantracebacktothefamousRussianmathematicianLyapunov(1857–1918),whofirstrewroteanonlinearequationN[u(r,t)]=L0[u(r,t)]+N0[u(r,t)]=f(r,t),(1.1)whererandtdenotethespatialandtemporalvariables,u(r,t)aunknownfunction,f(r,t)aknownfunction,L0andN0arelinearandnonlinearoperator,respectively,tosuchanewequationL0[u(r,t)]+qN0[u(r,t)]=f(r,t),(1.2)whereqhasnophysicalmeaning.Then,LyapunovregardedqasasmallparametertogainperturbationapproximationsX+∞23mu≈u0+u1q+u2q+u3q+···=u0+umq,(1.3)m=1andfinallygainedapproximationX+∞u≈u0+um,(1.4)m=1bysettingq=1,whereL0[u0(r,t)]=f(r,t),L0[u1(r,t)]=−N0[u0(r,t)],···(1.5)andsoon.ItshouldbeemphasizedthatonehasnofreedomtochoosethelinearoperatorL0inLyapunov’sartificialsmall-parametermethod:itisexactlythelinearpartofthewholeleft-handsideoftheoriginalequationN[u]=f,whereN=L0+N0.Thus,whenL0iscomplicatedor“singular”
12October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.14S.-J.Liao(forexample,itdoesnotcontainthehighestderivative),itisdifficult(orevenimpossible)tosolvethehigh-orderapproximationequation(1.5).Be-sides,theconvergenceoftheapproximationseries(1.4)isnotguaranteedingeneral.Evenso,Lyapunov’sexcellentworkisamilestoneofanalyticapproximationmethods,becauseitisindependentoftheexistenceofphys-icalsmallparameter,eventhoughitfirstregardsqasa“smallparameter”butfinallyenforcesittobe1thatishowevernot“small”strictlyfrommathematicalviewpoints.Theso-called“Adomiandecompositionmethod”(ADM)[6–8]wasde-velopedfromthe1970stothe1990sbyGeorgeAdomian,thechairoftheCenterforAppliedMathematicsattheUniversityofGeorgia,USA.Ado-mianrewrote(1.1)intheformN[u(r,t)]=LA[u(r,t)]+NA[u(r,t)]=f(r,t),(1.6)whereLAoftencorrespondstothehighestderivativeoftheequationunderconsideration,NA[u(r,t)]givestheleftpart,respectively.ApproximationsoftheADMarealsogivenby(1.4),too,whereLA[u0(r,t)]=f(r,t),LA[um(r,t)]=−Am−1(r,t),m≥1,(1.7)withtheso-calledAdomialpolynomial("#)1∂kX+∞nAk(r,t)=kNAun(r,t)q.(1.8)k!∂qn=0q=0SincethelinearoperatorLAissimplythehighestderivativeoftheconsid-eredequation,itisconvenienttosolvethehigh-orderapproximationequa-tions(1.7).ThisisanadvantageoftheADM,comparedto“Lyapunov’sartificialsmall-parametermethod”[5].However,theADMdoesnotpro-videsusfreedomtochoosethelinearoperatorLA,whichisrestrictedtoberelatedonlytothehighestderivative.Besides,like“Lyapunov’sartificialsmall-parametermethod”[5],theconvergenceoftheapproximationseries(1.4)givenbytheADMisstillnotguaranteed.Essentially,bothofthe“Lyapunov’sartificialsmallparametermethod”andthe“Adomiandecompositionmethod”transferanonlinearproblemintoaninfinitenumberoflinearsub-problems,withoutsmallphysicalpa-rameter.However,theyhavetwofundamentalrestrictions.First,onehasnofreedomandflexibilitytochoosethelinearoperatorsL0orLA,sinceL0isexactlythelinearpartofNandLAcorrespondstothehighestderivative,respectively.Second,thereisnowaytoguaranteetheconvergenceoftheapproximationseries(1.4).Thesecondonesismoreserious,sincedivergent
13October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod5approximationsaremostlyuseless.Thus,likeperturbationmethods,thetraditionalnon-perturbationmethods(suchasLyapunov’sartificialsmallparametermethodandtheADM)areoftenvalidforweaklynonlinearprob-lemsinmostcases.Intheory,itisveryvaluabletodevelopanewkindofanalyticapproxi-mationmethodwhichshouldhavethefollowingcharacteristics:(1)itisindependentofsmallphysicalparameter;(2)itprovidesusgreatfreedomandflexibilitytochoosetheequation-typeandsolutionexpressionofhigh-orderapproximationequations;(3)itprovidesusaconvenientwaytoguaranteetheconvergenceofapprox-imationseries.Oneofsuchkindofanalyticapproximationmethods,namelythe“homo-topyanalysismethod”(HAM)[9–17],wasdevelopedbyShijunLiaofrom1990sto2010s,togetherwithcontributionsofmanyotherresearchersintheoryandapplications.ThebasicideasoftheHAMwithitsbriefhistoryaredescribedbelow.1.2.AbriefhistoryoftheHAMThebasicideasof“Lyapunov’sartificialsmall-parametermethod”canbegeneralizedintheframeofthehomotopy,afundamentalconceptoftopol-ogy.ForanonlinearequationN[u(r,t)]=f(r,t),(1.9)Liao[9]proposetheso-called“homotopyanalysismethod”(HAM)byusingthehomotopy,abasicconceptintopology:(1−q)L[ϕ(r,t;q)−u0(r,t)]=c0qH(r,t){N[ϕ(r,t;q)]−f(r,t)},(1.10)whereLisanauxiliarylinearoperatorwiththepropertyL[0]=0,Nisthenonlinearoperatorrelatedtotheoriginalequation(1.9),q∈[0,1]istheembeddingparameterintopology(calledthehomotopyparameter),ϕ(r,t;q)isthesolutionof(1.10)forq∈[0,1],u0(r,t)isaninitialguess,c06=0istheso-called“convergence-controlparameter”,andH(r,t)isanauxiliaryfunctionthatisnon-zeroalmosteverywhere,respectively.Notethat,intheframeofthehomotopy,wehavegreatfreedomtochoosetheauxiliarylinearoperatorL,theinitialguessu0(r,t),theauxiliaryfunctionH(r,t),andthevalueoftheconvergence-controlparameterc0.
14October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.16S.-J.LiaoWhenq=0,duetothepropertyL[0]=0,wehavefrom(1.10)thesolutionϕ(r,t;0)=u0(r,t).(1.11)Whenq=1,sincec06=0andH(r,t)6=0almosteverywhere,Eq.(1.10)isequivalenttotheoriginalnonlinearequation(1.9)sothatwehaveϕ(r,t;1)=u(r,t),(1.12)whereu(r,t)isthesolutionoftheoriginalequation(1.9).Thus,astheho-motopyparameterqincreasesfrom0to1,thesolutionϕ(r,t;q)ofEq.(1.10)varies(ordeforms)continuouslyfromtheinitialguessu0(r,t)tothesolu-tionu(r,t)oftheoriginalequation(1.9).Forthissake,Eq.(1.10)iscalledthezeroth-orderdeformationequation.Here,itmustbeemphasizedonceagainthatwehavegreatfreedomandflexibilitytochoosetheauxiliarylinearoperatorL,theauxiliaryfunctionH(r,t),andespeciallythevalueoftheconvergencecontrolparameterc0inthezeroth-orderdeformationequation(1.10).Inotherwords,thesolutionϕ(r,t;q)ofthezeroth-orderdeformationequation(1.10)isalsodependentuponallaoftheauxiliarylinearoperatorL,theauxiliaryfunctionH(r,t)andtheconvergence-controlparameterc0asawhole,eventhoughtheyhavenophysicalmeanings.ThisisakeypointoftheHAM,whichwewilldiscussindetailslater.AssumethatL,H(r,t)andc0areproperlychosensothatthesolutionϕ(r,t;q)ofthezeroth-orderdeformationequation(1.10)alwaysexistsforq∈(0,1)andbesidesitisanalyticatq=0,andthattheMaclaurinseriesofϕ(r,t;q)withrespecttoq,i.e.X+∞ϕ(r,t;q)=u(r,t)+u(r,t)qm(1.13)0mm=1convergesatq=1.Then,dueto(1.12),wehavetheapproximationseriesX+∞u(r,t)=u0(r,t)+um(r,t).(1.14)m=1Substitutingtheseries(1.13)intothezeroth-orderdeformationequation(1.10)andequatingthelike-powerofq,wehavethehigh-orderapproxima-tionequationsforum(r,t),calledthemth-orderdeformationequationL[um(r,t)−χmum−1(r,t)]=c0H(r,t)Rm−1(r,t),(1.15)aMorestrictly,ϕ(r,t;q)shouldbereplacedbyϕ(r,t;q,L,H(r,t),c0).Onlyforthesakeofsimplicity,weusehereϕ(r,t;q),butshouldalwayskeepthispointinmind.
15October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod7where("#!)1∂kX+∞nRk(r,t)=N(r,t)q−f(r,t),(1.16)k!∂qkn=0q=0withthedefinition0,whenk≤1,χk=(1.17)1,whenk≥2.Forvarioustypesofnonlinearequations,itiseasyandstraightforwardtousethetheoremsprovedinChapter4ofLiao’sbook[11]tocalculatethetermRk(r,t)ofthehigh-orderdeformationequation(1.15).ItshouldbeemphasizedthattheHAMprovidesusgreatfreedomandflexibilitytochoosetheauxiliarylinearoperatorLandtheinitialguessu0.Thus,differentfromallotheranalyticmethods,theHAMprovidesusgreatfreedomandflexibilitytochoosetheequationtypeandsolutionexpressionofthehigh-orderdeformationequation(1.15)sothatitssolutioncanbeof-tengainedwithoutgreatdifficulty.Noticethat“theessenceofmathematicsliesentirelyinitsfreedom”,aspointedoutbyGeorgCantor(1845–1918).Moreimportantly,thehigh-orderdeformationequation(1.15)containstheconvergence-controlparameterc0,andtheHAMprovidesgreatfreedomtochoosethevalueofc0.Mathematically,ithasbeenprovedthattheconvergence-controlparameterc0canadjustandcontroltheconvergenceregionandratiooftheapproximationseries(1.14).Fordetails,pleaserefertoLiao[10,12,13]andespecially§5.2to§5.4ofhisbook[11].So,unlikeallotheranalyticapproximationmethods,theconvergence-controlparameterc0oftheHAMprovidesusaconvenientwaytoguaranteetheconvergenceoftheapproximationseries(1.14).Infact,itistheconvergence-controlparameterc0thatdifferstheHAMfromallotheranalyticmethods.Atthemth-orderofapproximation,theoptimalvalueoftheconvergence-controlparameterc0canbedeterminedbytheminimumofresidualsquareoftheoriginalgoverningequation,i.e.dEm=0,(1.18)dc0whereZ("m#)2XEm=Nun(r,t)−f(r,t)dΩ.(1.19)Ωn=0Besides,ithasbeenprovedbyLiao[16]thatahomotopyseriessolution(1.14)mustbeoneofsolutionsofconsideredequation,aslongasitis
16October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.18S.-J.Liaoconvergent.Inotherwords,foranarbitraryconvergence-controlparameterc0∈Rc,whereRc=c0:limEm(c0)→0(1.20)m→+∞isaninterval,thesolutionseries(1.14)isconvergenttothetruesolutionoftheoriginalequation(1.9).Fordetails,pleaserefertoLiao[16]andChapter3ofhisbook[11].Insummary,theHAMhasthefollowingadvantages:(a)itisindependentofanysmall/largephysicalparameters;(b)itprovidesusgreatfreedomandlargeflexibilitytochooseequationtypeandsolutionexpressionoflinearhigh-orderapproximationequations;(c)itprovidesusaconvenientwaytoguaranteetheconvergenceofapprox-imationseries.Inthisway,nearlyallrestrictionsandlimitationsofthetraditionalnon-perturbationmethods(suchasLyapunov’sartificialsmallparametermethod[5],theAdomiandecompositionmethod[6–8],theδ-expansionmethod[18]andsoon)canbeovercomebymeansoftheHAM.Besides,ithasbeengenerallyproved[10,12,13]thattheLyapunov’sartificialsmallparametermethod[5],theAdomiandecompositionmethod[6–8]andtheδ-expansionmethod[18]areonlyspecialcasesoftheHAMforsomespeciallychosenauxiliarylinearoperatorLandconvergence-controlparameterc0.Especially,theso-called“homotopyperturbationmethod”(HPM)[19]proposedbyJihuanHein1998(sixyearslaterafterLiao[9]proposedtheearlyHAMin1992)wasonlyaspecialcaseoftheHAMwhenc0=−1,andthushas“nothingnewexceptitsname”[20].SomeresultsgivenbytheHPMaredivergenteveninthewholeintervalexceptthegiveninitial/boundaryconditions,andthus“itisveryimportanttoinvestigatetheconvergenceofapproximationseries,otherwiseonemightgetuselessresults”,aspointedoutbyLiangandJeffrey[21].Fordetails,see§6.2ofLiao’sbook[11].Thus,theHAMismoregeneralintheoryandwidelyvalidinpracticeformoreofnonlinearproblemsthanotheranalyticapproximationtechniques.Incalculus,thefamousEulertransformisoftenusedtoacceleratecon-vergenceofaseriesortomakeadivergentseriesconvergent.Itisinter-estingthatonecanderivetheEulertransformintheframeoftheHAM,andgiveasimilarbutmoregeneraltransform(calledthegeneralizedEulertransform),asshowninChapter5ofLiao’sbook[11].Thisprovidesusa
17October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod9theoreticalcornerstoneforthevalidityandgeneralityoftheHAM.Theintroductionoftheso-called“convergence-controlparameter”c0inthezeroth-orderdeformationequation(1.10)isamilestonefortheHAM.Fromphysicalviewpoint,the“convergence-controlparameter”c0hasnophysicalmeaningssothatconvergentseriesofsolutiongivenbytheHAMmustbeindependentofc0.Thisisindeedtrue:thereexistssucharegionRcthat,forarbitraryc0∈Rc,theHAMseriesconvergestothetruesolutionoftheoriginalequation(1.9),asillustratedbyLiao[10,11].However,ifc06∈Rc,thesolutionseriesdiverges!So,fromamathematicalviewpoint,the“convergence-controlparameter”isakeypointoftheHAM,whichprovidesusaconvenientwaytoguaranteetheconvergenceofthesolutionseries.Infact,itistheso-called“convergence-controlparameter”thatdifferstheHAMfromallotheranalyticapproximationmethods.Theintroductionofthebasicconcepthomotopyintopologyisalsoamilestoneoftheanalyticapproximationmethodsfornonlinearproblems.ItisthehomotopythatprovidesusgreatfreedomandlargeflexibilitytochoosetheauxiliarylinearoperatorLandinitialguessu0inthezeroth-orderdeformationequation(1.10),whichdeterminetheequationtypeandsolutionexpressionofthehigh-orderdeformationequations(1.15).Besides,itisthehomotopythatprovidesusthefreedomtointroducetheso-called“convergence-controlparameter”c0in(1.10),whichbecomesnowacor-nerstoneoftheHAM.Notethatitisimpossibletointroducesuchkindof“convergence-controlparameter”intheframeofperturbationtechniquesandthetraditionalnon-perturbationmethods(suchasLyapunov’sartificialsmallparameter,Adomiandecompositionmethodandsoon).ThefreedomonthechoiceoftheauxiliarylinearoperatorLissolargethatthesecond-ordernonlinearGelfandequationcanbesolvedconveniently(withgoodagreementwithnumericalresults)intheframeoftheHAMevenbymeansofaforth-orderauxiliarylinearoperator(fortwodimen-sionalGelfandequation)orasixth-orderauxiliarylinearoperator(forthreedimensionalGelfandequation),respectively,asillustratedbyLiao[14].Al-thoughitistruethattheauxiliarylinearoperator(withthesamehighestorderofderivativeasthatofconsideredproblem)canbechosenstraight-forwardlyinmostcases,suchkindoffreedomoftheHAMshouldbetakenintoaccountsufficientlybytheHAMcommunitywhennecessary,especiallyforsomevaluablebutchallengingproblems(someofthemaresuggestedbelowin§1.5).Inaddition,bymeansoftheabove-mentionedfreedomoftheHAM,theconvergenceofapproximationsolutioncanbegreatlyacceleratedinthe
18October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.110S.-J.LiaoframeoftheHAMbymeansoftheiteration,theso-calledhomotopy-Pad´etechniqueandsoon.Fordetails,pleasereferto§2.3.5to§2.3.7ofLiao’sbook[11].Indeed,“theessenceofmathematicsliesentirelyinitsfreedom”,aspointedoutbyGeorgCantor(1845–1918).SuchkindofgreatfreedomoftheHAMshouldprovideusgreatpossi-bilitytosolvesomeopenquestions.Oneofthemisdescribedbelow.Thesolutionofthehigh-orderdeformationequation(1.15)canbeexpressedintheform−1um(r,t)=−χmum−1(r,t)+L[c0H(r,t)Rm−1(r,t)],(1.21)whereL−1istheinverseoperatorofL.ForafewauxiliarylinearoperatorL,itsinverseoperatorissimple.However,inmostcases,itisnotstraight-forwardtosolvetheabovelineardifferentialequation.Canwedirectlychoose(ordefine)theinverseauxiliarylinearoperatorL−1soastosolve(1.15)conveniently?ThisispossibleintheframeoftheHAM,sinceintheorytheHAMprovidesusgreatfreedomandlargeflexibilitytochoosetheauxiliarylinearoperatorL.Ifsuccessful,itwouldberatherefficientandconvenienttosolvethehigh-orderdeformationequation(1.15).ThisisaninterestingbutopenquestionfortheHAMcommunity,whichdeservestobestudiedindetails.Notethatsomeinterestingproblemsaresuggestedin§1.5.1.3.SomeadvancesoftheHAMSince1992whenLiao[9]proposedtheearlyHAM,theHAMhasbeendevelopinggreatlyintheoryandapplications,duetothecontributionsofmanyresearchersindozensofcountries.Unfortunately,itisimpossibletodescribealloftheseadvancesindetailsinthisbriefreview,andeveninthisbook.Infact,theHAMhasbeensuccessfullyappliedtonumerous,varioustypesofnonlinearproblemsinscience,engineeringandfinance.So,wehadtofocusonarathersmallportoftheseadvanceshere.1.3.1.Generalizedzeroth-orderdeformationequationThestartingpointoftheuseoftheHAMistoconstructtheso-calledzeroth-orderdeformationequation,whichbuildsaconnection(i.e.acontinuousmapping/deformation)betweenagivennonlinearproblemandarelativelymuchsimplerlinearones.So,thezeroth-orderdeformationequationisabaseoftheHAM.
19October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod11Givenanonlinearequation,wehavegreatfreedomandlargeflexibilityintheframeoftheHAMtoconstructtheso-calledzeroth-orderdefor-mationequationusingtheconcepthomotopyintopology.Especially,theconvergence-controlparameterc0playsanimportantroleintheframeoftheHAM.So,itisnaturaltoenhancetheabilityoftheso-called“convergencecontrol”bymeansofintroducingmoresuchkindofauxiliaryparameters.Duetotheabove-mentionedfreedomandflexibilityoftheHAM,therearenumerousapproachestodoso.Forexample,wecanconstructsuchakindofzeroth-orderdeformationequationwithK+1convergence-controlparameters:(1−q)L[ϕ(r,t;q)−u0(r,t)]!XK=cqn+1H(r,t){N[ϕ(r,t;q)]−f(r,t)},(1.22)nn=0whereϕ(r,t;q)isthesolution,NisanonlinearoperatorrelatedtoanoriginalproblemN[u(r,t)]=f(r,t),q∈[0,1]isthehomotopyparameter,u0isaninitialguess,Lisanauxiliarylinearoperator,H(r,t)isanauxiliaryfunctionwhichisnonzeroalmosteverywhere,andc={c0,c1,···,cK}isavectorof(K+1)non-zeroconvergence-controlparameters,respectively.Notethat,whenK=0,itgivesexactlythezeroth-orderdeformationequation(1.10).Thecorrespondinghigh-orderdeformationequationreadsmin{Xm−1,K}L[um(r,t)−χmum−1(r,t)]=H(r,t)cnRm−1−n(r,t),(1.23)n=0whereRn(r,t)andχnaredefinedbythesameformulas(1.16)and(1.17),respectively.WhenK=0,theabovehigh-orderdeformationequation(1.23)isexactlythesameas(1.15).Atthemth-orderofapproximation,theoptimalconvergence-controlparametersaredeterminedbytheminimumoftheresidualsquareoftheoriginalequation,i.e.dEm=0,0≤n≤min{m−1,K},(1.24)dcnwhereEmisdefinedby(1.19).Fordetails,pleaserefertoChapter4ofLiao’sbook[11].WhenK→+∞,itisexactlytheso-called“optimalhomotopyasymptoticmethod”[22].So,the“optimalhomotopyasymptotic
20October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.112S.-J.Liaomethod”[22]isalsoaspecialcaseoftheHAM,asshownin§3.2.2and§6.3ofLiao’sbook[11].Intheory,themoretheconvergence-controlparameters,thelargertheabilitytocontroltheconvergenceoftheHAMseries.However,itisfound[16]thatmuchmoreCPUtimesisneededinpracticewhenmoreconvergence-controlparametersareused.Inmostcases,oneoptimalconvergence-controlparameterisgoodenoughtogainconvergentresultsbymeansoftheHAM.Consideringthecomputationalefficiency,oneuptothreeconvergence-controlparametersaregenerallysuggestedintheframeoftheHAM.Fordetails,pleasereferto§2.3.3,§2.3.4,§4.6.1andChapter3ofLiao’sbook[16].Itshouldbeemphasizedonceagainthat,intheframeofthehomotopyintopology,wehaverathergreatfreedomandlargeflexibilitytoconstructtheso-calledzeroth-orderdeformationequation.Intheory,givenanonlinearequationN[u(r,t)]=f(r,t),wecanalwaysproperlychooseaninitialguessu0(r,t)andanauxiliarylinearoperatorLtoconstructsuchazeroth-orderdeformationequationinarathergeneralformA[u0(r,t),L,ϕ(r,t;q),c;q]=0(1.25)thatitholdsϕ(r,t;0)=u0(r,t),whenq=0,(1.26)andϕ(r,t;1)=u(r,t),whenq=1,(1.27)i.e.,whenq=1thezeroth-orderdeformationequation(1.25)isequivalenttotheoriginalnonlinearequationN[u(r,t)]=f(r,t).UsingthetheoremsgiveninChapter4ofLiao’sbook[11],itiseasytogainthecorrespondinghigh-orderdeformationequations.Here,c={c0,c1,···,cK}isavectorofconvergence-controlparameters,whoseoptimalvaluesarede-terminedbytheminimumofresidualsquareoftheoriginalequation.Notethat(1.25)israthergeneral:thezeroth-orderdeformationequations(1.10)and(1.22),andevenEq.(1.2)forLyapunov’sartificialsmallparametermethod,areonlyspecialcasesof(1.25).Somecommonlyusedzeroth-orderdeformationequationsaredescribedin§4.3ofLiao’sbook[11]asspecialcasesofthegeneralizedzeroth-orderdeformationequation(1.25).Intheory,thereareaninfinitenumberofdifferentwaystoconstructazeroth-orderdeformationequation(1.25).Therefore,intheframeofthe
21October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod13HAM,weindeedhavehugefreedomandflexibility.Suchkindoffreedomandflexibilitycomesfromthehomotopy,abasicconceptintopology.Intheory,thiskindoffreedomandflexibilityprovidesusgreatabilitytosolvesomeinterestingbutchallengingnonlinearproblems(someofthemaresuggestedbelowin§1.5ofthischapter),ifwecanclearlyknowhowtousetheminaproperwaywithstimulatedimagination!Inpractice,itissuggestedtofirstlyusethezeroth-orderdeformationequation(1.10),sinceitworksformostofnonlinearproblems,asillustratedbyLiao[10,11].Ifunsuccessful,onecanfurtherattemptalittlemorecomplicatedzeroth-orderdeformationequations,suchas(1.22).Finally,weemphasizeonceagainthat,intheory,onehashugefreedomtoconstructazeroth-orderdeformationequation(1.25)satisfyingbothof(1.26)and(1.27),aslongasoneclearlyknowshowtousesuchkindoffreedom.1.3.2.SpectralHAMandcomplicatedauxiliaryoperatorAlthoughtheHAMprovidesusgreatfreedomtochoosetheauxiliarylinearoperatorL,itmightbedifficulttosolvethelinearhigh-orderdeformationequation(1.15)or(1.23)exactly,ifLiscomplicated.Thisismainlybecausemostoflineardifferentialequationshavenoclosed-formsolutions,i.e.theirsolutionsaremostlyexpressedbyaninfiniteseries.So,inordertoexactlysolvehigh-orderdeformationequationsintheframeoftheHAM,weoftenshouldchooseareasonablebutsimpleenoughauxiliarylinearoperatorL.This,however,restrictstheapplicationsoftheHAM.Thisisthemainreasonwhyonlyafewsimpleauxiliarylinearoperators,suchasLu=u0,Lu=xu0+u,Lu=u0+u,Lu=u00+uandsoon,havebeenmostlyusedintheframeoftheHAM,wheretheprimedenotesthedifferentiationwithrespecttox.Theseauxiliarylinearoperatorscorrespondtosomefundamentalfunctionssuchaspolynomial,exponential,trigonometricfunctionsandtheircombination.Therearemanyspecialfunctionsgovernedbylineardifferentialequa-tions.Althoughmanysolutionscanexpressedbythesespecialfunctions,theyarehardlyusedintheframeoftheHAMuptonow,becausethecorrespondinghigh-orderdeformationequationsoftenbecomemoreandmoredifficulttosolve.Thisisapity,sinceintheorytheHAMindeedprovidesusfreedomtousespecialfunctionstoexpresssolutionsofmanynonlineardifferentialequations.Currently,VanGorder[23]madeaninter-
22October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.114S.-J.Liaoestingattemptinthisdirection.IntheframeoftheHAM,VanGorder[23]expressedanalyticapproximationsoftheFitzhugh–Nagumoequationbymeansoferrorfunction,Gaussianfunctionandsoon.ThekeyisthatVanGorder[23]chosesuchanauxiliarylinearoperator2z2−1Lu=u00+u0,zwheretheprimedenotesthedifferentiationwithrespecttoz,andespeciallysuchaproperauxiliaryfunctionH(z)=z|z|,thatthecorrespondinghigh-orderdeformationequationscanbesolvedeasily.Fordetails,pleaserefertoVanGorder[23],VajraveluandVanGorder[24]and§4.6ofthisbook.ThisexampleillustratesonceagainthattheHAMindeedprovidesusgreatfreedom,i.e.lotsofpossibilities.Thekeyishowtousesuchkindoffreedom!Generallyspeaking,solutionofacomplicatedlinearODE/PDEshouldbeexpressedinaserieswithaninfinitenumberofterms.Mathematically,suchaseriesleadstothelargerandlargerdifficultytogainhigher-orderanalyticapproximationsofanonlinearproblem.Fortunately,fromphysicalviewpoint,itisoftenaccurateenoughtohaveanalyticapproximationswithmanyenoughterms.Currently,usingtheSchmidt-Gramprocess,Zhao,LinandLiao[25]suggestedaneffectivetruncationtechniqueintheframeoftheHAM,whichcanbeusedtogreatlysimplifytheright-handsideofthehigh-orderdeformationequations,suchas(1.15)and(1.23),priortosolvingthem.Inthisway,muchCPUtimecanbesaved,evenwithoutlossofaccuracy.In2010,Motsaetal.[26,27]suggestedtheso-called“spectralhomotopyanalysismethod”(SHAM)usingtheChebyshevpseudospectralmethodtosolvethelinearhigh-orderdeformationequationsandchoosingtheauxiliarylinearoperatorLintermsoftheChebyshevspectralcollocationdifferentia-tionmatrix[28].Intheory,anyacontinuousfunctioninaboundedintervalcanbebestapproximatedbyChebyshevpolynomial.So,theSHAMpro-videsuslargerfreedomtochoosetheauxiliarylinearoperatorLandinitialguessintheframeoftheHAM.ItisvaluabletoexpandtheSHAMfornonlinearpartialdifferentialequations.Besides,itiseasytoemploytheoptimalconvergence-controlparameterintheframeoftheSHAM.Thus,theSHAMhasgreatpotentialtosolvemorecomplicatednonlinearprob-lems,althoughfurthermodificationsintheoryandmoreapplicationsareneeded.ForthedetailsabouttheSHAM,pleasereferto[26,27]andChap-ter3ofthisbook.
23October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod15Chebyshevpolynomialisjustoneofspecialfunctions.TherearemanyotherspecialfunctionssuchasHermitepolynomial,Legendrepolynomial,Airyfunction,Besselfunction,Riemannzetafunction,hypergeometricfunctions,errorfunction,Gaussianfunctionandsoon.SincetheHAMprovidesusextremelylargefreedomtochooseauxiliarylinearoperatorLandinitialguess,itshouldbepossibletodevelopa“generalizedspectralHAM”whichcanuseproperspecialfunctionsforsomenonlinearproblems.Especially,combinedtheSHAM[26,27]withtheabove-mentionedtrunca-tiontechniquesuggestedbyZhao,LinandLiao[25],itwouldbepossibletouse,whennecessary,morecomplicatedauxiliarylinearoperatorsintheframeoftheHAMsothatsomedifficultnonlinearproblemscanbesolved.1.3.3.PredictorHAMandmultiplesolutionsManynonlinearboundaryvalueproblemshavemultiplesolutions.Ingen-eral,itisdifficulttogainthesedualsolutionsbymeansofnumericaltech-niques,mainlybecausedualsolutionsareoftenstronglydependentuponinitialconditionsbutwedonotknowhowtochoosethemexactly.Compar-ativelyspeaking,itisalittlemoreconvenienttouseanalyticapproximationmethodstosearchformultiplesolutionsofnonlinearproblems,sinceana-lyticmethodsadmitunknownvariablesininitialguess.Forexample,letusconsiderasecond-ordernonlineardifferentialequa-tionofatwo-pointboundaryvalueproblem:N[u(x)]=0,u(0)=a,u(1)=b,(1.28)whereNisa2nd-ordernonlineardifferentialoperator,aandbareknownconstants,respectively.Assumethatthereexistmultiplesolutionsu(x).Thesemultiplesolutionsmusthavesomethingdifferent.Withoutlossofgenerality,assumethattheyhavedifferentfirst-orderderivativeu0(0)=σ,whereσisunknown.Obviously,differentinitialguessu0(x)mightleadtomultiplesolutions.Fortunately,theHAMprovidesusgreatfreedomtochooseinitialguessu0(x).Asmentionedbefore,suchkindoffreedomisonecornerstoneoftheHAM.So,intheframeoftheHAM,itisconvenientforustochoosesuchaninitialguessu0(x)thatitsatisfiesnotonlythetwoboundaryconditionsu(0)=a,u(1)=bbutalsotheadditionalconditionu0(0)=σ.Inthisway,theinitialguessu0(x)containsanunknownparameterσ,calledbyLiao(seeChapter8of[11])themultiple-solution-controlparameter.Then,theanalyticapproximationsgainedbytheHAMcontainatleasttwoun-
24October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.116S.-J.Liaoknownauxiliaryparameters:theconvergence-controlparameterc0andthemultiple-solution-controlparameterσ.AssuggestedbyLiao(seeChapter8of[11]),theoptimalvaluesofc0andσcanbedeterminedbytheminimumoftheresidualsquareofgoverningequations.Inthisway,multiplesolu-tionsofsomenonlineardifferentialequationscanbegained,asillustratedbyLiao(seeChapter8of[11]).IntheframeoftheHAM,AbbasbandyandShivanian[29,30]developedadifferentialbutratherinterestingapproachtogaindualsolutions,namelythePredictorHAM(PHAM).Forsimplicity,letususethesameequation(1.28)asanexampletodescribeitsbasicideas.Firstofall,anadditionalconditionsuchasu0(0)=σisintroducedwiththeunknownparameterσ.Then,intheframeoftheHAM,onesolvesthenonlineardifferentialequationN[u(x)]=0,butwiththetwoboundaryconditionsu0(0)=σandu(1)=b.Then,u(0),theHAMapproximationatx=0,containsatleasttwounknownparameters:oneistheso-calledconvergence-controlparameterc,theotherisσ=u0(0),calledthemultiple-solution-control0parameterbyLiao(seeChapter8of[11])intheabove-mentionedapproach.Substitutingtheexpressionofu(0)intotheboundaryconditionu(0)=agivesanonlinearalgebraicequationaboutc0andσ.Fromthephysicalviewpoint,σ=u0(0)hasphysicalmeanings,buttheconvergence-controlparameterc0doesnot.Iftheorderofapproximationishighenough,onecangainconvergent,accurateenoughmultiplevaluesofσforproperlychosenvaluesofc0inafiniteinterval,asillustratedin[29,30].Inthisway,onecanfindmultiplesolutionsofagivennonlinearproblem.Fordetails,pleaserefertoChapter2ofthisbook.IntheframeoftheHAM,somenewbranchesofsolutionsforviscousboundary-layerflowswerefound[31,32],andthemultipleequilibrium-statesofresonantwavesindeepwater[33]andinfinitewaterdepth[34]werediscoveredforthefirsttime,tothebestofauthor’sknowledge.Alloftheseillustratethepotential,noveltyandvalidityoftheHAMtogivesomethingnewanddifferent.ThisisasuperiorityoftheHAMtonumericalmethodsandsomeotheranalyticapproximationtechniques.Certainly,itisvaluabletoapplytheHAMtodiscoversomenewsolutionsofothernonlinearproblems!1.3.4.ConvergenceconditionandHAM-basedsoftwareTheoreticallyspeaking,theHAMindeedprovidesusgreatfreedomtochooseinitialguess,auxiliarylinearoperator,convergence-controlparame-
25October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod17ter,equation-typeandsolution-expressionofhigh-orderdeformationequa-tion,andsoon.However,itisstillnotveryclearhowtousethesefreedomintheframeoftheHAM,mainlybecauselittlemathematicaltheoremshavebeenprovedinanabstractway.Somestudiesonthestabilityofauxiliarylinearoperatorandconvergence-controlparameteroftheHAMaredescribedinChapter4ofthisbook.SomecurrentworksaboutconvergenceconditionoftheHAMseriesaredescribedinChapter5.ItshouldbeespeciallyemphasizedthatParkandKim[35,36]success-fullyappliedtheHAMtosolveafewclassicproblemsinfinance,andgaveconvergenceconditionsfortheiranalyticapproximations.Itisratherinter-estingthattheyevengaveanerrorestimationfortheiranalyticapproxima-tionsin[36].Currently,ParkandKimusedtheHAMtosolveanabstractlinearproblemwithrespecttoboundedlinearoperatorsfromaBanachspacetoaBanachspace,andrigorouslyprovedthatthehomotopysolutionexistsinthesensethataseriesoftheproblemconvergesinaBanachnormsenseifthelinearoperatorsatisfiessomemildconditions.Theirfantasticworksareveryimportant,andmightpioneeranewresearchdirectionandstyle(i.e.abstractproof)intheframeoftheHAM.SuchkindofabstractmathematicaltheoremsintheframeoftheHAMaremorevaluableanduseful,ifnonlineargoverningequationsandespeciallytheinfluenceoftheconvergence-controlparameterontheconvergencecouldbeconsidered.Ontheotherside,theHAMhasbeensuccessfullyappliedtonumerousnonlinearproblemsinvariousfieldsofscienceandengineering.Theseap-plicationsshowthegeneralvalidityandnoveltyoftheHAM.Unfortunately,itisimpossibletomentionallofthemhereindetails.Asexamplesamongthesenumerousapplications,aHAM-basedapproachaboutboundary-layerflowsofnanofluidisgiveninChapter6ofthisbook.Inaddition,anappli-cationoftheHAMfortime-fractionalboundary-valueproblemisillustratedinChapter7.TosimplifysomeapplicationsoftheHAM,twoHAM-basedsoftwareweredeveloped.TheHAM-basedMaplepackageNOPH(version1.0.2)forperiodicoscillationsandlimitcyclesofnonlineardynamicsystemsisdescribedinChapter8ofthisbookwithvariousapplications.Itisfreeavailableonlineathttp://numericaltank.sjtu.edu.cn/NOPH.htmwithasimpleuser’sguide.Besides,theHAM-basedMathematicapackageBVPh(version2.0)forcouplednonlinearordinarydifferentialequations
26October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.118S.-J.LiaowithboundaryconditionsatmultiplepointsaregiveninChapter9ofthisbook.Itisfreeavailableonlineathttp://numericaltank.sjtu.edu.cn/BVPh.htmwithasimpleuser’sguideandsomeexamplesofapplication.BothofthesetwoHAM-basedsoftwareareeasy-to-useanduser-friendly.TheygreatlysimplifysomeapplicationsoftheHAM,andareespeciallyhelpfulforthebeginnersoftheHAM.1.4.RelationshipstoothermethodsInpuremathematics,thehomotopyisafundamentalconceptintopologyanddifferentialgeometry.TheconceptofhomotopycanbetracedbacktoJulesHenriPoincar´e(1854–1912),aFrenchmathematician.Ahomotopydescribesakindofcontinuousvariation(ordeformation)inmathematics.Forexample,acirclecanbecontinuouslydeformedintoasquareoranellipse,theshapeofacoffeecupcandeformcontinuouslyintotheshapeofadoughnutbutcannotbedistortedcontinuouslyintotheshapeofafootball.Essentially,ahomotopydefinesaconnectionbetweendifferentthingsinmathematics,whichcontainsamecharacteristicsinsomeaspects.Inpuremathematics,thehomotopyiswidelyusedtoinvestigateexistenceanduniquenessofsolutionsofsomeequations,andsoon.Inappliedmathematics,theconceptofhomotopywasusedlongagotodevelopsomenumericaltechniquesfornonlinearalgebraicequations.Theso-called“differentialarclengthhomotopycontinuationmethod”werepro-posedin1970sbyKeller[37,38].However,theglobalhomotopymethodscanbetracedasfarbackastheworkofLahaye[39]in1934.Tosolveanon-linearalgebraicequationf(x)=0bymeansofthehomotopycontinuationmethod,onefirstconstructssuchahomotopyH(x,q)=qf(x)+(1−q)g(x),(1.29)whereq∈[0,1]isthehomotopyparameter,g(x)isafunctionforwhichazeroisknownorreadilyobtained.AsdiscussedbyWayburnandSeader[40],thechoiceofg(x)isarbitrary,butthetwomostwidelyusedfunctionsaretheNewtonhomotopyH(x,q)=qf(x)+(1−q)[f(x)−f(x0)],andthefixed-pointhomotopyH(x,q)=qf(x)+(1−q)(x−x0),
27October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod19wherex0isanarbitrarystartingpoint.Thelocusofsolutionsdefinesthehomotopypaththatistrackedwithsomecontinuationmethod.Conse-quently,homotopycontinuationmethodsconsistnotonlyofthehomotopyequationitself,butalsothehomotopypathtrackingmethod,i.e.ofsomecontinuationstrategy.Homotopycontinuationmethodsareusuallybasedupondifferentiatingthehomotopyequation(1.29)withrespecttothearclengths,whichgivestheequation∂Hdx∂Hdq+=0.(1.30)∂xds∂qdsTakingintoaccountthearc-lengthrelation22dxdq+=1dsdsandtheinitialconditionH(x0,0)=0,weobtainaninitialvalueprob-lem.Then,pathtrackingbasedontheinitialvalueproblemisnumericallycarriedoutwithapredictor-correctoralgorithmtogainasolutionoftheoriginalequationf(x)=0.Someeleganttheoremsofconvergenceareprovedintheframeofthehomotopycontinuationmethod.Fordetailsofthehomotopycontinuationmethod,pleasereferto[41–49].Unlikethehomotopycontinuationmethodthatisaglobalconvergentnumericalmethodmainlyfornonlinearalgebraicequations,theHAMisakindofanalyticapproximationmethodmainlyfornonlineardifferentialequations.So,theHAMisessentiallydifferentfromthehomotopycon-tinuationmethod,althoughbothofthemarebasedonthehomotopy,thebasicconceptofthetopology.NotethattheHAMusesmuchmorecom-plicatedhomotopyequation(1.10)or(1.22)than(1.29)forthehomotopycontinuationmethod.Furthermore,theHAMprovideslargerfreedomtochoosetheauxiliarylinearoperatorL.Mostimportantly,theso-calledconvergence-controlparameterc0isintroducedforthefirsttime,tothebestofourknowledge,inthehomotopyequation(1.10)or(1.22)sothattheHAMprovidesusaconvenientwaytoguaranteetheconvergenceofseriesseries.Notethatthehomotopyequation(1.29)ofthehomotopycontinuationmethoddoesnotcontainsuchkindofconvergence-controlpa-rameteratall.So,theconvergence-controlparameterc0isindeedanovel.Infact,itistheconvergence-controlparameterc0whichdifferstheHAMfromallotheranalyticapproximationmethods.Inaddition,theHAMlogicallycontainsmanyotheranalyticapproxima-tionmethodsandthusisrathergeneral.Forexample,ithasbeengenerally
28October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.120S.-J.Liaoproved[10–13]thattheLyapunov’sartificialsmallparametermethod[5],theAdomiandecompositionmethod[6–8],theδ-expansionmethod[18]areonlyspecialcasesoftheHAMforsomespeciallychosenauxiliarylinearoperatorLandconvergence-controlparameterc0.Furthermore,theso-called“optimalhomotopyasymptoticmethod”[22]developedin2008isalsoaspecialcaseofthehomotopyequation(1.22)whenK→+∞,too,aspointedoutbyLiao(see§6.3ofLiao’sbook[11]).Especially,theso-called“homotopyperturbationmethod”(HPM)[19]proposedbyJihuanHein1998(sixyearslaterafterLiao[9]proposedtheearlyHAMin1992)wasonlyaspecialcaseoftheHAMwhenc0=−1,asprovedin[20],andthusithas“nothingnewexceptitsname”[20].SomeresultsgivenbytheHPMaredivergenteveninthewholeintervalexceptthegiveninitial/boundaryconditions,andthus“itisveryimportanttoinvestigatetheconvergenceofapproximationseries,otherwiseonemightgetuselessresults”,aspointedoutbyLiangandJeffrey[21].Formoredetails,see§6.2ofLiao’sbook[11].Inaddition,eventhefamousEulertransformincalculuscanbederivedintheframeoftheHAM(seeChapter5ofLiao’sbook[11]).ThisprovidesusatheoreticalcornerstoneforthevalidityandgeneralityoftheHAM.Insummary,basedontheconceptofhomotopytopology,theHAMisanovelanalyticapproximationmethodforhighlynonlinearproblems,withgreatfreedomandflexibilitytochooseequation-typeandsolutionex-pressionofhigh-orderapproximationequationsandalsowithaconvenientwaytoguaranteetheconvergence,sothatitmightovercomerestrictionsofperturbationtechniquesandothernon-perturbationmethods.1.5.Chanceandchallenge:somesuggestedproblemsAnytrulynewmethodsshouldgivesomethingnoveland/ordifferent,orsolvesomedifficultproblemsthatcannotbesolvedwithsatisfactionbyothermethods.Unlikeotheranalyticapproximationmethods,theHAMprovidesusgreatfreedomandflexibilitytochooseequation-typeandsolutionexpres-sionofhigh-orderapproximationequations,andespeciallyasimplewaytoguaranteetheconvergenceofsolutionseries.Thus,theHAMprovidesusalargepossibilityandchancetogivesomethingnovelordifferent,andtoattacksomedifficultnonlinearproblems.Forexample,somenewso-lutions[31,32]ofboundary-layerflowshavebeenfoundbymeansoftheHAM,whichhadbeenneglectedevenbynumericaltechniquesandhad
29October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod21beenneverreported.Someanalyticapproximationsfortheoptimalex-erciseboundaryofAmericanputoptionweregiven,whicharevalidfromacoupleofyears(see[50,51])uptoeven20years(seeChapter13ofLiao’sbook[11])priortoexpiry,andthusmuchbetterthantheasymp-totic/perturbationapproximationsthatareoftenvalidonlyforacoupleofdaysorweeks.Besides,theHAMhasbeensuccessfullyemployedtosolvesomecomplicatednonlinearPDEs:themultipleequilibrium-statesofreso-nantwavesindeepwater[33]andinfinitewaterdepth[34]werediscoveredbymeansoftheHAMforthefirsttime,tothebestofourknowledge,whichgreatlydeepenandenrichourunderstandingsaboutresonantwaves.Allofthesesuccessfulapplicationsshowtheoriginality,validityandgeneralityoftheHAMfornonlinearproblems,andencourageustoapplytheHAMtoattacksomefamous,challengingnonlinearproblems.SomeoftheseproblemsaresuggestedbelowfortheHAMcommunity,especiallyforbrave,enterprising,youngresearchers.1.5.1.PeriodicsolutionsofchaoticdynamicsystemsItiswellknownthatchaoticdynamicsystemshavetheso-called“but-terflyeffect”[52,53],say,thecomputer-generatednumericalsimulationshavesensitivedependencetoinitialconditions(SDIC).Forexample,thenonlineardynamicsystemofLorenzequations[52]x˙=σ(y−x),(1.31)y˙=rx−y−xz,(1.32)z˙=xy−bz,(1.33)haschaoticsolutionincaseofr=28,b=8/3andσ=10formostofgiveninitialconditionsx0,y0,z0ofx,y,zatt=0.However,forsomespecialinitialconditionssuchasx0=−13.7636106821,y0=−19.5787519424,z0=27;x0=−9.1667531454,y0=−9.9743951128,z0=27;x0=−13.5683173175,y0=−19.1345751139,z0=27,theabovedynamicsystemofLorenzequationhasunstableperiodicsolu-tions,asreportedbyViswanath[54].Aperiodicsolutionu(t)withtheperiodThasthepropertyu(t)=u(t+nT)
30October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.122S.-J.Liaoforarbitrarytimet≥0andarbitraryintegersn,evenift→+∞andn→∞.Thispropertycannotbecheckedstrictlybymeansofnumericalapproaches,sinceallnumericalintegrationsimulationsaregainedinafi-niteintervaloftime.Naturally,aperiodicsolutionshouldbeexpressedanalyticallybyperiodicbasefunctionssuchastrigonometricfunctions.So,theoreticallyspeaking,itisinherentlybettertouseanalyticapproximationmethodstosearchforperiodicsolutionsofchaoticdynamicsystemsthannumericalones.Infact,asillustratedbyLiaoinChapter13ofhisbook[10],theHAMcanbeemployedtogainperiodicsolutionofnonlineardynamicsystems.CanweemploytheHAMtogaintheabove-mentionedunstableperiodicsolutionsofLorenzequationfoundbyViswanath[54]?Moreimportantly,itwouldbeveryinterestingiftheHAMcouldbeemployedtofindsomenewperiodicsolutionsofLorenzequationwithphysicalparametersleadingtochaos!ThisismainlybecauseLorenzequationisoneofthemostfamousonesinnonlineardynamicsandnonlinearscience.1.5.2.PeriodicorbitsofNewtonianthree-bodyproblemLetusconsideroneofthemostfamousprobleminmechanicsandappliedmathematics:theNewtonianthree-bodyproblem,say,themotionofthreecelestialbodiesundertheirmutualgravitationalattraction.Letx1,x2,x3denotethethreeorthogonalaxes.Thepositionvectoroftheithbodyisexpressedbyri=(x1,i,x2,i,x3,i),wherei=1,2,3.LetTandLdenotethecharacteristictimeandlengthscales,andmithemassoftheithbody,respectively.UsingNewtoniangravitationlaw,themotionofthethreebodiesaregovernedbythecorrespondingnon-dimensionalequationsX3(xk,j−xk,i)x¨k,i=ρj3,k=1,2,3,(1.34)Ri,jj=1,j6=iwhere"#1/2X3R=(x−x)2(1.35)i,jk,jk,ik=1andmiρi=,i=1,2,3(1.36)m1denotestheratioofthemass.
31October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod23AccordingtoH.Poincar´e,orbitsofthree-bodyproblemareunintegrableingeneral.Althoughchaoticorbitsofthree-bodyproblemswidelyexist,threefamiliesofperiodicorbitswerefound:(1)theLagrange–Eulerfamily,datingbacktotheanalyticalsolutionsinthe18thcentury(onerecentorbitwasgivenbyMoore[55]);(2)theBroucke–Hadjidemetriou–H´enonfamily,datingbacktothemid-1970s[56–61];(3)theFigure-8family,discoveredin1993byMoore[55]andextendedtotherotatingcases[62–65].Notethatnearlyallofthesereportedperiodicorbitsareplanar.In2013,SuvakovandDmitraˇsinovi´c[66]foundthatthereexistfourclassesofplanarˇperiodicorbitsofthree-bodyproblem,withtheabovethreefamiliesbelong-ingtooneclass.Besides,theyreportedthreenewclassesofplanarperiodicorbitsandgavethecorrespondinginitialconditionsforeachclass.Forthedetailsoftheir15planarperiodicorbits,pleaserefertothegallery[67].SuvakovandDmitraˇsinovi´c[66]foundthesenewclassesofplanarperi-ˇodicorbitsbymeansofaniterativenumericalintegrationapproachwithoutusingmultipleprecision.So,itisunknownwhetherornotthenumericalsimulationsdepartthecorrespondingperiodicorbitsforratherlargetime,i.e.t→∞.Asmentionedbefore,itisbetterandmorenaturaltoexpressaperiodicsolutionu(t)withtheperiodTinseriesofperiodicbasefunctions(withthesameperiodT)sothatu(t)=u(t+nT)canholdforarbitraryintegernandarbitrarytimetevenift→∞.Thus,itisvaluabletoapplytheHAMtodoublecheckallofthereportedperiodicorbitsin[66],andmoreimportantly,tofindsomecompletelynewperiodicorbits!Notethatnearlyalloftheperiodicorbitsofthree-bodyproblemre-porteduptonowareplanar.Therefore,itisvaluableandinterestingiftheHAMcanbeappliedtofindsomeperiodicorbitsofNewtonianthree-bodyproblems,whicharenotplanar,i.e.threedimensional.Mathematicallyspeaking,weshoulddeterminesuchunknowninitialpositionsr1,r2,r3,unknowninitialvelocitiesr˙1,r˙2,r˙3andunknowncorrespondingmass-ratiosρ1,ρ2,ρ3ofthreebodiesintheframeoftheHAMthatEqs.(1.34)haveperiodicsolutionxk,i(t)=xk,i(t+nT)forarbitrarytimetandintegern,whereTistheunknowncorrespondingperiodtobedetermined,andi,k=1,2,3.Thisisavaluable,interestingbutchallengingproblemfortheHAMcommunity.
32October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.124S.-J.Liao1.5.3.ViscousflowpastasphereOneofthemostfamous,classicalprobleminfluidmechanicsisthesteady-stateviscousflowpastasphere[68–75],governedbytheNavier-Stokesequation,i.e.asystemofnonlinearpartialdifferentialequations.Considerthesteady-stateviscousflowpastasphereinauniformstream.Howlargeisthedragofthesphereduetotheviscosityoffluid?Tostudythesteady-stateviscousflowpastasphere,thesphericalco-ordinates~r=(r,θ,φ)isoftenused.Sincetheproblemhasaxialsymmetry,onecanusetheStokesstreamfunctionψ(r,θ)definedthroughthefollowingrelations:11vr=2ψθ,vθ=−ψr,vφ=0.(1.37)rsin(θ)rsin(θ)Thestreamfunctionψ(r,θ)isgovernedbythedimensionlessequationR∂(ψ,D2ψ)D4ψ=+2D2ψLψ,(1.38)r2∂(r,µ)subjecttotheboundaryconditionsψ(1,µ)=0,(1.38a)∂ψ(r,µ)=0,(1.38b)∂rr=1ψ(r,µ)12lim=(1−µ),(1.38c)r→∞r22whereR=aU∞/νistheReynoldsnumberandµ≡cos(θ),(1.39)∂21−µ2∂2D2≡+,(1.40)∂r2r2∂µ2µ∂1∂L≡+.(1.41)1−µ2∂rr∂µHere,adenotestheradiusofthesphereandU∞theuniformstreamvelocityatinfinity,respectively,accordingtothenotationofProudmanandPearson[71].AsmentionedbyLiao[74],thedragcoefficientreadsZ124∂ψCD=−pµ+dµ,(1.42)R∂r2−1r=1wherethepressurepisgivenbyZ11∂3ψp=−dµ.(1.43)µ(1−µ2)∂r3r=1
33October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod25Unfortunately,neitherthelinearizationmethod[68–70]northepertur-bationtechniques[71,72]canprovideananalyticapproximationofdragcoefficientCDvalidforRd>3,whereRd=dU∞/ν=2Rforthediam-eterdofthesphere.Especially,the3rd-ordermultiple-scaleperturbationapproximationofCDgivenbyChesterandBreach[72]wasvalideveninasmallerintervalofReynoldsnumberthanthe2nd-ordermultiple-scaleperturbationresultofProudmanandPearson[71].Thisimpliestheinva-lidityofperturbationmethodsforthisfamousproblem.So,“theideaofusingcreepingflowtoexpandintothehighReynoldsnumberregionhasnotbeensuccessful”,aspointedoutbyWhiteinhistextbook[76].Besides,themethodofrenormalizationgroupcannotessentiallymodifytheseanalyticresults[75],either.In2002,Liao[74]employedtheHAMtosolvethesteady-stateviscousflowpastasphereandgainedaanalyticapproximationofdragcoefficientCD,whichagreewellwithexperimentaldatainaconsiderablylargerin-tervalRd≤30.However,thecorrespondingexperimentsindicatethatthesteady-stateviscousflowpastasphereexistsuntilRd≈100.So,strictlyspeaking,thisHAMresultgivenin[74]isnotsatisfactory.Theoreticallyspeaking,itisveryinterestingandvaluableifonecangiveanaccurateenoughanalyticresultofthedragcoefficientCDvalidforthesteady-stateviscousflowpastasphereuptoRd≈100,mainlybecauseitisoneofthemostfamous,classicalproblemsinfluidmechanicswithahistoryofmorethan150year!Canwesolvethisfamous,classicalproblembymeansoftheHAM?1.5.4.ViscousflowpastacylinderThesteady-stateviscousflowpastaninfinitecylinderisalsooneofthemostfamous,classicalproblemsinfluidmechanicswithalonghistory.Forthesteady-stateviscousflowpastaninfinitecylinder,itisnaturaltousecylindricalcoordinates~r=(r,θ,z).Sincetheproblemistwodimensional,itisconvenienttousetheLagrangianstreamfunctionψ(r,θ)definedbyProudmanandPearson[71]:1∂ψ∂ψur=,uθ=−,uz=0.(1.44)r∂θ∂rThestreamfunctionψ(r,θ)isgovernedbyR∂(ψ,∇2)∇4ψ(r,θ)=−r,(1.45)rr∂(r,θ)
34October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.126S.-J.Liaosubjecttotheboundaryconditionsψ(r=1,θ)=0,(1.45a)∂ψ(r,θ)=0,(1.45b)∂rr=1ψ(r,θ)lim=sin(θ),(1.45c)r→∞rwhere∂21∂1∂22422∇r≡∂r2+r∂r+r2∂θ2,∇r≡∇r∇r.Here,R=aU∞/νistheReynoldsnumber,aandU∞denotetheradiusofcylinderandtheuniformstreamvelocityatinfinity,respectively.Asreviewedin[75],neitherthelinearizationmethodnorperturbationtechniquecangivegoodanalyticapproximationofthedragcoefficientCDofacylinderforR≥3.Infact,eventhemethodofrenormalizationgroupcannotmodifytheseresultsgreatly[75].So,itisstillanopenquestion.Theoreticallyspeaking,itisvaluabletogainanaccurateanalyticexpres-sionofdragcoefficientCDvalidforlargeReynoldsnumberuptoR≈40,beyondwhichtheperiodicVonKarm´anvortexoccurs.Thisismainlybe-causeitisoneofthemostfamous,historicalprobleminfluidmechanics.Canthisfamous,classicalproblembesolvedbymeansoftheHAM?1.5.5.NonlinearwaterwavesTheHAMhasbeensuccessfullyappliedtosolvesomenonlinearwaveequa-tions.Especially,intheframeoftheHAM,themultipleequilibrium-statesofresonantwavesindeepwater[33]andinfinitewaterdepth[34]werediscoveredforthefirsttime,tothebestoftheauthor’sknowledge.Thus,theHAMprovidesusaconvenienttooltoinvestigatesomecomplicatedwaveproblems.Strictlyspeaking,waterwavesaregovernedbyEulerequationwithtwononlinearboundaryconditionssatisfiedonanunknownfreesurface,whichhoweverareratherdifficulttosolveingeneral.BasedontheexactEulerequation,somesimplifiedwavemodelsforshallowwaterwaves,suchastheKdVequation[77],Boussinesqequation[78],Camassa–Holm(CH)equa-tion[79],andsoon,arederivedbyassumingtheexistenceofsomesmallphysicalparametersinshallowwater.AlthoughtheseshallowwaterwaveequationsaremuchsimplerthantheexactEulerequation,theycanwellex-
35October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod27plainmanyphysicalphenomena,suchassolitonwaves,wavepropagationsandinteractionsinshallowwater,wavebreaking,andsoon.Forexample,thecelebratedCamassa–Holm(CH)equation[79]ut+2ωux−uxxt+3uux=2uxuxx+uuxxx,(1.46)subjecttotheboundaryconditionu=0,ux=0,uxx=0,asx→±∞,(1.47)canmodelbothphenomenaofsolitoninteractionandwavebreaking(see[80]),whereu(x,t)denotesthewaveelevation,x,tarethetemporalandspatialvariables,ωisaconstantrelatedtothecriticalshallowwaterwavespeed,thesubscriptdenotesthepartialdifferentiation,respectively.Mathe-matically,theCHequationisintegrableandbi-Hamiltonian,thuspossessesaninfinitenumberofconservationlawsininvolution[79].Inaddition,itisassociatedwiththegeodesicflowontheinfinitedimensionalHilbertmani-foldofdiffeomorphismsofline(see[80]).Thus,theCHequation(1.46)hasmanyintriguingphysicalandmathematicalproperties.AspointedoutbyFushssteiner[81],theCHequation(1.46)even“hasthepotentialtobecomethenewmasterequationforshallowwaterwavetheory”.Especially,whenω=0,theCHequation(1.46)hasthepeakedsolitarywaveu(x,t)=cexp(−|x−ct|),whichwasfoundfirstbyCamassaandHolm[79].Thefirstderivativeofthepeakedsolitarywaveisdiscontinuousatthecrestx=ct.LiketheCHequation,manyshallowwaterequationsadmitpeakedand/orcuspedsolitarywaves.Theseequationswithpeakedand/orcuspedsolitarywaveshavebeenwidelyinvestigatedmathematically,andthousandsofrelatedar-ticleshavebeenpublished.However,tothebestoftheauthor’sknowledge,peakedandcuspedsolitarywaveshaveneverbeengaineddirectlyfromtheexactEulerequation!Thisisverystrange.Logicallyspeaking,sincethesesimplifiedequations(liketheCHequation)aregoodenoughapprox-imationsoftheEulerequationinshallowwater,theexactEulerequationshouldalsoadmitthepeakedand/orcuspedsolitarywavesaswell.Canwegainsuchkindofpeakedand/orcuspedsolitarywavesoftheexactwaveequationbymeansoftheHAM,iftheyindeedexist?Eitherpositiveornegativeanswerstothisquestionhaveimportantscientificmean-ings.Ifsuchkindofpeakedsolutionsoftheexactwaveequationindeedexist,itcangreatlyenrichanddeepenourunderstandingsaboutpeaked
36October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.128S.-J.Liaosolitarywaves.Ifthepeakedsolitarywavesgivenbytheexactwaveequa-tionexistsmathematicallybutisimpossibleinphysics,wehadtocheckthephysicalvalidityofthepeakedsolitarywaves.So,thisisaninterestingandvaluablework,althoughwithgreatchallenge.Forsomeattemptsinthisdirection,pleaserefertoLiao[82],whoproposedageneralizedwavemodelbasedonthesymmetryandthefullynonlinearwaveequations,whichadmitsnotonlythetraditionalwaveswithsmoothcrestbutalsopeakedsolitarywaves.ItisfoundthatthepeakedsolitarywavessatisfyKelvin’stheoremeverywhere.Besides,thesepeakedsolitarywavesincludethefa-mouspeakedsolitarywavesoftheCamassa–Holmequation.So,thegener-alizedwavemodel[82]isconsistentwiththetraditionalwavetheories.Itisfound[82]thatthepeakedsolitarywaveshavesomeunusualcharacter-isticsquitedifferentfromthetraditionalones,althoughitisstillanopenquestionwhetherornottheyarereasonableinphysicsiftheviscosityoffluidandthesurfacetensionareconsidered.Inaddition,theso-called“roguewave”[83,84]isahottopicofnonlinearwaves.Certainly,itisvaluabletoapplytheHAMtodosomeinvestigationsinthisfield.Insummary,itistruethattheproblemssuggestedaboveareindeeddif-ficult,butveryvaluableandinterestingintheory.Infact,therearemanysuchkindofinterestingbutdifficultproblemsinscience,engineeringandfinance.Itshouldbeemphasizedthat,unlikeallotheranalyticapprox-imationmethods,theHAMprovidesusgreatfreedomandflexibilitytochooseequation-typeandsolutionexpressionofhigh-orderapproximationequations,andbesidesaconvenientwaytoguaranteetheconvergenceofsolutionseries.AspointedoutbyGeorgCantor(1845–1918),“theessenceofmathematicsliesentirelyinitsfreedom”.Hopefully,thegreatfreedomandflexibilityoftheHAMmightcreatesomenovelideasandinspiresomebrave,enterprising,youngresearcherswithstimulatedimaginationtoat-tackthemwithsatisfactory,muchbetterresults.Chancealwaysstayswithchallenges!AcknowledgmentThisworkispartlysupportedbyNationalNaturalScienceFoundationofChina(ApprovalNo.11272209),theFoundationforShanghaiLeadingScientists,andStateKeyLaboratoryofOceanEngineering(ApprovalNo.GKZD010056).
37October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod29References[1]A.H.Nayfeh,PerturbationMethods.JohnWiley&Sons,NewYork,1973.[2]M.VanDyke,PerturbationMethodsinFluidMechanics.Parabolic,Stan-ford,1975.[3]E.J.Hinch,PerturbationMethods.CambridgeUniversityPress,Cambridge,1991.[4]J.A.Murdock,Perturbations:TheoryandMethods.JohnWiley&Sons,NewYork,1991.[5]A.M.Lyapunov,GeneralProblemonStabilityofMotion(Englishtransla-tion).Taylor&Francis,London,1992(theoriginalonewaspublishedin1892inRussian).[6]G.Adomian,Nonlinearstochasticdifferentialequations,J.Math.Anal.Appl.55:441–452(1976).[7]G.Adomian,Areviewofthedecompositionmethodandsomerecentresultsfornonlinearequations,Comput.Math.Appl.21:101–127(1991).[8]G.Adomian,SolvingFrontierProblemsofPhysics:TheDecompositionMethod.KluwerAcademicPublishers,Boston,1994.[9]S.J.Liao,Ontheproposedhomotopyanalysistechniquesfornonlinearprob-lemsanditsapplication.Ph.D.dissertation,ShanghaiJiaoTongUniversity,1992.[10]S.J.Liao,BeyondPerturbation:IntroductiontotheHomotopyAnalysisMethod.Chapman&Hall/CRCPress,BocaRaton,2003.[11]S.J.Liao,HomotopyAnalysisMethodinNonlinearDifferentialEquations.Springer&HigherEducationPress,Heidelberg,2012.[12]S.J.Liao,Anexplicit,totallyanalyticapproximationofBlasius’viscousflowproblems,Int.J.Non-LinearMech.34:759–778(1999).[13]S.J.Liao,Onthehomotopyanalysismethodfornonlinearproblems,Appl.Math.Comput.147:499–513(2004).[14]S.J.LiaoandY.Tan,Ageneralapproachtoobtainseriessolutionsofnon-lineardifferentialequations,Stud.Appl.Math.119:297–354(2007).[15]S.J.Liao,Notesonthehomotopyanalysismethod:somedefinitionsandtheorems,Commun.NonlinearSci.Numer.Simulat.14:983–997(2009).[16]S.J.Liao,Anoptimalhomotopy-analysisapproachforstronglynonlineardifferentialequations,Commun.NonlinearSci.Numer.Simulat.15:2315–2332(2010).[17]S.J.Liao,OntherelationshipbetweenthehomotopyanalysismethodandEulertransform,Commun.NonlinearSci.Numer.Simulat.15:1421–1431,2010.[18]A.V.Karmishin,A.T.ZhukovandV.G.Kolosov,MethodsofDynamicsCal-culationandTestingforThin-walledStructures(inRussian).Mashinostroye-nie,Moscow,1990.[19]J.H.He,Anapproximatesolutiontechniquedependinguponanartificialparameter,Commun.NonlinearSci.Numer.Simulat.3:92–97(1998).[20]M.SajidandT.Hayat,ComparisonofHAMandHPMmethodsfornon-linearheatconductionandconvectionequations.Nonlin.Anal.B.9:2296–
38October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.130S.-J.Liao2301(2008).[21]S.X.LiangandD.J.Jeffrey,Comparisonofhomotopyanalysismethodandhomotopyperturbationmethodthroughanevaluationequation,Commun.NonlinearSci.Numer.Simulat.14:4057–4064(2009).[22]V.MarincaandN.Heri¸sanu:Applicationofoptimalhomotopyasymptoticmethodforsolvingnonlinearequationsarisinginheattransfer,Int.Com-mun.HeatMass.35:710–715(2008).[23]R.A.VanGorder:GaussianwavesintheFitzhugh-Nagumoequationdemon-strateoneroleoftheauxiliaryfunctionH(x)inthehomotopyanalysismethod,Commun.NonlinearSci.Numer.Simulat.17:1233–1240(2012).[24]K.VajraveluandR.A.VanGorder,NonlinearFlowPhenomenaandHomo-topyAnalysis:FluidFlowandHeatTransfer.Springer,Heidelberg,2012.[25]Y.L.Zhao,Z.L.LinandS.J.Liao,AniterativeHAMapproachfornon-linearboundaryvalueproblemsinasemi-infinitedomain,Comput.Phys.Commun.184:2136–2144(2013).[26]S.S.Motsa,P.SibandaandS.Shateyi,AnewspectralhomotopyanalysismethodforsolvinganonlinearsecondorderBVP,Commun.NonlinearSci.Numer.Simulat.15:2293–2302(2010).[27]S.S.Motsa,P.Sibanda,F.G.AuadandS.Shateyi,AnewspectralhomotopyanalysismethodfortheMHDJeffery-Hamelproblem,Computer&Fluids.39:1219–1225(2010).[28]W.S.DonandA.Solomonoff,AccuracyandspeedincomputingtheCheby-shevcollocationderivative.SIAMJ.Sci.Comput.16:1253–1268(1995).[29]S.AbbasbandyandE.Shivanian,Predictionofmultiplicityofsolutionsofnonlinearboundaryvalueproblems:Novelapplicationofhomotopyanalysismethod,Commun.NonlinearSci.Numer.Simulat.15:3830–3846(2010).[30]S.AbbasbandyandE.Shivanian,Predictorhomotopyanalysismethodanditsapplicationtosomenonlinearproblems,Commun.NonlinearSci.Numer.Simulat.16:2456–2468(2011).[31]S.J.Liao,Anewbranchofsolutionsofboundary-layerflowsoveranimper-meablestretchedplate,Int.J.HeatMassTran.48:2529–2539(2005).[32]S.J.LiaoandE.Magyari,Exponentiallydecayingboundarylayersaslimit-ingcasesoffamiliesofalgebraicallydecayingones,Z.angew.Math.Phys.57,777–792(2006).[33]S.J.Liao,Onthehomotopymultiple-variablemethodanditsapplicationsintheinteractionsofnonlineargravitywaves,Commun.NonlinearSci.Numer.Simulat.16:1274–1303(2011).[34]D.L.Xu,Z.L.Lin,S.J.LiaoandM.Stiassnie,Onthesteady-statefullyresonantprogressivewavesinwateroffinitedepth,J.FluidMech.710:379–418(2012).[35]S.H.ParkandJ.H.Kim,Homotopyanalysismethodforoptionpricingunderstochasticvolatility,Appl.Math.Lett.24:1740–1744(2011).[36]S.H.ParkandJ.H.Kim,Asemi-analyticpricingformulaforlookbackop-tionsunderageneralstochasticvolatilitymodel,StatisticsandProbabilityLetters,83:2537–2543(2013).[37]H.B.Keller,Numericalsolutionofbifurcationandnonlineareigenvalueprob-
39October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod31lems.InP.Rabinowitz,ApplicationofBifurcationTheory(p.359),AcademicPress,NewYork,1977.[38]H.B.Keller,GlobalhomotopiesandNewtonmethods.InC.deBoor,etal.,RecentAdvancesinNumericalAnalysis(p.73),AcademicPress,NewYork,1978.[39]E.Lahaye,Unem´ethodeder´esolutiond’unecat´egoried´equationstranscen-dantes.ComptesRendusdel’Acad´emiedesSciencesParis,198:1840–1842(1934).[40]T.L.WayburnandJ.D.Seader,Homotopycontinuationmethodsforcomputer-aidedprocessdesign,ComputersinChemicalEngineering.11:7(1987).[41]B.C.Eaves,Homotopiesforcomputationoffixedpoints,Math.Program-ming.3:1–22(1972).[42]S.N.Chow,J.Mallet-ParetandJ.A.Yorke,Findingzerosofmaps:Homo-topymethodsthatareconstructivewithprobabilityone,Math.Comp.32:887–899(1978).[43]L.T.Watson,Solvingthenonlinearcomplementarityproblembyahomotopymethod,SIAMJ.ControlOptim.17:36–46(1979).[44]J.C.Alexander,T.Y.LiandJ.A.Yorke,Piecewisesmoothhomotopies,inHomotopyMethodsandGlobalConvergence,pp.1–14,B.C.Eaves,F.J.Gould,H.O.PeitgenandM.J.Todd,eds.,Plenum,NewYork,1983.[45]L.T.Watson,M.Sosonkina,R.C.Melville,A.P.MorganandH.F.Walker,Algorithm777:HOMPACK90:AsuiteofFortran90codesforgloballyconvergenthomotopyalgorithms,ACMTrans.Math.Software.23:514–549(1997).[46]T.Y.Li,Numericalsolutionofmultivariatepolynomialsystemsbyhomotopycontinuationmethods,ActaNumerica.6:399–436(1997).[47]T.Y.Li,Solvingpolynomialsystemsbythehomotopycontinuationmethod,HandbookofNumericalAnalysis,Vol.XI,pp.209–304,editedbyP.G.Cia-rlet,North-Holland,Amsterdam,2003.[48]T.Lee,C.TsaiandT.Y.Li,HOM4PS-2.0:Asoftwarepackageforsolv-ingpolynomialsystemsbythepolyhedralhomotopycontinuationmethod,Computing.83:109–133(2008).[49]T.Y.LiandC.H.Tsai,HOM4PS-2.0para:ParallelizationofHOM4PS-2.0forsolvingpolynomialsystems,ParallelComputing.35:226–238(2009).[50]S.P.Zhu,AnexactandexplicitsolutionforthevaluationofAmericanputoptions,Quant.Financ.6:229–242(2006).[51]J.Cheng,S.P.ZhuandS.J.Liao,AnexplicitseriesapproximationtotheoptimalexerciseboundaryofAmericanputoptions,Commun.NonlinearSci.Numer.Simulat.15:1148–1158(2010).[52]E.N.Lorenz,Deterministicnon-periodicflows,J.Atmos.Sci.20:130–141(1963).[53]E.N.Lorenz,TheEssenceofChaos.UniversityofWashingtonPress,Seattle,1995.[54]D.Viswanath,ThefractalpropertyoftheLorenzattractor,PhysicaD.190:115–128(2004).
40October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.132S.-J.Liao[55]C.Moore,Braidsinclassicaldynamics,Phys.Rev.Lett.70:3675–3679(1993).[56]R.BrouckeandD.Boggs,PeriodicorbitsinthePlanarGeneralThree-BodyProblem,Celest.Mech.11:13–38(1975).[57]J.D.HadjidemetriouandT.Christides,Familiesofperiodicorbitsintheplanarthree-bodyproblem,Celest.Mech.12:175–187(1975).[58]J.D.Hadjidemetriou,Thestabilityofperiodicorbitsinthethree-bodyprob-lem,Celest.Mech.12:255–276(1975).[59]R.Broucke,Onrelativeperiodicsolutionsoftheplanargeneralthree-bodyproblem,Celest.Mech.12:439–462(1975).[60]M.H´enon,Afamilyofperiodicsolutionsoftheplanarthree-bodyproblem,andtheirstability,Celest.Mech.13:267–285(1976).[61]M.H´enon,Stabilityofinterplaymotions,Celest.Mech.15:243–261(1977).[62]M.Nauenberg,Periodicorbitsforthreeparticleswithfiniteangularmomen-tum,Phys.Lett.A.292:93–99(2001).[63]A.Chenciner,J.FejozandR.Montgomery,RotatingEights:I.thethreeΓifamilies,Nonlinearity.18:1407(2005).[64]R.Broucke,A.ElipeandA.Riaguas,Onthefigure-8periodicsolutionsinthethree-bodyproblem,Chaos,SolitonsFractals.30:513–520(2006).[65]M.Nauenberg,Continuityandstabilityoffamiliesoffigureeightorbitswithfiniteangularmomentum,Celest.Mech.97:1–15(2007).[66]M.SuvakovandV.Dmitraˇsinovi´c,Threeclassesofnewtonianthree-bodyˇplanarperiodicorbits,Phys.Rev.Lett.110:114301(2013).[67]http://suki.ipb.ac.rs/3body/.[68]G.G.Stokes,Ontheeffectoftheinternalfrictionoffluidsonthemotionofpendulums,Trans.CambridgePhilos.Soc.9:8(1851).[69]C.W.Oseen,UberdieStokesscheFormelunddieverwandteAufgabeinder¨Hydrodynamik,Arkiv.Mat.Astron.Phsik.6:29(1910).[70]S.Goldstein,ThesteadyflowofviscousfluidpastafixedsphericalobstacleatsamllReynoldsnumbers,Proc.R.Soc.London,Ser.A.123:225(1929).[71]I.ProudmanandJ.R.A.Pearson,ExpansionsatsmallReynoldsnumbersfortheflowpastasphereandacircularcylinder,J.FluidMech.2:237(1957).[72]W.ChesterandD.R.Breach,OntheflowpastasphereatlowReynoldsnumber,J.FluidMech.37:751–760(1969).[73]E.R.Lindgren,ThemotionofasphereinanincompressibleviscousfluidatReynoldsnumbersconsiderablylessthanone,Phys.Scr.60:97(1999).[74]S.J.Liao,Ananalyticapproximationofthedragcoefficientfortheviscousflowpastasphere,Int.J.Non-LinearMech.37:1–18(2002).[75]J.VeyseyandN.Goldenfeld,Simpleviscousflow:Fromboundarylayerstotherenormalizationgroup,Rev.Mod.Phys.79:883(2007).[76]F.M.White,ViscousFluidFlow.McGraw-Hill,NewYork,1991.[77]D.J.KortewegandG.deVries,Onthechangeofformoflongwavesad-vancinginarectangularcanal,andonanewtypeoflongstationarywaves,Phil.Mag.39:422–443(1895)[78]J.Boussinesq,Th´eoriedesondesetdesremousquisepropagentlelongd’un
41October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.1ChanceandChallenge:ABriefReviewofHomotopyAnalysisMethod33canalrectangulairehorizontal,encommuniquantauliquidecontenudanscecanaldesvitessessensiblementpareillesdelasurfaceaufond,JournaldeMath’ematiquesPuresetAppliqu´ees.17:55–108(1872).[79]R.CamassaandD.D.Holm,Anintegrableshallowwaterequationwithpeakedsolitons,Phys.Rev.Lett.71:1661–1664(1993).[80]A.Constantin,Existenceofpermanentandbreakingwavesforashallowwaterequation:ageometricapproach,Ann.Inst.Fourier,Grenoble.50:321–362(2000).[81]B.Fuchssteiner,Sometricksfromthesymmetry-toolboxfornonlinearequa-tions:generalizationsoftheCamassa-Holmequation,PhysicaD.95:296–343(1996).[82]S.J.Liao,Dopeakedsolitarywaterwavesindeedexist?ArXiv:1204.3354,2012.[83]I.DidenkulovaandE.Pelinovsky,Roguewavesinnonlinearhyperbolicsys-tems(shallow-waterframework),Nonlinearity.24:R1–R18(2011).[84]A.CaliniandC.M.Schober,DynamicalcriteriaforroguewavesinnonlinearSchr¨odingermodels,Nonlinearity.25:R99–R116(2012).
42May2,201314:6BC:8831-ProbabilityandStatisticalTheoryPST˙wsThispageintentionallyleftblank
43October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.2Chapter2PredictorHomotopyAnalysisMethod(PHAM)SaeidAbbasbandy∗andElyasShivanian†DepartmentofMathematicsImamKhomeiniInternationalUniversityGhazvin,Iran∗abbasbandy@yahoo.com,†shivanian@sci.ikiu.ac.irThischapterintroducesamethodtopredictthemultiplicityoftheso-lutionsofnonlinearboundaryvalueproblems.Thisprocedurecanbeeasilyappliedonnon-linearordinarydifferentialequationswithbound-aryconditionssothatitcalculateseffectivelytheallbranchesofthesolutions(ontheconditionthat,thereexistsuchsolutionsfortheprob-lem)analyticallyatthesametime.Inthismanner,forpracticaluseinscienceandengineering,thismethodmightgivenewunfamiliarclassofsolutionswhichisoffundamentalinterest.Contents2.1.Preliminaries....................................362.2.Descriptionofthemethod.............................372.2.1.Zeroth-orderdeformationequation.....................382.2.2.High-orderdeformationequation......................392.2.3.Predictionofthemultiplesolutions....................402.3.Convergenceanalysis................................482.4.Someillustrativemodels..............................542.4.1.Nonlinearproblemarisinginheattransfer................552.4.2.StronglynonlinearBratu’sequation....................602.4.3.Nonlinearreaction-diffusionmodel.....................662.4.4.Mixedconvectionflowsinaverticalchannel...............712.5.Concludingremarks.................................79References.........................................8135
44October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.236S.AbbasbandyandE.Shivanian2.1.PreliminariesManyofthemathematicalmodelingofthephysicalphenomenainscienceandengineeringoftenleadtononlineardifferentialequations.Therearealotofmethods,fromthepastuptonow,togivenumericallyapproximatesolutionsofnonlineardifferentialequationssuchasEulermethod,Runge-Kuttamethod,multistepmethod,Taylorseriesmethod,Hybridmethods,familyoffinitedifferencemethods[1,2],familyoffiniteelementmeth-ods[3],meshlessmethods,differentialquadrature,spectralmethods[4–6]etc.Therearealsomanyofmethodswhichgiveanalyticallyapproximatesolutionslikeforexampleperturbationmethods[7,8],theartificialsmallparametermethod[9],theδ-expansionmethod[10],andrecentlynotewor-thytechniquesconsistsoftheAdomiandecompositionmethod[11],thevariationaliterationmethod[12]andsoon.Inthisregard,onemayaskthequestion:Dotheapproximatemethodsenabletopredictmultiplicityofsolutionsofthenonlineardifferentialequations?Intheotherwords,canweforecastexistenceofmultiplesolutionsofnonlinearequationsbyap-proximatemethodsandatthesametimeobtainallbranchesofsolutions?Answertothisquestionissomehowdifficultbyconcerningthattheapprox-imatemethodsusuallyconvergetoonesolutionbyoneinitialguessthatisexactlymeaningof“convergence”.Nevertheless,thepresentchapterisgoingtoinfractthisconventionbyintellectualusinghomotopyanalysismethod.Thepurposeofthepresentchapteristointroduceamethod,probablyforthefirsttime,topredictthemultiplicityofthesolutionsofnonlinearboundaryvalueproblems.Thisprocedurecanbeeasilyappliedonnonlin-earordinarydifferentialequationswithboundaryconditions.Thismethod,aswillbeseen,besidesanticipatingofmultiplicityofthesolutionsofthenonlineardifferentialequations,calculateseffectivelyallbranchesoftheso-lutions(ontheconditionthat,thereexistsuchsolutionsfortheproblem)analyticallyatthesametime.Inthismanner,forpracticaluseinscienceandengineering,thismethodmightgivenewunfamiliarclassofsolutionswhichisoffundamentalinterestandfurthermore,theproposedapproachconvincestoapplyitonnonlinearequationsbytoday’spowerfulsymbolicsoftwareprogramssothatitdoesnotneedtediousstagesofevaluationandcanbeusedwithoutstudyingthewholetheory.Infact,thistechniquehasnewpointofviewtohomotopyanalysismethod.AsitiswellknownintheframeofHAM,theconvergence-controllerparameterplaysimportantroletoguaranteetheconvergenceofthesolutionsofnonlineardifferentialequa-
45October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.2PredictorHomotopyAnalysisMethod(PHAM)37tions.Itisshownthattheconvergence-controllerparameterplaysafunda-mentalroleinthepredictionofmultiplicityofsolutionsandallbranchesofsolutionsareobtainedsimultaneouslybyoneinitialapproximationguess,oneauxiliarylinearoperatorandoneauxiliaryfunction.Notwithstandingboththeperturbationandnon-perturbationmethodssuchastheartificialsmallparametermethod,theδ-expansionmethod,theAdomiandecompositionmethodandvariationaliterationmethodcannotprovideuswithasimplewaytoadjustandcontroltheconvergenceregionandrateofgivenapproximateseries,theHAMcanguaranteetheconver-genceoftheseriessolutionsbyconvergence-controllerparameter~.Infact,thistraitmakesHAMtobedifferentfromtheothersanalyticaltechniqueswhichareusedtoapproachtojustonesolutionandsothatpossiblytolosetheothersolutions.Thereforeinthisway,presentproceduremaybegener-atesnewclassofsolutionsforfurtherphysicalinterpretationsinengineeringandsciences.Thelegitimacyandreliabilityofthemethodischeckedbyitsapplica-tiontofourimportantnonlinearequationsnamelynonlinearheattransferequation[13,14],stronglynonlinearBratu’sequation[15,16],nonlinearreaction-diffusionequation[17]andthemodelofmixedconvectionflowsinaverticalchannel[18,19].Allthesenonlinearproblemsadmitmultiplesolutionsforsomevaluesofparametersoftheequations.2.2.DescriptionofthemethodToillustratetheprocedureconsiderthefollowingnonlineardifferentialequation:N[u(r)]=0,r∈Ω,(2.1)withboundaryconditions∂uBu,=0,r∈Γ,(2.2)∂nwhereNisgeneralnonlinearoperator,Bisaboundaryoperator,andΓistheboundaryofthedomainΩ.Thecrucialstepofthetechniqueisthattheboundaryvalueproblem(2.1)and(2.2)shouldbereplacedbyequivalentproblemsothattheconditions(2.2)involveanunknownparameterlikeδ(prescribedparameter)andaresplitto0∂uBu,δ,=0,r∈Γandu(α)=β,(2.3)∂n
46October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.238S.AbbasbandyandE.Shivanianwhereu(α)=βistheforcingconditionthatcomesfromoriginalconditions(2.2).Now,homotopyanalysismethodisappliedontheproblem(2.1)withtheconditions(2.3)exceptforcingconditionasfollows:N[u(r)]=0,r∈Ω,(2.4)0∂uBu,δ,=0,r∈Γ.(2.5)∂n2.2.1.Zeroth-orderdeformationequationWesupposethatallthesolutionsu=u(r)ofproblem(2.4)canbeexpressedbythesetofbasefunctions{ωi(r),i=0,1,2,...}intheformX+∞u=u(r)=anωn(r),(2.6)n=0whereanarecoefficientstobedetermined.Letu0(r,δ)denoteaninitialapproximationguessoftheexactsolutionu(r)whichsatisfiesboundaryconditions(2.5)automatically.Also,asthatiswellknownintheframeofHAM,assume~6=0denoteconvergence-controllerparameter,H(r)6=0anauxiliaryfunction,andLanauxiliarylinearoperator.Nowusingp∈[0,1]asanembeddingparam-eter,weconstructthegeneralzeroth-orderdeformationequationandthecorrespondingboundaryconditionsasfollow:(1−p)L[ϕ(r,δ;p)−u0(r,δ)]=p~H(r)N[ϕ(r,δ;p)],(2.7)0ϕ(r,δ;p)Bϕ(r,δ;p),δ,=0,r∈Γ,(2.8)∂nwhereϕ(r,δ;p)isanunknownfunctiontobedetermined.Whenp=0,thezeroth-orderdeformationequation(2.7)becomesL[ϕ(r,δ;0)−u0(r,δ)]=0,(2.9)whichgivesϕ(r,δ;0)=u0(r,δ).Whenp=1,Eq.(2.7)leadstoN[ϕ(r,δ;1)]=0,(2.10)whichisexactlythesameastheoriginalEq.(2.1)providedthatϕ(r,δ;1)=u(r,δ).
47October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.2PredictorHomotopyAnalysisMethod(PHAM)39Wenowexpandthefunctionϕ(r,δ;p)inaTaylorseriestotheembeddingparameterp.ThisTaylorexpansioncanbewrittenintheformX+∞ϕ(r,δ;p)=u(r,δ)+u(r,δ)pm,(2.11)0mm=1where1∂mϕ(r,δ;p)um(r,δ)=,m=0,1,2,...,+∞.(2.12)m!∂pmAsitiswellknowninduringtheframeofHAM[20,21],whenthelinearoperatorL,theinitialapproximationu0(r,δ),theauxiliaryparameter~6=0,andtheauxiliaryfunctionH(r)6=0arechosenproperly,theseries(2.11)convergesforp=1,andthusX+∞X+∞u(r,δ)=u0(r,δ)+um(r,δ)=anωn(r),(2.13)m=1n=0willbethesolutionofthenonlinearproblem(2.4)and(2.5)aswillbeprovedlater.2.2.2.High-orderdeformationequationAssumethatthelinearoperatorL,theinitialapproximationu0(r,δ),andtheauxiliaryfunctionH(r)6=0arechosenproperly(itisworthmen-tioningherethat~6=0so-calledconvergence-controllerparameterwillbedeterminedlater),theunknownfunctionsum(r,δ)inEq.(2.13)canbedeterminedwiththeaidofthehigh-orderdeformationequationsasfol-−→lows.Atfirstwedefinethevectorun={u0(r),u1(r),...,un(r)}then,differentiatingthezeroth-orderdeformationequation(2.7)mtimeswithrespecttotheembeddingparameterp,dividingitbym!,settingsubse-quentlyp=0andtakingintoaccounttheboundaryconditions(2.8),oneobtainsthemth-orderdeformationequation−→L[um(r,δ)−χmum−1(r,δ)]=~H(r)Rm(um−1,r,δ),(2.14)subjecttotheboundaryconditions∂mϕ(r,δ;p)B0ϕ(r,δ;p),δ,=0,r∈Γ,(2.15)∂pm∂np=0where0,m≤1,χm=(2.16)1,m>1,
48October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.240S.AbbasbandyandE.Shivanianand1∂m−1N[ϕ(r,δ;p)]Rm(−→um−1,r,δ)=(m−1)!∂pm−1p=0hPi∂m−1Nn=+∞u(r,δ)pn1n=0n=.(2.17)(m−1)!∂pm−1p=0Thehigh-orderdeformationequation(2.14)obviously,isjusttheordinarydifferentialequationwithboundarycondition(2.15)and,canbeeasilysolvedbyusingsomesymbolicsoftwareprogramssuchasMathematicaorMaple.Inthisway,startingbyu0(r,δ),weobtainthefunctionsum(r,δ)form=1,2,3,...fromEqs.(2.14)and(2.15)successively.Accordingly,theMth-orderapproximatesolutionoftheproblems(2.4)and(2.5)isgivenby0XMXMu(r,δ)≈UM(r,δ,~)=u0(r,δ)+um(r,δ)=anωn(r).(2.18)m=1n=02.2.3.PredictionofthemultiplesolutionsItisnoteworthytoindicatethatuptothisstage,thelinearoperatorL,theinitialapproximationu0(r,δ),andtheauxiliaryfunctionH(r)6=0havebeenchosenproperlysothattheseriessolutions(2.18)wouldbeconver-gence.However,therearestilltwounknownparametersinseries(2.18)namelyδ(prescribedparameter)and~(convergence-controllerparameter)whichshouldbedetermined.Itisessentialthatexistenceofuniqueormul-tiplesolutionsintermsofthebasicfunctions(2.6)fortheoriginalboundaryvalueproblem(2.1)dependsonthefactwhethertheforcingcondition(2.3)(u(α)=β),admitsuniqueormultiplevaluesfortheformallyintroducedparameterδintheboundaryconditions(2.3).Thisstageiscalledruleofmultiplicityofsolutionsthatisacriterioninordertoknowhowmanysolutionstheboundaryvalueproblem(2.1)admits.Theso-calledruleofmultiplicityofsolutionsisappliedasfollows:ConsidertheMth-orderapproximatesolution(2.18)andsetintheforcingcondition(2.3)(u(α)=β),sothefollowingequationisderivedu(α)≈Um(α,δ,~)=β.(2.19)Theaboveequationhastwounknownparametersnamelyδ,and~whichcontrolstheconvergenceoftheHAMseries(2.18).ItisabasicfeatureofHAMthattheseriessolution(2.18)convergesatr=αonlyinthatrange
49October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.2PredictorHomotopyAnalysisMethod(PHAM)41of~,wheretheparameterδdoesnotchangewiththevariationof~.Thismeansthatintheplotofδasfunctionof~accordingtoEq.(2.19)inanimplicitlyway,intheconvergencerangeoftheseriesu(α)aplateauoccurs.Thenumberofsuchhorizontalplateauswhereδ(~)becomesconstant,givesthemultiplicityofthesolutionsofproblems(2.1)and(2.2),(wecallthispointasruleofmultiplicityofsolutions).Inthefollowing,wegiveaprooftotheaboveassertion.Wenoticethatthereisadirectconnectionbetweentheexistencesofmultiplesolutionsandthenumberofhorizontalplateausintheplotofδ(~),soatfirst,letusdiscussaboutthefundamentalof~-curveanditsrelationtotheTaylor’sseries.Thefollowingtheoremisbasicandsimpletounderstand.Theorem2.1.Iff(x)iscontinuouson[a,b]andisdifferentiablein(a,b)andf(a)=f(b)=0,thenthereexistsatleastonepointξ∈(a,b)suchthat0f(ξ)=0.Proof.SeeRef.[22].Now,wediscussaboutoccurrenceof~-curveforthearbitrarybutsmoothfunctionf(x).Theorem2.2.Supposethatg(~)beacontinuousfunctionontointerval[a,b]andallderivativesoff:[a,b]→RexistandhaveacommonMsothatmax|f(k)(x)|≤M,forallk.(2.20)x∈[a,b]Furthermore,assumethatGn(x,α)betheTaylorpolynomialofdegreenforf(x)aboutsomeα∈(a,b),sayα=g(~),then∀ε>0andγ∈(a,b)thereexistsN∈Nandinterval(c,d)sothat∀~∈(c,d)andn≥N:|f(γ)−Gn(γ,g(~))|<ε.(2.21)Proof.Letβ∈[a,b]bethepointatwhichwewanttodeterminetheerror.Wesuppose(withoutlossofgenerality)thatβ>α.LetXn(β−x)k(k)s(x)=f(β)−f(x)−f(x),(2.22)k!k=10thens(x)existsforx∈(a,b)andn0(β−x)(n+1)s(x)=−f(x).(2.23)n!
50October28,201311:14WorldScientificReviewVolume-9inx6inAdvances/Chap.242S.AbbasbandyandE.ShivanianNow,considerthefunctionn+1β−xU(x)=s(x)−s(α),(2.24)β−αthenU(α)=U(β)=0.(2.25)n+1β−xFromthedifferentiabilityofs(x)and,itfollowsthatU(x)isβ−αdifferentiableonanysubintervalof(a,b).Now,weapplyTheorem2.1toU(x)ontheinterval[α,β].Thereforethereexistsξβ∈(α,β)sothat0U(ξβ)=0,(2.26)whichimpliesnn(β−ξβ)(n+1)(β−ξβ)−f(ξβ)+(n+1)n+1s(α)=0,(2.27)n!(β−α)then,sinceβ6=ξβ,n+1(β−α)(n+1)s(α)=f(ξβ).(2.28)(n+1)!Now,fromEq.(2.22),wehaven+1(β−α)(n+1)s(α)=f(β)−Gn(β,α)=f(ξβ).(2.29)(n+1)!Sincetheparameterβhasbeenchosenarbitrarythenn+1(x−α)(n+1)∀x∈[α,b],α∈[a,b]:f(x)−Gn(x,α)=f(ξx),(2.30)(n+1)!whereξx∈(α,x).Supposethatγ∈(a,b)andε>0,letα∈(a,γ)thenitisclearthatthereexistsNsuchthatn+1(γ−α)ε∀n≥N:<,(2.31)(n+1)!MthenfromEqs.(2.20)and(2.30),n+1(γ−α)(n+1)ε|f(γ)−Gn(γ,α)|=|f(ξγ)|<·M=ε.(2.32)(n+1)!MTherefore,wehaveproved∀γ∈(a,b),α∈(a,γ),ε>0,∃N⇒∀n≥N:|f(γ)−Gn(γ,α)|<ε.
51October24,201310:44WorldScientificReviewVolume-9inx6inAdvances/Chap.2PredictorHomotopyAnalysisMethod(PHAM)43Sinceg(~)iscontinuousfunctionontointerval[a,b]thenthereexistsin-terval(c,d)suchthatg{(c,d)}=(a,γ).Hencethestatement(2.32)isreadequivalentlyas∀γ∈(a,b),~∈(c,d),ε>0,∃N⇒∀n≥N:|f(γ)−Gn(γ,g(~))|<ε,andtheproofiscompleted.Corollary2.1.Supposethatf(x)issufficientlysmoothoninterval[a,b],g(x)beacontinuousfunctionontointerval[a,b]andGn(x,g(~))betheTaylorpolynomialofdegreenaboutg(~).Then,asngoesinfinity,intheplotofGn(γ,g(~)),a<γ