Numerical Methods for Hyperbolic Conservation Laws Lecture 4
|
|
- Linda Sparks
- 5 years ago
- Views:
Transcription
1 Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
2 Nonlinear scalar conservation law Traffic flow: q: density of cars, u: velocity of cars. u = u(q) = u m (1 q) f (q) = qu(q) = u m q(1 q): flux q t + f (q) x = 0 Quasilinear form characteristics: X (t) = f (q(x (t), t)) q t + f (q)q x = 0 d dt q(x (t), t) = X (t)q x + q t = 0 characteristics (if exist) are straight lines with slope f (q), along which q remains constant. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
3 Rarefaction Waves Assume q x (x, 0) < 0 in traffic flow. Then f (q) = u m (1 q) f (q) = u m < 0 f is increasing in x. The characteristics are spreading out. In case of a Riemann problem, with q l > q r, we get a centered rarefaction wave. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
4 Formation of shocks If q x (x, 0) > 0 then characteristic can cross each other. multi-valued solutions Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
5 Vanishing viscosity and Equal-area rule Viscous equation: q t + f (q) x = ɛq xx Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
6 Shock speed Apply conservation law: xq r xq l = tf (q l ) tf (q r ) With s t = x: s(q r q l ) = f (q l ) f (q r ) so shock speed: s = f (q l) f (q r ) q l q r Rankine-Hugoniot jump condition Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
7 Solution of Riemann problem is self-similar, i.e, solution is constant along rays of x/t = c q(x, t) = q(x/t) Solve for rarefaction waves: q t (x, t) = x t q (x/t) = f (q) x = 1 t f ( q(x/t)) q (x/t) so f ( q(x/t)) = x/t Traffic flow: so q(x/t) = 1 [1 u m [1 q(x/t)] = x/t, x u m t ], for f (q l ) x/t f (q r ) Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
8 Weak solution Definition: The function q(x, t) is a weak solution of the conservation law q t + f (q) x = 0 with initial data q(x, 0) if the following holds for all test function φ C0 1 (i.e. C 1 functions with compact support): [qφ t + f (q)φ x ] dx dt = q(x, 0)φ(x, 0) dx 0 Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
9 Manipulating conservation laws Burgers equation shock speed: u t + (u /) x = 0 s = 1 ul ur = 1 u l u r (u l + u r ). Multiply the equation by u and manipulate: (u ) t + ( 3 u3 ) x = 0 shock speed s 1 = ul 3 ur 3 3 ul ur s s 1 Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation LawsOxford, LectureSpring, / 49
10 Entropy conditions Lax entropy condition A discontinuity (q l, q r ) with speed s (satisfying Rankine-Hogoniot condition) must satisfy f (q l ) > s > f (q r ) characteristics enter the shock from both sides. Oleinik entropy condition Assume f (q) > 0. The q(x, t) is the entropy solution if there is a constant E > 0 such that for all a > 0, t > 0 and x R, q(x + a, t) q(x, t) a < E t. decay of positive waves Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
11 Read: entropy functions, ch N-Wave decay, ch Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
12 Non-convex flux function Example: Two-phase flow and the Buckley-Levrett equation q: saturation of water 1 q: saturation of oil fractional flow: Conservation laws: f (q) = q q + a(1 q) (a < 1) q t + f (q) x = 0 Plot of f : S-shaped with an inflection point, where f changes sign. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
13 Solution of Riemann problem with q l = 1, q r = 0. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
14 Entropy condition for non-convex flux Lax condition is NOT enough! O-condition (by Oleinik) A discontinuity (q l, q r ) with speed s (satisfying R-H condition) is entropy admissible if f (q) f (q l ) s f (q) f (q r ) q q l q q r for all q between q l and q r. Graphical interpretation: If q l < q r, then the graph of f on [q l, q r ] must lie above the secant line connecting (q l, f (q l )) and (q r, f (q r )). If q l > q r, then the graph of f on [q r, q l ] must lie below the secant line connecting (q l, f (q l )) and (q r, f (q r )). Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
15 Convex Hull construction Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
16 Convex Hull construction example Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
17 FVM for nonlinear scalar conservation laws Conservative form Q n+1 i = Qi n t [ ] F n x i+ F n 1 i 1 Godunov s method of first order: F n i = f (Q ) 1 i 1 Q : solution of Riemann problem with Q i 1 i 1, Q i along ξ = 0. Assume either f < 0 or f > 0 (i.e. genuinely nonlinear/convex flux). Riemann solution consists of one simple wave, either a shock or a rarefaction. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
18 5 cases: Only in (c) is Q i 1 different from Q i 1, Q i. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
19 In (c) the solution is a rarefaction wave, with f (Q i 1 ) < 0 < f (Q i ). called transonic rarefaction (in gas dynamics). Recall solution of a rarefaction wave f ( q(ξ)) = ξ At ξ = 0, one has f (q s ) = 0. q s : stagnation point q s is between Q i 1 and Q i. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
20 Summary: Godunov flux for f > 0: f (Q i 1 ), if Q i 1 > q s and s > 0 F n i = f (Q 1 i ), if Q i < q s and s < 0 f (q s ), if Q i 1 < q s < Q i where s = f (Q i) f (Q i 1 ) Q i Q i 1 RH speed (similar for f < 0... leave as exercise.) Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
21 Also, one can use: F n i = 1 min f (q), if Q i 1 < Q i Q i 1 q Q i ( ) max f (q), if Q i 1 > Q i Q i q Q i 1 (*) holds for both convex and non-convex flux! Verify it (as an exercise). Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
22 If a shock with s = 0 occurs: By RH condition, one has f (Q i 1 ) = f (Q i ) one can use either f (Q i 1 ) or f (Q i ). Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
23 Define the fluctuations: A + Q i 1 A Q i 1 = f (Q i ) f (Q i 1 ) = f (Q i 1 ) f (Q i 1 ) Godunov method takes the same form as for linear equation: Q n+1 i = Q n i dt x (A+ Q i 1 + A Q i+ 1 ) Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
24 A simplified version: Solving every Riemann problem with a jump. W i 1 s i 1 = Q i Q i 1 { [f (Q i ) f (Q i 1 )]/(Q i Q i 1 ), if Q i 1 Q i = f (Q i ) if Q i 1 = Q i (Useful in high resolution correction terms.) Use A + Q i 1 A Q i 1 = s + i 1 W i 1 = s i 1 W i 1 (+) Simple (-) Trouble with a transonic rarefaction Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
25 Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
26 Entropy fix for transonic rarefaction: if f (Q i 1 ) < 0 < f (Q i ), use instead A + Q i 1 A Q i 1 = f (Q i ) f (q s ) = f (q s ) f (Q i 1 ) Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
27 Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
28 High resolution methods: Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
29 The Engquist-Osher method Approach: Treat all the solution of Riemann problem as a rarefaction. Set A + Q i 1 A Q i 1 = = Qi Q i 1 (f (q)) + dq Qi Q i 1 (f (q)) dq Then F i 1 Qi = f (Q i 1 ) + (f (q)) dq Q i 1 (or) = f (Q i ) Qi Q i 1 (f (q)) + dq = 1 [f (Q i 1) + f (Q i )] 1 Qi Q i 1 f (q) dq Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
30 F i 1 = 1 [f (Q i 1) + f (Q i )] 1 Qi Q i 1 f (q) dq If f > 0: F i 1 = f (Q i 1) If f < 0: F i 1 = f (Q i) If in transonic rarefaction: F i 1 = f (q s) If in transonic shock, i.e., f (Q i 1 ) > 0 > f (Q i ), then consistent F i 1 = f (Q i 1) + f (Q i ) f (q s ) (+) avoided expansion shock, and entropy condition always holds (+) extension to system is possible Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
31 E-scheme Osher showed that, if (E-scheme) then sign(q i Q i 1 ) [F i 1 f (q)] 0 q between Q i 1, Q i TVD if CFL number sufficiently small (leave as exercise) convergent to entropy weak solutions at most first order Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
32 High-resolution TVD methods One can extend from the method for linear problems. where Q n+1 i = Qi n t [ ] A + Q x i 1 + A Q i+ 1 t [ Fi+ 1 x F ] i 1 s i 1 : R-H speed F i 1 = 1 s i 1 W i 1 : limited version of the wave [ 1 t x s i 1 ] W i 1 Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
33 If one uses a TVD limiter φ(θ), the the method is TVD. convergence (NB! Only valid for scalar conservation laws!) Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
34 A partial informal analysis: Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
35 The important of the conservative form Burger equation u t + 1 (u ) x = 0 with u > 0. Using Upwind method: U n+1 i = Ui n 1 t [ (U n x i ) (Ui 1) n ] Discretize directly the quasilinear form u t + uu x = 0 U n+1 i = U n i t x Un i (U n i U n i 1) Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
36 Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
37 Lax-Wendroff Theorem Lax-Wendroff Theorem Let Q ε (x, t) denote the numerical approximation computed by a consistent and conservative method on grid ( t, x). Suppose that, as ε 0 (i.e. t, x 0), Q ε (x, t) converges to q(x, t), i.e., T b 0 a Q ε (x, t) q(x, t) dx dt 0 TV(Q ε (, t)) < M 0 t T, ε > 0 Then q(x, t) is a weak solution of the conservation law, i.e. 0 (qφ t + f (q)φ x ) dx dt = φ(x, 0)q(x, 0) dx ( ) for all test function φ C 1 0. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
38 Proof: Denote by (Q n i, Φn i ) the discrete values of (Qε, φ) at a grid point. Conservative method: Multiply with Φ n i and sum over i, n: n=0 i= Φ n i (Q n+1 i Using summation-by-parts: m k=1 Q n+1 i = Qi n t x [F n i+ F n 1 i ] 1 Q n i ) = t x n=0 i= Φ n i [F n i+ F n 1 i ] 1 m 1 a k (b k b k 1 ) = a m b m a 1 b 0 (a k+1 a k )b k Since φ has compact support, we have k=1 i Φ 0 i Q 0 i (Φ n i Φ n 1 i )Qi n = t x n=1 i (Φ n i+1 Φ n i )F n i 1 n=0 i Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
39 i Φ 0 i Q 0 i (Φ n i Φ n 1 i )Qi n = t x n=1 i (Φ n i+1 Φ n i )F n i 1 n=0 i Can write as I + II = III where I = x t II = x t (Φ n i Φ n 1 n=1 i n=0 i t i ) Qi n (Φ n i+1 Φn i ) F n x i 1 III = x i Φ 0 i Q 0 i Since φ C0 1, the limit as x 0, t 0 I 0 φ t q dx dt, III φ(x, 0)q(x, 0) dx Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
40 For II: use the consistence of the numerical flux. II = x t For example, F n = F(Q i 1 i 1 n, Qn i ), Use q = Q n i : n=0 i (Φ n i+1 Φn i ) F n x i 1 F(Q n i 1, Q n i ) f ( q) L max( Q n i q, Q n i 1 q ) F(Q n i 1, Q n i ) f ( q) L Q n i Since TV(Q n ) < M, as ε 0, we must have Q n i so F i 1 f (Qn i ) as ε 0 for almost all i. Then II q(x, t) is a weak solution. 0 Qi 1 n, 0 for almost all i. Q n i 1 φ x (x, t)f (q(x, t)) dx dt Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
41 Entropy condition But: weak solutions are not unique. Need entropy conditions. η(q): convex entropy function ψ(q): entropy flux, with ψ (q) = η (q)f (q) Entropy condition: η(q(x, t)) + ψ(q(x, t)) 0 t x in the weak sense Discrete version: enough to require η(q n+1 i ) η(qi n ) t ( ) Ψ n x i+ Ψ n 1 i 1 where Ψ n = Ψ(Q i 1 i 1 n, Qn i ): numerical entropy flux Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
42 Entropy consistency of Godunov method Given q n (x, t n ) piecewise constant. Let q n (x, t) be the exact entropy solution on [t n, t n+1 ]. Integrating η(q(x, t)) t + ψ(q(x, t)) x 0 over (x i 1, x i+ 1 ) (t n, t n+1 ): xi+ 1 x i 1 tn+1 xi+ η( q n 1 (x, t n+1 )) dx η( q n (x, t n )) dx x i 1 tn+1 ψ( q n (x i 1, t)) dt ψ( q n (x i+ 1 t n t n = t ψ(q ) t ψ(q ) i 1 i+ 1, t)) dt Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
43 So 1 xi+ 1 η( q n (x, t n+1 )) dx η(qi n ) t x x x (ψ(q ) ψ(q )) i+ 1 i 1 i 1 Since η (q) > 0, Jensen s inequality gives η(q n+1 i ) = η 1 xi+ 1 q n (x, t n+1 ) dx 1 x x x i 1 xi+ 1 x i 1 η( q n (x, t n+1 )) dx This leads to the discrete entropy condition: η(q n+1 i ) η(qi n ) t x (ψ(q ) ψ(q )) i+ 1 i 1 Solutions of Godunov method converges to entropy weak solutions. Wen Shen (Penn State) Numerical Methods for Hyperbolic Conservation Laws Oxford, Lecture Spring, / 49
AMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
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 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 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 informationThe Riemann problem. The Riemann problem Rarefaction waves and shock waves
The Riemann problem Rarefaction waves and shock waves 1. An illuminating example A Heaviside function as initial datum Solving the Riemann problem for the Hopf equation consists in describing the solutions
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
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 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 informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University Abstract We study two systems of conservation laws for polymer flooding in secondary
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University. Email: wxs27@psu.edu November 5, 2015 Abstract We study two systems of conservation
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
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 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 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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
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 informationModeling and Numerical Approximation of Traffic Flow Problems
Modeling and Numerical Approximation of Traffic Flow Problems Prof. Dr. Ansgar Jüngel Universität Mainz Lecture Notes (preliminary version) Winter Contents Introduction Mathematical theory for scalar conservation
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 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 information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
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 informationIntroduction to hyperbolic PDEs and Shallow Water Equations
Introduction to hyperbolic PDEs and Shallow Water Equations Technische Universität München Fundamentals of Wave Simulation - Solving Hyperbolic Systems of PDEs Supervisor Leonhard Rannabauer 8/12/2017
More informationGlobal Riemann Solver and Front Tracking Approximation of Three-Component Gas Floods
Global Riemann Solver and Front Tracking Approximation of Three-Component Gas Floods Saeid Khorsandi (1), Wen Shen (2) and Russell T. Johns (3) (1) Department of Energy and Mineral Engineering, Penn State
More informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationNONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL
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 informationDefinition and Construction of Entropy Satisfying Multiresolution Analysis (MRA)
Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 6 Definition and Construction of Entropy Satisfying Multiresolution Analysis (MRA Ju Y. Yi Utah State University
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 informationSolutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness
Solutions in the sense of distributions Definition, non uniqueness 1. Notion of distributions In order to build weak solutions to the Hopf equation, we need to define derivatives of non smooth functions,
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 informationNumerical Schemes Applied to the Burgers and Buckley-Leverett Equations
University of Reading Numerical Schemes Applied to the Burgers and Buckley-Leverett Equations by September 4 Department of Mathematics Submitted to the Department of Mathematics, University of Reading,
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 informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
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 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 informationMATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS
K CHAPMAN & HALL/CRC Monographs and Surveys in Pure and Applied Mathematics I 18 MATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS ANDREI G. KULIKOVSKII NIKOLAI V. POGORELOV ANDREI YU. SEMENOV
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 informationTraffic Flow Problems
Traffic Flow Problems Nicodemus Banagaaya Supervisor : Dr. J.H.M. ten Thije Boonkkamp October 15, 2009 Outline Introduction Mathematical model derivation Godunov Scheme for the Greenberg Traffic model.
More informationGalilean invariance and the conservative difference schemes for scalar laws
RESEARCH Open Access Galilean invariance and the conservative difference schemes for scalar laws Zheng Ran Correspondence: zran@staff.shu. edu.cn Shanghai Institute of Applied Mathematics and Mechanics,
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationExistence Theory for Hyperbolic Systems of Conservation Laws with General Flux-Functions
Existence Theory for Hyperbolic Systems of Conservation Laws with General Flux-Functions Tatsuo Iguchi & Philippe G. LeFloch Abstract For the Cauchy problem associated with a nonlinear, strictly hyperbolic
More informationQuasi-linear first order equations. Consider the nonlinear transport equation
Quasi-linear first order equations Consider the nonlinear transport equation u t + c(u)u x = 0, u(x, 0) = f (x) < x < Quasi-linear first order equations Consider the nonlinear transport equation u t +
More informationNotes: Outline. Shallow water equations. Notes: Shallow water equations. Notes:
Outline Nonlinear hyperbolic systems Shallow water equations Shock waves and Hugoniot loci Integral curves in phase plane Compression and rarefaction R.J. LeVeque, University of Washington IPDE 2011, July
More informationSHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS
SHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE Abstract. The present paper deals with the following hyperbolic elliptic coupled system, modelling
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
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 informationModeling Rarefaction and Shock waves
April 30, 2013 Inroduction Inroduction Rarefaction and shock waves are combinations of two wave fronts created from the initial disturbance of the medium. Inroduction Rarefaction and shock waves are combinations
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 informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More informationHyperbolic Conservation Laws Past and Future
Hyperbolic Conservation Laws Past and Future Barbara Lee Keyfitz Fields Institute and University of Houston bkeyfitz@fields.utoronto.ca Research supported by the US Department of Energy, National Science
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 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 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 informationShock on the left: locus where cars break behind the light.
Review/recap of theory so far. Evolution of wave profile, as given by the characteristic solution. Graphical interpretation: Move each point on graph at velocity c(ρ). Evolution as sliding of horizontal
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 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 informationOn finite time BV blow-up for the p-system
On finite time BV blow-up for the p-system Alberto Bressan ( ), Geng Chen ( ), and Qingtian Zhang ( ) (*) Department of Mathematics, Penn State University, (**) Department of Mathematics, University of
More informationEducjatipnal. L a d ie s * COBNWALILI.S H IG H SCHOOL. I F O R G IR L S A n B k i n d e r g a r t e n.
- - - 0 x ] - ) ) -? - Q - - z 0 x 8 - #? ) 80 0 0 Q ) - 8-8 - ) x ) - ) -] ) Q x?- x - - / - - x - - - x / /- Q ] 8 Q x / / - 0-0 0 x 8 ] ) / - - /- - / /? x ) x x Q ) 8 x q q q )- 8-0 0? - Q - - x?-
More informationRegularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws
1 Regularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws Sébastien Blandin, Xavier Litrico, Maria Laura Delle Monache, Benedetto Piccoli and Alexandre Bayen Abstract
More informationCONVERGENCE OF THE LAX-FRIEDRICHS SCHEME AND STABILITY FOR CONSERVATION LAWS WITH A DISCONTINUOUS SPACE-TIME DEPENDENT FLUX
Chin. Ann. Math. 25B:3(24),287 318. CONVERGENCE OF THE LAX-FRIEDRICHS SCHEME AND STABILITY FOR CONSERVATION LAWS WITH A DISCONTINUOUS SPACE-TIME DEPENDENT FLUX K. H. KARLSEN J. D. TOWERS Abstract The authors
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 information1. Introduction. We study the convergence of vanishing viscosity solutions governed by the single conservation law
SIAM J. NUMER. ANAL. Vol. 36, No. 6, pp. 1739 1758 c 1999 Society for Industrial and Applied Mathematics POINTWISE ERROR ESTIMATES FOR SCALAR CONSERVATION LAWS WITH PIECEWISE SMOOTH SOLUTIONS EITAN TADMOR
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationGas Dynamics Equations: Computation
Title: Name: Affil./Addr.: Gas Dynamics Equations: Computation Gui-Qiang G. Chen Mathematical Institute, University of Oxford 24 29 St Giles, Oxford, OX1 3LB, United Kingdom Homepage: http://people.maths.ox.ac.uk/chengq/
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationIntroduction to nonlinear wave models
Introduction to nonlinear wave models Marko Nedeljkov Department of Mathematics and Informatics, University of Novi Sad Trg Dositeja Obradovića 4, 21000 Novi Sad, Serbia marko.nedeljkov@dmi.uns.ac.rs http://www.dmi.uns.ac.rs
More informationApproximate Solutions of Generalized Riemann Problems for Hyperbolic Conservation Laws and Their Application to High Order Finite Volume Schemes
Approximate Solutions of Generalized Riemann Problems for Hyperbolic Conservation Laws and Their Application to High Order Finite Volume Schemes Dissertation zur Erlangung des Doktorgrades der Fakultät
More informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
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 informationA re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws
A re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c,3 a Department of Applied Mathematics and National
More informationCLASSROOM NOTES PART II: SPECIAL TOPICS. APM526, Spring 2018 Last update: Apr 11
CLASSROOM NOTES PART II: SPECIAL TOPICS APM526, Spring 2018 Last update: Apr 11 1 Function Space Methods General Setting: Projection into finite dimensional subspaces t u = F (u), u(t = 0) = u I, F : B
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationCoupling conditions for transport problems on networks governed by conservation laws
Coupling conditions for transport problems on networks governed by conservation laws Michael Herty IPAM, LA, April 2009 (RWTH 2009) Transport Eq s on Networks 1 / 41 Outline of the Talk Scope: Boundary
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 informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 253 Monotonicity of Implicit Finite Difference Methods for Hyperbolic Conservation Laws Michael
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 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 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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationRegularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws
Regularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws Sébastien Blandin, Xavier Litrico, Maria Laura Delle Monache, Benedetto Piccoli and Alexandre Bayen Abstract
More informationA FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS
A FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS Ruben Juanes 1 and Knut-Andreas Lie 2 1 Stanford University, Dept. Petroleum Engineering, USA 2 SINTEF IKT, Dept., Norway ECMOR IX, August 30 September
More informationFinite Volume Methods for Hyperbolic Problems
Finite Volume Methods for Hyperbolic Problems This book contains an introduction to hyperbolic partial differential equations and a powerful class of numerical methods for approximating their solution,
More informationClassification of homogeneous quadratic conservation laws with viscous terms
Volume 26, N. 2, pp. 251 283, 2007 Copyright 2007 SBMAC ISSN 0101-8205 www.scielo.br/cam Classification of homogeneous quadratic conservation laws with viscous terms JANE HURLEY WENSTROM 1 and BRADLEY
More informationWAVE INTERACTIONS AND NUMERICAL APPROXIMATION FOR TWO DIMENSIONAL SCALAR CONSERVATION LAWS
Computational Fluid Dynamics JOURNAL 4(4):46 January 6 (pp.4 48) WAVE INTERACTIONS AND NUMERICAL APPROXIMATION FOR TWO DIMENSIONAL SCALAR CONSERVATION LAWS Matania BEN ARTZI Joseph FALCOVITZ Jiequan LI
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 informationAN ALTERNATING DESCENT METHOD FOR THE OPTIMAL CONTROL OF THE INVISCID BURGERS EQUATION IN THE PRESENCE OF SHOCKS
December 7, 2007 10:10 WSPC/INSTRUCTION FILE Castro- Palacios- Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company AN ALTERNATING DESCENT METHOD FOR THE OPTIMAL CONTROL
More informationAPPROXIMATE SOLUTIONS OF GENERALIZED RIEMANN PROBLEMS FOR NONLINEAR SYSTEMS OF HYPERBOLIC CONSERVATION LAWS
APPROXIMATE SOLUTIONS OF GENERALIZED RIEMANN PROBLEMS FOR NONLINEAR SYSTEMS OF HYPERBOLIC CONSERVATION LAWS CLAUS R. GOETZ AND ARMIN ISKE Abstract. We study analytical properties of the Toro-Titarev solver
More informationSolving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation
Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation W.L. Jin and H.M. Zhang August 3 Abstract: In this paper we study the Payne-Whitham (PW) model as
More informationGeneralized Riemann Problems: From the Scalar Equation to Multidimensional Fluid Dynamics
1 Generalized Riemann Problems: From the Scalar Equation to Multidimensional Fluid Dynamics Matania Ben-Artizi and Joseph Falcovitz Institute of Mathematics, the Hebrew University of Jerusalem, Jerusalem
More informationConvergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates
Convergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates Xin Wen and Shi Jin Abstract We study the L 1 -error estimates
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 information