Finite Volume Schemes: an introduction
|
|
- Cornelius Eaton
- 6 years ago
- Views:
Transcription
1 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 in Scienze dell Ingegneria Civile e Ambientale Corso di Metodi Numerici marzo 2008
2 Outline Hyperbolic Conservation Laws Classical and Weak solutions The Riemann problem Finite Volume Methods in 1D Finite Volume Methods in 2D The lectures about Finite Volume and Mixed Hybrid Finite Element will be on line at A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
3 Hyperbolic Conservation Laws Finite Volumes (FV) schemes are most useful for modelling hyperbolic conservation laws, which, in their simplest (scalar) form may be expressed by the PDE u t + f = 0 on a domain Ω, with u(x, t) specified on inflow boundaries. Example: the equation for conservation of mass in groundwater contamination problem. Assumptions: no diffusion, no mechanical dispersion, density of the contaminant and the water constant, velocity of the fluid mixture constant, Let ρ(x, t) be the density of the contaminant at point x and time t and v(x, t) be the velocity of the fluid mixture. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
4 The total pollutant mass in two sections 1 and 2 is given by x 2 x 1 ρ(x, t) dx. If no source or sinks are present, the mass can change only because of the mixture flowing across the endsections in x 1 or x 2. The rate of mass change is given by: d dt x2 x 1 ρ(x, t) dx = ρ(x 1, t)v(x 1, t) ρ(x 2, t)v(x 2, t) By integrating in time from t k to t k+1 we obtain the integral form of the conservation law. Assuming ρ and v differentiable, with simple calculations we obtain the differential form of the conservation law: ρ t + (ρv) x = 0 A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
5 Classical and Weak solutions Solution of hyberbolic conservation laws may be obtained if initial and boundary conditions are given. The simplest problem is the initial value problem, or Cauchy problem, defined for < x < and t 0. We must specify initial conditions only: u t + f = 0 u(x, 0) = u 0 (x) < x <. Classical solutions of this problem are constant along the characteristics curves defined by x(t) such that dx dt = f (u(x(t), t)) t 0 x(0) = x 0 A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
6 Shock waves u is constant along these characteristics and the characteristics travel at constant velocity which is equal to f (u 0 (x 0 )). Generally, classical solutions exist only for t [0, T ] where T 1 = inf x u 0 (x)f (u 0 (x)). At the time t = T the characteristics first cross, the function u(x, t) has an infinite slope the wave is said to break by analogy with waves on a beach and a shock forms. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
7 Weak solutions To allow discontinuities, which arise in a natural way in this situation, we define a weak solution of conservation law. Definition (Weak solution) A function u(x, t), bounded and measurable, is called a weak solution of the conservation law, if for each φ C (R R + ), the following equality holds: + 0 [φ t u + φ x f (u)] dx dt = + φ(x, 0)u(x, 0) dx. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
8 The Riemann problem Definition (The Riemann problem) A Riemann problem is the conservation law together with initial data consisting of two constant states separated by a single discontinuity, { ul x < 0, u 0 (x) = u r x > 0. Example: The Riemann problem in Burgers equation u t + ( 1 2 u2 ) x = 0 A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
9 If u l > u r, there is a unique weak solution, { ul x < st u(x, t) = x > st. s = (u l + u r ) is the shock speed 2 If u l < u r, there are infinitely many weak solutions, since between the points u l t < x < u r t, there is no information available from the characteristics. To determine the correct physical behavior we adopt the vanishing viscosity approach by adding the term ɛu xx depending on small parameter ɛ u r u t + ( 1 2 u2 ) x = ɛu xx. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
10 If the initial data is smooth and ɛ very small then, before the wave begins to break, the ɛu xx term is negligible compared to other terms and the solutions to the two PDEs look nearly identical. As the wave begins to break, the term u xx grows much faster than u x and at some point the ɛu xx term is comparable to the other terms and begins to play a role. This term keeps the solution smooth for all time, preventing the breakdown of solutions that occurs for the hyperbolic problem. The weak solution is called rarefaction wave or expansion fan: u l x < u l t u(x, t) = x/t u l t x u r t x > u r t u r A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
11 u l u r Shock wave Rarefaction wave u r u l 0 0 A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
12 Some common numerical schemes To solve numerically hyberbolic conservation laws some common numerical schemes are: Finite Differences Finite Volumes Finite Elements Discontinuous Galerkin We will see more in detail (but not too much...) the Finite Volume approach. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
13 Properties of the ideal scheme The ideal scheme for the resolution of hyperbolic conservation laws would reflect the properties of the underlying physical and mathematical system: Conservative: for correct capturing of discontinuities; Accurate: obvious...; Free of spurious oscillations: spurious oscillations have no physical meaning. This goal is remarkably difficult to achieve; Continuous: for convergence to the steady state solution. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
14 The Finite Volume method In the Finite Volume method the three main steps to follow are: Partition the computational domain into control volumes (or control cells) - wich are not necessarily the cells of the mesh. Discretize the integral formulation of the conservation laws over each control volume (by applying the divergence theorem). Solve the resulting set of algebraic equations or update the values of the dependent variables. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
15 FV in 1D We will see two FV schemes in 1D: Godunov s method and van Leer s method. Godunov s method uses the exact solution of the local Riemann problem and does not produce oscillations around discontinuities, but it is only first order accurate and the solutions display large numerical viscosity. van Leer s method can be viewed as a generalization of Godunov s method, it is a high resolution method, characterized by second order accuracy on smooth solutions and the absence of spurious oscillations. Both methods are conservative. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
16 Conservative method Definition Given a uniform grid with time step t and spatial mesh size x, a numerical method is said to be conservative if the corresponding scheme can be written as: v n+1 j = vj n λ(g n g n ), j Z n 0 j+ 1 j where vj n approximates u(x j, t n ) at the point (x j = j x, t n = n t), λ = t x and g : R2k R is a continuous function, called the numerical flux (function), that defines a (2k + 1)-point scheme. g n j+ 1 2 = g(v n j k+1,..., v n j+k). A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
17 Motivation of the conservative scheme v n j can be view as an approximation of u n j = 1 x xj+1/2 x j 1/2 u(x, t n ) dx. u n j is the average of u(, t n ) on the cell [x j 1/2, x j+1/2 ] (where x j±1/2 = x j ± x 2 ). From the integral form of the conservation law, we have: xj+1/2 x j 1/2 xj+1/2 u(x, t n+1 ) dx = u(x, t n ) dx [ t n+1 t n x j 1/2 f (u(x j+1/2, t)) dt t n+1 t n f (u(x j 1/2, t)) dt]. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
18 Dividing by x and using u n j we get u n+1 j = u n j 1 t n+1 x [ f (u(x j+1/2, t)) dt t n t n+1 t n f (u(x j 1/2, t)) dt]. The numerical flux function can be considered as an average flux through x j+1/2 over the time interval [t n, t n+1 ], g n j+1/2 = 1 t t n+1 t n f (u(x j+1/2, t)) dt. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
19 Godunov s scheme 1. Projection step Given data vj n at time t n, construct a piecewise constant function ˆv j n (x, t n ) defined by ˆv n j (x, t n ) = v n j x j 1/2 x x j+1/2. v_(j 1) v_(j+1) v_j x_(j 1) x_j x_(j+1) A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
20 2. Evolution step Solve the local Riemann problem at the cell interfaces, that is, on each subinterval [x j, x j+1 ] and for t t n, solve ˆv n j t + n f (ˆv j ) x{ ˆv n j (x, t n ) = = 0 vj n, x j < x < x j+1/2 vj+1 n x j+1/2 < x < x j+1 3. Projection step Define the approximation v n+1 j at time t n+1 by averaging the Riemann problem solution ˆv j n at the time t n+1, so that v n+1 j = 1 x xj+1/2 x j 1/2 ˆv n j (x, t n+1 ) dx. These values are then used to define new piecewise constant data ˆv n+1 j (x, t n+1 ) and the process repeats. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
21 Example In the simple case of a linear advection equation u t + au x = 0, with a > 0, the third step gives v n+1 j = v n j a t x (v j n vj 1 n ). u_(j 1) u_j < > a t u_(j 1) < > a t x_(j 1) x_j x_(j+1) A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
22 Stability condition The CFL (Courant-Friedrichs-Lewy) condition λ max f (u) 1 u 2, where λ = t is a necessary stability condition stating that x the domain of dependence of the method includes the domain of dependence of the PDE. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
23 van Leer s method First and third steps of Godunov s methods are of a numerical nature and can be modified without influencing the physics, in order to define spatially high order accurate scheme. van Leer s scheme is a spatially second order accurate scheme (known also as MUSCL (Monotone Upstream-centered Scheme for Conservation Laws) method). The straightforward replacement of the first-order scheme by a second-order accurate interpolation leads to the generation of oscillations around discontinuities. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
24 To overcome this limitations, non linear components are introduced. This concept was introduced by van Leer under the form of limiters, i.e. functions that control the gradient of the computed solution to prevent the appearance of unphysical over/under shoots. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
25 1 Reconstruction step: the dependent variable is interpolated using a piecewise linear function ˆv starting from the cell averages vj n, vj 1 n, v j+1 n. To this end define Sj n to be the slope on the jth cell calculated using vj 1 n and v j+1 n. Then v ˆv n j n + (x x j ) Sn j x j 1/2 < x < x j+1/2, x (x) = ˆv n (x j 1/2 ) x x j 1/2 ˆv n (x j+1/2 ) x x j+1/2. Note that taking Sj n = 0 for all j and n recovers Godunov s method; A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
26 2 Evolution step: the waves are propagated across cell interfaces according to an exact or approximate solution of a local Riemann problem that uses the interpolated values ˆv n (x) as initial conditions. One solves t w + x f (w) = 0 x R, t n t t n+1 w(x, t n ) = ˆv n (x). This step yields w(, t n+1 ). 3 Cell-averaging step: v n+1 j is obtained by projecting the solution w(x, t n+1 ) onto the piecewise constant functions v n+1 j = 1 xj+1/2 w(x, t n+1 ) dx. x x j 1/2 A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
27 a) b) x j-1 x j x j+1 x x x j-1 j j+1 v t c) d) x x x j-1 j j+1 x x x j-1 j j+1 a) reconstruction; b) limiting; c) evolution; d) cell-averaging step. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
28 Slope limiter The numerical scheme is completely defined after specification of the slope limiter S n j. Poor choice of slopes may give oscillations in the solution. Generation of oscillations can be prevented by acting on their production mechanism and introducing nonlinear correction factors, so called the limiters, that force the method to be total variation diminishing (TVD). Definition A numerical method to solve hyperbolic conservation laws is called total variation diminishing if + j= v n+1 j+1 v n+1 j + j= v n j+1 v n j A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
29 To obtain a slope limiter, let S n j = Ŝn j Φ n j where Ŝn j represents the actual slope approximation and Φ j = Φ(θj n ) is a limiter function, defined in such a way that the method is TVD and θj n is the ratio of two consecutive gradients, i.e.: θj n = v j n vj 1 n vj+1 n v. j n A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
30 Limiter functions Various limiter functions can be found in literature. For example: van Leer s limiter: Φ(θ) = θ + θ 1 + θ, minmod limiter that represents the lowest boundary of the second-order TVD region: { min(θ, 1) if θ > 0 Φ(θ) = 0 if θ 0. the superbee limiter, that represents the upper limit of the second-order TVD region and has been introduced by Roe: Φ(θ) = max[0, min(2θ, 1), min(θ, 2)]. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
31 Numerical example with periodic boundary conditions u t + u x = 0, with initial condition u 0 (x) periodic of period 2 defined on the interval [ 1, 1] as 1 1 x < x < 0.25 u 0 (x) = x < x < x 1 We solve this problem at time t = 2 s with λ = 0.5, by setting t = s and x = m. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
32 1.0 (a) Godunov scheme 1.0 (b) van Leer scheme van Leer limiter (c) van Leer s scheme minmod limiter 1.0 (d) van Leer s scheme superbee limiter A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
33 Numerical example with non periodic boundary conditions u t + u x = 0, with the following initial condition: { sin(πx) 0 x 2 u(x, 0) = u 0 (x) = 0 otherwise. We solve this problem at time t = 2 s with λ = 0.5 and x = m. A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
34 1.0 (a) Godunov s scheme 1.0 (b) van Leer s scheme van Leer limiter (c) van Leer s scheme minmod limiter 1.0 (d) van Leer s scheme superbee limiter A. Mazzia (DMMMSA) Finite Volume a.a. 2007/ / 34
Mixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
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 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 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 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 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 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 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 informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationThe 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 informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
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 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 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 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 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 informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
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 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 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 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 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 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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
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 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 informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
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 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 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 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 informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationA Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids
Diese Arbeit wurde vorgelegt am Lehrstuhl für Mathematik (MathCCES) A Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids CES Seminar Computational Engineering
More informationAdvection algorithms II. Flux conservation, subgrid models and flux limiters
Chapter 4 Advection algorithms II. Flu conservation, subgrid models and flu limiters In this chapter we will focus on the flu-conserving formalism of advection algorithms, and we shall discuss techniues
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 informationState of the Art MHD Methods for Astrophysical Applications p.1/32
State of the Art MHD Methods for Astrophysical Applications Scott C. Noble February 25, 2004 CTA, Physics Dept., UIUC State of the Art MHD Methods for Astrophysical Applications p.1/32 Plan of Attack Is
More informationSpatial Discretization
Spatial Discretization Shengtai Li and Hui Li Los Alamos National Laboratory 1 Overview Finite Difference Method Finite Difference for Linear Advection Equation 3 Conservation Laws Modern Finite Difference
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
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 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 informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
More 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 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 informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
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 informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More 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 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 (Some) Algorithms for Edge Gyrokinetics
Comparison of (Some) Algorithms for Edge Gyrokinetics Greg (G.W.) Hammett & Luc (J. L.) Peterson (PPPL) Gyrokinetic Turbulence Workshop, Wolfgang Pauli Institute, 15-19 Sep. 2008 w3.pppl.gov/~hammett Acknowledgments:
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 informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More 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 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 informationRecovery of high order accuracy in radial basis function approximation for discontinuous problems
Recovery of high order accuracy in radial basis function approximation for discontinuous problems Chris L. Bresten, Sigal Gottlieb 1, Daniel Higgs, Jae-Hun Jung* 2 Abstract Radial basis function(rbf) methods
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 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 informationThe Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics
The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
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 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 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 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 informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
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 informationLarge Time Step Scheme Behaviour with Different Entropy Fix
Proceedings of the Paistan Academy of Sciences: A. Physical and Computational Sciences 53 (): 13 4 (016) Copyright Paistan Academy of Sciences ISSN: 0377-969 (print), 306-1448 (online) Paistan Academy
More 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 informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
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 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 informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
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 informationAdvanced numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationTotal Variation Theory and Its Applications
Total Variation Theory and Its Applications 2nd UCC Annual Research Conference, Kingston, Jamaica Peter Ndajah University of the Commonwealth Caribbean, Kingston, Jamaica September 27, 2018 Peter Ndajah
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationSecond Order Positive Schemes by means of Flux Limiters for the Advection Equation
Second Order Positive Schemes by means of Flux Limiters for the Advection Equation Riccardo Fazio and Alessandra Jannelli Abstract In this paper, we study first and second order positive numerical methods
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 informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
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 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 informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More information