From Finite Volumes to Discontinuous Galerkin and Flux Reconstruction
|
|
- Philomena Jacobs
- 6 years ago
- Views:
Transcription
1 From Finite Volumes to Discontinuous Galerkin and Flux Reconstruction Sigrun Ortleb Department of Mathematics and Natural Sciences, University of Kassel GOFUN 2017 OpenFOAM user meeting March 21, 2017
2 Numerical simulation of fluid flow This includes flows of liquids and gases such as flow of air or flow of water src: NCTR src: NOAA Spreading of Tsunami waves Weather prediction src: NASA Goddard Flow through sea gates Requirements on numerical solvers CC BY 3.0 Flow around airplanes High accuracy of computation Detailed resolution of physical phenomena Stability and efficiency, robustness Compliance with physical laws (e.g. conservation) 1
3 Contents 1 The Finite Volume Method 2 The Discontinuous Galerkin Scheme 3 SBP Operators & Flux Reconstruction 4 Current High Performance DG / FR Schemes
4 Derivation of fluid equations Based on physical principles: conservation of quantities & balance of forces mathematical tools: Reynolds transport & Gauß divergence theorem Different formulations: Integral conservation law d u dx + F(u, u) n dσ = dt V V V s(u, x, t) dx Partial differential equation u t + F = s embodies the physical principles V n 2
5 Derivation of the continuity equation Based on Reynolds transport theorem d u(x, t) dx = dt V t rate of change in moving volume = u(x, t) V t t rate of change in fixed volume dx + u(x, t) v n dσ V t + convective transfer through surface (V t control volume of fluid particles) 3
6 Derivation of the continuity equation Based on Reynolds transport theorem d u(x, t) dx = dt V t rate of change in moving volume = u(x, t) V t t rate of change in fixed volume Physical principle: conservation of mass dx + u(x, t) v n dσ V t + convective transfer through surface (V t control volume of fluid particles) dm dt = d ρ ρ dx = dt V t V t t dx + ρv n dσ = 0 V t 3
7 Derivation of the continuity equation Based on Reynolds transport theorem d u(x, t) dx = dt V t rate of change in moving volume = u(x, t) V t t rate of change in fixed volume Physical principle: conservation of mass dx + u(x, t) v n dσ V t + convective transfer through surface (V t control volume of fluid particles) dm dt = d ρ ρ dx = dt V t V t t dx + ρv n dσ = 0 V t Divergence theorem yields [ ] ρ t + (ρv) dx = 0 ρ t + (ρv) = 0 V V t continuity equation 3
8 The compressible Navier-Stokes equations... are based on conservation of mass, momentum and energy u t + F = s [ ] u t + Finv + (A(u) u) = s inviscid & viscous fluxes Conservative variables u R 5, fluxes F R 3 5 and sources s R 5 ρ ρv 0 u = ρv, F = ρe ρv v + pi τ (ρe + p)v κ T τ v, s = ρg ρ(q + g v) convective fluxes, heat fluxes & surface forces p pressure, T temperature τ viscous stress tensor κ thermal conductivity body forces & heat sources g gravitational & electromag. forces q intern. heat sources simplified programming by representation in same generic form sufficient to develop discretization schemes for generic conservation law 4
9 General discretization techniques Finite differences / differential form approximation of nodal values and nodal derivatives easy to derive, efficient essentially limited to structured meshes Finite volumes / integral form approximation of cell means and integrals conservative by construction suitable for arbitrary meshes difficult to extend to higher order Finite elements / weak form weighted residual formulation quite flexible and general suitable for arbitrary meshes x i,j+1 x i 1,j x i,j x i+1,j x i,j 1 5
10 Finite volume schemes... based on the integral rather than the differential form Integral conservation enforced for small control volumes V i defined by computational mesh V = K i=1 V i Degrees of freedom: cell means u i (t) = 1 u(x, t) dx V i V i To be specified: concrete definition of control volumes type of approximation inside these cell-centered vs. vertex-centered possibly staggered for different variables numerical method for evaluation of integrals and fluxes 6
11 Why the integral form? Because this is the form directly obtained from physics. 1D scalar hyperbolic conservation law u (u) + f = 0 t x f (u) = 1 2 u2 u + u u = 0 t x (Burgers equation) 1, x < 0, init. cond.: u 0(x) = cos(πx), 0 x 1, 1, x > 1. PDE theory tells us: As long as the exact solution is smooth, it is constant along charactereristic curves 7
12 Why the integral form? Because this is the form directly obtained from physics. 1D scalar hyperbolic conservation law u (u) + f = 0 t x f (u) = 1 2 u2 u + u u = 0 t x (Burgers equation) 1, x < 0, init. cond.: u 0(x) = cos(πx), 0 x 1, 1, x > 1. PDE theory tells us: As long as the exact solution is smooth, it is constant along charactereristic curves Characteristic curves are straight lines and cross smooth solution breaks down integral form (time integrated) still holds 7
13 Why the integral form? Because this is the form directly obtained from physics. 1D scalar hyperbolic conservation law u (u) + f = 0 t x f (u) = 1 2 u2 u + u u = 0 t x (Burgers equation) 1, x < 0, init. cond.: u 0(x) = cos(πx), 0 x 1, 1, x > 1. PDE theory tells us: As long as the exact solution is smooth, it is constant along charactereristic curves Characteristic curves are straight lines and cross smooth solution breaks down integral form (time integrated) still holds It is important to ensure correct shock speed 7
14 Discretization in conservative form In 1D, the FV scheme can be regarded as a FD scheme in conservative form u i (t) = 1 xi+1/2 u(x, t) dx x i x i 1/2 8
15 Discretization in conservative form In 1D, the FV scheme can be regarded as a FD scheme in conservative form u i (t) = 1 xi+1/2 u(x, t) dx x i x i 1/2 On a control volume V = [x i 1/2, x i+1/2 ], the exact solution fulfills d dt xi+1/2 x i 1/2 u dx + [f (u)] x i+1/2 x i 1/2 = 0 u n+1 i ui n discretized: x i t + f (u n i+1/2 ) f (un i 1/2 ) = 0 8
16 Discretization in conservative form In 1D, the FV scheme can be regarded as a FD scheme in conservative form u i (t) = 1 xi+1/2 u(x, t) dx x i x i 1/2 On a control volume V = [x i 1/2, x i+1/2 ], the exact solution fulfills d dt xi+1/2 x i 1/2 u dx + [f (u)] x i+1/2 x i 1/2 = 0 u n+1 i ui n discretized: x i t + f (u n i+1/2 ) f (un i 1/2 ) = 0 flux values f (u i±1/2 ) depending on unknown face quantities u i 1/2, u i+1/2 reconstruction necessary from available data..., u i 1, u i, u i+1,... Introduction of numerical flux functions f u n+1 ( i f (ui n, ui+1) n f (ui 1, n ui n ) ) = u n i x i t 8
17 Discretization in conservative form In 1D, the FV scheme can be regarded as a FD scheme in conservative form u i (t) = 1 xi+1/2 u(x, t) dx x i x i 1/2 On a control volume V = [x i 1/2, x i+1/2 ], the exact solution fulfills d dt xi+1/2 x i 1/2 u dx + [f (u)] x i+1/2 x i 1/2 = 0 u n+1 i ui n discretized: x i t + f (u n i+1/2 ) f (un i 1/2 ) = 0 flux values f (u i±1/2 ) depending on unknown face quantities u i 1/2, u i+1/2 reconstruction necessary from available data..., u i 1, u i, u i+1,... Introduction of numerical flux functions f The heart of FV schemes u n+1 i = u n i x i t ( f (u n i, u n i+1) f (u n i 1, u n i ) ) 8
18 Classical numerical flux functions linked to Riemann problems & characteristic directions 9
19 Classical numerical flux functions linked to Riemann problems & characteristic directions Upwind methods Scalar linear equation a > 0 ( u t + a u x = 0) u n+1 i = u n i x t (aun i au n i 1) (f (u i, u i+1 ) = au i ) Linear system of equations f (u i, u i+1 ) = A + u i + A u i+1 u n+1 i Nonlinear systems = u n i t ( A + (u n i u n i 1) + A (u n i+1 u n i ) ) x Flux vector splitting f(u) = f + (u) + f (u) f (u i, u i+1 ) = f + (u i ) + f (u i+1 ) Steger & Warming, van Leer, AUSM and variants 9
20 Classical approximate Riemann solvers for Euler equations Roe scheme exact solution to linear Riemann problem u t + A LR(u L, u R ) u x = 0 10
21 Classical approximate Riemann solvers for Euler equations Roe scheme exact solution to linear Riemann problem u t + A LR(u L, u R ) u x = 0 Properties of Roe matrix A LR A LR A(u) = Df(u) A LR (u, u) = A(u) A LR is diagonalizable f(u R ) f(u L ) = A LR (u R u L ) (mean value property) entropy-fix needed 10
22 Classical approximate Riemann solvers for Euler equations Roe scheme exact solution to linear Riemann problem u t + A LR(u L, u R ) u x = 0 Properties of Roe matrix A LR A LR A(u) = Df(u) A LR (u, u) = A(u) A LR is diagonalizable f(u R ) f(u L ) = A LR (u R u L ) (mean value property) entropy-fix needed HLL scheme Godunov-type scheme 10
23 Classical approximate Riemann solvers for Euler equations Roe scheme exact solution to linear Riemann problem u t + A LR(u L, u R ) u x = 0 Properties of Roe matrix A LR A LR A(u) = Df(u) A LR (u, u) = A(u) A LR is diagonalizable f(u R ) f(u L ) = A LR (u R u L ) (mean value property) entropy-fix needed HLL scheme Godunov-type scheme approximates only one intermediate state based on integral conservation law 10
24 What about higher order schemes?
25 Challenges posed by hyperbolic conservation laws src: NOAA gust front CC-BY-SA-3.0 bow waves Computation of discontinuous solutions (shocks) Unphysical oscillations Oscillations of approximate solution Needs additional numerical dissipation 11
26 Difficulties regarding discontinuous solutions Example: Burgers equation u t + u u x = 0, u(x, 0) = sin(x) Exact solution 12
27 Difficulties regarding discontinuous solutions Example: Burgers equation u t + u u x = 0, u(x, 0) = sin(x) Exact solution Projection Goal: Controll oscillations 12
28 A first approach towards higher order FV schemes with linear reconstruction modify left and right states u n+1 i = u n i + x i t ( f (u n i, u n i+1) f (u n i 1, u n i ) ) = 0 13
29 A first approach towards higher order FV schemes with linear reconstruction modify left and right states u n+1 i = u n i + x i t ( f (u n i, u n i+1) f (u n i 1, u n i ) ) = 0 linear reconstruction within cells u(x) = u i + s(x x i ) preserve integral means how to compute slopes s? prevent creation of new max/min 13
30 A first approach towards higher order FV schemes with linear reconstruction modify left and right states u n+1 i = u n i + x i t ( f (u n i, u n i+1) f (u n i 1, u n i ) ) = 0 linear reconstruction within cells u(x) = u i + s(x x i ) preserve integral means how to compute slopes s? prevent creation of new max/min enforce TVD property (relates to properties of exact solution) i u n+1 i+1 u n+1 i i u n i+1 u n i sufficient condition { } x si x s i 0, 2 u i u i 1 u i+1 u i 13
31 A first approach towards higher order FV schemes with linear reconstruction modify left and right states u n+1 i = u n i + x i t ( f (u n i, u n i+1) f (u n i 1, u n i ) ) = 0 linear reconstruction within cells u(x) = u i + s(x x i ) preserve integral means how to compute slopes s? prevent creation of new max/min enforce TVD property (relates to properties of exact solution) i u n+1 i+1 u n+1 i i u n i+1 u n i sufficient condition { } x si x s i 0, 2 u i u i 1 u i+1 u i Typical example: s i = 1 x minmod(u i+1 u i, u i u i 1 ), minmod(a, b) = a a < b, ab > 0 b a b, ab > 0 0 otherwise 13
32 The ENO reconstruction ENO stands for essentially non-oscillatory higher order reconstruction via interpolation adaptively choose stencil avoid interpolation across shocks 14
33 The ENO reconstruction ENO stands for essentially non-oscillatory higher order reconstruction via interpolation adaptively choose stencil avoid interpolation across shocks ENO approach: successive increase of polynomial degree via Newton interpolation compare divided differences obtained by left or right extension 14
34 The ENO reconstruction ENO stands for essentially non-oscillatory higher order reconstruction via interpolation adaptively choose stencil avoid interpolation across shocks ENO approach: successive increase of polynomial degree via Newton interpolation compare divided differences obtained by left or right extension ENO properties: constructs only one reconstruction polynomial prone to round off errors 14
35 The WENO reconstruction WENO stands for weighted essentially non-oscillatory A priori: choice of main (central) stencil as well as secondary stencils Compute polynomial reconstruction on each stencil, conserve integral means Compute weights depending on oscillatory behaviour of reconstruction Evaluate weighted sum of reconstruction polynomials 15
36 The WENO reconstruction WENO stands for weighted essentially non-oscillatory A priori: choice of main (central) stencil as well as secondary stencils Compute polynomial reconstruction on each stencil, conserve integral means Compute weights depending on oscillatory behaviour of reconstruction Evaluate weighted sum of reconstruction polynomials WENO properties: The higher the discretization order the higher the number of required neighbor cells Unstructured grids: difficult to construct stencils High demand on resources (CPU time / memory requirement) 15
37 From Finite Volumes..... to Discontinuous Galerkin
38 Discontinuous Galerkin schemes discontinuous approximate solutions modern space discretization d u h Φ dx + F (u h dt, u+ h, n) Φ dσ F(u h ) Φ dx = q h Φ dx V i V i V i V i FV: Φ = Φ 0 DG: Φ = Φ 0, Φ 1,..., Φ N Closer link to given physical equations High accuracy & flexibility Compact domains of dependance Highly adapted to computations in parallel 16
39 High resolution of DG scheme Double Mach reflection: Shock hitting fixed wall u u + Inflow Outflow Fixed wall Initial shock Excellent shock resolution & Detailed representation of fine structures Density distribution: FV vs. DG scheme FV scheme DG scheme (N = 4) 17
40 The triangular grid DG scheme Hyperbolic conservation law in 2D u(x, t) + F(u(x, t)) = 0, (x, t) Ω R+ t Initial conditions: Boundary conditions: u(x, 0) = u 0 (x) inflow/outflow, reflecting walls 18
41 The triangular grid DG scheme Hyperbolic conservation law in 2D u(x, t) + F(u(x, t)) = 0, (x, t) Ω R+ t Initial conditions: Boundary conditions: u(x, 0) = u 0 (x) inflow/outflow, reflecting walls Approximation u h,n (x, t): piecewise polynomial in x, degree N T h = {τ 1, τ 2,..., τ #T h} Triangulation Ω u h,n (, t) 18
42 The triangular grid DG scheme Multiplication by test functions Φ P N (τ i ), Integration over τ i d u Φ dx + F(u) Φ dx = 0 dt τ i τ i Use divergence theorem d u Φ dx + F(u) n Φ dσ F(u) Φ dx = 0 dt τ i τ i τ i u h,n H(u h,n, u+ h,n, n) F(u h,n) Nutze orthogonale Polynombasis { Φ 1, Φ 2,..., Φ q(n) } von P N (τ i ) u h,n τi (x, t) = q(n) k=1 û i k (t)φ k(x), q(n) = (N + 1)(N + 2)/2 19
43 The triangular grid DG scheme Multiplication by test functions Φ P N (τ i ), Integration over τ i d u Φ dx + F(u) Φ dx = 0 dt τ i τ i Use divergence theorem d u Φ dx + F(u) n Φ dσ F(u) Φ dx = 0 dt τ i τ i τ i u h,n F (u h,n, u+ h,n, n) F(u h,n) 20
44 The triangular grid DG scheme Multiplication by test functions Φ P N (τ i ), Integration over τ i d u Φ dx + F(u) Φ dx = 0 dt τ i τ i Use divergence theorem d u Φ dx + F(u) n Φ dσ F(u) Φ dx = 0 dt τ i τ i τ i u h,n F (u h,n, u+ h,n, n) F(u h,n) Use orthogonal polynomial basis { Φ 1, Φ 2,..., Φ q(n) } of P N (τ i ) u h,n τi (x, t) = q(n) k=1 û i k (t)φ k(x), q(n) = (N + 1)(N + 2)/2 20
45 With respect to orthogonal basis Time Evolution of coefficients ( d dt ûi k = F (u h,n, u+ h,n, n)φ kdσ + τ i Quadrature rules ) F(u h,n ) Φ k dx / Φ k 2 L 2 τ i 21
46 With respect to orthogonal basis Time Evolution of coefficients ( d dt ûi k = F (u h,n, u+ h,n, n)φ kdσ + τ i Quadrature rules ) F(u h,n ) Φ k dx / Φ k 2 L 2 τ i System of ODEs for coefficients û i k d ) dt Û(t) = L h,n (Û(t), t, Û = [ û i ] k k=1,...,q(n), i=1,...,#t h 21
47 With respect to orthogonal basis Time Evolution of coefficients ( d dt ûi k = F (u h,n, u+ h,n, n)φ kdσ + τ i Quadrature rules ) F(u h,n ) Φ k dx / Φ k 2 L 2 τ i System of ODEs for coefficients û i k d ) dt Û(t) = L h,n (Û(t), t, Û = [ û i ] k k=1,...,q(n), i=1,...,#t h e.g. Runge Kutta time integration Cockburn and Shu ( , 1998) 21
48 With respect to orthogonal basis Time Evolution of coefficients ( d dt ûi k = F (u h,n, u+ h,n, n)φ kdσ + τ i Quadrature rules ) F(u h,n ) Φ k dx / Φ k 2 L 2 τ i System of ODEs for coefficients û i k d ) dt Û(t) = L h,n (Û(t), t, Û = [ û i ] k k=1,...,q(n), i=1,...,#t h e.g. Runge Kutta time integration Cockburn and Shu ( , 1998) Allows easy incorporation of modal filters Meister, Ortleb, Sonar 12 21
49 Damping strategies for DG Modify approximate solution Modify equation Limiters (H)WENO- Reconstruction Explicit dissipation u + F(u) = ɛ u t Use neighboring data at shocks Often N = 1 Cockburn,Shu 89 Krivodonova 07 Use stencil data Weighted interpol. polynomials Qiu,Shu 04/05 In conservation law or discretization Time step O(h 2 ) Jaffre,Johnson,Szepessy 95 Persson,Peraire 06 Feistauer,Kučera 07 22
50 Damping for spectral methods 1D periodic case: t u + x f (u) = 0, x [ π, π] Fourier method: t u N + x P Nf (u N ) = 0 u N (x, t) = û k (t)e ikx k N Modal Filtering Modify coefficients at times t n un(x,t σ n ) = ( ) k σ ûk n e ikx N k N σ : [0, 1] [0, 1] Gottlieb, Lustman, Orzag 81 23
51 Damping for spectral methods 1D periodic case: t u + x f (u) = 0, x [ π, π] Fourier method: t u N + x P Nf (u N ) = 0 u N (x, t) = û k (t)e ikx k N Modal Filtering Modify coefficients at times t n un(x,t σ n ) = ( ) k σ ûk n e ikx N k N Spectral viscosity Add special viscosity term [ p+1 p p u ] N ɛ N ( 1) x p Q N x p Gottlieb, Lustman, Orzag 81 Tadmor 89 24
52 Damping for spectral methods 1D periodic case: t u + x f (u) = 0, x [ π, π] Fourier method: t u N + x P Nf (u N ) = 0 u N (x, t) = û k (t)e ikx k N Modal Filtering Modify coefficients at times t n un(x,t σ n ) = ( ) k σ ûk n e ikx N k N Simple implementation Link Spectral viscosity Add special viscosity term [ p+1 p p u ] N ɛ N ( 1) x p Q N x p Convergence theory, Parameter choice 25
53 Shock-density interaction (Shu-Osher test case) Initial conditions (ρ, v 1, v 2, p) = { ( , , 0, ) if x < 4 ( sin(5x), 0, 0, 1) if x 4 DG with modal filtering approximate density solution (t = 1.8) Discretization: N = 5, K =
54 Further classical test cases forward facing step density N = 5, K = 6496 double Mach reflection density N = 4, K = shock-vortex interaction pressure N = 7, K =
55 Finite Volumes & Discontinuous Galerkin..... and beyond?
56 Contents 1 The Finite Volume Method 2 The Discontinuous Galerkin Scheme 3 SBP Operators & Flux Reconstruction 4 Current High Performance DG / FR Schemes
57 Conservative schemes in 1D Finite Volume Method Discontinuous Galerkin Method Cockburn/Shu 89 Spectral Difference Method Kopriva/Kolias 96, Liu et al. 06, Wang et al. 07 Flux Reconstruction Method Huynh 11 VCJH Energy Stable FR Method Vincent et al. 10 SBP-SAT schemes Originally: Kreiss/Scherer 74 28
58 The 1D DG scheme Scalar hyperbolic conservation law in 1D u(x, t) + f (u(x, t)) = 0, t > 0, x Ω = [α, β] t x 29
59 The 1D DG scheme Scalar hyperbolic conservation law in 1D u(x, t) + f (u(x, t)) = 0, t > 0, x Ω = [α, β] t x Subdivision of Ω in Ω = Ω i = [x i, x i+1 ] i i approximation of u on Ω i uh i N+1 (x, t) = u i l (t)φi l (x) k=1 u i = (u1 i,..., ui p+1 )T with basis functions Φi l solution vector The DG scheme in strong form Ω i uh i f t Φi h kdx+ Ωi i x Φi kdx = [fi 1,i f h i (x i )]Φ i k(x i ) [fi,i+1 f h i (x i+1 )]Φ i k(x i+1 ) 29
60 The 1D DG scheme Scalar hyperbolic conservation law in 1D u(x, t) + f (u(x, t)) = 0, t > 0, x Ω = [α, β] t x Subdivision of Ω in Ω = Ω i = [x i, x i+1 ] i i approximation of u on Ω i uh i N+1 (x, t) = u i l (t)φi l (x) k=1 u i = (u1 i,..., ui p+1 )T with basis functions Φi l solution vector The DG scheme in strong form Ω i uh i f t Φi h kdx+ Ωi i x Φi kdx = [fi 1,i f h i (x i )]Φ i k(x i ) [fi,i+1 f h i (x i+1 )]Φ i k(x i+1 ) equivalent to M i dui dt + S i f i = [(f h f )Φ i ] x i+1 x i M i kl = Ω i Φ i k Φi l dx S i kl = Ω i Φ i k x Φi l dx Φ i = (Φ i 1,..., Φi p+1 )T 29
61 SBP schemes Generalized definition of 1D SBP scheme Del Rey Fernández et al. 14 M symmetric positive definite D := M 1 S approximates x S + S T = B with (x µ ) T Bx ν = (x i+1 ) µ+ν (x i ) µ+ν SBP mimics integration by parts 30
62 SBP schemes Generalized definition of 1D SBP scheme Del Rey Fernández et al. 14 M symmetric positive definite D := M 1 S approximates x S + S T = B with (x µ ) T Bx ν = (x i+1 ) µ+ν (x i ) µ+ν SBP mimics integration by parts fulfilled by DG scheme M du dt + Sf = [(f h f )Φ] x i+1 x i M kl = Ω i Φ k Φ l dx S kl = Ω i Φ k x Φ ldx B kl = [Φ k Φ l ] x i+1 x i Gauss-Lobatto (GLL) and Gauss-Legendre (GL) DG schemes: B GLL = diag{ 1, 0,..., 0, 1}, B GL,N=1 = diag{ 3, 3}, B GL,N=2 = 1 1 ξ 2 0 ξ 3 ξ 3 1 ξ 2 ξ ξ 3 ξ 3 ξ ξ 3 ξ 3, ξ =
63 Advantages of SBP schemes Provable linear stability Energy stability w.r.t. 1 2 u 2 M 31
64 Advantages of SBP schemes Provable linear stability Energy stability w.r.t. 1 2 u 2 M Relation to quadrature formulae Corresponding QF preserve certain properties of functional e.g. discrete divergence theorem Hicken/Zingg 13 31
65 Advantages of SBP schemes Provable linear stability Energy stability w.r.t. 1 2 u 2 M Relation to quadrature formulae Corresponding QF preserve certain properties of functional e.g. discrete divergence theorem Hicken/Zingg 13 Correct discretization of split form conservation laws split forms (e.g. skew-symmetric) Better control of oscillations preservation of secondary quantities, e.g. kinetic energy 31
66 Advantages of SBP schemes Provable linear stability Energy stability w.r.t. 1 2 u 2 M Relation to quadrature formulae Corresponding QF preserve certain properties of functional e.g. discrete divergence theorem Hicken/Zingg 13 Correct discretization of split form conservation laws split forms (e.g. skew-symmetric) Better control of oscillations preservation of secondary quantities, e.g. kinetic energy possible lack of discrete conservation 31
67 Advantages of SBP schemes Provable linear stability Energy stability w.r.t. 1 2 u 2 M Relation to quadrature formulae Corresponding QF preserve certain properties of functional e.g. discrete divergence theorem Hicken/Zingg 13 Correct discretization of split form conservation laws split forms (e.g. skew-symmetric) Better control of oscillations preservation of secondary quantities, e.g. kinetic energy possible lack of discrete conservation SBP schemes: equivalent telescoping form Fisher et al. 12 if convergent, then weak solution (Lax-Wendroff) 31
68 Observation Discontinuous Galerkin Energy Stable Flux Reconstruction (VCJH) 32
69 Observation Discontinuous Galerkin Energy Stable Flux Reconstruction (VCJH) These methods meet as SBP schemes! 32
70 Spectral difference and flux reconstruction schemes The SD scheme [Wang et al. 07] SD u f + t x = 0 Construction of f SD f SD f h u h 33
71 Spectral difference and flux reconstruction schemes The SD scheme [Wang et al. 07] SD u f + t x = 0 Construction of f SD Generalized by FR scheme [Huynh 11] f SD f FR = fh i + [fi 1,i fh i (x i )] g L + [fi,i+1 fh i (x i+1 )] g R }{{}}{{} f CL { f CR where g L, g R P N+1 gl (x i ) = 1 g with R (x i ) = 0 g L (x i+1 ) = 0 g R (x i+1 ) = 1 f h u h DG for g L, g R right & left Radau polynomials 33
72 Spectral difference and flux reconstruction schemes The SD scheme [Wang et al. 07] SD u f + t x = 0 Construction of f SD Generalized by FR scheme [Huynh 11] f SD f FR = fh i + [fi 1,i fh i (x i )] g L + [fi,i+1 fh i (x i+1 )] g R }{{}}{{} f CL { f CR where g L, g R P N+1 gl (x i ) = 1 g with R (x i ) = 0 g L (x i+1 ) = 0 g R (x i+1 ) = 1 in matrix-vector form [Allaneau/Jameson 11] f h u h DG for g L, g R right & left Radau polynomials expand u h, f h and g L, g R in same basis {Φ k } multiply by M with M kl = Ω i Φ k Φ l dx M du dt + S f = f CLM g L f CRM g R 33
73 Spectral difference and flux reconstruction schemes The SD scheme [Wang et al. 07] SD u f + t x = 0 Construction of f SD Generalized by FR scheme [Huynh 11] f SD f FR = fh i + [fi 1,i fh i (x i )] g L + [fi,i+1 fh i (x i+1 )] g R }{{}}{{} f CL { f CR where g L, g R P N+1 gl (x i ) = 1 g with R (x i ) = 0 g L (x i+1 ) = 0 g R (x i+1 ) = 1 in matrix-vector form [Allaneau/Jameson 11] f h u h DG for g L, g R right & left Radau polynomials expand u h, f h and g L, g R in same basis {Φ k } multiply by M with M kl = Ω i Φ k Φ l dx M du dt + S f = f CLM g f CRM g L R [ du ] dt + D f = f CLg f CRg, D = L R M 1 S 33
74 Relation of FR framework to DG scheme Reformulate M du dt + Sf = f CLM g L f CRM g R f CL = f i 1,i f i h (x i ) f CR = f i,i+1 f i h (x i+1) 34
75 Relation of FR framework to DG scheme Reformulate f CL = f i 1,i f i h (x i ) M du dt + Sf = f CLM g L f CRM g f CR = fi,i+1 f h i (x i+1) R As M g = g Φdx = [gφ] x i+1 x i gφ dx ( g is representation of g ) Ω i Ω i 34
76 Relation of FR framework to DG scheme Reformulate f CL = f i 1,i f i h (x i ) M du dt + Sf = f CLM g L f CRM g f CR = fi,i+1 f h i (x i+1) R As M g = g Φdx = [gφ] x i+1 x i gφ dx ( g is representation of g ) Ω i Ω i we have (as g R (x i ) = g L (x i+1 ) = 0) RHS FR = f CL Φ(x i ) f CR Φ(x i+1 ) }{{} RHS DG + (f CL g L + f CR g R ) Φ dx Ω } i {{} Deviation from DG 34
77 Relation of FR framework to DG scheme Reformulate f CL = f i 1,i f i h (x i ) M du dt + Sf = f CLM g L f CRM g f CR = fi,i+1 f h i (x i+1) R As M g = g Φdx = [gφ] x i+1 x i gφ dx ( g is representation of g ) Ω i Ω i we have (as g R (x i ) = g L (x i+1 ) = 0) RHS FR = f CL Φ(x i ) f CR Φ(x i+1 ) }{{} RHS DG + Two equivalent formulations for FR scheme: M du dt + S f = RHS DG + (f CL g L + f CR g R ) Φ dx Ω i du dt + D f = f CL g L f CR g R (f CL g L + f CR g R ) Φ dx Ω } i {{} Deviation from DG 34
78 VCJH schemes: energy stable FR Derivation of energy stable FR schemes: M du dt + S f = RHS DG + (f CL g L + f CR g R ) Φ dx Ω i K du dt + K Df = f CLK g f CRK g, K pos. semidef. with K D = 0 L R 35
79 VCJH schemes: energy stable FR Derivation of energy stable FR schemes: M du dt + S f = RHS DG + (f CL g L + f CR g R ) Φ dx Ω i K du dt + K Df = f CLK g f CRK g, K pos. semidef. with K D = 0 L R Summing up yields ( ) ( ) ( ) du M + K dt +Sf = RHS DG +f CL g L Φ dx K g + f CR g R Φ dx K g L R Ω i Ω i 35
80 VCJH schemes: energy stable FR Derivation of energy stable FR schemes: M du dt + S f = RHS DG + (f CL g L + f CR g R ) Φ dx Ω i K du dt + K Df = f CLK g f CRK g, K pos. semidef. with K D = 0 L R Summing up yields ( ) ( ) ( ) du M + K dt +Sf = RHS DG +f CL g L Φ dx K g + f CR g R Φ dx K g L R Ω i Ω i VCJH schemes [Vincent/Castonguay/Jameson 10] Choose g L, g R such that red terms vanish for suitable K 35
81 VCJH schemes: energy stable FR Derivation of energy stable FR schemes: M du dt + S f = RHS DG + (f CL g L + f CR g R ) Φ dx Ω i K du dt + K Df = f CLK g f CRK g, K pos. semidef. with K D = 0 L R Summing up yields ( ) ( ) ( ) du M + K dt +Sf = RHS DG +f CL g L Φ dx K g + f CR g R Φ dx K g L R Ω i Ω i VCJH schemes [Vincent/Castonguay/Jameson 10] Choose g L, g R such that red terms vanish for suitable K Similar to DG: M M + K (modified mass matrix) filtered DG scheme [Allaneau/Jameson 11] 35
82 VCJH schemes: energy stable FR Derivation of energy stable FR schemes: M du dt + S f = RHS DG + (f CL g L + f CR g R ) Φ dx Ω i K du dt + K Df = f CLK g f CRK g, K pos. semidef. with K D = 0 L R Summing up yields ( ) ( ) ( ) du M + K dt +Sf = RHS DG +f CL g L Φ dx K g + f CR g R Φ dx K g L R Ω i Ω i VCJH schemes [Vincent/Castonguay/Jameson 10] Choose g L, g R such that red terms vanish for suitable K Similar to DG: M M + K (modified mass matrix) filtered DG scheme [Allaneau/Jameson 11] Fulfills SBP property! [D = M 1 S = (M + K) 1 S] 35
83 Comparison of low order DGSBP schemes & Use of kinetic energy preservation and skew-symmetric forms
84 Smooth solutions to 1D Navier-Stokes equations Non-linear acoustic pressure wave ρ(x, 0) = 1, v(x, 0) = 1, p(x, 0) = sin(2πx), x [0, 1] periodic BC, viscosity µ = 0.002, Prandtl number Pr = 0.72 Gauss-Legendre vs. Gauss-Lobatto nodes [N = 1 on 80 cells, KEP flux, T = 20; reference: N = 3 on 500 cells] 1.01 DG Gauss Legendre DG Gauss Lobatto Reference Solution DG Gauss Legendre DG Gauss Lobatto Reference Solution Pressure Kinetic Energy x x Higher accuracy of Gauss-Legendre DG scheme. 36
85 2D decaying homogeneous turbulence Computed on cartesian grid discretizing Ω = [0, 2π] 2, periodic b.c. T = 1 T = 2 T = 5 T = 10 T = 0: Initial energy spectrum given in Fourier space by ( ) [ 2s+1 ( E(k) = as 1 k exp s + 1 ) ( ) ] 2 k 2 k p k p 2 k p for wave number k = k 2 x + k 2 y (Parameters k p = 12, a s = ) 37
86 Comparison standard DG vs. DG-KEP scheme (I) Energy spectrum T = 10 Gauss nodes, N = 1 Re= E(k) E(k) DG standard 40 cells DG skew 40 cells Reference t=10 Initial energy spectrum k DG standard 80 cells DG skew 80 cells Reference t=10 Initial energy spectrum k cells cells SBP operators allow for conservative discretization of fluid equations in skew-symmetric form. 38
87 Comparison standard DG vs. DG-KEP scheme (II) Re= E(k) E(k) DG standard 80 cells DG skew 80 cells Reference t= DG standard 160 cells DG skew 160 cells Reference t=10 Initial energy spectrum Initial energy spectrum k k cells cells Better representation of energy spectrum for KEP scheme. Specifically for in underresolved case. 39
88 Current successful implementations of DG and FR schemes
89 The FLEXI Project DG space discretization explicit time stepping massive scalability F. Hindenlang, G. J. Gassner, C. Altmann, A. Beck, M. Staudenmaier, C. Munz, Explicit discontinuous Galerkin methods for unsteady problems, Computers & fluids 61, pp ,
90 Flux reconstrution with Python: PyFR high order flux reconstruction parallel platforms team named a 2016 finalist for Gordon Bell Prize F. D. Witherden, A. M. Farrington, P. E. Vincent, PyFR: An Open Source Framework for Solving Advection-Diffusion Type Problems on Streaming Architectures using the Flux Reconstruction Approach, Computer Physics Communications 185, pp ,
91 Spectral/hp Element Framework: Nektar++ CG / DG operators hierarchical and nodal expansion bases temporal and spatial adaption, MPI parallel communication C. D. Cantwell, D. Moxey, A. Comerford, A. Bolis, G. Rocco, G. Mengaldo, D. De Grazia, S. Yakovlev, J.-E. Lombard, D. Ekelschot, B. Jordi, H. Xu, Y. Mohamied, C. Eskilsson, B. Nelson, P. Vos, C. Biotto, R.M. Kirby, S.J. Sherwin, Nektar++: An open-source spectral/ element framework, Computer Physics Communications 192, pp ,
92 DG for OpenFOAM? 4/b/b05 1/index.de.jsp developement within DG framework BoSSS (Bounded Support Spectral Solver) to be successively implemented in OpenFOAM level set method for multiphase flow high order DG method arbitrarily high accuracy at phase interfaces with cut-cell method N. Müller, S. Krämer-Eis, F. Kummer, M. Oberlack, A high-order discontinuous Galerkin method for compressible flows with immersed boundaries, Int. J. Numer. Meth. Engng. 110, pp. 3 30,
93 Summary 1 The Finite Volume Method 2 The Discontinuous Galerkin Scheme 3 SBP Operators & Flux Reconstruction 4 Current High Performance DG / FR Schemes 44
94 Thank you for your attention!
ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY ON INTERIOR NODAL SETS
6th European Conference on Computational Mechanics (ECCM 6 7th European Conference on Computational Fluid Dynamics (ECFD 7 11 15 June 018, Glasgow, UK ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationEntropy stable high order discontinuous Galerkin methods. for hyperbolic conservation laws
Entropy stable high order discontinuous Galerkin methods for hyperbolic conservation laws Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Tianheng Chen, and with Yong Liu
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationChapter 1. Introduction and Background. 1.1 Introduction
Chapter 1 Introduction and Background 1.1 Introduction Over the past several years the numerical approximation of partial differential equations (PDEs) has made important progress because of the rapid
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationHigh-Order Methods for Diffusion Equation with Energy Stable Flux Reconstruction Scheme
49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposi 4-7 January 11, Orlando, Florida AIAA 11-46 High-Order Methods for Diffusion Equation with Energy Stable Flux
More informationarxiv: v2 [math.na] 9 Apr 2018
oname manuscript o. will be inserted by the editor) The BR1 Scheme is Stable for the Compressible avier-stokes Equations Gregor J. Gassner Andrew R. Winters Florian J. Hindenlang David A. Kopriva arxiv:1704.03646v
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationAntony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationRunge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter
Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter Jun Zhu, inghui Zhong, Chi-Wang Shu 3 and Jianxian Qiu 4 Abstract In this paper, we propose a new type of weighted
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationSub-Cell Shock Capturing for Discontinuous Galerkin Methods
Sub-Cell Shock Capturing for Discontinuous Galerkin Methods Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 39, U.S.A. A shock capturing strategy for higher order
More informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationA central discontinuous Galerkin method for Hamilton-Jacobi equations
A central discontinuous Galerkin method for Hamilton-Jacobi equations Fengyan Li and Sergey Yakovlev Abstract In this paper, a central discontinuous Galerkin method is proposed to solve for the viscosity
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationEmerging Numerical Methods for Atmospheric Modeling
Chapter 9 Emerging Numerical Methods for Atmospheric Modeling Ramachandran D. Nair, Michael N. Levy and Peter H. Lauritzen Abstract This chapter discusses the development of discontinuous Galerkin (DG)
More informationNumerische Mathematik
Numer. Math. ) :545 563 DOI.7/s--443-7 Numerische Mathematik A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang Chi-Wang Shu Received: 7 July / Revised: 5 September
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationAn efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes
An efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes Christian Klingenberg, Frank Pörner, Yinhua Xia Abstract
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationarxiv: v1 [math.na] 23 Jun 2017
On the use of kinetic energy preserving DG-schemes for large eddy simulation David Flad 1 arxiv:1706.07601v1 [math.na] 23 Jun 2017 Institute for Aerodynamics and Gas Dynamics, University of Stuttgart,
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationDISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS
DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS CHI-WANG SHU Abstract. In these lectures we give a general survey on discontinuous Galerkin methods for solving
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationOn the Non-linear Stability of Flux Reconstruction Schemes
DOI 10.1007/s10915-011-9490-6 TECHNICA NOTE On the Non-linear Stability of Flux econstruction Schemes A. Jameson P.E. Vincent P. Castonguay eceived: 9 December 010 / evised: 17 March 011 / Accepted: 14
More informationarxiv: v1 [physics.comp-ph] 21 Nov 2018
Implicit High-Order Flux Reconstruction Solver for High-Speed Compressible Flows Ray Vandenhoeck a,, Andrea Lani b a Department of Mechanical Engineering, KU Leuven, Celestijnenlaan 300, 3001, Leuven,
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationEssentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws
NASA/CR-97-0653 ICASE Report No. 97-65 Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws Chi-Wang Shu Brown University Institute for Computer
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws
ADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws Martin Käser and Armin Iske Abstract. ADER schemes are recent finite volume methods for hyperbolic conservation laws, which can be
More informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS
ALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS HAILIANG LIU AND MICHAEL POLLACK Abstract. In this work, we propose a high resolution Alternating Evolution Discontinuous
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationArt Checklist. Journal Code: Article No: 6959
Art Checklist Journal Code: JCPH Article No: 6959 Disk Recd Disk Disk Usable Art # Y/N Format Y/N Remarks Fig. 1 Y PS Y Fig. 2 Y PS Y Fig. 3 Y PS Y Fig. 4 Y PS Y Fig. 5 Y PS Y Fig. 6 Y PS Y Fig. 7 Y PS
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationThe construction of discretely conservative finite volume schemes that also globally conserve energy or entropy
The construction of discretely conservative finite volume schemes that also globally conserve energy or entropy Antony Jameson Stanford University Aerospace Computing Laboratory Report ACL 007 January
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More information