Solving the Euler Equations!
|
|
- Tamsyn Conley
- 6 years ago
- Views:
Transcription
1 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 ρ E = e + u / h = e + p / ρ; H = h + u / = E + p / ρ; p = ρt; e = e( T ); c v = de /dt h = h( T); c p = dh /dt = c p c v ; γ = c p /c v ; p = ( γ )ρe The Euler equations for D flow:! or! where! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p ρ t f + x F = 0 ρ ρu f = ρu ; F = ρu + p ρe ρu(e + p / ρ) Predictor step! Variables! f j +/ Predictor-corrector method: imiting the variables! Find:! Final step! n +/ f j = f n j Δt h F n n j +/ F j / n +/ = f j + Ψ n +/ f n +/ ( j f j ) n +/ f j +/ = f j + Ψ n +/ f n +/ j + f j F j +/ n +/ = F ( f ) n +/, f n +/ ( j +/ j +/ ) n f + j = f n j Δt h F n +/ n +/ ( j +/ F j / ) 5! 6! The Euler equations can be solved using the flux limited high order methods described earlier by finding the fluxes using solutions to the iemann problem! F j +/ n +/ = F ( f ) n +/, f n +/ ( j +/ j +/ ) In principle we can solve this problem, the iemann problem, exactly by assuming constant states and then integrate the fluxes over the time step. In the last lecture we did so for the case when the fluids are initially at rest.! Although the iemann problem can be solved for the general case when the fluids are not stationary, the solution is expensive (involving solving a nonlinear equation for the pressure ratio across the shock) and the only information that we need from the solution is the flux across the cell boundary.! Therefore, usually we use approximate iemann solvers. Different approximations are possible but generally the result is a complex process. Here we will only outline it briefly!
2 7! 8! The oe approximate iemann solver was one of the first method to compute the fluxes in a simpler way.! It is based on approximating the Euler equation by a linear equation! f t + A f x = 0 whose fluxes can be found analytically. The linearization is done in such a way that the correct wave speed is preserved. The jump in the solution across a shock can be written in terms of the eigenvectors:! 3 f f = α p p p= oe introduced the following matrix A! A = 0 0 γ 3 u ( 3 γ )u γ γ u 3 uh H ( γ )u γ u where the average variables are defined by! u = ρ u + ρ u ; H = ρ H + ρ H ρ + ρ ρ + ρ 9! The fluxes are given by! 3 F i+/ = F + F λ p α p where the eigenvalues are the same as for the original! and the scaled eigenvectors are given by! = u c H u c ; = p= u u p λ = u c ; λ = c ; λ 3 = u + c ; ; = u + c H + u c ; α α α 3 0! The alphas are given by! = u 4c + ( γ ) u c c + ( γ ) u c γ u γ c u 4c ( γ ) u c c ( γ ) u c Δρ j+/ = ρ j+/ ρ j+/ ; Δρu j+/ = ρu j+/ γ c u ( γ ) u c c ρu j+/ γ c Where the jumps across the cell boundary are:! and (again)! u = ρ u + ρ u ; H = ρ H + ρ H ; c = γ ρ + ρ ρ + ρ H u Δρ j+/ Δρu j+/ ΔρE j+/ ; ΔρE j+/ = ρe j+/ ρe j+/! The fluxes are given by! F i+/ F i+/ 3 F i+/ = ρu ( + ( ρu) j+/ j+/ ) λ α λ α λ α 3 3 ( j+/ ) = ( γ ) ( ρe) j+/ ( γ 3) ( ρu) ( ρ) j+/ j+/ + ( γ ) ( ρe) j+/ = ρu ρu γ 3 ( ρ) j+/ λ α ( u c ) λ α u λ α ( u + c ) 3 3 H ( + ( ρu) j+/ j+/ H ) λ α ( h u c ) λ α u λ α ( h + u c ) 3 3! rd order unge-kutta! f () = f n + Δt f n f () = f () + Δt f () f n+ = f n + f () 3 rd order unge-kutta! f () = f n + Δt f n f () = 3 4 f n + 4 f () + () Δt f 4 f n+ = 3 f n + 3 f () + () Δt( f ) 3 Time Integration! 4 rd order unge-kutta! f (0) = f n f () = f (0) + 4 Δt f n f () = f (0) + () Δt( f ) 3 f (3) = f (0) + () Δt f f n+ = f (0) + Δt f ()
3 3! 4! esults from:! PINCIPES OF COMPUTATIONA! FUID DYNAMICS by P. Wesseling,! Delft University of Technology,! The Netherlands! 48 gridpoints! Matlab codes available at the website! 5! 6! 48 gridpoints! The oe approximate iemann solver generally gives well behaved results but it does allow for expansion shocks in some cases. This can be corrected by the so-called entropy fix (see Wesseling)! 7! 8! The TVD property gives monotonic behavior for D scalar conservation laws but generally not for nonlinear systems or multidimensional flows. However, generally it is found that schemes that satisfy the TVD property when applied to D scalar equation do well for systems. When used for multidimensional flows by splitting, they are also generally found to be well behaved.! Several variants of similar methods have been proposed.! One of the better known is the AUSM (Advection Upstream Splitting Method) family of methods, originally due to iou and Steffen (993), where the fluxes are split into convective and a pressure part. The method and its derivatives have been used widely for many aerospace and other applications, including multiphase flows.!
4 The Euler Equations in twodimensions! For the Euler equations:! ρ ρu ρv ρu + ρu + p + ρuv = 0 t ρv x ρuv y ρv + p ρe ρu(e + p /ρ) ρv(e + p/ρ) where! E = e + (u + v ) /; p = (γ )ρe In vector form! f t + F x + G y = 0 Two-dimensional problems are almost always solved by splitting where we apply one-dimensional methods to each direction in sequence! Flux vector splitting! F ± = ± ρ 4c u ± c ( γ )u ± c γ v v [ + ( γ )u ± c ] ( γ ) f t + F + x + F x + G + y + G y = 0 G ± = ± ρ 4c v ± c u ( γ )v ± c γ u [ + ( γ )v ± c ] ( γ ) esults from a review paper:! P. Woodward and P. Colella:! Computational domain! Shock! h=/30 h The Numerical Simulation of two- Dimensional Fluid Flow with Strong Shocks.! J. Comput. Phys. 54 (984), 5-73.! Given inflow! 3! wall!! Outflow! h=/60 h=/0 0 Godunovʼs method! MacCormacʼs method!
5 6! h=/30 h=/60 h Higher Order and more recent methods! h=/0 0 MUSC! PPM! 7! 8! Beyond linear: econstruction of higher order approximations for the function in each cell (ENO and WENO).! The critical step in the methods discussed so far is the construction of a linear slope in each cell and the limitation of this slope to prevent oscillations. For higher order methods, a higher order profile needs to be constructed! Constructing an interpolation polynomial from the cell averages: For anything higher than second order (linear) the problem is that the average value in the cell is not equal to the value at the center.! To get around this we look at the primitive function:! x v( ξ) V x The lower bound is arbitrary and can be replaced! dξ j- j- j j+ j+ j-3/ j-/ j+/ j+3/ 9! 30! Since this is the integral over the cells, the discrete version is exact at the cell boundaries! i x j+/ = v( ξ) V x i+/ dξ = f i Δx j= x j / j- j- j j+ j+ j-3/ j-/ j+/ j+3/ i j= V ( x i 3/ ); V ( x i / ); V ( x i+/ ); V ( x i+3/ ); A polynomial interpolating the edge values is given by P(x) and we denote its derivative by p(x) Then it can be shown that! x j+/ p(x) = P'(x) p ξ dξ = P' ξ dξ = P(x j+/ ) P(x j / ) x j / x j+/ x j / x j+/ = V (x j+/ ) V (x j / ) = v ξ dξ = f i Δx x j / That is, the integral of p(x) over the cell is equal to the cell average f i
6 3! 3! Thus, p(x) gives the correct average value in each cell and the integrated value gives the exact values of the primitive function at the cell boundaries.! We need to write down a polynomial P(x) that interpolates the values of the primitive function of the cell boundaries and then differentiate this polynomial to get p(x), which lets us compute the variables at the cell boundary! j- j- j j+ j+ j-3/ j-/ j+/ j+3/ The interpolation polynomial is often taken to be the agrangian Polynomial! k P(x) = V (x i r + m / ) m=0 k l =0 l m x x i r +l / x i r / x i r +l / Where r determines where we start and k is the order! 33! Example: Second order ENO! j- j- j j+ j+ j-3/ j-/ j+/ j+3/ The question is now which point we select. We start by interpolating over one cell (linear). To add one point we can add either the point to the left or the right. In ENO we select the points based on the minimum absolute value of the divided differences of the function values!. Construct left and right slopes by connecting the average values in adjacent cells!. Select the downstream flux by using the smaller slope! Δf + j = f j + f j Δf j-! j! j+! j = f j f j f j +/ = U>0! f j f j f j + amin Δf + ( j, Δf j ), if u j + u j + f j amin Δf + ( j +, Δf j + ), if u + u j j + f j / f j +/ > 0 < 0 f j+ f t + u f x = 0 f * j = f n j Δt h u n j f j n+ = f j n Δt h Second order ENO scheme for the linear advection equation! f n n ( j +/ f j / ) f n j +/ f n * * * ( ( j / ) + u j ( f j +/ f j / )) u n j f j +/ = = amin a,b a, a < b b, b a f j + amin Δf + ( j,δf j ), if u + u j j + f j amin Δf + ( j +,Δf j + ), if u + u j j + > 0 < 0 Second order ENO scheme for the linear advection equation! f t + u f x = 0 Blue: nd order ENO! ed: st Upwind! Δf j + = f j + f j Δf j = f j f j
7 37! 38! j- j- j j+ j+ j-3/ j-/ j+/ j+3/ In the WENO (weighted essentially non-oscillating) scheme we use all the points but weigh the contribution of each according to a smoothness criteria. High-order WENO represents the current state-of-the-art in computing of flows with sharp interfaces! Increasingly we see methods developed for the inviscid Euler equation with shocks being used for the advection part of the Navier-Stokes solvers.! For a recent reference on WENO see, for example:! C-W Shu. High Order Weighted Essential Nonoscillatory Schemes for Convection Dominated Problems. SIAM eview, Vol. 5 (009), 8-6.! 39! 40! CIP-gradient augmentation! Other Approaches! The CIP (Constrained Interpolation Polynomial)! Method (Yabe)! In addition to advecting the marker function f, its derivative is advected by fitting a third order polynomial through the function and its derivatives.! Start with! Introduce! In D, the advection of the derivative is given by! f and g given! Therefore, the derivative is translated with velocity u, just as the function. In D splitting is used to separate translation and deformation! New f and g! 4! CIP-gradient augmentation! The CIP method results in very accurate advection and for a sharp interface it greatly reduces overshoots, but does not eliminate them completely! Enormous progress has been made in solution techniques for hyperbolic systems with shocks in the last twenty years. Advanced methods are now able to resolve complex shocks within a grid space or two, even in multidimensional situations for a large range of governing parameters and physical complexity.! Here, we have only examined relatively elementary aspects of methods for hyperbolic systems, but this short introduction should have taught you methods to solve such systems and introduced you to literature.!
8 43! A few words about the project:! I asked you to select your own parameters. In doing so keep in mind that while Nature usually presents us with problems that have multiphysics problems with a large range of scales and extreme values of the physical properties, in CFD we prefer single physics, small range of scales and O() values for the physical properties. In the project you are simply showing that the code works and looking at trends. Using Matlab you are limited to modest resolution and run-time so you should select you parameters commensurate with the resolution you are using!!
HFVS: 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationComparison of (Some) Algorithms for Edge Gyrokinetics
Comparison of (Some) Algorithms for Edge Gyrokinetics Greg (G.W.) Hammett & Luc (J. L.) Peterson (PPPL) Gyrokinetic Turbulence Workshop, Wolfgang Pauli Institute, 15-19 Sep. 2008 w3.pppl.gov/~hammett Acknowledgments:
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
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 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 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 informationNumerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand
Institute of Aerospace Thermodynamics Numerical modelling of phase change processes in clouds Challenges and Approaches Martin Reitzle Bernard Weigand Introduction Institute of Aerospace Thermodynamics
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 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 informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
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 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 informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More 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 informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
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 informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
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 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 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 informationSub-Cell Shock Capturing for Discontinuous Galerkin Methods
Sub-Cell Shock Capturing for Discontinuous Galerkin Methods Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 39, U.S.A. A shock capturing strategy for higher order
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 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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationActive Flux for Advection Diffusion
Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA
More informationEffect of Applied Magnetic Field on Shock Boundary Layer Interaction
University of Kentucky UKnowledge Mechanical Engineering Faculty Publications Mechanical Engineering 1-2012 Effect of Applied Magnetic Field on Shock Boundary Layer Interaction Ovais U. Khan University
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
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 informationAffordable, entropy-consistent, Euler flux functions
Affordable, entropy-consistent, Euler flux functions (with application to the carbuncle phenomenon) Phil Roe Aerospace Engineering University 0f Michigan Ann Arbor Presented at HYP 2006 1 Entropy pairs
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 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 informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
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 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 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 informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationLow-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln NASA Publications National Aeronautics and Space Administration 1999 Low-Dissipative High-Order Shock-Capturing Methods
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 informationHIGH ORDER FINITE VOLUME SCHEMES
HIGH ORDER FINITE VOLUME SCHEMES Jean-Pierre Croisille Laboratoire de mathématiques, UMR CNRS 71 22 Univ. Paul Verlaine-Metz LMD, Jan. 26, 2011 Joint work with B. Courbet, F. Haider, ONERA, DSNA (Numerical
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 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 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 informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
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 Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationChapter 2 Finite-Difference Discretization of the Advection-Diffusion Equation
Chapter Finite-Difference Discretization of the Advection-Diffusion Equation. Introduction Finite-difference methods are numerical methods that find solutions to differential equations using approximate
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 informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More 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 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 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 limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
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 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 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 informationThe Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1
The Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 94305
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 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 informationA Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes
A Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory SIAM Annual Meeting
More informationA Comparison of Artificial Viscosity Sensors for the Discontinuous Galerkin Method. Joshua Favors
A Comparison of Artificial Viscosity Sensors for the Discontinuous Galerkin Method by Joshua Favors A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
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 informationSpectral Difference Method for Unstructured Grids II: Extension to the Euler Equations
Journal of Scientific Computing, Vol. 32, No. 1, July 2007 ( 2006) DOI: 10.1007/s10915-006-9113-9 Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations Z. J. Wang, 1,4
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 Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationarxiv: v2 [math.na] 24 Mar 2016
arxiv:1504.04107v2 [math.na] 24 Mar 2016 Strong stability preserving explicit linear multistep methods with variable step size Yiannis Hadjimichael David I. Ketcheson Lajos Lóczi Adrián Németh March 21,
More informationON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY ON INTERIOR NODAL SETS
6th European Conference on Computational Mechanics (ECCM 6 7th European Conference on Computational Fluid Dynamics (ECFD 7 11 15 June 018, Glasgow, UK ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationIntroduction to Physical Acoustics
Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
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 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 informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
More informationWhat is a flux? The Things We Does Know
What is a flux? Finite Volume methods (and others) (are based on ensuring conservation by computing the flux through the surfaces of a polyhedral box. Either the normal component of the flux is evaluated
More information