A Pseudo-Wavelet Scheme for the Two-Dimensional Na

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

APseudo-WaveletSchemefortheTwo-DimensionalNavier-StokesEquationsPhilippeCharton1;3ValeriePerrier2;3March2,1996AbstractNumericalschemestakingadvantageofthewaveletbasescapabilitiestocom-pressbothfunctionsandoperatorsarepresented.GalerkinmethodsfortheheatandPoissonequationsindimensiontwoaredevelopedandnumericallytested.Thenapseudo-waveletscheme,usingcollocationforthenonlinearterm,isdenedforthetwo-dimensionalNavier-Stokesequations.Numericaltestsproveitsaccuracy.SubmittedtoMatematicaAplicadaeComputacional1IREMIA,UniversitedeLaReunion,15Av.ReneCassin,BP7151,97715SaintDenisMessag.Cedex9,France(Philippe.Charton@univ-reunion.fr)2Laboratoired’Analyse,GeometrieetApplications,UA742InstitutGalilee,UniversiteParisNord,AvJ.B.Clement,93430Villetaneuse,France3LaboratoiredeMeteorologieDynamique,E.N.S.,24rueLhomond,75231Pariscedex05,France(perrier@lmd.ens.fr).11IntroductionDespitealargenumberofworksandcomputationsdonesincecomputationaluiddy-namicswasborn25yearsago,thedirectnumericalsimulationofturbulenceisstillagreatchallenge.Itconcernspracticalpurposessuchasindustrialormeteorologicalcom-putations,aswellasfundamentalstudies.DirectnumericalsimulationofturbulencerequirestheintegrationintimeofthenonlinearNavier-Stokesequations.However,atlargeReynoldsnumber,thatiswhennonlinearinteractionsarefardominentuponviscouseects,turbulentowsgenerateincreasinglysmallscales.Inconsequence,toberealistic,thediscretizationinspace(andcorrelativelyintime)oughttohandleahugenumberofdegreesoffreedom.Indimensiontwo,directnumericalsimulationofhomogeneousturbulentowsintheincompressiblecasecanbeperformeduptoquitelargeReynoldsnumbersbywayofspectralFouriertechniques;howevertheseReynoldsnumbersaretoolowtocomparetolargescaleatmosphericdynamicsforexample.ThesituationismuchmoredicultindimensionthreewhereonlymoderateReynoldsnumbersareathandevenwiththepresentlargestsupercomputers.Manytentativeshavebeendoneorareunderwaytoovercomethisproblem;LargeEddySimulations(LES)andsubgrid-scalestechniquestrytoseparatethelargestscales,thatareexplicitlycomputed,fromthesmallscales,thatareparametrizedorstatisticallyrepresented.Vortexmethodsorcontoursurgeryareableindimensiontwo,togenerateverythinscaleshowevertheycan’tadapttoverygeneralproblems.Anotherpropertyofturbulentowsistheirstrongintermittency.Smallscalesmaybelocalizedonlyinaverysmallfractionofspace.Fromthispropertyonemaydreamonnewbasisfunctionsmoresuitabletorepresentthisintermittentspatialstructurewithonlyafewnumberofdegreesoffreedom.Theproblemishowtoconcentratethedegreesoffreedomattheplaceswheresmallscalesarepresent?Thisasksrsttodetectwheresmallscalesareorwillhappen,andthentorepresentthemwithonlyafewparameters.Inthatcontext,waveletsappearedtobeverypromising.Indeedweknow,fromthemathematicalpropertiesofwaveletbases,thatwhereaeldhasasingularityofagivenorder,itswaveletcoecientsatsmallscalesgrowatarelatedpowerofthescale[26,28].Thismeansthatatagivensmallscalethenumberofdegreesoffreedomneedednolongerdependsonthescaleitselfbutonthenumberofactivesingularitiesoftheeldatthisscale.Thispropertywhichcomesfromthedoublelocalization,inspaceandscale,ofthewavelets,hasbeenusedforanalysisandcompressionofturbulentelds[18].Itisalsothebasisofseveralwaveletmethodsdesignedtosolvepartialdierentialequations(PDEs).SchematicallythewaveletbasedmethodsforPDEscanbeseparatedintothreeclasses.Inarstclassofmethods,waveletsareused,intheframeworkofaclassicalgrid-adaptivenumericalcode,todetectwherethegridhastoberenedorcoarsedtooptimallyrepresentthesolution[37,29].Multiresolutionanalysisandtheirassociatedscalefunctionbasesmaybeusedasal-ternativebasesinGalerkinmethods[30,24,42,41].Suchmethodshavethusconvergencepropertiessimilartotheonesofspectralmethods,andsimultaneouslypartialderivativeoperatorsdiscretizesimilarlyasinnite-dierencemethods.However,asthesemethodsdon’tusewaveletsbutratherscalefunctionsasbasisfunctions,theycan’tbeadaptivemethodsandcan’treducesignicantlythenumberofdegreesoffreedominanumericalcode.Thethirdclass,theonlyonewhichusesthecompressionpropertiesofwaveletbases,2containsspecicwaveletmethodsforPDEs.Intheliterature,manytentativeshavebeenperformed,oftenbasedonGalerkinorPetrov-Galerkinmethods.Someofthemtakead-vantageofthewaveletcompressionofthesolution[32,33,3],othersuseinsteadthewaveletcompressionoftheoperator[5,25],theaimofourstudyistotakebenetsimultaneouslyofthesetwoproperties.Thus,ourmethodtakesplaceinthespecic-waveletGalerkinmethodcategoryforsolvingthe2DNavier-Stokesequations.Thealgorithmisbasedonfastmatrix-vectorproducts.Thenonlineartermistreatedbymeansofcollocation.A2Dschemefortheheatequation,alreadyannouncedin[8],isdeveloped.TherstresultsforthewaveletNavier-Stokescodewerepresentedin[10]andpublishedinChartons’PhDthesis[7].Nearlyatthesametimeawavelet-algorithmforthe2DNavier-Stokesequation,basedonthebiorthogonalapproachandaFFT-Poissonsolver,wastestedbyFrohlichandSchneider[22].Thepaperwillbeorganizedasfollows:Section2describesthechosenwaveletapproachforsolvingtheone-dimensionalheatequation,itiscomparedtoclassicalm

1 / 21
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功