Notes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
|
|
- Percival Harrell
- 5 years ago
- Views:
Transcription
1 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 Reading: Chapter 3 Recall: Some slides have section numbers on footer. $CLAW/book Examples from the book. Gallery of applications. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 Diffusive flux q(x, t) concentration β diffusion coefficient (β > 0) diffusive flux βq x (x, t) q t + f x 0 diffusion equation: q t (βq x ) x βq xx (if β const). Heat equation: Same form, where q(x, t) density of thermal energy κt (x, t), T (x, t) temperature, κ heat capacity, flux βt (x, t) (β/κ)q(x, t) q t (x, t) (β/κ)q xx (x, t). R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.2 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.2 Advection-diffusion q(x, t) concentration that advects with velocity u and diffuses with coefficient β: Advection-diffusion equation: flux uq βq x. q t + uq x βq xx. If β > 0 then this is a parabolic equation. Advection dominated if u/β (the Péclet number) is large. Fluid dynamics: parabolic terms arise from thermal diffusion and diffusion of momentum, where the diffusion parameter is the viscosity. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.2 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.2
2 The Riemann problem The Riemann problem consists of the hyperbolic equation under study together with initial data of the form { ql if x < 0 q(x, 0) q r if x 0 Piecewise constant with a single jump discontinuity from q l to q r. The Riemann problem is fundamental to understanding The mathematical theory of hyperbolic problems, Godunov-type finite volume methods Why? Even for nonlinear systems of conservation laws, the Riemann problem can often be solved for general q l and q r, and consists of a set of waves propagating at constant speeds. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8 The Riemann problem for advection The Riemann problem for the advection equation q t + uq x 0 with { ql if x < 0 q(x, 0) q r if x 0 has solution { ql if x < ut q(x, t) q(x ut, 0) if x ut consisting of a single wave of strength W q r q l propagating with speed s u. q r R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8 Riemann solution for advection q(x, T ) x t plane q(x, 0) R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8
3 Discontinuous solutions Note: The Riemann solution is not a classical solution of the PDE q t + uq x 0, since q t and q x blow up at the discontinuity. Integral form: d x2 q(x, t) dx uq(x, t) uq(x 2, t) dt x Integrate in time from t to t 2 to obtain x2 x q(x, t 2 ) dx t2 t x2 x uq(x, t) dt q(x, t ) dx t2 t uq(x 2, t) dt. The Riemann solution satisfies the given initial conditions and this integral form for all x 2 > x and t 2 > t 0. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.7 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.7 Discontinuous solutions Vanishing Viscosity solution: The Riemann solution q(x, t) is the limit as ɛ 0 of the solution q ɛ (x, t) of the parabolic advection-diffusion equation q t + uq x ɛq xx. For any ɛ > 0 this has a classical smooth solution: R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec..6 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec..6 Diagonalization of linear system Consider constant coefficient linear system q t + Aq x 0. Suppose hyperbolic: Real eigenvalues λ λ 2 λ m, Linearly independent eigenvalues r, r 2,..., r m. Let R r r 2 r m m m matrix of eigenvectors. Then Ar p λ p r p means that AR RΛ where λ λ 2 Λ... diag(λ, λ 2,..., λ m ). λ m AR RΛ A RΛR and R AR Λ. Similarity transformation with R diagonalizes A. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.9 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.9
4 Diagonalization of linear system Consider constant coefficient linear system q t + Aq x 0. Multiply system by R : R q t (x, t) + R Aq x (x, t) 0. Introduce RR I: R q t (x, t) + R ARR q x (x, t) 0. Use R AR Λ and define w(x, t) R q(x, t): w t (x, t) + Λw x (x, t) 0. Since R is constant! This decouples to m independent scalar advection equations: w p t (x, t) + λp w p x(x, t) 0. p, 2,..., m. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.9, 3. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.9, 3. Solution to Cauchy problem Suppose q(x, 0) q(x) for < x <. From this initial data we can compute data w (x) R q(x) The solution to the decoupled equation w p t + λp wx p 0 is w p (x, t) w p (x λ p t, 0) w p (x λ p t). Putting these together in vector gives w(x, t) and finally q(x, t) Rw(x, t). We can rewrite this as m m q(x, t) w p (x, t) r p w p (x λ p t) r p p p R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3. Linear acoustics Example: Linear acoustics in a d gas tube p p(x, t) pressure perturbation q u u(x, t) velocity Equations: p t + κu x 0 Change in pressure due to compression ρu t + p x 0 Newton s second law, F ma where K bulk modulus, and ρ unperturbed density of gas. Hyperbolic system: p 0 κ p + 0. u /ρ 0 u t x R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.9.
5 Eigenvectors for acoustics u0 K A 0 /ρ 0 u 0 (acoustics relative to flow with speed u 0 ) Eigenvectors: r ρ0 c 0, r 2 ρ0 c 0. Check that Ar p λ p r p, where λ u 0 c 0, λ 2 u 0 + c 0. with c 0 K 0 /ρ 0 K 0 ρ 0 c 2 0. Note: Eigenvectors are independent of u 0. Let Z 0 ρ 0 c 0 K 0 ρ 0 impedance. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.8 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 2.8 Physical meaning of eigenvectors Eigenvectors for acoustics: r ρ0 c 0 Z0, r 2 ρ0 c 0 Z0. Consider a pure -wave (simple wave), at speed λ c 0, If q (x) q + w (x)r then q(x, t) q + w (x λ t)r Variation of q, as measured by q x or q q(x + x) q(x) is proportional to eigenvector r, e.g. q x (x, t) w x(x λ t)r R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.4, 3.5 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.4, 3.5 Physical meaning of eigenvectors Eigenvectors for acoustics: r ρ0 c 0 Z0, r 2 ρ0 c 0 Z0. In a simple -wave (propagating at speed λ c 0 ), px Z0 β(x) u x The pressure variation is Z 0 times the velocity variation. Similarly, in a simple 2-wave (λ 2 c 0 ), px Z0 β(x) u x The pressure variation is Z 0 times the velocity variation. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5
6 Acoustic waves q(x, 0) p (x) 0 p (x) 2Z 0 Z0 + p (x) Z0 2Z 0 w (x, 0)r + w 2 (x, 0)r 2 p (x)/2 p (x)/(2z 0 ) + p (x)/2 p (x)/(2z0 ). R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5 Solution by tracing back on characteristics The general solution for acoustics: q(x, t) w (x λ t, 0)r + w 2 (x λ 2 t, 0)r 2 w (x + c 0 t, 0)r + w 2 (x c 0 t, 0)r 2 Recall that w(x, 0) R q(x, 0), i.e. w (x, 0) l q(x, 0), w 2 (x, 0) l 2 q(x, 0) where l and l 2 are rows of R. R l l 2 Note: l and l 2 are left-eigenvectors of A: l p A λ p l p since R A ΛR. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5 Solution by tracing back on characteristics The general solution for acoustics: q(x, t) w (x λ t, 0)r + w 2 (x λ 2 t, 0)r 2 w (x + c 0 t, 0)r + w 2 (x c 0 t, 0)r 2 (x, t) x λ 2 t x c 0 t x λ t x + c 0 t R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.6 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.6
7 Solution by tracing back on characteristics The general solution for acoustics: q(x, t) w (x λ t, 0)r + w 2 (x λ 2 t, 0)r 2 q(x, t) w 2 constant w constant w 2 (x λ 2 t, 0) l 2 q(x λ 2 t, 0) w (x λ t, 0) l q(x λ t, 0) R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.5 Riemann Problem Special initial data: q(x, 0) { ql if x < 0 q r if x > 0 Example: Acoustics with bursting diaphram (u l u r 0) Pressure: Acoustic waves propagate with speeds ±c. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec Riemann Problem for acoustics Waves propagating in x t space: Left-going wave W q m q l and right-going wave W 2 q r q m are eigenvectors of A. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.8
8 Riemann Problem for acoustics In x t plane: q(x, t) w (x + ct, 0)r + w 2 (x ct, 0)r 2 Decompose q l and q r into eigenvectors: q l wl r + wl 2 r2 q r wrr + wrr 2 2 Then q m w rr + w 2 l r2 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.9 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.9 Riemann Problem for acoustics q l wl r + wl 2 r2 q r wrr + wrr 2 2 Then q m wrr + wl 2 r2 So the waves W and W 2 are eigenvectors of A: W q m q l (wr wl )r W 2 q r q m (wr 2 wl 2 )r2. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.9 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.9 Riemann solution for a linear system Linear hyperbolic system: q t + Aq x 0 with A RΛR. General Riemann problem data q l, q r lr m. Decompose jump in q into eigenvectors: q r q l m α p r p p Note: the vector α of eigen-coefficients is α R (q r q l ) R q r R q l w r w l. Riemann solution consists of m waves W p lr m : W p α p r p, propagating with speed s p λ p. R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.9 R.J. LeVeque, University of Washington IPDE 20, June 22, 20 FVMHP Sec. 3.9
Conservation 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 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 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 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 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 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 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 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 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 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 informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
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 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 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 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 informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationTraffic Flow Problems
Traffic Flow Problems Nicodemus Banagaaya Supervisor : Dr. J.H.M. ten Thije Boonkkamp October 15, 2009 Outline Introduction Mathematical model derivation Godunov Scheme for the Greenberg Traffic model.
More 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 informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
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 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 informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationThe 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 informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
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 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 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 informationPhase Plane Behavior of Solitary Waves in Nonlinear Layered Media
Phase Plane Behavior of Solitary Waves in Nonlinear Layered Media Randall J. LeVeque 1 and Darryl H. Yong 2 1 Department of Applied Mathematics, University of Washington, Box 35242, Seattle, WA 98195-242
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 informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
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 informationAdvection, Conservation, Conserved Physical Quantities, Wave Equations
EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary
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 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 informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
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 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 informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More 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 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 informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More 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 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 informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
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 informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
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 informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
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 information25 Self-Assessment. 26 Eigenvalue Stability for a Linear ODE. Chapter 7 Eigenvalue Stability
Chapter 7 Eigenvalue Stability As we have seen, while numerical methods can be convergent, they can still exhibit instabilities as the number of timesteps n increases for finite t. For example, when applying
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
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 informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationDUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs
DUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs MATH 4354 Fall 2005 December 5, 2005 1 Duhamel s Principle
More informationMath 220A - Fall 2002 Homework 5 Solutions
Math 0A - Fall 00 Homework 5 Solutions. Consider the initial-value problem for the hyperbolic equation u tt + u xt 0u xx 0 < x 0 u t (x, 0) ψ(x). Use energy methods to show that the domain of dependence
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationLECTURE 1 INTRODUCTION TO NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS I. NONLINEAR DIFFUSION EQUATION
LECTURE 1 INTRODUCTION TO NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS I. NONLINEAR DIFFUSION EQUATION UGUR G. ABDULLA I am going to start a series of lectures on nonlinear partial differential equations.
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationLu u. µ (, ) the equation. has the non-zero solution
MODULE 18 Topics: Eigenvalues and eigenvectors for differential operators In analogy to the matrix eigenvalue problem Ax = λx we shall consider the eigenvalue problem Lu = µu where µ is a real or complex
More informationMath 124A October 11, 2011
Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
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 informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationReading: P1-P20 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering)
Chapter 1. Partial Differential Equations Reading: P1-P0 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering) Before even looking
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 informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More 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 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 informationContents of lecture 2b. Lectures 2a & 2b. Physical vs. computational coordinates [2] Physical vs. computational coordinates [1]
Contents of lecture b Lectures a & b P. A. Zegeling Mathematical Institute Utrecht University The Netherlands Parameter-free non-singular moving grids in D: Theory & properties Application to resistive
More informationVibrating Strings and Heat Flow
Vibrating Strings and Heat Flow Consider an infinite vibrating string Assume that the -ais is the equilibrium position of the string and that the tension in the string at rest in equilibrium is τ Let u(,
More informationOPAC102. The Acoustic Wave Equation
OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating
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 informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationConvergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates
Convergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates Xin Wen and Shi Jin Abstract We study the L 1 -error estimates
More informationNotes on theory and numerical methods for hyperbolic conservation laws
Notes on theory and numerical methods for hyperbolic conservation laws Mario Putti Department of Mathematics University of Padua, Italy e-mail: mario.putti@unipd.it January 19, 2017 Contents 1 Partial
More informationNUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY
NUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY Stefan Johansson Department of Information Technology Scientific Computing
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 informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
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 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 informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationApproximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD
Approximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD Andrea Mignone Collaborators: G. Bodo, M. Ugliano Dipartimento di Fisica Generale, Universita di Torino (Italy)
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
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 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 informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
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 informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
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 information