E = where γ > 1 is a constant spesific to the gas. For air, γ 1.4. Solving for p, we get. 2 ρv2 + (γ 1)E t
|
|
- Holly Shelton
- 5 years ago
- Views:
Transcription
1 . The Euler equations The Euler equations are often used as a simplification of the Navier-Stokes equations as a model of the flow of a gas. In one space dimension these represent the conservation of mass, momentum and energy, and read ) v E t + v v + p ve + p) Here denotes the density of the gas, v the velocity, p the pressure and E the energy. To close this system, i.e., to reduce the number of unknowns to the number of equations, one can add a constitutive law relating these. Such laws are often called equations of state and are deduced from thermodynamics. For a so called ideal polytropic gas the equation of state takes the form E = p γ + v, where γ > is a constant spesific to the gas. For air, γ.4. Solving for p, we get ) p = γ ) E γ =. v = γ ) E γ q, where the momentum q = v. Inserting this in the Euler equations v v + γ 3 v + γ )E E v =. γe γ t v) In the conserved variables, q and E, this system of conservation laws read q 3) q 3 γ ) q + + γ ) E =. ) E q 3 Set u = Then the Jacobian dfu) reads t γ Eq γ q and fu) = E dfu) = Introducing the enthalpy as the Jacobian can be rewritten as dfu) = γ Eq γ γ 3 3 γ ) q γ Eq q + γ ) E. ) γ q 3 ) q 3 γ) q γ. + γ ) q3 γ E 3 3γ ) q γ q H = E + p = γ E ) γ q = γ ) p + γ v, ) v 3 γ) v γ ). v 3 vh H γ ) v γv γ 3 To find its eigenvalues, we compute the determinant det λi dfu)) = λ [ λ 3 γ)v) λ γv) + γ ) γ ) v H )]
2 + 3 γ v λ γv) + γ ) vh γ ) v 3 = λ [ λ 3vλ + γ3 γ)v + γ ) v + γ )H ] + 3 γ v λ γ + )v3 + γ )Hv [ = λ λ 3vλ + v + ] γ + ) v γ )H γ + )v3 + γ )vh [ = λ λ v)λ v) + ] γ + )v γ )H γ + )v3 + γ )vh [ = λ v) λλ v) + ] γ + )v γ )H [ = λ v) λ v) v )] γ + )v + γ )H [ )] γ = λ v) λ v) H v ). This can be simplified further by introducing the sound speed c, by c = γp. We then calculate H v = γ E γ )v v Therefore Thus the eigenvalues of the Jacobian are = γ E γv ) E = γ v = γ E v ) = γ p γ. detλi dfu)) = λ v) [ λ v) c ]. 4) λ u) = v c, λ u) = v, λ 3 u) = v + c. As for the corresponding eigenvectors, we write these as r i =, y i, z i ), and we see that y i = λ i, and z i = λ i ) γ γ 3)v + λ i γ 3)v. For i = we find that z = v ) γ γ 3)v + vγ 3)v + c cv γ 3)cv ) γ Recall that this is in, q, E) coordinates. = v + c γ cv
3 For i = 3 we similarly calculate ) γp = v + cv γ ) = H cv. z 3 = H + cv, and for i = it is straightforward to see that z = v /. Summing up we have the following eigenvalues and eigenvectors λ u) = v c, r u) = v c, H cv 5) λ u) = v, r u) = v, v λ 3 u) = v + c, r 3 u) = v + c. H + cv At this point it is convenient to introduce the concept of an i Riemann invariant. An i Riemann invariant is a function R = R, q, E) such that R is constant along the integral curves of r i. In other words, an i Riemann invariant satisfies Ru) r i ) =. The usefulness of this is that if we can find for each of the three eigenvectors, two Riemann invariants Ru) and Ru), then we can possibly solve the equations R, q, E) = R l, q l, E l ), R, q, E) = Rl, q l, E l ) to obtain a formula for the rarefaction waves. This is the equivalent to finding an implicit solution of the ordinary differential equation, u = ru), defining the rarefaction curves. It turns out that we have the following Riemann invariants { S, i =, Riemann invariants: v + c γ {, v, i =, Riemann invariants: p, { S, i = 3, Riemann invariants: v c γ, where we have introduced the entropy S S = log ) p γ. Verify these statements! Now we can try to obtain solution formulas for the rarefaction curves, for i =, this curve is given by ) γ p = p l, v = v l + c ) ) γ )/ l. l γ This curve is parameterized by. We must check which half of the curve to use, this will be the part where λ = v c is increasing. On the curve we have v) c) = v l + c ) ) γ )/ ) / l γp) γ l l 3
4 4 = v l + c l γ = v l + c l γ = v l + c l γ ) ) γ )/ ) / γpl l ) ) γ )/ ) γ )/ c l l l γ + l ) ) γ )/. l l ) γ )/ Since γ >, we see that v) c) is decreasing in, and for the -rarefaction wave we must use < l. Since p) is increasing in, this also means that we use the part where p < p l. Therefore we can use p as a parameter in the curve for v and write the -rarefaction curve as v p) = v l + c l γ ) ) γ )/γ) p, p p l. General theory tells us that at least for p close to p l ) this curve can be continued smoothly as a -shock curve. To find the rarefaction curve of the third family, we adopt the viewpoint that u r is fied, and we wish to find u as a function of u r. In the same way as for v this leads to the formula v 3 p) = v r + c r γ p l ) ) γ )/γ) p, p p r. To find how the density varies along the rarefaction curves we can use that the entropy s is constant, leading to ) /γ p =. l p l Now we turn to the computation of the Hugoniot loci. We view the left state u l as fied, and try to find the right state u, recall the notation u = u u l. The Rankine Hugoniot relations for ) are 6) p r s = v s v = v + p s E = ve + p), where s no longer denotes the entropy, but the speed of the discontinuity. Now we introduce new variables by w = v s and m = w. Then the first equation in 6) reads s s l = w + s l w l s l, which implies that m =. Similarly the second equation reads sw + s sw l s l = w + s) l w l + s) + p s m + s = w + w + s l w l l w l s l + p s = w + p + s. Hence mw + p =. Finally the third equation in 6) reads which implies se se l = Ew + Es + pw + ps E l w l E l s p l w l p l s E = E ) l m + pw p l w l + s p l E = E l + p l p ) l m sm w l
5 E + p = m sw c = m γ + w + s) sw c = m γ + w Hence the Rankine Hugoniot conditions are equivalent to 7) m = mw + p = c m γ + w =. We immediately find one solution by setting m =, which implies p = and v =. This is the contact discontinuity. Hence we assume that m to find the other Hugoniot loci. Now we introduce auiliary parameters Using these we have that 8) Then the third equation in 7) reads which can be rearranged as so that 9) π = p p l, and z = l. c c l = π z c l γ + w l = c l and w w l = z. c l π γ z + w l z π ) ) = wl γ z z, wl c l ) = zz π) γ z. Net recall that p = c /γ. Using this, the second equation in 7) reads or which again can be rearranged as z c γ + w = lc l γ ) c z γ + w = c l γ + w l + lw l ) c l π γz + w l z = c l γ + w l. Dividing by c l, we can solve for w l/c l ), viz. ) wl ) = zπ ) γ z. Equating ) and 9) and solving for z yields ) z = βπ + π + β, c l 5
6 6 where Using this epression for z in 9) we get wl c l β = γ + γ. ) = πβ+ π+β γ πβ+ π+β π ) πβ+) π+β) = πβ + ) π ) γ π ) β ) = πβ + γ β. Note that γ > implies β > so that this is always well defined. Since w l = v l s, we can use this to get an epression for the shock speed, βπ + ) s = v l c l γ β, where we use for the first family and + for the third. Net, using 8), v s v l s = z which can be used to epress v as a function of π. βπ + v = v l c l γ β ) ± π + β πβ + c πβ + l γ β ) ) β )π ) = v l c l γ β πβ + ) ) πβ + = v l c l γγ ) π πβ + ) /, where we take the sign for the first family and the + sign for the third. To see how the density varies along the Hugoniot loci, we use that = l z, or 3) = l πβ + π + β, which holds for both the first and the third family. Recall that for the first family we obtained a rarefaction solution if p p l, hence we have a shock solution if p p l. This means that the whole solution curve for waves of the first family is given by ) ) γ )/γ) γ p p l, p p l, 4) v p) = v l + c l ) ) / p γγ ) p l + β p p l, p pl. To find the density along this solution curve we have the formula ) /γ p p l, p pl, 5) p) = l +β p p l β+, p p p l. p l Since the second family is linearly degenerate, we can use the whole integral curve of r. Using the Riemann invariants, this is given simply as 6) v = v l, p = p l.
7 7 For the third family, we take the same point of view as for the shallow water equations; we keep u r fied and look for states u such that the Riemann problem { u <, u, ) = u r >, is solved by a wave shock or rarefaction) of the third family. By much the same calculations as for the first family we end up with ) ) γ )/γ) γ p p r, p p r rarefaction part), 7) v 3 p) = v r c r ) ) / p γγ ) p r + β p p r, p pr shock part). Regarding the density along this curve, it will change according to ) /γ p p r, p pr, 8) 3 p) = r +β p ), p p r. pr β+ p pr Now for any l, the curve v p) is a strictly decreasing function of p for p in the range [, ) taking values in the set, v l + c l /γ )]. Similarly, for any r v p) is a strictly increasing function of p, taking values in the set [v r c r /γ ), ). It follows that these curves will intersect in one point p m, v m ) if v r c r γ v l + c l γ. If this does not hold, then v does not intersect v 3 and we have a solution with vacuum. However if they intersect, then we can find the solution of the Riemann problem by finding p m, v m ).... Sod s shock tube problem. We consider an initial value problem similar to the dam break problem for shallow water. The initial velocity is everywhere zero, but the pressure to the left is higher than the pressure on the right. Specifically we set { <, 9) p, ) = v, ) =,, ) =., We have used γ =.4. In Figure we show the solution to this Riemann problem in the p, v) plane and in the, t) plane. We see that the solution consists of a leftward moving rarefaction wave of the first family, followed by a contact discontinuity and a shock wave of the third family. In Figure we show the 8 Solution in the p, v)plane Solution in the, t) plane v p).8.7 v 4 3 p m,v m ) v p) t p r,v r ) p l,v l ) p 3 Figure. The solution of the Riemann problem 9).
8 8 pressure, velocity, density and the Mach number as functions of /t. The Mach number is defined to be v /c, so that if this is larger than, the flow is called supersonic. The solution found here is actually supersonic between the contact discontinuity and the shock wave. Pressure.4 Velocity. p 8 6 v Density Mach number Figure. Pressure, velocity, density and the Mach number for the solution of 9)... The Euler equations and entropy. We shall show that the physical entropy is in fact also a mathematical entropy for the Euler equations, in the sense that ) S) t + vs), weakly for any weak solution u =, q, E) which is the limit of solutions to the viscous approimation. To this end, it is convenient to introduce the internal specific energy, defined by e = E ) v. Then the Euler equations read ) e + )) v + t t + v) = v) t + v + p ) ) = v + ev + pv =.
9 9 For classical solutions, this is equivalent to the nonconservative form t + v + v = ) v t + vv + p = Verify this! We have that e t + ve + p v =. ) p S = log γ ) γ )e = log γ = γ ) log) loge) logγ ). Thus we see that S = γ > and S e = e <. These inequalities are general, and thermodynamical mumbo-jumbo implies that they hold for any equation of state, not only for polytropic gases. For classical solutions we can compute Therefore S t = S t + S e e t = γ v + v ) + ve + p ) e v = γ ) p ) v γ ) e e ) v e = vs. S t + vs = for smooth solutions to the Euler equations. This states that the entropy of a particle of the gas remains constant as the particle is transported with the velocity v. Furthermore, S) t = t S + S t = v) S vs = vs). Thus for smooth solutions the specific entropy ηu) = Su) is conserved 3) S) t + vs) =. Of course, combining this with ) and viewing the entropy as an independent unknown, we have four equations for three unknowns, so we cannot automatically epect to have a solution. Sometimes one considers models where the energy is not conserved but the entropy is, so called isentropic flow. In models of isentropic flow the third equation in ) is replaced by the conservation of entropy 3). To show that ) holds for viscous limits, we must show that the map u ηu) = u)su), eu)) is conve. We have that η is conve if its Hessian d η is a positive definite matri. For the moment we use the convention that all vectors are column vectors, and for a vector a, a T denotes its transpose. Then we compute d η = d u)su)) = S) T + S ) T + d S.
10 Now we view Su) as Su), eu)) and compute d Su), eu)) = S ) T + S e e ) T + e) T ) + S ee e e) T + S e d e. If we use this in the previous equation, we end up with where C is given by d ηu) = S + S ) ) T + S e e) T + e ) T ) + S ee e e) T S e C, ) C = d e + e) T + e ) T. Trivially we have that =,, ) T. Furthermore so we have E q + q eu) = E, ) T = e + T v, v, ). e = 3, q, The Hessian of e is given by E d 3 q 3 q 4 3 e = q 3 = e v v v. Net, Then e) T + e ) T = e + v v ) + e + v v ) = e + v v v. Now introduce the matri D by C = e v v v e + v v v = v v v. D = v v + e v. We have that D is invertible, and thus d η is positive definite if and only if Dd ηd T is positive definite. Now Dd ηu)d T = S + S ) D D ) T + S e D D e) T + D e D ) T ) + S ee D e D e) T S e DCD T. We compute D =,
11 D e =, DCD T =, and using this Dd ηu)d T = S + S ) + S e + S ee S e S + S S e = S e S e S ee = e. e γ Hence A has three positive eigenvalues, and is positive definite, therefore also d η is positive definite, and η is conve..5 Entropy 5 Specific entropy S η Figure 3. The entropy and specific entropy for the solution of the Riemann problem 9). In Figure 3 we show the entropy and the specific entropy for the solution of Riemann problem 9). The entropy decreases as the shock and the contact discontinuity passes, while it is constant across the rarefaction wave. Analogously to the shallow water equations, we can also check whether ) holds for the solution of the Riemann problem. We know that this will hold if and only if s S + vs. Using the epression giving the shock speed, ), we calculate s S + vs = S s + v ) + l s S + v l S ) βπ + = ± l c l γ β S,
12 where we use the + sign for the first family and the sign for the second. Hence the entropy will decrease if and only if S < for the first family, and S > for the third family. Note in passing that for the contact discontinuity, s = v, and thus s S + vs = v S + v S =. Therefore, as epected, entropy is conserved across a contact discontinuity. We consider shocks of the first family, and view S as a function of π = p/p l. Recall that for these shocks π >. S = S S l p = log To check if hπ) < = h) we differentiate l ) log ) p = γ logz) logπ) ) βπ + = γ log logπ) π + β =: hπ). h β π) = γ π + β)βπ + ) π = γβ )π π + β)βπ + ) ) ππ + β)βπ + ) β = γ )πβ ) π ) ππ + β)βπ + ) β = ππ + β)βπ + ) π π ) β = ππ + β)βπ + ) π + ) <. Thus S is monotonously decreasing along the Hugoniot locus of the first family. We also see that ) holds only if p p l for waves of the first family. For shocks of the third family, an identical computation shows that ) holds only if p p l. Eercise. Show that the La inequalities v l c l > s > v p) c p), if and only if p p l, where v is given by 4), c by 5) and s is the speed of the shock, given by ) using the sign. Similarly we have v 3 p) + c 3 p) > s > v r + c r, where v 3, c 3 and s are given by 7), 8) and ) with the + sign) respectively.
Lecture 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 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 informationShock Waves. 1 Steepening of sound waves. We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: kˆ.
Shock Waves Steepening of sound waves We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: v u kˆ c s kˆ where u is the velocity of the fluid and k is the wave
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 informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationK. Ambika and R. Radha
Indian J. Pure Appl. Math., 473: 501-521, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0200-9 RIEMANN PROBLEM IN NON-IDEAL GAS DYNAMICS K. Ambika and R. Radha School of Mathematics
More informationIV. Compressible flow of inviscid fluids
IV. Compressible flow of inviscid fluids Governing equations for n = 0, r const: + (u )=0 t u + ( u ) u= p t De e = + ( u ) e= p u+ ( k T ) Dt t p= p(, T ), e=e (,T ) Alternate forms of energy equation
More information15. Eigenvalues, Eigenvectors
5 Eigenvalues, Eigenvectors Matri of a Linear Transformation Consider a linear ( transformation ) L : a b R 2 R 2 Suppose we know that L and L Then c d because of linearit, we can determine what L does
More informationChapter 5: The Hydrodynamical Riemann Problem
Chapter 5: The Hydrodynamical Riemann Problem 5.) Introduction 5..) General Introduction to the Riemann Problem We have seen in Chapter 4 that even Burgers equation, the simplest non-linear scalar conservation
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 informationShock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids
Shock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids Andrea Mentrelli Department of Mathematics & Research Center of Applied Mathematics (CIRAM) University of Bologna, Italy Summary
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 informationSimple waves and a characteristic decomposition of the two dimensional compressible Euler equations
Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute
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 informationNotes: Outline. Shallow water equations. Notes: Shallow water equations. Notes:
Outline Nonlinear hyperbolic systems Shallow water equations Shock waves and Hugoniot loci Integral curves in phase plane Compression and rarefaction R.J. LeVeque, University of Washington IPDE 2011, July
More 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 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 informationx a(x) H(U), , H(U) =
EXACT RIEMANN SOLUTIONS TO COMPRESSIBLE EULER EQUATIONS IN DUCTS WITH DISCONTINUOUS CROSS SECTION EE HAN, MAREN HANTKE, AND GERALD WARNECKE Abstract. We determine completely the exact Riemann solutions
More informationGlobal Existence of Large BV Solutions in a Model of Granular Flow
This article was downloaded by: [Pennsylvania State University] On: 08 February 2012, At: 09:55 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationShock and Expansion Waves
Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More 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 informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
More 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 informationClassification of homogeneous quadratic conservation laws with viscous terms
Volume 26, N. 2, pp. 251 283, 2007 Copyright 2007 SBMAC ISSN 0101-8205 www.scielo.br/cam Classification of homogeneous quadratic conservation laws with viscous terms JANE HURLEY WENSTROM 1 and BRADLEY
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationLecture Notes on Hyperbolic Conservation Laws
Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide
More informationVortex stretching in incompressible and compressible fluids
Vortex stretching in incompressible and compressible fluids Esteban G. Tabak, Fluid Dynamics II, Spring 00 1 Introduction The primitive form of the incompressible Euler equations is given by ( ) du P =
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 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 informationModeling of two-phase flow > Direct steam generation in solar thermal power plants
Modeling of two-phase flow > Direct steam generation in solar thermal power plants Pascal Richter RWTH Aachen University EGRIN Pirirac-sur-Mer June 3, 2015 Novatec Solar GmbH Direct steam generation Solar
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationNumerische Mathematik
Numer. Math. (2007 106:369 425 DOI 10.1007/s00211-007-0069-y Numerische Mathematik Hyperbolic balance laws: Riemann invariants and the generalized Riemann problem Matania Ben-Artzi Jiequan Li Received:
More informationAerothermodynamics of high speed flows
Aerothermodynamics of high speed flows AERO 0033 1 Lecture 6: D potential flow, method of characteristics Thierry Magin, Greg Dimitriadis, and Johan Boutet Thierry.Magin@vki.ac.be Aeronautics and Aerospace
More informationCompressible Euler equations with second sound: asymptotics of discontinuous solutions. Beixiang Fang Reinhard Racke
Universität Konstanz Compressible Euler equations with second sound: asymptotics of discontinuous solutions Beixiang Fang Reinhard Racke Konstanzer Schriften in Mathematik Nr. 306, August 0 ISSN 430-3558
More informationA Hybrid Algorithm for the Baer-Nunziato Model Using the Riemann Invariants
A Hybrid Algorithm for the Baer-Nunziato Model Using the Riemann Invariants Smadar Karni Gerardo Hernández-Dueñas Abstract The paper considers the Baer-Nunziato model for two-phase flow in porous media,
More informationComputational Fluid Dynamics. PHY 688: Numerical Methods for (Astro)Physics
Computational Fluid Dynamics Hydrodynamics When we discussed PDEs, we focused so far on scalar PDEs Often we wish to study systems of PDEs. Here we'll look at the equations of hydrodynamics Nonlinear system
More informationarxiv: v3 [math.ap] 26 May 2015
Delta shock wave for a 3 3 hyperbolic system of conservation laws arxiv:153.6693v3 [math.ap] 26 May 215 Richard De la cruz Juan Galvis Juan Carlos Juajibioy Leonardo Rendón August 27, 218 XV International
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 informationLESSON 35: EIGENVALUES AND EIGENVECTORS APRIL 21, (1) We might also write v as v. Both notations refer to a vector.
LESSON 5: EIGENVALUES AND EIGENVECTORS APRIL 2, 27 In this contet, a vector is a column matri E Note 2 v 2, v 4 5 6 () We might also write v as v Both notations refer to a vector (2) A vector can be man
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 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 informationTo subsonic flow around the wing profile with shock waves in supersonic zones. Equation for stream function.
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 1, Number 4 (016), pp. 347-351 Research India Publications http://www.ripublication.com/gjpam.htm To subsonic flow around the wing
More informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More informationWELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL
WELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL VIVIEN DESVEAUX, MARKUS ZENK, CHRISTOPHE BERTHON, AND CHRISTIAN KLINGENBERG Abstract. The present paper concerns the derivation
More information= H(W)z x, hu hu 2 + gh 2 /2., H(W) = gh
THE EXACT RIEMANN SOUTIONS TO SHAOW WATER EQUATIONS EE HAN AND GERAD WARNECKE Abstract. We determine completely the eact Riemann solutions for the shallow water equations with a bottom step including the
More informationJordan Canonical Form Homework Solutions
Jordan Canonical Form Homework Solutions For each of the following, put the matrix in Jordan canonical form and find the matrix S such that S AS = J. [ ]. A = A λi = λ λ = ( λ) = λ λ = λ =, Since we have
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 informationEigenvectors and Hermitian Operators
7 71 Eigenvalues and Eigenvectors Basic Definitions Let L be a linear operator on some given vector space V A scalar λ and a nonzero vector v are referred to, respectively, as an eigenvalue and corresponding
More informationNotes #4a MAE 533, Fluid Mechanics
Notes #4a MAE 533, Fluid Mechanics S. H. Lam lam@princeton.edu http://www.princeton.edu/ lam October 23, 1998 1 The One-dimensional Continuity Equation The one-dimensional steady flow continuity equation
More informationComputationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity:
Diagonalization We have seen that diagonal and triangular matrices are much easier to work with than are most matrices For example, determinants and eigenvalues are easy to compute, and multiplication
More informationThe Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow
The Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow D.W. Schwendeman, C.W. Wahle 2, and A.K. Kapila Department of Mathematical Sciences, Rensselaer Polytechnic
More informationUnsteady Wave Motion Shock Tube Problem - Shock Reflection
Unsteady Wave Motion Shock Tube Problem - Shock Reflection Niklas Andersson Division of Fluid Dynamics Department of Applied Mechanics Chalmers University of Tecnology 8 februari 09 Shock Reflection When
More informationFirst- and Second-order finite volume methods for the one-dimensional nonconservative Euler system
First- and Second-order finite volume methods for the one-dimensional nonconservative Euler system Stephane Clain, David Rochette To cite this version: Stephane Clain, David Rochette. First- and Second-order
More informationSimple waves and characteristic decompositions of quasilinear hyperbolic systems in two independent variables
s and characteristic decompositions of quasilinear hyperbolic systems in two independent variables Wancheng Sheng Department of Mathematics, Shanghai University (Joint with Yanbo Hu) Joint Workshop on
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationA = 3 1. We conclude that the algebraic multiplicity of the eigenvalues are both one, that is,
65 Diagonalizable Matrices It is useful to introduce few more concepts, that are common in the literature Definition 65 The characteristic polynomial of an n n matrix A is the function p(λ) det(a λi) Example
More informationA MULTIWAVE APPROXIMATE RIEMANN SOLVER FOR IDEAL MHD BASED ON RELAXATION II - NUMERICAL IMPLEMENTATION WITH 3 AND 5 WAVES
A MULTIWAVE APPROXIMATE RIEMANN SOLVER FOR IDEAL MHD BASED ON RELAXATION II - NUMERICAL IMPLEMENTATION WITH 3 AND 5 WAVES FRANÇOIS BOUCHUT CHRISTIAN KLINGENBERG 2 KNUT WAAGAN 3 Abstract. In the first part
More informationRayleigh processes in single-phase fluids
Rayleigh processes in single-phase fluids M. S. Cramer Citation: Physics of Fluids (1994-present) 18, 016101 (2006); doi: 10.1063/1.2166627 View online: http://dx.doi.org/10.1063/1.2166627 View Table of
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationHervé Guillard, INRIA Projet Smash, B.P. 93, Sophia-Antipolis Cedex, France,
TRAVELING WAVE ANALYSIS OF TWO-PHASE DISSIPATIVE MODELS Hervé Guillard, INRIA Projet Smash, B.P. 93, 06902 Sophia-Antipolis Cedex, France, Herve.Guillard@sophia.inria.fr Joint work with : Mathieu Labois,
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationarxiv: v1 [math.ap] 29 May 2018
Non-uniqueness of admissible weak solution to the Riemann problem for the full Euler system in D arxiv:805.354v [math.ap] 9 May 08 Hind Al Baba Christian Klingenberg Ondřej Kreml Václav Mácha Simon Markfelder
More informationCouplage de modèles en thermohydraulique. Troisième rapport de synthèse
Couplage de modèles en thermohydraulique. Troisième rapport de synthèse Annalisa Ambroso, Benjamin Boutin, Christophe Chalons, Frédéric Coquel, Thomas Galié Edwige Godlewski, Frédéric Lagoutière, Pierre-Arnaud
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationLinear Algebra Final Exam Solutions, December 13, 2008
Linear Algebra Final Exam Solutions, December 13, 2008 Write clearly, with complete sentences, explaining your work. You will be graded on clarity, style, and brevity. If you add false statements to a
More informationAnnalisa Ambroso 1, Christophe Chalons 2, Frédéric Coquel 3 and Thomas. Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher RELAXATION AND NUMERICAL APPROXIMATION OF A TWO-FLUID TWO-PRESSURE DIPHASIC MODEL
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 informationChapter 4: Non-Linear Conservation Laws; the Scalar Case. The conservation law in eqn. (4.1) will be hyperbolic if its eigenvalue ( ) ( )
Chapter 4: Non-inear Conservation aws; the Scalar Case 4.1) Introduction In the previous chapter we developed an understanding of monotonicity preserving advection schemes and iemann solvers for linear
More 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 informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationCharbel Farhat, Arthur Rallu, Sriram Shankaran
A Higher-Order Generalized Ghost Fluid Method for the Poor for the Three-Dimensional Two-Phase Flow Computation of Underwater Eplosions and Implosions Charbel Farhat, Arthur Rallu, Sriram Shankaran Department
More informationHyperbolicity singularities in rarefaction waves
Hyperbolicity singularities in rarefaction waves Alexei A. Mailybaev and Dan Marchesin Abstract For mixed-type systems of conservation laws, rarefaction waves may contain states at the boundary of the
More informationThe Euler Equations! Advection! λ 1. λ 2. λ 3. ρ ρu. c 2 = γp/ρ. E = e + u 2 /2 H = h + u 2 /2; h = e + p/ρ. 0 u 1/ρ. u p. t + A f.
http://www.nd.edu/~gtryggva/cfd-course/! Advection! Grétar Tryggvason! Spring! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = x ( / ) ρe ρu E + p ρ E = e + u / H = h + u
More informationThe RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
More informationDesigning Information Devices and Systems I Spring 2016 Official Lecture Notes Note 21
EECS 6A Designing Information Devices and Systems I Spring 26 Official Lecture Notes Note 2 Introduction In this lecture note, we will introduce the last topics of this semester, change of basis and diagonalization.
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More 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 informationCHAPTER 1 The Compressible Euler System in Two Space Dimensions. Introduction
1 CHAPTER 1 The Compressible Euler System in Two Space Dimensions Yuxi Zheng Department of Mathematics The Pennsylvania State University E-mail: yzheng@math.psu.edu Abstract We analyze the self-similar
More informationSingularity formation for compressible Euler equations
Singularity formation for compressible Euler equations Geng Chen Ronghua Pan Shengguo Zhu Abstract In this paper, for the p-system and full compressible Euler equations in one space dimension, we provide
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationIntroduction to Turbomachinery
1. Coordinate System Introduction to Turbomachinery Since there are stationary and rotating blades in turbomachines, they tend to form a cylindrical form, represented in three directions; 1. Axial 2. Radial
More informationStability of Patterns
Stability of Patterns Steve Schecter Mathematics Department North Carolina State University 1 2 Overview Motivation Systems of viscous conservation laws and their variants are partial differential equations
More informationPhysics Dec Time Independent Solutions of the Diffusion Equation
Physics 301 10-Dec-2004 33-1 Time Independent Solutions of the Diffusion Equation In some cases we ll be interested in the time independent solution of the diffusion equation Why would be interested in
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationA New Relativistic Hydrodynamic Code
A New Relativistic Hydrodynamic Code Dongsu Ryu, Indranil Chattopadhyay Department of Astronomy and Space Science, Chungnam National University, Daejeon 305-764, Korea: ryu@canopus.cnu.ac.kr, indra@canopus.cnu.ac.kr
More informationRepeated Eigenvalues and Symmetric Matrices
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationSTRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION
Electronic Journal of Differential Equations, Vol. 216 (216, No. 126, pp. 1 16. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN
More informationLU Factorization. A m x n matrix A admits an LU factorization if it can be written in the form of A = LU
LU Factorization A m n matri A admits an LU factorization if it can be written in the form of Where, A = LU L : is a m m lower triangular matri with s on the diagonal. The matri L is invertible and is
More informationMAS. Modelling, Analysis and Simulation. Modelling, Analysis and Simulation. Building your own shock tube. J. Naber
C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a MAS Modelling, Analysis and Simulation Modelling, Analysis and Simulation Building your own shock tube J. Naber REPORT MAS-E5 JANUARY 5 CWI
More informationA Linearised Riemann Solver for Godunov-Type Methods. E.F.Toro
Cranfield»nun «32 College of Aeronautics Report No.9116 September 1991 TECHNISCHE UNIVERSITEIT DELFT LUCHTVAART- EN RUIMTEVAARnECHNIEK BIBLIOTHEEK Kluyverweg 1-2629 HS DELFT A Linearised Riemann Solver
More information2 Eigenvectors and Eigenvalues in abstract spaces.
MA322 Sathaye Notes on Eigenvalues Spring 27 Introduction In these notes, we start with the definition of eigenvectors in abstract vector spaces and follow with the more common definition of eigenvectors
More informationIn this chapter we will be concerned with compressible gas flows. This forms the main focus of this lecture. Name Symbol Unit (CGS) Unit (SI)
Chapter 1 Equations of hydrodynamics In this chapter we will be concerned with compressible gas flows. This forms the main focus of this lecture. 1.1 Basic quantities The basic quantities that describe
More informationOn Fluid Maxwell Equations
On Fluid Maxwell Equations Tsutomu Kambe, (Former Professor, Physics), University of Tokyo, Abstract Fluid mechanics is a field theory of Newtonian mechanics of Galilean symmetry, concerned with fluid
More informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 1 I. DERIVATION OF THE DISCRETIZED PATH INTEGRAL We begin our discussion of the Feynman path integral with the canonical ensemble. The epressions for the
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More information