AMath 574 February 11, 2011
|
|
- Evelyn Wilkerson
- 5 years ago
- Views:
Transcription
1 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 14, 2011
2 Entropy-violating numerical solutions Riemann problem for Burgers equation at t = 1 with u l = 1 and u r = 2: 2.5 Godunov with no entropy fix 2.5 Godunov with entropy fix High resolution with no entropy fix 2.5 High resolution with entropy fix R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 12.3]
3 Vanishing viscosity solution We want q(x, t) to be the limit as ɛ 0 of solution to q t + f(q) x = ɛq xx. This selects a unique weak solution: Shock if f (q l ) > f (q r ), Rarefaction if f (q l ) < f (q r ). R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
4 Vanishing viscosity solution We want q(x, t) to be the limit as ɛ 0 of solution to q t + f(q) x = ɛq xx. This selects a unique weak solution: Shock if f (q l ) > f (q r ), Rarefaction if f (q l ) < f (q r ). Lax Entropy Condition: A discontinuity propagating with speed s in the solution of a convex scalar conservation law is admissible only if f (q l ) > s > f (q r ), where s = (f(q r ) f(q l ))/(q r q l ). R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
5 Vanishing viscosity solution We want q(x, t) to be the limit as ɛ 0 of solution to q t + f(q) x = ɛq xx. This selects a unique weak solution: Shock if f (q l ) > f (q r ), Rarefaction if f (q l ) < f (q r ). Lax Entropy Condition: A discontinuity propagating with speed s in the solution of a convex scalar conservation law is admissible only if f (q l ) > s > f (q r ), where s = (f(q r ) f(q l ))/(q r q l ). Note: This means characteristics must approach shock from both sides as t advances, not move away from shock! R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
6 Approximate Riemann solvers For nonlinear problems, computing the exact solution to each Riemann problem may not be possible, or too expensive. Often the nonlinear problem q t + f(q) x = 0 is approximated by q t + A i 1/2 q x = 0, q l = Q i 1, q r = Q i for some choice of A i 1/2 f (q) based on data Q i 1, Q i. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
7 Approximate Riemann solvers For nonlinear problems, computing the exact solution to each Riemann problem may not be possible, or too expensive. Often the nonlinear problem q t + f(q) x = 0 is approximated by q t + A i 1/2 q x = 0, q l = Q i 1, q r = Q i for some choice of A i 1/2 f (q) based on data Q i 1, Q i. Solve linear system for α i 1/2 : Q i Q i 1 = p αp i 1/2 rp i 1/2. Waves W p i 1/2 = αp i 1/2 rp i 1/2 propagate with speeds sp i 1/2, r p i 1/2 are eigenvectors of A i 1/2, s p i 1/2 are eigenvalues of A i 1/2. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
8 Approximate Riemann solvers q t + Âi 1/2q x = 0, q l = Q i 1, q r = Q i Often Âi 1/2 = f (Q i 1/2 ) for some choice of Q i 1/2. In general  i 1/2 = Â(q l, q r ). R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
9 Approximate Riemann solvers q t + Âi 1/2q x = 0, q l = Q i 1, q r = Q i Often Âi 1/2 = f (Q i 1/2 ) for some choice of Q i 1/2. In general  i 1/2 = Â(q l, q r ). Roe conditions for consistency and conservation: Â(q l, q r ) f (q ) as q l, q r q,  diagonalizable with real eigenvalues, For conservation in wave-propagation form,  i 1/2 (Q i Q i 1 ) = f(q i ) f(q i 1 ). R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
10 Approximate Riemann solvers For a scalar problem, we can easily satisfy the Roe condition  i 1/2 (Q i Q i 1 ) = f(q i ) f(q i 1 ). by choosing  i 1/2 = f(q i) f(q i 1 ) Q i Q i 1. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
11 Approximate Riemann solvers For a scalar problem, we can easily satisfy the Roe condition  i 1/2 (Q i Q i 1 ) = f(q i ) f(q i 1 ). by choosing  i 1/2 = f(q i) f(q i 1 ) Q i Q i 1. Then r 1 i 1/2 = 1 and s1 i 1/2 = Âi 1/2 (scalar!). Note: This is the Rankine-Hugoniot shock speed. = shock waves are correct, rarefactions replaced by entropy-violating shocks. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
12 Weak solutions to Burgers equation u t + ( 1 2 u2) x = 0, u l = 1, u r = 2 Characteristic speed: u Rankine-Hugoniot speed: 1 2 (u l + u r ). Physically correct rarefaction wave solution: R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
13 Weak solutions to Burgers equation u t + ( 1 2 u2) x = 0, u l = 1, u r = 2 Characteristic speed: u Rankine-Hugoniot speed: 1 2 (u l + u r ). Entropy violating weak solution: R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
14 Weak solutions to Burgers equation u t + ( 1 2 u2) x = 0, u l = 1, u r = 2 Characteristic speed: u Rankine-Hugoniot speed: 1 2 (u l + u r ). Another Entropy violating weak solution: R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
15 Transonic rarefactions Sonic point: u s = 0 for Burgers since f (0) = 0. Consider Riemann problem data u l = 0.5 < 0 < u r = 1.5. In this case wave should spread in both directions: R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
16 Transonic rarefactions Entropy-violating approximate Riemann solution: s = 1 2 (u l + u r ) = 0.5. Wave goes only to right, no update to cell average on left. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
17 Transonic rarefactions If u l = u r then Rankine-Hugoniot speed is 0: Similar solution will be observed with Godunov s method if entropy-violating approximate Riemann solver used. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
18 Entropy-violating numerical solutions Riemann problem for Burgers equation at t = 1 with u l = 1 and u r = 2: 2.5 Godunov with no entropy fix 2.5 Godunov with entropy fix High resolution with no entropy fix 2.5 High resolution with entropy fix R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 12.3]
19 Approximate Riemann solver Q n+1 i = Q n i t [ A + Q x i 1/2 + A ] Q i+1/2. For scalar advection m = 1, only one wave. W i 1/2 = Q i 1/2 = Q i Q i 1 and s i 1/2 = u, A Q i 1/2 = s i 1/2 W i 1/2, A + Q i 1/2 = s + i 1/2 W i 1/2. For scalar nonlinear: Use same formulas with W i 1/2 = Q i 1/2 and s i 1/2 = F i 1/2 / Q i 1/2. Need to modify these by an entropy fix in the trans-sonic rarefaction case. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 12.3]
20 Entropy fix Q n+1 i Revert to the formulas = Q n i t [ A + Q x i 1/2 + A ] Q i+1/2. A Q i 1/2 = f(q s ) f(q i 1 ) A + Q i 1/2 = f(q i ) f(q s ) if f (Q i 1 ) < 0 < f (Q i ). left-going fluctuation right-going fluctuation High-resolution method: still define wave W and speed s by W i 1/2 = Q i Q i 1, { (f(qi ) f(q s i 1/2 = i 1 ))/(Q i Q i 1 ) if Q i 1 Q i f (Q i ) if Q i 1 = Q i. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 12.3]
21 Godunov flux for scalar problem The Godunov flux function for the case f (q) > 0 is F n i 1/2 = = f(q i 1 ) if Q i 1 > q s and s > 0 f(q i ) if Q i < q s and s < 0 f(q s ) if Q i 1 < q s < Q i. min f(q) Q i 1 q Q i if Q i 1 Q i max f(q) if Q i Q i 1, Q i q Q i 1 Here s = f(q i) f(q i 1 ) Q i Q i 1 is the Rankine-Hugoniot shock speed. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 12.1]
22 Entropy-violating numerical solutions Riemann problem for Burgers equation with q l = 1 and q r = 2: 2.5 Godunov with no entropy fix 2.5 Godunov with entropy fix High resolution with no entropy fix 2.5 High resolution with entropy fix R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 12.3]
23 Entropy (admissibility) conditions We generally require additional conditions on a weak solution to a conservation law, to pick out the unique solution that is physically relevant. In gas dynamics: entropy is constant along particle paths for smooth solutions, entropy can only increase as a particle goes through a shock. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
24 Entropy (admissibility) conditions We generally require additional conditions on a weak solution to a conservation law, to pick out the unique solution that is physically relevant. In gas dynamics: entropy is constant along particle paths for smooth solutions, entropy can only increase as a particle goes through a shock. Entropy functions: Function of q that behaves like physical entropy for the conservation law being studied. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
25 Entropy (admissibility) conditions We generally require additional conditions on a weak solution to a conservation law, to pick out the unique solution that is physically relevant. In gas dynamics: entropy is constant along particle paths for smooth solutions, entropy can only increase as a particle goes through a shock. Entropy functions: Function of q that behaves like physical entropy for the conservation law being studied. NOTE: Mathematical entropy functions generally chosen to decrease for admissible solutions, increase for entropy-violating solutions. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
26 Entropy functions A scalar-valued function η : lr m lr if the Hessian matrix η (q) with (i, j) element η ij(q) = 2 η q i q j is positive definite for all q, i.e., satisfies is a convex function of q v T η (q)v > 0 for all q, v lr m. Scalar case: reduces to η (q) > 0. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
27 Entropy functions Entropy function: η : lr m lr Entropy flux: ψ : lr m lr chosen so that η(q) is convex and: η(q) is conserved wherever the solution is smooth, η(q) t + ψ(q) x = 0. Entropy decreases across an admissible shock wave. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
28 Entropy functions Entropy function: η : lr m lr Entropy flux: ψ : lr m lr chosen so that η(q) is convex and: η(q) is conserved wherever the solution is smooth, η(q) t + ψ(q) x = 0. Entropy decreases across an admissible shock wave. Weak form: x2 x 1 η(q(x, t 2 )) dx + x2 x 1 t2 t 1 η(q(x, t 1 )) dx with equality where solution is smooth. t2 ψ(q(x 1, t)) dt t 1 ψ(q(x 2, t)) dt R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
29 Entropy functions How to find η and ψ satisfying this? η(q) t + ψ(q) x = 0 For smooth solutions gives η (q)q t + ψ (q)q x = 0. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
30 Entropy functions How to find η and ψ satisfying this? η(q) t + ψ(q) x = 0 For smooth solutions gives η (q)q t + ψ (q)q x = 0. Since q t = f (q)q x this is satisfied provided ψ (q) = η (q)f (q) R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
31 Entropy functions How to find η and ψ satisfying this? η(q) t + ψ(q) x = 0 For smooth solutions gives η (q)q t + ψ (q)q x = 0. Since q t = f (q)q x this is satisfied provided ψ (q) = η (q)f (q) Scalar: Can choose any convex η(q) and integrate. Example: Burgers equation, f (u) = u and take η(u) = u 2. Then ψ (u) = 2u 2 = Entropy function: ψ(u) = 2 3 u3. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
32 Weak solutions and entropy functions The conservation laws ( ) 1 u t + 2 u2 = 0 and both have the same quasilinear form x u t + uu x = 0 ( u 2 ) t + ( 2 3 u3 ) x = 0 but have different weak solutions, different shock speeds! R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
33 Weak solutions and entropy functions The conservation laws ( ) 1 u t + 2 u2 = 0 and both have the same quasilinear form x u t + uu x = 0 ( u 2 ) t + ( 2 3 u3 ) x = 0 but have different weak solutions, different shock speeds! Entropy function: η(u) = u 2. A correct Burgers shock at speed s = 1 2 (u l + u r ) will have total mass of η(u) decreasing. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
34 Entropy functions x2 x 1 η(q(x, t 2 )) dx + x2 x 1 t2 t 1 η(q(x, t 1 )) dx t2 ψ(q(x 1, t)) dt ψ(q(x 2, t)) dt t 1 comes from considering the vanishing viscosity solution: q ɛ t + f(q ɛ ) x = ɛq ɛ xx R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
35 Entropy functions x2 x 1 η(q(x, t 2 )) dx + x2 x 1 t2 t 1 η(q(x, t 1 )) dx t2 ψ(q(x 1, t)) dt ψ(q(x 2, t)) dt t 1 comes from considering the vanishing viscosity solution: Multiply by η (q ɛ ) to obtain: q ɛ t + f(q ɛ ) x = ɛq ɛ xx η(q ɛ ) t + ψ(q ɛ ) x = ɛη (q ɛ )q ɛ xx. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
36 Entropy functions x2 x 1 η(q(x, t 2 )) dx + x2 x 1 t2 t 1 η(q(x, t 1 )) dx t2 ψ(q(x 1, t)) dt ψ(q(x 2, t)) dt t 1 comes from considering the vanishing viscosity solution: qt ɛ + f(q ɛ ) x = ɛqxx ɛ Multiply by η (q ɛ ) to obtain: η(q ɛ ) t + ψ(q ɛ ) x = ɛη (q ɛ )qxx. ɛ Manipulate further to get η(q ɛ ) t + ψ(q ɛ ) x = ɛ ( η (q ɛ )qx ɛ ) x ɛη (q ɛ ) (qx) ɛ 2. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
37 Entropy functions Smooth solution to viscous equation satisfies η(q ɛ ) t + ψ(q ɛ ) x = ɛ ( η (q ɛ )qx ɛ ) x ɛη (q ɛ ) (qx) ɛ 2. Integrating over rectangle [x 1, x 2 ] [t 1, t 2 ] gives x2 x 1 η(q ɛ (x, t 2 )) dx = + ɛ ɛ ( t2 t 1 t2 t 1 t2 x2 t 1 x2 x 1 ψ(q ɛ (x 2, t)) dt η(q ɛ (x, t 1 )) dx t2 t 1 ) ψ(q ɛ (x 1, t)) dt [ η (q ɛ (x 2, t)) q ɛ x(x 2, t) η (q ɛ (x 1, t)) q ɛ x(x 1, t) ] dt x 1 η (q ɛ ) (q ɛ x) 2 dx dt. Let ɛ 0 to get result: Term on third line goes to 0, Term of fourth line is always 0. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
38 Entropy functions Weak form of entropy condition: 0 [ φt η(q) + φ x ψ(q) ] dx dt + φ(x, 0)η(q(x, 0)) dx 0 for all φ C0 1 (lr lr) with φ(x, t) 0 for all x, t. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
39 Entropy functions Weak form of entropy condition: 0 [ φt η(q) + φ x ψ(q) ] dx dt + φ(x, 0)η(q(x, 0)) dx 0 for all φ C0 1 (lr lr) with φ(x, t) 0 for all x, t. Informally we may write η(q) t + ψ(q) x 0. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec. 11.4]
40 Lax-Wendroff Theorem Suppose the method is conservative and consistent with q t + f(q) x = 0, F i 1/2 = F(Q i 1, Q i ) with F( q, q) = f( q) and Lipschitz continuity of F. If a sequence of discrete approximations converge to a function q(x, t) as the grid is refined, then this function is a weak solution of the conservation law. Note: Does not guarantee a sequence converges (need stability). Two sequences might converge to different weak solutions. Also need to satisfy an entropy condition. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
41 Sketch of proof of Lax-Wendroff Theorem Multiply the conservative numerical method by Φ n i to obtain Q n+1 i Φ n i Q n+1 i = Q n i t x (F n i+1/2 F n i 1/2 ) = Φ n i Q n i t x Φn i (F n i+1/2 F n i 1/2 ). R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
42 Sketch of proof of Lax-Wendroff Theorem Multiply the conservative numerical method by Φ n i to obtain Q n+1 i Φ n i Q n+1 i = Q n i t x (F n i+1/2 F n i 1/2 ) = Φ n i Q n i t x Φn i (F n i+1/2 F n i 1/2 ). This is true for all values of i and n on each grid. Now sum over all i and n 0 to obtain n=0 i= Φ n i (Q n+1 i Q n i ) = t x n=0 i= Φ n i (F n i+1/2 F n i 1/2 ). Use summation by parts to transfer differences to Φ terms. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
43 Sketch of proof of Lax-Wendroff Theorem Obtain analog of weak form of conservation law: x t [ n=1 i= + ( Φ n i Φ n 1 i t n=0 i= ) Q n i ( Φ n i+1 Φ n i x ) ] Fi 1/2 n = x i= Φ 0 i Q 0 i. Consider on a sequence of grids with x, t 0. Show that any limiting function must satisfy weak form of conservation law. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
44 Analog of Lax-Wendroff proof for entropy Show that the numerical flux function F leads to a numerical entropy flux Ψ such that the following discrete entropy inequality holds: η(q n+1 i ) η(q n i ) t [ ] Ψ n i+1/2 x Ψn i 1/2. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
45 Analog of Lax-Wendroff proof for entropy Show that the numerical flux function F leads to a numerical entropy flux Ψ such that the following discrete entropy inequality holds: η(q n+1 i ) η(q n i ) t [ ] Ψ n i+1/2 x Ψn i 1/2. Then multiply by test function Φ n i, sum and use summation by parts to get discrete form of integral form of entropy condition. = If numerical approximations converge to some function, then the limiting function satisfies the entropy condition. R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
46 Entropy consistency of Godunov s method For Godunov s method, F (Q i 1, Q i ) = f(q i 1/2 ) where Q i 1/2 is the constant value along x i 1/2 in the Riemann solution. Let Ψ n i 1/2 = ψ(q i 1/2 ) Discrete entropy inequality follows from Jensen s inequality: The value of η evaluated at the average value of q n is less than or equal to the average value of η( q n ), i.e., ( ) 1 xi+1/2 η q n (x, t n+1 ) dx 1 xi+1/2 η ( q n (x, t n+1 ) ) dx. x x i 1/2 x x i 1/2 R.J. LeVeque, University of Washington AMath 574, February 14, 2011 [FVMHP Sec ]
Notes: 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 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 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 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 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 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 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 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 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 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 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 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 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 informationCharacteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:
AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,
More informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
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 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 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 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 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 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 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 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 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 informationClawpack Tutorial Part I
Clawpack Tutorial Part I Randall J. LeVeque Applied Mathematics University of Washington Conservation Laws Package www.clawpack.org (pdf s will be posted and green links can be clicked) Some collaborators
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 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 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 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 informationNumerical and mathematical analysis of a five-equation model for two-phase flow
SINTEF Energy Research Master s Thesis PFE-Master Numerical and mathematical analysis of a five-equation model for two-phase flow Author Pedro José Martínez Ferrer Supervisors Svend Tollak Munkejord Tore
More informationWave Propagation Software, Computational Science, and Reproducible Research
Wave Propagation Software, Computational Science, and Reproducible Research Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF and DOE Outline Hyperbolic
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 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 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 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 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 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 informationNumerical simulation of cylindrical solitary waves in periodic media
Numerical simulation of cylindrical solitary waves in periodic media Manuel Quezada de Luna David I. Ketcheson June, 3 Abstract We study the behavior of nonlinear waves in a two-dimensional medium with
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More 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 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 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 informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
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 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 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 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 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 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 informationSOME PRELIMINARY NOTES ON HYPERBOLIC CONSERVATION LAWS
SOME PRELIMINARY NOTES ON HYPERBOLIC CONSERVATION LAWS KAYYUNNAPARA T. JOSEPH CENTRE FOR APPLICABLE MATHEMATICS TATA INSTITUTE OF FUNDAMENTAL RESEARCH BANGALORE 5665 1. Introduction Many of the physical
More informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationKINETIC INDUCED MOMENT SYSTEMS FOR THE SAINT VENANT EQUATIONS
KINETIC INDUCED MOMENT SYSTEMS FOR THE SAINT VENANT EQUATIONS DIANA CRISTINA GIL MONTOYA Faculty of Applied Physics and Mathematics Gdans University of Technology, Poland Department of Mathematics University
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 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 informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationc Copyright 2014 Jihwan Kim
c Copyright 2014 Jihwan Kim Finite Volume Methods for Tsunamis Generated by Submarine Landslides Jihwan Kim A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor
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 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 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 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 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 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 informationTwo-layer shallow water system and its applications
Proceedings of Symposia in Applied Mathematics Two-layer shallow water system and its applications Jihwan Kim and Randall J. LeVeque Abstract. The multi-layer shallow water system is derived by depth averaging
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 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 informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationLarge Time Step Scheme Behaviour with Different Entropy Fix
Proceedings of the Paistan Academy of Sciences: A. Physical and Computational Sciences 53 (): 13 4 (016) Copyright Paistan Academy of Sciences ISSN: 0377-969 (print), 306-1448 (online) Paistan Academy
More informationHIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF STATES FOR SOLVING HYPERBOLIC SYSTEMS WITH NONCONSERVATIVE PRODUCTS
MATHEMATICS OF COMPUTATION Volume 75, Number 255, July 26, Pages 113 1134 S 25-5718(6)1851-5 Article electronically published on March 21, 26 HIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF
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 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 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 CONSERVATION PROPERTIES OF H(DIV )-CONFORMING LEAST-SQUARES FINITE ELEMENT METHODS FOR THE BURGERS EQUATION
SIAM JOURNAL ON SCIENTIFIC COMPUTING, SUBMITTED JUNE 23, ACCEPTED JUNE 24 NUMERICAL CONSERVATION PROPERTIES OF H(DIV )-CONFORMING LEAST-SQUARES FINITE ELEMENT METHODS FOR THE BURGERS EQUATION H. DE STERCK,
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 informationA Well-Balanced Path-Integral f-wave Method for Hyperbolic Problems with Source Terms
A Well-Balanced Path-Integral f-wave Method for Hyperbolic Problems with Source Terms Randall J. LeVeque 1 Abstract. Systems of hyperbolic partial differential equations with source terms (balance laws)
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 informationEntropy and the numerical integration of conservation laws
Pysics Procedia Pysics Procedia 00 2011) 1 28 Entropy and te numerical integration of conservation laws Gabriella Puppo Dipartimento di Matematica, Politecnico di Torino Italy) Matteo Semplice Dipartimento
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 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 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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationHeuristical and numerical considerations for the carbuncle phenomenon
Heuristical and numerical considerations for the carbuncle phenomenon arxiv:57.666v2 [physics.comp-ph] 6 Mar 27 Friedemann Kemm March 7, 27 Keywords: Carbuncle phenomenon, High speed flow, Shock instability,
More informationA Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws
A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws S. Jaisankar, S. V. Raghurama Rao CFD Centre, Department of Aerospace Engineering, Indian Institute of Science, Bangalore, India-56 Abstract
More informationEntropy-based moment closure for kinetic equations: Riemann problem and invariant regions
Entropy-based moment closure for kinetic equations: Riemann problem and invariant regions Jean-François Coulombel and Thierry Goudon CNRS & Université Lille, Laboratoire Paul Painlevé, UMR CNRS 854 Cité
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 informationNumerical Solution of the Six-Equation Two-Phase Single Velocity Model with Finite Relaxation
Numerical Solution of the Six-Equation Two-Phase Single Velocity Model with Finite Relaxation Sigbjørn Løland Bore Applied Physics, NTNU December 19, 214 Abstract In this report we investigate a recently
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 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 informationCONVERGENCE RATE OF APPROXIMATE SOLUTIONS TO CONSERVATION LAWS WITH INITIAL RAREFACTIONS
CONVERGENCE RATE OF APPROXIMATE SOLUTIONS TO CONSERVATION LAWS WITH INITIAL RAREFACTIONS HAIM NESSYAHU AND TAMIR TASSA Abstract. We address the question of local convergence rate of conservative Lip +
More informationThe RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
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 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 informationShock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids
Shock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids Andrea Mentrelli Department of Mathematics & Research Center of Applied Mathematics (CIRAM) University of Bologna, Italy Summary
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 informationScalar conservation laws with moving density constraints arising in traffic flow modeling
Scalar conservation laws with moving density constraints arising in traffic flow modeling Maria Laura Delle Monache Email: maria-laura.delle monache@inria.fr. Joint work with Paola Goatin 14th International
More informationSTRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION
Electronic Journal of Differential Equations, Vol. 216 (216, No. 126, pp. 1 16. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN
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 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 informationClawpack Tutorial Part I
Clawpack Tutorial Part I Randall J. LeVeque Applied Mathematics University of Washington Conservation Laws Package www.clawpack.org Slides posted at http://www.clawpack.org/links/tutorials http://faculty.washington.edu/rjl/tutorials
More information