COMPUTATIONAL METHODS AND ALGORITHMS Vol II - Computational Methods for Compressible Flow Problems - Remi Abgrall
|
|
- Noah Norton
- 5 years ago
- Views:
Transcription
1 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, numerical examples Contents 1. Introduction. A Brief Description of the Solutions.1 On the Notion of Solutions. Hyperbolicity: Characteristic Form of the Equations 3. Numerical Schemes for 1-D Problems 3.1 The Principle of Conservative Schemes 3. The Godunov Scheme 3.3 The Roe Scheme 3.4 Kinetic Schemes 3.5 High Order Accuracy Scheme 4. Schemes for Multidimensional Problems 5. Numerical Examples 5.1 A One-dimensional Example: The Sod Problem 5. A Two-dimensional Example 6. Conclusions Bibliography Summary In this chapter, we provide an introduction to the numerical schemes devoted to the simulation of compressible flow problems. We first review a few elements of modeling of compressible flow problems and then provide a brief description of the solutions to one dimensional problems. We identify a class of partial differential equations (PDEs) which share common features with those describing flow problems. In particular, we show that it is not possible to expect globally smooth solutions. Then we describe a class of numerical schemes that are able to approximate the solutions, first for onedimensional problems and for the multidimensional problems. Finally, numerical examples are provided. 1. Introduction In many industrial applications, fluid flow cannot be approximated by the incompressible Navier Stokes equations: either the speed is too high (e.g. flows around aircraft or reentry vehicles), or the density is not uniform (e.g. in car engines), or acoustic effects are not negligible (e.g. when an aircraft is landing). In these cases and
2 many others, one has to rely on the Navier Stokes equations for compressible flow., For simplicity, in this chapter we focus on flows where viscous effects can be neglected. The state of flow is represented by a vector of conserved variables W = ( ρρu,, E) T where ρ is the density, u is the flow velocity, and 1 E = ρε + ρu is the total energy. Here, ε represents the specific internal energy. in the light of the basic principles of continuum mechanics, W satisfies the following: q div Fq ( ) 0 t + = for x Ω and t > 0 q( x,0) = q ( x ) for x Ω (1) 0 In (1), Ω is an open subset of ( E+ p) u N, N = 1,,3, the flux F is given by Fq ( ) = ( + ρ u ρu u pid) () The system is closed by a relation between the pressure p and the state variables. In this paper, we consider the simplest example of an ideal gas for which 1 p ( 1) = γ E ρu. (3) In (3), γ is the ratio of specific heats and is γ =1.4 in all applications we present in the paper The first relation given by (1), namely, ρ div( ) 0 t + ρu = (4) signifies conservation of mass. The second relation, ρ ( u) + div( ρu u ) + div p = 0, t is for conservation of momentum and application of the fundamental principles of mechanics. The nest relation is for conservation of total energy, E div(( E p) ) 0, t + + u = and represents the first law of thermodynamics.
3 There is yet another relation, namely the second law of thermodynamics. It states that the entropy of a fluid particle increases with time, that is s s 0. t + u. (5) By multiplying (5) by ρ and (4) by s, and then summing up, we get ρs + div( ρ su ) 0. (6) t In (5) and (6), s represents specific entropy. In the case of an ideal gas, we have γ s= s0 + cυ log( pρ ). (7) One of its important properties is that the function 1 s = s( ρ, E ρu ) is a concave function of 1 ( ρ, E ρu ). The problem of relation (5) is that it is written in nonconservative form. In particular, as we see later, it is natural to admit solution for which neither u nor s are continuous, so that the product u. s may be undefined in general. Interest in relation (7) is two fold. First, it is written in conservative form so that, as we see later, the problem of non-smooth solution can be solved. Second, the mathematical entropy, S = ρs( ρ, u, E) is a convex function of its arguments.. A Brief Description of the Solutions It is very natural to consider non continuous solutions. Take for example the situation of a shock tube. This is an experimental device which is used to generate shock waves. It is a tube filled with gas. It is initially made of two sections separated by a diaphragm, see Figure 1. The left and right chambers are filled with gas with different velocity and pressure. In this experimental setup, the initial velocity is set to zero; the 1D-Riemann problem, Figure 1: the shock tube problem q F + = 0 x, t > 0 t x ql if x < 0 q( x,0) = qr if x > 0 (8)
4 plays a very important role in the numerical simulation of compressible flows. In the Riemann problem, the initial condition is discontinuous. An analysis of the solution of the shock tube problem shows that it is equivalent to the superposition of two piston problems, see Figure. Figure : The piston problems In one case, the piston is moving to the left, in the other case, it is moving to the right. If one considers an experiment with a tube filled with gas where a piston, at time t = 0 moves with a velocity V to the left, then the density and the pressure in the vicinity of the piston will decrease, and a wave is created from the piston. If we sit far enough from the piston, we do not feel the effect of changing the pressure: in general, this wave connects the area where the solution is still uniform and the fan. The connection is continuous but not continuously differentiable. If one now pushes the piston to the right, one can show that the state flow parameters are uniform in the neighborhood of the piston, also uniform far enough, and the two uniform states are connected by a shock wave, that is a moving discontinuity where all the flow parameters are discontinuous. If the piston is moved smoothly, depending on the conditions, the flow may be continuous at any time or discontinuities can be generated. In the case where the solution is not smooth, what is the meaning of the system (1)? The answer has been provided by P.Lax..1 On the Notion of Solutions The idea is to integrate (1) against compactly supported smooth test functions 1 N + ϕ C c ( Ω [0, T ]). Assuming that q is smooth, if we integrate (1) multiplied by ϕ, we get
5 Ω[0,Τ] ϕ q + F( q) ϕ dxdt q(,0). (,0) d = 0 t x ϕ x x Ω (9) The advantage of (9) w.r.t. (1) is that q does not appear anymore with derivatives: they have been transferred to the terms containing ϕ. Thanks to this, we can define the weak 1 solutions of (1) as the functions in L ( Ω [0, T]) L ( Ω [0, T]) such that for any 1 ( [0, ]) N + ϕ Cc Ω T, (9) holds. Similarly, the inequality (7) has to be understood as the solutions q such that for any ϕ C 1 ( Ω [0, T]), ϕ 0, one has c Ω[0,Τ] ϕ ρs + ρsu. ϕ dxdt s(,0) (,0) d 0. t x ϕ x x Ω This relation shows the advantage of (5) with respect to (6). Solutions that satisfies the inequalities (6) for all the convex entropy U associated to an entropy flux G, that is G = U. F and Ω[0,Τ] q W ϕ U + G. ϕ dxdt U(,0) (,0) d 0. t x ϕ x x Ω are called entropy solutions. For a characterization of the couple entropy-flux for the Euler equations, ( UG,, ) see [15]. The solutions we are interested in are piecewise smooth. The discontinuities of the solutions lie on a piecewise smooth curve, see Figure 3. By using test functions whose supports are localized in the neighborhood of, we get the Rankine-Hugoniot relations (10) [ F n ( q)] = σ[ q] (11) + where [ f ] = (.) (.) is the jump of f at M, F speed of the discontinuity. The inequality (7) becomes n is the normal flux and σ is the [ ρu.ns ] σ[ s] (1) In (1), the definition of [ ] f is important because of the inequality. Equation (1) indicates that the specific entropy of a fluid particle increases across a discontinuity.
6 Figure 3: Discontinuous solutions on a piecewise smooth curve Tadmor [3] showed that the solution (if it is bounded) satisfies the following minimum principle s( x, t) min s( y,0) y x t u where x is the Euclidean norm of x and. Hyperbolicity: Characteristic Form of the Equations (13) u is the L norm of the velocity field A straightforward calculation shows that if nn (,n ) is any vector, then the matrix n xa + n yb, where A and B are the Jacobian matrices of the flux, is diagonalizable in. The system (1) is said to be hyperbolic. More precisely, set c the speed of sound p p c : = = γ. ρ ρ s x The eigenvalues of n xa + n yb are λ 0 = u.n, λ ± = u.n± c. The eigenvalue λ 0 is double (triple in 3D). If t = ( ny, n ), the eigenvectors are respectively x y
7 1 0 1 u c r0 r -ny 0 r u ± n n t x = υ = ±, n υ ± n c + x y u υ. t H ± n u.n n we skip the left eigenvectors 0, 0, ±, they can be found in [13]. t It is possible to find special solutions to (1) of the form q= q( σ), σ such that dq rq ( ( )) d σ = σ where r is any of the eigenvectors of n A + n B and (. σ = σ x ω, t ), ω being a fixed x direction. We are looking for planar waves. Using (1), we see that ( χ, t) σ( χ, t) satisfies an equation of the type σ ( q(, t)) 0 t + σ λ σω, χ =. These three equations are non linear scalar equations. These particular solutions are characterized as follows 1. waves associated with λ = u ± c, these are fans where pρ γ and u c remain γ 1 constant. waves associated with λ = uu, and p remain constant. The functions pρ γ, u c, u and p are the Riemann invariants. γ 1 y When analyzing the structure of the Rankine Hugoniot relations, it can be easily seen that the piecewise constant solutions must be of the form 1 either u and p remain constant, ρ and s can be arbitrary: these are contact discontinuities; or all the variables are discontinuous, they are shock waves. Let us denote with a superscript + the post-shock conditions and with a superscript the pre-shock conditions and denoting by σ is the shock speed and υ = u σ, and M = υ the c relative match number of the shock, we have
8 ρ ρ + 1 γ 1 = + γ+ 1 M γ+ 1 υ υ p p γ 1 = + γ+ 1 M γ 1 γ γ 1 = M γ+ 1 γ+ 1 From this we see that the density and the pressure increase through a shock. The opposite is true for the velocity. We observe that u and p are invariant parameters for the wave associated with This remark is useful in the effective solution of the Riemann problem. λ = u. These particular solutions are building blocks of the solution of the Riemann problem 8. They can be observed in Figures 7, 8, and 9. More on the structure of the solution can be found in [13, 8] Bibliography TO ACCESS ALL THE 36 PAGES OF THIS CHAPTER, Visit: [1] R. Abgrall.( 1988) Generalisation of the Roe scheme for the computation of mixture of perfect gases. Rech. Aérospatiale (English edition), : [] R. Abgrall.(1994) On Essentially Non Oscillatory Schemes on Unstructured Meshes: Analysis and Implementation. J. Comp. Phys., 114(1): [3] R. Abgrall and S. Karni.() Compressible multifluid flows. J. Comput. Physics, submitted. [4] R. Abgrall, S. Lantéri, and Th. Sonar. (1996) Eno schemes for compressible fluid dynamics. ZAMM, à paraitre. Rapport de l Institut fur Mathematik der Universitaet Hamburg, num 116. [5] J.P. Boris and D.L. Book.(1973) Flux-corrected transport i: SHASTA a fluid transport algorithm that works. J. Comput. Phys, 11: [6] P. Cargo and G.Gallice. (1996) Construction of a roe matrix for the ideal magnetohydrodynamics equations. Z. Angew. Math. Mech., 76:
9 [7] P. Colella, H.M. Glaz, and R.E. Ferguson. (1989) Multifluid algorithms for Eulerian finite difference methods. unpublished. [8] F. Coquel and P.G. LeFloch. (1996) An entropy satisfying muscl scheme for systems of conservation laws. Numer. Math., 74(1):1-33. [9] H. Deconinck, P.L. Roe, and R. Struijs. (1993) A multidimensional generalisation of Roe s differnce splitter for the Euler equations. Computer and Fluids :15-. [10] H. Deconinck, R. Struijs, G. Bourgeois, and P.L. Roe. (1993) Compact advection schemes on unstructured meshes. VKI Lecture Series , Computational Fluid Dynamics. [11] B. Dubroca, Communication privée. [1] E. Godlewski and P.A. Raviart. (1991) Hyperbolic systems of conservation laws. Mathématiques et Applications. Ellipses, Paris. [13] E. Godlewski and P.A. Raviart. (1996) Numerical approximation of systems of conservation laws, volume 118 of Applied Mathematical Sciences, Springer. [14] A. Harten, S. Osher, B. Engquist, and S.R. Chakravarthy. (1987) Uniformly high order accurate nonoscilaltory schemes iii. J. Comput. Phys., 71: [15] Amiram Harten. (1983) On the symmetric form of systems of conservation laws with entropy. J. Comput. Phys.. [16] P. Jenny, B. Mueller, and H. Thomann. (1997) Correction of conservative euler solvers for gas mixtures. J. Comp. Phys, 13: [17] S. Karni. (1996) Multi-component flow calculations by a consistent primitive algorithm. J. Comp. Phys., 11: [18] S. Karni. (1998) (E.F. Toro and J.F. Clarke eds.). A level-set scheme for compressible interfaces. Pages [19] B. Khobalatte and B. Perthame. (1994) A maximum principle on the entropy and minimal limitations for kinetic schemes. Math. Comp., 6: [0] P. Lax and B. Wendroff. (1960) Systems of conservation laws. Comm. Pure Appl. Math., 13: [1] W.F. Noh and P. Woodward. (1976) Slic (simple line interface calculation). Lawrence Livermore National Laboratory Rep. UCRL [] B. Perthame. (1990) Boltzman type equations for gas dynamics and the entropy property. SIAM J. Numer. Anal., 7(6): [3] B. Perthame. (1999) An introduction to kinetic schemes for gas dynamics, volume 5 of Lecture Notes in Computational Sciences and Engineering, chapter I. Springer Verlag. [4] K.G. Powell. (1994) An approximate Riemann solvers for magnetohydrodynamics (that works in more than one dimension). Technical Report 94-4, ICASE. [5] R. Abgrall. () Toward the ultimate conservative scheme: Following the quest. J. Comput. Phys, submitted.
10 [6] P.L. Roe. (1981) Approximate Riemann solvers, parameter vectors and difference schemes. J. Comp. Phys., 43: [7] R. Saurel and R. Abgrall. (1999) Some models and methods for compressible multifluid and multiphase flows. J. Comp. Phys., 150(): [8] D. Serre. (1996) Systéme de lois de conservation, volume 1 et. Diderot. [9] C.W.Shu and S. Osher. (1988) Efficient implementation of essentially non oscillatory shock capturing schemes. J. Comput. Phys., 77: [30] R. Struijs, H. Deconinck, and P.L. Roe. (1991) Fluctuation Splitting Schemes for the D Euler equations. VKI LS , Computational Fluid dynamics. [31] P.K. Sweby. (1984) High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM J. Numer. Anal., 1: [3] E. Tadmor. (1987) The numerical viscosity of entropy stable schemes for systems of conservation laws. Math. Comp., 49: [33] B.van Leer. (1974) Towards the ultimate conservative scheme. ii. monotonicity and conservation combined in a second order scheme. J. Comput. Phys., 14: [34] B. van Leer. (1979) Towards the ultimate conservative scheme. v.a second order sequel of Godunov s scheme. J. Comput. Phys., 3: [35] P.R. Woodward and P. Collela. (1984) The numerical simulation of two dimensional flows with strong shocks. J. Comput Phys., 54: [36] P.R. Woodward and P. Collela. (1984) The piece-wise parabolic method (ppm) for gas dynamic simulations. J. Comput. Phys., 54: [37] H.C. Yee and A. Harten. (1987) Implicit tvd schemes for hyperbolic conservation laws in curvilinear coordinates. AIAA Journal, 5():66-74.
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 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 informationGas 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 informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationProjection 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 informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
More informationPOSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS. Cheng Wang. Jian-Guo Liu
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 3 Number May003 pp. 0 8 POSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS
More informationFUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC PROBLEMS WITH DISCONTINUITIES
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK FUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC
More informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
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 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 informationAn Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term
An Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term Alexander Kurganov Department of Mathematics, Tulane University, 683 Saint Charles Avenue, New Orleans, LA 78,
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More 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 informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationA 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 informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More 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 informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationAmiram HARTEN LIST OF PUBLICATIONS
1 Amiram HARTEN LIST OF PUBLICATIONS ARTICLES 1. S. Cuperman and A. Harten The Solution of One-Fluid Equations with Modified Thermal Conductivity for the Solar Wind Cosmic-Electrodynamics I, 205-217 (1970).
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 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 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 informationc 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 informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationTwo-Dimensional Riemann Solver for Euler Equations of Gas Dynamics
Journal of Computational Physics 167, 177 195 (2001) doi:10.1006/jcph.2000.6666, available online at http://www.idealibrary.com on Two-Dimensional Riemann Solver for Euler Equations of Gas Dynamics M.
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationReal Gas Computation Using an Energy Relaxation Method and High-order WENO Schemes
NASA/CR-998-87 ICASE Report No. 98-4 Real Gas Computation Using an Energy Relaation Method and High-order WENO Schemes Philippe Montarnal and Chi-Wang Shu Brown University, Providence, Rhode Island Institute
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
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 informationApplications 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 informationANALYSIS AND IMPLEMENTATION OF THE GAS-KINETIC BGK SCHEME FOR COMPUTATIONAL GAS DYNAMICS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 25, 21±49 (1997) ANALYSIS AND IMPLEMENTATION OF THE GAS-KINETIC BGK SCHEME FOR COMPUTATIONAL GAS DYNAMICS CHONGAM KIM*, KUN XU, LUIGI MARTINELLI
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 informationOn 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 informationRESEARCH 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 informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationLarge Time Step Scheme Behaviour with Different Entropy Fix
Proceedings of the Paistan Academy of Sciences: A. Physical and Computational Sciences 53 (): 13 4 (016) Copyright Paistan Academy of Sciences ISSN: 0377-969 (print), 306-1448 (online) Paistan Academy
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More 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 informationResearch Article A New Flux Splitting Scheme Based on Toro-Vazquez and HLL Schemes for the Euler Equations
Computational Methods in Physics Volume 4, Article ID 8734, 3 pages http://dx.doi.org/.55/4/8734 Research Article A New Flux Splitting Scheme Based on Toro-Vazquez and HLL Schemes for the Euler Equations
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 informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
More informationJournal of Computational Physics
Journal of Computational Physics 9 () 759 763 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Short Note A comment on the computation
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 informationArt Checklist. Journal Code: Article No: 6959
Art Checklist Journal Code: JCPH Article No: 6959 Disk Recd Disk Disk Usable Art # Y/N Format Y/N Remarks Fig. 1 Y PS Y Fig. 2 Y PS Y Fig. 3 Y PS Y Fig. 4 Y PS Y Fig. 5 Y PS Y Fig. 6 Y PS Y Fig. 7 Y PS
More informationComputations of Compressible Multifluids
Journal of Computational Physics 69, 594 623 (200) doi:0.006/jcph.2000.6685, available online at http://www.idealibrary.com on Computations of Compressible Multifluids Rémi Abgrall and Smadar Karni, Mathématiques
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 informationThe 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 informationOn the Modeling and Simulation of a Laser-Induced Caviation Bubble
1 On the Modeling and Simulation of a Laser-Induced Caviation Bubble Ali Zein Supervisor: Prof. Dr. Gerald Warnecke Co-supervisor: Dr. Maren Hantke Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More information2 Introduction Standard conservative schemes will admit non-physical oscillations near material interfaces, which can be avoided by the application of
A Conservative Approach to the Multiphase Euler Equations without Spurious Pressure Oscillations Xu-Dong Liu Ronald P. Fedkiw Stanley Osher June 8, 998 Abstract In [4], the pressure evolution equation
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
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 informationDepartment of Mathematics University of California Santa Barbara, Santa Barbara, California,
The Ghost Fluid Method for Viscous Flows 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 9435 Email:fedkiw@cs.stanford.edu Xu-Dong Liu Department of Mathematics
More informationHyperbolic Problems: Theory, Numerics, Applications
AIMS on Applied Mathematics Vol.8 Hyperbolic Problems: Theory, Numerics, Applications Fabio Ancona Alberto Bressan Pierangelo Marcati Andrea Marson Editors American Institute of Mathematical Sciences AIMS
More informationHigh 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 informationGas Evolution Dynamics in Godunov-type Schemes and Analysis of Numerical Shock Instability
NASA/CR-1999-208985 ICASE Report No. 99-6 Gas Evolution Dynamics in Godunov-type Schemes and Analysis of Numerical Shock Instability Kun Xu ICASE, Hampton, Virginia Institute for Computer Applications
More informationExtremum-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 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 informationOn 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 informationThe Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics
The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationA Simple Model for Cavitation with Non-condensable Gases
A Simple Model for Cavitation with Non-condensable Gases Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène Mathis To cite this version: Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène
More informationAn interface treating technique for compressible multi-medium flow with. Runge-Kutta discontinuous Galerkin method. Abstract
An interface treating technique for compressible multi-medium flow with Runge-Kutta discontinuous Galerkin method Chunwu Wang 1 and Chi-Wang Shu 2 Abstract The high-order accurate Runge-Kutta discontinuous
More informationA 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 informationComputational 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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationMeasure-valued - strong uniqueness for hyperbolic systems
Measure-valued - strong uniqueness for hyperbolic systems joint work with Tomasz Debiec, Eduard Feireisl, Ondřej Kreml, Agnieszka Świerczewska-Gwiazda and Emil Wiedemann Institute of Mathematics Polish
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationNumerische Mathematik
Numer. Math. ) :545 563 DOI.7/s--443-7 Numerische Mathematik A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang Chi-Wang Shu Received: 7 July / Revised: 5 September
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationCausal 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 informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationFinite 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 informationSolution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers
Solution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers Alexander Kurganov, 1, * Eitan Tadmor 2 1 Department of Mathematics, University of Michigan, Ann Arbor, Michigan
More informationMATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS
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
More informationSecond order entropy diminishing scheme for the euler equations
Second order entropy diminishing scheme for the euler equations Frédéric Coquel, Philippe Helluy, Jacques Schneider To cite this version: Frédéric Coquel, Philippe Helluy, Jacques Schneider. Second order
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
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 informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationA 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 informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
More informationFlux corrected transport solver for solving 1D multiphase equations for drilling applications
Computational Methods in Multiphase Flow VII 169 Flu corrected transport solver for solving 1D multiphase equations for drilling applications P. A. Varadarajan & P. S. Hammond Schlumberger Gould Research,
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 informationLinear high-resolution schemes for hyperbolic conservation laws: TVB numerical evidence 1
Linear high-resolution schemes for hyperbolic conservation laws: TVB numerical evidence C. Bona a,2, C. Bona-Casas a,3, J. Terradas b,4 arxiv:885v2 [gr-qc] 2 Dec 28 a Departament de Física. Universitat
More information