The 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.
|
|
- Doreen Higgins
- 5 years ago
- Views:
Transcription
1 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 /; h = e + p/ρ p = ρrt; e = e( T ); c v = de /dt h = h( T); c p = dh /dt R = c p c v ; = c p /c v ; p = ( )ρe ρ ρu ρu + ρu + p = x ( / ) ρe ρu E + p ρ p = ( )ρe c = p/ρ Expanding the derivative and rearranging the equations! ρ u ρ ρ u + u /ρ u p ρc x = + A or! x = u p The Characteristics for the Euler Equations are found by finding the eigenvalues for A T! det(a T λi) = Find! det(a T λi) = or:! or! (( ) c ) = ( u λ) u λ ( u λ) = λ = u (( u λ) c ) = u λ = ±c λ = u ± c u λ det ρ u λ ρc = /ρ u λ Therefore! λ = u c; λ = u + c; λ 3 = u The Rankine-Hugoniot conditions! For a hyperbolic system! + F x = The speed of a discontinuity (shock) is found by moving to a frame where a shock moving with speed s is stationary!! s x + F x = Integrating across the shock yields! [ ] = [ F] s f
2 For the Euler equations:! ρ ρu ρu + ρu + p = x ( / ) ρe ρu E + p ρ The Rankine-Hugoniot conditions are:! s( ρ L ρ R ) = (( ρu) L ( ρu) R ) s( ( ρu) L ( ρu) R ) = (( ρu + p) ρu + p L ( ) R ) s( ( ρe ) L ρe ) = (( ρu(e + p/ρ) ) L ρu(e + p/ρ) ) ( ) R ( ) R The Shock-Tube Problem! Exact Solution! The shock tube problem! u L,, ρ L! u R,, ρ R! Density! L! 5! 3!! R! Pressure! Expansion Fan! Contact! Shock! Expansion Fan! Contact! Shock!! 3!! given:!, ρ L, u L,, ρ R, u R α = +! 3!! given:!, ρ L, u L,, ρ R, u R c L = ρ L c R = ρ R Consider the case > :!!Shock separates R and!!contact discontinuity separates and 3!!Expansion fan separates 3 and L! The Rankine-Hugoniot conditions give a nonlinear P = p relation for the pressure ump across the shock! P = p L ( ) ( c R / c L )(P ) ( + ( + ) ( P ) ) which can be solved by iteration! /( )
3 The speed of the shock is! + ( +)P s shock = u R + c R /! 3!! x shock = x + s shock t The speed of the contact is! s contact = u 3 = u = u L + c L P p R x contact = x + s contact t The left hand side of the fan moves with speed! s fl = c L x fl = x + s fl t The right hand side of the fan moves with speed! s fr = u c L x fr = x + s fr t Left uniform state! given:!, ρ L, u L, c L = ρ L! 3!! In the expansion fan ()! + (s x) ρ = ρ c L α t u = u L + s x + t + (s x) p = p c L α t Behind the contact (3)! p 3 = p u 3 = u ρ 3 = ρ L P p R / Behind the contact (3)! p 3 = p u 3 = u ρ 3 = ρ L P p R /! 3!! Behind the shock ()! ρ = + αp α + P ρ R p = P u = u L + c L P p R α = + Right uniform state! given:!, ρ R, u R c R = ρ R The speed of sound is given by! c = p ρ The Mach number is defined as the ratio of the local velocity over the speed of sound! Ma = u c Ma < subsonic! Ma > supersonic! Test case:! Shocktube problem of G.A. Sod, JCP 7:, 97! = 5 ; ρ L =.; u L = = ; ρ R =.5; u R = t final =.5 x =.5 Subsonic case!.... x 3 5 velocity t=.5! Velocity! Pressure! Density! Mach number 3 5 Mach Number! velocity x pressure Mach number
4 Test case:! Pressure! Shocktube problem of G.A. Sod, JCP 7:, 97! = 5 ; ρ L =.; u L = = ; ρ R =.; u R = t final =.5 x =.5 Supersonic case! Velocity! Mach Number! Density! Solve using second order Lax-Wendroff! For the fluid-dynamic system of equations (Euler equations):! ρ ρu ρu + ρu + p = x ( / ) ρe ρu E + p ρ where! E = e + u / ; Add the artificial viscosity to RHS:! = αh x ρ u u x x u p = ( ) ρe * f + Solutions of the D Euler equation using Lax-Wendroff! ( ).5 Δt n n ( F ) =.5 f n n + f + ( ) n f + = f n Δt h F * * + F F' = F αh ρ u h F + With an artificial viscosity term added to the corrector step! u u x x Test case:! Effect of α! Shocktube problem of G.A. Sod, JCP 7:, 97! = 5 ; ρ L =.; u L = = ; ρ R =.5; u R = Final time:.5!.... nx=; α = nx=; α =.5 α =. α =.5 α =.5 x pressure α =.5 α =. α =.5 α =
5 nx=; nx=; nx=5; time=.5 x pressure.... Effect of resolution! α = Outline of L-W program! for istep=:! for i=:nx,p(i)=..; for i=:nx-!!!%prediction step! rh(i)=..! ruh(i)=..! reh(i)=..! for i=:nx,ph(i)=..; for i=:nx-!!!%correction step! r(i)=..! ru(i)=..! re(i)=..! for i=:nx,u(i)=ru(i)/r(i); for i=:nx-!!!%artificial viscosity! ru(i)=..! re(i)=..! time=time+dt,istep! nx=9; artvisc=.5;! hold off! %gg=.;p_left=;p_right=;r_left=;r_right=.;! gg=.;p_left=;p_right=;r_left=;r_right=.5;! xl=.;h=xl/(nx-);time=;! r=zeros(,nx);ru=zeros(,nx);re=zeros(,nx);p=zeros(,nx);! rh=zeros(,nx);ruh=zeros(,nx);reh=zeros(,nx);ph=zeros(,nx);! for i=:nx,r(i)=r_right;ru(i)=.;re(i)=p_right/(gg-); for i=:nx/; r(i)=r_left; re(i)=p_left/(gg-); rh=r;ruh=ru;reh=re;ph=p;! dt=.5*h/sqrt(.*max([7+p_right/r_right,7+p_left/r_left]) )! for istep=:! for i=:nx,p(i)=(gg-)*(re(i)-.5*(ru(i)*ru(i)/r(i))); for i=:nx-!!%prediction step! rh(i)=.5*(r(i)+r(i+))-(.5*dt/h)*(ru(i+)-ru(i));! ruh(i)=.5*(ru(i)+ru(i+))-(.5*dt/h)*((ru(i+)^/r(i+))+p(i+)-(ru(i)^/r(i))-p(i));! reh(i)=.5*(re(i)+re(i+))-...!!(.5*dt/h)*((re(i+)*ru(i+)/r(i+))+(ru(i+)*p(i+)/r(i+))...!!! -(re(i)* ru(i) /r(i))- (ru(i) *p(i) /r(i)));! for i=:nx,ph(i)=(gg-)*(reh(i)-.5*(ruh(i)*ruh(i)/rh(i))); for i=:nx-!!%correction step! r(i)=r(i)-(dt/h)*(ruh(i)-ruh(i-));! ru(i)=ru(i)-(dt/h)*((ruh(i)^/rh(i))-(ruh(i-)^/rh(i-))+ph(i)-ph(i-));! re(i)=re(i)-(dt/h)*((reh(i)*ruh(i)/rh(i))-(reh(i-)*ruh(i-)/rh(i-))+...!!(ruh(i)*ph(i)/rh(i))-(ruh(i-)*ph(i-)/rh(i-)));! for i=:nx,u(i)=ru(i)/r(i); Flux Vector Splitting! for i=:nx-!!%artificial viscosity! ru(i)=ru(i)+artvisc*(.5*dt/h)*( (r(i)+r(i+))*(u(i+)-u(i))*abs(u(i+)-u(i))...! -(r(i)+r(i-))*(u(i)-u(i-))*abs(u(i)-u(i-)) );! re(i)=re(i)+artvisc*(.5*dt/h)*...! ( (u(i+)+u(i))*(r(i)+r(i+))*(u(i+)-u(i))*abs(u(i+)-u(i))...! -(u(i)+u(i-))*(r(i)+r(i-))*(u(i)-u(i-))*abs(u(i)-u(i-)) );! time=time+dt,istep! plot(p,'b','linewidth',),title('pressure'); pause! if(time >.5)break, For upwind schemes, it is necessary to determine the upstream direction. For systems with many characteristics running both left and right, there is not one upstream direction! Upwind Scheme - Revisited! Generalized Upwind Scheme (for both U > and U < )! Define! f f UΔt n n ( f f ), h U > UΔt n n ( f f ), h U < n+ n = f n+ n = f + U + = ( U + U ), U = ( U U ) The two cases can be combined into a single expression:! [ ] n f + = f n Δt h U + ( f n f n ) + U ( f n + f n )
6 A system of hyperbolic equations! For a system of equations there are generally waves running in both directions. To apply upwinding, the fluxes must be decomposed into left and right running waves! + F x = can be written in the form! + [ A] x = ; The system is hyperbolic if! Steger-Warming (979)! [ A] = F [ T ] [ A] [ T] = [ λ]; [ T ] = T q T q N F = [ A]f = T The matrix of eigenvalues Hence! [ A] = A + [ ][ λ] [ T] f [ λ] [ λ] = [ λ + ] + [ λ ] [ ] + [ A ] = [ T ][ λ + ] T is divided into two matrices! [ ] T [ ] + [ T] λ [ ] ( ) Define! F = F + + F F + = [ A + ]f, F = [ A ]f Conservation law becomes! + F + x + F x = Example: -D Hyperbolic Equation! Leading to! v t w t + c f c f x = v x w x = f = v w ; F = c w v ; v = ; w = x [ A] = F = c [ T ] = q T = c T λ + q c ; [ ] = c [ A + ] = [ T ] [ λ + ][ T] ; [ A ] = [ T] [ λ ][ T ] F + = [ A + ]f; F = [ A ]f Leading to:! F + = cv c w ; F = cv c w v + cw v cw c [ λ ] = Solve using first order upwinding with flux splitting! For a nonlinear system of equations such as the Euler equations, there is some arbitrariness in how the system is split!
7 Here we solve the one-dimensional Euler equation using the van Leer vector flux splitting! Van Leer! F + = ρ (u + c) c ( )u + c ( )u + c ( ) F = ρ (u c) c ( )u c c ( )u ( ) van Leer vector flux splitting! ρu ρu + p = ρ ρu(e + p c ρ ) (u + c) For example, the mass flux:! ρ c (u + c) ρ (u c) c ( )u + c ( )u + c ( ) ρ (u c) c ( ) = ρ c = ρ c u + uc + c u + uc c ( )u c c ( )u ( ) uc = ρu Steger-Warming! F + = ρ Several other splitting schemes are possible, such as:! ( )u + c ( )u + (u + c) ( )u 3 + (u + c)3 + 3 (u + c)c ( ) F = ρ (u c) u c (u c) + 3 c F + = ρ (u + c) c Rewrite the flux terms in terms of Mach number:! ( )u + c ( )u + c ( ) = ρc u c + c + u c c + u c = ρc ( M + ) c + M c + M + F + x + F - x = ρ ρu + ρc ρe x ( M + c ) + M + x ρc ( M c ) + M = c + M c M + F+ x + F- x = Solve by! f n+ = f n Δt h F + + F - Where! ( ) n Δt ( ) n - - h F + F F + = ρc ( M + c ) + M ; F = ρc ( M c ) + M c + M c M
8 for istep=:maxstep for i=:nx,c(i)=sqrt( gg*(gg-)*(re(i)-.5*(ru(i)^/r(i)))/r(i) );end for i=:nx,u(i)=ru(i)/r(i);end; for i=:nx,m(i)=u(i)/c(i);end for i=:nx- %upwind rn(i)=r(i)-(dt/h)*(... (.5*r(i)*c(i)*(m(i)+)^) - (.5*r(i-)*c(i-)*(m(i-)+)^)+... (-.5*r(i+)*c(i+)*(m(i+)-)^) - (-.5*r(i)*c(i)*(m(i)-)^) ); run(i)=ru(i)-(dt/h)*(... (.5*r(i)*c(i) *(m(i)+)^) *(( +.5*(gg-)*m(i)) **c(i) /gg) -... (.5*r(i-)*c(i-)*(m(i-)+)^)*(( +.5*(gg-)*m(i-))**c(i-)/gg) +... (-.5*r(i+)*c(i+)*(m(i+)-)^)*((-+.5*(gg-)*m(i+))**c(i+)/gg) -... (-.5*r(i)*c(i) *(m(i)-)^) *((-+.5*(gg-)*m(i)) **c(i) /gg) ); ren(i)=re(i)-(dt/h)*(... (.5*r(i)*c(i) *(m(i)+)^) *((+.5*(gg-)*m(i))^ **c(i)^ /(gg^-)) -... (.5*r(i-)*c(i-)*(m(i-)+)^)*((+.5*(gg-)*m(i-))^**c(i-)^/(gg^-)) +... (-.5*r(i+)*c(i+)*(m(i+)-)^)*((-.5*(gg-)*m(i+))^**c(i+)^/(gg^-)) -... (-.5*r(i)*c(i) *(m(i)-)^) *((-.5*(gg-)*m(i))^ **c(i)^ /(gg^-)) ); end Effect of resolution! Shocktube problem of G.A. = 5 ; ρ L =.; u L = Sod, JCP 7:, 97! = ; ρ R =.5; u R = x Pressure nx=; maxstep=3! nx=5; maxstep=! Final time:.5! Density... nx=; maxstep=3! nx=; maxstep=!. nx=5; maxstep=! nx=5; maxstep=5! Summary! Numerical solutions of the onedimensional Euler equations! Upwind/flux splitting! Lax-Wendroff/artificial viscosity!
( ) κ ( f j +1. ( ) 1 + κ ( f j +2. ( ) j +1/ 2 ) = Uf j +1/ 2. = r j +1/ 2. Ψ( r) = r + r 1 + r. Ψ( r) = minmod(r,1) j 1/ 2 ( F n +1/ 2 F n
ttp://users.wpi.edu/~gretar/me6.tml! Numerical Metods! or Hyperbolic! Equatios VIII! Grétar Tryggvaso! Sprig! Upwid or oliear systems:!!flux vector splittig!!goduovʼs metod! Higer Order Metod!!Higer order
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 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 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 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 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 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 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 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 informationShock Waves. 1 Steepening of sound waves. We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: kˆ.
Shock Waves Steepening of sound waves We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: v u kˆ c s kˆ where u is the velocity of the fluid and k is the wave
More informationThe 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 informationAA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow
AA210A Fundamentals of Compressible Flow Chapter 13 - Unsteady Waves in Compressible Flow The Shock Tube - Wave Diagram 13.1 Equations for irrotational, homentropic, unsteady flow ρ t + x k ρ U i t (
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 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 information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More 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 informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More 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 informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
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 informationDrag of a thin wing and optimal shape to minimize it
Drag of a thin wing and optimal shape to minimize it Alejandro Pozo December 21 st, 211 Outline 1 Statement of the problem 2 Inviscid compressible flows 3 Drag for supersonic case 4 Example of optimal
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More 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 informationRadiative & Magnetohydrodynamic Shocks
Chapter 4 Radiative & Magnetohydrodynamic Shocks I have been dealing, so far, with non-radiative shocks. Since, as we have seen, a shock raises the density and temperature of the gas, it is quite likely,
More informationNumerical Hydrodynamics: A Primer. Wilhelm Kley Institut für Astronomie & Astrophysik & Kepler Center for Astro and Particle Physics Tübingen
A Primer Wilhelm Kley Institut für Astronomie & Astrophysik & Kepler Center for Astro and Particle Physics Tübingen Bad Honnef 2016 Hydrodynamics: Hydrodynamic Equations The Euler-Equations in conservative
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 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 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 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 informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
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 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 informationIsentropic Duct Flows
An Internet Book on Fluid Dynamics Isentropic Duct Flows In this section we examine the behavior of isentropic flows, continuing the development of the relations in section (Bob). First it is important
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
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 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 informationIntroduction to Gas Dynamics All Lecture Slides
Introduction to Gas Dynamics All Lecture Slides Teknillinen Korkeakoulu / Helsinki University of Technology Autumn 009 1 Compressible flow Zeroth law of thermodynamics 3 First law of thermodynamics 4 Equation
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 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 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 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 informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationSELF-SIMILAR SOLUTIONS FOR THE TRIPLE POINT PARADOX IN GASDYNAMICS
SIAM J. APPL. MATH. Vol. 0, No. 0, pp. 000 000 c XXXX Society for Industrial and Applied Mathematics SELF-SIMILAR SOLUTIONS FOR THE TRIPLE POINT PARADOX IN GASDYNAMICS ALLEN M. TESDALL, RICHARD SANDERS,
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 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 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 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 informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More 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 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 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 informationFundamentals of compressible and viscous flow analysis - Part II
Fundamentals of compressible and viscous flow analysis - Part II Lectures 3, 4, 5 Instantaneous and averaged temperature contours in a shock-boundary layer interaction. Taken from (Pasquariello et al.,
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 informationComputational Astrophysics 7 Hydrodynamics with source terms
Computational Astrophysics 7 Hydrodynamics with source terms Oscar Agertz Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term
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 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 information1. (20 pts total 2pts each) - Circle the most correct answer for the following questions.
ME 50 Gas Dynamics Spring 009 Final Exam NME:. (0 pts total pts each) - Circle the most correct answer for the following questions. i. normal shock propagated into still air travels with a speed (a) equal
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
More informationShock and Expansion Waves
Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid
More 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 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 informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationGravity Waves Gravity Waves
Gravity Waves Gravity Waves 1 Gravity Waves Gravity Waves Kayak Surfing on ocean gravity waves Oregon Coast Waves: sea & ocean waves 3 Sound Waves Sound Waves: 4 Sound Waves Sound Waves Linear Waves compression
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
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 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 informationUnsteady Wave Motion Shock Tube Problem - Shock Reflection
Unsteady Wave Motion Shock Tube Problem - Shock Reflection Niklas Andersson Division of Fluid Dynamics Department of Applied Mechanics Chalmers University of Tecnology 8 februari 09 Shock Reflection When
More informationDetonation Simulations with a Fifth-Order TENO Scheme
Commun. Comput. Phys. doi: 10.4208/cicp.OA-2018-0008 Vol. x, No. x, pp. 1-37 xxx 201x Detonation Simulations with a Fifth-Order TENO Scheme Haibo Dong 1, Lin Fu 2, Fan Zhang 3,, Yu Liu 4 and Jun Liu 1
More informationNumerical resolution of a two-component compressible fluid model with interfaces
Numerical resolution of a two-component compressible fluid model with interfaces Bruno Després and Frédéric Lagoutière February, 25 Abstract We study a totally conservative algorithm for moving interfaces
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 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 information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
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 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 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 informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
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 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 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 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 informationShock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids
Shock and Rarefaction Waves in a Hyperbolic Model of Incompressible Fluids Andrea Mentrelli Department of Mathematics & Research Center of Applied Mathematics (CIRAM) University of Bologna, Italy Summary
More 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 finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More 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 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 informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationEulerian interface-sharpening methods for hyperbolic problems
Eulerian interface-sharpening methods for hyperbolic problems Application to compressible multiphase flow Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan 11:55-12:25, March 05,
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 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 informationQuasi-linear first order equations. Consider the nonlinear transport equation
Quasi-linear first order equations Consider the nonlinear transport equation u t + c(u)u x = 0, u(x, 0) = f (x) < x < Quasi-linear first order equations Consider the nonlinear transport equation u t +
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationAngular momentum preserving CFD on general grids
B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration
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 information