Computer Methods in Applied Mechanics and Engineering

一些ME专业提升的论文。

ThisarticleappearedinajournalpublishedbyElsevier.Theattachedcopyisfurnishedtotheauthorforinternalnon-commercialresearchandeducationuse,includingforinstructionattheauthorsinstitution

andsharingwithcolleagues.Otheruses,includingreproductionanddistribution,orsellingorlicensingcopies,orpostingtopersonal,institutionalorthirdparty

websitesareprohibited.Inmostcasesauthorsarepermittedtoposttheirversionofthearticle(e.g.inWordorTexform)totheirpersonalwebsiteorinstitutionalrepository.AuthorsrequiringfurtherinformationregardingElsevier’sarchivingandmanuscriptpoliciesare

encouragedtovisit:

http://doc.guandang.net/copyright

一些ME专业提升的论文。

Author's personal copy

Available online at http://doc.guandang.net

Comput.MethodsAppl.Mech.Engrg.197(2008)

2131–2146

http://doc.guandang.net/locate/cma

Afastimmersedboundarymethodusinganullspaceapproach

andmulti-domainfar- eldboundaryconditions

TimColonius*,KunihikoTaira

DivisionofEngineeringandAppliedScience,CaliforniaInstituteofTechnology,CA91125,USAReceived21March2007;receivedinrevisedform3August2007;accepted6August2007

Availableonline12September2007

Abstract

Wereportonthecontinueddevelopmentofaprojectionapproachforimplementingtheimmersedboundarymethodforincompress-ible owsintwoandthreedimensions.BoundaryforcesandpressureareregardedasLagrangemultipliersthatenabletheno-slipandpergence-freeconstraintstobeimplicitlydeterminedtoarbitraryprecisionwithnoassociatedtime-steprestrictions.Inordertoaccel-eratethemethod,wefurtherimplementanullspace(discretestreamfunction)methodthatallowsthepergence-freeconstrainttobeautomaticallysatis edtomachineroundo .Byemployingafastsinetransformtechnique,thelinearsystemtodeterminetheforcescanbesolvede cientlywithdirectoriterativetechniques.Amulti-domaintechniqueisdevelopedinordertoimprovefar- eldboundaryconditionsthatarecompatiblewiththefastsinetransformandaccountfortheextensivepotential owinducedbythebodyaswellasvorticitythatadvects/di usestolargedistancefromthebody.Themulti-domainandfasttechniquesarevalidatedbycomparingtotheexactsolutionsforthepotential owinducedbystationaryandpropagatingOseenvorticesandbyanimpulsively-startedcircularcyl-inder.Speed-upsofmorethananorder-of-magnitudeareachievedwiththenewmethod.Ó2007ElsevierB.V.Allrightsreserved.

1991MSC:76D05;76M12PACS:47.11.+j

Keywords:Immersedboundarymethod;Fractionalstepmethod;Projectionmethod;Nullspacemethod;Vorticity/streamfunctionformulation;Far- eldboundaryconditions;Multi-domainmethod;FastPoissonsolver;Finitevolumemethod;Incompressibleviscous ow

1.Introduction

Intheimmersedboundarymethod(IBmethod),immersedsurfacesaregeneratedbyforcesatasetofLagrangianpoints[29,20,19].The owissolvedonanEuleriangridthatdoesnotconformtothebodygeometry–typicallyauniformCartesiangridisused.Theboundaryforcesthatexistassingularfunctionsalongthesurfaceinthecontinuousequationsaredescribedbydiscretedeltafunctionsthatsmear(regularize)theforcinge ectovertheneighboringEuleriancells.

Correspondingauthor.

E-mailaddresses:colonius@caltech.edu(T.Colonius),kunihiko@cal-tech.edu(K.Taira).

0045-7825/$-seefrontmatterÓ2007ElsevierB.V.Allrightsreserved.doi:10.1016/j.cma.2007.08.014

*

IntheoriginalIBmethod,surfaceswereviewedas ex-ibleelasticmembraneswithaconstitutiverelation(e.g.Hooke’slaw)relatingtheforcestothemotionoftheLagrangianpoints[28].Thistechniquewaslaterextendedtosurfaceswithprescribedmotion(andinparticularrigidbodies)bytakingthespringconstanttobelarge[2,16].Goldsteinetal.[9]appliedtheconceptoffeedbackcontroltocomputetheforceontherigidimmersedsurface.Thedi erencebetweenthevelocitysolutionandtheboundaryvelocityisusedinaproportional-integralcontroller.Con-stitutivelawsareeliminatedinthedirectforcingmethod[21,7];forcinginthemomentumequationisdeterminedbypenalizingtheslipatthe(interpolated)surface.Fortheaforementionedtechniquesthatutilizeconstitutiverela-tions,thechoiceofgain(sti ness)http://doc.guandang.netrgegain

一些ME专业提升的论文。

2132T.Colonius,K.Taira/Comput.MethodsAppl.Mech.Engrg.197(2008)2131–2146

resultsinrestrictionsonthetimestep,whilesmallgainresultsinsliperror.1Directforcingmethodssimilarlyresultinasliperroratthesurface.Whilethesliperrorisreportedtobesmall[7],themagnitudecannotbeestimatedinadeductivemanner.FurtherinformationregardingtheIBmethodandhigher-orderextensionsaregiveninarecentreview[20].

AnalternativeistoregardtheboundaryforcesasLagrangemultiplierswhosevaluesarechosentosatisfytheno-slipconstraint[8,36].Byintroducingappropriateregularizationandinterpolationoperatorsandgroupingthepressureandforceunknownstogether,thediscretizedincompressibleNavier–Stokesequationscanbeformulatedwithastructurealgebraicallyidenticaltothetraditionalfractionalstepmethod[36].Thepressureandforceunknownsarefoundbysolvinga(modi ed)Poissonequa-tion.Inwhatfollows,werefertothismethodastheimmersedboundaryprojectionmethod(IBPM).

TheprincipleadvantagesoftheIBPMtechniquearethatthecontinuityandno-slipconstraintscanbesatis ed(toarbitraryaccuracy)implicitlyatthenexttimelevel,andthattheCourantnumberisonlylimitedbythechoiceoftimemarchingschemesfortheviscousandadvectiontermsinthemomentumequation.Further,itispossibletoarrangealloperationssothatthemethodisuniformlysec-ond-orderaccurateintime,andsothatthematrixarisingfromimplicittreatmentoftheviscoustermsinthemomen-tumequationaswellasthemodi edPoissonmatrixarebothsymmetricandpositivede nite.Consequentlytheconjugate-gradientmethodcanbeusedtosolvethelinearsystems.However,iterativesolutionofthelinearsystemsresultsinaconvergenceerror.Thispresentsnodi cultyinthemomentumequationwherethesolutionneedonlybeconvergedtotheextentthatitissmallerthanotherdis-cretizationerrors.Butinthemodi edPoissonequation,convergenceerrorsdirectlyimpacttheaccuracytowithwhichthepergence-freeandno-slipconstraintsaresatis- ed.Whiletheerrorscanbemadearbi …… 此处隐藏：55722字，全部文档内容请下载后查看。喜欢就下载吧 ……