Well-balanced DG scheme for Euler equations with gravity
|
|
- Marianna Tyler
- 6 years ago
- Views:
Transcription
1 Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore Dept. of Mathematics, IIT Delhi 26 October, 2016 Supported by Airbus Foundation Chair at TIFR-CAM, Bangalore 1 / 59
2 Euler equations with gravity q t + f x = where ρ ρu q = ρu, f = p + ρu 2, E (E + p)u Calorically ideal gas p = (γ 1) 0 ρ Φ x ρu Φ = gravitational potential [E 12 ρu2 ], γ = c p c v > 1 2 / 59
3 Hydrostatic solutions Fluid at rest u e = 0 Mass and energy equation satisfied t ρ + x (ρu) = 0, t E + x [(E + p)u] = ρu x Φ Momentum equation dp e dx = ρ dφ e dx Need additional assumptions to solve this equation (1) 3 / 59
4 Hydrostatic solutions: Isothermal case Thermally ideal gas model p = ρrt, R = gas constant Isothermal hydrostatic state, i.e., T e (x) = T e = const, then ( ) Φ(x) p e (x) exp = const (2) RT e Density ρ e (x) = p e(x) RT e 4 / 59
5 Hydrostatic solutions: Polytropic case Polytropic hydrostatic state p e ρ ν e = α = const, ν > 1 constant (3) From (1) and (3), we obtain ναρe ν 1 (x) + Φ(x) = β = const (4) ν 1 5 / 59
6 Scope of present work Nodal DG scheme Gauss-Lobatto-Legendre nodes Arbitrary quadrilateral cells in 2-D Well-balanced for hydrostatic solutions isothermal solutions under ideal gas polytropic solutions (with Markus Zenk) Any consistent numerical flux Discontinuous density: contact preserving flux 6 / 59
7 Finite element method Conservation law with source term stationary solution q e q t + f(q) x = s(q) f(q e ) x = s(q e ) Weak formulation: Find q(t) V such that d (q, ϕ) + a(q, ϕ) = (s(q), ϕ), dt ϕ V FEM with quadrature: Find q h (t) V h such that d dt (q h, ϕ h ) h + a h (q h, ϕ h ) = (s h (q h ), ϕ h ) h, ϕ h V h 7 / 59
8 Finite element and well-balanced Stationary solution is not a polynomial, q e / V h. Let q e,h = Π h (q e ), Π h : V V h, interpolation or projection FEM is well-balanced if a h (q e,h, ϕ h ) = (s h (q e,h ), ϕ h ) h, ϕ h V h because if q h (0) = q e,h = q h (t) = q e,h t How to construct a h, s h to achieve well-balancing? 8 / 59
9 Mesh and basis functions Partition domain into disjoint cells C i = (x i 1, x 2 i+ 1 ), x i = x 2 i+ 1 x 2 i 1 2 Approximate solution inside each cell by a polynomial of degree N I i 1 I i I i 1 9 / 59
10 Mesh and basis functions Map C i to a reference cell, say [0, 1] x = ξ x i + x i 1 2 (5) On reference cell, ξ j, 0 j N are Gauss-Lobatto-Legendre nodes, roots of (1 ξ 2 )P N (ξ) in [ 1, +1] l j (ξ) = Lagrange polynomials using GLL points l j (ξ k ) = δ jk, 0 j, k N Basis functions in physical coordinates ϕ j (x) = l j (ξ), 0 j N Derivative of ϕ j : apply the chain rule of differentiation d dx ϕ j(x) = l j(ξ) dξ dx = 1 x i l j(ξ) 10 / 59
11 Mesh and basis functions x j C i denote the physical locations of the GLL points x j = ξ j x i + x i 1, 0 j N 2 11 / 59
12 Discontinuous Galerkin Scheme Consider the single conservation law with source term q t + f = s(q, x) x Solution inside cell C i is polynomial of degree N q h (x, t) = Approximate the flux N q j (t)ϕ j (x), q j (t) = q h (x j, t) j=0 f(q h ) f h (x, t) = N f(q h (x j, t))ϕ j (x) = j=0 N f j (t)ϕ j (x) j=0 12 / 59
13 Discontinuous Galerkin Scheme Gauss-Lobatto-Legendre quadrature N φ(x)ψ(x)dx (φ, ψ) h = x i ω q φ(x q )ψ(x q ) C i q=0 Semi-discrete DG: For 0 j N d dt (q h, ϕ j ) h + ( x f h, ϕ j ) h +[ ˆf i+ 1 f h (x )]ϕ 2 i+ 1 j (x ) i [ ˆf (6) i 1 f h (x + )]ϕ 2 i 1 j (x + ) = (s i 1 h, ϕ j ) h 2 2 where ˆf i+ 1 = ˆf(q, q + ) is a numerical flux function. This is also 2 i+ 1 i called a DG Spectral Element Method. 13 / 59
14 Numerical flux Consistency of numerical flux ˆf(q, q) = f(q) Def: Contact property The numerical flux ˆf is said to satisfy contact property if for any two states w L = [ρ L, 0, p] and w R = [ρ R, 0, p] we have ˆf(q(w L ), q(w R )) = [0, p, 0] states w L, w R in the above definition correspond to a stationary contact discontinuity. Contact Property = exactly supports stationary contact Examples: Roe, HLLC, etc. 14 / 59
15 Approximation of source term: isothermal case Let T i = temperature corresponding to the cell average value in cell C i Rewrite the source term in the momentum equation as (Xing & Shu) s = ρ Φ ( ) ( x = ρr T Φ i exp R T i x exp Φ ) R T i Source term approximation: For x C i ( ) s h (x) = ρ h (x)r T Φ(x) N ( i exp R T exp Φ(x ) j) i x R T ϕ j (x) (7) i Source term in the energy equation 1 ρ h (ρu) h s h j=0 15 / 59
16 Approximation of source term: polytropic case Define H(x) inside each cell C i as H(x) = ν [ ] ν 1 ν 1 ln (β i Φ(x)), x C i να i α i, β i : constants to be chosen. Rewrite source term s(x) = ρ Φ x = ν 1 ρ(β i Φ(x)) exp( H(x)) ν x exp(h(x)), The source term is approximated as x C i s h (x) = ν 1 ρ h (x)(β i Φ(x)) exp( H(x)) ν x [ ν β i = max 0 j N ν 1 ] p j + Φ(x j ) ρ j N exp(h(x j ))ϕ j (x) j=0, α i = p j ρ ν j 16 / 59
17 Well-balanced property Well-balanced property Let the initial condition to the DG scheme (6), (7) be obtained by interpolating the isothermal/polytropic hydrostatic solution corresponding to a continuous gravitational potential Φ. Then the scheme (6), (7) preserves the initial condition under any time integration scheme. Proof: For continuous hydrostatic solution, q h (0) is continuous. By flux consistency ˆf i+ 1 f h (x ) = 0, ˆfi 2 i+ 1 1 f h (x + ) = i 1 2 Above is true even if density is discontinuous, provided flux satisfies contact property. = density and energy equations are well-balanced 17 / 59
18 Well-balanced property Momentum eqn: flux f h has the form N f h (x, t) = p j (t)ϕ j (x), j=0 p j = pressure at the GLL point x j Isothermal initial condition, T i = T e = const. The source term evaluated 18 / 59
19 Well-balanced property at any GLL node x k is given by s h (x k ) = ( ) Φ(xk ) N ρ h (x k )RT }{{} e exp RT e p j=0 k = ( ) Φ(xk ) N p k exp RT e = = = j=0 N ( Φ(xk ) p k exp RT e j=0 j=0 ( exp Φ(x ) j) RT e x ϕ j(x k ) ( exp Φ(x ) j) RT e x ϕ j(x k ) ) exp ( Φ(x ) j) RT e x ϕ j(x k ) N ( ) ( Φ(xj ) p j exp exp Φ(x ) j) RT e RT e x ϕ j(x k ) N j=0 p j x ϕ j(x k ) = x f h(x k ) 19 / 59
20 Well-balanced property Since x f h (x k ) = s h (x k ) at all the GLL nodes x k we can conclude that ( x f h, ϕ j ) h = (s h, ϕ j ) h, 0 j N = scheme is well-balanced for the momentum equation The proof for polytropic case is similar. 20 / 59
21 2-D Euler equations with gravity q t + f x + g y = s q = vector of conserved variables, (f, g) = flux vector and s = source term, given by ρ q = ρu ρv, f = E s = ρu p + ρu 2 ρuv (E + p)u 0 ρ Φ x ρ Φ y ( u Φ, g = ) x + v Φ y ρv ρuv p + ρv 2 (E + p)v 21 / 59
22 2-D hydrostatic solution: Isothermal case Momentum equation p e = ρ e Φ Assuming the ideal gas equation of state p = ρrt and a constant temperature T = T e = const, we get ( ) Φ(x, y) p e (x, y) exp = const (8) RT e We will exploit the above property of the hydrostatic state to construct the well-balanced scheme. 22 / 59
23 Mesh and basis functions A quadrilateral cell K and reference cell ˆK = [0, 1] [0, 1] y 1 4 K 2 3 T K η 4 ˆK 3 x 1 2 ξ 1-D GLL points ξ r [0, 1], 0 r N 23 / 59
24 Mesh and basis functions Tensor product of GLL points N = 1 N = 2 N = 3 Basis functions in cell K: i (r, s) ϕ K i (x, y) = l r (ξ)l s (η), (x, y) = T K (ξ, η) Computation of x f h requires derivatives of the basis functions x ϕk i (x, y) = l r(ξ)l s (η) ξ x + l r(ξ)l s(η) η x 24 / 59
25 DG scheme in 2-D Consider a single conservation law with source term Solution inside cell K q h (x, y, t) = q t + f x + g x = s M qi K (t)ϕ K i (x, y), M = (N + 1) 2 i=1 Approximate fluxes inside the cell by interpolation at the same GLL nodes f h (x, y) = M f(q h (x i, y i ))ϕ K i (x, y) = i=1 M f i ϕ K i (x, y) i=1 g h (x, y) = M g(q h (x i, y i ))ϕ K i (x, y) = i=1 M g i ϕ K i (x, y) i=1 25 / 59
26 DG scheme in 2-D Quadrature on element K using GLL points (φ, ψ) K = N r=0 s=0 Quadrature on the faces of the cell K N φ(ξ r, ξ s )ψ(ξ r, ξ s )ω r ω s J K (ξ r, ξ s ) (φ, ψ) K = e K (φ, ψ) e where (, ) e are one dimensional quadrature rules using the subset of GLL nodes located on the boundary of the cell. 26 / 59
27 DG scheme in 2-D Semi-discrete DG scheme: For 1 i M d ( qh, ϕ K ) i dt K + ( x f h + y g h, ϕ K ) i K ( ) + ˆFh F h, ϕk i = ( s h, ϕ K i K ) (9) Numerical flux function ˆF h = ˆF (q h, q+ h, n), n = (n x, n y ), unit outward normal to K Trace of flux on K from cell K F h = f h n x + g h n y 27 / 59
28 Approximation of source term: Isothermal case Let T K = temperature corresponding to the cell average value in cell K Rewrite source term in the x momentum equation s = ρ Φ ( ) ( x = ρr T Φ K exp R T K x exp Φ ) R T K Approximation of source term for (x, y) K ( ) s h (x, y) = ρ h (x, y)r T Φ(x, y) K exp R T K x M j=1 ( exp Φ(x ) j, y j ) R T ϕ K j (x, y) K (10) 28 / 59
29 Well-balanced property Let the initial condition to the DG scheme (9), (10) be obtained by interpolating the isothermal/polytropic hydrostatic solution corresponding to a continuous gravitational potential Φ. Then the scheme preserves the initial condition under any time integration scheme. 29 / 59
30 Limiter TVB limiter of Cockburn-Shu How to preserve hydrostatic solution? If residual in cell K is zero, then dont apply limiter in that cell. 30 / 59
31 Time integration scheme: dq dt = R(t, q) Second order accurate SSP Runge-Kutta scheme q (1) = q n + tr(t n, q n ) q (2) = 1 2 qn [q(1) + tr(t n + t, q (1) )] q n+1 = q (2) Third order accurate SSP Runge-Kutta scheme q (1) = q n + tr(t n, q n ) q (2) = 3 4 qn [q(1) + tr(t n + t, q (1) )] q (3) = 1 3 qn [q(2) + tr(t n + t/2, q (2) )] q n+1 = q (3) 31 / 59
32 Numerical Results 32 / 59
33 1-D hydrostatic test: Well-balanced property Potential Initial state u = 0, Φ = x, Φ = sin(2πx) ρ = p = exp( Φ(x)) Mesh ρu ρ E 25x e e e-14 50x e e e x e e e x e e e-14 Table: Well-balanced test on Cartesian mesh using Q1 and potential Φ = y 33 / 59
34 1-D hydrostatic test: Well-balanced property Mesh ρu ρ E 25x e e e-14 50x e e e x e e e x e e e-14 Table: Well-balanced test on Cartesian mesh using Q2 and potential Φ = x Mesh ρu ρ E 25x e e e-13 50x e e e x e e e x e e e-13 Table: Well-balanced test on Cartesian mesh using Q1 and Φ = sin(2πx) 34 / 59
35 1-D hydrostatic test: Well-balanced property Mesh ρu ρ E 25x e e e-13 50x e e e x e e e x e e e-13 Table: Well-balanced test on Cartesian mesh using Q2 and Φ = sin(2πx) 35 / 59
36 1-D hydrostatic test: Evolution of perturbations Potential Φ(x) = x Initial condition u = 0, ρ = exp( x), p = exp( x) + η exp( 100(x 1/2) 2 ) 36 / 59
37 1-D hydrostatic test: Evolution of perturbations Initial 100 cells 1000 cells Initial 50 cells 500 cells Pressure perturbation Pressure perturbation x (a) x (b) Figure: Evolution of perturbations for η = 10 2 : (a) Q1 (b) Q2 37 / 59
38 1-D hydrostatic test: Evolution of perturbations Initial 100 cells 1000 cells Initial 50 cells 500 cells Pressure perturbation Pressure perturbation x (a) x (b) Figure: Evolution of perturbations for η = 10 4 : (a) Q1 (b) Q2 38 / 59
39 Shock tube problem Potential Φ(x) = x and initial condition { (1, 0, 1) x < 1 (ρ, u, p) = 2 (0.125, 0, 0.1) x > Q1, 100 cells Q2, 100 cells FVM, 2000 cells Q1, 200 cells Q2, 200 cells FVM, 2000 cells Density 0.6 Density x x 39 / 59
40 2-D hydrostatic test: Well-balanced test Potential Φ = x + y Hydrostatic solution is given by ( ρ = ρ 0 exp ρ ) ( 0g (x + y), p = p 0 exp ρ ) 0g (x + y), u = v = 0 p 0 p 0 ρu ρv ρ E Q1, e e e e-13 Q1, e e e e-13 Q1, e e e e-13 Q2, e e e e-13 Q2, e e e e-13 Q2, e e e e / 59
41 2-D hydrostatic test: Evolution of perturbation ( p = p 0 exp ρ ) ( 0g (x + y) + η exp 100 ρ ) 0g [(x 0.3) 2 + (y 0.3) 2 ] p 0 p 0 41 / 59
42 2-D hydrostatic test: Evolution of perturbation y y x (a) x Figure: Pressure perturbation on mesh at time t = 0.15 (a) Q 1 (b) Q 2. Showing 20 contours lines between to (b) 42 / 59
43 2-D hydrostatic test: Evolution of perturbation y y x (a) x Figure: Pressure perturbation on mesh at time t = 0.15 (a) Q 1 (b) Q 2. Showing 20 contours lines between to (b) 43 / 59
44 2-D hydrostatic test: Discontinuous density Potential Temperature is discontinuous T = Hydrostatic pressure and density p = Φ(x, y) = y { T l, y < 0 T u, y > 0 { p 0 e y/rt l, y 0 p 0 e y/rtu, y > 0, ρ = Density is discontinuous at y = 0. { p/rt l, y < 0 p/rt u, y > 0 44 / 59
45 2-D hydrostatic test: Discontinuous density Case 1, T l = 1, T u = 2 Q1, 25x e e e e-13 Q1, 50x e e e e-13 Q2, 25x e e e e-13 Q2, 50x e e e e-13 Table: Well-balanced test for Rayleigh-Taylor problem Case 2, T l = 2, T u = 1 Q1, 25x e e e e-13 Q1, 50x e e e e-13 Q2, 25x e e e e-13 Q2, 50x e e e e-13 Table: Well-balanced test for Rayleigh-Taylor problem 45 / 59
46 2-D hydrostatic test: Discontinuous density Ligther fluid on top 46 / 59
47 2-D hydrostatic test: Discontinuous density Heavier fluid on top 47 / 59
48 Order of accuracy Exact solution of Euler equations with gravity given by (Xing & Shu) ρ = sin(π(x + y t(u 0 + v 0 ))), u = u 0, v = v 0 p = p 0 + t(u 0 + v 0 ) x y cos(π(x + y t(u 0 + v 0 )))/π Compute solution error in L 2 norm at time t = 0.1 1/h ρu ρv ρ E Error Rate Error Rate Error Rate Error Rate e e e e e e e Table: Convergence of error for degree N = 1 1/h ρu ρv ρ E Error Rate Error Rate Error Rate Error Rate e e e e e e e e e e e e e e e e Table: Convergence of error for degree N = 2 48 / 59
49 Radial Rayleigh-Taylor problem cells 956 cells 49 / 59
50 Radial Rayleigh-Taylor problem: Well-balanced test Radial potential Hydrostatic solution Φ(x, y) = r = x 2 + y 2 ρ = p = exp( r) ρu ρv ρ E Q1, 30x e e e e-14 Q1, 50x e e e e-14 Q1, 100x e e e e-14 Q2, 30x e e e e-14 Q2, 50x e e e e-14 Q2, 100x e e e e-14 Table: Well balanced test for radial Rayleigh-Taylor problem on Cartesian mesh 50 / 59
51 Radial Rayleigh-Taylor problem: Well-balanced test ρu ρv ρ E Q1, e e e e-16 Q1, e e e e-16 Q1, e e e e-15 Q2, e e e e-13 Q2, e e e e-13 Q2, e e e e-13 Table: Well balanced test for radial Rayleigh-Taylor problem on unstructured mesh 51 / 59
52 Radial Rayleigh-Taylor problem: Perturbations Initial pressure and density where { { e r r r 0 e r r r i p = e r α +r 0 (1 α), ρ = α r > r 0 r > r i 1 α e r α +r (1 α) 0 α r i = r 0 (1 + η cos(kθ)), α = exp( r 0 )/(exp( r 0 ) + ρ ) The density jumps by an amount ρ > 0 at the interface defined by r = r i whereas the pressure is continuous. Following [?], we take ρ = 0.1, η = 0.02, k = / 59
53 Radial Rayleigh-Taylor problem: Perturbations Cartesian mesh of and Q1 basis 53 / 59
54 Well-balanced versus non well-balanced (a) (b) (c) Well balanced test for radial Rayleigh-Taylor problem on mesh. (a) Initial density, (b) well-balanced scheme, density at t = 1.5, (c) non well-balanced scheme, density at t = / 59
55 Summary Well-balanced DG scheme for Euler with gravity Preserves isothermal hydrostatic solutions for ideal gas model Preserves polytropic hydrostatic solution Any numerical flux can be used Cartesian, unstructured meshes 55 / 59
56 Extensions Define variables v = [ρ/α, u, p/β] Interpolate v onto V h v h = Π h v(q h ) Source DG scheme in 1-D s h (x, t) = ρ h(x, t) p 0 α h (x) ρ 0 x d dt (q h, ϕ j ) h + ( x f h (v h ), ϕ j ) h N β j ϕ j (x) j=0 +[{ ˆf(q(v h ), q(v+ h )) f(q(v h ))}ϕ j ] i+ 1 2 [{ ˆf(q(v h ), q(v+ h )) f(q(v+ h ))}ϕ+ j ] i 1 2 = (s h (q h ), ϕ j ) h 56 / 59
57 Extensions GL nodes can be used Well-balanced for any equation of state Well-balanced on adaptive meshes with hanging nodes 57 / 59
58 Thank You 58 / 59
59 References 59 / 59
Well-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 informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
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 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 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 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 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 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 informationFinite volume scheme for Euler equations with gravity which is well-balanced for general equations of state and grid systems
Finite volume scheme for Euler equations with gravity which is well-balanced for general equations of state and grid systems Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg,, and Friedrich
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 informationOn divergence-free reconstruction schemes for CED and MHD
On divergence-free reconstruction schemes for CED and MHD Praveen Chandrashekar praveen@math.tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore-560065, India
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 informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
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 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 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 informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
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 informationON 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 informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationModeling Shock Waves Using Exponential Interpolation Functions with the Least-Squares Finite Element Method
Old Dominion University ODU Digital Commons Mechanical & Aerospace Engineering Theses & Dissertations Mechanical & Aerospace Engineering Spring 26 Modeling Shock Waves Using Exponential Interpolation Functions
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 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 informationA Comparison of Artificial Viscosity Sensors for the Discontinuous Galerkin Method. Joshua Favors
A Comparison of Artificial Viscosity Sensors for the Discontinuous Galerkin Method by Joshua Favors A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationA Robust CFL Condition for the Discontinuous Galerkin Method on Triangular Meshes
A Robust CFL Condition for the Discontinuous Galerkin Method on Triangular Meshes N. Chalmers and L. Krivodonova Abstract When the discontinuous Galerkin (DG) method is applied to hyperbolic problems in
More informationHybrid DG/FV schemes for magnetohydrodynamics and relativistic hydrodynamics
Hybrid DG/FV schemes for magnetohydrodynamics and relativistic hydrodynamics Jonatan Núñez-De La Rosa, Claus-Dieter Munz To cite this version: Jonatan Núñez-De La Rosa, Claus-Dieter Munz. Hybrid DG/FV
More informationEarth System Modeling Domain decomposition
Earth System Modeling Domain decomposition Graziano Giuliani International Centre for Theorethical Physics Earth System Physics Section Advanced School on Regional Climate Modeling over South America February
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 informationDiscontinuous Galerkin methods Lecture 2
y y RMMC 2008 Discontinuous Galerkin methods Lecture 2 1 Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu y 1 0.75 0.5 0.25 0-0.25-0.5-0.75 y 0.75-0.0028-0.0072-0.0117 0.5-0.0162-0.0207-0.0252
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
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 informationK. Ambika and R. Radha
Indian J. Pure Appl. Math., 473: 501-521, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0200-9 RIEMANN PROBLEM IN NON-IDEAL GAS DYNAMICS K. Ambika and R. Radha School of Mathematics
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 informationTides in Higher-Dimensional Newtonian Gravity
Tides in Higher-Dimensional Newtonian Gravity Philippe Landry Department of Physics University of Guelph 23 rd Midwest Relativity Meeting October 25, 2013 Tides: A Familiar Example Gravitational interactions
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 New Class of High-Order Energy Stable Flux Reconstruction Schemes for Triangular Elements
DOI 10.1007/s10915-011-9505-3 A New Class of High-Order Energy Stable Flux Reconstruction Schemes for Triangular Elements P. Castonguay P.E. Vincent A. Jameson Received: December 010 / Revised: 7 April
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 informationEfficient conservative ADER schemes based on WENO reconstruction and space-time predictor in primitive variables
Zanotti and Dumbser RESEARCH arxiv:1511.4728v2 [math.na] 23 Dec 215 Full list of author information is available at the end of the article Efficient conservative ADER schemes based on WENO reconstruction
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 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 informationTectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop
.. Tectonic evolution at mid-ocean ridges: geodynamics and numerical modeling. Second HPC-GA Workshop Marco Cuffaro, Edie Miglio, Mattia Penati, Marco Viganò Politecnico di Milano - MOX, Dipartimento di
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 informationA Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations
A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations Institute for Mathematics Applied to the Geosciences (IMAGe) National Center for Atmospheric Research (NCAR) Boulder CO
More informationA second-order asymptotic-preserving and positive-preserving discontinuous. Galerkin scheme for the Kerr-Debye model. Abstract
A second-order asymptotic-preserving and positive-preserving discontinuous Galerkin scheme for the Kerr-Debye model Juntao Huang 1 and Chi-Wang Shu Abstract In this paper, we develop a second-order asymptotic-preserving
More informationA lumped nodal DGTD-PIC method to ensure charge conservation for the 3D Vlasov-Maxwell system on nonconforming grids
A lumped nodal DGTD-PIC method to ensure charge conservation for the 3D Vlasov-Maxwell system on nonconforming grids B. Fornet Nuclétudes Wednesday, October 13 1 Some Physics 2 DGTD-PIC : why and how 3
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationA positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations
A positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations Sashank Srinivasan a, Jonathan Poggie a, Xiangxiong Zhang b, a School of Aeronautics and Astronautics,
More informationSpectral Difference Method for Unstructured Grids II: Extension to the Euler Equations
Journal of Scientific Computing, Vol. 32, No. 1, July 2007 ( 2006) DOI: 10.1007/s10915-006-9113-9 Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations Z. J. Wang, 1,4
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 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 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 informationarxiv: v1 [math.na] 23 Sep 2015
An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry Niklas Wintermeyer a, Andrew R. Winters
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationMaximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments
Maximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments By Xiangxiong Zhang 1 and Chi-Wang Shu 1 Department of Mathematics, Brown University,
More informationIFE for Stokes interface problem
IFE for Stokes interface problem Nabil Chaabane Slimane Adjerid, Tao Lin Virginia Tech SIAM chapter February 4, 24 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2 Problem statement
More informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationImproving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on
Improving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on Matthew R. Norman Scientific Computing Group National Center for Computational Sciences Oak Ridge
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationOn Numerical Simulation of Flow Problems by Discontinuous Galerkin and Finite Volume Techniques
On Numerical Simulation of Flow Problems by Discontinuous Galerkin and Finite Volume Techniques Ondřej Winter 1, Petr Sváček 1 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Technical
More informationFrom Finite Volumes to Discontinuous Galerkin and Flux Reconstruction
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
More informationA high-order triangular discontinuous Galerkin oceanic shallow water model
Calhoun: The NPS Institutional Archive Faculty and Researcher Publications Faculty and Researcher Publications 8 A high-order triangular discontinuous Galerkin oceanic shallow water model Giraldo, F. X.
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 informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
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 informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
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 informationLocally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials
Locally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials Harald Klimach, Michael Gaida, Sabine Roller harald.klimach@uni-siegen.de 26th WSSP 2017 Motivation Fluid-Dynamic
More informationA spectral finite volume transport scheme on the cubed-sphere
Applied Numerical Mathematics 57 (2007) 1021 1032 www.elsevier.com/locate/apnum A spectral finite volume transport scheme on the cubed-sphere Vani Cheruvu a,, Ramachandran D. Nair b, Henry M. Tufo a,b
More informationPositivity-preserving Lagrangian scheme for multi-material compressible flow. Abstract
Positivity-preserving Lagrangian scheme for multi-material compressible flow Juan Cheng and Chi-Wang Shu 2 Abstract Robustness of numerical methods has attracted an increasing interest in the community
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 informationMixed Mimetic Spectral Elements for Geophysical Fluid Dynamics
for Geophysical Fluid Dynamics Dave Lee Los Alamos National Laboratory Outline Connection of finite volumes to differential forms Key ideas of differential forms Differential forms for discrete data Construction
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
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 informationWell-balanced schemes for the Euler equations with gravitation
Well-balanced schemes for the Euler equations with gravitation R. Käppeli and S. Mishra Research Report No. 213-5 Seminar für Angewandte Mathematik Eidgenössische Technische Hochschule CH-892 Zürich Switzerland
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationLecture - Implementation and numerical aspects of DG-FEM in 1D
Lecture - Implementation and numerical aspects of DG-FEM in 1D Ph.D. Course: An Introduction to DG-FEM for solving partial differential equations Allan P. Engsig-Karup Scientific Computing Section DTU
More informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
More informationA gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions
Acta Mech. Sin. 2017) 334):733 741 DOI 10.1007/s10409-017-0695-2 RESEARCH PAPER A gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions Xiaodong Ren 1
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
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 informationCell-centered discontinuous Galerkin discretization for two-dimensional Lagrangian hydrodynamics
Cell-centered discontinuous Galerkin discretization for two-dimensional Lagrangian hydrodynamics François Vilar To cite this version: François Vilar. Cell-centered discontinuous Galerkin discretization
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationA New Class of High-Order Energy Stable Flux Reconstruction Schemes
DOI 10.1007/s10915-010-940-z A New Class of High-Order Energy Stable Flux Reconstruction Schemes P.E. Vincent P. Castonguay A. Jameson Received: May 010 / Revised: September 010 / Accepted: 5 September
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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 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 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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationA novel discontinuous Galerkin method using the principle of discrete least squares
A novel discontinuous Galerkin method using the principle of discrete least squares Jan Glaubitz and Philipp Öffner October 4, 2017 In this work, a novel discontinuous Galerkin DG) method is introduced
More information