Introduction to hyperbolic PDEs and Shallow Water Equations
|
|
- Cameron Wilson
- 5 years ago
- Views:
Transcription
1 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
2 Presentation Topics Hyperbolic PDEs Conservation Laws Solution of Advection Equation The Riemann Problem linear and nonlinear characteristics The weak solution Rankine Hugoniot Condition Admissibility Conditions Shock and rarefaction waves Analytical solution of the Dam-Break Problem 2/57
3 Hyperbolic PDEs Definition: An nxn system of equations is strictly hyperbolic if the Jacobian matrix has n real distinct eigenvalues. 3/57
4 Conservation Laws Definition: The overall particular measurable property of an isolated physical system does not change as the system evolves over time. 4/57
5 Conservation Laws formulation [time change of quantity]+[inflow]-[outflow]+[generation]-[destruction]=0 The generation and destruction terms can be omitted for the homogeneous case The mathematical formulation is simplified to ut+f(u)x=0 If f(u)x is smooth ut+fx(u)ux=0 5/57
6 The advection equation ut+ūux=0 u(x,t0)=u0(x,0) (or in conservation form ut+fx(u)ux=0 with fx(u)=ū) ubar=constant Solve with method of characteristics: introduce a new variable s. This is the gradient of the line in the x-t plane along which u remains constant. du ds =0 And by using the chain rule we get dt =1 ds and dx =ū dt Consequently: s=t (t0=0) and x=x0+ūt With these data, solving du =0 ds we get u(x,t)= u0 (x0)=u0(x-ūt) 6/57
7 Characteristics of the advection equation Interpretation of the advection equation solution Image from 1 The lines are parallel as the system is linear 7/57
8 The Riemann problem A conservation system with piecewise initial conditions e.g. ut+ūux=0 With initial conditions q ( x, 0)= { q l for x <0 q r for x > 0 It follows that the solution of the 1d advection equation of the Riemann problem is Image from 1 8/57
9 Generalization for nxn linear conservation system Briefly, since the system is hyperbolic it can be decomposed as n independent equations and the summation of their solutions is the general solution qt+aqx=0 where q is a vector of conserved quantities A is the flux jacobian matrix. In the linear case it is constant Since the system is hyperbolic it has a complete set of eigenvectors and real eigenvalues A R-1ΛR where R ={r1,r2,,rn} and Λ=diag{λ1,λ2,.λn} 9/57
10 System decoupling qt+rλr-1qx=0 R-1qt+ΛR-1qx=0 For ease, a new variable ω is assigned ω=r-1q ωtt+λωx=0 and ω(x,0)= R-1q(x,0) these leads to n scalar advection equations ωtp+λpωxp=0 where ωp is the pth element of the ω vector (scalar) Solution : ωp(x,t)=ωp(x-λpt,0) n q(x,t)=rω(x,t)= ω p ( x λ p t, 0) r p 1 A linear weighted combination of eigenvectors Image from 1 10/57
11 Solution for linear Riemann problem(a) Reminder: q ( x, 0)= { q l for x <0 q r for x > 0 Initial values can be decomposed into a linear combination of the eigenvectors n ql= ω r 1 p l p n and qr= p p ω rr 1 For any point, the solution would be a combination of the eigenvalue decomposition of ql and qr. 11/57
12 Solution for linear Riemann problem(b) E.g. for a 2x2 system there are the following possibilities, depending on the position of x at time t. 1) q(x,t)=ωl1r1+ωl2r2 2)q(x,t)=ωl1r1+ωr2r2 3)q(x,t)=ωr1r1+ωr2r2 Images from [1] Consequently, the overall solution can be divided into three areas 12/57
13 Solution for linear Riemann problem(c) The previous solution can be interpreted as a series of jumps that connect the ql and qr states Image from [1] It can be seen that any two neighboring states differ only on the coefficient of one eigenvector. Also, for the difference between the leftmost and rightmost state the following must hold ql-qr=ωl1r1+ωl2r2 -ωr1r1+ωr2r2 = a1r1+a2r2 where a1=ω1l-ω1r and a2=ω2l-ω2r 13/57
14 Solution for linear Riemann problem(d) Generalizing ql-qr=ωl1r1+ωl2r2 -ωr1r1+ωr2r2 = a1r1+a2r2 ql-qr=ra a can be found initially by decomposing ql-qr Each solution domain can be found by starting from either qr or ql and adding the appropriate number of jumps depending on point coordinates and characteristic speed. q(x,t)=ql+ p a r p p p : λ < x/ t or equivalently q(x,t)=qr+ a pr p p p : λ > x/ t 14/57
15 Non linear system Shallow water equations. Conserved quantities: h (hρ:mass) and hu (huρ:momentum) They are derived from the mass and momentum equilibrium accounting for hydrostatic pressure due to wave height. q q= h = 1 uh q2 f (q)= hu 1 hu2 + gh2 2 for smooth solutions the normal conservation expression qt+f(q)x=0 can be re written as qt+f (q)qx=0 With f ' (q)= u + gh 2u 15/57
16 Solution of a non linear hyperbolic system Question: are there still characteristic curves on which the solution remains constant? d q( x (t ), t )=q t + X ' (t )q x =0 dt For smooth solutions the conservation law takes the form qt+f (q)qx=0 Combining the 2 it follows that f (q)qx=x (t)qx This is valid if X (t) is an eigenvalue of the flux Jacobian The slope of the characteristic curve is the eigenvalue BUT the eigenvalue is now q dependent! 16/57
17 Eigenvalues of the shallow water equation The eigenvectors and eigenvalues pairs of the flux jacobian can found 1 λ =u gh 2 λ =u+ gh 1 r= 2 r = 1 u- gh 1 u+ gh Which are q dependent. Therefore, lines exist on the x-t plane along which q is constant, BUT they are not parallel! 17/57
18 Linear and Non Linear Characteristics(b) In this case seeing the gradient we can deduce that λ1(ql)<λ1(qr) This leads to a rarefaction A continuous solution to the Riemann problem In this case seeing the gradient we can deduce that λ2(ql)>λ2(qr) This leads to a shock A discontinuous solution to the Riemann problem Image from [1] 18/57
19 Linear and Non Linear Characteristics(a) Images from [1] 19/57
20 The weak solution Discontinuities cant be handled by using the smooth conservation law Analytical solution is not always feasible (especially for larger systems) Solution: weak solution Idea :Instead of solving the equation exactly and require the solution to be continuously differentiable, loosen the constrain by solving for u(t,x) for which it holds that ( uφ t + f ( u) φ x dxdt )= 0 For every continuously differentiable function with compact support φ Cc1 20/57
21 Rankine Hugoniot Jump condition(a) As seen, the Riemann problem can result in discontinuities which must be resolved. { ul for x < λt ur for λt > 0 Definition: φ=φ(t,x) any continuously differentiable function with compact support Let Ω be the domain. Compact support means φ(t)=0 on Ω boundaries U ( x, 0)= ( Uφ t + f (U ) φ x dxdt )= 0 ( + ) ( Uφ t + f (U ) φ x dxdt )=0 Apply divergence theorem Ω - Ω + ( + ) (div(uφ,f(u)φ))dxdt =0 Ω - Ω + 21/57
22 Rankine Hugoniot Jump condition(b) Remark. Due to compact support, φ(x,0) 0 only on the interface The interface is the line x=λt ds the differential of the interface Image from [4] (along the line x=λt) n+ds=(λ,-1)dt and n-ds=(-λ,1)dt 22/57
23 Rankine Hugoniot Jump condition(c) Splitting the integrals along the interface normal n vds + n vds=0 + Ω - + Ω - Substituting [ λu u f (u )+f (u )]φ(t, λt )dt =0 For the above to hold for every φ it follows that λ (u u )=f (u ) f (u ) Rankine Hugoniot Condition 23/57
24 Admissibility conditions Weak solution doesn't provide unique solutions, therefore the following constrains can be adopted Vanishing viscosity ute+f(ue)=euexx which converges to u as ε 0 i.e. the second order term acts as a smoother Entropy inequality Call entropy a continuously differentiable function η R n R And entropy flux a q Rn R So that it holds Dη(u)Df(u)=Dq(u) It follows that if u=u(t,x) us a C1 solution of the system of conservation laws η(u)t+q(u)x=0 i.e. we require an additional conservation law to hold thus decreasing the weak solution s possible breach of conservation. Stability condition As seen before, a shock is split into two minor shocks connected by an intermediate state. In order for the shock to be preserved and not dissipate the hint shock must be faster than the front one f ( u* ) f (u- ) * u u - f (u + ) f (u * ) u+ u * 24/57
25 Addressing the strange formations The solution can t be decomposed as the summation of weighted eigenvectors (due to nonlinearity) Idea: Find an intermediate state. Image from [1] Transition:depending on the speed difference of the characteristics choose either Shock transition Rarefaction transition 25/57
26 Resolving the shock(a) Idea: use the Rankine-Hugoniot condition f(ql)-f(qr)=s(ql-qr) Consequently, the speed of the shock is s and can be found solving the above equation 26/57
27 Resolving the shock(b) Image from [1] The shock is resolved by separating the crashing regions by a line with slope=s as calculated by the Rankine Hugoniot condition 27/57
28 Resolving the shock(c) The family of possible q* that provide a shock can be calculated as follows: First, determine that at the interface between q and q* there is a shock. This can be done by taking the λp(qr) and comparing it to λp(ql) if at any characteristic it holds that λp(qr)<λp(ql), or vice versa, there is a shock formation. The Rankine Hugoniot condition must hold, so: s(q-q*)=f(q*)-f(q) where q* is the intermediate state and q the extreme state Specifically for the shallow water equation we have s(h-h*)=hu-h*u* s(hu-h*u*)=hu2-h*u*u*+g(h2-h*h*)/2 Solving for s and substituting we get a quadratic equation that has the following solution 28/57
29 Resolving the shock(d) g (h) h* u (h)=u± (h h *) 2 h* h Each sign choice corresponds to one family of characteristics. In our case the sign corresponds to the λ1 family and the + to the λ2 (the - family approaches q* tangent to the vector r1 and the + tangent to the vector r2) It can be easily determined how many, if any, are physically correct. (entropy condition) 29/57
30 The entropy condition (Lax) For the shock to be physically correct, it must hold that The characteristics of the shock must satisfy λp(ql)>s>λp(qr)>0 or λp(ql)<s<λp(qr)<0 λj(ql)<s λj(qr)<s for j<p λj(ql)>s λj(qr)>s for j>p with the eigenvalues being ordered as λ1<λ2<.λn 30/57
31 Resolving the rarefaction(a) An integral curve is a parametric curve that represents a specific solution to an ordinary differential equation or system of equations. If the differential equation is represented as a vector field or slope field, then the corresponding integral curves are tangent to the field at each point. 31/57
32 Resolving the rarefaction(b) It follows that if q(ξ) is an integral curve of the vector field rp, it holds that q (ξ)=a(ξ)rp(q(ξ)) Take the shallow water equation as example q1(ξ) =α(ξ) which is satisfied with q1(ξ)=α(ξ)ξ We will examine the parameter value a=1 for ease and readability of the formulas 32/57
33 Resolving the rarefaction(c) q2(ξ) =q2(ξ)/q1(ξ)- (gq1(ξ) Fixing a point on the curve (any point) and solving for q2(ξ) and r1 q2(ξ)=ξul+2ξ( (gq1*- (gq1(ξ)) 33/57
34 Resolving the rarefaction(d) Or in terms of actual physical quantities hu=hlul+2h( (gq1*- (gq1(ξ)) similar procedure is followed for r2 Plotting the curves for the SWE gives (with fixed point being ql or qr depending on the problem configuration Image from [2] 34/57
35 Resolving the rarefaction(e) hu=h*u*+2h( (gh*- (gh) this equation can be rewritten as hu+2h (gh=h*u*+2h (gh*) Remark: since the points can be anywhere on the curve w1(q)=u+2h (gh) is constant along the curve equally w2(q)=u-2h (gh) Constant functions on the integral curve of rp are called p-riemann invariants 35/57
36 Resolving the rarefaction(f) Since wp is constant along the rp integral curve, d p it follows that w (q (ξ ))=0 dξ Expanding and substituting q (ξ)=q(ξ)rp wprp=0 wp is the gradient with respect to q 36/57
37 Resolving the rarefaction(g) The rarefaction is a continuous solution, therefore a strong one. The conservation law must be fulfilled qt=q (ξ(x,t))ξt and qx=q (ξ(x,t))ξx Substituting into conservation law ξt+λp(q(ξ))ξx=0 d p p λ = λ q ' (ξ )=0 dξ 37/57
38 Resolving the rarefaction(h) The rarefaction is a continuous solution, therefore a strong one. The conservation law must be fulfilled qt=q (ξ(x,t))ξt and qx=q (ξ(x,t))ξx Substituting into conservation law ξt+λp(q(ξ))ξx=0 The variation of λ along the curve is given by d p p λ = λ q ' (ξ )=0 dξ 38/57
39 Resolving the rarefaction(i) A centered rarefaction wave is a special wave, were the solution is constant along rays through the origin. The solution for this wave is self similar i.e. q(x,t)=q(x/t) in this case, x/t=ξ (ξ speed of propagation) consequently ξ=λp(q(ξ)) 39/57
40 Resolving the rarefaction(almost done) ξt+λp(q(ξ))ξx=0 in this case becomes ξ=λp(q(ξ)) Differentiate w.r.t. ξ 1= λ q ' (ξ )=0 Substitute q (ξ)=a(ξ)rp(q(ξ)) p p 1= λ ar p Solving for a and substituting we get rp q ' (ξ )= p p λ r 40/57
41 Resolving the rarefaction(last one I swear) rp q ' (ξ )= p p λ r Can be solved. As boundary condition we use the known state which the rarefaction connects to the intermediate state Image from [1] 41/57
42 Finding the intermediate state As discussed above, we can find a series of intermediate points that connect the left and right state with shocks and rarefactions. To find the correct intermediate point 2 conditions must hold It must the states in a physically correct way It must connect both states Solution: Find the common point that fulfills both those conditions (analytically or numerically (or graphically)) 42/57
43 Hugoniot loci The intermediate state must satisfy both the conditions coming from the left and from the right. All possible states that connect ql q* and qr are plotted. Admissibility conditions must be satisfied or the state discarded (this includes expected transition form or unphysical behavior, eg negative height) q* is chosen as as the only state satisfying all the above conditions Image from [1] 43/57
44 Structure of the rarefaction The rarefaction can be seen as a set of intermediate points from the intermediate point to the state it connects to. Therefore, structure of the rarefaction is given by the the intermediate states between q* and the side it connects to. Images from [1] 44/57
45 General strategy for solving the non linear 2x2 Riemann problem Determine the nature of the wave (shock or rarefaction) either by direct eigenvalue comparison or by an appropriate entropy condition Determine the intermediate state between the two waves Determine the structure of any rarefaction wave 45/57
46 Remark It might happen that a shock and rarefaction collide Images from [5] It is treated precisely as if two characteristic lines collided, i.e. with a shock or a rarefaction (in this example shock) that fulfills all the necessary conditions (Rankine Hugoniot, entropy) 46/57
47 Image from [5] 47/57
48 Dam break problem (a) Following the above procedure we can approach the dam breaking problem A system of shallow water equations with initial conditions { q ( x, 0)= h=hl u=0 for x <0 h=h r u=0 for x > 0 Assuming hl=6 hr=4 As discussed previously Consequently λ2(ql)>λ2(qr)>0>λ1(qr)>λ1(ql) 1 λ =u gh 2 and λ =u+ gh λ1(ql)<λ1(qr)<0 the left state is connected to the intermediate via a rarefaction 0<λ2(ql)<λ2(qr) the right state is connected to the intermediate via a shock 48/57
49 Dam break problem (b) The solution has the form As discussed above, we can use Image from [1] the shock and rarefaction equations to create families of intermediate states. From the resulting curves we only need to take into account the 1 rarefaction that connects the intermediate with ql and the 2 shock that connects qr with the intermediate 49/57
50 Dam break problem (c) Rarefaction connects the 1 characteristics and ql with q* therefore we use the function q ' (ξ )= with 1 r= 1 u- gh 1 λ =u gh 1 r 1 1 λ r 1 λ = u/h- 0.5 g/h 1/hu 50/57
51 Dam break problem (d) 1 r q ' (ξ )= 1 1 λ r Substituting into h=(a-ξ)2/9g with A the integration constant The left side of the rarefaction is known (q=ql) at ξ=ul- ghl A=ul+ ghl Equally, all the intermediate points must satisfy the same conditions This way using h=(ul+ ghl-ξ)2/9g we can find the intermediate points Alternatively, we can use the Riemann invariants and determine the points from A=ul+ ghl=u*+ gh* 51/57
52 Dam break problem (e) Shock connects the 2 characteristics and q* with qr therefore we use this function g (h) h * u (h)=u + (h h* ) 2 h* h the entropy condition excludes all the unphysical points Finding the common intermediate we get Images from [1] 52/57
53 Dam break problem (f) And the physical behavior Images from [1] 53/57
54 Different initial conditions Depending on the initial conditions the possible solutions are qualitatively Two shocks Two rarefactions Shock-rarefaction Rarefaction-shock Determined with the procedure outlined in the previous slides 54/57
55 Literature [1] Lectures of subject tsunami simulation at [2] LeVeque, R. (2002). Finite Volume Methods for Hyperbolic Problems (Cambridge Texts in Applied Mathematics) [3] Hyperbolic Conservation Laws An Illustrated Tutorial, Lectures Notes [4] Modelling and Optimization of Flows on Networks Cetraro School 2007 [5] Riemann Problems of the Shallow Water Equations at 55/57
56 56/57
57 57/57
Notes: 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 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 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 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 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 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 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 informationAMath 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, < 1, 1 +, 1 < < 0, ψ() = 1, 0 < < 1, 0, > 1, so that it verifies ψ 0, ψ() = 0 if 1 and ψ()d = 1. Consider (ψ j ) j 1 constructed as
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More 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 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 informationThe shallow water equations Lecture 8. (photo due to Clark Little /SWNS)
The shallow water equations Lecture 8 (photo due to Clark Little /SWNS) The shallow water equations This lecture: 1) Derive the shallow water equations 2) Their mathematical structure 3) Some consequences
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems 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, 2014 1 / 56 We study some basic
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 informationHyperconcentrated 1D Shallow Flows on Fixed Bed with Geometrical Source Term Due to a Bottom Step
J Sci Comput (2011) 48:319 332 DOI 10.1007/s10915-010-9457-z Hyperconcentrated 1D Shallow Flows on Fixed Bed with Geometrical Source Term Due to a Bottom Step L. Begnudelli G. Rosatti Received: 24 November
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 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 informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationA 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 informationComputational Fluid Dynamics. PHY 688: Numerical Methods for (Astro)Physics
Computational Fluid Dynamics Hydrodynamics When we discussed PDEs, we focused so far on scalar PDEs Often we wish to study systems of PDEs. Here we'll look at the equations of hydrodynamics Nonlinear system
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 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 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 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 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 information1 Basic Second-Order PDEs
Partial Differential Equations A. Visintin a.a. 2011-12 These pages are in progress. They contain: an abstract of the classes; notes on some (few) specific issues. These notes are far from providing a
More informationHyperbolicity singularities in rarefaction waves
Hyperbolicity singularities in rarefaction waves Alexei A. Mailybaev and Dan Marchesin Abstract For mixed-type systems of conservation laws, rarefaction waves may contain states at the boundary of the
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 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 informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationNumerische Mathematik
Numer. Math. (2007 106:369 425 DOI 10.1007/s00211-007-0069-y Numerische Mathematik Hyperbolic balance laws: Riemann invariants and the generalized Riemann problem Matania Ben-Artzi Jiequan Li Received:
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 informationApplication of the relaxat ion met hod to model hydraulic jumps
Application of the relaxat ion met hod to model hydraulic jumps P. J. Montgomery Mathematics and Computer Science Program, University of Northern British Columbia, Prince George, Canada. Abstract A finite
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 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 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 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 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 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 informationOn a simple model of isothermal phase transition
On a simple model of isothermal phase transition Nicolas Seguin Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie Paris 6 France Micro-Macro Modelling and Simulation of Liquid-Vapour Flows
More informationE = where γ > 1 is a constant spesific to the gas. For air, γ 1.4. Solving for p, we get. 2 ρv2 + (γ 1)E t
. The Euler equations The Euler equations are often used as a simplification of the Navier-Stokes equations as a model of the flow of a gas. In one space dimension these represent the conservation of mass,
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 informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationAdaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations
Adaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations Thida Pongsanguansin, Khamron Mekchay and Montri Maleewong Abstract The adaptive Discontinuous Galerkin DG method for solving
More 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 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 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 informationWELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL
WELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL VIVIEN DESVEAUX, MARKUS ZENK, CHRISTOPHE BERTHON, AND CHRISTIAN KLINGENBERG Abstract. The present paper concerns the derivation
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 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 informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationConstrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics
Project work at the Department of Mathematics, TUHH Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Katja Baumbach April 4, 005 Supervisor:
More informationIV. Compressible flow of inviscid fluids
IV. Compressible flow of inviscid fluids Governing equations for n = 0, r const: + (u )=0 t u + ( u ) u= p t De e = + ( u ) e= p u+ ( k T ) Dt t p= p(, T ), e=e (,T ) Alternate forms of energy equation
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 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 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 informationChapter 4: Non-Linear Conservation Laws; the Scalar Case. The conservation law in eqn. (4.1) will be hyperbolic if its eigenvalue ( ) ( )
Chapter 4: Non-inear Conservation aws; the Scalar Case 4.1) Introduction In the previous chapter we developed an understanding of monotonicity preserving advection schemes and iemann solvers for linear
More informationCranfield ^91. College of Aeronautics Report No.9007 March The Dry-Bed Problem in Shallow-Water Flows. E F Toro
Cranfield ^91 College of Aeronautics Report No.9007 March 1990 The Dry-Bed Problem in Shallow-Water Flows E F Toro College of Aeronautics Cranfield Institute of Technology Cranfield. Bedford MK43 OAL.
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 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 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 informationDIMINISHING FUNCTIONALS FOR NONCLASSICAL ENTROPY SOLUTIONS SELECTED BY KINETIC RELATIONS
DIMINISHING FUNCTIONALS FOR NONCLASSICAL ENTROPY SOLUTIONS SELECTED BY KINETIC RELATIONS MARC LAFOREST AND PHILIPPE G. LEFLOCH Abstract. We consider nonclassical entropy solutions to scalar conservation
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
More informationS. N. Kruzhkov s lectures on first-order quasilinear PDEs
to appear in: Analytical and Numerical Aspects of PDEs de Gruyter 2009 S. N. Kruzhkov s lectures on first-order quasilinear PDEs Gregory A. Chechkin and Andrey Yu. Goritsky Abstract. The present contribution
More informationBefore we look at numerical methods, it is important to understand the types of equations we will be dealing with.
Chapter 1. Partial Differential Equations (PDEs) Required Readings: Chapter of Tannehill et al (text book) Chapter 1 of Lapidus and Pinder (Numerical Solution of Partial Differential Equations in Science
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 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 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 informationHysteresis rarefaction in the Riemann problem
Hysteresis rarefaction in the Riemann problem Pavel Krejčí 1 Institute of Mathematics, Czech Academy of Sciences, Žitná 25, 11567 Praha 1, Czech Republic E-mail: krejci@math.cas.cz Abstract. We consider
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 informationSalmon: Lectures on partial differential equations
4 Burger s equation In Lecture 2 we remarked that if the coefficients in u x, y,! "! "x + v x,y,! "! "y = 0 depend not only on x,y but also on!, then the characteristics may cross and the solutions become
More informationRotation & nonlinear effects in shallow water: the Rossby adjustment problem
Rotation & nonlinear effects in shallow water: the Rossby adjustment problem Ross Tulloch April 29, 2004 1 1 Shallow water equations Consider, as in [6], a homogeneous incompressible 2D fluid in a channel
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 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 informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationCausality, hyperbolicity, and shock formation in Lovelock theories
Causality, hyperbolicity, and shock formation in Lovelock theories Harvey Reall DAMTP, Cambridge University HSR, N. Tanahashi and B. Way, arxiv:1406.3379, 1409.3874 G. Papallo, HSR arxiv:1508.05303 Lovelock
More informationArchimedes Center for Modeling, Analysis & Computation. Singular solutions in elastodynamics
Archimedes Center for Modeling, Analysis & Computation Singular solutions in elastodynamics Jan Giesselmann joint work with A. Tzavaras (University of Crete and FORTH) Supported by the ACMAC project -
More informationThe Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow
The Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow D.W. Schwendeman, C.W. Wahle 2, and A.K. Kapila Department of Mathematical Sciences, Rensselaer Polytechnic
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More information