International Engineering Research Journal
|
|
- Cory Stevens
- 5 years ago
- Views:
Transcription
1 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 Pune University, Matoshri college of Engineering, Nashik, India. Mechanical Department, SavitribaiPhule Pune University, Matoshri college of Engineering, Nashik, India. ϯ Mechanical Department, SavitribaiPhule Pune University, Sandip Institute of Engineering and Management, Nashik, India. Abstract In this work high resolution schemes are developed to solve gas dynamics problems. Two dimensional Euler equations are used as governing equations. One of the method to solve the Euler equations that are hyperbolic in nature, is to treat them as initial value problem. The Riemann solvers are adopted for solving the initial value problems. Riemann solvers are mainly of two types, exact and approximate solvers. The Godunov solver provides an exact solution whereas HLL and HLLC an approximate solution. These methods are developed to solve gas dynamic problems. In this work Godunov solver, HLLC solvers are considered for obtaining solutions to several D Sod tests. First, these methods are benchmarked with the solutions to exact Riemann problems. Later, a qualitative comparison is presented for understanding the performance of these methods. Keywords:High Resolution Methods, Riemann Solvers, Exact Riemann Solvers, Approximate Riemann Solvers, HLLC.. Introduction In this work high resolution schemes are developed to solve gas dynamics problems. One of the method to solve the Euler equations that are hyperbolic in nature, is to treat them as initial value problem. Riemann problem offers a solution strategy to solve a system of hyperbolic conservation equations. Following are the equations that are considered for discretization, t ρ ρu ρv ρe + x ρu ρu + p ρuv ρu(e + p/ρ + y ρv ρuv ρv + p ρv(e + p/ρ = () In case of Euler equation the Riemann contain so called shock tube problem, the basic problem of gas dynamics. The Riemann problem solution of two types Exact and Approximate. The Godunov's method provides exact solution whereas, HLL (Harten, Lax and Van Leer) approach is an approximation for inter-cell numerical flux. It offers a very efficient and robust approximation to Godunov type methods. The shortcomings of HLL approach are overcome in a modified approach called the HLLC (C stands for contact) Riemann solver. First, all the three Riemann solvers the Godunov, HLL and HLLC are benchmarked with the exact solutions of the five Riemann problems, the Sod tests. The conditions in these tests are too hypothetical to exist in a real gas dynamic problem. An agreement to the results of these problems ensures robustness of the solver. Later these approximate Riemann solvers are used for solving another five -D tests In 988 S. F. Davis have done a number of algorithms for obtaining the bounds on smallest and the largest signal velocities in the solution to Riemann problem where bounds are used to construct the approximate Riemann solution of Harten, Lax, and van Leer which contains only one intermediate state. This approximate Riemann solution is used to construct first- and second-order Godunov type methods. The methods described in this paper require only the characteristic velocities of a hyperbolic system. Therefore, these methods are easily applied to any system of conservations laws. He has applied the method of this paper to the steady supersonic Euler equations, the unsteady Euler equations with a nonconvex equation of state on the nonlinear shallow water equations. In 988 B. Enfield.has describe a new approximate Riemann solver for compressible gas flow. In contrast to previous Riemann solvers, where a numerical approximation for the pressure and the velocity at the contact discontinuity is computed, he derived a numerical approximation for the largest and smallest signal velocity in the Riemann problem. Having obtained the numerical signal velocities, he used theoretical results by Harten, Lax and van Leer to obtain the full approximation. A stability condition for the numerical signal velocities is derived. He has also demonstrate a relation between the signal velocities and the dissipation contained in the corresponding Godunov-type method. Also given numerical results for the one- and two-dimensional
2 Special Edition PGCON-MECH-7 compressible gas dynamics equations. In 3 Liska, R. and Wendroff, B. presented results of computations with eight explicit finite difference schemes on a suite of one-dimensional and two-dimensional test problems for the Euler equations. Both dimensionally split and two-dimensional schemes are represented, as are central and upwind-biased methods, and all are at least second-order. For D tests they had chosen five D (in x) Riemann problems. In 99 E. F. Toro investigated Riemann solvers and numerical methods for fluid dynamics. He considered more difficult set of D problems has been who describes in detail several popular methods and shows their behavior on his tests, all of which have easily computed exact solutions. In 994 E. F. Toro, M. Spruce and W. Spears presented new ways of obtaining wave speed estimates also implemented improved Riemann solver in the second order WAF method and tested for D problem.. Methodology. Understand shock tube problem and Initial value problems.. Use of Euler equations to solve IVPs and discretize these equations. 3.Obtain Exact solutions of IVPs (Five Sod tests), solving Riemann problems. 4. Use of these Riemann problems to obtain exact solutions, develop the Godunov's solver.. Benchmarking with Sod test D. 6. Develop approximate solution strategy, using HLLC approach for D. 7. Benchmarking with Sod test D. 8. Development of high resolution method for solving D U= [ρ, ρu, ρv, ρw,e] T, F(U)=[ ρu, ρu +p, ρuv, ρuw,u(e+p)] T, G(U)=[ ρv, ρuv, ρv +p. ρvw,v(e+p)] T, H(U)=[ ρw, ρuw, ρvw, ρw +p,w(e+p)] T. Here U = is a vector of conserved variables, F(U), G(U), and H(U) are flux vectors. ρ is density of fluid, u, v, w are velocity components in x, y, z directions. p is pressure and E is total energy per unit volume of the fluid. For one-dimensional case, Eqn.() reduces to U t+f(u) x= (3) Here, U = [ρ, ρu, E] T, F(U)=[ ρu, ρu +p,u(e+p)] T. Actually U, F, G, H are vectors but they are not be highlighted henceforth. 3. Basic Riemann Problem For one-dimensional non stationary flows, the state U of a perfect gas is specified completely by three independent variables ρ, u, p and two constants γ, R for a particular gas. From these variables and constants, all of the other state properties may be obtained, and the state variable U may completely be defined over a domain. The Riemann problem for one-dimensional time dependent Euler equations is represented as below: If the vectors of conserved variables U and flux vector F(U) have their usual form as in (3), solve U t+f(u) x= (3) Subjected to U x, = UL if x < ; UR if x <. Here, subscripts L and R represent the two discrete initial left and right states respectively.. Governing Equations The fluid, through which the wave propagating, is compressible and inviscid,. The fluid obeys the ideal gas equation,. No external forces act on the fluid, Thus reducing the Navier-Stokes equation to Euler equations which are written in the Cartesian vector form as, U +F(U) +G(U) Y+H(U) Z= () Where, Fig.Schematic representation of General Riemann problem To simplify the solution procedure, a vector of primitive variables W is often used instead of vector U: W = [ρ,u,p] T (4)
3 Special Edition PGCON-MECH-7 As shown in Fig.., the domain of interest is the (x-t) plane with x ranging from - to.as time progresses the discontinuities between the two initial states are break into leftward and rightward moving waves which are separated by the contact discontinuity. Each wave may either be a shock wave or a rarefaction wave (expansion wave) depending on the initial data, giving rise to four possible wave patterns as shown in Fig.. The region in between the left and right waves is called as the star region and the flow properties in this region are known as star properties.. On both side of contact discontinuity same gas exists.. The gas in the problem under consideration is a perfect gas. Properties of fluid flow across the shock are related by Rankine-Hugoniot relations, which are noting but applications of conservation equations across the shock wave, i.e., conservation of mass, momentum and energy. These are listed below. ρ u =ρ u () p +ρ u =p + ρ u (6) u (E +p )= u (E +p ) (7) 4.. Shock Wave As stated earlier, Rankine-Hugoniot conditions are applied across theshock. After certain mathematical manipulations one gets the relation for star properties as below. QL= ρ *Lu * (8) P *L= ρ L γ γ + +p pl γ p γ +pl + (9) u * = u L ƒ L (p *, W L) () ƒ L(p*,W L)=(p*p L) A L () p + B L A L= () γ+ ρ L Fig. Possible wave structures of Riemann problem The problem of determining the types of waves, their strengths and the flow properties in each region for a given set of initial conditions is called a Riemann problem. 4. Exact Solution of Riemann Solver A number of different exact Riemann solvers have been developed by different investigators. Godunov [] [], Hart Lax Leer pioneered the work in this field. He implicitly related the star pressure on both sides of the contact discontinuity to the mass fluxes passing through left and right waves. Thus a single equation is terms of star pressure are obtained. He used a fixed point iteration scheme for solving the pressure equation. In Godunov s second method, a higher order Newton s iterative method is used. Though this method is computationally accurate than the first one, it is time consuming due to the use of higher order derivatives. 4. Exact Riemann Solver Toro [] presented an exact Riemann solver for ideal gases. This method uses first order Newton- Raphson iterative method for solving the implicit pressure equation. For this analysis, the following assumptions are made: B L= γ γ+ p L(3) Thus the expression of f Lcontains all the known terms, i.e., the conditions to the left of the left nonlinear wave, and also the unknown star pressure. Above equations represents a case when the left wave is shock wave. Similar relations writtenfor the right moving wave is the short wave Expansion Wave The following equations are written assuming that right moving wave structure is an expansion wave. Similar equations exist for a case when the left moving wave is an expansion wave. ρ *=ρ R p p R C = γp ρ γ C *R= a R p p R γ γ () (6) (4) u *- c R = ur- c R (7) γ γ u *= u R+ ƒ R (p *,W R) (8) f R p, W R = C R γ p p R γ γ (9) Thus two expressions for star velocity () and (8). Eliminating the star velocity following expression is obtained: 3
4 Special Edition PGCON-MECH-7 f R p, W R f L p, W L + f R p, W R + u = () u = u R -u L () This () is the required expression in terms of star pressure. HLLC is the improved version of Hart-Lax-Leer Riemann solver by incorporating the contact surface in wave pattern. The wave estimates SL, S*, SR are assumed. The procedure to be adopted for obtaining the HLLC fluxes is as follows, Compute the wave speeds SL, S*, SR. 4. Godunov s Method Godunov's finite volume method for the specific case of time dependent one-dimensional Euler equations are be discussed []. In the discretized form Godunov's Scheme is written as, U n+ i = U n i + F i + F i+ () with inter cell numerical flux is given as F i+ = F U i+ (3) During the time marching one needs to obey the Courant Number (CFL) criterion. t = C cfl x S max (4) Ccfl is known as Courant number which must be in the range of to. Transmissive or transparent boundary arises from need to define finite or very small computational domain. 4..Solution algorithm for one dimensional Godunov's method is as follows:. Physical one dimensional space is discretized into the number of cells.. For the given unsteady problem, the initial flow variable values are assigned as the cell center values. 3. Outside the physical boundary, two fictitious cells, one outside left boundary and the other outside the right boundary are taken. The boundary conditions are taken to be trans- missive as per the requirement. 4. For each of the cell faces, an exact Riemann problem is created with the left cell values as left data set and the right cell values as right data set.. The solution of these Riemann problems gives the corresponding cell face values. 6. The solutions of the Riemann problems also give a maximum wave velocity fir that time step. 7. This maximum velocity is used in Courant number condition (6.4) of stability for deciding the time step size. 8. To calculate variable values at next time level, (6.) is used. 9. The code is be made to run for the desired number of time steps. 4.3 HLLC Riemann Solver Fig.3Solution in the star region for HLLC Riemann solver. Compute the approximate States. Compute the HLLC fluxes F L if S L Fi+ = F L = F L + S L U L U L ifs L S F R if S R F R = F R + S R U R U R ifs S R () Where, F L = F(U L) and F R = F(U R) 3. Compute new values of the states using the formula. U i n+ = U i n + Δt Δx (F i - F i + )(6). Results and Disscusions In order to illustrate the performance of solvers numerical results were presented for one dimensional test problems.. Sod s shock tube problem The test consists of solving one dimensional Euler s equations (), here the value of gamma is taken as.4, the left and right states are separated by diaphragm positioned in the middle of the tube. The domain is discretized with computing cells and the Courant number of computations is.8... One-dimensional Tests The performance of each high-resolution solver is assessed with ten one-dimensional (-D) tests. These include five standard Sod tests. The initial conditions for these nine tests are summarized in Table. A short description of the five sod tests respectively is given below. Test-: The solution of the mild test consists of shock, a contact and rarefaction, all moving towards right. Test-: The solution of this test consist of two symmetric rarefaction and a contact moving with zero velocity. The region between the rarefaction is near vacuum. 4
5 Velocity Pressure Special Edition PGCON-MECH-7 Test-3: The solution for this test consists of a strong shock followed by a moving contact and refraction. It is designed for testing the robustness of solver. Test-4: The solution for this test consists of three strong discontinuities moving to right this test is also designed for assessing the robustness of a solver. Test-: This test similar to test 3 but with the nonlinear velocities, and it assesses a solver s ability to resolve slowly-moving contact discontinuities. Table Initial condition for -D Riemann problems D L U L P L D R U R P R Five test cases are checked using this code whose solution were available in the literature. The results obtained were exactly in agreement with literature []. Table represents the input condition for four test cases tested with the code. Table. Shows the exact solution of the star properties for the four test cases. These results match with the actual calculations and those published in the literature. Table Exact solution Test Ρ* L Ρ* R P* u* E E E E Fig. 4 to 8 present a selection of results for test to. The density, Pressure Velocity and internal energy distributions are plotted and compared with the corresponding exact solutions
6 velocity Pressure Special Edition PGCON-MECH HLLC GODUNOV EACT Fig.4 Property distribution for SOD Test Fig. Property distribution for SOD Test HLLC GODUNOV EACT 6
7 Velocity Pressure Velocity Pressure Special Edition PGCON-MECH Fig.6 Property distribution for SOD Test
8 Special Edition PGCON-MECH Fig.7 Property distribution for SOD Test HLLC GODUNOV EACT Fig.8 Property distribution for SOD Test Conclusion In case of Euler equation the Riemann contain the shock tube problem which is studied. Riemann solvers are used for solving five -D tests and allover five -D tests. Results obtained are exactly in agreement with literature. These results match with the actual calculations and those published in the literature. The conditions in these tests are too hypothetical to exist in a real gas dynamic problem. An agreement to the results of these problems ensures robustness of the solver. The Godunov's method and HLLC (Harten, Lax and Van Leer contact) are used further for benchmarking with SOD -D tests and remaining -D tests. References 8
9 Special Edition PGCON-MECH-7 S. F. Davis,(988), Simplfied second order Godunov Type methods, SIAM J. Sci. Stat. Comput., 9: B.Enfield.,(988),On Godunov Type methods for gas dynamics, SIAM J. Number Anal., (): Harten, A.,(983),High resolution schemes for hyperbolic conversion laws,j. compute. Physso:3-69. Harten, A and Hyman, J.M.,(983), Self-adjusting grid method for one-dimensional hyperbolic conversion laws, s. Compute. Phys., SO: C. B. Laney,(), Computational Gas dynamics, Cambridge University Press, First edition, New York, U.S., 998. R. J. Leveque, finite volume method for hyperbolic problems, Cambridge University Press, First edition, Cambridge, U.K. Liska,R.,Wendroff, B.,(3), Comparison of several differences schemes on -D and -D test problems fortheeuler equations, SIAM. Lax, P.D. and liu,. D.,(998), Solution of two dimensional Riemann Problems of gas dynamics by positive schemes. SIAMJ. scicomput., 9: Takayama,K. and Inoue,O.,(99),Shock wave diffraction over go degree sharp comes.posters presented at 8th ISSW. Shock waves, :3-3. E. F. Toro,(99), Riemann solvers and Numerical Methods for Fluid Dynamics: A Practical Approach, Springer Verlag, Berlin. E. F. Toro, M. Spruce and W. Spears,(994),Restoration of the contact surface in HLL- Riemann solver, shock waves, 4:-34. E. F. Toro and A. Chakraborty.,(994),Development of an Approximate Riemann solver for steady supersonic Euler equations, Aeronoutical Journal, 98: Woodword, P., colella,p.,(984),the numerical simulation of two dimensional fluid flow with strong shocks. J. Compute.Phys., 4:-73 9
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 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 informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More 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 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 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 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 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 informationA Linearised Riemann Solver for Godunov-Type Methods. E.F.Toro
Cranfield»nun «32 College of Aeronautics Report No.9116 September 1991 TECHNISCHE UNIVERSITEIT DELFT LUCHTVAART- EN RUIMTEVAARnECHNIEK BIBLIOTHEEK Kluyverweg 1-2629 HS DELFT A Linearised Riemann Solver
More 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 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 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 informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More 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 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 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 informationNumerical investigation of shock wave dense particles cloud interaction
25 th ICDERS August 2 7, 215 Leeds, UK Numerical investigation of shock wave dense particles cloud interaction Pavel S. Utkin Moscow Institute of Physics and Technology Dolgoprudny, Moscow region, Russia
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 informationThe Euler Equations! Advection! λ 1. λ 2. λ 3. ρ ρu. c 2 = γp/ρ. E = e + u 2 /2 H = h + u 2 /2; h = e + p/ρ. 0 u 1/ρ. u p. t + A f.
http://www.nd.edu/~gtryggva/cfd-course/! Advection! Grétar Tryggvason! Spring! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = x ( / ) ρe ρu E + p ρ E = e + u / H = h + u
More 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 informationAdvection / 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 informationIntroduction to Riemann Solvers
CO 5 BOLD WORKSHOP 2, 2 4 June Introduction to Riemann Sovers Oskar Steiner . Three major advancements in the numerica treatment of the hydrodynamic equations Three major progresses in computationa fuid
More informationA 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 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 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 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 informationApproximate 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 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 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 informationGODUNOV-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 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 informationDEVELOPMENT 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 informationSlip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics
Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics J.J.W. van der Vegt and H. van der Ven Nationaal Lucht- en Ruimtevaartlaboratorium National
More informationChapter 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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More 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 informationA Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws
J Sci Comput (011) 46: 359 378 DOI 10.1007/s10915-010-9407-9 A Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws Fei Teng Li Yuan Tao Tang Received: 3 February 010 / Revised:
More informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
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 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 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 simple numerical scheme for the 2D shallow-water system
A simple numerical scheme for the D shallow-water system Jie Hu, Graduate Research Student, The National Hydraulics and Environmental Laboratory LNHE (Laboratoire National d Hydraulique et Environnement),
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More 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 informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More 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 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 informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
More informationCOMPUTATIONAL 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 informationModeling full-mach-range cavitating flow with sharp interface model
Center for Turbulence Research Proceedings of the Summer Program 28 83 Modeling full-mach-range cavitating flow with sharp interface model By X. Y. Hu, N. A. Adams, E. Johnsen AND G. Iaccarino In this
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 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 informationSimple waves and a characteristic decomposition of the two dimensional compressible Euler equations
Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
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 informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
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 informationA wave propagation method for compressible multicomponent problems on quadrilateral grids
A wave propagation method for compressible multicomponent problems on quadrilateral grids K.-M. Shyue Department of Mathematics, National Taiwan University, Taipei 16, Taiwan Abstract We describe a simple
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 informationThe PLUTO code for astrophysical gasdynamics
Mem. S.A.It. Suppl. Vol. 13, 67 c SAIt 009 Memorie della Supplementi The PLUTO code for astrophysical gasdynamics A. Mignone 1, 1 Dipartimento di Fisica Generale Amedeo Avogadro, Università degli Studi
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 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 informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationMATANIA BEN-ARTZI, JIEQUAN LI AND GERALD WARNECKE
A DIRECT EULERIAN GRP SCHEME FOR COMPRESSIBLE FLUID FLOWS MATANIA BEN-ARTZI JIEQUAN LI AND GERALD WARNECKE Abstract. A direct Eulerian generalized Riemann problem (GRP) scheme is derived for compressible
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
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 informationAntony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
More informationHFVS: An Arbitrary High Order Flux Vector Splitting Method
HFVS: An Arbitrary High Order Flu Vector Splitting Method Yibing Chen, Song Jiang and Na Liu Institute of Applied Physics and Computational Mathematics, P.O. Bo 8009, Beijing 00088, P.R. China E-mail:
More informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
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 informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationDirect Numerical Solution of the Steady 1D Compressible Euler Equations for Transonic Flow Profiles with Shocks
Direct Numerical Solution of the Steady 1D Compressible Euler Equations for Transonic Flow Profiles with Shocks Hans De Sterck, Scott Rostrup Department of Applied Mathematics, University of Waterloo,
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 informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationAn Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations. Thesis by Cheng Chi
An Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations Thesis by Cheng Chi In Partial Fulfillment of the Requirements For the Degree of Masters of Science King Abdullah
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationA WIMF Scheme for the Drift-Flux Two-Phase Flow Model
A WIMF Scheme for the Drift-Flux Two-Phase Flow Model Steinar Eve A,B,D, Tore Flåtten B,A and Svend Tollak Munkeord C A Centre of Mathematics for Applications (CMA), 153 Blindern, NO-316 Oslo, Norway B
More informationJournal of Computational Physics
Journal of Computational Physics 8 (9) 77 798 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A central Rankine Hugoniot solver
More informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationE = 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 informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationu-= (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 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 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 informationIMPORTANCE OF NUMERICAL EFFICIENCY FOR REAL TIME CONTROL OF TRANSIENT GRAVITY FLOWS IN SEWERS
1106 September 11~16, 2005, Seoul, Korea IMPORTANCE OF NUMERICAL EFFICIENCY FOR REAL TIME CONTROL OF TRANSIENT GRAVITY FLOWS IN SEWERS ARTURO S. LEON 1, MOHAMED S. GHIDAOUI 2, ARTHUR R. SCHMIDT 3 and MARCELO
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationThe Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow
The Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow D.W. Schwendeman, C.W. Wahle 2, and A.K. Kapila Department of Mathematical Sciences, Rensselaer Polytechnic
More informationChapter 5: The Hydrodynamical Riemann Problem
Chapter 5: The Hydrodynamical Riemann Problem 5.) Introduction 5..) General Introduction to the Riemann Problem We have seen in Chapter 4 that even Burgers equation, the simplest non-linear scalar conservation
More informationNotes: Outline. Shallow water equations. Notes: Shallow water equations. Notes:
Outline Nonlinear hyperbolic systems Shallow water equations Shock waves and Hugoniot loci Integral curves in phase plane Compression and rarefaction R.J. LeVeque, University of Washington IPDE 2011, July
More 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 Well-Balanced Gas-Kinetic Scheme for the Shallow-Water Equations with Source Terms
Journal of Computational Physics 178, 533 562 (2002) doi:10.1006/jcph.2002.7040, available online at http://www.idealibrary.com on A Well-Balanced Gas-Kinetic Scheme for the Shallow-Water Equations with
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 informationA Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems
A Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems Dia Zeidan,a) and Eric Goncalves 2 School of Basic Sciences and Humanities, German Jordanian University, Amman, Jordan
More information