Engineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical

Size: px
Start display at page:

Download "Engineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical"

Transcription

1 Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013

2 Outline Outline Geurts Book 1 Geurts Book

3 Origin This lecture is strongly based on the book: Outline Geurts Book Bernard J. Geurts : Elements of Direct and Simulation books/072/index.html

4 Part I Second Lecture

5 The idea Reduced flow description Low-pass filtering (in wave number space), with characteristic length of non-linear terms Closure problems Important to consider the ratio /η. As /η O(1), LES converge to DNS. /η is the main parameter of the LES. It characterise the resolution/quality of various results of the simulation.

6 The governing Incompressible, isochoric fluid ρ = const. Constant viscosity ν = const. Continuity equation: i u i = 0 (1) t u i + u j j u i = p i ρ + ν j j u i (2) Let us assume that ρ = 1 so we can skip it from the.

7 Transformation properties of the Galilean invariance etc.

8 Notation of the filtering Let us define an operator L which reduces the spatial scales of the flow-field: L : ϕ ϕ 1D example (ϕ(x, t) or e.g. u(x, t)): ϕ = L(ϕ)(x, t) = x+ /2 x /2 ϕ(ξ, t) dξ (3) 2) = l/16 Two filtering is shown: 1) = l/4,

9 Preliminary comment on filter width We are interested in a quantitative variable q of a specific flow (for example flow around cylinder, channel flow). The filtered flow-field can be used also for the evaluation. Depending on the interesting quantity (drag coefficient, turbulence level), the accuracy of the estimation q err will clearly depend on the applied filter length i.e. how many details of the flow we keep. For some quantity the very smoothed signal ( is big) gives already a good prediction (q err is small), for other quantities fines scales are also responsible and only filtering at smaller scales can be applied.

10 Notation of the filtering In 3D a generic homogeneous filter can be defined: ϕ (x, t) = L(ϕ)(x, t) = + G (x x)ϕ(x, t) dx = G ϕ (4) This a convolution filter (multiplied by something shifted over)! No temporal filtering is carried out! Causality: future does not influence the past Temporal averaging would lead us to RANS Normalization of the kernel: + G (x ) dx = 1 (5)

11 Commutation of the filter This type of filtering is a linear integral operator. The above defined homogeneous convolution filter commutes both with temporal and spatial derivation: t f = t f ; i f = i f (6)

12 Commuting of the mom. eq. Because of the commutation property: (ρ = 1 and everything is rearranged to the LHS) t u i + j ( u i u j ) + i p ν j j u i = 0 (7) If we want to see the same terms on LHS as in the orig. NS. we have to add the required term j ( u i u j ) and subtract the one available j ( u i u j ) on both side of the equation: t u i + j ( u i u j ) + i p ν j j u i = j τ ij (8) where: τ ij def = u i u j u i u j (9)

13 SGS tensor τ ij is called the SubGrid Scale stress tensor. It is responsible for the effect of the unresolved (removed by the filter) scales. Its actual effect is j τ ij which needs to be led. It is called historically sub grid scale, since the filtering was originally (but not very logically) related to the numerical grid size (h).

14 SGS tensor as a commutator Let us define the product operator: Π(u i, u j ) def = u i u j (10) The SGS stress can be written using this as a commutator: τ ij = u i u j u i u j = = L(Π(u i, u j )) Π(L(u i ), L(u i )) = = [L, Π](u i, u j ) (11) where [a, b] def = ab ba is the commutator symbol for spaces where addition and multiplication are defined. It measures the strength of not commuting.

15 Identities of the commutator operator Wikipedia section 1 [A, A] = 0 2 [A, B] = [B, A], anti-commutativity 3 [A, [B, C]] + [B, [C, A]] + [C, [A, B]] = 0, Jacobi identity 4 [A, BC] = [A, B]C + B[A, C], Leibnitz rule Example for Leibnitz rule: A B C def = d x ( ) = ( ) def = f ( ) def = g( ) (12) I.e.: (fg) = f g + fg

16 Non-uniform filtering In convolution filtering the filter width ( ) is constant. Practically the mesh size is changing (e.g.: see channel in wall normal direction)! 1D example (ϕ(x, t) or e.g. u(x, t)): ϕ = L(ϕ)(x, t) = x+ +(x) x (x) g(x, ξ) ϕ(ξ, t) dξ (13) (x) where (x) = + + is the spatially varying filter width.

17 Commutation error at non-uniform filtering i ϕ i ϕ because of the varying filter width. This can be written more formally: L( i ϕ) i L(ϕ) = [L, i ](ϕ) 0 (14)

18 Non-uniformly filtered set of Continuity The continuity equation: j u j = [L, j ](u j ) (15) i.e. the filtered flow-field is not any more divergence free (solenoidal). Variation of the filter width is related to local mass production or destruction.

19 Non-uniformly filtered set of Momentum The momentum equation: t u i + j ( u i u j ) + i p ν j j u i = j ([L, Π](u i, u j )) ([L, j ](p) ν[l, jj ](u i ) + [L, j ](Π(u i, u j )) (16) The RHS is much more complicated now. The new terms even destroy the conservation property of the original. Flow structures travelling in the direction of filter width change experience a scale reduction or occurrence simply by the change in the filter.

20 Non-uniformly filtered set of Momentum E.g. a vortex arriving to a region of higher filter width will be smoothed corresponding to local kinetic energy reduction, which is against conservation. The effect can be quantified by material derivative of the filter width: u j j

21 Summary We arrive at the LES-template : NS( U ) = R(U, U ) (17) where R is the residual because of the filtering. R(U, U) = 0, since in this case L = Id identity operator, and DNS should be recovered. In the ling exercise R(U, U ) M( U ), i.e. the residual has to be approximated using the filtered flow-field.

22 How to develop the? How to separate between large and small scales? Spatial filtering, smoothing using a kernel function ϕ (x j, t) def = G (r i ; x j ) ϕ(x j r i, t)dr i (18) V

23 kernel G is the filtering kernel with a typical size of. G has a compact support (its definition set where the value is non-zero is closed) in its first variable To be the filtered value of a constant itself it has to be true: G (r i ; x j )dr i = 1 (19) V If G (r i ; x j ) is homogeneous in its second variable and isotropic in its first variable than G ( r i ) is a function of only one variable

24 kernel Examples

25 Physical space Fluctuation: ϕ def = ϕ ϕ (20) ϕ 0, a difference compared to Reynolds averaging

26 Spectral space Recall: the cutting wavenumber (κ c ), below which ling is needed

27 If using the previously defined (homogeneous, isotropic) filter Averaging and the derivatives commute (exchangeable) i u i = 0 (21) t u i + u j j u i = 1 ρ p + ν j j u i j τ ij (22) 3D (because turbulence is 3D) unsteady (because the large eddies are unsteady)

28 Sub Grid Scale stress τ ij is called Sub-Grid Scale stress SGS from the times when filtering was directly associated to the grid τ ij def = u i u j u i u j (23) It represents the effect of the filtered scales It is in a form of a stress tensor Should be dissipative to represent the dissipation on the filtered small scale

29 Same as in RANS τ ij 1 3 τ kkδ ij = 2ν t s ij (24) Relatively a better approach (compared to the application for RANS) since the small scales are more universal Dissipative if ν t > 0.

30 Smagorinsky ν t = (C s ) 2 S (25) S def = 2s ij s ij (26) C s Smagorinsky constant to be determined using spectral theory of turbulence using validations on real flow computations to be prescribed determine the computational cost (if too small) determine the accuracy (if too big) resolving 80% of the energy is a good compromise

31 Scale Similarity Let us assume that the cuted small scales are similar to the kept large scales! A logical : τ ij def = u i u j u i u j (27) The unfiltered quantities are led by the filtered ones. Properties It is not dissipative enough It gives feasible shear stresses (from experience) Logical to combine with Smag.!

32 Dynamic approach The idea is the same as in the scale similarity The theory is more complicated Any can be made dynamic Dynamic Smagorinsky is widely used (combining the two advantages)

33 Periodicity Periodicity is used to infinite long domain The length of periodicity is given by the length scales of turbulence

34 Inlet Much more difficult than in RANS Turbulent structures should be represented Vortices should be synthesized Separate precursor simulation to provide real turbulence

35 Wall y + 1 (28) x + 50 (29) z (30)

Turbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.

Turbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics. Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena

More information

Modelling of turbulent flows: RANS and LES

Modelling of turbulent flows: RANS and LES Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012

More information

Reliability of LES in complex applications

Reliability of LES in complex applications Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow

More information

Turbulence Modeling I!

Turbulence Modeling I! Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!

More information

Engineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical

Engineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical Outline Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2011 Outline Outline Part I First Lecture Connection between time and ensemble average Ergodicity1 Ergodicity

More information

Turbulence: Basic Physics and Engineering Modeling

Turbulence: Basic Physics and Engineering Modeling DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois

More information

An Introduction to Theories of Turbulence. James Glimm Stony Brook University

An Introduction to Theories of Turbulence. James Glimm Stony Brook University An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing

More information

An evaluation of a conservative fourth order DNS code in turbulent channel flow

An evaluation of a conservative fourth order DNS code in turbulent channel flow Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical

More information

Computational Fluid Dynamics 2

Computational Fluid Dynamics 2 Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2

More information

Regularization modeling of turbulent mixing; sweeping the scales

Regularization 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 information

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline

More information

Database analysis of errors in large-eddy simulation

Database analysis of errors in large-eddy simulation PHYSICS OF FLUIDS VOLUME 15, NUMBER 9 SEPTEMBER 2003 Johan Meyers a) Department of Mechanical Engineering, Katholieke Universiteit Leuven, Celestijnenlaan 300, B3001 Leuven, Belgium Bernard J. Geurts b)

More information

A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries

A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation

More information

Turbulent Boundary Layers & Turbulence Models. Lecture 09

Turbulent Boundary Layers & Turbulence Models. Lecture 09 Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects

More information

Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh

Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh 40th Fluid Dynamics Conference and Exhibit 8 June - July 00, Chicago, Illinois AIAA 00-506 Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows Xiaochuan Chai and Krishnan Mahesh University

More information

A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS

A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS June 30 - July 3, 2015 Melbourne, Australia 9 7B-4 A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS C.-Y. Chang, S. Jakirlić, B. Krumbein and C. Tropea Institute of Fluid Mechanics and Aerodynamics /

More information

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and

More information

2. Conservation Equations for Turbulent Flows

2. Conservation Equations for Turbulent Flows 2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging

More information

Lecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.

Lecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of

More information

Eulerian models. 2.1 Basic equations

Eulerian models. 2.1 Basic equations 2 Eulerian models In this chapter we give a short overview of the Eulerian techniques for modelling turbulent flows, transport and chemical reactions. We first present the basic Eulerian equations describing

More information

COMMUTATION ERRORS IN PITM SIMULATION

COMMUTATION ERRORS IN PITM SIMULATION COMMUTATION ERRORS IN PITM SIMULATION B. Chaouat ONERA, 93 Châtillon, France Bruno.Chaouat@onera.fr Introduction Large eddy simulation is a promising route. This approach has been largely developed in

More information

HYBRID LES-RANS: Inlet Boundary Conditions for Flows With Recirculation

HYBRID LES-RANS: Inlet Boundary Conditions for Flows With Recirculation 1 HYBRID LES-RANS: Inlet Boundary Conditions for Flows With Recirculation Lars Davidson Div. of Fluid Dynamics Dept. of Applied Mechanics Chalmers University of Technology, Göteborg, Sweden E-mail lada@chalmers.se

More information

CHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.

CHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1. 1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy

More information

Introduction to Turbulence and Turbulence Modeling

Introduction 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 information

Multiscale Computation of Isotropic Homogeneous Turbulent Flow

Multiscale Computation of Isotropic Homogeneous Turbulent Flow Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible

More information

Before we consider two canonical turbulent flows we need a general description of turbulence.

Before we consider two canonical turbulent flows we need a general description of turbulence. Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent

More information

Nonequilibrium Dynamics in Astrophysics and Material Science YITP, Kyoto

Nonequilibrium Dynamics in Astrophysics and Material Science YITP, Kyoto Nonequilibrium Dynamics in Astrophysics and Material Science 2011-11-02 @ YITP, Kyoto Multi-scale coherent structures and their role in the Richardson cascade of turbulence Susumu Goto (Okayama Univ.)

More information

Three-dimensional wall filtering formulation for large-eddy simulation

Three-dimensional wall filtering formulation for large-eddy simulation Center for Turbulence Research Annual Research Briefs 6 55 Three-dimensional wall filtering formulation for large-eddy simulation By M. Shoeybi AND J. A. Templeton 1. Motivation and objectives Large-eddy

More information

Chapter 13. Eddy Diffusivity

Chapter 13. Eddy Diffusivity Chapter 13 Eddy Diffusivity Glenn introduced the mean field approximation of turbulence in two-layer quasigesotrophic turbulence. In that approximation one must solve the zonally averaged equations for

More information

Some remarks on grad-div stabilization of incompressible flow simulations

Some 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

Estimation of Turbulent Dissipation Rate Using 2D Data in Channel Flows

Estimation of Turbulent Dissipation Rate Using 2D Data in Channel Flows Proceedings of the 3 rd World Congress on Mechanical, Chemical, and Material Engineering (MCM'17) Rome, Italy June 8 10, 2017 Paper No. HTFF 140 ISSN: 2369-8136 DOI: 10.11159/htff17.140 Estimation of Turbulent

More information

Anisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray

Anisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale

More information

Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions

Simulating 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 information

SUBGRID MODELS FOR LARGE EDDY SIMULATION: SCALAR FLUX, SCALAR DISSIPATION AND ENERGY DISSIPATION

SUBGRID MODELS FOR LARGE EDDY SIMULATION: SCALAR FLUX, SCALAR DISSIPATION AND ENERGY DISSIPATION SUBGRID MODELS FOR LARGE EDDY SIMULATION: SCALAR FLUX, SCALAR DISSIPATION AND ENERGY DISSIPATION By Sergei G. Chumakov A dissertation submitted in partial fulfillment of the requirements for the degree

More information

LES 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. 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 information

On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows

On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows Center for Turbulence Research Annual Research Briefs 1998 267 On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows By Jeffrey S. Baggett 1. Motivation and objectives

More information

Optimizing calculation costs of tubulent flows with RANS/LES methods

Optimizing calculation costs of tubulent flows with RANS/LES methods Optimizing calculation costs of tubulent flows with RANS/LES methods Investigation in separated flows C. Friess, R. Manceau Dpt. Fluid Flow, Heat Transfer, Combustion Institute PPrime, CNRS University

More information

Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows

Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations

More information

LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS

LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,

More information

Turbulence modelling. Sørensen, Niels N. Publication date: Link back to DTU Orbit

Turbulence modelling. Sørensen, Niels N. Publication date: Link back to DTU Orbit Downloaded from orbit.dtu.dk on: Dec 19, 2017 Turbulence modelling Sørensen, Niels N. Publication date: 2010 Link back to DTU Orbit Citation (APA): Sørensen, N. N. (2010). Turbulence modelling. Paper presented

More information

Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions

Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering

More information

Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step

Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step Copyright c 2004 Tech Science Press CMC, vol.1, no.3, pp.275-288, 2004 Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step B. Wang 1, H.Q. Zhang 1, C.K. Chan 2 and X.L.

More information

L.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133, ORSAY CEDEX, FRANCE fax number:

L.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133, ORSAY CEDEX, FRANCE  fax number: Large Eddy Simulations of a spatially developing incompressible 3D mixing layer using the v-ω formulation. C. TENAUD, S. PELLERIN, A. DULIEU and L. TA PHUOC L.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133,

More information

Green s Functions and Distributions

Green s Functions and Distributions CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where

More information

Turbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i

Turbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)

More information

Boundary 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 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 information

Hybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model

Hybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de

More information

Mass Transfer in Turbulent Flow

Mass Transfer in Turbulent Flow Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.

More information

RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD

RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,

More information

meters, we can re-arrange this expression to give

meters, we can re-arrange this expression to give Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.

More information

Flow Structure Investigations in a "Tornado" Combustor

Flow Structure Investigations in a Tornado Combustor Flow Structure Investigations in a "Tornado" Combustor Igor Matveev Applied Plasma Technologies, Falls Church, Virginia, 46 Serhiy Serbin National University of Shipbuilding, Mikolayiv, Ukraine, 545 Thomas

More information

model and its application to channel ow By K. B. Shah AND J. H. Ferziger

model and its application to channel ow By K. B. Shah AND J. H. Ferziger Center for Turbulence Research Annual Research Briefs 1995 73 A new non-eddy viscosity subgrid-scale model and its application to channel ow 1. Motivation and objectives By K. B. Shah AND J. H. Ferziger

More information

Predicting natural transition using large eddy simulation

Predicting natural transition using large eddy simulation Center for Turbulence Research Annual Research Briefs 2011 97 Predicting natural transition using large eddy simulation By T. Sayadi AND P. Moin 1. Motivation and objectives Transition has a big impact

More information

Basic Fluid Mechanics

Basic Fluid Mechanics Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible

More information

COMPARISON OF DIFFERENT SUBGRID TURBULENCE MODELS AND BOUNDARY CONDITIONS FOR LARGE-EDDY-SIMULATIONS OF ROOM AIR FLOWS.

COMPARISON OF DIFFERENT SUBGRID TURBULENCE MODELS AND BOUNDARY CONDITIONS FOR LARGE-EDDY-SIMULATIONS OF ROOM AIR FLOWS. 7 TH INTRNATINAL CNFRNC N AIR DISTRIBTIN IN RMS, RMVNT 2 pp. 31-36 CMPARISN F DIFFRNT SBGRID TRBLNC MDLS AND BNDARY CNDITINS FR LARG-DDY-SIMLATINS F RM AIR FLWS. D. Müller 1, L. Davidson 2 1 Lehrstuhl

More information

The lattice Boltzmann equation (LBE) has become an alternative method for solving various fluid dynamic

The lattice Boltzmann equation (LBE) has become an alternative method for solving various fluid dynamic 36th AIAA Fluid Dynamics Conference and Exhibit 5-8 June 2006, San Francisco, California AIAA 2006-3904 Direct and Large-Eddy Simulation of Decaying and Forced Isotropic Turbulence Using Lattice Boltzmann

More information

A Finite-Element based Navier-Stokes Solver for LES

A Finite-Element based Navier-Stokes Solver for LES A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,

More information

Explicit algebraic Reynolds stress models for boundary layer flows

Explicit algebraic Reynolds stress models for boundary layer flows 1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both

More information

Turbulence - Theory and Modelling GROUP-STUDIES:

Turbulence - Theory and Modelling GROUP-STUDIES: Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence

More information

Automatic Eddy Viscosity Assignment for 2-D Hydrodynamic Model of Szczecin Bay

Automatic Eddy Viscosity Assignment for 2-D Hydrodynamic Model of Szczecin Bay PUBLS. INST. GEOPHYS. POL. ACAD. SC., E-6 (390), 006 Automatic Eddy Viscosity Assignment for -D Hydrodynamic Model of Szczecin Bay Ryszard EWERTOWSKI 1, 1 Technical University of Szczecin, Building and

More information

Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions

Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions Journal of Physics: Conference Series PAPER OPEN ACCESS Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions To cite this article: A Rosiak and

More information

2 GOVERNING EQUATIONS

2 GOVERNING EQUATIONS 2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates

More information

Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence

Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal and Vijaya Rani PI: Sarma L. Rani Department of Mechanical and Aerospace

More information

A Low Reynolds Number Variant of Partially-Averaged Navier-Stokes Model for Turbulence

A Low Reynolds Number Variant of Partially-Averaged Navier-Stokes Model for Turbulence Int. J. Heat Fluid Flow, Vol., pp. 65-669 (), doi:.6/j.ijheatfluidflow... A Low Reynolds Number Variant of Partially-Averaged Navier-Stokes Model for Turbulence J.M. Ma,, S.-H. Peng,, L. Davidson, and

More information

Module 3: Velocity Measurement Lecture 15: Processing velocity vectors. The Lecture Contains: Data Analysis from Velocity Vectors

Module 3: Velocity Measurement Lecture 15: Processing velocity vectors. The Lecture Contains: Data Analysis from Velocity Vectors The Lecture Contains: Data Analysis from Velocity Vectors Velocity Differentials Vorticity and Circulation RMS Velocity Drag Coefficient Streamlines Turbulent Kinetic Energy Budget file:///g /optical_measurement/lecture15/15_1.htm[5/7/2012

More information

MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran

MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER

More information

Accretion Disks I. High Energy Astrophysics: Accretion Disks I 1/60

Accretion Disks I. High Energy Astrophysics: Accretion Disks I 1/60 Accretion Disks I References: Accretion Power in Astrophysics, J. Frank, A. King and D. Raine. High Energy Astrophysics, Vol. 2, M.S. Longair, Cambridge University Press Active Galactic Nuclei, J.H. Krolik,

More information

Modeling of turbulence in stirred vessels using large eddy simulation

Modeling of turbulence in stirred vessels using large eddy simulation Modeling of turbulence in stirred vessels using large eddy simulation André Bakker (presenter), Kumar Dhanasekharan, Ahmad Haidari, and Sung-Eun Kim Fluent Inc. Presented at CHISA 2002 August 25-29, Prague,

More information

Feedback Control of Turbulent Wall Flows

Feedback Control of Turbulent Wall Flows Feedback Control of Turbulent Wall Flows Dipartimento di Ingegneria Aerospaziale Politecnico di Milano Outline Introduction Standard approach Wiener-Hopf approach Conclusions Drag reduction A model problem:

More information

RANS-LES inlet boundary condition for aerodynamic and aero-acoustic. acoustic applications. Fabrice Mathey Davor Cokljat Fluent Inc.

RANS-LES inlet boundary condition for aerodynamic and aero-acoustic. acoustic applications. Fabrice Mathey Davor Cokljat Fluent Inc. RANS-LES inlet boundary condition for aerodynamic and aero-acoustic acoustic applications Fabrice Mathey Davor Cokljat Fluent Inc. Presented by Fredrik Carlsson Fluent Sweden ZONAL MULTI-DOMAIN RANS/LES

More information

Multi-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah

Multi-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary

More information

Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles

Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &

More information

On the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows

On the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows PHYSICS OF FLUIDS VOLUME 11, NUMBER 5 MAY 1999 On the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows L. Shao a) Laboratoire de Mécanique des Fluides

More information

1 Mean Flow Equations

1 Mean Flow Equations ME 543 Averaged Equations 23 January 2018 These notes are meant as a supplement to the text, and as an aid for the class notes, but not as an independent document. 1 Mean Flow Equations 1.1 Mean flow equations

More information

Turbulence models and excitation of solar oscillation modes

Turbulence models and excitation of solar oscillation modes Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives

More information

Publication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada

Publication 97/2. An Introduction to Turbulence Models. Lars Davidson,   lada ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature

More information

Modeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!

Modeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics! http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations

More information

A G-equation formulation for large-eddy simulation of premixed turbulent combustion

A G-equation formulation for large-eddy simulation of premixed turbulent combustion Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent

More information

cfl Copyright by Fotini V. Katopodes 2000 All rights reserved.

cfl Copyright by Fotini V. Katopodes 2000 All rights reserved. A theory for the subfilter-scale model in large-eddy simulation Fotini V. Katopodes, Robert L. Street, Joel H. Ferziger March, 2000 Technical Report 2000-K1 Environmental Fluid Mechanics Laboratory Stanford,

More information

LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL

LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals 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 information

Interaction(s) fluide-structure & modélisation de la turbulence

Interaction(s) fluide-structure & modélisation de la turbulence Interaction(s) fluide-structure & modélisation de la turbulence Pierre Sagaut Institut Jean Le Rond d Alembert Université Pierre et Marie Curie- Paris 6, France http://www.ida.upmc.fr/~sagaut GDR Turbulence

More information

36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson

36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson 36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through

More information

6 VORTICITY DYNAMICS 41

6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS As mentioned in the introduction, turbulence is rotational and characterized by large uctuations in vorticity. In this section we would like to identify some

More information

A NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS

A NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS June 3 - July 3, 25 Melbourne, Australia 9 B-3 A NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS Ryan N. King Department of Mechanical Engineering University of Colorado Boulder, CO, 839, USA ryan.n.king@colorado.edu

More information

Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena

Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena 29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied

More information

Application of Compact Schemes to Large Eddy Simulation of Turbulent Jets

Application of Compact Schemes to Large Eddy Simulation of Turbulent Jets Journal of Scientific Computing, Vol. 21, No. 3, December 2004 ( 2004) Application of Compact Schemes to Large Eddy Simulation of Turbulent Jets Ali Uzun, 1 Gregory A. Blaisdell, 2 and Anastasios S. Lyrintzis

More information

Introduction to Turbulence, its Numerical Modeling and Downscaling

Introduction to Turbulence, its Numerical Modeling and Downscaling 560.700: Applications of Science-Based Coupling of Models October, 2008 Introduction to Turbulence, its Numerical Modeling and Downscaling Charles Meneveau Department of Mechanical Engineering Center for

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

Lecture 3: The Navier-Stokes Equations: Topological aspects

Lecture 3: The Navier-Stokes Equations: Topological aspects Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without

More information

Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue

Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue 11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates

More information

Direct Numerical Simulations of converging-diverging channel flow

Direct Numerical Simulations of converging-diverging channel flow Intro Numerical code Results Conclusion Direct Numerical Simulations of converging-diverging channel flow J.-P. Laval (1), M. Marquillie (1) Jean-Philippe.Laval@univ-lille1.fr (1) Laboratoire de Me canique

More information

AER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly

AER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly 1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence

More information

Validation of an Entropy-Viscosity Model for Large Eddy Simulation

Validation of an Entropy-Viscosity Model for Large Eddy Simulation Validation of an Entropy-Viscosity Model for Large Eddy Simulation J.-L. Guermond, A. Larios and T. Thompson 1 Introduction A primary mainstay of difficulty when working with problems of very high Reynolds

More information

LES ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT

LES ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT 2th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT Wang T.*, Gao S.F., Liu Y.W., Lu Z.H. and Hu H.P. *Author

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

AA214B: 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 information

A Simple Turbulence Closure Model

A Simple Turbulence Closure Model A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: Mean velocity: Turbulent velocity: Gradient operator: Advection operator: V =

More information

Inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows. UPMC Univ Paris 06, CNRS, UMR 7190, Paris, F-75005, France a)

Inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows. UPMC Univ Paris 06, CNRS, UMR 7190, Paris, F-75005, France a) Inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows Y. Ling, 1 S. Balachandar, 2 and M. Parmar 2 1) Institut Jean Le Rond d Alembert, Sorbonne Universités, UPMC Univ Paris

More information

Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows

Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of

More information

Modeling and simulation of bedload transport with viscous effects

Modeling 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 information