On a class of numerical schemes. for compressible flows
|
|
- Tobias Shaw
- 5 years ago
- Views:
Transcription
1 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. Cadarache C.E.A Cadarache
2 Some key words: Compressible Euler and Navier Stokes Implicit or semi-implicit (or explicit) schemes Internal energy balance formulation (rather than total energy formulation) even in the presence of shocks All Mach scheme Finite volume / non conforming finite element Upwinding with respect to the material velocity (no Riemann solver) Colocated or staggered grids Stability via a kinetic energy inequality Consistency in the Lax-Wendroff sense if a conservative numerical scheme for a hyperbolic system of conservation laws converges, then it converges towards a weak solution.
3 The Euler (NS) equations: total energy vs. internal energy Euler (Navier-Stokes) equations: t ϱ + div(ϱu) = 0, t (ϱu) + div(ϱu u) divτ + p = 0, t (ϱe) + div [ (ϱe + p)u ] = div(τ u), (mass) (mom) (tot.en) p = (γ 1) ϱe, E = 1 2 u 2 + e. For regular functions, (mom) u & (mass) (kin.en): 1 2 t (ϱ u 2 ) div(ϱ u 2 u) + p u = div(τ ) u. (kin.en) Subtracting from (tot.en) yields the internal energy balance: t (ϱe) + div(ϱeu) + p divu = τ : u, (int.en) which implies e 0. Incompressible schemes use the internal energy (or temperature) equation.
4 Internal energy Dealing with the internal energy: positive internal energy convenient for incompressible problems ρe is not a conserved variable conserved variables : ρ, ρu, ρe Formulations based on variables other than the conserved variables (non-conservative variables) fail at shock waves. They give the wrong jump conditions; consequently they give the wrong shock strength, the wrong shock speed and thus the wrong shock position. (Toro, 1999) density Test 5 of [ Toro chapter 4] - Density at t = 0.035, n = 2000 cells, with and without corrective source terms, and analytical solution. 5 0 with without exact x
5 Right and wrong shock speed for Bürgers Bürgers equation: for regular positive solutions (B) : t u + x (u 2 ) = 0 (BS) : t u x u3 = 0. No longer true with irregular solutions: Rankine Hugoniot gives σ = u l + u r and σ = 4 3 (u l + u r ). Weak solutions of (B) weak solutions of (BS). Figure: Explicit upwind Scheme for (B) (left) and (BS) (right) with different mesh sizes, CFL = 1.
6 Bürgers, numerical diffusion Scalar conservation law: t u + x (f (u)) = 0, f 0 Upwind scheme: h δt (un+1 i u n i ) + f (u n i ) f (u n i 1 ) = 0. Upwinding is formarly similar to add a numerical diffusion. t u + x (f (u)) x ( hf (u) δt f 2 (u) x u) = 0 2 The CFL condition states that hf (u) δt f 2 (u) 0 (i.e. δt f (u) h) In the case of the Bürgers equation it gives t u + x (u 2 ) x ((hu 2δt u 2 ) x u) = 0, x R, t R +
7 Bürgers, non conservative numerical diffusion In the case of the equivalent equation t (u 2 ) + (4/3) x (u 3 ) = 0, u > 0 upwinding is formally consistent with Dividing by 2u, this leads to t (u 2 ) x (u3 ) x ((2hu 2 4δt u 3 ) x u) = 0, t u + x (u 2 ) 1 u x ((hu2 2δt u 3 ) x u) = 0, x R, t R + This is a numerical diffusion (thanks to the CFL condition) but not in conservative form. A non conservative diffusion may lead to wrong discontinuities
8 Benefit from a non conservative numerical diffusion? t u + x (u 2 ) = 0 + initial condition (B) t (u 2 ) x (u3 ) = 0 + initial condition (BS) Explicit upwind scheme on (BS) formally equivalent to: t u + x (u 2 ) 1 u x ((hu2 2δt u 3 ) x u) = 0. }{{} non conservative numerical diffusion. Negative result for a non conservative diffusion :-( Non conservative numerical diffusion on (B) yields { wrong shock velocity for (B) correct shock velocity for (BS) Positive result for a non conservative diffusion? :-) Non conservative numerical diffusion on (BS) yields { wrong shock velocity for (BS) correct shock velocity for (B)? How do we choose the non conservative numerical diffusion?
9 Non conservative numerical diffusion on (BS) Viscous Bürgers: Multiplying by 2u: t u + x (u 2 ) ε xx u = 0. (B) ε t (u 2 ) x (u3 ) 2uε xx u = 0. (BS) ε Discretize (BS) ε instead of (BS): t (u 2 ) x (u3 ) uε 0 h xx u = 0, (BS) ε with 2ε = ε 0 h.
10 Non conservative numerical diffusion on (BS) t (u 2 ) x (u3 ) uε 0 h xx u = 0, (BS) ε with 2ε = ε 0 h. Centered finite volume with non conservative diffusion ( (n)) 2 ( (n 1)) 2 4δt u i = u i + u(n 1) + u (n 1) 3 i 1 i 3h 2 + δt [ h 2 ε 0h u (n 1) i u(n 1) i u (n 1) i u (n 1) i+1 2 2u (n 1) i ] + u (n 1). i+1 Centered Scheme for (BS) ε0 h, ] uε 0 h xx u = ε 0 h [ x (u x u) + x ( u2 2 )
11 From Bürgers to Euler For regular solutions, Bürgers: t u + x (u 2 ) = 0 t (u 2 ) x (u3 ) = 0. Euler: t ϱ + div(ϱu) = 0, t (ϱu) + div(ϱu u) + p = 0, t (ϱe) + div [ (ϱe + p)u ] = 0. t ϱ + div(ϱu) = 0, t (ϱu) + div(ϱu u) + p = 0, t (ϱe) + div(ϱeu) + p divu = 0. we had an equation, we now have a system... Idea: add a non conservative corrective term to the internal energy equation. Which term? Inspiration comes from copying the formal derivation of the internal energy equation at the discrete level.
12 Total energy= kinetic energy + internal energy Total energy (Euler): t (ϱe) + div [ (ϱe + p)u ] = 0, E = e u 2 = t (ρe) + div(ρeu) + p divu t (ρ u 2 ) div(ρ u 2 u) + u p = 0. From mass balance, for regular z: t (ρz) + div(ρzu) =( t ϱ + div(ϱu)) }{{} z + ρ t z + ρu z, = 0. = 1 2 t (ρui 2 ) + 1 [ ] 2 div(ρu2 i u) = ρu i t u i + ρu i u u i = u i ρ t u i + ρu u i = [ u i t (ρu i ) + div(ρu i u) ] = u i i p, 1 i 3. From momentum balance: = 1 2 t (ρ u 2 ) div(ρ u 2 u) = u [ t (ρu) + div(ρu u) ] = u p. and t (ρe) + div(ρeu) + p div(u) = 0.
13 Strategy for a numerical scheme for Euler How to get correct weak solutions of Euler equations while solving the internal energy balance? More precisely: General idea: (int.en) d + (kin.en) d (tot.en) d 1- (mom) d and (mass) d discrete kinetic energy (E k ) balance 2- (int.en) d : ð t E k + div d (E k u) + u d p + R = 0 R : non conservative residual term ð t (ρe) + div d (eu) + pdiv d u = R R 0 3- Perform Lax-Wendroff consistency analysis of the scheme: (a) Suppose bounds and convergence for a sequence of discrete solutions control in BV and L, convergence in L p, for p 1. (b) Let ϕ be a regular function, interpolate, test the kinetic energy balance, test the internal energy balance, and pass to the limit in the scheme. The corrective term in the internal energy balance is such that, at the limit, the weak form of the total energy equation is recovered.
14 Required discrete properties Discrete transport property, i.e. discrete equivalent of t (ρz) + div(ρzu) = ρ t z + ρu z, z = u i. = Compatible discretization of mass and momentum balance equation Discrete duality i.e. discrete equivalent of div(pu) = pdivu + u p. Positivity of the residual R 0 in the discrete kinetic energy balance equation (to ensure the positivity of the internal energy). Points to be taken care of when designing the scheme(s).
15 Time discretization : I - Decoupled (explicit) choice Decoupled 1: natural ordering, bad idea..., ρ u e p 1 δt (ϱn+1 ϱ n ) + div(ϱ n u n) ) = 0, (mass) d ϱ n+1 1 δt (ϱn+1 u n+1 ϱ n u n ) + div(ϱ n u n u n ) + p n = 0, (mom) d u n+1. 1 δt (ϱn+1 e n+1 ϱ n e n ) + div(ϱ n e n u n ) + p n divu n+1 = R n, (int.en) d e n+1 p n+1 = (ϱ n+1, e n+1 ) = (γ 1) ϱ n+1 e n+1, (eos) d p n+1 Decoupled 2: better idea..., ρ e p u 1 δt (ϱn+1 ϱ n ) + div(ϱ n u n) ) = 0, (mass) d ϱ n+1 1 δt (ϱn+1 e n+1 ϱ n e n ) + div(ϱ n e n u n ) + p n divu n = R n, (int.en) d e n+1 p n+1 = (ϱ n+1, e n+1 ) = (γ 1) ϱ n+1 e n+1, (eos) d p n+1 1 δt (ϱn+1 u n+1 ϱ n u n ) + div(ϱ n u n u n ) + p n+1 = 0, (mom) d u n+1.
16 Time discretization : I -Decoupled choice 1000 pressure exact 1000 pressure exact pressure pressure ρ u e p x ρ e p u x [E. Toro, Riemann solvers and numerical methods for fluid dynamics, third edition, test 3 of chapter 4].
17 Decoupled scheme : Why choose ρ e p u? Derivation of the entropy balance at the continuous level. For regular functions, an entropy s satisfies t s + u s = 0 or equivalently t (ρs) + div(ρus) = 0; for Euler perfect gas s = ψ(ρ, e) = ln(ρ) 1 γ 1 lne. t ϱ + u ϱ + ϱ divu = 0 ϱψ(ρ, e) t e + u e + p divu = 0 eψ(ρ, e) ϱ [ t s + u s + ϱ ϱs + p ] ϱ es divu = 0. }{{} =0 We need div(ϱu) and p divu at the same time level to mimick this computation at the discrete level. Choice 1 ρ u e p div(ϱ n u n ) and p n divu n+1 Choice 2 ρ e p u div(ϱ n u n ) and p n divu n
18 Time discretization : II - Implicit or semi-implicit choice Implicit ρ n+1 ρ n + div(ρ n+1 u n+1 ) = 0, δt (mass) d 1 δt (ϱn+1 u n+1 ϱ n u n ) + div(ϱ n+1 u n+1 u n+1 ) + p n+1 = 0, (mom) d 1 δt (ϱn+1 e n+1 ϱ n e n ) + div(ϱ n+1 e n+1 u n+1 ) + p n+1 divu n+1 = R n+1, (int.en) d p n+1 = (ϱ n+1, e n+1 ) = (γ 1) ϱ n+1 e n+1, Semi-implicit ( ρ Pressure gradient scaling step: ( p) n+1 n ) 1/2 = ( p n ρ n 1 ) (eos) d Prediction step Solve for ũ n+1 : 1 ( ρnũ n+1 ρ n 1 u n) + div(ρ n ũ n+1 u n ) + ( p) n+1 = 0, (mom) δt d Correction step Solve for p n+1, e n+1, ρ n+1 and u n+1 : 1 δt ρn (u n+1 ũ n+1 ) + ( p n+1 ) ( p) n+1 = 0, ρ n+1 ρ n + div(ρ n+1 u n+1 ) = 0, (mass) δt d 1 δt (ρn+1 e n+1 ρ n e n ) + div(ρ n+1 e n+1 u n+1 ) + p n+1 (div(u n+1 )) = R n+1, (int.en) d ρ n+1 = ϱ ( e n+1, p n+1). (eos)
19 Meshes K p K, ρ K, e K, u K Colocated Advantage: Easier Data structure Easily refined Pressure correction scheme studied for the Euler equations (C. Zaza s thesis). u σw u σn K p K, ρ K, e K u σs u σe u (2) σ N K u σ (1) W u σ (1) p K, ρ K, e E K Staggered: Crouzeix-Raviart (on simplices) Rannacher-Turek (on quadrangles) full velocities on the edges (faces) MAC: normal velocities on the edges (faces) Advantage: Stable in the incompressible limit u (2) σ S
20 Space discretization: Finite volume discretization of the mass equation K ρ n+1 ρ n + div(ρ n+1 u n+1 ) = 0, (mass) δt (mass) discretization of the fluxes: K δt K (ρ n+1 K ρ n K ) + ρ n+1 ρ n σ E(K ) δt F n+1 = σ ρn+1 K,σ σ un+1 σ n K,σ, numerical flux through σ. F n+1 K,σ = 0, + (ρ n+1 u n+1 n K ) = 0. K σ F ρ K,σ K u σ K ρ n+1 σ upwind approximation of ρ n+1 at the face σ with respect to u n+1 σ n K,σ. Positive density: ρ n+1 > 0 if (ρ n > 0 and ρ > 0 at inflow boundary) σ
21 Discrete internal energy equation and E.O.S. 1 δt (ρn+1 e n+1 ρ n e n ) + div(ρ n+1 e n+1 u n+1 ) + p n+1 (div(u n+1 )) = R n+1 Discretization by upwind finite volume of the discrete internal energy K δt (ρn+1 e n+1 ρ n K K K en K ) + σ E(K ) eσ n+1 upwind choice positivity of e (if R n+1 0) K K (divu) K = σ u σ n K,σ. σ E(K ) F n+1 K,σ en+1 σ + K p n+1 K (divu n+1 ) K = R n+1 K, discrete E.O.S. p n+1 K = (γ 1)ρ n+1 K e n+1 K. (eos) d
22 Discretization of the momentum equation 1 ( ρnũ n+1 ρ n 1 u n) + div(ρ n ũ n+1 u n ) + ( p) n+1 (mom) n δt (mom) n D σ + ρ n ũ n+1 ρ n 1 u n ρ n ũ n+1 u n n K D σ δt D σ }{{} C(ρ n, u n ) + D σ ( p) n+1 = 0. Space discretization D σ (ρ n D δt σ ũ n+1 σ ρ n 1 u n D σ σ) + Fσ,ɛ n ũn+1 ɛ + D σ ρn D σ ρ n 1 ( p n ) σ = 0, ɛ E(D σ) D σ }{{} C d (ρ n, u n ) Grad-div duality : K p K (divu) K + D σ u σ ( p) σ = 0 K T σ E D σ ( p n ) σ = σ (pl n pn K ) n K,σ for σ = K L. ρ n D σ? F n σ,ɛ? K ρ Dσ F K,σ L u σ D σ σ F σ,ɛ
23 Discretization of the convection operator Choice of ρ n D σ, ρ n 1 D σ and F n σ,ɛ in C d (ρ n, u n, ũ n+1 )? Discretize C d (ρ n, u n, ũ n+1 ) so as to obtain a discrete kinetic energy balance. Copy the continuous kinetic energy balance: (mom) u & (mass) (kin.en) with some formal algebra... using t ρ + div(ρ u) = 0. Do the same at the discrete level? Momentum on dual cells, mass on primal cells... Idea: reconstruct a mass balance on the the dual cells Choose ρ Dσ = 1 ( ) DK,σ ρ K + D L,σ ρ L Dσ F σ,ɛ : linear combination of the primal fluxes (F K,σ ) σ E(K ). so that a discrete mass balance holds on the dual cells D σ: σ E int, D σ δt (ρ n+1 D σ ρ n D σ ) + ɛ E(D σ) F n+1 σ,ɛ = 0. Then take C d (ρ n, u n, ũ n+1 ) = Dσ δt with u n+1 ɛ (ρ n D σ ũ n+1 σ ρd Fσ,ɛ n ũn+1 ɛ n 1 σ u n σ) + ɛ E(D σ) = 1 2 (un+1 σ + u n+1 σ )
24 Discrete kinetic energy balance: computation of R σ Continuous setting: Multiply continuous momentum by u: ( ) t (ρ u) + div(ρ u u) + p = 0 u... with some formal algebra... using t ρ + div(ρ u) = 0, continuous kinetic energy balance: t ( 1 2 ρ u 2 ) + div ( ( 1 2 ρ u 2 ) u ) + p u = 0 (kin.en) Discrete setting: Similarly, multiply discrete momentum by ũ n+1 σ : ( Dσ (ρ n D δt σ ũ n+1 σ ρ n 1 σ u n σ) + ) F σ,ɛũn+1 n ɛ + D σ ( p n+1 ) σ = 0 ũ n+1 σ... with some algebra... using ɛ E(D σ) Mass on dual cell: Dσ (ρ n δt Dσ ρn 1 σ ) + F n σ,ɛ = 0. ɛ E(Dσ ) Scaling of the pressure ρ n 1 ( p) n+1 2 = ρ n p n 2. Correction equation 1 δt ρn u n+1 + p n+1 = 1 δt ρn ũ n+1 + ( p) n+1 discrete kinetic energy balance: 1 2 D [ σ ρ n δt σ u n+1 σ 2 ρ n 1 D σ u n σ 2] D σ ( p n+1 ) σ u n+1 σ ɛ=d σ D σ Fσ,ɛ n ũn+1 σ ũ n+1 σ +Rn+1 σ +Pn+1 σ = 0 with Rσ n+1 0, (kin.en) σ
25 Choice of R K R σ = Dσ ρ n 1 δt functions. ũ n+1 D σ σ u n σ 2 0 for regular functions, but NOT for discontinuous By definition of ρ Dσ, for σ = K L, R σ = D K,σ ρ n 1 δt ũ n+1 K σ u n σ 2 + D L,σ δt ρ n 1 ũ n+1 L σ u n 2 σ R K = 1 2 D K,σ ρ n 1 ũ n+1 δt K σ σ E(K ) u n 2 σ R K R σ = 0 K T σ E
26 Passage to the limit: total energy recovered Kinetic energy (kin) n σ = Dσ 2δt (ϱn+1 σ un+1 σ 2 ϱ n σ un σ 2 ) + 1 F n 2 σ,ɛũn+1 σ ũn+1 σ + ( p) n+1 σ un+1 σ = R n+1 σ + Pn+1 σ, ɛ=dσ D σ Internal energy (int) n K = K δt (ρn+1 K e n+1 K ρ n K en K ) + σ E(K ) F n K,σ en σ + K pn+1 K (divu n+1 ) K = R n+1 K, ϕ : test function Multiply (kin) σ by interpolate ϕ n σ and (int) K by interpolate ϕ n K n σ E(kin) n σ ϕn σ + (int) n K ϕn K = δtr K ϕ n K δtr σϕ n σ + δtp σϕ n σ. n K T n K T n σ E n σ E }{{}}{{}}{{} T [ ] ρ E t ϕ + (ρ E + p) u ϕ Ω ρ 0 (x) E 0 (x) ϕ(x, 0) Ω In particular, the pressure terms combine themselves to converge to pu ϕ.
27 Numerical tests - I Euler, high Mach ρ p u (1) u (2) Mach 3 facing step (Woodward Collela) MAC space discretization, uniform grid, δt = h/4 = 0.001, (u1 + c = 4 at the inlet boundary). Artificial viscosity µ = 0.001, which roughly corresponds to the numerical viscosity associated to an upwinding of the convection term µupw ' ρ u h/2 divided by 5.
28 Numerical tests - II Flow past cylinder, low Mach Flow past a cylinder, benchmark Schäfer and S. Turek, Mach 0.003, Re 100. Pressure correction scheme, Rannacher-Turek FE. coarse mesh fine mesh Mesh Space unks c d,max c l,max St m m m m Reference range Table: Drag and lift coefficients and Strouhal number.
29 Numerical tests - II Flow past cylinder, high Mach Flow past a cylinder, Mach ' 3, Re ' 100. pext = γ /10 ρ, c = 0.1. mes = impermeability and and perfect slip condition at the upper and lower boundaries u = (1, 0)t. at inlet. δt = Rannacher-Turek FE. ρ p u (1) u (2) t = 5, mesh of 106
30 Summary: main features of the class of schemes Colocated or staggered discretization Upwind choice for ρ, positivity of the density, Compatible discretization of (mass) and (mom) & careful choice of ρ Dσ and fluxes in (mom) to recover a mass conservation on the dual cells discrete kinetic energy inequality Compatible discretization of (mass) and (int.en) & upwind choice for e K under CFL positivity of e Conservation of total mass Existence of a solution to the scheme (topological degree argument) Velocity and pressure are constant through the contact discontinuity: p n and u n constant p n+1 and u n+1 constant. Consistency : under compactness assumptions, the discrete solution tends to a weak solution of the Euler systems. Numerical discontinuous solutions have correct shocks.
31 Recent and on going work Higher order schemes for the explicit MUSCL and additional viscosity. Convergence of the scheme to an entropy weak solution Low Mach limit for barotropic NS (with K. Saleh, M.-H. Vignal). Convergence for barotropic NS (γ 3) and strong-weak error estimates for barotropic NS (with D. Maltese, and A. Novotny). Reactive flows. Analysis of the MAC scheme in primitive variables for incompressible flows (K. Mallem s thesis)
Convergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More informationA second order scheme on staggered grids for detonation in gases
A second order scheme on staggered grids for detonation in gases Chady Zaza IRSN PSN-RES/SA2I/LIE August 1, 2016 HYP 2016, RWTH Aachen Joint work with Raphaèle Herbin, Jean-Claude Latché and Nicolas Therme
More informationAn adaptive projection scheme for all-mach flows applied to shock diffraction
An adaptive projection scheme for all-mach flows applied to shock diffraction CEA DEN/DM2S/STMF/LMEC AMU I2M UMR 7373/AA Group January 20, 2015 All-Mach solver? Approximate Riemann solvers for M 1 (Guillard-Viozat,
More informationFinite difference MAC scheme for compressible Navier-Stokes equations
Finite difference MAC scheme for compressible Navier-Stokes equations Bangwei She with R. Hošek, H. Mizerova Workshop on Mathematical Fluid Dynamics, Bad Boll May. 07 May. 11, 2018 Compressible barotropic
More informationConsistency analysis of a 1D Finite Volume scheme for barotropic Euler models
Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models F Berthelin 1,3, T Goudon 1,3, and S Mineaud,3 1 Inria, Sophia Antipolis Méditerranée Research Centre, Proect COFFEE Inria,
More informationEntropy estimates for a class of schemes for the euler equations
Entropy estimates for a class of schemes for the euler equations Thierry Gallouët, Raphaele Herbin, J.-C Latché, N Therme To cite this version: Thierry Gallouët, Raphaele Herbin, J.-C Latché, N Therme.
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationarxiv: v1 [math.na] 16 Oct 2007
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher arxiv:070.2987v [math.na] 6 Oct 2007 AN UNCONDITIONNALLY STABLE PRESSURE CORRECTION
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 informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationarxiv: v1 [math.na] 17 Jan 2017
CONVERGENCE OF THE MARKER-AND-CELL SCHEME FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON NON-UNIFORM GRIDS T. Gallouët 1, R. Herbin 2, J.-C. Latché 3 and K. Mallem 4 arxiv:1701.04553v1 [math.na] 17
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 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 informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
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 informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
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 informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
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 information.u= 0 ρ( u t +(u. )u)= ρ g p+.[µ( u+ t u)]
THETIS is a numerical simulation tool developed by University of Bordeaux. It is a versatile code to solve different problems: fluid flows, heat transfers, scalar transports or porous mediums. The potential
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 informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
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 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 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 informationModeling and simulation of bedload transport with viscous effects
Introduction Modeling and simulation of bedload transport with viscous effects E. Audusse, L. Boittin, M. Parisot, J. Sainte-Marie Project-team ANGE, Inria; CEREMA; LJLL, UPMC Université Paris VI; UMR
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 informationQUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES
1 QUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES O. ALLAIN(*), A. DERVIEUX(**), I. ABALAKIN(***), S. CAMARRI(****), H. GUILLARD(**), B. KOOBUS(*****),
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 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 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 informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationIntroduction to numerical simulation of fluid flows
Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural
More informationCFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki
CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing Lakhdar Remaki Outline What we doing in CFD? CFD in Industry Shock-filter model for mesh adaptation
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More 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 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 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 informationarxiv: v1 [physics.comp-ph] 10 Aug 2015
Numerical experiments on the efficiency of local grid refinement based on truncation error estimates Alexandros Syrakos a,, Georgios Efthimiou a, John G. Bartzis a, Apostolos Goulas b arxiv:1508.02345v1
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
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 informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
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 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 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 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 informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
16th International Conference on Finite Elements in Flow Problems DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology
More informationNumerical investigation of shock wave dense particles cloud interaction
25 th ICDERS August 2 7, 215 Leeds, UK Numerical investigation of shock wave dense particles cloud interaction Pavel S. Utkin Moscow Institute of Physics and Technology Dolgoprudny, Moscow region, Russia
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
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 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 informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationStudy of a numerical scheme for miscible two-phase flow in porous media. R. Eymard, Université Paris-Est
Study of a numerical scheme for miscible two-phase flow in porous media R. Eymard, Université Paris-Est in collaboration with T. Gallouët and R. Herbin, Aix-Marseille Université C. Guichard, Université
More informationHervé Guillard, INRIA Projet Smash, B.P. 93, Sophia-Antipolis Cedex, France,
TRAVELING WAVE ANALYSIS OF TWO-PHASE DISSIPATIVE MODELS Hervé Guillard, INRIA Projet Smash, B.P. 93, 06902 Sophia-Antipolis Cedex, France, Herve.Guillard@sophia.inria.fr Joint work with : Mathieu Labois,
More informationHigh-order accurate Lagrange-remap hydrodynamic schemes on staggered Cartesian grids
High-order accurate Lagrange-remap hydrodynamic schemes on staggered Cartesian grids Gautier Dakin 1, Hervé Jourdren 1 Université Pierre et Marie Curie, LJLL, Paris, France CEA, DAM, DIF, F-9197 Arpajon,
More informationOn two fractional step finite volume and finite element schemes for reactive low Mach number flows
Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes
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 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 informationRobust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems
Robust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems M. Aaqib Afaq Institute for Applied Mathematics and Numerics (LSIII) TU Dortmund 13 July 2017
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationFinite volumes schemes preserving the low Mach number limit for the Euler system
Finite volumes schemes preserving the low Mach number limit for the Euler system M.-H. Vignal Low Velocity Flows, Paris, Nov. 205 Giacomo Dimarco, Univ. de Ferrara, Italie Raphael Loubere, IMT, CNRS, France
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More 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 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationA Higher-Order Generalized Ghost Fluid Method for the Poor for the Three-Dimensional Two-Phase Flow Computation of Underwater Implosions
A Higher-Order Generalized Ghost Fluid Method for the Poor for the Three-Dimensional Two-Phase Flow Computation of Underwater Implosions Charbel Farhat, Arthur Rallu, Sriram Shankaran Department of Mechanical
More informationFinite 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 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 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 informationThe Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics
The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
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 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 informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
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 informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
More 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 informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationIntroduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués
Introduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués Contents Preliminary considerations 3. Study of the movement of a fluid................ 3.
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
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 information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
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 informationStability and consistency of kinetic upwinding for advection diffusion equations
IMA Journal of Numerical Analysis (006 6, 686 7 doi:0.093/imanum/drl005 Advance Access publication on March, 006 Stability and consistency of kinetic upwinding for advection diffusion equations MANUEL
More informationA variational approach to design a numerical scheme on an arbitrary moving grid for N-fluid flow with thermodynamic consistency
A variational approach to design a numerical scheme on an arbitrary moving grid for N-fluid flow with thermodynamic consistency T. VAZQUEZ-GONZALEZ a, A. LLOR b, C. FOCHESATO c a. CEA, DAM, DIF, 91297
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
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 information