Approximation des systemes hyperboliques par elements finis continus non uniformes en dimension quelconque
|
|
- Alicia Craig
- 6 years ago
- Views:
Transcription
1 Approximation des systemes hyperboliques par elements finis continus non uniformes en dimension quelconque Jean-Luc Guermond and Bojan Popov Department of Mathematics Texas A&M University Séminaire du Laboratoire Jacques-Louis Lions UPMC 04 Nov 2016
2 Acknowledgments Collaborators: Murtazo Nazarov (Uppsala University) Vladimir Tomov (LLNL) Young Yang (Penn State University) Laura Saavedra (Universidad Politécnica de Madrid) Support:
3 Hyperbolic systems 1 Hyperbolic systems 2 FE approximation 3 Hyperbolic systems + ALE 4 Maximum wave speed Hyperbolic systems
4 Hyperbolic systems The PDEs Hyperbolic tu + r f(u) =0, (x, t) 2 D R +. u(x, 0) = u 0 (x), x 2 D. D open polyhedral domain in R d. f 2C 1 (R m ; R m d ), the flux. u 0,admissibleinitialdata. Periodic BCs or u 0 has compact support (to simplify BCs)
5 Hyperbolic systems Examples Linear transport: f(u) = Burgers equation: f (u) = 1 2 u2 Tra c flow equation: f (u) =v max(1 Bucley-Leverett: f (u) = Shallow water: p-system Euler equations... u = 2 R q 1+d, f(u) = u u u max )u u 2 u 2 +a(1 u) 2, a is constant 1! q T 1 q q R 2 g 2 I (1+d) d d
6 Formulation of the problem Assumptions 9 admissible set A s.t. for all (u l, u r ) 2Athe 1D Riemann problem ( u l if x < tv x (n f(v)) = 0, v(x, 0) = u r if x > 0. has a unique entropy solution u(u l, u r )(x, t) foralln 2 R d, knk`2 =1. There exists an invariant set A A,i.e., A is convex. u(u l, u r )(x, t) 2 A, 8t 0, 8x 2 R, 8u l, u r 2 A.
7 Formulation of the problem Examples of invariant sets Invariant domains are convex for genuinely nonlinear systems (Ho (1979, 1985), Chueh, Conley, Smoller (1973)). Scalar conservation in R d : A =[a, b], 8a apple b 2 R. Euler: A = { >0, e > 0, s a}, 8a 2 R,wheres is the specific entropy. p-system (1D): etc. U =(v, u) T A := {U 2 R + R a apple W 2 (U) apple W 1 (U) apple b}, 8a apple b 2 R where W 1 and W 2 are the Riemann invariants Z 1 p W 1 (U) =u + p 0 (s) ds, and W 2 (U) =u v Z 1 p v p 0 (s) ds.
8 FE approximation 1 Hyperbolic systems 2 FE approximation 3 Hyperbolic systems + ALE 4 Maximum wave speed Hyperbolic systems
9 Approximation (time and space) FE space/shape functions {T h } h>0 shape regular conforming mesh sequence {' 1,...,' I },positive+partitionofunity( P j2{1:i } ' j =1) Ex: P 1, Q 1,Bernsteinpolynomials(anydegree) m i := R D ' i dx, lumpedmassmatrix(m i = P j2i(s i ) R D ' i ' j dx)
10 Approximation (time and space) Algorithm: Galerkin Set u h (x, t) = P I j=1 U j (t)' j (x). Galerkin + lumped mass matrix m Z + r (f(u h ))' i dx =0 D Algorithm: Galerkin + First-order viscosity + Explicit Euler by Un+1 i t Approximate f(u h )by P j2i(s i ) (f(un j ))' j U n+1 i m i t U n i Z + D U n i 0 1 X (f(u n j ))' j A ' i dx + X dij n (Un i U n j ) =0. j2i(s i ) j2i(s i ) How should we choose artificial viscosity d n ij?
11 Approximation (time and space) Algorithm: Galerkin + First-order viscosity + Explicit Euler Introduce Then U n+1 i m i Z c ij = ' i (x)r' j (x) dx. D t U n i = X j c ij f(u j )+d n ij U j Observe that conservation implies P j c ij =0,(partitionofunity) We define d n ii such that P j d n ij =0,(conservation).. Remark Rest of the talk applies to any method that can be formalized as above. (FV, DG, FD, etc.)
12 Approximation (time and space) Algorithm: Galerkin + First-order viscosity + Explicit Euler Observe that conservation implies P j c ij =0and P j d n ij =0. U n+1 i m i t U n i = X j c ij (f(u i ) f(u j )) + dij n (U i + U j ). Try to construct convex combination... U n+1 i = U n i (1 + 2 t m i D ii )+ X j6=i = U n i (1 X j6=i 2 t m i d n ij )+X j6=i t c ij (f(u i ) f(u j )) + dij n m (U i + U j ) i 2 t dij n m i Introduce intermediate states U(U i, U j ) 1 2 (U i + U j )+ c! ij 2dij n (f(u i ) f(u j )) U(U i, U j ):= 1 2 (U i + U j )+ c ij 2dij n (f(u i ) f(u j )).
13 Approximation (time and space) Algorithm: Galerkin + First-order viscosity + Explicit Euler Now construct convex combination U n+1 i = U n i (1 X 2 t dij n m )+X j6=i i j6=i Are the states U(U i, U j )goodobjects? 2 t m i d n ij U(U i, U j )
14 Approximation (time and space) Algorithm: Galerkin + First-order viscosity + Explicit Euler Define n ij = c ij /kc ij k`2 2 R d,(unitvector). f ij (U) :=n ij f(u) isanhyperbolicfluxbydefinitionofhyperbolicity! Then U(U i, U j ):= 1 2 (U i + U j )+ kc ij k`2 2dij n (f ij (U i ) f ij (U j )).
15 Approximation (time and space) Lemma (GP (2015)) Consider the fake 1D Riemann tv x (n ij f(v)) = 0, v(x, 0) = ( U i if x < 0 U j if x > 0. Let max(f, n ij, U i, U j ) be maximum wave speed in 1D Riemann problem Then U(U i, U j )= R v(x, t) dx with fake time t = kc ij k`2 2d ij n,provided kc ij k`2 2d n ij max(f, n ij, U i, U j )= t max(f, n ij, U i, U j ) apple 1 2 Define viscosity coe cient d n ij := max(f, n ij, U i, U j )kc ij k`2, j 6= i.
16 Approximation (time and space) Theorem (GP (2015)) Provided CFL condition, (1 2 t m i D ii ) 0. Local invariance: U n+1 i 2 Conv{U(U n i, Un j ) j 2I(S i )}. Global invariance: The scheme preserves all the convex invariant sets. (Let A be a convex invariant set, assume U 0 2 A, thenu n+1 i 2 A for all n 0.) Discrete entropy inequality for all the entropy pairs (, q): Z m i t ( (Un+1 i ) (U n i )) + D r ( h q(u n h ))' i dx + X i6=j2i(s i ) d ij (U n j ) apple 0.
17 Approximation (time and space) Is it new? Loose extension of non-staggered Lax-Friedrichs to FE. Similar results proved by Ho (1979, 1985), Perthame-Shu (1996), Frid (2001) in FV context and compressible Euler. Some relation with flux vector splitting theory of Bouchut-Frid (2006). Not aware of similar results for arbitrary hyperbolic systems and continuous FE.
18 A priori error estimate for scalar equations: Definition of mollifiers Let >0and = kfk Lip Consider mollifiers! and! 8 1 >< t apple, 3 2 t! (t) := 3 >: 2 apple t apple2, 0 otherwise,! (x) := d l=1! (x l ), x := (x 1,...,x d ). Following Kruskov (1970), define (x, y, t, s) :=! (x y)! (t s), 8(y, s) 2 D [0, T ]. Following Cockburn-Gremaud (1996,1998), define Z t (t) :=! (s) ds. 0
19 A priori error estimate for scalar equations: A useful lemma Lemma (Guermond, Popov ( )) Assume u 0 2 BV ( ). Leteu h : D [0, T ]! R be any approximate solution. Assume that there is aboundedfunctionalonlipschitzfunctionssothat8k 2 [u min, u max], 8 2 Wc 1,1 (D [0, T ]; R + ): Z T Z eu h t + sgn(eu h k)(f(eu h ) f(k)) r dx dt 0 D + k h (eu h (T ) k) h (, T h ) k`1 k h (eu h (0) k) h (, h h) k`1 apple ( ), h where k k`1 h is the discrete L 1 -norm and T T h apple t, 0 h apple t, >0 is a uniform constant. Then the following estimate holds ku(, T ) ũ h (, T )k L 1 ( ) apple c (( + h) u 0 BV ( ) + ) R t0 R where D ( ) dy ds := sup 0appletappleT. (t) Generalization of results by Cockburn-Gremaud (1996) and Bouchut-Perthame (1998) based on Kruskov (1970), Kuznecov (1976).
20 A priori error estimate for scalar equations English translation Control on all the Kruskov entropies ) Convergence estimate. Theorem (Guermond, Popov ( )) Assume u 0 2 BV and f Lipschitz. Let u h be the first-order viscosity solution. Then there is c 0,uniform,suchthatthefollowingholdsifCFLapple c 0 : (i) ku(t ) u h (T )k L 1 ((0,T );L 1 ) apple ch 1 2 if a priori BV estimate on u h. (ii) ku(t ) u h (T )k L 1 ((0,T );L 1 ) apple ch 1 4 otherwise. BV estimate is trivial in 1D (Harten s lemma). BV estimate can be proved in nd on special meshes. Similar results for FV Chainais-Hillairet (1999), Eymard et al (1998) First error estimates for explicit continuous FE method (as far as we know).
21 High-order extension Higher-order in time Use SSP method to get higher-order in time. Strong Stability Preserving methods (SSP), Kraaijevanger (1991) (amazing paper), Gottlieb-Shu-Tadmor (2001), Spiteri-Ruuth (2002) Ferracina-Spijker (2005), Higueras (2005), etc.:
22 High-order extension Higher-order in time The midpoint rule is not SSP Heun s method is SSP w (1) = u n + tl(t n, u n ) w (2) = w (1) + tl(t n + t, w (1) ) u n+1 = 1 2 un w (2). SSPRK(3, 3): w (1) = u n + tl(t n, u n ), z (1) = w (1) + tl(t n + t, w (1) ), w (2) = 3 4 un z(1), z (2) = w (2) + tl(t n t, w (2) ), u n+1 = 1 3 un z(2).
23 High-order extension Remark on SSP SSP is not about positivity, it is about convexity. Let U! S t (U) besspschemebasedoneulerstepu! E t (U), t apple t 0, Let A be a convex set then If Euler step E t (U) isinvariantdomainpreservingina then SSP step S t (U) isinvariantdomainpreservingina SSP methods preserve convex domains that are invariant for forward Euler time stepping (It s all about convexity).
24 High-order extension Higher-order in space: Entropy viscosity Use entropy viscosity (or something else) FCT or other limitation (work in progress)
25 Strong explosion; ent. vis. sol. 1.5 million P 2 nodes (author: Murtazo Nazarov; 1.5 million P 2 nodes)
26 Mach 10 ramp, ent. vis. sol. 1.2 million P 2 nodes (author: Murtazo Nazarov; 1.2 millions P 2 nodes)
27 Mach 3, Step, ent. vis. sol P 1 nodes and viscous sol P 1 nodes
28 Hyperbolic systems + ALE 1 Hyperbolic systems 2 FE approximation 3 Hyperbolic systems + ALE 4 Maximum wave speed Hyperbolic systems
29 ALE formulation Instead of tracking the characteristics (there are too many), we want to move the mesh. ALE formulation Let : R d R +! R d be a uniformly Lipschitz mapping (R d 3 7! (, t) 2 R d invertible on [0, t ]) Let v A (x, t) =@ t ( 1 t (x), t) Arbitrary Lagrangian Eulerian velocity We are going to use v A to move the mesh. Lemma The following holds in the distribution sense (in time) over [0, t ] for every function 2 C0 0(Rd 1 ; R) (with the notation '(x, t) := ( t (x))): t u(x, t)'(x, t) dx = r (u v A f(u))'(x, t) dx. ZR d R d
30 Finite elements Geometric Finite elements Let (T 0 h ) h>0 be a shape-regular sequence of matching meshes. Reference Lagrange finite element ( b K, b P geo, b geo )forgeometry Lagrange nodes {ba i } i2{1:n geo sh } and Lagrange shape functions {b geo i } geo i2{1:n sh } {a n i } i2{1:i geo } collection of all the Lagrange nodes in the mesh Th n j geo : Th n {1:ngeo sh }!{1:I geo } geometric connectivity array Geometric transformation T n K : b K T n K (bx) =! K defined by X i2{1:n geo sh n a } b j geo (i,k) geo i (bx). ) Mesh motion controlled by motion of Lagrange nodes
31 Finite elements Approximating Finite elements Reference finite element ( b K, b P, b )} Shape functions b i (x) 0, P i2{1:n sh } b i (bx) =1 Finite element spaces P(T n h ):={v 2C0 (D n ; R); v K P d (Th n ):=[P(T h n )]d, P m(th n ):=[P(T h n )]m. { n i } i2{1:i } global shape functions in P(T n h ). T n K 2 b P, 8K 2T n h },
32 Finite elements The algorithm Initialization: m 0 i := R R d n i (x) dx u h0 := P i2{1:i } U0 i 0 i 2 P m(t 0 h ) ALE velocity field given: w n = P i2{1:i } Wn i n i Mesh motion: a n+1 i = a n i + tw n (a n i ). Mass update: (do not use m n+1 i = R n+1 D i dx!) Z m n+1 i = m n i + t Update approximation field u n+1 h S n i 2 P d (T n h ), n i (x)r wn (x) dx. m n+1 i U n+1 i t m n i Un i X j2i(s i n ) Z + r R d d n ij Un j X j2{1:i } (f(u n j ) Un j W n j ) n j (x) n i (x) dx =0,
33 Finite elements Definition of d n ij Consider flux g n j (v) :=f(v) v Wn j, j 2{1:I } Consider one-dimensional Riemann tv x (g n j (v) nn ij )=0, (x, t) 2 R R+, v(x, 0) = ( U n i if x < 0 U n j if x > 0. Define d n ij by dij n =max( max(gj n, nn ij, Un i, Un j )kcn ij k`2, max(gn i, nn ji, Un j, Un i )kcn ji k`2 ). Note that max(g n j, nn ij, Un i, Un j )=max( L(f, n n ij, Un i, Un j ) Wn j nn ij, R(f, n n ij, Un i, Un j ) Wn j nn ij ).
34 Conservation and invariant domain property Theorem (GPSY (2015)) The total mass P i2{1:i } mn i Un i is conserved. Provided CFL condition, (1 2 t m n D ii ) 0. i Local invariance: U n+1 i 2 Conv{U(U n i, Un j ) j 2I(S i )}. Global invariance. Let A be a convex invariant set, assume U 0 2 A, then U n+1 i 2 A for all n 0. The scheme preserves all the convex invariant sets. Discrete entropy inequality for any entropy pair (, q) 1t mn+1 i (U n+1 i ) m n i (Un i ) apple X Z r R d j2i(s i n ) X j2i(s i n ) (q(u n j ) d n ij (Un j ) (Un j )Wn j ) n j (x) n i (x) dx Corollary (GPSY (2015)) The scheme preserves constant states (Discrete Global Conservation Law (DGCL))
35 Hyperbolic systems FE approximation Hyperbolic systems + ALE Maximum wave speed 2D u + r ( 12 u 2 ) = 0, u0 (x) = 1S, with := (1, 1)T, S := (0, 1)2 Figure: Burgers equation, mesh. Left: Q1 FEM with 25 contours; Center left: Final Q1 mesh; Center right: P1 FEM with 25 contours; Right: Final P1 mesh.
36 Hyperbolic systems FE approximation Hyperbolic systems + ALE Maximum wave speed Nonconvex flux (KPP u + r f(u) = 0, u0 (x) = kxk `2 < , with f(u) = (sin u, cos u)t Figure: KPP problem, mesh. Left: Q1 FEM with 25 contours; Center left: Final Q1 mesh; Center right: P1 FEM with 25 contours; Right: Final P1 mesh.
37 Euler Compressible Euler, 2D Noh problem, = 5 3 Initial data x 0 (x) =1.0, u 0 (x) = 1 x6=0, p 0 (x) = kxk`2 Q 1 P 1 #dofs L 2 -norm L 1 -norm L 2 -norm L 1 -norm E E E E E E E E E E E Table: Noh problem, convergence test, T = 0.6, CFL = 0.2
38 Compressible Euler, 2D Noh problem, = 5 3 Figure: Noh problem at t =0.6, mesh. From left to right: density field with Q 1 approximation (25 contour lines); mesh with Q 1 approximation; density field with P 1 approximation (25 contour lines); mesh with P 1 approximation.
39 Compressible Euler, 3D Noh problem, = 5 3 Figure: Density cuts for the 3D Noh problem at t =0.6. Figure: 3D Noh problem at t = MPI tasks division.
40 Maximum wave speed 1 Hyperbolic systems 2 FE approximation 3 Hyperbolic systems + ALE 4 Maximum wave speed Hyperbolic systems
41 How to compute local viscosity? d n ij := 2 max(f, n ij, U i, U j )kc ij k`2,forj 6= i. max(f, n ij, U i, U j )ismaxwavespeedforriemannproblem
42 Riemann fan for Euler, p =( 1) e Structure of the Riemann problem (Lax (1957), Bressan (2000), Toro (2009)). Waves 1 and 3 are genuinely nonlinear (either shock or rarefaction) Wave 2 is linearly degenerate (contact) w L =( L, u L, p L ), wl =( L, u, p ), wr =( R, u, p ), w R =( R, u R, p R ),
43 Maximum wave speed bound Euler system, p =( 1) e Given the states U L and U R,wehave 1 = u L a L 1+ (p 1 p L ) < 3 = u R +a R 1+ (p 1 p R ) p L 2 p R 2 where p is the pressure of the intermediate state. Then and define max(u L, U R )=max( 1, 3 ). In practice we just need a good upper bound of p : p p.then max(u L, U R )=max( 1(p ), 3(p ) ).
44 Maximum wave speed bound To avoid computing p,itisacommonpracticetoestimate max( u L + a L, u R + a R ) This estimate is inaccurate and can be wrong. max by
45 Maximum wave speed bound Counter-example 1: 1-wave and the 3-wave are both shocks Toro 2009, L R u L u R p L p R max but max( u L + a L, u R + a R ) 29.97, large overestimation
46 Maximum wave speed bound Counter-example 2: 1-wave is a shock and the 3-wave is an expansion L R u L u R p L p R max but max( u L + a L, u R + a R ) 1.183, large underestimation
47 Definition of p Let p be the zero of R, then 0 1 p B a 1 L + a R 2 (u R u L ) C 1 1 A 2 2 a L pl + a R pr 2 1 Lemma (GP (2016)) We have p < p in the physical range of, 1 < apple 5 3. p is an upper bound on p. min(p L, p R ) apple p apple p (starting guess for cubic Newton alg., GP (2016))
48 Conclusions Continuous finite elements Continuous FE are viable tools to solve hyperbolic systems. Continuous FE are viable alternatives to DG and FV. Continuous FE are easy to implement and parallelize. Exa-scale computing will need simple, robust, methods. Current and future work Convergence analysis, error estimates beyond first-order. Extension to DG. Extension of BBZ to higher-order polynomials (order 3 and higher). Extension of BBZ to systems (Shallow water, Euler). Extension to equations with source terms (Radiative transport, Radiative hydrodynamics).
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
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 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 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 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 informationHyperbolic Problems: Theory, Numerics, Applications
AIMS on Applied Mathematics Vol.8 Hyperbolic Problems: Theory, Numerics, Applications Fabio Ancona Alberto Bressan Pierangelo Marcati Andrea Marson Editors American Institute of Mathematical Sciences AIMS
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 informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
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 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 Viscosity and L1-based Approximations of PDEs: Exploiting Sparsity
https://livelink.ebs.afrl.af.mil/livelink/llisapi.dll Page 1 of 2 10/28/2015 AFRL-AFOSR-VA-TR-2015-0337 Entropy Viscosity and L1-based Approximations of PDEs: Exploiting Sparsity Jean-Luc Guermond TEXAS
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
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 information1. Introduction. We are interested in the scalar hyperbolic conservation law
SIAM J. NUMER. ANAL. Vol. 4, No. 5, pp. 1978 1997 c 005 Society for Industrial and Applied Mathematics ON CONVERGENCE OF MINMOD-TYPE SCHEMES SERGEI KONYAGIN, BOJAN POPOV, AND OGNIAN TRIFONOV Abstract.
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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
More informationOn Convergence of Minmod-Type Schemes
On Convergence of Minmod-Type Schemes Sergei Konyagin, Boan Popov, and Ognian Trifonov April 16, 2004 Abstract A class of non-oscillatory numerical methods for solving nonlinear scalar conservation laws
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 informationDo Finite Volume Methods Need a Mesh?
Do Finite Volume Methods Need a Mesh? Michael Junk Fachbereich Mathematik, Universität Kaiserslautern, 67663 Kaiserslautern, Germany Abstract. In this article, finite volume discretizations of hyperbolic
More informationOn High Order Strong Stability Preserving Runge Kutta and Multi Step Time Discretizations
Journal of Scientific Computing, Vol. 5, Nos. /, November 005 ( 005) DOI: 0.007/s095-00-65-5 On High Order Strong Stability Preserving Runge Kutta and Multi Step Time Discretizations Sigal Gottlieb Received
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 informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More 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 informationModeling uncertainties with kinetic equations
B. Després (LJLL- University Paris VI) thanks to G. Poette (CEA) and D. Lucor (Orsay) B. Perthame (LJLL) E. Trélat (LJLL) Modeling uncertainties with kinetic equations B. Després (LJLL-University Paris
More informationARTICLE IN PRESS Mathematical and Computer Modelling ( )
Mathematical and Computer Modelling Contents lists available at ScienceDirect Mathematical and Computer Modelling ournal homepage: wwwelseviercom/locate/mcm Total variation diminishing nonstandard finite
More informationEntropy Viscosity Method for the Single Material Euler Equations in Lagrangian Frame
Entropy Viscosity Method for the Single Material Euler Equations in Lagrangian Frame Jean-Luc Guermond a, Bojan Popov a, Vladimir Tomov b, a Department of Mathematics, Texas A&M University 3368 TAMU, College
More informationStrong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators
Journal of Scientific Computing, Vol. 8, No., February 3 ( 3) Strong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators Sigal Gottlieb
More informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
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 informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
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 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 informationWEAK ASYMPTOTIC SOLUTION FOR A NON-STRICTLY HYPERBOLIC SYSTEM OF CONSERVATION LAWS-II
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 94, pp. 1 14. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu WEAK ASYMPTOTIC
More informationGlobal Existence of Large BV Solutions in a Model of Granular Flow
This article was downloaded by: [Pennsylvania State University] On: 08 February 2012, At: 09:55 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
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 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 informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationThe inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method
The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas
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 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 informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationLARGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMERICS FOR 1-D SCALAR CONSERVATION LAWS
LAGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMEICS FO 1-D SCALA CONSEVATION LAWS L. I. IGNAT, A. POZO, E. ZUAZUA Abstract. In this paper we analyze the large time asymptotic behavior of the discrete
More informationFinite element quasi-interpolation, best approximation and applications
Finite element quasi-interpolation, best approximation and applications + some recent results in MHD and hyperbolic systems ICES, UT Austin, Feb 4 2016 Acknowledgments Collaborators: Alexandre Ern (University
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 informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationMath 5588 Final Exam Solutions
Math 5588 Final Exam Solutions Prof. Jeff Calder May 9, 2017 1. Find the function u : [0, 1] R that minimizes I(u) = subject to u(0) = 0 and u(1) = 1. 1 0 e u(x) u (x) + u (x) 2 dx, Solution. Since the
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 informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
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 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 informationM. HERTY, CH. JÖRRES, AND B. PICCOLI
EXISTENCE OF SOLUTION TO SUPPLY CHAIN MODELS BASED ON PARTIAL DIFFERENTIAL EQUATION WITH DISCONTINUOUS FLUX FUNCTION M. HERTY, CH. JÖRRES, AND B. PICCOLI Abstract. We consider a recently [2] proposed model
More informationLecture Notes on Hyperbolic Conservation Laws
Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide
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 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 informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationNumerical methods for conservation laws with a stochastically driven flux
Numerical methods for conservation laws with a stochastically driven flux Håkon Hoel, Kenneth Karlsen, Nils Henrik Risebro, Erlend Briseid Storrøsten Department of Mathematics, University of Oslo, Norway
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 informationEntropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing
Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Siddhartha Mishra Seminar for Applied Mathematics ETH Zurich Goal Find a numerical scheme for conservation
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 informationOscillating waves and optimal smoothing effect for one-dimensional nonlinear scalar conservation laws
arxiv:1302.1345v1 [math.ap] 6 Feb 2013 Oscillating waves and optimal smoothing effect for one-dimensional nonlinear scalar conservation laws Pierre Castelli and Stéphane Junca February 7, 2013 Pierre Castelli
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationStrong Stability of Singly-Diagonally-Implicit Runge-Kutta Methods
Strong Stability of Singly-Diagonally-Implicit Runge-Kutta Methods L. Ferracina and M. N. Spijker 2007, June 4 Abstract. This paper deals with the numerical solution of initial value problems, for systems
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
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 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 informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationNonlinear Regularizing Effects for Hyperbolic Conservation Laws
Nonlinear Regularizing Effects for Hyperbolic Conservation Laws Ecole Polytechnique Centre de Mathématiques Laurent Schwartz Collège de France, séminaire EDP, 4 mai 2012 Motivation Consider Cauchy problem
More informationOn universality of critical behaviour in Hamiltonian PDEs
Riemann - Hilbert Problems, Integrability and Asymptotics Trieste, September 23, 2005 On universality of critical behaviour in Hamiltonian PDEs Boris DUBROVIN SISSA (Trieste) 1 Main subject: Hamiltonian
More informationA model for a network of conveyor belts with discontinuous speed and capacity
A model for a network of conveyor belts with discontinuous speed and capacity Adriano FESTA Seminario di Modellistica differenziale Numerica - 6.03.2018 work in collaboration with M. Pfirsching, S. Goettlich
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 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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationDownloaded 10/11/16 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 54, No. 4, pp. 2466 2489 c 2016 Society for Industrial and Applied Matematics INVARIANT DOMAINS AND FIRST-ORDER CONTINUOUS FINITE ELEMENT APPROXIMATION FOR HYPERBOLIC SYSTEMS
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 informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationFifth order multi-moment WENO schemes for hyperbolic conservation laws
J. Sci. Comput. manuscript No. (will be inserted by the editor) Fifth order multi-moment WENO schemes for hyperbolic conservation laws Chieh-Sen Huang Feng Xiao Todd Arbogast Received: May 24 / Accepted:
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
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 informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More 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 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 informationDiscontinuous Galerkin methods for first-order PDEs
Discontinuous Galerkin methods for first-order PDEs (ENPC) CNAM, Paris, 19.06.2015 Outline Friedrichs systems (steady linear PDEs) design of dg methods convergence analysis for smooth solutions unified
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
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 informationMAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS
MATHEMATICS OF COMPUTATION S 005-5780054-7 Article electronically published on May 6, 0 MAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS ORHAN MEHMETOGLU AND BOJAN POPOV Abstract.
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationFinite volume method for conservation laws V
Finite volume method for conservation laws V Schemes satisfying entropy condition Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
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 informationCalculus of Variations. Final Examination
Université Paris-Saclay M AMS and Optimization January 18th, 018 Calculus of Variations Final Examination Duration : 3h ; all kind of paper documents (notes, books...) are authorized. The total score of
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More information