MATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS

Size: px
Start display at page:

Download "MATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS"

Transcription

1 K CHAPMAN & HALL/CRC Monographs and Surveys in Pure and Applied Mathematics I 18 MATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS ANDREI G. KULIKOVSKII NIKOLAI V. POGORELOV ANDREI YU. SEMENOV CHAPMAN & HALL/CRC Boca Raton London New York Washington, D.C.

2 Contents 1 Hyperbolic Systems of Partial Differential Equations Quasilinear systems Hyperbolic systems of quasilinear differential equations Definitions Systems of conservation laws Mechanical examples Nonstationary equations of gas dynamics Stationary Euler equations Shallow water equations Equations of ideal magnetohydrodynamics Elasticity equations Properties of solutions Classical solutions Generalized solutions Small-amplitude shocks Evolutionary conditions for shocks Entropy behavior on discontinuities Disintegration of a small arbitrary discontinuity 31 2 Numerical Solution of Quasilinear Hyperbolic Systems "introduction Methods based on the exact solution of the Riemann problem The Godunov method of the first order Exact solution of the Riemann problem Methods based on approximate Riemann problem solvers Courant-Isaacson-Rees-type methods Roe's scheme The Osher numerical scheme Generalized Riemann problem The Godunov method of the second order Multidimensional scheihes and their stability conditions Reconstruction procedures and slope limiters Preliminary remarks TVD schemes Monotone and limiting reconstructions 80

3 2.7.4 Genuine TVD and TVD limiting reconstructions TVD limiters of nonsymmetric stencil Multidimensional reconstruction Boundary conditions for hyperbolic systems General notions Nonreflecting boundary conditions Evolutionary boundary conditions Shock-fitting methods Floating shock fitting Shock fitting on moving grids Entropy correction procedures Final remarks 119 Gas Dynamic Equations Systems of governing equations Two-temperature gas dynamic equations The mixture of ideal gases in chemical nonequilibrium The Godunov method for gas dynamic equations Exact solution of the Riemann problem Elementary solution 1: Shockwave Elementary solution 2: Contact discontinuity Elementary solution 3: Rarefaction wave General exact solution An arbitrary EOS Approximate Riemann problem solvers The Courant-Isaacson-Rees method for an arbitrary EOS Computation of shock-induced phenomena by the CIR method The CIR-simulation of jet-like structures in laser plasma Roe's method Roe's Riemann problem solver for an arbitrary EOS Osher-Solomon numerical scheme Shock-fitting methods Discontinuities as boundaries of the computational region Floating shock-fitting procedures Shock-fitting on moving grids Self-adjusting grids Stationary gas dynamics Systems of governing equations The Godunov method. The CIR and Roe's schemes Exact solution of the Riemann problem General exact solution Solar wind - interstellar medium interaction Physical formulation of the problem Nonreflecting boundary conditions Numerical results 221

4 3.7.4 A note on Godunov-type methods for relativistic hydrodynamics Shallow Water Equations System of governing equations The Godunov method for shallow water equations Exact solution of the Riemann problem v Elementary solution 1: Hydraulic jump Elementary solution 2: Tangential discontinuity Elementary solution 3: Riemann wave General exact solution Results of numerical analysis Approximate Riemann problem solvers The CIR method Roe's method The Osher-Solomon solver Stationary shallow water equations System of governing equations The Godunov method. The CIR and Roe's schemes Exact solution of the Riemann problem General exact solution 275 Magnetohydrodynamic Equations MHD system in the conservation-law form Classification of MHD discontinuities Evolutionary MHD shocks Evolutionary diagram Convenient relations on MHD shocks Evolutionarity of perpendicular, parallel, and singular shocks Jouget points High-resolution numerical schemes for MHD equations The Osher-type method Piecewise-parabolic method Roe's characteristic decomposition method Numerical tests with the Roe-type scheme Modified MHD system Shock-capturing approach and nonevolutionary solutions in MHD Preliminary remarks Simplified MHD equations and related discontinuities Shock structure in solutions of the simplified system Nonstationary processes in the structure of nonevolutionary shock waves Numerical experiments based on the full set of MHD equations Numerical disintegration of a compound wave Strong background magnetic field Elimination of numerical magnetic charge 348

5 Xll Preliminary remarks Application of the vector potential The use of an artificial scalar potential _ Application of the modified MHD system Application of staggered grids Solar wind interaction with the magnetized interstellar medium Statement of the problem Numerical algorithm Numerical results: axisymmetric case Numerical results: rotationally perturbed flow A note on the MHDflowover an infinitely conducting cylinder Numerical results: three-dimensional modelling Solid Dynamics Equations System of governing equations Solid dynamics with an arbitrary EOS Conservative form of elastoviscoplastic solid dynamics Dynamics of thin shells CIR method for the calculation of solid dynamics problems Numerical simulation of spallation phenomena CIR method for studying the dynamics of thin shells The Klein-Gordon equation Dynamics equations of cylindrical shells.... x Dynamics equations of orthotropic shells Selection of rapidly oscillating components Nonclassical Discontinuities and Solutions of Hyperbolic Systems Evolutionary conditions in nonclassical cases Structure of fronts. Additional boundary conditions on the fronts Equations describing the discontinuity structure Formulation of the structure problem and additional assumptions :2.3 Behavior of the solution as -+ ±oo Additional relations on discontinuities Main result and its discussion A remark on deriving additional relations when condition (7.2.7) is not satisfied ' Hugoniot manifold Behavior of the Hugoniot curve in the vicinity of Jouget points and nonuniqueness of solutions of self-similar problems Nonlinear small-amplitude waves in anisotropic elastic media Basic equations Quasilongitudinal waves Quasitransverse waves Riemann waves Shockwaves 452

6 7.4.6 Self-similar problems and nonuniqueness of solutions Waves in viscoelastic media, vanishing viscosity Role of the wave anisotropy and passage to the limit g -> Final conclusions Electromagnetic shock waves in ferromagnets Long-wave approximation. Elastic analogy Structure of electromagnetic shock waves The set of admissible discontinuities Nonuniqueness of solutions Shock waves in composite materials Basic equations and the discontinuity structure Discontinuity structure; admissible discontinuities Casefc> Casefc< Longitudinal nonlinear waves in elastic rods Large-scale model Model for moderate-scale motions Equations describing the discontinuity structure Admissible discontinuities More precise large-scale model. Nonuniqueness Ionization fronts in a magnetic field Large-scale model Moderate-scale model The set of admissible discontinuities The simplest self-similar problem Variation of the gas velocity across ionization fronts Constructing the solution of the piston problem Discussion 501 Bibliography 503 Index 535 xm

NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS

NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction

More information

A Finite Volume Code for 1D Gas Dynamics

A Finite Volume Code for 1D Gas Dynamics A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler

More 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

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

EqWorld INDEX.

EqWorld INDEX. EqWorld http://eqworld.ipmnet.ru Exact Solutions > Basic Handbooks > A. D. Polyanin and V. F. Zaitsev, Handbook of Nonlinear Partial Differential Equations, Chapman & Hall/CRC, Boca Raton, 2004 INDEX A

More information

High-resolution finite volume methods for hyperbolic PDEs on manifolds

High-resolution finite volume methods for hyperbolic PDEs on manifolds High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution

More information

Approximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD

Approximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD Approximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD Andrea Mignone Collaborators: G. Bodo, M. Ugliano Dipartimento di Fisica Generale, Universita di Torino (Italy)

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

Chapter 1. Introduction

Chapter 1. Introduction Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear

More information

Finite Volume Schemes: an introduction

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

DEVELOPMENT AND APPLICATION OF GENERALIZED MUSTA SCHEMES

DEVELOPMENT AND APPLICATION OF GENERALIZED MUSTA SCHEMES European Conference on Computational Fluid Dynamics ECCOMAS CFD 26 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 26 DEVELOPMENT AND APPLICATION OF GENERALIZED MUSTA SCHEMES V.

More information

The Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics

The Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation

More 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

Projection Dynamics in Godunov-Type Schemes

Projection Dynamics in Godunov-Type Schemes JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and

More 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

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

Recapitulation: Questions on Chaps. 1 and 2 #A

Recapitulation: Questions on Chaps. 1 and 2 #A Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?

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

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

Finite Volume Methods for Hyperbolic Problems

Finite Volume Methods for Hyperbolic Problems Finite Volume Methods for Hyperbolic Problems This book contains an introduction to hyperbolic partial differential equations and a powerful class of numerical methods for approximating their solution,

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

Part 1: Numerical Modeling for Compressible Plasma Flows

Part 1: Numerical Modeling for Compressible Plasma Flows Part 1: Numerical Modeling for Compressible Plasma Flows Dongwook Lee Applied Mathematics & Statistics University of California, Santa Cruz AMS 280C Seminar October 17, 2014 MIRA, BG/Q, Argonne National

More information

State of the Art MHD Methods for Astrophysical Applications p.1/32

State of the Art MHD Methods for Astrophysical Applications p.1/32 State of the Art MHD Methods for Astrophysical Applications Scott C. Noble February 25, 2004 CTA, Physics Dept., UIUC State of the Art MHD Methods for Astrophysical Applications p.1/32 Plan of Attack Is

More information

RESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method

RESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method RESEARCH HIGHLIGHTS (Last update: 3 rd April 2013) Here I briefly describe my contributions to research on numerical methods for hyperbolic balance laws that, in my view, have made an impact in the scientific

More information

A note on the carbuncle in shallow water simulations

A note on the carbuncle in shallow water simulations A note on the carbuncle in shallow water simulations Friedemann Kemm Institute for Applied Mathematics and Scientific Computing, Brandenburgische Technische Universität Cottbus, Platz der Deutschen Einheit

More information

Multi-D MHD and B = 0

Multi-D MHD and B = 0 CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B

More information

Target Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu

Target Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Muon Collider/Neutrino Factory Collaboration Meeting May 26 28, CERN, Geneva U.S. Department of Energy Target Simulations Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Center for Data Intensive

More information

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

R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant

R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant CONTENTS I. Introductory Remarks S1. General Information about the Variety of Solutions.

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 Abstract We study two systems of conservation laws for polymer flooding in secondary

More information

Godunov methods in GANDALF

Godunov methods in GANDALF Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least

More information

A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS

A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor

More information

FDM for wave equations

FDM 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 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 Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws

A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University

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

Heuristical and numerical considerations for the carbuncle phenomenon

Heuristical and numerical considerations for the carbuncle phenomenon Heuristical and numerical considerations for the carbuncle phenomenon arxiv:57.666v2 [physics.comp-ph] 6 Mar 27 Friedemann Kemm March 7, 27 Keywords: Carbuncle phenomenon, High speed flow, Shock instability,

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

APPLIED PARTIAL DIFFERENTIAL EQUATIONS

APPLIED PARTIAL DIFFERENTIAL EQUATIONS APPLIED PARTIAL DIFFERENTIAL EQUATIONS AN I N T R O D U C T I O N ALAN JEFFREY University of Newcastle-upon-Tyne ACADEMIC PRESS An imprint of Elsevier Science Amsterdam Boston London New York Oxford Paris

More information

A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS

A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods

More 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

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

Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws

Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang

More 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 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

Partial Differential Equations

Partial Differential Equations Partial Differential Equations Analytical Solution Techniques J. Kevorkian University of Washington Wadsworth & Brooks/Cole Advanced Books & Software Pacific Grove, California C H A P T E R 1 The Diffusion

More information

The Hopf equation. The Hopf equation A toy model of fluid mechanics

The Hopf equation. The Hopf equation A toy model of fluid mechanics The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van

More information

Advection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II

Advection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises

More information

arxiv: v2 [math.ap] 1 Jul 2011

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

A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws

A 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

NUMERICAL METHODS IN ASTROPHYSICS An Introduction

NUMERICAL METHODS IN ASTROPHYSICS An Introduction -1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,

More information

Introduction to Magnetism and Magnetic Materials

Introduction to Magnetism and Magnetic Materials Introduction to Magnetism and Magnetic Materials Second edition David Jiles Ames Laboratory, US Department of Energy Department of Materials Science and Engineering and Department of Electrical and Computer

More information

Introduction to Magnetohydrodynamics (MHD)

Introduction to Magnetohydrodynamics (MHD) Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality

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

Shock-capturing Schemes for a Collisionless Two-fluid Plasma Model

Shock-capturing Schemes for a Collisionless Two-fluid Plasma Model Shock-capturing Schemes for a Collisionless Two-fluid Plasma Model E. Alec Johnson Department of Mathematics, UW-Madison Presented on August 29, 2 at Sandia National Laboratories c 23, Kenneth Lang, Tufts

More information

Nonlinear Problems of Elasticity

Nonlinear Problems of Elasticity Stuart S. Antman Nonlinear Problems of Elasticity With 105 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo Hong Kong Barcelona Budapest Contents Preface vn Chapter I. Background

More information

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

A Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes )

A Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes ) A Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes ) Takahiro MIYOSHI and Kanya KUSANO 1) Hiroshima University, Higashi-Hiroshima 739-856, Japan 1) Nagoya University,

More information

High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation

High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with

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

Comparison of Approximate Riemann Solvers

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

Various Hydro Solvers in FLASH3

Various Hydro Solvers in FLASH3 The Center for Astrophysical Thermonuclear Flashes Various Hydro Solvers in FLASH3 Dongwook Lee FLASH3 Tutorial June 22-23, 2009 An Advanced Simulation and Computing (ASC) Academic Strategic Alliances

More information

Propagation of discontinuities in solutions of First Order Partial Differential Equations

Propagation of discontinuities in solutions of First Order Partial Differential Equations Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in

More information

Workshop on Compressible Navier-Stokes Systems and Related Problems (I) March 5-10, 2018 TITLE & ABSTRACT

Workshop on Compressible Navier-Stokes Systems and Related Problems (I) March 5-10, 2018 TITLE & ABSTRACT Workshop on Compressible Navier-Stokes Systems and Related Problems (I) March 5-10, 2018 TITLE & ABSTRACT (Last updated: 6 March 2018) Classification of asymptotic states for radially symmetric solutions

More information

METHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS

METHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS METHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS V.I. Agoshkov, P.B. Dubovski, V.P. Shutyaev CAMBRIDGE INTERNATIONAL SCIENCE PUBLISHING Contents PREFACE 1. MAIN PROBLEMS OF MATHEMATICAL PHYSICS 1 Main

More information

c 1999 Society for Industrial and Applied Mathematics

c 1999 Society for Industrial and Applied Mathematics SIAM J. SCI. COMPUT. Vol. 21, No. 3, pp. 1115 1145 c 1999 Society for Industrial and Applied Mathematics A SIMPLE METHOD FOR COMPRESSIBLE MULTIFLUID FLOWS RICHARD SAUREL AND RÉMI ABGRALL Abstract. A simple

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

On the entropy stability of Roe-type finite volume methods

On the entropy stability of Roe-type finite volume methods Proceedings of Symposia in Applied Mathematics On the entropy stability of Roe-type finite volume methods Mária Lukáčová - Medvid ová and Eitan Tadmor Abstract. We study the entropy stability of a class

More information

Extremum-Preserving Limiters for MUSCL and PPM

Extremum-Preserving Limiters for MUSCL and PPM arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA

More information

A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws

A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws S. Jaisankar, S. V. Raghurama Rao CFD Centre, Department of Aerospace Engineering, Indian Institute of Science, Bangalore, India-56 Abstract

More information

Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2

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

Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column

Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of

More information

Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws

Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang

More information

Generalized Functions Theory and Technique Second Edition

Generalized Functions Theory and Technique Second Edition Ram P. Kanwal Generalized Functions Theory and Technique Second Edition Birkhauser Boston Basel Berlin Contents Preface to the Second Edition x Chapter 1. The Dirac Delta Function and Delta Sequences 1

More information

THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK

THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK G.P. Zank, I. Kryukov, N. Pogorelov, S. Borovikov, Dastgeer Shaikh, and X. Ao CSPAR, University of Alabama in Huntsville Heliospheric observations

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

A simple robust and accurate a posteriori subcell finite

A simple robust and accurate a posteriori subcell finite A simple robust and accurate a posteriori subcell finite volume limiter for the discontinuous Galerkin method Michael Dumbser ERC Grant STiMulUs FET-HPC Grant ExaHyPE 1/ 69 Séminaire LJLL, Paris, France,

More information

ELECTRODYNAMICS OF CONTINUOUS MEDIA

ELECTRODYNAMICS OF CONTINUOUS MEDIA ELECTRODYNAMICS OF CONTINUOUS MEDIA by L. D. LANDAU and E. M. LIFSHITZ Institute of Physical Problems, USSR Academy of Sciences Volume 8 of Course of Theoretical Physics Translated from the Russian by

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

Numerical Hydrodynamics in Special Relativity

Numerical Hydrodynamics in Special Relativity Numerical Hydrodynamics in Special Relativity J. M. Martí Departamento de Astronomía y Astrofísica Universidad de Valencia 46100 Burjassot (Valencia), Spain martij@godunov.daa.uv.es E. Müller Max-Planck-Institut

More information

Chapter 4: Non-Linear Conservation Laws; the Scalar Case. The conservation law in eqn. (4.1) will be hyperbolic if its eigenvalue ( ) ( )

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

Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models

Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models 0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,

More information

Comparison of (Some) Algorithms for Edge Gyrokinetics

Comparison of (Some) Algorithms for Edge Gyrokinetics Comparison of (Some) Algorithms for Edge Gyrokinetics Greg (G.W.) Hammett & Luc (J. L.) Peterson (PPPL) Gyrokinetic Turbulence Workshop, Wolfgang Pauli Institute, 15-19 Sep. 2008 w3.pppl.gov/~hammett Acknowledgments:

More 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

An efficient shock-capturing central-type scheme for multidimensional relativistic flows. II. Magnetohydrodynamics

An efficient shock-capturing central-type scheme for multidimensional relativistic flows. II. Magnetohydrodynamics A&A 400, 397 413 (2003) DOI: 10.1051/0004-6361:20021641 c ESO 2003 Astronomy & Astrophysics An efficient shock-capturing central-type scheme for multidimensional relativistic flows II. Magnetohydrodynamics

More information

The Physics of Fluids and Plasmas

The Physics of Fluids and Plasmas The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the

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

Sung-Ik Sohn and Jun Yong Shin

Sung-Ik Sohn and Jun Yong Shin Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic

More information

x a(x) H(U), , H(U) =

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

Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu

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

c 2003 Society for Industrial and Applied Mathematics

c 2003 Society for Industrial and Applied Mathematics SIAM J. SCI. COMPUT. Vol., No., pp. 57 7 c 00 Society for Industrial and Applied Mathematics ON THE ARTIFICIAL COMPRESSION METHOD FOR SECOND-ORDER NONOSCILLATORY CENTRAL DIFFERENCE SCHEMES FOR SYSTEMS

More information

Numerical Solution to Strong Cylindrical Shock Wave in the Presence of Magnetic Field

Numerical Solution to Strong Cylindrical Shock Wave in the Presence of Magnetic Field Numerical Solution to Strong Cylindrical Shock Wave in the Presence of Magnetic Field Narsimhulu Dunna, Addepalli Ramu and Dipak Kumar Satpathi Abstract One-dimensional flow of converging shock waves in

More information

VISCOUS FLUX LIMITERS

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

Shock Reflection-Diffraction, Nonlinear Conservation Laws of Mixed Type, and von Neumann s Conjectures 1

Shock Reflection-Diffraction, Nonlinear Conservation Laws of Mixed Type, and von Neumann s Conjectures 1 Contents Preface xi I Shock Reflection-Diffraction, Nonlinear Conservation Laws of Mixed Type, and von Neumann s Conjectures 1 1 Shock Reflection-Diffraction, Nonlinear Partial Differential Equations of

More information

GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS

GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,

More information

A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws

A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace

More 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

A 5-wave relaxation solver for the shallow water MHD system

A 5-wave relaxation solver for the shallow water MHD system A 5-wave relaxation solver for the shallow water MHD system François Bouchut, Xavier Lhébrard Abstract The shallow water magnetohydrodynamic system describes the thin layer evolution of the solar tachocline.

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

arxiv:astro-ph/ v1 17 May 2002

arxiv:astro-ph/ v1 17 May 2002 Astronomy & Astrophysics manuscript no. (will be inserted by hand later) An efficient shock-capturing central-type scheme for multidimensional relativistic flows arxiv:astro-ph/0205290v1 17 May 2002 I.

More information