Discontinuous Galerkin methods Lecture 2
|
|
- Ethan Rich
- 5 years ago
- Views:
Transcription
1 y y RMMC 2008 Discontinuous Galerkin methods Lecture 2 1 Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu y y Scattering about a vertical cylinder in a finite-width channel cal model based on the linear Padé (2,2) ro tational velocity version we set up a test-0.25 for 0.01 open-channel flow. We seek to model the scattering of an incident wave field which is prop agating toward a bottom-mounted rigid cylin der positioned in the middle of a finite-width x x y x 8.4 Scattering about a vertical cylinder in a finitewidth channel A numerical study is carried out to both test the consistency of the imposed boundary conditions and investigate how the geometric representation of the the domain may 0.25 affect the computed solution. In a numeri- channel Due to the symmetry, the solution is mathe E E E E E E E E E E E E E E E x Figure 8.12: Scattering x of waves about a Darmstadt International Workshop, October 2004 p.79 cylinder in a finite-width channel.
2 A brief overview of what s to come Lecture 1: Introduction, Motivation, History Lecture 2: Basic elements of DG-FEM Lecture 3: Linear systems and some theory Lecture 4: A bit more theory and discrete stability Lecture 5: Attention to implementations Lecture 6: Nonlinear problems and properties Lecture 7: Problems with discontinuities and shocks Lecture 8: Higher order/global problems
3 Lecture 2 Lets briefly recall the issues A bit of notation A second look at the scheme Returning to the example Looking forward
4 Let us recall Local flexibility to achieve high-order and geometric flexibility in the spirit of FEM Explicit scheme and problem control in the spirit of FVM... but many answers remain unanswered How do we achieve high-order accuracy? How do we choose the numerical flux? Is the scheme stable? What is the price?
5 A bit of notation It is the multi-element component of FEM/FVM which gives the geometric flexibility Ω Ω h = K k=1 D k, D k-1 D k D k+1 x 1 l =L x k-1 x K r =xk l x k r =xk+1 l r =R and the solution is approximated as Local norms K u(x, t) u h (x, t) = u k h(x, t), k=1 (u, v) D k = uv dx, u 2 D =(u, u) k D, k D k n inner product and norm
6 A bit of notation Global norms (u, v) Ω,h = K (u, v) D k, u 2 Ω,h =(u, u) Ω,h. k=1 To deal with the solutions across an interface, we define The jump: The average: u [[u]] = ˆn u + ˆn + u +, {{u}} = u + u +, 2 u + = local/interior solution = neighbor/exterior solution outward pointing normal
7 The basics of DG-FEM To understand the role of the different choices, let us consider an example u t + f(u) =0, x [L, R] =Ω, x as ( )=. This is subject to the ap First recall d dt u h 2 Ω,h = a ( u 2 (R) u 2 (L) ), t x f(u) =au. by multiplication of Eq. (2.1) by ( Leading to suitable boundary conditions as u(l, t) =g(t) if a 0, u(r, t) =g(t) if a and energy conservation when nonoverlapping elements u(r) =u(l)
8 The basics of DG-FEM Let us now assume that the local solution is x D k : u k h(x, t) = N p n=1 Where we have Some local modal basis A local nodal basis û k n(t)ψ n (x) = N p i=1 u k h(x k i,t)l k i (x). have introduced two complementary expressions for the lo l k i (x) ψ n (x)
9 The basics of DG-FEM Let us now assume that the local solution is x D k : u k h(x, t) = We form the local residual h and form the local residual x D k : R h (x, t) = uk h t + f k h x and require this to vanish locally in a Galerkin sense N p n=1 Where we have Some local modal basis A local nodal basis û k n(t)ψ n (x) = D k R h (x, t)l k j (x) dx =0, N p i=1 u k h(x k i,t)l k i (x). have introduced two complementary expressions for the lo l k i (x) ψ n (x) ack to the finite element scheme, we
10 The first DG schemes The lack of solution uniqueness at the interface is addressed as in FVM by a numerical flux f = f (u h,u+ h ), it is consistent [i.e u k h dl k j D k t lk j fh k dx glk j dx = [ f l k j g Gauss theorem once again, [ ] [ ] with the corresponding strong form being orem once again, D k R h (x, t)l k j (x) dx = [ (f k h f )l k j ] x k+1 x k. ] x k+1 x k, Naturally, the choice of the flux is important!
11 The basics of DG-FEM To simplify the notation, introduce Weak: M k ij = D k l k i (x)l k j (x) dx, S k ij = Dk l ki (x) dlk j dx dx. to obtain the two basic forms of DG-FEM M k duk h dt Strong: (S k ) T f k h M k g k h = f (x k+1 )l k (x k+1 )+f (x k )l k (x k ) M k duk h dt + S k f k h M k g k h =(f k h (x k+1 ) f (x k+1 ))l k (x k+1 ) (f k h (x k ) f (x k ))l k (x k ).
12 The basics of DG-FEM Let us consider the strong form ( ) M k d ( dt uk h + S k ( ) au k ) [ ] x k h = l k (x)(au k h (au h ) r ). ( ) x k l ced the nodal versions of the local operators ( ) N p N p u T h M k u h = u k h(x k D k i )l k i (x) u k h(x k j )l k j (x) dx = u k h 2 D ; k i=1 j=1 First note that Then also note that hat is, it recovers the local energy. Furthermore, consider u T h S k u h = D k N p i=1 = 1 2 [(uk h) 2 ] xk r x k l u k h(x k i )l k i (x). N p j=1 u k h(x k j ) dlk j dx dx = D k u k h(x)(u k h(x)) dx
13 The basics of DG-FEM Combining this, we achieve from d dt M k d dt uk h + S k ( au k) [ = h l k (x)(au k h (au h ) ) ced the nodal versions of the local operators ( ) the local energy estimate u k h 2 D k h = a[(uk h) 2 ] xk r [ x k l +2 [ u k h(au k h (au) ) ] x k r x k l as for the original equation, require that which should behave in such a way that ] x k r x k l.. K k=1 d dt u k h 2 D k = d dt u h 2 Ω,h 0. Stability
14 The basics of DG-FEM Let us consider the contribution from one interface each of which looks like ˆn (au 2 h(x ) 2u h (x )(au) (x ) ) + ˆn + (au 2 h(x + ) 2u h (x + )(au) (x + ) ) 0 at every interface. Here ( ) refers to the left (e.g., ), and right side We have the freedom ) to ( choose numerical flux to guarantee this. Consider f α = 1 ) ( f =(au) = {{au}} + a 1 α 2 α = 0, - upwind flux pwind α = 1, fl- central/average flux α = 0, w [[u]]. This yields this the gives local a contribution term from each interface a (1 α)[[u h ]] 2, 0, for 0 α 1
15 The basics of DG-FEM What about the boundary conditions? a>0:u(l, t) =g(t) Approach #1: f L = ag(t),f R = au k h (xk r ) Approach #2: f L = au k h (x1 l )+2g(t),f R = au k h (xk r )
16 The basics of DG-FEM Approach #1: K 1 d dt u h 2 Ω,h = a (1 α) [[u k h(x k r)]] 2 (1 α)a(u 1 h(x 1 l )) 2 a(u K h (x K r )) 2. For 0 Approach #2: k=1 1 we have global stability. tflow. This yields a global energy estimate as K 1 d dt u h 2 Ω,h = a (1 α) [[u k h(x k r)]] 2 a(u 1 h(x 1 l )) 2 a(u K h (x K r )) 2. k=1 gain, we recover stability for 0 a 0, 0 α 1 1. A particularly interesting case Stability Energy dissipation for α 1 Energy conservation for α =1,u(L) =u(r)
17 The basics of DG-FEM What did we learn from this? Stability is enforced through the flux choice. No restrictions on the local basis, e.g., it need not be polynomial, and is chosen to provide accuracy. The numerical solution is discontinuous between elements. Boundary conditions and interface conditions are imposed weakly. All operators are local. Due to the weak interface based coupling, there are no restrictions on element size and local approximation.
18 Back to the example Let us consider a simple example u u 2π =0, x [0, 2π], t x u(x, 0) = sin(lx), l = 2π λ, ons and initial condition as N\ K Convergence rate 1 4.0E E E E E E E E E E E E E E E E E E E E E E E The error clearly behaves as u u h Ω,h Ch N+1. proximation that gives the
19 Back to the example What about time dependence Final time (T) π 10π 100π 1000π 2000π (N,K)=(2,4) 4.3E E E-01 >1 >1 (N,K)=(4,2) 3.3E E E E E-01 (N,K)=(4,4) 3.1E E E E E-03 The error behaves as u u h Ω,h C(T )h N+1 (c 1 + c 2 T )h N+1,
20 Back to the example u(x,2π) x/2π u(x,2π) x/2π Central flux Upwind flux
21 Back to the example What about cost? N\K dominates and destroys the expected convergence rate. Higher order is cheaper Time C(T )K(N + 1) 2, N\K Convergence rate with the final time in a linear 1 4.0E E E E E E E E E E E E E E E E E E E E E E E
22 A few remarks We know now It is the flux that gives stability It is the local basis that gives accuracy. The scheme is very (VERY) flexible.
23 A few remarks We know now It is the flux that gives stability It is the local basis that gives accuracy. The scheme is very (VERY) flexible. BUT - we have doubled the number of degrees of freedom along the interfaces In 1D not a big deal -- penalty is N+1 N
24 Does it generalize? Let us first consider the scalar conservation law u t + f(u) x =0, x [L, R], where the e initial boundary conditions conditions are u(l, t) =g 1 (t) when f u (u(l, t)) 0, u(r, t) =g 2 (t) when f u (u(r, t)) 0. tly as for the linear case discussed previousl Assume as usual that x D k : u k h(x, t) = N p n=1 û k n(t)ψ n (x) = N p i=1 u k h(x k i,t)l k i (x). have introduced two complementary expressions for the lo
25 Does it generalize? From this we directly recover ( ) u k h t φk h fh k (u k h) dφk h dx = ( dx ) D k D k ˆn f φ k h dx, strong form and the ( corresponding ) strong form D k ( u k h t + f h k(uk h ) ) φ k h dx = x D k ˆn (f k h (u k h) f ) φ k h dx, Here we have the general local test functions φ k h Vk h. x D k : φ k h(x) = N p ˆφ k nψ n (x), n=1 tions as
26 Does it generalize? This yields exactly local unknowns N p D k N p equations for the and the corresponding strong form D k Here we have also introduced ( u k h t ψ n fh k (u k h) dψ ) n ( dx ) = dx ( )) ( u k h t + f h k(uk h ) ) ψ n dx = x D k ˆn f ψ n dx, ( ( ) ) D k ˆn (f k h (u k h) f ) ψ n dx, the weak and strong semi-discrete form, respectively. In both c x D k : f k h (u k h)= N p n=1 ˆf k nψ n (x) = N p i=1 The details of this will return! f h (x k i )l k i (x).
27 Does it generalize? The only thing that remains unknown is the flux f = f (u h,u+ h ) We rely on the hugely successful theory of finite volume monotone schemes
28 Does it generalize? The only thing that remains unknown is the flux f = f (u h,u+ h ) We rely on the hugely successful theory of finite volume monotone schemes The flux is consistent The flux is monotone f(u h )=f (u h,u h ) increasing in a f (a, b){ decreasing in b
29 Does it generalize? There are many choices for this The Lax-Friedrich flux where the global LF flux is given by and the local LF flux is obtained for NOTE: for f LF (a, b) = f(a)+f(b) 2 choice. The definition of C C + C 2 ˆn (a b), allows for some v max f u(s) inf u h (x) s sup u h (x) f = au max f u(s). min(a,b) s max(a,b) we recover the flux seen earlier
30 Does it generalize?... but the FV literature is filled with alternatives Exact Riemann solvers Godunov fluxes Engquist-Osher fluxes Approximate Riemann fluxes (Roe, van Leer, HLLC etc)
31 Does it generalize?... but the FV literature is filled with alternatives Exact Riemann solvers Godunov fluxes Engquist-Osher fluxes Approximate Riemann fluxes (Roe, van Leer, HLLC etc) To keep things simple we shall mainly focus on the LF flux which generally works very well, but is also the most dissipative flux.
32 Does it generalize? Let us now consider the system of conservation laws u t + f(u) =0, x [L, R], x omponent is introduced to o s u =[u 1 (x, t),...,u m (x, t)] T m m. where the boundary conditions are B L u(l, t) =g 1 (t) at x = L, B R u(r, t) =g 2 (t) at x = R, B The only essential difference is that C in the LF flux depends on the eigenvalues of f u = f u.
33 Does it generalize? The last frontier -- multidimensional problem No essential difference - the weak form ( u k ( ) ( h t φk h f k h(u k h) φ k ) h dx = D k The strong form D k ˆn f φ k h dx, the strong form ( ) ( ) ) u k ( h + f k D k t h(u k h) φ k h dx = ˆn f k h(u k h) f ) φ k D k h dx. with the LF-flux f = {{f h (u h )}} + C 2 [[u h]]. C = max u ( λ ˆn f ), u eigenvalue of the matrix.
34 Lets summarize We already know a lot about the basic DG-FEM Stability is provided by carefully choosing the numerical flux. Accuracy appear to be given by the local solution representation. We can utilize major advances on monotone schemes to design fluxes. The scheme generalizes with very few changes to very general problems -- multidimensional systems of conservation laws.
35 Lets summarize We already know a lot about the basic DG-FEM Stability is provided by carefully choosing the numerical flux. Accuracy appear to be given by the local solution representation. We can utilize major advances on monotone schemes to design fluxes. The scheme generalizes with very few changes to very general problems -- multidimensional systems of conservation laws. At least in principle -- but what can we actually prove?
Lecture - 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 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 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 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 informationResidual Distribution. basics, recents developments, relations with other techniques
Residual Distribution basics, recents developments, relations with other techniques MARIO RICCHIUTO July 19, 2012 INTRODUCTION with historical perspective t u+ F(u) = 0 In the 80 s numerical techniques
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 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 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 informationDG-FEM for PDE s Lecture 4
Univ Coruna/USC/Univ Vigo 2010 DG-FEM for PDE s Lecture 4 Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu A brief overview of what s to come Lecture 1: Introduction and DG in 1D Lecture 2: Implementation
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 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 informationApplication of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling
1 Application of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling Xin Wang ABSTRACT This work will explore the discontinuous Galerkin finite element method (DG-FEM) for solving acoustic wave
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 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 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 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 informationFinite Element methods for hyperbolic systems
Finite Element methods for hyperbolic systems Eric Sonnendrücker Max-Planck-Institut für Plasmaphysik und Zentrum Mathematik, TU München Lecture notes Wintersemester 14-15 January 19, 15 Contents 1 Introduction
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 informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationC e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a
C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a Modelling, Analysis and Simulation Modelling, Analysis and Simulation Bilinear forms for the recovery-based discontinuous Galerkin method
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
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 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 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 information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
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 informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationDiscontinuous Galerkin Method for 1D Shallow Water Flow with Water Surface Slope Limiter
Clemson University TigerPrints All Theses Theses 1-010 Discontinuous Galerkin Method for 1D Shallow Water Flow with Water Surface Slope Limiter Wencong Lai Clemson University, wlai@clemson.edu Follow this
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 informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
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 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 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 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 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 informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
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 informationInvestigation of Godunov Flux Against Lax Friedrichs' Flux for the RKDG Methods on the Scalar Nonlinear Conservation Laws Using Smoothness Indicator
American Review of Mathematics and Statistics December 2014, Vol. 2, No. 2, pp. 43-53 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s). 2014. All Rights Reserved. Published by American
More informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationFirst order BSSN formulation of Einstein s field equations
David Brown 1 Peter Diener 2 3 Jan Hesthaven 4 Frank Herrmann 3 Abdul Mroué 5 Olivier Sarbach 6 Erik Schnetter 7 Manuel Tiglio 3 Michael Wagman 4 1 North Carolina State University 2 Louisiana State University
More informationRECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES
RECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES FOR INVISCID AND VISCOUS PROBLEMS. R. Abgrall Team Bacchus INRIA Bordeaux Sud Ouest and UniversitédeBordeauxTalence,France Roscoff,
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 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 informationDistances, volumes, and integration
Distances, volumes, and integration Scribe: Aric Bartle 1 Local Shape of a Surface A question that we may ask ourselves is what significance does the second fundamental form play in the geometric characteristics
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 informationComparison of Approximate Riemann Solvers
Comparison of Approximate Riemann Solvers Charlotte Kong May 0 Department of Mathematics University of Reading Supervisor: Dr P Sweby A dissertation submitted in partial fulfilment of the requirement for
More informationNodal Discontinuous Galerkin Methods:
Jan S Hesthaven and Tim Warburton Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications LIST OF CORRECTIONS AND CLARIFICATIONS September 6, 2010 Springer List of corrections and
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 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 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 informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 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 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 information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationUncertainty Quantification in Computational Science
DTU 2010 - Lecture I Uncertainty Quantification in Computational Science Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu Objective of lectures The main objective of these lectures are To offer
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 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 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 informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationSuperconvergence of the Direct Discontinuous Galerkin Method for Convection-Diffusion Equations
Superconvergence of the Direct Discontinuous Galerkin Method for Convection-Diffusion Equations Waixiang Cao, Hailiang Liu, Zhimin Zhang,3 Division of Applied and Computational Mathematics, Beijing Computational
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 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 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 informationNodal Discontinuous Galerkin Methods:
Jan S Hesthaven and Tim Warburton Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications LIST OF CORRECTIONS AND CLARIFICATIONS August 12, 2009 Springer List of corrections and clarifications
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
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 informationOn the entropy stability of Roe-type finite volume methods
Proceedings of Symposia in Applied Mathematics On the entropy stability of Roe-type finite volume methods Mária Lukáčová - Medvid ová and Eitan Tadmor Abstract. We study the entropy stability of a class
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
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 informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationSolution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers
Solution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers Alexander Kurganov, 1, * Eitan Tadmor 2 1 Department of Mathematics, University of Michigan, Ann Arbor, Michigan
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 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 informationElliptic Kirchhoff equations
Elliptic Kirchhoff equations David ARCOYA Universidad de Granada Sevilla, 8-IX-2015 Workshop on Recent Advances in PDEs: Analysis, Numerics and Control In honor of Enrique Fernández-Cara for his 60th birthday
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 informationA DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR TIME DEPENDENT PARTIAL DIFFERENTIAL EQUATIONS WITH HIGHER ORDER DERIVATIVES
MATHEMATCS OF COMPUTATON Volume 77, Number 6, April 008, Pages 699 730 S 005-5718(07)0045-5 Article electronically published on September 6, 007 A DSCONTNUOUS GALERKN FNTE ELEMENT METHOD FOR TME DEPENDENT
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
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 informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationRunge Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems
Journal of Scientific Computing, Vol. 16, No. 3, September 2001 ( 2001) Review Article Runge Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems Bernardo Cockburn 1 and Chi-Wang Shu
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
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 informationFROM GODUNOV TO A UNIFIED HYBRIDIZED DISCONTINUOUS GALERKIN FRAMEWORK FOR PARTIAL DIFFERENTIAL EQUATIONS
FROM GODUNOV TO A UNIFIED HYBRIDIZED DISCONTINUOUS GALERKIN FRAMEWORK FOR PARTIAL DIFFERENTIAL EQUATIONS TAN BUI-THANH Abstract. By revisiting the basic Godunov approach for linear system of hyperbolic
More information