Well-balanced DG scheme for Euler equations with gravity
|
|
- Leslie Freeman
- 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 Higher Order Numerical Methods for Evolutionary PDEs: Applied Mathematics Meets Astrophysical Applications (15w5134) BIRS, Banff Centre, May, 2015 Supported by Airbus Foundation Chair at TIFR-CAM, Bangalore 1 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 8 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
32 Numerical Results 32 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
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 / 55
46 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 46 / 55
47 Radial Rayleigh-Taylor problem cells 956 cells 47 / 55
48 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 48 / 55
49 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 49 / 55
50 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 [2], we take ρ = 0.1, η = 0.02, k = / 55
51 Radial Rayleigh-Taylor problem: Perturbations Cartesian mesh of and Q1 basis 51 / 55
52 Extensions: Gauss-Legendre nodes Define variables (isothermal case) Interpolate onto GLL nodes DG scheme in 1-D v = [ρ exp(φ/r T i ), u, p exp(φ/r T i )] v h = Π h v(q h ) d dt (q h, ϕ j ) h + ( x f h (q h ), ϕ j ) h +[{ ˆ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 Well-balanced on adaptive meshes with hanging nodes 52 / 55
53 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 53 / 55
54 Thank You 54 / 55
55 References [1] R. Käppeli and S. Mishra, Well-balanced schemes for the Euler equations with gravitation, J. Comput. Phys., 259 (2014), pp [2] Randall. J. LeVeque and Derek. S. Bale, Wave propagation methods for conservation laws with source terms, in Hyperbolic Problems: Theory, Numerics, Applications, Rolf Jeltsch and Michael Fey, eds., vol. 130 of International Series of Numerical Mathematics, Birkhäuser Basel, 1999, pp [3] Yulong Xing and Chi-Wang Shu, High order well-balanced WENO scheme for the gas dynamics equations under gravitational fields, J. Sci. Comput., 54 (2013), pp / 55
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 Dept. of
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationReceived 6 August 2005; Accepted (in revised version) 22 September 2005
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. -34 Commun. Comput. Phys. February 6 A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
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 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 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 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 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 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 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 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 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 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 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 informationarxiv: v1 [math.na] 10 Sep 2018
A global divergence conforming DG method for hyperbolic conservation laws with divergence constraint Praveen handrashekar arxiv:809.0394v [math.na] 0 Sep 08 Abstract We propose a globally divergence conforming
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 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 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 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 informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
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 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 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 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 informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
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 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 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 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 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 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 informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
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 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 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 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 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 informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
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 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 informationCubed-Sphere Grids and Galerkin Approaches
2008 Global Atmospheric Core Workshop National Center for Atmospheric Research (NCAR) September 23-24, Boulder CO, USA. Institute for Mathematics Applied to the Geosciences (IMAGe) Cubed-Sphere: Central
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 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 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 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 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 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 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 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 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 informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
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 informationSurprising Computations
.... Surprising Computations Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca www.cs.ubc.ca/ ascher/ Uri Ascher (UBC) Surprising Computations Fall 2012 1 / 67 Motivation.
More informationAdaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations
Adaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations Thida Pongsanguansin, Khamron Mekchay and Montri Maleewong Abstract The adaptive Discontinuous Galerkin DG method for solving
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 informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
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 informationDiscontinuous Galerkin methods for Hyperbolic PDEs. Lecture 1
Discontinuous Galerkin methods for Hyperbolic PDEs. Lecture 1 Olindo Zanotti olindo.zanotti@unitn.it University of Trento Laboratory of Applied Mathematics PiTP 2016 Computational Plasma Astrophysics 18
More informationAngular momentum preserving CFD on general grids
B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
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 informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
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 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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
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 information