Info. No lecture on Thursday in a week (March 17) PSet back tonight
|
|
- Shon Fisher
- 5 years ago
- Views:
Transcription
1 Lecture
2 Info No lecture on Thursday in a week (March 7) PSet back tonight
3 Nonlinear transport & conservation laws What if transport becomes nonlinear?
4 Remember: Nonlinear transport A first attempt at understanding what is going on comes from analyzing the characteristic lines Simplest example: lnviscid Burgers' equation We can use the implicit ansatz du dx u(x, t) =u(x ut, 0) Characteristics for initial conditions all meet in (x,t)=(,)! u(x, 0) = x Solution is not defined at (x,t)=(,) => so what do we do?
5 Nonlinear transport & conservation laws A more famous example of issues with characteristics is the Riemann problem: Take again Burger s equation: Riemann problem: How does solution evolve in time for initial conditions du dx u(x, 0) = A, x < 0, and u(x, 0) = B, x 0 t t A > B: shock A < B: fan
6 How to find shock wave or fan solution? Physically, we could add a bit of diffusion to avoid discontinuities, then let D -> 0 Not very feasible numerically see Lect. Note: There is an example in the book (analytical) for u t + uu x = u x x
7 Conservation laws Burger s equation creates discontinuities Mathematical description using integral form (=conservation law) In general, conservation laws read du dx u t + f(u, u x,...) x =0() d dt Z b a udx + f(u(b),u x (b),...) f(u(a),u x (a),...)=0 Special form above: Diffusion is also a conservation law
8 Weak solutions Integral form is not very handy, because we have to show it holds for all times and all intervals Instead: Use the weak form: Multiply PDE with a test function ɸ and integrate over time and space. Assuming ɸ has compact support (almost everywhere zero), integration by parts gives (for Burger s equation): u t + f x (u) =0, Z 0 Z [ tu + x f(u)]dxdt = Z Lecture (x, 0)u(x, 0)dx For Burgers eq: One can show (PSet 2) that the shock in the Riemann problem with initial condition u(x, 0) = ( u L x<0 u R x>0 with ul>ur ul t ur here: ul=0, ur =- is solved by: u(x, t) = ( u L u R x<st x>st with shock speed s = u L + u R 2 ul ur
9 RH jump condition The previous example was specific for Burger s equation. But there is a general condition that determines the shock speed. Rankine-Hugoniot (RH) jump condition: The shock speed s of a conservation law with flux f is s = f(u L) f(u R ) = [f] u L u R [u] Lecture u L, u R : value of u on the left & right of the shock Also valid for initial conditions that are not piecewise linear. Then u L, u R : value of u infinitesimally close to shock
10 Rarefaction waves How does the Riemann problem look like when there are no shocks? ul t ur Correct (weak) solution (for Burgers eq): Rarefaction wave: u(x, t) = 8 >< x t >: u L u R x<u L t u L t apple x apple u R t x>u R t ul ur A < B: fan ul ur Think of cars starting to drive at green light Unfortunately, a shock would also be a weak solution in this case!!! We need another criterium to decide for a shock vs. a fan
11 Entropy condition Note that characteristics must collide for a shock to occur, while they separate for fans. This can be formulated in the entropy condition: If f 0 (u L ) >s>f 0 (u R ) => shocks otherwise: fans with u(x,t) = x/t
12 A first numerical try Simple idea: Discretize flux, not PDE! du d du du Start from f (u) = 0 and not from - + f '(u) - = 0 at dx at dx E.g. upwind finite differences (assuming f > 0) Upwind (for f '>o) Uj,n+l - Uj,n + f (Uj+l,n)- f (Uj,n) = 0 At Ax The advantage of this approach (apart from allowing discont. u) is that the conservation law is automatically fulfilled: B Discrete conservation uj,n+l - uj, - FA,^) = 0 j=a Lecture
13 Numerical flux function In fact, we have just seen that any scheme of the form U j,n+ t U j,n = x [F (U j p,n,u j p+,n,...,u j+q,n ) F (U j p+,n,u j p+2,n,...,u j+q+,n )] is a conservative scheme (meaning total U inside a domain is determined only by flux F at boundaries) Of course we cannot freely combine values of U into a flux: Consistency The flux F (Uj+,,..., Uj-,) satisfies F (u,..., u) = f (u) Question: How to choose flux F in the best way
14 A first numerical try Upwind finite differences (assuming f > 0) Upwind (for f '>o) Uj,n+l - Uj,n + f (Uj+l,n)- f (Uj,n) = 0 At Ax Using F j,n = f(u j,n ), this can be written as U j,n+ t U j,n = x [F j+,n F j,n ] We can extend this scheme to capture both flux directions: ( f(u j,n ) f 0 (U j,n ) > 0 F j,n = f(u j,n ) f 0 (U j,n ) < 0 First-order method, but fails for certain initial conditions (see LeVeque. p. 34)
15 Lecture 8.086
16 Last time: Conservation laws Burger s equation creates discontinuities Mathematical description using integral form (=conservation law) In general, conservation laws read du dx u t + f(u, u x,...) x =0() d dt Z b a udx + f(u(b),u x (b),...) f(u(a),u x (a),...)=0 Special form above: Diffusion is also a conservation law
17 Last time: RH jump condition The previous example was specific for Burger s equation. But there is a general condition that determines the shock speed. Rankine-Hugoniot (RH) jump condition: The shock speed s of a conservation law with flux f is s = f(u L) f(u R ) = [f] u L u R [u] Lecture u L, u R : value of u on the left & right of the shock Also valid for initial conditions that are not piecewise linear. Then u L, u R : value of u infinitesimally close to shock
18 Last time: Conservation laws Simple idea: Discretize flux, not PDE! du d du du Start from f (u) = 0 and not from - + f '(u) - = 0 at dx at dx E.g. upwind finite differences (assuming f > 0) Upwind (for f '>o) Uj,n+l - Uj,n + f (Uj+l,n)- f (Uj,n) = 0 At Ax Using F j,n = f(u j,n ), this can be written as U j,n+ t U j,n = x [F j+,n F j,n ] Advantage: conservation law is automatically fulfilled! We can extend this scheme to capture both flux directions: F j,n = ( f(u j,n ) f 0 (U j,n ) > 0 f(u j,n ) f 0 (U j,n ) < 0
19 Lax-Friedrichs for nonlinear conservation laws Remember Lax-Friedrichs for linear -way wave eq: Lax- Friedrichs uj,n+l - 5(uj+l,n uj-,n) Uj+l,n - Uj-l,n = C At 2Ax The nonlinear, conservative formulation is: This can be written in the standard discrete conservative form (see before) U j,n+ t U j,n = x [F j+,n F j,n ] by introducing the LF flux Lax-Friedrichs flux Ax F~~ 5 = [f(u~) - 2at(U3+l - UJ) j+ + f (u3+l)]
20 Lax-Friedrichs for nonlinear conservation laws As in the linear case, Lax-Friedrichs is only first-order accurate. It captures both directions, but the central difference flux term creates extra diffusion! We could just try to use Lax-Wendroff (which was second order), but this creates oscillations Schemes that avoid oscillations are called total variation diminishing (TVD) methods Define total variation (TV): TV(u) = / dx TV(U) = x IUJ+I - U,I O" du A TVD method has the property TV(Un+l) < TV(Un)
21 Godunov s method Some problems so far: Upwind is good, but we need to know direction (i.e. velocity) of wave. Problem: System of equations! Lax-Friedrichs has no upwind direction problems, but high diffusion and thus low accuracy 2nd order methods (e.g. LW) are prone to oscillations. How can we find schemes (in particular numerical flux functions) that are high order, direction-independent and TVD? Godunov s method: A systematic way to construct TVD numerical flux functions (more details in Strang s book and in R.J. LeVeque, Numerical Methods for Conservation Laws) But: Godunov s theorem: Linear numerical schemes for solving partial differential equations (PDE s) that are TVD can be at most first-order accurate.
22 Flux-limiting methods Linear in Godunov s theorem means a scheme where the scheme itself (flux function, coefficients etc.) are constant, i.e do not change their form depending on the current solution. => Use nonlinear numerical scheme, e.g. flux-limiter methods Idea: Combine high-order flux F H (good for smooth part) with low-order flux F L (TVD, good for shock) F(U;j) = FL(U;j) + 4'(U;j) [FH(U;j) - FL(U;j)] Choose Φ close to 0 near shock, near away from shock
23 Flux-limiter for Lax-Wendroff Remember Lax-Wendroff for -way wave eq. u t +au x =0 U j,n+ t U j,n = a U j+,n U j,n 2 x Can write it like an upwind plus correction method: U j,n+ = U j,n r(u j,n U j,n ) + t U 2 a2 j+,n 2U j,n + U j,n x 2 2 r( r)(u j+,n 2U j,n + U j,n ) The correction term increases accuracy (next term in series expansion), but we have seen it leads to oscillations Lecture Corresponding flux: F (U, j) =au j + 2 a( r)(u j+ U j ) Then, LW can be written as U j,n+ t U j,n = x [F j+,n F j,n ] (F j+,n = F(U,j+) etc.)
24 Flux-limiter for Lax-Wendroff Lax-Wendroff flux (prev. slide): F (U, j) =au j + 2 a( r)(u j+ U j ) This flux is an upwind flux (au j ) corrected with a high-order flux! We replace F(U,j) by F (U, j) =au j + 2 a( r)(u j+ U j ) j How should we choose? j Naturally, it should depend on the quantity It can be shown: No oscillations for j = U j U j U j+ U j close to for U smooth far from for shock in U 0 apple ( ) apple 2 and 0 apple ( ) apple 2 for all Θ For nonlinear equations: replace speed a by the local speed (see LeVeque, p. 8)
25 Flux-limiter for Lax-Wendroff LW-flux with flux-limiter: F (U, j) =au j + 2 a( r)(u j+ U j ) j Choices for ( ) φ(θ) 2t LW is not total vari tion diminis k, Θ
26 Finite volume method Another interpretation of the discrete conservation law: Finite volume method Directly approximates the integral form: - (j+i)~x Cell average Uj,n approximates uj,, = - / u(x, nat) dx AX (,-+)AX Change in I cell average '' U(X, (j++)~x (3++)Ax dt (J-;)*x ) dx + [f(uj++,n) - i(uj-+,n)] = 0 Integrate in time from n to n+ to obtain the finite volume conservation law /' U(X, (j-+)ax LA, tn+~) (n+l)at dx - l-;)ax (3+4)Ax U(X, tn) dx + (n+l)at f (uj+;,t) dt- LA, f (u,-;, t) dt = 0
27 Finite volume method - /(j+i)~x Using uj,, = - u(x, nat) dx the integral form (j++)~x tn+~) /' U(X, (j-+)ax LA, (n+l)at can be written as t (ū j,n+ AX (,-+)AX dx - l-;)ax (3+4)Ax U(X, tn) dx + (n+l)at f (uj+;,t) dt- LA, f (u,-;, t) dt = 0 ū j,n )+ t x Z (n+) n t t f(u j+ 2,t)dt t x Z (n+) n t t f(u j,t)dt =0 2 Using the time-averaged flux f j+ 2,n = t Z (n+) n t t f(u j+ 2,t)dt this becomes t (ū j,n+ ū j,n )+ x fj+ 2,n fj 2,n =0
28 Finite volume method Treat ū j,n as independent variables U j,n, approximate average fluxes by F j+ : t (U j,n+ U j,n )+ x F j+ 2,n F j 2,n =0 2,n Easiest choice for numerical flux F: Assume f is constant over one time step, ie. F j+ 2,n = f(u j+ 2,n t ) Note: We need to interpolate U to half-grid points! Note: Other choices of F can lead to implicit schemes
29 Why finite volumes? Grid can be arbitrary, conservation law always fulfilled y j+ NW N NE y j WW W nw w n P ne e n e E y EE y j- sw s se y SW x S SE j i x x i- x i x i+
30 AREPO finite volume code (Volker Springel) Example
31 Navier-Stokes / Flow problems Tricky equation, tricky numerics Incompressible Navier-Stokes equation: du -+(u-v)u=-vp+-nu+ dt Re divu=v=u=o f momentum balance continuity/incompr. equation In this dimensionless form, the only physical parameter is the Reynolds number Re inertial forces (a velocity U) (a length L) Reynoldsnumber Re= % viscous forces (kinematic viscosity v) - Pressure p is a Lagrange multiplier, adjusted such that the incompressibility condition is fulfilled.
32 Example: 2D flow in a lid-driven cavity Consider 2D box, with no-slip boundary conditions at three boundaries: = [0,l x ] [0,l y ]. The domain is In 2D, we can write the NS equations in component form as u t + p x = (u 2 ) x (uv) y + Re (u xx + u yy ) v t + p y = (uv) x (v 2 ) y + Re (v xx + v yy ) Lecture u x + v y = 0 The boundary conditions will be u(x, l y )=u N (x) u(x, 0) = u S (x) u(0,y)=0 u(l x,y)=0 v(x, l y ) = 0 v(x, 0) = 0 v(0,y)=v W (y) v(l x,y)=v E (y)
Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
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 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 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 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 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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More 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 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 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 informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More 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 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 information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
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 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 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 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 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 informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
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 informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
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 informationCLASSROOM NOTES PART II: SPECIAL TOPICS. APM526, Spring 2018 Last update: Apr 11
CLASSROOM NOTES PART II: SPECIAL TOPICS APM526, Spring 2018 Last update: Apr 11 1 Function Space Methods General Setting: Projection into finite dimensional subspaces t u = F (u), u(t = 0) = u I, F : B
More 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 informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
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 informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
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 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 informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
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 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 informationNumerical Schemes Applied to the Burgers and Buckley-Leverett Equations
University of Reading Numerical Schemes Applied to the Burgers and Buckley-Leverett Equations by September 4 Department of Mathematics Submitted to the Department of Mathematics, University of Reading,
More 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 informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More 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 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 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 Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
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 informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
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 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 informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationα 2 )(k x ũ(t, η) + k z W (t, η)
1 3D Poiseuille Flow Over the next two lectures we will be going over the stabilization of the 3-D Poiseuille flow. For those of you who havent had fluids, that means we have a channel of fluid that is
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 informationChapter 3 Burgers Equation
Chapter 3 Burgers Equation One of the major challenges in the field of comple systems is a thorough understanding of the phenomenon of turbulence. Direct numerical simulations (DNS) have substantially
More informationGalilean invariance and the conservative difference schemes for scalar laws
RESEARCH Open Access Galilean invariance and the conservative difference schemes for scalar laws Zheng Ran Correspondence: zran@staff.shu. edu.cn Shanghai Institute of Applied Mathematics and Mechanics,
More 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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
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 informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationA stochastic particle system for the Burgers equation.
A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u
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 informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
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 informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationNon-Oscillatory Central Schemes for a Traffic Flow Model with Arrhenius Look-Ahead Dynamics
Non-Oscillatory Central Schemes for a Traffic Flow Model with Arrhenius Look-Ahead Dynamics Alexander Kurganov and Anthony Polizzi Abstract We develop non-oscillatory central schemes for a traffic flow
More informationTMA 4195 Mathematical Modelling December 12, 2007 Solution with additional comments
Problem TMA 495 Mathematical Modelling December, 007 Solution with additional comments The air resistance, F, of a car depends on its length, L, cross sectional area, A, its speed relative to the air,
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
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 informationIntroduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués
Introduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués Contents Preliminary considerations 3. Study of the movement of a fluid................ 3.
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 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 informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
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 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 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 informationMultiscale Hydrodynamic Phenomena
M2, Fluid mechanics 2014/2015 Friday, December 5th, 2014 Multiscale Hydrodynamic Phenomena Part I. : 90 minutes, NO documents 1. Quick Questions In few words : 1.1 What is dominant balance? 1.2 What is
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 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 informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationAntony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
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 informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationHomotopy Perturbation Method for Solving Systems of Nonlinear Coupled Equations
Applied Mathematical Sciences, Vol 6, 2012, no 96, 4787-4800 Homotopy Perturbation Method for Solving Systems of Nonlinear Coupled Equations A A Hemeda Department of Mathematics, Faculty of Science Tanta
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
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 informationAdvanced Methods for Numerical Fluid Dynamics and Heat Transfer (MVKN70)
2017-08-24 Advanced Methods for Numerical Fluid Dynamics and Heat Transfer (MVKN70) 1 Lectures and course plans Week 1, 2: RY: Dr. Rixin Yu: tel: 222 3814; e-mail: rixin.yu@energy.lth.se Week 3; BS: Prof.
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 informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
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 informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
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 informationHydrodynamics. l. mahadevan harvard university
Thermodynamics / statistical mechanics : Hydrodynamics N>>1; large number of interacting particles! Hydrodynamics : L >> l; T >> t; long wavelength, slow time - average over (some) microscopic length and
More information