Entropy stable schemes for compressible flows on unstructured meshes
|
|
- Cornelius Gaines
- 5 years ago
- Views:
Transcription
1 Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep SCPDE-2014, Hong Kong 9th November 2014
2 Work done with: Praveen Chandrashekar, TIFR-CAM, Bangalore Siddhartha Mishra, Seminar for Applied Mathematics, ETH Zurich Ulrik S. Fjordholm, NTNU, Trondheim Funded by: AIRBUS Group Corporate Foundation Chair in Mathematics of Complex Systems, established in TIFR/ICTS, Bangalore
3 Conservation laws Consider the (hyperbolic) system Examples U t + f 1 x + f 2 y = 0 (x, t) R2 R + U(x, 0) = U 0 (x) x R 2 Shallow water equations (Geophysics) Euler equations (Aerodynamics) MHD equations (Plasma physics) Non linearities = disc. soln. = weak (distributional) soln.
4 Entropy framework Entropy-entropy flux pair (η(u), q(u) = (q 1 (U), q 2 (U)) t η(u) + div (q(u)) 0 where V = η (U) entropy variables. Single out a physically relevant solution d dt η(u)dx 0 = U(., t) L 2 C No global existence, uniqueness results for generic multidimensional systems
5 Entropy stable schemes Cartesian meshes Tadmor (1987): Entropy conservative flux. Lefloch et al. (2002): Higher order entropy conservative schemes. Explicit entropy stable schemes by Roe and Ismail (2009), PC (2013). Fjordholm et al. (2011): Entropy stable ENO schemes (TeCNO). Madrane et al. (2012): First order entropy stable schemes on unstructured meshes.
6 Outline 2D conservation laws: Discretisation (unstructured) Entropy conservative/stable schemes High-order schemes (sign property) Numerical results 2D Navier-Stokes equations: Global entropy relations Discrete analogue Boundary conditions Numerical results
7 Discretization Discretize domain into triangles T, with nodes i, j, k,etc Primary cell: k n T i T l n T e i T e n T e j j e n e i
8 Discretization Discretize domain into triangles T, with nodes i, j, k,etc Primary cell: Dual control volumes C i Median dual Voronoi dual
9 Discretization Discretize domain into triangles T, with nodes i, j, k,etc Primary cell: Dual control volumes C i Median dual Voronoi dual Suitable for complex geometries!!
10 Semi-discrete scheme du i dt + 1 C i F ij = 0 j i U i cell average over C i F ij = F(U i, U j, n ij ) is the numerical flux satisfying 1 Consistency: 2 Conservation: F(U, U, n) = F(U, n) := f 1 (U)n 1 + f 2 (U)n 2 F(U 1, U 2, n) = F(U 2, U 1, n) U 1, U 2, n
11 First order entropy stable flux 1 Entropy conservative flux F ij = F ij s.t. Sufficient condition (Tadmor) where dη(u i ) dt + 1 C i qij = 0 j i V ijf ij = ψ(u j, n ij ) ψ(u i, n ij ) ψ(u, n) := V(U) F(U, n) q(u, n) }{{} q 1n 1+q 2n 2 (entropy potential) Entropy variable based dissipation F ij = F ij 1 2 D ij V ij, D ij = D ij 0 = dη(u i) dt + 1 C i q ij 0 1 A. Madrane, UF, SM, and E. Tadmor. Entropy conservative and entropy stable finite volume schemes for multi-dimensional conservation laws on unstructured meshes. (in review) j i
12 First order entropy stable flux Kinetic energy and entropy conservative flux (PC) for Euler equations. Dissipation operator D ij = R ij Λ ij R ij R ij scaled eigenvectors of F U (Barth) Λ ij = diag[ λ 1,..., λ n ] Roe type KEPES Kinetic energy and entropy conservative flux + Roe dissipation
13 High-order diffusion V ij O( x ij ) = first order flux Idea: For each x ij, reconstruct V in C i, C j with polynomials V i (x), V j (x) respectively. V ij = V i (x ij ), V ji = V j (x ij ), V ij = V ji V ij Sign property (Fjordholm et al.) For each x ij, define the scaled entropy variables Z = Rij V. Then numerical flux F ij = F ij 1 2 R ijλ ij Z ij = F ij 1 2 D ij V ij, V ij = (R T ij) 1 Z ij is entropy stable if the sign property holds for Z componentwise. sign( Z ij ) = sign( Z ij )
14 Second order (limited) reconstruction i 1 i j + 1 j Define Z i = RijV i, Z j = RijV j The componentwise reconstructed scaled variables are Z ij =Z i + 1 ) ( 2 minmod f ij, } b ij Z ji =Z j 1 ) need h Z i = Rij ( h V i 2 minmod f ji, b ji
15 Transonic flow past NACA-0012 airfoil angle of attack = 2 degrees, freestream M = 0.85 KEPES KEPES-TeCNO Mach number, 20 equally spaced contours between 0.5 and 1.5 ROE (MUSCL)
16 Subsonic flow past a cylinder freestream M = 0.3, symmetric, isentropic flow KEPES KEPES-TeCNO KEPES2 Scheme Minimum Maximum Percent deviation from s KEPES % % KEPES-TeCNO % % KEPES % % Table : Physical entropy bounds, with freestream s =
17 2D Navier-Stokes Equations Initial boundary valued problem U t + f 1 x + f 2 y = g 1 x + g 2 x = (x, y) Ω R 2, t R + y U(x, 0) = U 0 (x) x Ω B(U(x, t)) = h(x, t) x Ω
18 2D Navier-Stokes Equations Initial boundary valued problem U t + f 1 x + f 2 y = g 1 x + g 2 x = (x, y) Ω R 2, t R + y U(x, 0) = U 0 (x) x Ω B(U(x, t)) = h(x, t) x Ω Specific entropy pair for symmetrization η(u) = ρs γ 1, q 1(U) = ρus γ 1, q 2(U) = ρvs γ 1
19 2D Navier-Stokes Equations Initial boundary valued problem U t + f 1 x + f 2 y = g 1 x + g 2 x = (x, y) Ω R 2, t R + y U(x, 0) = U 0 (x) x Ω B(U(x, t)) = h(x, t) x Ω Specific entropy pair for symmetrization η(u) = ρs γ 1, q 1(U) = ρus γ 1, q 2(U) = ρvs γ 1 Viscous fluxes in terms of entropy variables g 1 = K 11 (V) V x + K 12(V) V y, where [ ] K11 K K = 12 0 K 21 K 22 g 2 = K 21 (V) V x + K 22(V) V y
20 Global entropy relation Integrating NSE against V gives us d η = q(u, n) K V, V + V G(U, n) dt Ω Ω Ω Ω q(u, n) + V G(U, n) where and ( ) x V V = y V Ω Ω G(U, n) = g 1 (U)n 1 + g 2 (U)n 2 R 8 obtained from V = ( x V, y V ) R 4 2
21 Semi-discrete scheme Notations j i = { all vertices j neighbouring vertex i } i T = { all vertices i belonging to triangle T } T i = { all triangles T having vertex i } Γ = { all boundary edges of the primary mesh } Γ i = { all boundary edges of the primary mesh having vertex i } C i du i dt = j i F ij + T i G T nt i 2 F ie + ne G e 2 e Γ i e Γ i
22 Choosing flux terms Interior inviscid flux F ij is an entropy stable flux
23 Choosing flux terms Interior inviscid flux F ij is an entropy stable flux Interior viscous flux G T = (G T 1, G T 2 ) is chosen as G T α = K T α1 h xv T + K T α2 h y V T, α = 1, 2
24 Choosing flux terms Interior inviscid flux F ij is an entropy stable flux Interior viscous flux G T = (G T 1, G T 2 ) is chosen as G T α = K T α1 h xv T + K T α2 h y V T, α = 1, 2 Boundary inviscid flux F ie is chosen as F ie = F E ie = 1 F ie = F ie ( U, U b, n e ) F ρ ie 2 (pb i n e) + u i F ρ ie F E ie [ ( ρ u γp γ 1, F ρ ie (ρ = ub n e 2 ) u b ne 2 ] BC implemented weakly i ) [ + (p b p) u n e 2 i ] i
25 Choosing flux terms Interior inviscid flux F ij is an entropy stable flux Interior viscous flux G T = (G T 1, G T 2 ) is chosen as G T α = K T α1 h xv T + K T α2 h y V T, α = 1, 2 Boundary inviscid flux F ie is chosen as F ie = F ie ( U, U b, n e ) Boundary viscous flux G e = (G e 1, G e 2) is chosen as G e n e = G Te n e + C e
26 Discrete global entropy relation Pre-multiplying the scheme by V i and summing over all nodes. [ d ( ) ( ) ] ρs u b n η i C i e ρs u b n e dt γ 1 2 γ 1 2 i e Γ i j [ ] (V b i + V b j) + G Te (Vi + Vj) n e + C e 2 2 e Γ which is consistent with d η dt Ω Ω q(u, n) + Ω V G(U, n)
27 Homogeneous boundary conditions Slip BC (Euler) Continuous setup d η 0 dt Ω Discrete setup d dt η i C i 0 i
28 Homogeneous boundary conditions Slip BC (Euler) Continuous setup d η 0 dt Ω Discrete setup d dt η i C i 0 i No-slip BC with zero heat flux Continuous setup G e n e = 0 τ n e 0 d dt Ω η 0 ( C e = Discrete setup G e 0 n e = τ h n e, 0 ( ( κ h V Te,(4) ) )) n R(V e,(4) ) 2 e d dt η i C i 0 i
29 Numerical results: Lid driven cavity Top lid moving with velocity (u = 1, v = 0). No-slip BC on other 3 faces. Comparing with data of Ghia et al y Referrence Roe KEPES TECNO v Referrence Roe KEPES TECNO u u on vert. line through center x v on hor. line through center
30 Flow past a cylinder Inflow: ( ) y(h y) u = 4u m H 2, 0, u m = 1.5, Re = 100 back
31 Flow past a cylinder Lift Inflow: 0.6 ( ) y(h y) u = 4u m H 2, 0, u m = 1.5, Re = Time
32 Conclusion Second order (limited) entropy stable finite volume scheme for system of conservation laws, on unstructured meshes. Discretized viscous fluxes for Navier-Stokes equations with appropriate boundary (homogeneous) fluxes to obtain global entropy stability. Next steps: Higher order finite volume schemes? Entropy stable boundary conditions for general BC? Parellelisation and 3D implementation.
33 Thank you.
34 2D Euler Equations where Choose ρ ρu ρv U = ρu ρv, f 1(U) = ρu 2 + p ρuv, f 2(U) = ρuv ρv 2 + p E u(e + p) v(e + p) η(u) = ρs γ 1, ( ) 1 p E = ρ 2 u2 + (γ 1)ρ q 1(U) = ρus γ 1, q 2(U) = ρvs γ 1 s = ln (p) γ ln (ρ) (physical entropy) γ s γ 1 β u 2 V = 2βu, β = ρ 2p = 1 RT 2β
35 Entropy conservative flux for Euler Kinetic energy and entropy conservative flux F,ρ ρu n F = F,m1 F,m2 = pn 1 + uf,ρ pn 2 + vf,ρ F,e F,e where u n = un 1 + vn 2, p = ρ 2β, [ ] F,e 1 = 2(γ 1) β 1 2 u 2 F,ρ + u F,m The crucial property for kinetic energy preservation (Jameson) F m = pn + uf ρ for any consistent approximations of p and F ρ.
36 Entropy stable KEPES flux For Euler equations, we choose Roe type dissipation where D ij = R ij Λ ij R ij R = u añ 1 u ñ 2 u + añ 1 v añ 2 v ñ 1 v + añ 2 S H auñ 2 u 2 uñ 2 vñ 1 H + auñ [ ] ρ S = diag 2γ, (γ 1)ρ ρ, p, γ 2γ Λ = diag [u n a, u n, u n, u n + a] ñ unit face normal, uñ = u ñ γp a = speed of sound in air ρ H = a2 γ 1 + u 2 specific enthalpy 2
37 Second order (limited) reconstruction i 1 i j + 1 j The forward differences f ij = Z j Z i, f ji = Z j+1 Z j = [Z i +2(x j x i ) h Z j ] Z j The backward differences b ij = Z i Z i 1 = Z i [Z j 2(x j x i ) h Z i ], b ji = Z j Z i h Z i = R ij h V i
38 Entropy stability is important 1D modified shocktube test: (ρ L, u L, p L ) = (1.0, 0.75, 1.0), (ρ R, u R, p R ) = (0.125, 0.0, 0.1) Exact KEPES ROE Exact KEPES ROE density x velocity 0.5 Exact KEPES pressure ROE x Roe s solver (without fix) gives entropy violating shock. x x
39 Gradients n1,2 ni,2 n T 1 n T 2 i 1 T ni,1 n5,i n5, h Z i = R ij h V i n T 3 h V T = 1 T h V Te = 1 T (V i + V j ) nk T } 2 {{} trapezoidal rule [ (V b i + V b j ) 2 + (V j + V k ) 2 n e + (V j + V k ) 2 n T i + (V k + V i ) 2 n T j ] n Te i + (V k + V i ) n Te j 2
40 Gradients n1,2 ni,2 n T 1 n T 2 i 1 T ni,1 n5,i n5, h Z i = R ij h V i n T 3 h V i = T h V T T i T T i
41 Time integration: SSP-RK3: Explicit strong stability preserving RK3 du dt = L(u) v (0) = v n v (1) = v (0) + tl(v (0) v (2) = 3 4 v(0) + 1 [v (1) + tl(v (1)] 4 v (3) = 1 3 v(0) + 2 [v (2) + tl(v (2)] 3 v n+1 = v (3)
42 Step in wind tunnel Supersonic flow past forward step, M=3. Corner is the center of a rarefaction fan = singular point. Shocks undergo reflections.
43 2D Navier-Stokes Equations g 1 (U) = U t + f 1 x + f 2 y = g 1 x + g 2 y 0 τ 11 τ 21 uτ 11 + vτ 21 Q 1, g 2(U) = 0 τ 21 τ 22 uτ 21 + vτ 22 Q 2 where τ 11 = 4µ u 3 x 2µ ( v u 3 y, τ 12 = τ 21 = µ y v ), x τ 22 = 4µ 3 v y 2µ 3 u x, Q = (Q 1, Q 2 ) = κ T µ dynamic coeff. of viscosity κ coeff. of heat conductance
44 Inviscid boundary flux Boundary inviscid flux F ie is chosen as F ρ ie F ie = 1 2 (pb i n e) + u i F ρ ie, F ρ Fie E ie (ρ = ub n e 2 [ ( Fie E = ρ u γp ) ] u b [ ne + (p b p) u n e γ i ) i ] i
45 Numerical results: Viscous shock tube problem Balancing effects of physical viscosity, the heat flux and artificial numerical viscosity (ρ L, u L, p L ) = (1.0, 0.0, 1.0) (ρ R, u R, p R ) = (0.125, 0.0, 0.1) P r = 1.0, h κ = µc p P r We consider three flows regimes depending on µ. ND Numerical dissipation present Visc Physical viscosity present Heat Heat flux present
46 Numerical results: Viscous shock tube problem Regime 1: µ = Visc+Heat+ND Visc+ND Visc Visc+Heat Visc+Heat+ND Visc+ND Visc Visc+Heat density x vleocity x x Severe oscillations in the absence of ND.
47 Numerical results: Viscous shock tube problem Regime 2: µ = Visc+Heat+ND Visc+ND Visc Visc+Heat Visc+Heat+ND Visc+ND Visc Visc+Heat density x x vleocity x
48 Numerical results: Viscous shock tube problem Regime 2: µ = Visc+Heat+ND Visc+ND Visc Visc+Heat density x
49 Numerical results: Viscous shock tube problem Regime 3: µ = Visc+Heat+ND Visc+ND Visc Visc+Heat Visc+Heat+ND Visc+ND Visc Visc+Heat density x vleocity x x
50 Numerical results: Viscous shock tube problem Regime 2: µ = Visc+Heat+ND Visc+ND Visc Visc+Heat density x
51 Numerical results: Laminar flat-plate boundary layer Re = 10 5, M = 0.1 Slip BC on inlet portion Farfield BC on In-flow Adiabatic no-slip BC on flat plate Freestream pressure BC on top and outlet Inlet portion Flat plate
52 Numerical results: Laminar flat-plate boundary layer Comparing with Blasius semi-analytical solution on the vertical line through plate center 7 6 Blasius ROE KEPES-TECNO 7 6 Blasius ROE KEPES-TECNO eta eta u/uinf Streamwise velocity v*sqrt(2*rex)/uinf Vertical velocity
Finite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 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 informationarxiv: v1 [cs.na] 22 Sep 2012
Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier- Stokes equations Praveen Chandrashekar TIFR Center for Applicable Mathematics, Bangalore 5665, India
More informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationComputation of entropic measure-valued solutions for Euler equations
Computation of entropic measure-valued solutions for Euler equations Eitan Tadmor Center for Scientific Computation and Mathematical Modeling (CSCAMM) Department of Mathematics, Institute for Physical
More informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
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 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 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 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 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 Dept. of
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
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 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 informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
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 informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 5, No., pp. 544 573 c Society for Industrial and Applied Mathematics ARBITRARILY HIGH-ORDER ACCURATE ENTROPY STABLE ESSENTIALLY NONOSCILLATORY SCHEMES FOR SYSTEMS OF CONSERVATION
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationACCURATE NUMERICAL SCHEMES FOR APPROXIMATING INITIAL-BOUNDARY VALUE PROBLEMS FOR SYSTEMS OF CONSERVATION LAWS.
ACCURATE NUMERICAL SCHEMES FOR APPROXIMATING INITIAL-BOUNDARY VALUE PROBLEMS FOR SYSTEMS OF CONSERVATION LAWS. SIDDHARTHA MISHRA AND LAURA V. SPINOLO Abstract. Solutions of initial-boundary value problems
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 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 informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
More informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More information1. Introduction Some Basic Concepts
1. Introduction Some Basic Concepts 1.What is a fluid? A substance that will go on deforming in the presence of a deforming force, however small 2. What Properties Do Fluids Have? Density ( ) Pressure
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 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 informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
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 informationEntropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing
Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Siddhartha Mishra Seminar for Applied Mathematics ETH Zurich Goal Find a numerical scheme for conservation
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
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 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 informationNonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50
Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations Matthew McMullen Advisor: Antony Jameson Co-advisor: Juan Alonso Sponsor: Accelerated Strategic Computing Initiative
More informationComputation of NACA0012 Airfoil Transonic Buffet Phenomenon with Unsteady Navier-Stokes Equations
5th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 9-2 January 22, Nashville, Tennessee AIAA 22-699 Computation of NACA2 Airfoil Transonic Buffet Phenomenon with
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationRANS Solutions Using High Order Discontinuous Galerkin Methods
RANS Solutions Using High Order Discontinuous Galerkin Methods Ngoc Cuong Nguyen, Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 2139, U.S.A. We present a practical
More informationSteady waves in compressible flow
Chapter Steady waves in compressible flow. Oblique shock waves Figure. shows an oblique shock wave produced when a supersonic flow is deflected by an angle. Figure.: Flow geometry near a plane oblique
More informationBoundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity
Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with
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 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 informationAdjoint-based aerodynamic shape optimization
IT Licentiate theses 2003-012 Adjoint-based aerodynamic shape optimization OLIVIER AMOIGNON UPPSALA UNIVERSITY Department of Information Technology Adjoint-based aerodynamic shape optimization BY OLIVIER
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
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 informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationFoundations of Computational Mathematics: Hong Kong 2008
London Mathematical Society Lecture Note Series: 363 Foundations of Computational Mathematics: Hong Kong 008 Edited by Felipe Cucker City University of Hong Kong Allan Pinkus Technion Mike Todd Cornell
More informationCompressible Potential Flow: The Full Potential Equation. Copyright 2009 Narayanan Komerath
Compressible Potential Flow: The Full Potential Equation 1 Introduction Recall that for incompressible flow conditions, velocity is not large enough to cause density changes, so density is known. Thus
More informationSimulation of Condensing Compressible Flows
Simulation of Condensing Compressible Flows Maximilian Wendenburg Outline Physical Aspects Transonic Flows and Experiments Condensation Fundamentals Practical Effects Modeling and Simulation Equations,
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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
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 informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
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 informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationFundamentals of Aerodynamics
Fundamentals of Aerodynamics Fourth Edition John D. Anderson, Jr. Curator of Aerodynamics National Air and Space Museum Smithsonian Institution and Professor Emeritus University of Maryland Me Graw Hill
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
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 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 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 informationParallel Computation of Forced Vibration for A Compressor Cascade
44th AIAA Aerospace Sciences Meeting and Exhibit 9-2 January 26, Reno, Nevada AIAA 26-628 AIAA Paper 26-628 Parallel Computation of Forced Vibration for A Compressor Cascade Zongjun Hu and Gecheng Zha
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 informationModule-5: Hypersonic Boundary Layer theory. Lecture-20: Hypersonic boundary equation
Module-5: Hypersonic Boundary Layer theory Lecture-0: Hypersonic boundary equation 0.1 Governing Equations for Viscous Flows The Navier-Stokes (NS) equaadtions are the governing equations for the viscous
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 informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
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 informationSimulation of low Mach number flows
Simulation of low Mach number flows Improvement of accuracy and convergence of the TAU code using time derivative preconditioning Ralf Heinrich, Braunschweig, 9th February 008 Simulation of low Mach number
More informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
More informationFundamentals of Aerodynamits
Fundamentals of Aerodynamits Fifth Edition in SI Units John D. Anderson, Jr. Curator of Aerodynamics National Air and Space Museum Smithsonian Institution and Professor Emeritus University of Maryland
More informationHyperbolic Conservation Laws Past and Future
Hyperbolic Conservation Laws Past and Future Barbara Lee Keyfitz Fields Institute and University of Houston bkeyfitz@fields.utoronto.ca Research supported by the US Department of Energy, National Science
More informationPEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru
Governing Equations of Fluid Flow Session delivered by: M. Sivapragasam 1 Session Objectives -- At the end of this session the delegate would have understood The principle of conservation laws Different
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationAn Unstructured Grid Incompressible Navier-Stokes Algorithm for Convective Heat. Transfer Based on Artificial Compressibility. Jessica Elaine Kress
An Unstructured Grid Incompressible Navier-Stokes Algorithm for Convective Heat Transfer Based on Artificial Compressibility By Jessica Elaine Kress Approved: Lafayette Taylor Professor of Computational
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 informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationFLOW OVER A GIVEN PROFILE IN A CHANNEL WITH DYNAMICAL EFFECTS
Proceedings of Czech Japanese Seminar in Applied Mathematics 2004 August 4-7, 2004, Czech Technical University in Prague http://geraldine.fjfi.cvut.cz pp. 63 72 FLOW OVER A GIVEN PROFILE IN A CHANNEL WITH
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
More informationMasters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16
Masters in Mechanical Engineering Aerodynamics st Semester 05/6 Exam st season, 8 January 06 Name : Time : 8:30 Number: Duration : 3 hours st Part : No textbooks/notes allowed nd Part : Textbooks allowed
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 informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
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 informationHybrid RANS/LES simulations of a cavitating flow in Venturi
Hybrid RANS/LES simulations of a cavitating flow in Venturi TSFP-8, 28-31 August 2013 - Poitiers - France Jean Decaix jean.decaix@hevs.ch University of Applied Sciences, Western Switzerland Eric Goncalves
More informationAdjoint code development and optimization using automatic differentiation (AD)
Adjoint code development and optimization using automatic differentiation (AD) Praveen. C Computational and Theoretical Fluid Dynamics Division National Aerospace Laboratories Bangalore - 560 037 CTFD
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 information2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations
. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations We need to review the governing equations of fluid mechanics before examining the methods of computational aerodynamics in detail.
More informationarxiv: v1 [math.na] 30 Jan 2018
ENTROPY CONSERVATIVE SCHEMES AND THE RECEDING FLOW PROBLEM AYOUB GOUASMI, SCOTT MURMAN, AND KARTHIK DURAISAMY arxiv:1801.1013v1 [math.na] 30 Jan 018 Abstract. This work delves into the family of entropy
More informationOn divergence-free reconstruction schemes for CED and MHD
On divergence-free reconstruction schemes for CED and MHD Praveen Chandrashekar praveen@math.tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore-560065, India
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 informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More information