Statistical turbulence modelling in the oceanic mixed layer
|
|
- Gavin Davidson
- 5 years ago
- Views:
Transcription
1 Statistical turbulence modelling in the oceanic mixed layer Hans Burchard Baltic Sea Research Institute Warnemünde, Germany Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 1/44
2 Program of presentation NSE Two-equation closure models Two-equation closure models General Ocean Turbulence Model (GOTM) Examples (observations versus simulations) General Estuarine Transport Model (GETM) Three-dimensional model study Conclusions Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 2/44
3 Basic approaches Two approaches to averaged turbulence modelling: Statistical turbulence modelling: Convert NSE to Friedmann-Keller series, cut-off where suitable and parameterise unknown terms. Empirical turbulence modelling: Close equations on lowest order and parameterise relevant processes. Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 3/44
4 NSE Reynolds equation Momentum Equation: t v i + v j j v i ν jj v i + 2ε ijl Ω j v l = ip ρ g i ρ ρ. Reynolds averaging: Reynolds Equation: v i = v i + ṽ i, ρ = ρ + ρ,... t v i + v j j v i j (ν j v i ṽ j ṽ i )+2ε ijl Ω j v l = i p ρ g i ρ ρ. Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 4/44
5 Reynolds Stress Equation ) t ṽ i ṽ j + l ( v l ṽ i ṽ j + ṽ l ṽ i ṽ j ν l ṽ i ṽ j = l v i ṽ l ṽ j l v j ṽ l ṽ i }{{} P ij 2Ω l (ε ilm ṽ j ṽ m + ε jlm ṽ i ṽ m ) }{{} Ω ij 1 ρ {g i ṽ j ρ + g j ṽ i ρ } }{{} B ij 2ν ( l ṽ j )( l ṽ i ). }{{} ε ij 1 ρ ( ṽ i j p + ṽ j i p ) }{{} Π ij Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 5/44
6 Algebraic SMCs The following steps lead to different types of second-moment closures: Empirical closures of pressure-strain correlators. Neglect or simplification of advective and diffusive fluxes of second-moments. Neglect of rotational terms in the second-moment equations. Boundary layer assumption (neglect of horizontal gradients and non-hydrostatic effects).... and many more details... Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 6/44
7 Algebraic SMCs Turbulent Fluxes: ũ w = ν t z ū, w T = ν t z T Eddy Viscosity / Eddy Diffusivity: ν t = c µ (α M, α N ) k2 ε, ν t = c µ(α M, α N ) k2 ε. Shear Number, Buoyancy Number: α M = k2 ε 2 M 2, α N = k2 ε 2 N 2. Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 7/44
8 .22.2 Stability Functions Canuto et al. [21]: c µ c µ α M α M α N α N Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 8/44
9 Exact TKE-Equation ) 1 t k + j ( v j k + ṽ j 2ṽ2 i ν jk + 1 ρ ṽ j p = ṽ j ṽ i i v j }{{} P g ρ ṽ 3 ρ }{{} B ν ( j ṽ i ) 2, }{{} ε This TKE equation will be modelled as it is given above, the only parameterisations needed are for the turbulent f lux terms, for which usually the down-gradient approximation is used. Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 9/44
10 Dissipation equation Exact form (e.g. Wilcox [1998]): ) t ε + j ( v j ε + ṽ j ν( j ṽ i ) 2 ν j ε + 2 ν ρ i ṽ j i p ) = 2ν j v i ( i ṽ k j ṽ k + k ṽ i k ṽ j 2ν jk v i ṽ k j ṽ i }{{} 2ν g ρ j ṽ 3 j ρ }{{} B ε 2ν P ε ( j ṽ i k ṽ i j ṽ k + ν ( ij ṽ k ) 2 ) } {{ } ε ε (1) k-ε model (Launder and Spalding [1972]): ( ) νt t ε z z ε = ε σ ε k (c ε1p + c ε3 B c ε2 ε). Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 1/44
11 Mellor-Yamada model General relation between k, ε and L: L = c 3/4 µ k 3/2 k-kl model (Mellor and Yamada [1982]): t (kl) z (S l z (kl)) = ε L 2 [ E 1 P + E 3 B (1 + E 2 ( LLz ) 2 ) ε ]. Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 11/44
12 Length scale equations (cont d) Other approaches are using equations for ω = ε/k (k-ω model), k/ε, k 2 /ε,..., so why not using the generalised approach of a k n ε m equation? Generic length scale equation (Umlauf and Burchard [23]): t (k n ε m ) z ( νt σ nm k n ε m ) = k n 1 ε m (c nm1 P + c nm3 B c nm2 ε) This works without correction term only for m. Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 12/44
13 Total equilibrium (k-ε) R st i R i = R st i = c µ c µ c2ε c 1ε c 2ε c 3ε..25: Steady-state Richardson number. PSfrag replacements R i st KC RH CA CB c 3ε Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 13/44
14 Kato-Phillips experiment 45 MLD for Canuto et al. version A MLD / m Time / h Empirical Ri_st =.2 Ri_st =.3 Ri_st =.4 Ri_st =.6 Ri_st =.8 Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 14/44
15 Total equilibrium (k-kl) R st i R i = R st i = K M K H 1 E 1 1 E 3..25: Steady-state Richardson number R st i B/ε G H.2.15 PSfrag replacements E 3 Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 15/44
16 GOTM, Challenge Aim The Idea Key features Software Fortran code Test cases Forcing How to run? Information What's New Publications list FAQ User Group Hot Links Who's Who? Guestbook GOTM is a one-dimensional numerical model developed and supported by a core team of ocean modellers. GOTM aims at simulating accurately vertical exchange processes in the marine environment where mixing is known to play a key role. GOTM is freely available under the GPL (Gnu Public License). The interested user can download the source code, a set of test cases (Papa, November, Flex,...) and a comprehensive report. You are warmly invited to join the GOTM mailing list and send any comments/questions to the GOTM team or become a GOTM contributor. The GOTM developers are grateful to their sponsors. Page " maintained by webmaster. Last update: 1/28/ 18:1:2 PPM: 1/26/ 2:48:54 Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 16/44
17 Turbulence under waves How to compare these observations? How to simulate them??? Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 17/44
18 Breaking waves: Simulations eplacements 1 (z z)/hs 1 Terray et al. [1996] Drennan et al. [1996] Anis and Moum et al. [1995] Numerical, z /H s = 2. Numerical, z /H s = 1. Numerical, z /H s =.5 Log-Law No shear production ε / c w u 3 /H s Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 18/44
19 ũ 2 /(w ) 2 ũ 2 /(w ) 2 w 2 Free /(w ) 2 Convection w 2 /(w ) 2 T 2 /(T 2 ) T 2 /(T 2 ) Temperature ε/b KC ε/b RH -.2 CA CB Temperature ( LES -.4 emperature Flux T T max )/T -.6 Variance of w u T /(w T Variance ) of u -.8 Variance of w ũ 2 /(w -1 Variance ) 2 of w Variance of T w /(wvariance ) 2 of T Dissipation Rate Dissipation ( T T max )/T T 2 /(T 2 Rate ) z/dm Temperature Temperature Flux Variance of u Variance of w Variance of T z/dm z/dm Dissipation Rate Temperature Flux -.8 KC RH -1 CA CB LES w T /(w T ) -.8 KC RH -1 CA CB LES ε/b Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 19/44
20 Geographical overview Y Lago Maggiore Northern North Sea Liverpool Bay OWS Papa Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 2/44
21 Lago Maggiore, Italy Observations and simulations of T and ε (Stips et al. [22]) Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 21/44
22 Open Ocean: OWS Papa Temperature Observations Depth [m] Measurements OWS Papa 1961/ Julian Day 1961/62 T [ o C] Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 22/44
23 Open Ocean: OWS Papa Heat budget OWS Papa 1961/62 2e+9 Heat from T-profiles Heat from surface fluxes eplacements Q / J m 2 1.5e+9 1e+9 5e Julian Day 1961/62 Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 23/44
24 Open Ocean: OWS Papa Temperature Simulations: k-ε model with alg. SMC Depth [m] Simulation OWS Papa 1961/ Julian Day 1961/62 T [ o C] Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 24/
25 Open Ocean: OWS Papa Temperature profiles T-profiles, day 12 T-profiles, day 15 T-profiles, day Observed Simulated, KC Simulated, CA Observed Simulated, KC Simulated, CA Observed Simulated, KC Simulated, CA T-profiles, day 21 T-profiles, day 24 T-profiles, day Observed Simulated, KC Simulated, CA Observed Simulated, KC Simulated, CA Observed Simulated, KC Simulated, CA Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 25/44
26 Open Ocean: OWS Papa SST in deg C SST in deg C Sea Surface Temperature OWS Papa 1961/62 Observation Simulation KC Julian Day 1961/62 Observation Simulation CA Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 26/44
27 Northern North Sea Bathymetry and station map ' 45' 1 ' 1 15' 1 3' 1 45' 2 ' 59 45' 59 45' 61 1 Norway ' 59 3' Scotland ' 59 15' ' 59 ' 3' 45' 1 ' 1 15' 1 3' 1 45' 2 ' Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 27/44
28 Northern North Sea Observed temperature Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 28/44
29 Northern North Sea Simulated temperature Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 29/44
30 Northern North Sea Wind and Tides rag replacements Q / (W m 2 ) ( x ζ) 2 + ( y ζ) 2 rag replacements Q / (W m 2 ) ( x ζ) 2 + ( y ζ) 2 τ / (N m 2 ) τ / (N m 2 ) /1 12: ADCP Model Surface stress at station NNS Bed stress at station NNS 22/1 : 22/1 12: Date in /1 : 23/1 12: Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 3/44
31 Northern North Sea Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 31/44
32 NPZD model in GOTM Temperature Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 32/44
33 NPZD model in GOTM Nutrients Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 33/44
34 NPZD model in GOTM Phytoplankton Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 34/44
35 NPZD model in GOTM Zooplankton Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 35/44
36 NPZD model in GOTM Detritus Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 36/44
37 Liverpool Bay Section of Temperature and Salinity Fig 2.a Temperature (Degrees C) 1 6 Depth (m) -2-4 LB Fig 2.b Salinity (PSU) Depth (m) Rippeth, Fisher, Simpson [21] Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 37/44
38 Liverpool Bay Observed and simulated temperature and salinity Simpson, Burchard, Fisher, Rippeth [22] Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 38/44
39 Liverpool Bay Observed and simulated current velocity Simpson, Burchard, Fisher, Rippeth [22] Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 39/44
40 Liverpool Bay Observed and simulated dissipation rates Simpson, Burchard, Fisher, Rippeth [22] Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 4/44
41 GETM General Estuarine Transport Model (Burchard & Bolding [22]) General vertical coordinates Horizontal curvilinear coordinates Uses GOTM as turbulence model Drying & fl ooding Public Domain Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 41/44
42 North Sea simulation PSfrag Sea surface temperature replacements on Aug. 1, PSfrag replacements surface temperature 4 on Aug. 1, 1997 surface T / replacements temperature on Aug. 1, 1997 C 2 3 replacements 2 Aug. 1, T / C Sea surface y / nmsalinity on Aug. 1, x / nm T / C S/psu Temperature section on Aug. 1, 1997 Aug. 1, 1997 T / -1 2 S/psu C Temperature section -2 on Aug. 1, T / C5 y / nm z / m -5-6 y / nm y / nm T / C z / m Sea surface salinity on Aug. 1, x / nm Salinity section on Aug. 1, 1997 S/psu S/psu x / nm x / nm Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 42/44
43 North Sea / Baltic Sea SSS on June 1, 1997, x = y = 3 nm g replacements y / nm S / psu x / nm Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 43/44
44 Conclusions Reynolds decomposition provides a physically sound framework for developing turbulence models. Empirical parameters allow for calibrating such models to the real world. With such models, turbulence is described but not better understood. For mixed layers, two-equation models with algebraic SMCs are sufficient for reproducing observations. These models are economic enough for using them in 3D models for oceanic applications. Seminar at Bjerknes Centre, October 16, 23, Bergen, Norway p. 44/44
GOTM goes biogeochemical: Application to Gotland Sea
GOTM goes biogeochemical: Application to Gotland Sea Hans Burchard 1,2, Karsten Bolding 3, and Thomas Neumann 2 hans.burchard@io-warnemuende.de 1. Bolding & Burchard Hydrodynamics, Rostock, Germany 2.
More informationSCIENCE CHINA Earth Sciences. Turbulence and mixing in a freshwater-influenced tidal bay: Observations and numerical modeling
SCIENCE CHINA Earth Sciences RESEARCH PAPER November 2015 Vol.58 No.11: 2049 2058 doi: 10.1007/s11430-015-5093-7 Turbulence and mixing in a freshwater-influenced tidal bay: Observations and numerical modeling
More informationWQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
More informationImpact of Typhoons on the Western Pacific Ocean (ITOP) DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Impact of Typhoons on the Western Pacific Ocean (ITOP) DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment
More information6A.3 Stably stratified boundary layer simulations with a non-local closure model
6A.3 Stably stratified boundary layer simulations with a non-local closure model N. M. Colonna, E. Ferrero*, Dipartimento di Scienze e Tecnologie Avanzate, University of Piemonte Orientale, Alessandria,
More informationComparison 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 informationOn 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 informationComparison between vertical shear mixing and surface wave-induced mixing in the global ocean
Comparison between vertical shear mixing and surface wave-induced mixing in the global ocean Fangli Qiao and Chuanjiang Huang Key Laboratory of Marine Science and Numerical Modeling First Institute of
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationSurface wind mixing in the Regional Ocean Modeling System (ROMS)
DOI 10.1186/s40562-017-0090-7 RESEARCH LETTER Open Access Surface wind mixing in the Regional Ocean Modeling System (ROMS) Robin Robertson * and Paul Hartlipp Abstract Mixing at the ocean surface is key
More informationmeters, 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 informationEddy 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 informationModelling turbulent vertical mixing sensitivity using a 1-D version of NEMO
doi:10.5194/gmd-8-69-2015 Author(s) 2015. CC Attribution 3.0 License. Modelling turbulent vertical mixing sensitivity using a 1-D version of NEMO G. Reffray 1, R. Bourdalle-Badie 1, and C. Calone 2 1 Mercator
More informationMellor-Yamada Level 2.5 Turbulence Closure in RAMS. Nick Parazoo AT 730 April 26, 2006
Mellor-Yamada Level 2.5 Turbulence Closure in RAMS Nick Parazoo AT 730 April 26, 2006 Overview Derive level 2.5 model from basic equations Review modifications of model for RAMS Assess sensitivity of vertical
More informationTutorial 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 informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More informationA 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 informationA 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 informationModeling the influence of wave-enhanced turbulence in a shallow tide- and wind-driven water column
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jc004246, 2008 Modeling the influence of wave-enhanced turbulence in a shallow tide- and wind-driven water column
More informationESS Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1
ESS5203.03 - Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1 Text: J.R.Garratt, The Atmospheric Boundary Layer, 1994. Cambridge Also some material from J.C. Kaimal and
More informationA two-fluid model of turbulent two-phase flow for simulating turbulent stratified flows
Ocean Engineering 30 (2003) 153 161 www.elsevier.com/locate/oceaneng A two-fluid model of turbulent two-phase flow for simulating turbulent stratified flows Y.M. Shen a,, C.-O. Ng b, A.T. Chwang b a State
More informationAsymmetric Tidal Mixing due to the Horizontal Density Gradient*
418 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 38 Asymmetric Tidal Mixing due to the Horizontal Density Gradient* MING LI Horn Point Laboratory, University of Maryland Center for
More informationAnalysis of Near-Surface Oceanic Measurements Obtained During CBLAST-Low
Analysis of Near-Surface Oceanic Measurements Obtained During CBLAST-Low John H. Trowbridge Woods Hole Oceanographic Institution, MS#12, Woods Hole, MA 02543 phone: (508) 289-2296 fax: (508) 457-2194 email:
More informationModelling the Channel and the Bay of Biscay using the MARS model. Contributions: Dyneco Physed / Pelagos / Benthos EMH
Modelling the Channel and the Bay of Biscay using the MARS model Contributions: Dyneco Physed / Pelagos / Benthos EMH Overview 1. The MANGA Configuration 2. Model validation: Sea Surface Temperature 3.
More informationESCI 485 Air/sea Interaction Lesson 2 Turbulence Dr. DeCaria
ESCI 485 Air/sea Interaction Lesson Turbulence Dr. DeCaria References: Air-sea Interaction: Laws and Mechanisms, Csanady An Introduction to Dynamic Meteorology ( rd edition), J.R. Holton An Introduction
More informationTurbulence - 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 informationPUBLICATIONS. Journal of Geophysical Research: Oceans. Upper-ocean mixing due to surface gravity waves RESEARCH ARTICLE 10.
PUBLICATIONS Journal of Geophysical Research: Oceans RESEARCH ARTICLE Key Points: Nonbreaking waves and Langmuir turbulence are the most important terms from waves on ocean mixing 2-D wave spectrum should
More informationFrictional 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ρ 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 informationSimulations 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 informationPOLCOMS Metadata for the ARCoES project Keywords: POLCOMS, WAM, residual circulation, waves, Liverpool Bay, UK shelf
POLCOMS Metadata for the ARCoES project Keywords: POLCOMS, WAM, residual circulation, waves, Liverpool Bay, UK shelf POLCOMS is the Proudman Oceanographic Laboratory Coastal Ocean Modelling System. It
More informationIntroduc)on to modeling turbulence in oceanic flows
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 The biggest
More informationTurbulence 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 informationOCEAN HYDRODYNAMIC MODEL
Jurnal Teknologi Pengelolaan Limbah (Journal of Waste Management Technology), ISSN 1410-9565 Volume 10 Nomor 1 Juli 2007 (Volume 10, Number 1, July, 2007) Pusat Teknologi Limbah Radioaktif (Radioactive
More informationarxiv: v1 [physics.flu-dyn] 11 Oct 2012
Low-Order Modelling of Blade-Induced Turbulence for RANS Actuator Disk Computations of Wind and Tidal Turbines Takafumi Nishino and Richard H. J. Willden ariv:20.373v [physics.flu-dyn] Oct 202 Abstract
More informationThe mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient):
8. TURBULENCE MODELLING 1 SPRING 2019 8.1 Eddy-viscosity models 8.2 Advanced turbulence models 8.3 Wall boundary conditions Summary References Appendix: Derivation of the turbulent kinetic energy equation
More informationValidation of the Hydrodynamic Module of a 3D Ecological Model for the Ria de Aveiro Lagoon (Portugal)
International Congress on Environmental Modelling and Software Brigham Young University BYU ScholarsArchive th International Congress on Environmental Modelling and Software - Barcelona, Catalonia, Spain
More informationToward Accurate Coastal Ocean Modeling
Toward Accurate Coastal Ocean Modeling Peter C. Chu Naval Postgraduate School Monterey, CA 93943, USA Email: pcchu@nps.edu http://www.oc.nps.navy.mil/~chu International Council for Sciences, Scientific
More informationVertical resolution of numerical models. Atm S 547 Lecture 8, Slide 1
Vertical resolution of numerical models Atm S 547 Lecture 8, Slide 1 M-O and Galperin stability factors Atm S 547 Lecture 8, Slide 2 Profile vs. forcing-driven turbulence parameterization Mellor-Yamada
More informationContents. 1 Motivation 1
Contents I Contents 1 Motivation 1 2 Stratification in the Limfjord 2 2.1 About the Limfjord............................ 2 2.2 Production of stratification....................... 3 2.2.1 Solar heating...........................
More informationLecture 9: Tidal Rectification, Stratification and Mixing
Lecture 9: Tidal Rectification, Stratification and Mixing Chris Garrett 1 Additional Notes on Tidal Rectification This lecture continues the discussion of long-wavelength tidal flow over comparatively
More informationcentrifugal 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 informationPotential impact of offshore wind farms on water renewal in the Baltic Sea
Potential impact of offshore wind farms on water renewal in the Baltic Sea Hans Burchard hans.burchard@io-warnemuende.de Baltic Sea Research Institute Warnemünde, Germany Summer School, AWI-GKSS-IOW, Büsum,
More informationAppendix: Nomenclature
Appendix: Nomenclature 209 Appendix: Nomenclature Chapter 2: Roman letters A amplitude of surface wave m A j (complex) amplitude of tidal component [various] a right ascension rad or a magnitude of tidal
More informationColloquium 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 informationA 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 informationFRICTION-DOMINATED WATER EXCHANGE IN A FLORIDA ESTUARY
FRICTION-DOMINATED WATER EXCHANGE IN A FLORIDA ESTUARY By KIMBERLY ARNOTT A THESIS PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE
More informationOCN/ATM/ESS 587. Ocean circulation, dynamics and thermodynamics.
OCN/ATM/ESS 587 Ocean circulation, dynamics and thermodynamics. Equation of state for seawater General T/S properties of the upper ocean Heat balance of the upper ocean Upper ocean circulation Deep circulation
More informationNONLINEAR 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 informationChapter (3) TURBULENCE KINETIC ENERGY
Chapter (3) TURBULENCE KINETIC ENERGY 3.1 The TKE budget Derivation : The definition of TKE presented is TKE/m= e = 0.5 ( u 2 + v 2 + w 2 ). we recognize immediately that TKE/m is nothing more than the
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 informationAtm 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 informationThe applicability of Monin Obukhov scaling for sloped cooled flows in the context of Boundary Layer parameterization
Julia Palamarchuk Odessa State Environmental University, Ukraine The applicability of Monin Obukhov scaling for sloped cooled flows in the context of Boundary Layer parameterization The low-level katabatic
More informationBefore 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 informationNumerical 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 informationChapter 24. Turbulence and heat flux formulations Test case pycno Theory
Chapter 24 Turbulence and heat flux formulations 24.1 Test case pycno 24.1.1 Theory This test case is a 1-D problem describing the evolution of a wind-driven surface mixed layer in the absence of rotation.
More informationTURBULENT KINETIC ENERGY
TURBULENT KINETIC ENERGY THE CLOSURE PROBLEM Prognostic Moment Equation Number Number of Ea. fg[i Q! Ilial.!.IokoQlI!!ol Ui au. First = at au.'u.' '_J_ ax j 3 6 ui'u/ au.'u.' a u.'u.'u k ' Second ' J =
More informationSMS 303: Integrative Marine
SMS 303: Integrative Marine Sciences III Instructor: E. Boss, TA: A. Palacz emmanuel.boss@maine.edu, 581-4378 5 weeks & topics: diffusion, mixing, tides, Coriolis, and waves. Pre-class quiz. Mixing: What
More informationExplicit 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 informationOCEAN MODELING. Joseph K. Ansong. (University of Ghana/Michigan)
OCEAN MODELING Joseph K. Ansong Email: jkansong@umich.edu (University of Ghana/Michigan) OUTLINE INTRODUCTION MOTIVATION EQUATIONS OF MOTION INTRODUCTION TO ROMS (Regional Ocean Modeling System) EXAMPLES
More informationAssessment of the performance of a turbulence closure model: along the tidally-influenced Kaipara River to the estuary, NZ
Assessment of the performance of a turbulence closure model: along the tidally-influenced Kaipara River to the estuary, NZ Berengere S. Dejeans 1, Julia C. Mullarney 2, Iain T. MacDonald 3 and Glen M.
More informationOcean Dynamics. The Great Wave off Kanagawa Hokusai
Ocean Dynamics The Great Wave off Kanagawa Hokusai LO: integrate relevant oceanographic processes with factors influencing survival and growth of fish larvae Physics Determining Ocean Dynamics 1. Conservation
More informationEXAMPLE SHEET FOR TOPIC 3 AUTUMN 2013
EXAMPLE SHEET FOR TOPIC ATMN 01 Q1. se dimensional analysis to investigate how the capillary rise h of a liquid in a tube varies with tube diameter d, gravity g, fluid density ρ, surface tension σ and
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationNeeds work : define boundary conditions and fluxes before, change slides Useful definitions and conservation equations
Needs work : define boundary conditions and fluxes before, change slides 1-2-3 Useful definitions and conservation equations Turbulent Kinetic energy The fluxes are crucial to define our boundary conditions,
More informationLogarithmic velocity profile in the atmospheric (rough wall) boundary layer
Logarithmic velocity profile in the atmospheric (rough wall) boundary layer P =< u w > U z = u 2 U z ~ ε = u 3 /kz Mean velocity profile in the Atmospheric Boundary layer Experimentally it was found that
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationInteraction(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 informationFundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.
Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer
More informationTesting 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 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 informationOn the Velocity Gradient in Stably Stratified Sheared Flows. Part 2: Observations and Models
Boundary-Layer Meteorol (2010) 135:513 517 DOI 10.1007/s10546-010-9487-y RESEARCH NOTE On the Velocity Gradient in Stably Stratified Sheared Flows. Part 2: Observations and Models Rostislav D. Kouznetsov
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationGeneral Curvilinear Ocean Model (GCOM): Enabling Thermodynamics
General Curvilinear Ocean Model (GCOM): Enabling Thermodynamics M. Abouali, C. Torres, R. Walls, G. Larrazabal, M. Stramska, D. Decchis, and J.E. Castillo AP0901 09 General Curvilinear Ocean Model (GCOM):
More informationLecture 2. Turbulent Flow
Lecture 2. Turbulent Flow Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of this turbulent water jet. If L is the size of the largest eddies, only very small
More informationSynthesis of TKE Dissipation Rate Observations in the Ocean's Wave Zone
Synthesis of TKE Dissipation Rate Observations in the Ocean's Wave Zone Ayal Anis Department of Oceanography Texas A&M University Galveston, TX 77551 phone: (409) 740-4987 fax: (409) 740-4786 email: anisa@tamug.tamu.edu
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 informationNumerical 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 informationA suite of boundary conditions for the simulation of atmospheric ows using a modied RANS/k ε closure
6th OpenFOAM Workshop Penn State, State College, USA A suite of boundary conditions for the simulation of atmospheric ows using a modied RANS/k ε closure J. Sumner, C. Masson École de technologie supérieure
More informationWind driven mixing below the oceanic mixed layer
Wind driven mixing below the oceanic mixed layer Article Published Version Grant, A. L. M. and Belcher, S. (2011) Wind driven mixing below the oceanic mixed layer. Journal of Physical Oceanography, 41
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More informationNorthern Arabian Sea Circulation Autonomous Research (NASCar) DRI: A Study of Vertical Mixing Processes in the Northern Arabian Sea
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Northern Arabian Sea Circulation Autonomous Research (NASCar) DRI: A Study of Vertical Mixing Processes in the Northern
More informationAccretion 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 informationGoverning Equations for Turbulent Flow
Governing Equations for Turbulent Flow (i) Boundary Layer on a Flat Plate ρu x Re x = = Reynolds Number µ Re Re x =5(10) 5 Re x =10 6 x =0 u/ U = 0.99 層流區域 過渡區域 紊流區域 Thickness of boundary layer The Origin
More informationTurbulent 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 informationDynamics of the Ems Estuary
Dynamics of the Ems Estuary Physics of coastal systems Jerker Menninga 0439738 Utrecht University Institute for Marine and Atmospheric research Utrecht Lecturer: Prof. dr. H.E. de Swart Abstract During
More informationMass 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 informationLecture 2 Atmospheric Boundary Layer
Lecture Atmospheric Boundary Layer H. J. Fernando Ariona State niversity Region of the lower atmosphere where effects of the Earth surface are felt Surface fluxes of momentum, buoyancy.. Neutral, Convective,
More informationTurbulence 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( ) = 1005 J kg 1 K 1 ;
Problem Set 3 1. A parcel of water is added to the ocean surface that is denser (heavier) than any of the waters in the ocean. Suppose the parcel sinks to the ocean bottom; estimate the change in temperature
More informationBuoyancy 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 informationThe atmospheric boundary layer: Where the atmosphere meets the surface. The atmospheric boundary layer:
The atmospheric boundary layer: Utrecht Summer School on Physics of the Climate System Carleen Tijm-Reijmer IMAU The atmospheric boundary layer: Where the atmosphere meets the surface Photo: Mark Wolvenne:
More informationOverview of HYCOM activities at SHOM
Overview of HYCOM activities at SHOM Stéphanie Louazel, Stéphanie Corréard, Rémy Baraille, Annick Pichon, Cyril Lathuilière, Audrey Pasquet, Emeric Baquet LOM2015 Copenhagen 2 nd 4 th June 2015 French
More informationA numerical model of the circulation in the NW Mediterranean
A numerical model of the circulation in the NW Mediterranean A.M.Doglioli A.A.Petrenko, Z.Hu, M.Kersalé, F.Nencioli, I.Dekeyser, C.Estournel, P.Marsaleix Workshop MEUST Marseille, 8-9 November 2010 Marsaleix
More informationTHE EFFECT OF STRATIFICATION ON THE ROUGHNESS LENGTH AN DISPLACEMENT HEIGHT
THE EFFECT OF STRATIFICATION ON THE ROUGHNESS LENGTH AN DISPLACEMENT HEIGHT S. S. Zilitinkevich 1,2,3, I. Mammarella 1,2, A. Baklanov 4, and S. M. Joffre 2 1. Atmospheric Sciences,, Finland 2. Finnish
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationPhysics of Aquatic Systems
Physics of Aquatic Systems 4. Turbulent Transport in Surface Waters Contents of Session 4: Transport 4.1 Turbulent (eddy) fusion 4. Measurement of eddy fusivity in lakes Tracer methods (artificial tracer)
More informationNumerical modeling of an estuary: a comprehensive skill assessment
Numerical modeling of an estuary: a comprehensive skill assessment John C. Warner U.S. Geological Survey, Woods Hole, Massachusetts, U.S.A. W. Rockwell Geyer Woods Hole Oceanographic Institution, Woods
More informationCFD analysis of the transient flow in a low-oil concentration hydrocyclone
CFD analysis of the transient flow in a low-oil concentration hydrocyclone Paladino, E. E. (1), Nunes, G. C. () and Schwenk, L. (1) (1) ESSS Engineering Simulation and Scientific Software CELTA - Rod SC-41,
More informationNumerical 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