High order time space implicit discontinous Galerkin method with reduced evaluation costs

Florent Renac

Description

Niveau: Supérieur
THE FRENCH AEROSPACE LAB Equations and Numerical Approach Numerical experiments Summary and outlook Time implicit high-order discontinuous Galerkin method with reduced evaluation costs Florent Renac 1 , Frederic Coquel 2 & Claude Marmignon 1 1 ONERA/DSNA/NUMF 2 CMAP, Ecole Polytechnique JSO High Fidelity Flow Simulations October 8th, 2010 ONERA Châtillon Implicit DG with reduced evaluation costs 1/30

• jso high


• low order


• implicit procedure


• nonlinear di?usion


• high-order discontinuous


• model problem


• numerical experiments


• ?? ·

Sujets

• École polytechnique

Informations

• Publié par: chaeh
• Nombre de lectures: 25
• Langue: English

Extrait

THE FRENCH AEROSPACE LAB
w1EquationscostsandteNumericalCh?tillonAppClaudero,aHighc8th,[e-mail]cit&hMaigh-o1rderONERA/DSNA/NUMFdiscontinuousCMAPGalerkinEcolemethoolydchniquewithreducedFidelitevaluationFlocostsSimulationsFloerrent201ReONERAnacflorent.1Implicit,withFevaluationrederic1/30CoTHE FRENCH AEROSPACE LAB
steadyecientEquationsDGandoNumericalschemeApp.roDGaofcrabhAntNumericalinexper-StokecedureruncouplingimentsinSummaDOFryofandexpoutlo,ok1997.Motivationcostsandracticalgoal3DDGeq.)metphrodd:locompactDOFstencilelement:owapplicationell-suitedoforrequationsunstructuredwithmeshes,[e-mail],pBCappapplication.(e.g.timeNaviexplicitesdiscretization:Aim:strongestimplicitCroFLforestrictiontheassomethoci:atedoftowparderrabsoliceachtermreconstructiontimehigherimplicitrderdiscretizatisoton:solutihighncompnonlineautatipaoolicnalnumericalcosterimentsiBR2nBassiducedal.b2ndyTC,thewlaen,rflorent.geImplicitnumbwitherevaluationof2/30DOFTHE FRENCH AEROSPACE LAB
ImplicitdiusionEquationscedureandSummaNumericalnApperimentsroraoutlocevaluationhpNumericalNumericalexpnonlineae2DrequationimentsanSummaflorent.rywithand3/30outloImplicitokroOutline21expEquations1DandrNumericalequation[e-mail]sSpacediscretizatioTHE FRENCH AEROSPACE LAB
2DSpaceEquationswithandNumericalNumericalanAppEntropropanonlineacequationhflorent.Numerical4/30expsolutioenrcedureimentserimentsSummadiusionryrandSummaoutlooutlookImplicitNonlineaevaluationrequationdiusionyequationnEntropdiscretizatioyImplicitsorolution2Spaceexpdiscretization1DImplicitrpequatio[e-mail]neacostsrdiusion−∇·( ( , )∇ ) = ( ) Ω⊂R
( , ) ∈R ∂Ω = Γ ∪Γ
= Γ
∇ · = Γ
−∇·( ( , )∇ ) = ( ) Ω×( ,∞)
( , ) = ( ) Ω
THE FRENCH AEROSPACE LAB
fastuNonlinearellipticeinsokexppx0Numericalsolveind:hroblemcConsideraequationdu(1)0withWBwithromaxtApprupNumericaluaroblemnonlineadelrdiusionfunctionandofxuImplicitandxEquationsNandeBC(1)saontimeimentsrchingmetholutionusooutloyNonlineaDBEntropxequationuNe:thSummasflorent.xImplicitinwithmoevaluation0urD(2)[e-mail]5/30reducedDcostsonuU( )
U >
U
( , ) =U ( )
×∈R
( )
U =
( ) −∇·( ( )∇ ) = ( ) Ω×( ,∞)
THE FRENCH AEROSPACE LAB
folloisNumericalanrentropuuypafunctionryfoinr1(2)(inifpsolutionoutlouuCys00entropvsatisesobtainsBlineacedurecxabrofu:pandImplicitimentsdiscretizationxuuhCxSpaceNumericalxNonlinealutionuwhereokConesotheywingdrEntropspade)equationrflorent.olicImplicitroblemwithoevaluationv6/30udenitevmatrixtIntroSummaducert[e-mail]eΩ
[
Ω = Ω
=
V ={ϕ∈ (Ω ): ϕ| ∈P (Ω ), = , }Ω
X
P (Ω )={ϕ∈ (Ω ): ϕ( , )= α φ ( , )}
=
∈V
X
( , )= φ ( ) ( ) ∀ ∈Ω , ∈( ,∞)
=

+ =
= ( + )( + )
=
THE FRENCH AEROSPACE LAB
tiroljfoprtitionImplicitwp:discretizationjjoSpaceplutionvjNsop1spacey1Nk:EntroppwithNequationdphdiusioneakjhrrmNonlineaolynomialsokdiscontinuousoutloofandxLP20ryRemajeSummaDOFimentsdiscretizreeofx1expdiscreteypNumerical:hhNxptcplpal01rollxAppVljNumericaltxunctionandFyjEquationstcedureaSpacejdiscretNirhnumbzationr2fifsLerflorent.aImplicitonwithlement:evaluationp7/30hthe:[e-mail]ationof∀ = ,...,
Z Z I Z
φ ( ) + ( θ )·∇φ − φ = φ ( )
Ω Ω ∂Ω Ω
θ , ∇ +
- +v vh hσ
= {∇ + }·
+Ω
-Ω
X
nσ
, ∀ ∈Ω
σ∈∂Ω

+ −Z Z +
{ }=σ + −
φ =− {φ}( − )
+ −Ω ∪Ω σ
THE FRENCH AEROSPACE LAB
jaRSpacefocedurejjdiscretroNumericalkdiusionhsovhdSdVAppwNumericalqroNjEntropand1kksuEquationsexpxSummacxrcostsaxrliftingBR2opwitheratohrsrpequationhyImplicitrdiscretizationhvxhlutionSpaceandCkhNumericalimentserrxhdrydSkopdtjevaluationvhhhvjhhnvmeanflorent.eratoImplicitrmulation:reduc[e-mail]8/30d+ ( ) =
Z
= φ φ ( )
Ω
( + ) ( ) ( )
( − )+ ( ) =
( ) = ( Δ )
∂
= +
Δ ∂
= ×
O( )
THE FRENCH AEROSPACE LAB
rdMisdmatrix:V.dttSummadiscretizationRsizeimentsyVVrtheep0RemaMasssquamatriNx:rixMDOFklaexpnNumericalandhljrocMkSpaceak:lunsymmetricumatrixvDOFroNvmhisAppNdVxnBackwbackwadiscretizationrdTimeEulecedurerAscheme:impAicitequationAV1NumericaltnImplicitandR1VEquationsrdiusionAVanequation:[e-mail]luationry9/30V+ −ˆ ( , )
≤ ≤
< ≤
= =
+ +
+ +
v vh h
=
=
x xx xx j+ ½ x j+ ½j- ½ j- ½
+ - + -
THE FRENCH AEROSPACE LAB
moefsdiusionrepcedureresentecouplingpbImplicitetSummawNumericaleensimpDOFmplessimpcouplingsoappokea(a)rscthroughWnumericals.t.uxqmatrixefrvandh0:implicitSpacevEntropBasics:Whand(1/2)WvcoureductionNumericalhrocostEquationsLohigherwdeseringpcouplingpbImplicitet1Dwxaeenfomopdes2evaluationp(b)[e-mail]reducedycostsequationprroblemNonlinearesolution:outlomorydesimentss.t.r0fullcedureplingqexprohpasimpApp2andreconstructionWcoofbreducedy:florent.1Implicitlowithwevaluationo10/30rder