Shock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids

Size: px
Start display at page:

Download "Shock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids"

Transcription

1 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

2 Summary Compressible vs. incompressible fluids Thermodynamic restrictions Models of incompressibility (perfect incompressibility, QTI model, EQTI model) An example of EQTI equation of Shock and rarefaction waves in EQTI materials The Riemann problem in EQTI materials

3 Introduction Incompressibility is a useful idealization when modeling materials characterized by extreme resistance to volume changes Incompressible fluids are treated as the limit case of compressible fluids Purely mechanical problems The limit is fine :) Solutions of incompressible model equations are obtained as limit of compressible ones (P.-L Lions, N. Masmoudi, J. Math. Pures Appl. 77, , 1998; S. Klainerman, S. Majda, Comm. Pure Appl. Math. 35, , 1982; etc. etc.) Thermo-mechanical problems The limit is ambiguous :( There are several definitions of incompressibility (perfect incompressibility; quasithermal-incompressibility; extended-quasi-thermal-incompressibility)

4 Introduction From a mathematical point of view, compressible and incompressible materials have different treatments: compressible materials the pressure is a constitutive function incompressible materials the pressure is a Lagrange multiplier (associated to the constraint of incompressibility) In order to compare the solutions of compressible/incompressible materials, it is convenient to choose the pressure (instead of the density) as a thermodynamic variable along with temperature

5 Introduction Euler equations (perfect fluid in the absence of external body forces) (conservation of mass) (conservation of momentum) (conservation of energy) Independent variables: pressure temperature velocity components Equations of : density (or specific volume) specific internal energy

6 Thermodynamic restrictions The constitutive equations of a material must satisfy the entropy principle and the conditions for the thermodynamic stability 1. Entropy principle Gibbs equation (entropy) Chemical potential Once the thermal equation of is known: with or

7 Thermodynamic restrictions 2. Thermodynamic stability Enthalpy Specific heat The chemical potential must be a concave function of pressure and temperature: Thus, thermodinamic stability requires:

8 Thermodynamic restrictions After introducing - as usual when studying the compressibility of a material - the two parameters thermal expansion coefficient compressibility factor The concavity of the chemical potential requires Thus, the thermodinamic stability conditions are Adiabatic sound velocity:

9 Models of incompressibility Perfectly incompressible material According to this model, a material is incompressible when its constitutive functions do not depend on the pressure It is easily seen that: the specific volume must be constant! (Müller's paradox) (I. Müller, Thermodynamics, Pitman, London, 1985)

10 Models of incompressibility Quasi-Thermal-Incompressible (QTI) Material According to this model, only the density is assumed to be indipendent on the pressure In this case: If then The Müller's paradox is removed, and perfectly incompressible materials are obtained as a limit of QTI materials The chemical potential is not concave: QTI materials are thermodynamically unstable (K. R. Rajagopal et al., Math. Model. Meth. Appl. Sci. 6, , 1996; H. Gouin, A. Muracchini, T. Ruggeri, Continuum Mech. Thermodyn. 24, , 2012)

11 Models of incompressibility Extended-Quasi-Thermal-Incompressible (EQTI) Material A compressible material satisfying the thermodynamic conditions is called an extendedquasi-thermal-incompressible (EQTI) material if: with An EQTI material is a stable compressible material that approximates a Müller incompressible material to the second order Assuming reference It is easily obtained:

12 Models of incompressibility Extended-Quasi-Thermal-Incompressible (EQTI) Material Thermodynamic stability is guaranteed if If And EQTI materials approximate a perfectly incompressible material (if the pressure is smaller than a critical value) Thermodynamically stability is guaranteed (H. Gouin, T. Ruggeri, Internat. J. Non-Linear Mech. 47, , 2012; A. M., T. Ruggeri, Internat. J. Non-Linear Mech. 51, 87-96, 2013)

13 An example of EQTI fluid We perform a linear expansion of the specific volume near the reference Modified Bousinnesq equation is an equation of of a EQTI material with:

14 An example of EQTI fluid Modified Bousinnesq equation Thermal expansion coefficient and compressibility factor: Condition for thermodynamic stability

15 Shock waves in EQTI fluids Shock wave, perturbed shock velocity, s unperturbed Euler equations (1D) Eigenvalues and eigenvectors

16 Shock waves in EQTI fluids Shock wave, perturbed shock velocity, s unperturbed Euler equations (1D) Rankine-Hugoniot conditions Hugoniot locus One-parameter family of solutions (parameter: ) Selection rule?

17 Shock waves in EQTI fluids Shock wave, perturbed shock velocity, s unperturbed Euler equations (1D) Genuinely non-linear waves Hugoniot locus Lax condition: Linearly degenerate waves Lax condition: ble ssi i dm ina

18 Shock waves in EQTI fluids Shock wave, perturbed shock velocity, s unperturbed Euler equations (1D) Locally linearly degenerate waves for some Hugoniot locus ib iss m ad in le Liu condition: ina d mis sib le through

19 Shock waves in EQTI fluids Shock wave, perturbed shock velocity, s unperturbed Euler equations (1D) Rankine-Hugoniot conditions One-parameter family of solutions (parameter: )

20 Shock waves in EQTI fluids Shock waves in water specific volume fluid velocity temperature shock velocity

21 Shock waves in EQTI fluids Shock waves in the incompressible limit specific volume fluid velocity temperature shock velocity

22 Rarefaction waves in EQTI fluids Rarefaction wave, perturbed unperturbed Euler equations (1D) The Rarefaction wave is a similarity solution

23 Rarefaction waves in EQTI fluids Rarefaction wave, perturbed unperturbed Euler equations (1D)

24 Rarefaction waves in EQTI fluids Rarefaction wave, perturbed unperturbed Euler equations (1D)

25 Rarefaction waves in EQTI fluids Rarefaction wave, perturbed unperturbed Euler equations (1D) Taking into account that for a EQTI fluid After some manipulations...

26 Rarefaction waves in EQTI fluids Rarefaction wave, perturbed unperturbed Euler equations (1D) The integral curves are defined by

27 Rarefaction waves in EQTI fluids Rarefaction waves in water specific volume fluid velocity temperature 'development' velocity

28 Rarefaction waves in EQTI fluids Rarefaction waves in the incompressible limit specific volume fluid velocity temperature 'development' velocity

29 The Riemann Problem Euler equations (1D) Typical solution: Two nonlinear waves (shock/rarefaction waves) Contact discontinuity

30 The Riemann Problem The knowledge of the Rankine-Hugoniot curves of the shock waves and of the integral curves of the rarefaction waves allow to completely solve the Riemann problem Numerical results

31 References & Acknowledgments References I. Müller, Thermodynamics, Pitman, London (1985) H. Gouin, A. Muracchini, T. Ruggeri, Continuum Mech. Thermodyn. 24, (2012) H. Gouin, T. Ruggeri, Internat. J. Non-Linear Mech. 47, (2012) A. M., T. Ruggeri, Internat. J. Non-Linear Mech. 51, (2013) Acknowledgments Partially supported by GNFM/INdAM Young Researchers Project 2012 Hyperbolic Models for Incompressible Materials

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3

Non-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 information

Lecture 5.7 Compressible Euler Equations

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 information

Shock 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. 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 information

Various lecture notes for

Various 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 information

Simple 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 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 information

On the Müller paradox for thermal-incompressible media

On the Müller paradox for thermal-incompressible media On the Müller paradox for thermal-incompressible media Henri Gouin, Augusto Muracchini, Tommaso Ruggeri To cite this version: Henri Gouin, Augusto Muracchini, Tommaso Ruggeri. On the Müller paradox for

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

AA214B: 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 information

Nonlinear stability of compressible vortex sheets in two space dimensions

Nonlinear stability of compressible vortex sheets in two space dimensions of compressible vortex sheets in two space dimensions J.-F. Coulombel (Lille) P. Secchi (Brescia) CNRS, and Team SIMPAF of INRIA Futurs Evolution Equations 2006, Mons, August 29th Plan 1 2 3 Related problems

More information

E = where γ > 1 is a constant spesific to the gas. For air, γ 1.4. Solving for p, we get. 2 ρv2 + (γ 1)E t

E = where γ > 1 is a constant spesific to the gas. For air, γ 1.4. Solving for p, we get. 2 ρv2 + (γ 1)E t . 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,

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

AA214B: 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 information

Notes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow

Notes: 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 information

STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION

STRUCTURAL 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 information

Rarefaction wave interaction for the unsteady transonic small disturbance equations

Rarefaction wave interaction for the unsteady transonic small disturbance equations Rarefaction wave interaction for the unsteady transonic small disturbance equations Jun Chen University of Houston Department of Mathematics 4800 Calhoun Road Houston, TX 77204, USA chenjun@math.uh.edu

More information

The RAMSES code and related techniques I. Hydro solvers

The 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 information

CapSel Roe Roe solver.

CapSel 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 information

Notes: Outline. Shallow water equations. Notes: Shallow water equations. Notes:

Notes: 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 information

On the Dependence of Euler Equations on Physical Parameters

On 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 information

On the Cauchy Problems for Polymer Flooding with Gravitation

On the Cauchy Problems for Polymer Flooding with Gravitation On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University. Email: wxs27@psu.edu November 5, 2015 Abstract We study two systems of conservation

More information

A Very Brief Introduction to Conservation Laws

A 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 information

Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation

Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation W.L. Jin and H.M. Zhang August 3 Abstract: In this paper we study the Payne-Whitham (PW) model as

More information

The 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy

The 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University,

More information

The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method

The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas

More information

CapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x

CapSel 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 information

u-= (u, v), x>o, j u, (u,, v,), x<o, U(X 0) (1) (1), A A2 only when u 0, in which case A 0. THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY II*

u-= (u, v), x>o, j u, (u,, v,), x<o, U(X 0) (1) (1), A A2 only when u 0, in which case A 0. THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY II* SIAM J. APPL. MATH. Vol. 48, No. 6, December 1988 (C) 1988 Society for Industrial and Applied Mathematics 006 THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY II* E. ISAACSONS" AN[) B. TEMPLE:I: Abstract.

More information

Simple waves and characteristic decompositions of quasilinear hyperbolic systems in two independent variables

Simple 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 information

Low Mach number limit of the full Navier-Stokes equations

Low Mach number limit of the full Navier-Stokes equations Low Mach number limit of the full NavierStokes equations Thomas Alazard Mathématiques Appliquées de Bordeaux Université Bordeaux 1 351 cours de la Libération 33405 Talence FRANCE thomas.alazard@math.ubordeaux.fr

More information

Global Existence of Large BV Solutions in a Model of Granular Flow

Global 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 information

AMath 574 February 11, 2011

AMath 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 information

Gas Dynamics Equations: Computation

Gas Dynamics Equations: Computation Title: Name: Affil./Addr.: Gas Dynamics Equations: Computation Gui-Qiang G. Chen Mathematical Institute, University of Oxford 24 29 St Giles, Oxford, OX1 3LB, United Kingdom Homepage: http://people.maths.ox.ac.uk/chengq/

More information

Low Froude Number Limit of the Rotating Shallow Water and Euler Equations

Low Froude Number Limit of the Rotating Shallow Water and Euler Equations Low Froude Number Limit of the Rotating Shallow Water and Euler Equations Kung-Chien Wu Department of Pure Mathematics and Mathematical Statistics University of Cambridge, Wilberforce Road Cambridge, CB3

More information

Hyperbolic 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 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 information

IV. Compressible flow of inviscid fluids

IV. 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 information

Fluid flows through unsaturated porous media: An alternative simulation procedure

Fluid flows through unsaturated porous media: An alternative simulation procedure Engineering Conferences International ECI Digital Archives 5th International Conference on Porous Media and Their Applications in Science, Engineering and Industry Refereed Proceedings Summer 6-24-2014

More information

A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS

A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS KATARINA JEGDIĆ, BARBARA LEE KEYFITZ, AND SUN CICA ČANIĆ We study a Riemann problem for the two-dimensional isentropic gas dynamics equations

More information

Shock and Expansion Waves

Shock 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 information

Hervé Guillard, INRIA Projet Smash, B.P. 93, Sophia-Antipolis Cedex, France,

Hervé 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 information

Non-linear Scalar Equations

Non-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 information

Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases

Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases Heinrich Freistühler and Blake Temple Proceedings of the Royal Society-A May 2017 Culmination of a 15 year project: In this we propose:

More information

APPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS

APPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS 1 APPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS BIN CHENG Department of Mathematics University of Michigan Ann Arbor, MI 48109 E-mail:bincheng@umich.edu EITAN

More information

Chapter Two. Basic Thermodynamics, Fluid Mechanics: Definitions of Efficiency. Laith Batarseh

Chapter Two. Basic Thermodynamics, Fluid Mechanics: Definitions of Efficiency. Laith Batarseh Chapter Two Basic Thermodynamics, Fluid Mechanics: Definitions of Efficiency Laith Batarseh The equation of continuity Most analyses in this book are limited to one-dimensional steady flows where the velocity

More information

Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory

Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Gennady El 1, Roger Grimshaw 1 and Noel Smyth 2 1 Loughborough University, UK, 2 University of Edinburgh, UK

More information

Global Riemann Solver and Front Tracking Approximation of Three-Component Gas Floods

Global Riemann Solver and Front Tracking Approximation of Three-Component Gas Floods Global Riemann Solver and Front Tracking Approximation of Three-Component Gas Floods Saeid Khorsandi (1), Wen Shen (2) and Russell T. Johns (3) (1) Department of Energy and Mineral Engineering, Penn State

More information

Singularity formation for compressible Euler equations

Singularity 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 information

Numerical 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 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 information

Conical Shock Waves for Isentropic Euler System

Conical Shock Waves for Isentropic Euler System Conical Shock Waves for Isentropic Euler System Shuxing Chen Institute of Mathematical Research, Fudan University, Shanghai, China E-mail: sxchen@public8.sta.net.cn Dening Li Department of Mathematics,

More information

Computational Fluid Dynamics. PHY 688: Numerical Methods for (Astro)Physics

Computational 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 information

u-- (u, v), x < 0, v,+1/2{bu:+2uv},,=o u--(u,v), (1) (1) and note that A1 A2 only when u 0 in which case A 0.

u-- (u, v), x < 0, v,+1/2{bu:+2uv},,=o u--(u,v), (1) (1) and note that A1 A2 only when u 0 in which case A 0. SIAM J. APPL. MATH. Vol. 48, No. 6, December 1988 1988 Society for Industrial and Applied Mathematics OO7 THE RIEMANN PROBLEM NEAR A HYPERBOLIC SINGULARITY III* E. ISAACSON AND B. TEMPLE Abstract. This

More information

NONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS

NONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL

More information

Quasi-linear first order equations. Consider the nonlinear transport equation

Quasi-linear first order equations. Consider the nonlinear transport equation Quasi-linear first order equations Consider the nonlinear transport equation u t + c(u)u x = 0, u(x, 0) = f (x) < x < Quasi-linear first order equations Consider the nonlinear transport equation u t +

More information

Perturbation Theory 1

Perturbation Theory 1 Perturbation Theory 1 1 Expansion of Complete System Let s take a look of an expansion for the function in terms of the complete system : (1) In general, this expansion is possible for any complete set.

More information

0.3.4 Burgers Equation and Nonlinear Wave

0.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 information

On the Front-Tracking Algorithm

On the Front-Tracking Algorithm JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 7, 395404 998 ARTICLE NO. AY97575 On the Front-Tracking Algorithm Paolo Baiti S.I.S.S.A., Via Beirut 4, Trieste 3404, Italy and Helge Kristian Jenssen

More information

The Riemann problem. The Riemann problem Rarefaction waves and shock waves

The 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 information

Computational Astrophysics

Computational Astrophysics 16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing

More information

Thermodynamic Third class Dr. Arkan J. Hadi

Thermodynamic Third class Dr. Arkan J. Hadi 5.5 ENTROPY CHANGES OF AN IDEAL GAS For one mole or a unit mass of fluid undergoing a mechanically reversible process in a closed system, the first law, Eq. (2.8), becomes: Differentiation of the defining

More information

Rayleigh processes in single-phase fluids

Rayleigh 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 information

Answers to Problem Set Number 04 for MIT (Spring 2008)

Answers 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 information

On a simple model of isothermal phase transition

On a simple model of isothermal phase transition On a simple model of isothermal phase transition Nicolas Seguin Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie Paris 6 France Micro-Macro Modelling and Simulation of Liquid-Vapour Flows

More information

International Engineering Research Journal

International 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 information

The Euler Equation of Gas-Dynamics

The 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 information

K. Ambika and R. Radha

K. 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 information

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

CHAPTER 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 information

Shock on the left: locus where cars break behind the light.

Shock on the left: locus where cars break behind the light. Review/recap of theory so far. Evolution of wave profile, as given by the characteristic solution. Graphical interpretation: Move each point on graph at velocity c(ρ). Evolution as sliding of horizontal

More information

Numerical Methods for Hyperbolic Conservation Laws Lecture 4

Numerical 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 information

A 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 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 information

Relativistic Gases. 1 Relativistic gases in astronomy. Optical. Radio. The inner part of M87. Astrophysical Gas Dynamics: Relativistic Gases 1/73

Relativistic Gases. 1 Relativistic gases in astronomy. Optical. Radio. The inner part of M87. Astrophysical Gas Dynamics: Relativistic Gases 1/73 Relativistic Gases 1 Relativistic gases in astronomy Optical The inner part of M87 Radio Astrophysical Gas Dynamics: Relativistic Gases 1/73 Evidence for relativistic motion Motion of the knots in the

More information

Stability of Thick Spherical Shells

Stability of Thick Spherical Shells Continuum Mech. Thermodyn. (1995) 7: 249-258 Stability of Thick Spherical Shells I-Shih Liu 1 Instituto de Matemática, Universidade Federal do Rio de Janeiro Caixa Postal 68530, Rio de Janeiro 21945-970,

More information

2 In this paper, we give a brief summary of the front tracking algorithm for axisymmetric ow. We report simulations of spherical shock refraction by a

2 In this paper, we give a brief summary of the front tracking algorithm for axisymmetric ow. We report simulations of spherical shock refraction by a A Fast Algorithm for Moving Interface Problems S. Dutta 1,J.Glimm 1 2,J.W.Grove 3,D.H.Sharp 4, and Y. Zhang 1 1 Department of Applied Mathematics and Statistics, University at Stony Brook, Stony Brook,

More information

Can constitutive relations be represented by non-local equations?

Can constitutive relations be represented by non-local equations? Can constitutive relations be represented by non-local equations? Tommaso Ruggeri Dipartimento di Matematica & Centro di Ricerca per le Applicazioni della Matematica (CIRAM) Universitá di Bologna Fractional

More information

Author(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)

Author(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1) Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL

More information

Introduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)

Introduction 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 information

THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS

THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS Journal of Hyperbolic Differential Equations Vol., No. 4 005 909 917 c World Scientific Publishing Company THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS VOLKER ELLING, and TAI-PING LIU, Dept.

More information

Applications of the compensated compactness method on hyperbolic conservation systems

Applications of the compensated compactness method on hyperbolic conservation systems Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,

More information

A PARADIGM FOR TIME-PERIODIC SOUND WAVE PROPAGATION IN THE COMPRESSIBLE EULER EQUATIONS

A PARADIGM FOR TIME-PERIODIC SOUND WAVE PROPAGATION IN THE COMPRESSIBLE EULER EQUATIONS METHODS AND APPLICATIONS OF ANALYSIS. c 2009 International Press Vol. 16, No. 3, pp. 341 364, September 2009 005 A PARADIGM FOR TIME-PERIODIC SOUND WAVE PROPAGATION IN THE COMPRESSIBLE EULER EQUATIONS

More information

Shock Reflection-Diffraction, Nonlinear Partial Differential Equations of Mixed Type, and Free Boundary Problems

Shock Reflection-Diffraction, Nonlinear Partial Differential Equations of Mixed Type, and Free Boundary Problems Chapter One Shock Reflection-Diffraction, Nonlinear Partial Differential Equations of Mixed Type, and Free Boundary Problems Shock waves are steep fronts that propagate in compressible fluids when convection

More information

Anomalous wave structure in magnetized materials described by non-convex equations of state

Anomalous wave structure in magnetized materials described by non-convex equations of state This is a preprint of: Anomalous wave structure in magnetized materials described by non-convex equations of state, Susana Serna, Antonio Marquina, Phys. Fluids, vol. 6, 6, 4. DOI: [.63/.48545] Anomalous

More information

The RAMSES code and related techniques 2- MHD solvers

The 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 information

Chapter 5: The Hydrodynamical Riemann Problem

Chapter 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 information

Modelling and numerical methods for the diffusion of impurities in a gas

Modelling and numerical methods for the diffusion of impurities in a gas INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi

More information

Waves in a Shock Tube

Waves 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 information

L 1 Stability for Systems of Hyperbolic Conservation Laws

L 1 Stability for Systems of Hyperbolic Conservation Laws Contemporary Mathematics Volume 238, 1999 B 0-8218-1196-7-03547-7 L 1 Stability for Systems of Hyperbolic Conservation Laws Tai-Ping Liu and Tong Yang Abstract. In this paper, we summarize our results

More information

Central Schemes for Systems of Balance Laws Salvatore Fabio Liotta, Vittorio Romano, Giovanni Russo Abstract. Several models in mathematical physics a

Central Schemes for Systems of Balance Laws Salvatore Fabio Liotta, Vittorio Romano, Giovanni Russo Abstract. Several models in mathematical physics a Central Schemes for Systems of Balance Laws Salvatore Fabio Liotta, Vittorio Romano, Giovanni Russo Abstract. Several models in mathematical physics are described by quasilinear hyperbolic systems with

More information

A maximum principle for optimally controlled systems of conservation laws

A maximum principle for optimally controlled systems of conservation laws RENDICONTI del SEMINARIO MATEMATICO della UNIVERSITÀ DI PADOVA ALBERTO BRESSAN ANDREA MARSON A maximum principle for optimally controlled systems of conservation laws Rendiconti del Seminario Matematico

More information

The Riemann problem for a class of resonant hyperbolic systems of balance laws

The Riemann problem for a class of resonant hyperbolic systems of balance laws The Riemann problem for a class of resonant hyperbolic systems of balance laws Paola Goatin and Philippe G. LeFloch Abstract. We solve the Riemann problem for a class of resonant hyperbolic systems of

More information

Initial Boundary Value Problems for Scalar and Vector Burgers Equations

Initial Boundary Value Problems for Scalar and Vector Burgers Equations Initial Boundary Value Problems for Scalar and Vector Burgers Equations By K. T. Joseph and P. L. Sachdev In this article we stu Burgers equation and vector Burgers equation with initial and boundary conditions.

More information

Shock formation in the compressible Euler equations and related systems

Shock formation in the compressible Euler equations and related systems Shock formation in the compressible Euler equations and related systems Geng Chen Robin Young Qingtian Zhang Abstract We prove shock formation results for the compressible Euler equations and related systems

More information

Numerical Solutions to Partial Differential Equations

Numerical 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 information

A LIAPUNOV-SCHMIDT REDUCTION FOR TIME-PERIODIC SOLUTIONS OF THE COMPRESSIBLE EULER EQUATIONS

A LIAPUNOV-SCHMIDT REDUCTION FOR TIME-PERIODIC SOLUTIONS OF THE COMPRESSIBLE EULER EQUATIONS A LIAPUNOV-SCHMIDT REDUCTION FOR TIME-PERIODIC SOLUTIONS OF THE COMPRESSIBLE EULER EQUATIONS BLAKE TEMPLE AND ROBIN YOUNG Abstract. Following the authors earlier work in [9, 10], we show that the nonlinear

More information

Shock reflection and oblique shock waves

Shock reflection and oblique shock waves JOURNAL OF MATHEMATICAL PHYSICS 48, 12312 27 Shock reflection and oblique shock waves Dening Li a Department of Mathematics, West Virginia University, Morgantown, West Virginia 2656, USA Received 6 March

More information

Introduction to Fluid Mechanics. Chapter 13 Compressible Flow. Fox, Pritchard, & McDonald

Introduction to Fluid Mechanics. Chapter 13 Compressible Flow. Fox, Pritchard, & McDonald Introduction to Fluid Mechanics Chapter 13 Compressible Flow Main Topics Basic Equations for One-Dimensional Compressible Flow Isentropic Flow of an Ideal Gas Area Variation Flow in a Constant Area Duct

More information

Stability of Mach Configuration

Stability of Mach Configuration Stability of Mach Configuration Suxing CHEN Fudan University sxchen@public8.sta.net.cn We prove the stability of Mach configuration, which occurs in moving shock reflection by obstacle or shock interaction

More information

Math Partial Differential Equations 1

Math 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 information

Ray equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions

Ray equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions Proc. Indian Acad. Sci. (Math. Sci.) Vol. 126, No. 2, May 2016, pp. 199 206. c Indian Academy of Sciences Ray equations of a weak in a hyperbolic system of conservation laws in multi-dimensions PHOOLAN

More information

THE EULER EQUATIONS OF COMPRESSIBLE FLUID FLOW

THE EULER EQUATIONS OF COMPRESSIBLE FLUID FLOW BULLETIN (New Series) OF THE AMERICAN MATHEMATICAL SOCIETY Volume 44, Number 4, October 2007, Pages 581 602 S 0273-0979(07)01181-0 Article electronically published on June 18, 2007 THE EULER EQUATIONS

More information

Math 660-Lecture 23: Gudonov s method and some theories for FVM schemes

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 information

Finite Volume simulation of cavitating flows

Finite Volume simulation of cavitating flows Finite Volume simulation of cavitating flows Thomas Barberon, Philippe Helluy ISITV, Laboratoire MNC, BP 56, 83162 La Valette, France Abstract We propose a numerical method adapted to the modelling of

More information

Riemann Solvers and Numerical Methods for Fluid Dynamics

Riemann 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 information

COMPUTATIONAL METHODS AND ALGORITHMS Vol II - Computational Methods for Compressible Flow Problems - Remi Abgrall

COMPUTATIONAL METHODS AND ALGORITHMS Vol II - Computational Methods for Compressible Flow Problems - Remi Abgrall COMPUTATIONAL METHODS FOR COMPRESSIBLE FLOW PROBLEMS Rémi Abgrall Université Bordeaux I, Talence Cedex, France Keywords: computational methods, numerical schemes, 1- D problems, multidimensional problems,

More information

Notes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion

Notes: 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 information

Hyperbolic Conservation Laws Past and Future

Hyperbolic Conservation Laws Past and Future Hyperbolic Conservation Laws Past and Future Barbara Lee Keyfitz Fields Institute and University of Houston bkeyfitz@fields.utoronto.ca Research supported by the US Department of Energy, National Science

More information