Finite volume method for conservation laws V
|
|
- Willis Erick Arnold
- 5 years ago
- Views:
Transcription
1 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 February 1, / 18
2 Entropy condition Let U(u) be any convex entropy of the scalar conservation law with associated entropy flux F (u), i.e., u t + x f(u) = 0 U (u)f (u) = F (u) Then the unique entropy solution of the conservation law satisfies in the sense of distributions U t + x F (u) 0 / 18
3 Definition (Entropy consistent numerical scheme) The difference scheme H is said to be consistent with the entropy condition if there exists a continuous function G : R k R which satisfies the following requirements (1) consistency with the entropy flux F () cell entropy inequality G(v,..., v) = F (v) U(v n+1 j ) U(vj n) + G(vn j k+1,..., vn j+k ) G(vn j k,..., vn j+k 1 ) 0 t x The function G is called the numerical entropy flux. If we set U n j = U(v n j ), G n j+ 1 = G(v n j k+1,..., v n j+k) then we have to verify that U n+1 j Uj n λ(g n j+ G n 1 j ) 1 3 / 18
4 When the scheme is consistent with any entropy function, we say that it is an entropy scheme. Proposition 4.1 (GR1, Chap. 3) Let H be a 3-point difference scheme with C numerical flux which is consistent with any entropy condition. Then it is at most first order accurate. Theorem Assume that the hypothesis of Lax-Wendroff theorem hold. Assume moreover that the scheme is consistent with any entropy condition. Then the limit u is the unique entropy solution of the conservation law. Theorem A monotone consistent scheme is consistent with any entropy condition. Proof: It is enough to check the entropy condition for the Kruzkov s entropy functions: for any l R U(u) = u l, F (u) = sign(u l)(f(u) f(l)) 4 / 18
5 Set a b = min(a, b), a b = max(a, b) and define the numerical entropy flux G as follows G(v k+1,..., v k ) = g(v k+1 l,..., v k l) g(v k+1 l,..., v k l) We first prove the following identity v j l λ(g j+ 1 G j 1 ) = H(v j k l,..., v j+k l) H(v j k l,..., v j+k l) By definition of G j+ 1 G j+ 1 G j 1 = [g(v j k+1 l,..., v j+k l) g(v j k+1 l,..., v j+k l)] [g(v j k l,..., v j+k 1 l) g(v j k l,..., v j+k 1 l)] = [g(v j k+1 l,..., v j+k l) g(v j k l,..., v j+k 1 l)] [g(v j k+1 l,..., v j+k l) g(v j k l,..., v j+k 1 l)] and using the finite volume scheme, we get λ[g j+ 1 j 1 ] = [H(v j k l,..., v j+k l) v j l] [H(v j k l,..., v j+k l) v j l] 5 / 18
6 which gives the identity we wanted to show since v j l v j l = v j l Now since H is monotone and consistent H(v j k l,..., v j+k l) H(v j k,..., v j+k ) H(l,..., l) = v n+1 j In the same way we obtain H(v j k l,..., v j+k l) v n+1 j Combining the above results we verify the entropy condition l v n+1 j l vj n l + λ(g n j+ G n 1 j ) 0 1 The consistency condition of G is readily satisfied since G(v,..., v) = g(v l,..., v l) g(v l,..., v l) = f(v l) f(v l) = sign(v l)(f(v) f(l)) = F (u) l 6 / 18
7 Definition (E-scheme) A consistent, conservative scheme is called an E-scheme if its numerical flux satisfies sign(v j+1 v j )(g j+ 1 f(u)) 0 for all u between v j and v j+1. Remark: Note that an E-scheme is essentially 3-point. Indeed letting v j+1 v j with first v j+1 v j and then with v j+1 v j shows that g is essentially 3-point. Remark: A 3-point monotone scheme is an E-scheme. Since g(u, v) is non-decreasing in u and non-increasing in v, we obtain and therefore g(u, v) g(u, w) g(w, w) = f(w) if u w v g(u, v) g(w, v) g(w, w) = f(w) if u w v sign(v u)(g(u, v) f(w)) 0, for all w between u and v In particular, the Godunov scheme is an E-scheme under CFl 1. 7 / 18
8 Lemma Assume that CFL 1. Then E-fluxes are characterized by { gj+ 1 gg j+ 1 g j+ 1 gg j+ 1 if v j < v j+1 if v j > v j+1 where g G stands for Godunov numerical flux. Proof: Under CFL 1, the Godunov flux is given by { g G min u [vj,v j+ = j+1] f(u) if v j < v j+1 1 max u [vj+1,v j] f(u) if v j v j+1 Assume v j < v j+1. Then E-flux satisfies g j+ 1 f(u), v j u v j+1 Assume v j > v j+1. Then E-flux satisfies v j+1 u v j g j+ 1 f(u), 8 / 18
9 Lemma Assume that CFL 1. E-schemes are characterized by 0 Q G j+ 1 Q j+ 1, j Z Proof: Assume v j < v j+1. Then for E-scheme g j+ 1 gg j+ 1 1 (f j + f j+1 ) 1 Q j+ 1 (v j+1 v j ) 1 (f j + f j+1 ) 1 QG j+ (v 1 j+1 v j ) Proposition An E-scheme with differentiable numerical flux is at most first order accurate. 9 / 18
10 Theorem (Viscous form and entropy condition) Assume that the CFL condition λ max a(u) 1 holds. An E-scheme whose coefficient of numerical viscosity satisfies Q G j+ 1 Q j+ 1 1 is consistent with any entropy condition. The proof requires two lemmas which we first prove. The basic idea is to write any E-scheme as a convex combination of the Godunov scheme and a modified Lax-Friedrichs scheme, both of which satisfy entropy condition. 10 / 18
11 Lemma (Godunov scheme) Assume that CFL 1. The Godunov scheme can be written in the following way where H G (v j 1, v j, v j+1 ) = 1 (vg + v G+ ) j 1 j+ 1 v G j 1 v G+ j+ 1 = = x x x 0 0 x w R (x/ t; v j 1, v j )dx = v j λ(f j g G j 1 ) w R (x/ t; v j, v j+1 )dx = v j + λ(f j g G j+ 1 ) Moreover if (U, F ) is any entropy pair, we have U(v G± ) U(v j± 1 j ) ± λ[f (v j ) G G j± ] 1 where G G (u, v) = F (w R (0; u, v)) 11 / 18
12 Proof: The formulae for v G± follow from the derivation of the Godunov j± 1 scheme. Now since U is a convex function, we obtain by Jensen s inequality U ( x x 0 w R (x/ t; v j 1, v j )dx ) x U(w R (x/ t; v j 1, v j ))dx x 0 Since w R (x/ t; v j 1, v j ) is by definition the entropy solution of the Riemann problem u t + x u(x, 0) = it satisfies an entropy inequality U(u) t f(u) = 0, { v j 1 x < 0 v j x > 0 + F (u) x t [0, t] 0 1 / 18
13 Integrating this last inequality on the domain (x j 1, x j) (0, t) we get In the same way we have x U(v G ) U(w j 1 R (x/ t; v j 1, v j ))dx x 0 U(v j ) λ[f (v j ) G G j ] 1 U(v G+ 0 ) U(w j+ 1 R (x/ t; v j, v j+1 ))dx x x U(v j ) + λ[f (v j ) G G j+ ] 1 Remark: This also shows that Godunov scheme satisfies the entropy condition associated with (U, F ) since U(H G (v j 1, v j, v j+1 )) 1 [U(vG ) + U(v G+ )] j 1 j+ 1 U(v j ) λ[g G j+ G G 1 j ] 1 13 / 18
14 Lemma (Lax-Friedrichs modified) Consider the 3-point scheme H M (v j 1, v j, v j+1 ) = 1 4 (v j 1 + v j + v j+1 ) 1 λ(f j+1 f j 1 ) Then under the condition CFL 1, we can write where v M j 1 v M+ j+ 1 = = H M (v j 1, v j, v j+1 ) = 1 (vm j 1 + v M+ j+ 1 ) x 1 w R (x/ t; v j 1, v j )dx = 1 x x (v j + v j 1 ) λ(f j f j 1 ) 1 x x x Moreover if (U, F ) is any entropy pair, we have where w R (x/ t; v j, v j+1 )dx = 1 (v j + v j+1 ) λ(f j+1 f j ) U(v M± j± 1 ) U(v j ) ± λ[f (v j ) G M j± 1 ] G M (u, v) = 1 1 (F (u) + F (v)) (U(v) U(u)) 4λ 14 / 18
15 Proof: Under the condition CFL 1, we obtain the desired results by integrating the conservation law over the domains (x j 1, x j ) (t n, t n+1 ) and (x j, x j+1 ) (t n, t n+1 ). We get 0 = = tn+1 xj t n x x t 0 x j 1 ( wr t + ) x f(w R) dxdt w R (x/ t; v j 1, v j )dx x (v j + v j 1 ) + [f(w R ( x/(t); v j 1, v j )) f(w R ( x/(t); v j 1, v j ))]dt But the waves from x j 1 do not reach x j 1, x j so that w R ( x/(t); v j 1, v j ) = v j, w R ( x/(t); v j 1, v j ) = v j 1 This proves the formulae for v M±. j± 1 15 / 18
16 Since w R is the entropy solution, we integrate the entropy inequality over (x j 1, x j ) (t n, t n+1 ) and using Jensen s inequality U(v M ) j 1 x 1 U(w R (x/ t; v j 1, v j ))dx x x 1 (U(v j) + U(v j 1 )) λ(f (v j ) F (v j 1 )) = U(v j ) λ(f (v j ) G M j 1 ) Remark: The modified Lax-Friedrichs scheme is consistent with any entropy condition with numerical flux G M since U(H M (v j 1, v j, v j+1 )) 1 (U(vM ) + U(v M+ )) j 1 j+ 1 U(v j ) λ[g M j+ G M 1 j ] 1 Note that the numerical viscosity of this scheme is Q M j+ 1 = 1 16 / 18
17 Proof of Theorem (Viscous form and entropy condition): Let us write the scheme in viscous form v n+1 j = vj n λ (f j+1 n fj 1) n + 1 (Qn j+ v n 1 j+ Q n 1 j v n 1 j ) 1 and in averaged form where v n+1 j = 1 (v + v + ) j 1 j+ 1 v j 1 v + j+ 1 = v j λ(f j f j 1 ) Q j 1 v j 1 = v j λ(f j+1 f j ) + Q j+ 1 v j+ 1 We can write the Godunov and modified Lax-Friedrich schemes in the same form with superscript G and M. Now since Q M = 1 and Q G Q Q M = 1 we can write, with some 0 θ j / 18
18 It follows that Q j+ 1 = θ j+ 1 QG j+ + (1 θ 1 j+ 1 )QM j+, j Z 1 v ± j± 1 = θ j± 1 vg± + (1 θ j± 1 j± 1 )vm± j± 1 If (U, F ) is any entropy pair, then due to convexity of U where U(v n+1 j ) 1 U(v j 1 ) + 1 U(v+ j+ 1 ) 1 θ j 1 U(vG ) + 1 j 1 (1 θ j 1 )U(vM ) + j 1 1 θ j+ 1 U(vG+ j+ 1 U(v j ) λ(g j+ 1 G j 1 ) ) + 1 (1 θ j+ 1 )U(vM+ ) j+ 1 G j+ 1 = θ j+ 1 GG j+ + (1 θ 1 j+ 1 )GM j+ 1 is a consistent entropy flux associated with the E-scheme under consideration. Remark: Under the conditions of the above theorem Q G Q 1, the E-scheme is also TVD and L stable. 18 / 18
Math 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 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 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 informationCONVERGENCE OF THE LAX-FRIEDRICHS SCHEME AND STABILITY FOR CONSERVATION LAWS WITH A DISCONTINUOUS SPACE-TIME DEPENDENT FLUX
Chin. Ann. Math. 25B:3(24),287 318. CONVERGENCE OF THE LAX-FRIEDRICHS SCHEME AND STABILITY FOR CONSERVATION LAWS WITH A DISCONTINUOUS SPACE-TIME DEPENDENT FLUX K. H. KARLSEN J. D. TOWERS Abstract The authors
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 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 informationLARGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMERICS FOR 1-D SCALAR CONSERVATION LAWS
LAGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMEICS FO 1-D SCALA CONSEVATION LAWS L. I. IGNAT, A. POZO, E. ZUAZUA Abstract. In this paper we analyze the large time asymptotic behavior of the discrete
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 informationON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL
ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,
More informationLOCAL NEWS. PAVBBS any at this offioe. TAKE NOT IUE. AH accounts due the. firm qf MORRIS A' III HE mutt be paid to
U Q -2 U- VU V G DDY Y 2 (87 U U UD VY D D Y G UY- D * (* ) * * D U D U q F D G** D D * * * G UX UUV ; 5 6 87 V* " * - j ; j $ Q F X * * «* F U 25 ](«* 7» * * 75! j j U8F j» ; F DVG j * * F DY U» *»q*
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 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 informationFLUX IDENTIFICATION IN CONSERVATION LAWS 2 It is quite natural to formulate this problem more or less like an optimal control problem: for any functio
CONVERGENCE RESULTS FOR THE FLUX IDENTIFICATION IN A SCALAR CONSERVATION LAW FRANCOIS JAMES y AND MAURICIO SEP ULVEDA z Abstract. Here we study an inverse problem for a quasilinear hyperbolic equation.
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 informationOn Convergence of Minmod-Type Schemes
On Convergence of Minmod-Type Schemes Sergei Konyagin, Boan Popov, and Ognian Trifonov April 16, 2004 Abstract A class of non-oscillatory numerical methods for solving nonlinear scalar conservation laws
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 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 information1. Introduction. We are interested in the scalar hyperbolic conservation law
SIAM J. NUMER. ANAL. Vol. 4, No. 5, pp. 1978 1997 c 005 Society for Industrial and Applied Mathematics ON CONVERGENCE OF MINMOD-TYPE SCHEMES SERGEI KONYAGIN, BOJAN POPOV, AND OGNIAN TRIFONOV Abstract.
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 informationNumerische Mathematik
Numer. Math. ) :545 563 DOI.7/s--443-7 Numerische Mathematik A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang Chi-Wang Shu Received: 7 July / Revised: 5 September
More informationVariational formulation of entropy solutions for nonlinear conservation laws
Variational formulation of entropy solutions for nonlinear conservation laws Eitan Tadmor 1 Center for Scientific Computation and Mathematical Modeling CSCAMM) Department of Mathematics, Institute for
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 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 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 informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
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 informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 253 Monotonicity of Implicit Finite Difference Methods for Hyperbolic Conservation Laws Michael
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 informationA posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model
A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationLois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique
Lois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique Université Aix-Marseille Travail en collaboration avec C.Bauzet, V.Castel
More informationEntropy and the numerical integration of conservation laws
Pysics Procedia Pysics Procedia 00 2011) 1 28 Entropy and te numerical integration of conservation laws Gabriella Puppo Dipartimento di Matematica, Politecnico di Torino Italy) Matteo Semplice Dipartimento
More informationGeneralized Riemann Problems: From the Scalar Equation to Multidimensional Fluid Dynamics
1 Generalized Riemann Problems: From the Scalar Equation to Multidimensional Fluid Dynamics Matania Ben-Artizi and Joseph Falcovitz Institute of Mathematics, the Hebrew University of Jerusalem, Jerusalem
More informationCONVERGENCE OF RELAXATION SCHEMES FOR HYPERBOLIC CONSERVATION LAWS WITH STIFF SOURCE TERMS
MATHEMATICS OF COMPUTATION Volume 68, Number 227, Pages 955 970 S 0025-5718(99)01089-3 Article electronically published on February 10, 1999 CONVERGENCE OF RELAXATION SCHEMES FOR HYPERBOLIC CONSERVATION
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
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 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 informationFinite Di erence Approach for the Measure Valued Vanishing Dispersion Limit of Burgers Equation
Finite Di erence Approach for the Measure Valued Vanishing Dispersion Limit of Burgers Equation Master s Thesis of Simon Laumer Supervisor: Prof. Dr. Siddhartha Mishra July 0, 014 Abstract We are concerned
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationAttribution-NonCommercial-NoDerivs 2.0 KOREA. Share copy and redistribute the material in any medium or format
Attribution-NonCommercial-NoDerivs. KOREA You are free to : Share copy and redistribute the material in any medium or format Under the follwing terms : Attribution You must give appropriate credit, provide
More informationLOWELL JOURNAL NEWS FROM HAWAII.
U # 2 3 Z7 B 2 898 X G V ; U» VU UU U z z 5U G 75 G } 2 B G G G G X G U 8 Y 22 U B B B G G G Y «j / U q 89 G 9 j YB 8 U z U - V - 2 z j _ ( UU 9 -! V G - U GB 8 j - G () V B V Y 6 G B j U U q - U z «B
More informationLocal EleveD Wins Spectacular Game
G uu «w w, u b g b bu g B vw, bu v g u,, v, ub bu gu u w gv g u w, k wu u gw u k bu g w u, u v wu w, v u wu g g bu u, u u u w, z u w u u b g wu w k bu u u!, g u w u u w u Bv Vu «J G B 2 938 FY-X Y u u
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 informationPOSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS. Cheng Wang. Jian-Guo Liu
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 3 Number May003 pp. 0 8 POSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS
More informationOVER 150 J-HOP To Q JCATS S(MJ) BY NOON TUESDAY CO-ED S LEAGUE HOLDS MEETING
F F V à y y > * y y y! F! * * F k y è 2 3 U D Y F B U Y 3 * 93 G U P B PU FF; JH H D V 50 JHP QJ (J) BY UDY P G y D; P U F P 0000 GUU VD PU F B U q F' yyy " D Py D x " By ; B x; x " P 93 ; Py y y F H j
More informationMAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS
MATHEMATICS OF COMPUTATION S 005-5780054-7 Article electronically published on May 6, 0 MAXIMUM PRINCIPLE AND CONVERGENCE OF CENTRAL SCHEMES BASED ON SLOPE LIMITERS ORHAN MEHMETOGLU AND BOJAN POPOV Abstract.
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 informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
More informationA positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations
A positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations Sashank Srinivasan a, Jonathan Poggie a, Xiangxiong Zhang b, a School of Aeronautics and Astronautics,
More informationTHE ENTROlm~ DISSIPATION BY NUMERICAL VISCOSITY IN NONLINeAR CONSERVATIVE DIFFERENCE SC~I~ES Eitan Tadmor* School of Mathematical Sciences, Tel-Aviv University and Institute for Computer Applications in
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 informationFractional Roman Domination
Chapter 6 Fractional Roman Domination It is important to discuss minimality of Roman domination functions before we get into the details of fractional version of Roman domination. Minimality of domination
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 informationu,+f(u)x=o, A GEOMETRIC APPROACH TO HIGH RESOLUTION TVD SCHEMES*
SIAM J. NUMER. ANAL. Vol. 25, No. 2, April 1988 (C) 1988 Society for Industrial and Applied Mathematics OO2 A GEOMETRIC APPROACH TO HIGH RESOLUTION TVD SCHEMES* JONATHAN B. GOODMAN? AND RANDALL J. LeVEQUE$
More information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
More informationGrand Bapids Business Directory.
G 2 by U- b 5 y VU W HG WY UGU 2787 U 8 W UH VY Y G W Y- H W H G F U: yy 7 x V by 75 F VQ ) > q Fy y bb - W? G W W : 8 $8 6 2 8 8 36 5 35 8 5 7 bq U / Yy U qu b b b } GU H 7 W b U 8HUH Q F U y y? 8 x Q
More informationH r# W im FR ID A Y, :Q q ro B E R 1 7,.,1 0 1 S. NEPTUNE TH RHE. Chancelor-Sherlll Act. on Ballot at ^yisii/
( # Y Q q 7 G Y G K G Q ( ) _ ( ) x z \ G 9 [ 895 G $ K K G x U Y 6 / q Y x 7 K 3 G? K x z x Y Y G U UY ( x G 26 $ q QUY K X Y 92 G& j x x ]( ] q ] x 2 ] 22 (? ] Y Y $ G x j 88 89 $5 Y ] x U $ 852 $ ((
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 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 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 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 informationConvergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions
Convergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions Marie Hélène Vignal UMPA, E.N.S. Lyon 46 Allée d Italie 69364 Lyon, Cedex 07, France abstract. We are interested
More informationEntropy stable high order discontinuous Galerkin methods. for hyperbolic conservation laws
Entropy stable high order discontinuous Galerkin methods for hyperbolic conservation laws Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Tianheng Chen, and with Yong Liu
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 informationCONVERGENCE RATE OF APPROXIMATE SOLUTIONS TO CONSERVATION LAWS WITH INITIAL RAREFACTIONS
CONVERGENCE RATE OF APPROXIMATE SOLUTIONS TO CONSERVATION LAWS WITH INITIAL RAREFACTIONS HAIM NESSYAHU AND TAMIR TASSA Abstract. We address the question of local convergence rate of conservative Lip +
More 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 informationRearrangement on Conditionally Convergent Integrals in Analogy to Series
Electronic Proceedings of Undergraduate Mathematics Day, Vol. 3 (8), No. 6 Rearrangement on Conditionally Convergent Integrals in Analogy to Series Edward J Timko University of Dayton Dayton OH 45469-36
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 informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More informationConsistency analysis of a 1D Finite Volume scheme for barotropic Euler models
Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models F Berthelin 1,3, T Goudon 1,3, and S Mineaud,3 1 Inria, Sophia Antipolis Méditerranée Research Centre, Proect COFFEE Inria,
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 informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationModeling and Numerical Approximation of Traffic Flow Problems
Modeling and Numerical Approximation of Traffic Flow Problems Prof. Dr. Ansgar Jüngel Universität Mainz Lecture Notes (preliminary version) Winter Contents Introduction Mathematical theory for scalar conservation
More informationRunge Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems
Journal of Scientific Computing, Vol. 16, No. 3, September 2001 ( 2001) Review Article Runge Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems Bernardo Cockburn 1 and Chi-Wang Shu
More informationCS205B / CME306 Homework 3. R n+1 = R n + tω R n. (b) Show that the updated rotation matrix computed from this update is not orthogonal.
CS205B / CME306 Homework 3 Rotation Matrices 1. The ODE that describes rigid body evolution is given by R = ω R. (a) Write down the forward Euler update for this equation. R n+1 = R n + tω R n (b) Show
More informationAbstract. A front tracking method is used to construct weak solutions to
A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar
More informationMATH 217A HOMEWORK. P (A i A j ). First, the basis case. We make a union disjoint as follows: P (A B) = P (A) + P (A c B)
MATH 217A HOMEWOK EIN PEASE 1. (Chap. 1, Problem 2. (a Let (, Σ, P be a probability space and {A i, 1 i n} Σ, n 2. Prove that P A i n P (A i P (A i A j + P (A i A j A k... + ( 1 n 1 P A i n P (A i P (A
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 informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationCoupling algorithms for hyperbolic systems
Coupling algorithms for hyperbolic systems Edwige Godlewski Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie (Paris) Cemracs, July 2011 Collaborations with Frédéric Coquel ( CMAP, Cnrs),
More informationNumerical methods for conservation laws with a stochastically driven flux
Numerical methods for conservation laws with a stochastically driven flux Håkon Hoel, Kenneth Karlsen, Nils Henrik Risebro, Erlend Briseid Storrøsten Department of Mathematics, University of Oslo, Norway
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationNumerische Mathematik
Numer. Math. (001) 88: 513 541 Digital Object Identifier (DOI) 10.1007/s0011000033 Numerische Mathematik Convergence rates to discrete shocks for nonconvex conservation laws Hailiang Liu 1,, Jinghua Wang,,
More informationM545 Homework 8. Mark Blumstein. December 9, [f (x)] 2 dx = 1. Find such a function if you can. If it cannot be found, explain why not.
M545 Homework 8 Mark Blumstein ecember 9, 5..) Let fx) be a function such that f) f3), 3 [fx)] dx, and 3 [f x)] dx. Find such a function if you can. If it cannot be found, explain why not. Proof. We will
More informationThe intersection and the union of the asynchronous systems
BULETINUL ACADEMIEI DE ŞTIINŢE A REPUBLICII MOLDOVA. MATEMATICA Number 3(55), 2007, Pages 3 22 ISSN 1024 7696 The intersection and the union of the asynchronous systems Serban E. Vlad Abstract. The asynchronous
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 informationSkr. K. Nor. Vidensk. Selsk., 2003, no. 3, 49 pp
STABILITY FOR ENTROPY SOLUTIONS OF NONLINEAR DEGENERATE PARABOLIC CONVECTION-DIFFUSION EQUATIONS WITH DISCONTINUOUS COEFFICIENTS L 1 K. H. KARLSEN, N. H. RISEBRO, AND J. D. TOWERS Skr. K. Nor. Vidensk.
More informationNotes on Computational Mathematics
Notes on Computational Mathematics James Brunner August 19, 2013 This document is a summary of the notes from UW-Madison s Math 714, fall 2012 (taught by Prof. Shi Jin) and Math 715, spring 2013 (taught
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
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 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 informationSpectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem
Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v
More informationEvaluation of The Necessary of Agriculture Public Expenditure for Poverty Reduction and Food Security in Benin
MPRA Mu P RPE Av Evu N Aguu Pu Exu Pv Ru F Su B A C L Zg Uv E Lw O 010 O ://.u.u-u./447/ MPRA P N. 447, 7. Augu 010 0:17 UC qwugjkzxvqw ugjkzxvqwu gjkzxvqwug [ Evu N jkzxvqwugjkzx vqwugjkzxv qwugjkzxvqw
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 informationExistence result for the coupling problem of two scalar conservation laws with Riemann initial data
Existence result for the coupling problem of two scalar conservation laws with Riemann initial data Benjamin Boutin, Christophe Chalons, Pierre-Arnaud Raviart To cite this version: Benjamin Boutin, Christophe
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 information