关于particle-method-的个人学习心得.pptVIP

SolvingtheEPDiffEquationsbyParticleMethodbyHuMinganAndHuangJinziAndComparedwiththeresultsofFiniteVolumeMethod

Problem:EPDiffEquationsGeneralForm:1DForm:

Thesystemarisesindiversescientificapplications.Inparticular,itappliesinavarietyoffluidsituationsrangingfromsolitonstoturbulence.TheEPDiffequationcoincideswiththedispersionlesscaseoftheCamassa-Holm(CH)equationforshallowwaterinone-dimension(1-D)andtwo-dimensions(2-D).The1-DEPDiffequation(thatis,thedispersionlessCHequation),iscompletelyintegrableandthecorrespondingsolitarywavesthatemergefromtheinitialdata,exhibitaremarkablestability–theiridentityispreservedthroughnonlinearinteractions.Background:Property:

Method:ParticleMethodEquations:InitialData:

WeakSolution:Suchthat:where

TosolvetheEPDiffEquations,weneedtosolvetheODEs.Inthiscase,weneedtosolvethefollowingODEs:

WeusetheSSPRKmethodtosolvetheODEstohavep(t)andx(t).

Results:Inthenumericalexamples,weseta=1andplotthepicturesofuandx(t)p(t)toseetheirevolution.

N=1,P0=[2],X0=[0]Evolutionofu

Movieofu

ResultofFVmethodEvolutionofu

N=2,P0=[3,1],X0=[0,18]Evolutionofu

EvolutionofxEvolutionofpThebluelineisparticle1,theredlineisparticle2Theblacklineisthetotalmomentum,thesumofp.

Movieofu

EvolutionofuN=3,P0=[3,2,1],X0=[0,4,18]

EvolutionofpEvolutionofxThebluelineisparticle1,theredlineisparticle2andthegreenlineisparticle3.Theblacklineisthetotalmomentum,thesumofp.

Movieofu:

Anti-peakon:N=2,P0=[2,-2],X0=[0,17]Evolutionofu

EvolutionofxEvolutionofp

N=2,P0=[2,-2],X0=[0,17]Anotherpossibleresultunderthesameinitialconditions:Evolutionofu

EvolutionofxEvolutionofp

Movieofu

ComparedwithCUmethod:Inthiscase,wecanseethatPMismoreefficientandhavebetterresults.1.PMkeeptheformofthesolutionwellwhileCUmakeittoos