and discussed in where a spectrally accurate approach is applied to the direct numerical simulation DNS of high Reynolds number flows

Richard Pasquetti

Description

Abstract 1. Introduction (and discussed in [16]) where a spectrally accurate approach is applied to the direct numerical simulation (DNS) of high Reynolds number flows. Here we essentially focus on the spectral vanishing viscosity (SVV) method, which appears to be an e?cient stabilization technique possessing the property to preserve the spectral accuracy. It was initially * High Reynolds number flows are di?cult to compute, especially when using spectrally accurate nu- merical schemes. This directly results from the fact that spectral approximations are much less numerically di?usive than low-order ones, so that the non-artificially dissipated energy accumulates at the high spatial frequencies and finally leads to the divergence of the computations. One way to overcome this di?culty is to use stabilization techniques, but then the spectral accuracy of the algorithm is generally destroyed. This is, e.g., particularly obvious for approaches which essentially consist in adding some O?hr? hyper-viscous term, in the spirit of [14]. For a long time filtering techniques have also been proposed to overcome the stability problem. In the frame of spectral element approximations it is however essential to preserve the inter-el- ement continuity, as discussed in [2]. One of the most recent advances in this field has been proposed in [6] The possibility of using the spectral vanishing viscosity method for the spectral element computation of high Reynolds number incompressible flows is investigated.

• i?0


• viscosity method


• stabilized spectral


• vn ?


• svv method


• q?oxun ? goxg


• ?ox


• q?run ?

Sujets

• Antipolis
• Xiamen University

Informations

• Publié par: pefav
• Nombre de lectures: 13
• Langue: English

Extrait

JournalofComputationalPhysics196(2004)680–704 www.elsevier.com/locate/jcp Stabilizedspectralelementcomputationsofhigh Reynoldsnumberincompressibleflows ChuanjuXu a ,RichardPasquetti b, * a DepartmentofMathematics,XiamenUniversity,Xiamen,PRChina b Lab.J.A.Dieudonn
e,UMRCNRS6621,Universite
deNice-SophiaAntipolis,Nice,France Received16June2003;receivedinrevisedform14November2003;accepted15November2003 Abstract Thepossibilityofusingthespectralvanishingviscositymethodforthespectralelementcomputationofhigh Reynoldsnumberincompressibleflowsisinvestigated.Anexponentiallyaccuratestabilizedformulationisproposed andthenappliedtothecomputationofthe2Dwakeofacylinder.Suchaformulationcanbeeasilyimplementedin existingspectralelementsolvers,sinceonlymodifyingthecomputationoftheviscoustermwhilepreservingthe symmetryofthecorrespondingbilinearform.
2003ElsevierInc.Allrightsreserved. 1.Introduction HighReynoldsnumberflowsaredifficulttocompute,especiallywhenusingspectrallyaccuratenu- mericalschemes.Thisdirectlyresultsfromthefactthatspectralapproximationsaremuchlessnumerically diffusivethanlow-orderones,sothatthenon-artificiallydissipatedenergyaccumulatesatthehighspatial frequenciesandfinallyleadstothedivergenceofthecomputations.Onewaytoovercomethisdifficultyisto usestabilizationtechniques,butthenthespectralaccuracyofthealgorithmisgenerallydestroyed.Thisis, e.g.,particularlyobviousforapproacheswhichessentiallyconsistinaddingsomeOð h r Þhyper-viscousterm, inthespiritof[14].Foralongtimefilteringtechniqueshavealsobeenproposedtoovercomethestability problem.Intheframeofspectralelementapproximationsitishoweveressentialtopreservetheinter-el- ementcontinuity,asdiscussedin[2].Oneofthemostrecentadvancesinthisfieldhasbeenproposedin[6] (anddiscussedin[16])whereaspectrallyaccurateapproachisappliedtothedirectnumericalsimulation (DNS)ofhighReynoldsnumberflows. Hereweessentiallyfocusonthespectralvanishingviscosity(SVV)method,whichappearstobean efficientstabilizationtechniquepossessingthepropertytopreservethespectralaccuracy.Itwasinitially * Correspondingauthor. E-mailaddress:[e-mail](R.Pasquetti). 0021-9991/$-seefrontmatter
2003ElsevierInc.Allrightsreserved. doi:10.1016/j.jcp.2003.11.009

C.Xu,R.Pasquetti/JournalofComputationalPhysics196(2004)680–704681 developedfortheresolutionofhyperbolicequationsusingstandardFourierspectralmethods[21].The non-periodiccasewasthenconsideredintheframeofthespectralLegendreapproximation[13].Further refinementshavebeenrecentlysuggested,throughtheuseofaspectralhyper-viscous(ratherthanaviscous) termorthroughtheredefinitionofthestabilizationterm[8].Recently,ithasalsobeensuggestedtousethe SVVmethod,initsfirstformulation,forthelarge-eddysimulationofturbulentflows[10,15]. Inthispaper,ourgoalistocheckthecapabilitiesoftheSVVmethod,intermsofaccuracyandstability, whenitisimplementedinaNavier–Stokesspectralelementsolver. FirstweshowhowtoimplementtheSVVmethodintheframeofaspectralelementapproximation.The factthatcomplexmultidimensionalgeometriesandvectorvaluedfunctionsareconcernedmakethispoint non-trivial,sothatonecannotyetconsiderthatastandardwaytoimplementtheSVVmethodalready exists.Moreover,wesuggestusinganapproximateformwhichcanbeefficientlyimplementedinany spectralelementsolver.Theadvantageofsuchanapproximateformisthatthecomputationalcostper iteration(time-stepifadirectsolverisused)isroughlythesamewithandwithoutSVVstabilization. Secondweconsideranellipticequationsolvedbyasteepanalyticalsolutionandshowthatthecon- vergenceresultsobtainedwiththespectralelementapproximationarecoherentwiththoseobtainedinthe 1Dperiodiccase,whenusingFourierexpansions.AdetailedstudyoftheinfluenceoftheSVVtuning parametersontheconvergenceratesisprovided.Thenweconsidertheso-called‘‘Kovasznayflow’’,which isanexactsolutionoftheincompressibleNavier–Stokesequations,andagaincheckthecapabilitiesofthe SVVmethod. Third,inordertonumericallydemonstratethestabilizationpropertyofthemethod,wecompute2D wakesofacylinder,atReynoldsnumbersuptoRe ¼ 1000,i.e.,muchhigherthanthecriticalvalueas- sociatedwiththe2D–3Dtransition(Re
190). Finally,weconcludebyemphasizingtheinterestintheSVV-stabilizedspectralelementmethod(SEM) forthelarge-eddysimulation(LES)ofturbulentflows. 2.Stabilizedspectralelementformulation TheflowofanincompressibleNewtonianfluidisgovernedbythe‘‘incompressibleNavier–Stokes equations’’.ForanunsteadyflowinadomainXtheyread:
D t u  m r 2 u þr p ¼ sinX  R þ ; r u ¼ 0inX  R þ ; ð1Þ whereu, p andsdenotethevelocity,pressureandsourceterm,respectively, D t uthematerial(Lagrangian) derivativeofuwithrespecttotimetandmthedimensionlessviscosity(theinverseoftheReynoldsnumber). TheunsteadyNavier–Stokesequationsmustbeassociatedtoappropriateinitialandboundaryconditions, e.g.u ð t ¼ 0 Þ¼ 0(fluidatrest)andu j C ¼ 0(no-slipconditionattheboundaryCofX),inordertosetupa well-posedproblemthatonecanthentrytosolvenumerically. If(i)complexgeometriesareconsideredand(ii)highaccuracyisdesired,thentheSEMiswellsuited (see,e.g.[11]).ThespectralelementapproximationoftheweakformoftheincompressibleNavier–Stokes equationsyieldsthefollowingsemi-discretevariationalproblem,tobesolvedateachtime-stepafterthe time-discretization:Find u N 2 X N and p N 2 M N suchthat
ð D t u N ; v N Þþ m ðr u N ; r v N Þðr v N ;p N Þ¼ð s N ; v N Þ8 v N 2 X N ; ð 2 Þ ðr u N ;q N Þ¼ 0 8 q N 2 M N ; whereu N ,p N ands N denotethespectralelementapproximationsofu,pandsandwhere ð ; Þ isusedto denotethestandardL 2 ð X Þ innerproduct,withoutdifferenceifscalar,vectorialortensorialfunctionsare