on a failsafe flux limiting approach for the euler equations
TRANSCRIPT
![Page 1: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/1.jpg)
On a failsafe flux limiting approach for the Euler equations
Dmitri Kuzmin1 Matthias Möller2
1Chair of Applied Mathematics III, University Erlangen-Nuremberg, Germany
2Institute of Applied Mathematics (LS3), TU Dortmund, Germany
![Page 2: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/2.jpg)
Motivation
High-resolution schemes yield accurate results. But . . .
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0
0.2
0.4
0.6
0.8
1
Objective: to design flux corrected FEM with failsafe feature
![Page 3: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/3.jpg)
Outline
1 High-resolution schemeFinite element approximationFlux-correction algorithmFailsafe post-processing
2 ApplicationsConstrained initializationIdealized Z-pinch implosion modelIdeal MHD equations
3 Conclusions
![Page 4: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/4.jpg)
Finite element approximation
∂U
∂t+∇ · F(U) = 0, U =
ρρvρE
, F =
ρvρv ⊗ v + pI(ρE + p)v
Weak formulation
(using integration by parts)
∫Ω
W
[∂U
∂t+∇ · F(U)
]dx = 0, ∀W ∈ W
Group representation1 Uj(t) = U(xj , t), Fj(t) = F(Uj(t))
U(x, t) ≈∑
jϕj(x)Uj(t), F(U) ≈
∑jϕj(x)Fj(t)
1C.A.J. Fletcher, CMAME 1983, 37(2), pp. 225–244
![Page 5: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/5.jpg)
Finite element approximation
∂U
∂t+∇ · F(U) = 0, U =
ρρvρE
, F =
ρvρv ⊗ v + pI(ρE + p)v
Weak formulation (using integration by parts)∫
Ω
W∂U
∂tdx =
∫Ω
∇W · F(U) dx−∫
Γ
W n · F(U) ds, ∀W ∈ W
Group representation1 Uj(t) = U(xj , t), Fj(t) = F(Uj(t))
U(x, t) ≈∑
jϕj(x)Uj(t), F(U) ≈
∑jϕj(x)Fj(t)
1C.A.J. Fletcher, CMAME 1983, 37(2), pp. 225–244
![Page 6: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/6.jpg)
Finite element approximation
∂U
∂t+∇ · F(U) = 0, U =
ρρvρE
, F =
ρvρv ⊗ v + pI(ρE + p)v
Weak formulation (using integration by parts)∫
Ω
W∂U
∂tdx =
∫Ω
∇W · F(U) dx−∫
Γ
W n · F(U) ds, ∀W ∈ W
Group representation1 Uj(t) = U(xj , t), Fj(t) = F(Uj(t))
U(x, t) ≈∑
jϕj(x)Uj(t), F(U) ≈
∑jϕj(x)Fj(t)
1C.A.J. Fletcher, CMAME 1983, 37(2), pp. 225–244
![Page 7: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/7.jpg)
Finite element approximation, cont’d
Semi-discrete high-order scheme∑jmij
dUjdt
=∑
j(cji − sij) · Fj
=: RHi
∀i
mij =
∫Ω
ϕiϕj dx, cji =
∫Ω
∇ϕiϕj dx, sij =
∫Γ
ϕiϕjn ds
Semi-discrete low-order scheme
midUidt
=∑j
(cji − sij) · Fj +∑j 6=i
Dij(Uj − Ui)
=: RLi
∀i
mi =∑
jmij =
∫Ω
ϕi dx, Dij artificial viscosity
![Page 8: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/8.jpg)
Finite element approximation, cont’d
Semi-discrete high-order scheme∑jmij
dUjdt
=∑
j(cji − sij) · Fj
=: RHi
∀i
mij =
∫Ω
ϕiϕj dx, cji =
∫Ω
∇ϕiϕj dx, sij =
∫Γ
ϕiϕjn ds
Semi-discrete low-order scheme
midUidt
=∑j
(cji − sij) · Fj +∑j 6=i
Dij(Uj − Ui)
=: RLi
∀i
mi =∑
jmij =
∫Ω
ϕi dx, Dij artificial viscosity
![Page 9: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/9.jpg)
Finite element approximation, cont’d
Semi-discrete high-order scheme∑jmij
dUjdt
=∑
j(cji − sij) · Fj =: RHi ∀i
mij =
∫Ω
ϕiϕj dx, cji =
∫Ω
∇ϕiϕj dx, sij =
∫Γ
ϕiϕjn ds
Semi-discrete low-order scheme
midUidt
=∑j
(cji − sij) · Fj +∑j 6=i
Dij(Uj − Ui) =: RLi ∀i
mi =∑
jmij =
∫Ω
ϕi dx, Dij artificial viscosity
![Page 10: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/10.jpg)
Finite element approximation, cont’d
Relation between high- and low-order schemes∑j
mijdUjdt
= RHi ⇔ midUidt
= RLi +∑j 6=i
Fij
Raw antidiffusive flux
Fij = mij
(dUidt− dUj
dt
)+Dij(Ui − Uj), Fji = −Fij
Objective: linearize the raw antidiffusive fluxes and limit them toprevent the generation of nonphysical under-/overshoots
![Page 11: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/11.jpg)
Linearized FCT algorithm2
1 Compute the low-order solution at tn+1 = tn + ∆t
miULi − Uni
∆t= θRLi (UL) + (1− θ)RLi (Un), 0 < θ ≤ 1
2 Approximate the time derivative
midUidt
= RLi ⇒ dUidt≈ ULi =
1
miRLi (UL)
3 Linearize the raw antidiffusive fluxes
FLij = mij(ULi − ULj ) +Dij(U
Li − ULj ), FLji = −FLij
4 Apply the limited antidiffusive fluxes
miUn+1i = miU
Li + ∆t
∑j 6=i
αijFLij , 0 ≤ αij = αji ≤ 1
2D. Kuzmin, JCP 2009, 228(7), pp. 2517–2534
![Page 12: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/12.jpg)
Linearized FCT algorithm2
1 Compute the low-order solution at tn+1 = tn + ∆t
miULi − Uni
∆t= θRLi (UL) + (1− θ)RLi (Un), 0 < θ ≤ 1
2 Approximate the time derivative
midUidt
= RLi ⇒ dUidt≈ ULi =
1
miRLi (UL)
3 Linearize the raw antidiffusive fluxes
FLij = mij(ULi − ULj ) +Dij(U
Li − ULj ), FLji = −FLij
4 Apply the limited antidiffusive fluxes
miUn+1i = miU
Li + ∆t
∑j 6=i
αijFLij , 0 ≤ αij = αji ≤ 1
2D. Kuzmin, JCP 2009, 228(7), pp. 2517–2534
![Page 13: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/13.jpg)
Linearized FCT algorithm2
1 Compute the low-order solution at tn+1 = tn + ∆t
miULi − Uni
∆t= θRLi (UL) + (1− θ)RLi (Un), 0 < θ ≤ 1
2 Approximate the time derivative
midUidt
= RLi ⇒ dUidt≈ ULi =
1
miRLi (UL)
3 Linearize the raw antidiffusive fluxes
FLij = mij(ULi − ULj ) +Dij(U
Li − ULj ), FLji = −FLij
4 Apply the limited antidiffusive fluxes
miUn+1i = miU
Li + ∆t
∑j 6=i
αijFLij , 0 ≤ αij = αji ≤ 1
2D. Kuzmin, JCP 2009, 228(7), pp. 2517–2534
![Page 14: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/14.jpg)
Linearized FCT algorithm2
1 Compute the low-order solution at tn+1 = tn + ∆t
miULi − Uni
∆t= θRLi (UL) + (1− θ)RLi (Un), 0 < θ ≤ 1
2 Approximate the time derivative
midUidt
= RLi ⇒ dUidt≈ ULi =
1
miRLi (UL)
3 Linearize the raw antidiffusive fluxes
FLij = mij(ULi − ULj ) +Dij(U
Li − ULj ), FLji = −FLij
4 Apply the limited antidiffusive fluxes
miUn+1i = miU
Li + ∆t
∑j 6=i
αijFLij , 0 ≤ αij = αji ≤ 1
2D. Kuzmin, JCP 2009, 228(7), pp. 2517–2534
![Page 15: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/15.jpg)
Flux limiting for scalar equations
miun+1i = miu
Li + ∆t
∑j 6=i
αijfLij , fLji = −fLij , αij = αji
Zalesak’s limiter3 yields αij ’s such that the nodal values of the corrected
solution are bounded by the local extrema of the low-order solution
xi−1 xi xi+1
umini
uL umaxi
un+1i = uLi +
∆tmi
∑j 6=iαijf
Lij
3S. Zalesak, JCP 1979, 31(3), pp. 335–362
![Page 16: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/16.jpg)
Flux limiting for systems4
Apply flux limiter to a set of control variables simultaneously
Nodal transformation
fρij = Tρi F
Lij
fpij = Tpi F
Lij
fvij = T vi FLij
fρij
FCT
fpij
FCT
αij := minαρij, αpij, α
vij
FCT
f vij
Apply correction factors to the conservative antidiffusive fluxes
miUn+1i = miU
Li + ∆t
∑j 6=i
αijFLij , FLji = −FLij , αij = αji
4D. Kuzmin, M. M, J.N. Shadid, M. Shashkov, JCP 2010, 229(23), pp. 8766–8779
![Page 17: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/17.jpg)
Flux limiting for systems4, cont’d
Apply flux limiter to a set of control variables one after the other
Nodal transformation
fρij = Tρi F
Lij
fpij = Tpi (α
ρijF
Lij )
fvij = T vi (αpijαρijF
Lij )
FCTFCT FCT
fρij
αρij
fpij f vij
αpijαρij αvijα
pijα
ρij =: αij
Apply correction factors to the conservative antidiffusive fluxes
miUn+1i = miU
Li + ∆t
∑j 6=i
αijFLij , FLji = −FLij , αij = αji
4D. Kuzmin, M. M, J.N. Shadid, M. Shashkov, JCP 2010, 229(23), pp. 8766–8779
![Page 18: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/18.jpg)
Failsafe flux correction algorithm4
Flux correction machinery may fail to produce physically correct solutions:
∃i : umini ≤ uFCTi ≤ umax
i is violated for control variable u
Remedy: enforce physically-motivated constraints by post-processing
miUFCTi = miU
Li + ∆t
∑j 6=i
αijFLij
⇔ miUiL = miU
FCTi −∆t
∑j 6=i
αijFLij
Example: β(0)ij ≡ 0, β
(k)ij :=
k/K if failure is detected at i, jβ
(k−1)ij otherwise
4D. Kuzmin, M. M, J.N. Shadid, M. Shashkov, JCP 2010, 229(23), pp. 8766–8779
![Page 19: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/19.jpg)
Failsafe flux correction algorithm4
Flux correction machinery may fail to produce physically correct solutions:
∃i : umini ≤ uFCTi ≤ umax
i is violated for control variable u
Remedy: enforce physically-motivated constraints by post-processing
miUFCTi = miU
Li + ∆t
∑j 6=i
αijFLij
⇔
miUi(k) = miU
FCTi −∆t
∑j 6=i
β(k)ij αijF
Lij , k = 1, ...,K
Example: β(0)ij ≡ 0, β
(k)ij :=
k/K if failure is detected at i, jβ
(k−1)ij otherwise
4D. Kuzmin, M. M, J.N. Shadid, M. Shashkov, JCP 2010, 229(23), pp. 8766–8779
![Page 20: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/20.jpg)
Failsafe flux correction algorithm4
Flux correction machinery may fail to produce physically correct solutions:
∃i : umini ≤ uFCTi ≤ umax
i is violated for control variable u
Remedy: enforce physically-motivated constraints by post-processing
miUFCTi = miU
Li + ∆t
∑j 6=i
αijFLij
⇔
miUi(k) = miU
FCTi −∆t
∑j 6=i
β(k)ij αijF
Lij , k = 1, ...,K
Example: β(0)ij ≡ 0, β
(k)ij :=
k/K if failure is detected at i, jβ
(k−1)ij otherwise
4D. Kuzmin, M. M, J.N. Shadid, M. Shashkov, JCP 2010, 229(23), pp. 8766–8779
![Page 21: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/21.jpg)
Double Mach reflection5
Test: Roe-linearization + FCT, structured grid, Q1 finite elements
T = 0.2, Crank Nicolson time stepping (θ = 0.5), ∆t = 64h · 10−4
60
Γwall
ΓoutUL =
8.0
8.25 cos 30
−8.25 sin 30
116.5
ΓL ΓR
UR =
1.40.00.01.0
5P.R. Woodward, P. Colella, JCP 54, 115 (1984), pp. 115–173
![Page 22: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/22.jpg)
Double Mach reflection
Low-order, h = 1/64, ∆t = 10−4
FCT, αij(ρ, p), h = 1/64, ∆t = 10−4
![Page 23: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/23.jpg)
Double Mach reflection
Low-order, h = 1/128, ∆t = 5 · 10−5
FCT, αij(ρ, p), h = 1/128, ∆t = 5 · 10−5
![Page 24: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/24.jpg)
Double Mach reflection
Low-order, h = 1/256, ∆t = 2.5 · 10−5
FCT, αij(ρ, p), h = 1/256, ∆t = 2.5 · 10−5
![Page 25: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/25.jpg)
Double Mach reflection
Low-order, h = 1/512, ∆t = 1.25 · 10−5
FCT, αij(ρ, p), h = 1/512, ∆t = 1.25 · 10−5
![Page 26: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/26.jpg)
Double Mach reflection
FCT, αij(ρ, ρE), h = 1/512, ∆t = 1.25 · 10−5
FCT, αij(ρ, p), h = 1/512, ∆t = 1.25 · 10−5
![Page 27: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/27.jpg)
Sod’s shock tube problem
Test: Rusanov-type dissipation + FCT, αij(ρ, p), 10n× n grid, Q1 FEs
κu = log‖u2h − u4h‖1‖uh − u2h‖1
/ log 2
Time: 0.231
Crank Nicolson time steppingFCT Low-order
nfine κρ κp κρ κp20 0.624 1.027 0.193 0.62340 0.970 1.003 0.421 0.67180 1.079 1.005 0.575 0.701160 1.073 1.005 0.624 0.730
Backward Euler time steppingFCT Low-order
nfine κρ κp κρ κp20 0.671 0.982 0.190 0.61940 0.980 0.950 0.416 0.66980 0.977 0.947 0.575 0.701160 0.981 0.945 0.624 0.730
![Page 28: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/28.jpg)
Sod’s shock tube problem
Coarse mesh with contour plot of density variable at time T = 0.231
Crank Nicolson time stepping Backward Euler time steppingFCT Low-order FCT Low-order
#trias κρ κp κρ κp κρ κp κρ κp18,176 0.925 0.876 0.364 0.665 0.955 0.841 0.357 0.66272,704 0.874 0.800 0.539 0.679 0.820 0.732 0.536 0.679290,816 0.806 0.934 0.614 0.718 0.765 0.875 0.616 0.719
1,163,264 0.948 0.966 0.641 0.739 0.889 0.905 0.642 0.740
Linearized FCT algorithm yields accurate and non-oscillatory solutions using P1
and Q1 finite elements on structured and unstructured meshes, respectively.
![Page 29: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/29.jpg)
Outline
1 High-resolution schemeFinite element approximationFlux-correction algorithmFailsafe post-processing
2 ApplicationsConstrained initializationIdealized Z-pinch implosion modelIdeal MHD equations
3 Conclusions
![Page 30: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/30.jpg)
Constrained initialization
Test: discontinuous initial data with P1 FEs on unstructured meshΩ = (−0.5, 0.5)2, Ωin = (x, y) ∈ Ω|r =
√x2 + y2 < 0.13
Ωin Ω \ Ωin
ρ0 2.0 1.0v0 0.0 0.0p0 15.0 1.0
Pointwise initialization
Uh(xi, yi) := U0(xi, yi)
is not conservative!
∫Ω
ρ0 dx
∫Ω
(ρE)0 dx
∫Ω
ρh dx
∫Ω
(ρE)h dx
1.05309 4.35825 1.04799 4.17949
![Page 31: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/31.jpg)
Constrained initialization, cont’d
Conservative initialization by L2 projection
Uh =∑j
ϕjUj :
∫Ω
WhUhdx =
∫Ω
WhU0dx ∀Wh ∈ Wh
Consistent mass matrix
∑jmijU
Cj =
∫ΩϕiU0dx
Lumped mass matrix
miULi =
∫ΩϕiU0dx
![Page 32: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/32.jpg)
Constrained initialization, cont’d
Relation between consistent and lumped L2 projection
miUCi = miU
Li +
∑j 6=i
Fij , Fij = mij(UCi − UCj )
Constrained L2 projection
miUi = miULi +
∑j 6=i
αijFij
Apply failsafe flux correctionmachinery to compute
0 ≤ αij = αji ≤ 1
![Page 33: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/33.jpg)
Outline
1 High-resolution schemeFinite element approximationFlux-correction algorithmFailsafe post-processing
2 ApplicationsConstrained initializationIdealized Z-pinch implosion modelIdeal MHD equations
3 Conclusions
![Page 34: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/34.jpg)
Idealized Z-pinch implosion model6
Generalized Euler system coupled with scalar tracer equation
∂
∂t
ρρvρEρλ
+∇ ·
ρv
ρv ⊗ v + pIρEv + pvρλv
=
0f
f · v0
Equation of state
p = (γ − 1)ρ(E − 0.5|v|2)
Non-dimensional Lorentz force
f = (ρλ)(I(t)Imax
)2erreff
, 0 ≤ λ ≤ 1
λ = 1.0ρ = 106
ρ = 1.0λ = 0.0
ρ = 0.5λ = 0.0
v = 0.0, p = 1.0 everywhere
6J.W. Banks, J.N. Shadid, IJNMF 2009, 61(7), pp. 725–751
![Page 35: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/35.jpg)
Idealized Z-pinch implosion
![Page 36: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/36.jpg)
Idealized Z-pinch implosion
![Page 37: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/37.jpg)
Idealized Z-pinch implosion
![Page 38: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/38.jpg)
Idealized Z-pinch implosion
![Page 39: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/39.jpg)
Idealized Z-pinch implosion
![Page 40: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/40.jpg)
Idealized Z-pinch implosion
![Page 41: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/41.jpg)
Idealized Z-pinch implosion
![Page 42: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/42.jpg)
Idealized Z-pinch implosion
![Page 43: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/43.jpg)
Idealized Z-pinch implosion
FCT
Low-order
@@
@@
@@
@@
@I r = 0.5
![Page 44: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/44.jpg)
Outline
1 High-resolution schemeFinite element approximationFlux-correction algorithmFailsafe post-processing
2 ApplicationsConstrained initializationIdealized Z-pinch implosion modelIdeal MHD equations
3 Conclusions
![Page 45: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/45.jpg)
Ideal MHD equations
Idealized MHD equations
∂
∂t
ρρvBρE
+∇ ·
ρv
ρv ⊗ v + pI−B⊗Bv ⊗B−B⊗ v
ρEv + pv−(B · v)B
= 0
subject to ∇ ·B = 0
Divergence involution in 1D: ∂xBx = 0 ⇒ Bx = const
Hyperbolic conservation laws for 7 variables: ρ, v, By, Bz, ρE
Roe matrix (for arbitrary γ) by Cargo and Gallice7
FCT limiter is applied to control variables ρ, p, By and Bz
7P. Cargo, G. Gallice, JCP 1997, 136(2), pp.446–466
![Page 46: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/46.jpg)
Shock tube problem
γ = 1.4, Bx = 0.75, tfin = 0.1, 800 grid points
(ρ,v, By, Bz, p)T =
(1.0 , 0.0, 1.0, 0.0, 1.0)T if x ≤ 0.5
(0.125, 0.0,−1.0, 0.0, 0.1)T if x > 0.5
ρ
p
vx
vy
By
![Page 47: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/47.jpg)
Outline
1 High-resolution schemeFinite element approximationFlux-correction algorithmFailsafe post-processing
2 ApplicationsConstrained initializationIdealized Z-pinch implosion modelIdeal MHD equations
3 Conclusions
![Page 48: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/48.jpg)
Conclusions
Failsafe flux correction algorithm
ensures boundedness of physical quantities
preserves symmetry on unstructured grids
is applicable to ‘challenging’ applications
can be turned into a constrained L2 projection
Future research
extension to multidimensional MHD equations
treatment of the ∇ ·B = 0 involution
consider non-conforming FEs
Thank you for your attention!
![Page 49: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/49.jpg)
Conclusions
Failsafe flux correction algorithm
ensures boundedness of physical quantities
preserves symmetry on unstructured grids
is applicable to ‘challenging’ applications
can be turned into a constrained L2 projection
Future research
extension to multidimensional MHD equations
treatment of the ∇ ·B = 0 involution
consider non-conforming FEs
Thank you for your attention!
![Page 50: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/50.jpg)
References
1 C.A.J. Fletcher, The group finite element formulation. CMAME 1983,37(2), pp. 225–244.
2 D. Kuzmin, Explicit and implicit FEM-FCT algorithms with fluxlinearization. JCP 2009, 228(7), pp. 2517–2534.
3 S.T. Zalesak, Fully multidimensional flux-corrected transport algorithmsfor fluids. JCP 1979, 31(3), pp. 335–362.
4 D. Kuzmin, M. M, J.N. Shadid, M. Shashkov, Failsafe flux limiting andconstrained data projections for equations of gas dynamics. JCP 2010,229(23), pp. 8766–8779.
5 P.R. Woodward, P. Colella, The numerical simulation of two-dimensionalfluid flow with strong shocks. JCP 54, 115 (1984), pp. 115–173
6 J.W. Banks, J.N. Shadid, An Euler system source term that developsprototype Z-pinch implosions intended for the evaluation of shock-hydromethods. IJNMF 2009, 61(7), pp. 725–751.
7 P. Cargo, G. Gallice, Roe matrices for ideal MHD and systematicconstruction of Roe matrices for systems of conservation laws. JCP 1997,136(2), pp.446–466.
![Page 51: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/51.jpg)
Zalesak’s flux limiter5
P+i
P−i
ui
Consider sums of positive/negativeantidiffusive fluxes into each node
Limit antidiffusive flux if it exceedsthe distance to upper/lower bounds
Compute nodal correction factors
R+i = min1, Q+
i /P+i and R−i = min1, Q−i /P
−i
Limit antidiffusive flux for edge ij by
αij =
minR+
i , R−j for positive fluxes
minR−i , R+j for negative fluxes
5S. Zalesak, JCP 1979, 31(3), pp. 335–362
![Page 52: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/52.jpg)
Zalesak’s flux limiter5
Q+i
P+i Q−i
P−i
Consider sums of positive/negativeantidiffusive fluxes into each node
Limit antidiffusive flux if it exceedsthe distance to upper/lower bounds
Compute nodal correction factors
R+i = min1, Q+
i /P+i and R−i = min1, Q−i /P
−i
Limit antidiffusive flux for edge ij by
αij =
minR+
i , R−j for positive fluxes
minR−i , R+j for negative fluxes
5S. Zalesak, JCP 1979, 31(3), pp. 335–362
![Page 53: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/53.jpg)
Zalesak’s flux limiter5
Q+i
P+i Q−i
P−i
Consider sums of positive/negativeantidiffusive fluxes into each node
Limit antidiffusive flux if it exceedsthe distance to upper/lower bounds
Compute nodal correction factors
R+i = min1, Q+
i /P+i and R−i = min1, Q−i /P
−i
Limit antidiffusive flux for edge ij by
αij =
minR+
i , R−j for positive fluxes
minR−i , R+j for negative fluxes
5S. Zalesak, JCP 1979, 31(3), pp. 335–362
![Page 54: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/54.jpg)
Zalesak’s flux limiter5
Q+i
P+i Q−i
P−i
Consider sums of positive/negativeantidiffusive fluxes into each node
Limit antidiffusive flux if it exceedsthe distance to upper/lower bounds
Compute nodal correction factors
R+i = min1, Q+
i /P+i and R−i = min1, Q−i /P
−i
Limit antidiffusive flux for edge ij by
αij =
minR+
i , R−j for positive fluxes
minR−i , R+j for negative fluxes
5S. Zalesak, JCP 1979, 31(3), pp. 335–362
![Page 55: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/55.jpg)
Extended version of Zalesak’s FCT limiter
Input: auxiliary solution uL and antidiffusive fluxes fuij , where fuji 6= fuij
1 Sums of positive/negative antidiffusive fluxes into node i
P+i =
∑j 6=i max0, fuij, P−i =
∑j 6=i min0, fuij
2 Upper/lower bounds based on the local extrema of uL
Q+i = mi(u
maxi − uLi ), Q−i = mi(u
mini − uLi )
3 Correction factors αuij = αuji to satisfy the FCT constraints
αuij = minRij , Rji, Rij =
min1, Q+
i /P+i if fuij ≥ 0
min1, Q−i /P−i if fuij < 0
![Page 56: On a failsafe flux limiting approach for the Euler equations](https://reader030.vdocument.in/reader030/viewer/2022012700/61a367aa28f108792a48e9cb/html5/thumbnails/56.jpg)
Node-based transformation of control variables
Conservative variables: density, momentum, total energy
Ui = [ρi, (ρv)i, (ρE)i] , Fij =[fρij , f
ρvij , f
ρEij
], Fji = −Fij
Primitive variables V = TU : density, velocity, pressure
Vi = [ρi,vi, pi] , vi = (ρv)iρi
, pi = (γ − 1)[(ρE)i − |(ρv)i|2
2ρi
]Gij =
[fρij , f
vij , f
pij
]= T (Ui)Fij , T (Uj)Fji = Gji 6= −Gij
Raw antidiffusive fluxes for the velocity and pressure
fvij =fρvij −vif
ρij
ρi, fpij = (γ − 1)
[|vi|2
2 fρij − vi · fρvij + fρEij
]