Introduc)on to modeling turbulence in oceanic flows

Size: px
Start display at page:

Download "Introduc)on to modeling turbulence in oceanic flows"

Transcription

1 Introducon to modeling turbulence in oceanic flows Mehmet Ilıcak Uni Research, Bjerknes Centre for Climate Research, Bergen Norway. School on Ocean Climate Modelling, 05, Ankara, Turkey 09/8/05

2 The biggest challenge in geophysical flows Large scale motion, L ~ O0-000 km. High Reynolds number, Re UL/ν ~ O 0. Sampling/observation is a challenge: Lab. experiments are limited. Difficult to get χx,y,z,t in the real ocean/atmosphere. Rotation of Earth. Stratification: effect of density differences Vertical displacement are restricted. Mixing between different density layers. Presence of internal waves. 09/8/05

3 Geophysical flows Gulf Stream Great Red Spot gulfstream0image.gif 09/8/05

4 Turbulent coherent structures Lab experiment: Re400 Oil discharge from a ship: Re~0 7 from p. 00 of Van Dyke 98 09/8/05 4

5 Navier- Stokes equaons!u!t u!u!x v!u!y w!u!z "!p!x! #!x "!u $!x!v!t u!v!x v!v!y w!v!z "!p!y! #!x "!v $!x!w!t u!w!x v!w!y w!w!z "!p!z " g! '! #!x "!w $!x!u i!t u j!u i "!p! ' " g i!!x j!x i! # '!y "!u! # $!y '!z "!u $!z '! # '!y "!v! # $!y '!z "!v $!z '! # '!y "!w! # $!y '!z "!w $!z ' # "!u i!x j $!x j ' i, j,, 09/8/05 5

6 Modeling turbulent spectrum!u i!t u j!u i "!p! ' " g i!!x j!x i # "!u i!x j $!x j ' Re UL! Re [m / s]!00!0 [m] 0 " [m / s] 0 logek Energy containing range P ~k -5/ Inertial range fre Dissipation range ε logk 09/8/05

7 Possible modeling approaches Direct Numerical Simulations Resolved all scales of motion i.e. molecular viscosity, but how? DNS!u i!t u j RANS!u i "!p! ' " g i!!x j!x i LES # "!u i!x j $!x j ' Ilıcak et al /8/05 7

8 Possible modeling approaches DNS Resolved all scales N ~ Re 9/4 RANS Re! U!x RANS Resolved! mean state!u LES i!t u!u i j "!p! ' " g i! # "!u i!x Resolved large eddies j!x i!x j $!x j DNS LES ' Ilıcak et al /8/05 8

9 Kelvin-Helmholtz instability from p. 89 of Van Dyke 98 Upper fluid column is moving to the right faster then the lower one. 09/8/05 9

10 Direct numerical simulation DNS z D Non-hydrostatic model grid points in x,y and z directions Run in 8 processors x 09/8/05 0

11 Reynolds Averaging Navier- Stokes RANS u i U i uʹ i!u i!t!u U i j! u" j u i "!x j!x j #!p! # g i! $ "!U i!x i!x j!x j ' uʹ uʹ i t j U k uʹ uʹ i x k j uʹ uʹ uʹ k x i k j P ij ij ε ij x k uʹ ʹ iu j ν xk # unknowns > # equaons System is NOT closed The famous Turbulence Closure Problem 09/8/05

12 Turbulence closure models!u i!t!u U i j! u" j u i "!x j!x j #!p! # g i! $ "!U i!x i!x j!x j ' Eddy viscosity approach U uʹ wʹ K M z We need to find a way to compute eddy viscosity/diffusivity. Ex: KPP, k- epsilon, k- omega, MY, standard Smagorinsky. 09/8/05

13 Large Eddy Simulaons LES Large Eddy Simulation LES: Mixing by the large, energy-containing, eddies is handled through computation, while the effect of small, dissipative, turbulent eddies is modeled analytically.! < u > i!t u i < u > i u i!! < u > < u > i j!! ij "!p! " g i!!x j!x j!x i # "! < u > i!x j $!x j '! ij K M! < u > i!x j! < u > j!x i 09/8/05

14 Possible modeling approaches DNS N ~ Re 9/ 4 Resolved all scales RANS Resolved mean state LES Resolved large eddies DNS RANS LES Ilıcak et al /8/05 4

15 Possible modeling closures SGS turbulence closures; RANS Constant eddy viscosity/diffusivity Prandtl s 95 eq. model KPP Large et al. 994 eq. turbulence closures k-epsilon, k-omega, k- kl, MY LES Constant eddy viscosity/diffusivity Standard Smagorinsky Dynamical Smagorinsky #U u! w! "K M #z K M K M C! # $ kl "U i "x j "U j "x i ' 09/8/05 5

16 Ocean model equaons!u!t u!u!x v!u!y w!u!z "!p!x! #!x A M $!u!x! # '!y A M $!u!y! # '!z K M $!u!z '!v!t u!v!x v!v!y w!v!z "!p!y! #!x A M $!v!x! # '!y A M $!v!y! # '!z K M $!v!z '!T!t u!t!x v!t!y w!t!z! " $!x A H # Smagorinsky Biharmonic viscosity Leith Numerical i.e. ROMS!T!x '! " $!y A H #!T!y '! " $!z K H # Constant eddy viscosity Ad- hoc i.e. Price- Turner, KPP One/two eq. turbulence schemes.!t!z ' 09/8/05

17 One-equation turbulence closure K M kl! k / l where l is prescribed. 09/8/05 7

18 Two-equation turbulence closures K K M H c c kls kls M H 09/8/05 8

19 Dynamical Smagorinsky c < M L > ij ij < M M > ds kl M ij!!!! s!u! s!u ij "!! s u! s uij! "! s!u! s!u "!s ij u! s u! M ij! M ij /! kl c ds! < M Lij ij >!! < M kl M kl > 09/8/05 9

20 Eddy viscosity diffusivity #U u! w! "K M #z #T T 'w' "K H #z K K M H c c kls kls M H Pr t K K M H S S M H K M Eddy viscosity K H Eddy diffusivity kturbulent kinec energy lturbulent length scale Pr t Turbulent Prandtl Number S M and S H Very complex stability funcons 09/8/05 0

21 Pr t Ri g Pr K Pr t M 0 Pr π 0 0.; Ri exp Pr0 Ri k ; ε Ri f K H g f 0. K Pr M t Ri Ri g f Ilicak et al /8/05

22 Red Sea overflow T0 T 09/8/05

23 Experiments KPP K-Profile Parameterization PB07 Peters and Baumert 007 CA Canuto-A, CB Canuto-B G88 Galperin et al. 988, KC Kantha-Clayson 09/8/05

24 Propagaon of overflow 09/8/05 4

25 Eddy diffusivity Northern Channel 09/8/05 5

26 Summary There are different approaches to tackle turbulence in ocean models; DNS, RANS, LES Eddy viscosity/diffusivity approach is a modelling strategy. Different complexity of parameterizaons out there. Be aware of shortcomings and advantages of these models. 09/8/05

27 Upwelling test case C P P C 09/8/05 7

28 Upwelling test case 09/8/05 8

29 Upwelling test case 09/8/05 9

30 Stability funcons / ; / m m h h m h m h M m h H G f b G G f b G f b G f b G f b b cff cff G f s G f s s S cff G f s G f s s S Galperin et al. 988 h h H M h H A A G G S A A A A B S B A A G B A A S / ; / Canuto A and B / 4 / ; 4 4 / ; 7 ; 4 / / / ; ; /8 / / ; / µ c f b b b b b b s s s s s s 09/8/05 0

31 Hydrostatic vs. Non-hydrostatic Hydrostac model Non- hydrostac model Cannot resolve shear instabilities such as KH rolls. Depend on full closures such as nd order turbulence models Algebraic models P z Can resolve shear instabilities. Depend on resolution and subgrid scale parameterization LES VLES dw P ρg κ w dt z ρ ρg ρ 09/8/05 0 0

32 Vercal coordinate in OGCMs ρ z * σ 09/8/05

33 The turbulent Prandtl number Venayagamoorthy and Stretch 00 09/8/05

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

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

LES tutorial for Ocean 3D+1 MURI

LES tutorial for Ocean 3D+1 MURI LES tutorial for Ocean 3D+1 MURI Tamay Özgökmen, University of Miami Contents 1 A Discussion of Ocean Modeling Approaches 2 2 Large Eddy Simulation 4 2.1 The Concept of Spatial filtering.........................

More information

Quarter degree resolution ocean component of the Norwegian Earth System Model

Quarter degree resolution ocean component of the Norwegian Earth System Model Quarter degree resolution ocean component of the Norwegian Earth System Model Mats Bentsen 1,2, Mehmet Ilicak 1,2, and Helge Drange 3,2 1 Uni, Uni Research Ltd 2 Bjerknes Centre for Research 3 Geophysical

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

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

The Generalized Boundary Layer Equations

The Generalized Boundary Layer Equations The Generalized Boundary Layer Equations Gareth H. McKinley (MIT-HML) November 004 We have seen that in general high Reynolds number flow past a slender body such as an airfoil can be considered as an

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

Frictional boundary layers

Frictional boundary layers Chapter 12 Frictional boundary layers Turbulent convective boundary layers were an example of boundary layers generated by buoyancy forcing at a boundary. Another example of anisotropic, inhomogeneous

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

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

2.3 The Turbulent Flat Plate Boundary Layer

2.3 The Turbulent Flat Plate Boundary Layer Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The

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 Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:

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

Understanding and modeling dense overflows. Sonya Legg Princeton University AOMIP/FAMOS school for young scientists 2012

Understanding and modeling dense overflows. Sonya Legg Princeton University AOMIP/FAMOS school for young scientists 2012 Understanding and modeling dense overflows Sonya Legg Princeton University AOMIP/FAMOS school for young scientists 2012 What is an overflow? Dense water formation on shelf or marginal sea Dense water accelerates

More information

Buoyancy Fluxes in a Stratified Fluid

Buoyancy Fluxes in a Stratified Fluid 27 Buoyancy Fluxes in a Stratified Fluid G. N. Ivey, J. Imberger and J. R. Koseff Abstract Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed

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

Simulations for Enhancing Aerodynamic Designs

Simulations for Enhancing Aerodynamic Designs Simulations for Enhancing Aerodynamic Designs 2. Governing Equations and Turbulence Models by Dr. KANNAN B T, M.E (Aero), M.B.A (Airline & Airport), PhD (Aerospace Engg), Grad.Ae.S.I, M.I.E, M.I.A.Eng,

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

Diego, La Jolla, CA, USA Published online: 14 May 2014.

Diego, La Jolla, CA, USA Published online: 14 May 2014. This article was downloaded by: [University of California, San Diego] On: 17 November 2014, At: 16:21 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954

More information

NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS

NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44

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

(1) Transition from one to another laminar flow. (a) Thermal instability: Bernard Problem

(1) Transition from one to another laminar flow. (a) Thermal instability: Bernard Problem Professor Fred Stern Fall 2014 1 Chapter 6: Viscous Flow in Ducts 6.2 Stability and Transition Stability: can a physical state withstand a disturbance and still return to its original state. In fluid mechanics,

More information

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

Statistical turbulence modelling in the oceanic mixed layer

Statistical turbulence modelling in the oceanic mixed layer Statistical turbulence modelling in the oceanic mixed layer Hans Burchard hans.burchard@io-warnemuende.de Baltic Sea Research Institute Warnemünde, Germany Seminar at Bjerknes Centre, October 16, 23, Bergen,

More information

Ensemble averaged dynamic modeling. By D. Carati 1,A.Wray 2 AND W. Cabot 3

Ensemble averaged dynamic modeling. By D. Carati 1,A.Wray 2 AND W. Cabot 3 Center for Turbulence Research Proceedings of the Summer Program 1996 237 Ensemble averaged dynamic modeling By D. Carati 1,A.Wray 2 AND W. Cabot 3 The possibility of using the information from simultaneous

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

What is Turbulence? Fabian Waleffe. Depts of Mathematics and Engineering Physics University of Wisconsin, Madison

What is Turbulence? Fabian Waleffe. Depts of Mathematics and Engineering Physics University of Wisconsin, Madison What is Turbulence? Fabian Waleffe Depts of Mathematics and Engineering Physics University of Wisconsin, Madison it s all around,... and inside us! Leonardo da Vinci (c. 1500) River flow, pipe flow, flow

More information

Rayleigh-Taylor Unstable Flames

Rayleigh-Taylor Unstable Flames Rayleigh-Taylor Unstable Flames Elizabeth P. Hicks 1,2 and Robert Rosner 2 CIERA, Northwestern University 1 University of Chicago 2 CIERA Conference: September 2, 2011 1 Type Ia Supernovae Image: NASA

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

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

Eddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):

Eddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr): AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress

More information

centrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration

centrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit

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

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

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

Colloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.

Colloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p. Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš

More information

Numerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model

Numerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model Numerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model Hitoshi Sugiyama 1),Naoto Kato 1), Masahiro Ouchi 1) 1) Graduate School of Engineering,Utsunomiya University

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

Turbulence (January 7, 2005)

Turbulence (January 7, 2005) http://www.tfd.chalmers.se/gr-kurs/mtf071 70 Turbulence (January 7, 2005) The literature for this lecture (denoted by LD) and the following on turbulence models is: L. Davidson. An Introduction to Turbulence

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

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

Dynamics of turbulence under the effect of stratification and internal waves

Dynamics of turbulence under the effect of stratification and internal waves doi:10.5194/npg-22-337-2015 Author(s) 2015. CC Attribution 3.0 License. Dynamics of turbulence under the effect of stratification and internal waves O. A. Druzhinin 1 and L. A. Ostrovsky 2 1 The Institute

More information

Turbulence and Energy Transfer in Strongly-Stratified Flows

Turbulence and Energy Transfer in Strongly-Stratified Flows Turbulence and Energy Transfer in Strongly-Stratified Flows James J. Riley University of Washington Collaborators: Steve debruynkops (UMass) Kraig Winters (Scripps IO) Erik Lindborg (KTH) First IMS Turbulence

More information

AOSC 614. Chapter 4 Junjie Liu Oct. 26, 2006

AOSC 614. Chapter 4 Junjie Liu Oct. 26, 2006 AOSC 614 Chapter 4 Junjie Liu Oct. 26, 2006 The role of parameterization in the numerical model Why do we need parameterization process? Chapter 2 discussed the governing equations. Chapter 3 discussed

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

Dynamics of Atmospheres and Oceans

Dynamics of Atmospheres and Oceans Dynamics of Atmospheres and Oceans 51 (2011) 137 150 Contents lists available at ScienceDirect Dynamics of Atmospheres and Oceans journal homepage: www.elsevier.com/locate/dynatmoce Evaluation of turbulent

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

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

Topics in Fluid Dynamics: Classical physics and recent mathematics

Topics in Fluid Dynamics: Classical physics and recent mathematics Topics in Fluid Dynamics: Classical physics and recent mathematics Toan T. Nguyen 1,2 Penn State University Graduate Student Seminar @ PSU Jan 18th, 2018 1 Homepage: http://math.psu.edu/nguyen 2 Math blog:

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. Atmospheric Sciences 6150

A Simple Turbulence Closure Model. Atmospheric Sciences 6150 A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: V = V + v V = U i + u i Mean velocity: V = Ui + V j + W k =(U, V, W ) U i =(U

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

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

Mostafa Momen. Project Report Numerical Investigation of Turbulence Models. 2.29: Numerical Fluid Mechanics

Mostafa Momen. Project Report Numerical Investigation of Turbulence Models. 2.29: Numerical Fluid Mechanics 2.29: Numerical Fluid Mechanics Project Report Numerical Investigation of Turbulence Models Mostafa Momen May 2015 Massachusetts Institute of Technology 1 Numerical Investigation of Turbulence Models Term

More information

Chapter 6: Rational Expr., Eq., and Functions Lecture notes Math 1010

Chapter 6: Rational Expr., Eq., and Functions Lecture notes Math 1010 Section 6.1: Rational Expressions and Functions Definition of a rational expression Let u and v be polynomials. The algebraic expression u v is a rational expression. The domain of this rational expression

More information

A Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation by Sullivan

A Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation by Sullivan 耶鲁 - 南京信息工程大学大气环境中心 Yale-NUIST Center on Atmospheric Environment A Discussion on The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation

More information

Direct and Large Eddy Simulation of stably stratified turbulent Ekman layers

Direct and Large Eddy Simulation of stably stratified turbulent Ekman layers Direct and Large Eddy Simulation of stably stratified turbulent Ekman layers Stimit Shah, Elie Bou-Zeid Princeton University 64 th APS DFD Baltimore, Maryland Nov 21, 211 Effect of Stability on Atmospheric

More information

CFD Modelling of Compartment Fires

CFD Modelling of Compartment Fires CFD Modelling of Compartment Fires A validation study of a Large Eddy Simulation CFD code using compartment fire test data from Cardington By Nicholas Pope Presentation summary Introduction Fire modelling

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

Exercise 5: Exact Solutions to the Navier-Stokes Equations I

Exercise 5: Exact Solutions to the Navier-Stokes Equations I Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel

More 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

Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein

Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems Institute for Mathematics and Scientific Computing, Austria DISC Summerschool 5 Outline of the talk POD and singular value decomposition

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

Comparison of the simulated upper-ocean vertical structure using 1-dimensional mixed-layer models

Comparison of the simulated upper-ocean vertical structure using 1-dimensional mixed-layer models Comparison of the simulated upper-ocean vertical structure using 1-dimensional mixed-layer models Sonaljit Mukherjee 1 and Amit Tandon 1 1 Department of Mechanical Engineering, University of Massachusetts

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

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

Numerical Simulations of a Stratified Oceanic Bottom Boundary Layer. John R. Taylor - MIT Advisor: Sutanu Sarkar - UCSD

Numerical Simulations of a Stratified Oceanic Bottom Boundary Layer. John R. Taylor - MIT Advisor: Sutanu Sarkar - UCSD Numerical Simulations of a Stratified Oceanic Bottom Boundary Layer John R. Taylor - MIT Advisor: Sutanu Sarkar - UCSD Motivation Objective I: Assess and improve parameterizations of the bottom boundary

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

On stably stratified homogeneous shear flows subjected to rotation

On stably stratified homogeneous shear flows subjected to rotation Center for Turbulence Research Proceedings of the Summer Program 2000 241 On stably stratified homogeneous shear flows subjected to rotation By B. A. Pettersson Reif, A.Ooi AND P. A. Durbin Theoretical

More information

Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.

Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations

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

10. Buoyancy-driven flow

10. Buoyancy-driven flow 10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation

More information

Outline: Phenomenology of Wall Bounded Turbulence: Minimal Models First Story: Large Re Neutrally-Stratified (NS) Channel Flow.

Outline: Phenomenology of Wall Bounded Turbulence: Minimal Models First Story: Large Re Neutrally-Stratified (NS) Channel Flow. Victor S. L vov: Israeli-Italian March 2007 Meeting Phenomenology of Wall Bounded Turbulence: Minimal Models In collaboration with : Anna Pomyalov, Itamar Procaccia, Oleksii Rudenko (WIS), Said Elgobashi

More information

DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT

DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT 10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department

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

Tilting Shear Layers in Coastal Flows

Tilting Shear Layers in Coastal Flows DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Tilting Shear Layers in Coastal Flows Karl R. Helfrich Department of Physical Oceanography, MS-21 Woods Hole Oceanographic

More information

Testing Turbulence Closure Models Against Oceanic Turbulence Measurements

Testing Turbulence Closure Models Against Oceanic Turbulence Measurements Testing Turbulence Closure Models Against Oceanic Turbulence Measurements J. H. Trowbridge Woods Hole Oceanographic Institution Woods Hole, MA 02543 phone: 508-289-2296 fax: 508-457-2194 e-mail: jtrowbridge@whoi.edu

More information

Role of Turbulent Mixing in Plankton Dynamics Simulated by Large Eddy Simulation (LES) YIGN NOH Yonsei University, S. Korea

Role of Turbulent Mixing in Plankton Dynamics Simulated by Large Eddy Simulation (LES) YIGN NOH Yonsei University, S. Korea Role of Turbulent Mixing in Plankton Dynamics Simulated by Large Eddy Simulation (LES) YIGN NOH Yonsei University, S. Korea Application of LES of the ocean mixed layer to 1. the parameterization of the

More information

Atm S 547 Boundary Layer Meteorology

Atm S 547 Boundary Layer Meteorology Lecture 8. Parameterization of BL Turbulence I In this lecture Fundamental challenges and grid resolution constraints for BL parameterization Turbulence closure (e. g. first-order closure and TKE) parameterizations

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

Numerical studies of dispersion due to tidal flow through Moskstraumen, northern Norway

Numerical studies of dispersion due to tidal flow through Moskstraumen, northern Norway Numerical studies of dispersion due to tidal flow through Moskstraumen, northern Norway Lynge B.K. 1,2 Berntsen J. 3 Gjevik B. 1 Department of Mathematics, University of Oslo 1 Norwegian Hydrographic Service

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

Ocean Modeling II. Parameterized Physics. Peter Gent. Oceanography Section NCAR

Ocean Modeling II. Parameterized Physics. Peter Gent. Oceanography Section NCAR Ocean Modeling II Parameterized Physics Peter Gent Oceanography Section NCAR Ocean Modelling Challenges LO-RES (3 o ) O(100+ years/day) WORKHORSE (1 o ) O(10-100 years/day) HI-RES (0.1 o ) O(1 year/day)

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

A simple subgrid-scale model for astrophysical turbulence

A simple subgrid-scale model for astrophysical turbulence CfCA User Meeting NAOJ, 29-30 November 2016 A simple subgrid-scale model for astrophysical turbulence Nobumitsu Yokoi Institute of Industrial Science (IIS), Univ. of Tokyo Collaborators Axel Brandenburg

More information

+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2):

+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2): International Conference on Methods of Aerophysical Research, ICMAR 008 LARGE EDDY SIMULATION OF TURBULENT ROUND IMPINGING JET B.B. Ilyushin, D.V. Krasinsky Kutateladze Institute of Thermophysics SB RAS

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

Probability density function (PDF) methods 1,2 belong to the broader family of statistical approaches

Probability density function (PDF) methods 1,2 belong to the broader family of statistical approaches Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,

More information

The Simulation of Wraparound Fins Aerodynamic Characteristics

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

Periodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.

Periodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain. Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable

More information

Turbulence. 2. Reynolds number is an indicator for turbulence in a fluid stream

Turbulence. 2. Reynolds number is an indicator for turbulence in a fluid stream Turbulence injection of a water jet into a water tank Reynolds number EF$ 1. There is no clear definition and range of turbulence (multi-scale phenomena) 2. Reynolds number is an indicator for turbulence

More information

Introduction to Turbulence AEEM Why study turbulent flows?

Introduction to Turbulence AEEM Why study turbulent flows? Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and

More information

Large Eddy Simulation of Stratified Mixing in a Three-Dimensional Lock-Exchange System

Large Eddy Simulation of Stratified Mixing in a Three-Dimensional Lock-Exchange System Large Eddy Simulation of Stratified Mixing in a Three-Dimensional Lock-Exchange System by Tamay M. Özgökmen1, Traian Iliescu 2 and Paul F. Fischer 3 1. RSMAS/MPO, University of Miami, Miami, Florida 2.

More information

Computers and Mathematics with Applications. Investigation of the LES WALE turbulence model within the lattice Boltzmann framework

Computers and Mathematics with Applications. Investigation of the LES WALE turbulence model within the lattice Boltzmann framework Computers and Mathematics with Applications 59 (2010) 2200 2214 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Investigation

More information

Numerical simulations of heat transfer in plane channel flow

Numerical simulations of heat transfer in plane channel flow Numerical simulations of heat transfer in plane channel flow Najla EL GHARBI 1, 3, a, Rafik ABSI 2, b and Ahmed BENZAOUI 3, c 1 Renewable Energy Development Center, BP 62 Bouzareah 163 Algiers, Algeria

More information

Lattice Boltzmann Method for Fluid Simulations

Lattice Boltzmann Method for Fluid Simulations Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool

More 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