QUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES
|
|
- Myron Mathews
- 6 years ago
- Views:
Transcription
1 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(*****), T. KOZUBSKAYA(***), A.-C. LESAGE(*), A. MURRONE(**), M.-V. SALVETTI(****), S. WORNOM(**) (*)Societe Lemma, La Roquette sur Siagne, France, (**)INRIA, route des Lucioles, Sophia-Antipolis, France, (***)Institute for athematical Modeling, Moscou, Russie, (****)U. de Pise, Italie, (*****)U. de Montpellier IMFT, Toulouse, 9 mars
2 2 Overview Ingredients of the Mixed Element Volume (MEV) method: Finite Element: Finite Volume: P1 exactness, H1 consistency. conservations, positivity. IMFT, Toulouse, 9 mars
3 3 Outline of the talk 1. Basic options of the numerical method 2. Multiple levels/scales 3. Conservations 4. Positivity IMFT, Toulouse, 9 mars
4 4 1. BASIC OPTIONS OF THE MEV METHOD degrees of freedom at vertices i median dual cells C i variational or integral formulation with two types of test functions: φ i, continuous-p 1, and χ i, characteristic of C i. Finite-Element evaluation of second derivatives Finite-Element/Finite-Volume evaluation of first derivatives IMFT, Toulouse, 9 mars
5 5 First-order MEV W t + divf = divr Vol(i) n+1 (Wi Wi n ) = Σ j N(i) Φ ij t R(W ). φ i dx Upwind numerical integration (Roe): Φ ij = 1 2 (F(W i).n ij + F(W j ).n ij ) A (W j W i ) A = d dw F(W ).n ij IMFT, Toulouse, 9 mars
6 6 First extension: MUSCL Edge-based reconstructions using upwind elements of each edge. ( W ) ij = 1 W 3 F EM (T ij ) + 2 W 3 F EM (T ijk ) W ij = W i + ( W ) ij. ij Φ ij = 0.5(Φ(W ij ) + Φ(W ji )) δ A (W ji (W ij ). IMFT, Toulouse, 9 mars
7 7 Upwind-element reconstruction: 3D Edge-upwind tetrahedra IMFT, Toulouse, 9 mars
8 8 Arbitrary Lagrangian-Eulerian formulation t (V (x,t)w(t)) + F c (w(t),x,ẋ) = R(w(t),x) Ω n+1 i W n+1 i = Ω n i W n i t j V (i) Ω ij Φ ( W n+1 i,w n+1 j, ν ij, xn+1 ij t x n ij ) Φ(W,W,ν,κ) = F (W ) ν x + G(W ) ν y + H(W ) ν z κw. IMFT, Toulouse, 9 mars
9 9 Back to the fluid numerics : V6 ( W ) V6 ij. ij = 1 3 ( W ) Tij. ij (W j W i ) + ξ a ( ( W ) Tij. ij 2(W j W i ) + ( W ) Tji. ij) + ξ b ( ( W ) D ij. ij 2 ( W ) i. ij + ( W ) j. ij ) ( W ) D ij : linear interpolation of nodal gradients in nodes m and n. - Acoustics: (Abalakin-Dervieux-Kozubskaya, IJ Acoustics, 3:2, ,2004) IMFT, Toulouse, 9 mars
10 10 V6 : Tam s case for acoustics IMFT, Toulouse, 9 mars
11 11 2. MULTIPLES LEVELS/SCALES Agglomeration: several neighboring cells a coarse one. inconsistency for second-order problems. Agglomeration MG: corrector terms in the assembly of coarse diffusion terms Additive Multi-Level: smoothed basis functions Variational Multi-Scale (VMS) models IMFT, Toulouse, 9 mars
12 12 MEV implementation of Navier-Stokes + + Ω Ω Ω ρ t X i dω + SuppX i ρu.n dγ = 0 ρu t X i dω + ρ uu.n dγ + SuppX i σ φ i dω = 0 Ω E t X i dω + (E + P )u.n dγ + SuppX i λ T. φ i dω = 0 Ω SuppX i P n dγ Ω σu. φ i dω IMFT, Toulouse, 9 mars
13 13 VMS projection C M m(k) : macro-cell containing the cell C k. (V ol(c j ): volume of cell C j ) I k = {j / C j C M m(k) }. φ k = V ol(c k) φ j. V ol(c j ) j I k j I k W = k φ k W k ; W = W W IMFT, Toulouse, 9 mars
14 14 Smagorinski Large Eddy Simulation model S ij = 1 2 ( u i x j + u j x i ), u = (u 1,u 2,u 3), small scale velocity, S = 2S ij S ij, C s = 0.1, l = V ol(t l) 1 3,. µ t = ρ(c s ) 2 S, τ ij = µ t(2s ij 2 3 S kk δ ij). IMFT, Toulouse, 9 mars
15 15 MEV implementation of the VMS LES model (ended) + + Ω Ω Ω ρ t X i dω + SuppX i ρu.nx i dγ = 0 ρu t X i dω + ρ u u nx i dγ + SuppX i σ φ i dω + τ φ i dω = 0 Ω Ω E t X i dω + (E + P )u.nx i dγ + SuppX i C p µ t λ T. φ i dω + T. φ i dω = 0 P r t Ω Ω SuppX i P nx i dγ Ω σu. φ i dω IMFT, Toulouse, 9 mars
16 16 Results: flow past a square cylinder Reynolds=20,000. From Farhat-Koobus, IMFT, Toulouse, 9 mars
17 17 Results: flow past a square cylinder, end d LES C d C d C l S t l r C pb Classical LES VMS-LES Experiments C d C d C l S t l r C pb Lyn et al Luo et al Bearman et al Tab. 1 Bulk coefficients (*): LES simulations and experimental data (*) l r is the recirculation length behind the square cylinder, and C pb is the mean pressure coefficient on rear face at y = 0. IMFT, Toulouse, 9 mars
18 18 AND NEXT... IMFT, Toulouse, 9 mars
19 19 3. CONSERVATIONS IN ALE FORMULATION - Conservation of extensive quantities. - Geometric Conservation Law (GCL). - Energy budget. IMFT, Toulouse, 9 mars
20 20 Geometric Conservation law A ALE-GCL scheme computes exactly a uniform flow field t Ω n+1 i U n+1 i = Ω n i U n i j V (i) Ω ij Φ ( U n+1 i,u n+1 j, ν ij, xn+1 ij t x n ij ν ij = 0.5 (ν ij (x(t 1 + α 1 (t 2 t 1 ))) + ν ij (x(t 1 + α 2 (t 2 t 1 )))) Ω n+1 i Ω n i = ẋ i n i dγ Ω h (t) - Sufficient condition for 1st-order accuracy (Guillard-Farhat). - Practical stability and accuracy improvements. ) IMFT, Toulouse, 9 mars
21 21 Energy budget of a deforming fluid (1) Work transfers : pressure work in moment equation: M t 2. (x n+1 ij x n ij) = t 1 t ( ) Ω h,i Φ M Ω W n+1 i, ν ij, xn+1 ij x n ij. (x n+1 ij t i Ω h W ork t 2 = t Ω h,i p i ν i. (x n+1 ij x n ij) t 1 i Ω h x n ij) IMFT, Toulouse, 9 mars
22 22 Energy budget of a deforming fluid (2) The energy equation must satisfy the Geometric Conservation Law and this is obtained by an adhoc time integration: Ω n i E n i t Ω n+1 i E n+1 i = j V (i) Ω ij and ν ij specified by GCL. ( Ω ij Φ E U n+1 i,u n+1 j, ν ij, xn+1 ij t x n ij ) IMFT, Toulouse, 9 mars
23 23 Energy budget of a deforming fluid (3) Total energy variation: E 2 t 1 = t Ω h,i Φ E Ω i Ω h Slip condition: E 2 t 1 = t Ω h,i i Ω h ( ( W n+1 i, ν ij, xn+1 ij t x n ij Ω h,i (pu) i. ν i dγ ) ) IMFT, Toulouse, 9 mars
24 24 Energy budget of a deforming fluid (4) A way of integrating the above is: E t 2 = t Ω h,i p i ν i. (x n+1 ij t 1 i Ω h x n ij) Lemma: By replacing the energy flux by a product of boundary pressure times the GCL-preserving integration of mesh motion, we can derive a scheme that is conservative, satisfies GCL and have an exact energy budget: Work of pressure = Loss of total energy. Vàzquez-Koobus-Dervieux-Farhat, Spatial discretization issues for the energy conservation in compressible flow problems on moving grids, INRIA-RR4742 IMFT, Toulouse, 9 mars
25 25 4. POSITIVITY t (V (x,t)u(t)) + F c (w(t),x,ẋ) = 0 a n+1 i U n+1 i a n i U i n t + j N (i) φ(u n ij,u n ji,ν ij,κ ij ) = 0. Φ(u,u,ν,κ) = F (u) ν x + G(u) ν y + H(u) ν z κu. with φ(u,v,ν,κ) monotone: φ u 0 and φ v 0. IMFT, Toulouse, 9 mars
26 26 Limiters U ij = W F EM (T ij ). ij ; 0 U ij = W F EM (T ijk ). ij ( U) ij. ij = L( U ij, 0 U ij, higher U ij ) From usual conditions on limiteurs, and from the position of upwind-elements : U ij = Σp k (U i U k ) ; U ij = p j (U j U i ) where all p s are positive. IMFT, Toulouse, 9 mars
27 27 POSITIVITY, cont d a n+1 i U n+1 i a n i U i n t + j N (i) φ ij (U n ij,u n ji, ν ij, κ ij ) = 0. a n+1 i U n+1 i a n i U i n t = e i + j N (i) g ij (U n ji U n i ) + h ij (U n ij U n i ) IMFT, Toulouse, 9 mars
28 28 POSITIVITY, cont d g ij = t a i φ ij (U n ij,u n ji,ν ij, κ ij ) φ ij (U n ij,u n i,ν ij, κ ij ) U n ji U n i h ij = t a i φ ij (U n ij,u n i,ν ij, κ ij ) φ ij (U n i,u n i,ν ij, κ ij ) U n ij U n i e i = an i an+1 i a n i + a n+1 i j N (i) Ω ij κ ij IMFT, Toulouse, 9 mars
29 29 POSITIVITY, cont d Lemma: Assuming the discrete GCL condition and usual limiter properties, under a CFL condition, the above scheme satisfies the maximum/minimum principle. Cournde-Koobus-Dervieux, soumis Revue EEF, 2005 Farhat-Geuzaine-Grandmont, JCP 2001 IMFT, Toulouse, 9 mars
30 30 POSITIVITY, cont d F (U) = U = (ρ,ρu,ρv,ρw,e,ρy ), ρu ρu 2 + P ρuv ρuw (e + P )u ρuy G(U) = P = (γ 1) U F (U) + + G(U) + H(U) t x y z ρv ρuv ρv 2 + P ρvw H(U) = (e + P )v ρvy ( e 1 2 (u2 + v 2 + w 2 ) ) = 0 ρw ρuw ρvw ρw 2 + P (e + P )w ρwy Lemma: Assume that the Riemann solver is ρ-positive, then under a CFL condition, the limited ALE-MEV scheme preserves ρ-positivity and satisfies the maximum/minimum principle for Y = ρy/ρ. IMFT, Toulouse, 9 mars
31 31 Cournde-Koobus-Dervieux, soumis Revue EEF, 2005 IMFT, Toulouse, 9 mars
32 32 Application to a two-phase flow Five-equation quasi conservative reduced model t α kρ k + div(α k ρ k u) = 0 ρu + div(ρu u) + p = 0 t ρe + div(ρe + p)u = 0 t t α 2 + u. α 2 = α 1 α 2 ρ 1 a 2 1 ρ 2 a k=1 α k ρ ka 2 kdivu with e = ε + u 2 /2 and ρε = 2 k=1 α kρ k ε k (p,ρ k ). IMFT, Toulouse, 9 mars
33 33 NUMERICS FOR TWO-PHASE FLOW - Acoustic Riemann solver (Murrone-Guillard, Comp.Fluids 2003, JCP 2004)) - Second order limited MEV. - Source term: upwind integration of divu. - 3D parallel extension with MEV (Wornom et al.) IMFT, Toulouse, 9 mars
34 34 Falling drop (2D) IMFT, Toulouse, 9 mars
35 35 Shock on a drop (3D) IMFT, Toulouse, 9 mars
36 36 Conclusions The choice of a very simple scheme permits the proof of a very large set of properties: - many conservation statements can be made exactly satisfied, - Maximum principle and positivity, including the moving mesh case. Today s challenge are: - showing positivity statements for the new techniques of thick interfaces. - showing the above properties for new schemes. IMFT, Toulouse, 9 mars
Extension to moving grids
Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical
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 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 informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationA 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 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 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 informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More 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 informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
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 informationAssessment of VMS-LES and hybrid RANS VMS-LES models
Assessment of VMS-LES and hybrid RANS VMS-LES models A. Belme 1, S. Wornom 2, A. Dervieux 1, B. Koobus 3, M.V. Salvetti 4 1 INRIA, 2004 Route des Lucioles, Sophia-Antipolis, France E-mail: Anca.Belme@inria.fr,
More informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
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 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 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 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 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 informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
More informationSAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD
1 SAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD Bruno Koobus, 1,2 Laurent Hascoët, 1 Frédéric Alauzet, 3 Adrien Loseille, 3 Youssef Mesri, 1 Alain Dervieux 1 1 INRIA
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
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 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 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 informationCurved grid generation and DG computation for the DLR-F11 high lift configuration
ECCOMAS Congress 2016, Crete Island, Greece, 5-10 June 2016 Curved grid generation and DG computation for the DLR-F11 high lift configuration Ralf Hartmann 1, Harlan McMorris 2, Tobias Leicht 1 1 DLR (German
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 informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
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 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 informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More 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 informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
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 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 informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationA conservative Embedded Boundary method for an inviscid compressible flow coupled with a fragmenting structure
A conservative Embedded Boundary method for an inviscid compressible flow coupled with a fragmenting structure Maria Adela Puscas, Laurent Monasse, Alexandre Ern, Christian Tenaud, Christian Mariotti To
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 informationVariational multiscale large-eddy simulations of the flow past a circular cylinder : Reynolds number effects
Variational multiscale large-eddy simulations of the flow past a circular cylinder : Reynolds number effects Stephen Wornom a, Hilde Ouvrard b, Maria Vittoria Salvetti c, Bruno Koobus b, Alain Dervieux
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 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 informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
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 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 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 informationAdmissibility and asymptotic-preserving scheme
Admissibility and asymptotic-preserving scheme F. Blachère 1, R. Turpault 2 1 Laboratoire de Mathématiques Jean Leray (LMJL), Université de Nantes, 2 Institut de Mathématiques de Bordeaux (IMB), Bordeaux-INP
More informationarxiv: v1 [math.na] 23 Jun 2017
On the use of kinetic energy preserving DG-schemes for large eddy simulation David Flad 1 arxiv:1706.07601v1 [math.na] 23 Jun 2017 Institute for Aerodynamics and Gas Dynamics, University of Stuttgart,
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 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 informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationA High-Order Implicit-Explicit Fluid-Structure Interaction Method for Flapping Flight
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2690 A High-Order Implicit-Explicit Fluid-Structure Interaction Method
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 informationFlow simulation on moving boundary-fitted grids and application to fluid-structure interaction problems
Flow simulation on moving boundary-fitted grids and application to fluid-structure interaction problems Martin Engel and Michael Griebel Institute of Numerical Simulation, University of Bonn, Wegelerstr.
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 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 informationFrom Navier-Stokes to Saint-Venant
From Navier-Stokes to Saint-Venant Course 1 E. Godlewski 1 & J. Sainte-Marie 1 1 ANGE team LJLL - January 2017 Euler Eulerian vs. Lagrangian description Lagrange x(t M(t = y(t, z(t flow of the trajectories
More informationLarge Eddy Simulation of a turbulent premixed flame stabilized by a backward facing step.
Large Eddy Simulation of a turbulent premixed flame stabilized by a backward facing step. Angelo Murrone and Dominique Scherrer ONERA, 29 Avenue de la Division Leclerc, BP. 72, 92322 Châtillon, France
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 informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationρ t + (ρu j ) = 0 (2.1) x j +U j = 0 (2.3) ρ +ρ U j ρ
Chapter 2 Mathematical Models The following sections present the equations which are used in the numerical simulations documented in this thesis. For clarity, equations have been presented in Cartesian
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 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 informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
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 informationPositivity-preserving Lagrangian scheme for multi-material compressible flow. Abstract
Positivity-preserving Lagrangian scheme for multi-material compressible flow Juan Cheng and Chi-Wang Shu 2 Abstract Robustness of numerical methods has attracted an increasing interest in the community
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 informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
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 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 informationOn the behavior of upwind schemes in the low Mach number limit. IV : P0 Approximation on triangular and tetrahedral cells
On the behavior of upwind schemes in the low Mach number limit. IV : P0 Approximation on triangular and tetrahedral cells Hervé Guillard To cite this version: Hervé Guillard. On the behavior of upwind
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More 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 informationFailsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model
Failsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model Matthias Möller 1 joint work with: Dmitri Kuzmin 2 1 Institute of Applied Mathematics
More informationY. Abe, N. Iizuka, T. Nonomura, K. Fujii Corresponding author:
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2012-2801 ICCFD7-2801 Symmetric-conservative metric evaluations for higher-order finite
More informationPart I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
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 informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
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 informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
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 informationEfficient and Flexible Solution Strategies for Large- Scale, Strongly Coupled Multi-Physics Analysis and Optimization Problems
University of Colorado, Boulder CU Scholar Aerospace Engineering Sciences Graduate Theses & Dissertations Aerospace Engineering Sciences Spring 1-1-2016 Efficient and Flexible Solution Strategies for Large-
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 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 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 informationModelling and numerical simulation of bi-temperature Euler Equations in toroidal geometry
Modelling and numerical simulation of bi-temperature Euler Equations in toroidal geometry E. Estibals H. Guillard A. Sangam elise.estibals@inria.fr Inria Sophia Antipolis Méditerranée June 9, 2015 1 /
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 informationHigh-Order Hyperbolic Navier-Stokes Reconstructed Discontinuous Galerkin Method
AIAA SciTech Forum 7- January 29 San Diego California AIAA Scitech 29 Forum.254/6.29-5 High-Order Hyperbolic Navier-Stokes Reconstructed Discontinuous Galerkin Method Lingquan Li Jialin Lou and Hong Luo
More informationEfficient BDF time discretization of the Navier Stokes equations with VMS LES modeling in a High Performance Computing framework
Introduction VS LES/BDF Implementation Results Conclusions Swiss Numerical Analysis Day 2015, Geneva, Switzerland, 17 April 2015. Efficient BDF time discretization of the Navier Stokes equations with VS
More informationComparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows
Comparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows Dr Argiris I. Delis Dr Ioannis K. Nikolos (TUC) Maria
More informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
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 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 informationConvergence of QuickFlow
Convergence of QuickFlow A Steady State Solver for the Shallow Water Equations Femke van Wageningen-Kessels Report literature study December 2006 PSfrag replacements Supervisors Dr.ir. C. Vuik TUDelft
More informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
More information