A suite of boundary conditions for the simulation of atmospheric ows using a modied RANS/k ε closure
|
|
- Paula Hines
- 6 years ago
- Views:
Transcription
1 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 June 14th / 38
2 Introduction Wind farm power performance 2 / 38
3 Introduction Wind farm power performance 3 / 38
4 Introduction Oshore wind farms 4 / 38
5 Objective Derive and implement a general k ε surface layer model valid for both stable and neutral thermal stratication 5 / 38
6 Similarity theory Potential temperature proles Neutral Stable z [m] Θ [ C] Figure: Comparison neutral and stable surface layer 6 / 38
7 Similarity theory Velocity proles Neutral Stable z [m] 40 l m = κz/φ m l m = κz U [m/s] Figure: Comparison neutral and stable surface layer 7 / 38
8 Model development Imposing a new length scale The turbulent length scale is determined by the ε equation ( ) d νt dε ε 2 + Pε Cε2 = 0 dz σ ε dz k How to modify ε equation for desired mixing length and turbulence properties? 8 / 38
9 Model development The Apsley & Castro model Recast the ε equation as [ Pε = When l m << l max, When l m l max, Cε1 + (Cε2 Cε1) l m l max Π k ε Pε Cε1 k ε 2 Pε Cε2 0 k ] Πk ε k 9 / 38
10 Model development The Apsley & Castro model Recast the ε equation as [ Pε = When l m << l max, When l m l max, Cε1 + (Cε2 Cε1) l m l max Π k ε Pε Cε1 k ε 2 Pε Cε2 0 k ] Πk ε k 10 / 38
11 Model development The Apsley & Castro model Recast the ε equation as [ Pε = When l m << l max, When l m l max, Cε1 + (Cε2 Cε1) l m l max Π k ε Pε Cε1 k ε 2 Pε Cε2 0 k ] Πk ε k 11 / 38
12 Similarity theory Analytical proles Considering an arbitrary desired l m : Turbulence properties: U(z) = u z z 0 1 l m dz z 1 Θ(z) = σ θ θ dz z 0 l m k(z) = u2 Cµ ε(z) = u3 l ε 12 / 38
13 Model development A revised Apsley & Castro model Assuming equilibrium and writing correction as F d dz ( ) νt dε + F (Cε2 Cε1) ε2 = 0 σ ε dz k Substituting ν t = u l m, ε = u 3 /l ε and k = u 2 / Cµ and solving for F [ F = 1 dl m dl ε d 2 ( l ε κ 2 + l m dz dz dz 2 2l m dl ε l ε dz ) 2 ] 13 / 38
14 Model development A revised Apsley & Castro model Assuming equilibrium and writing correction as F d dz ( ) νt dε + F (Cε2 Cε1) ε2 = 0 σ ε dz k Substituting ν t = u l m, ε = u 3 /l ε and k = u 2 / Cµ and solving for F [ F = 1 dl m dl ε d 2 ( l ε κ 2 + l m dz dz dz 2 2l m dl ε l ε dz ) 2 ] 14 / 38
15 Model development Specic form Taking the mixing length as κz l m = 1 + κz/l max and relating the length scales with a simplied TKE budget equation 1 = 1 1 l ε l m κl yields F = ( 1 + R f ( )( ) 2 lm l ) 3 m l max l max where l m = u 3 0 /P k and R f G k/π k 15 / 38
16 Model development A caveat Turbulent viscosity For analytical proles be solutions, ν t = u l ε u l m k 2 l m k 2 ν t = Cµ = Cµ ε l ε ε (1 R f ) Note: (1 Rf ) is exactly equal to the ratio of length scales for the simplied TKE budget equation 16 / 38
17 Model development Turbulence transport equations 1D k ε model d dz d dz ( ν t dk ( νt dε σ ε dz dz ) + Π k ε = 0 ) ε 2 + Pε Cε2 = 0 k where Π k = P k + G k ( = ν t P k u w du dz ) 2 du dz G k αgθ w = ν t dθ αg σ θ dz 17 / 38
18 OpenFOAM implementation kepsilonlengthlimited // Shear production Gk_ = nut_*2*magsqr(symm(fvc::grad(u_))); // Link to theta field volscalarfield& theta = const_cast<volscalarfield&> (db().lookupobject<volscalarfield>("theta")); // Buoyancy production Gb_ = (alpha_*nut_/sigmatheta_) * g_ & fvc::grad(theta); // Total production volscalarfield G("RASModel::G", Gk_ + Gb_); 18 / 38
19 OpenFOAM implementation kepsilonlengthlimited // Calculate velocity scale u0_ = Foam::pow(Cmu_,0.25) * Foam::sqrt(k_); // Calculate mixing length lm_ = Foam::pow(u0_,3.0) / max(gk_, epsilonsmall_); // Update epsilon, G, Gk, Gb, and lm at the wall epsilon_.boundaryfield().updatecoeffs(); // Flux Richardson number volscalarfield Rf = min(1.0, -Gb_ / max(gk_, epsilonsmall_)); volscalarfield Rfprime = min(1.0, -Gb_ / max(g, epsilonsmall_)); // Weighting function f_ = Foam::sqr(1.0+Rfprime) * (lm_/lmax_+1.0) * Foam::pow((1.0-lm_/lMax_),3.0); 19 / 38
20 OpenFOAM implementation kepsilonlengthlimited // Dissipation equation tmp<fvscalarmatrix> epseqn ( fvm::ddt(epsilon_) + fvm::div(phi_, epsilon_) - fvm::sp(fvc::div(phi_), epsilon_) - fvm::laplacian(depsiloneff(), epsilon_) == (C1_ + (C2_ - C1_)*f_)*G*epsilon_/k_ - fvm::sp(c2_*epsilon_/k_, epsilon_) );... // Re-calculate viscosity nut_ = Cmu_*sqr(k_)/epsilon_ * max(1.0 - Rf, 1e-10); nut_.correctboundaryconditions(); 20 / 38
21 Boundary conditions Generalized Richards & Hoxey Outlet Inlet 21 / 38
22 Boundary conditions Generalized Richards & Hoxey Outlet Inlet 22 / 38
23 Boundary conditions Generalized Richards & Hoxey τ = ρu 2,sl k z = 0 ε z = u 3,sl κz 2 ν t z = κu,sl (1 + l κz max ) 2 dθ dz = q w σ θ c p ν t Outlet Zero gradients Inlet Equilibrium proles τ = ρu 2,w k z = 0 Wall functions dθ dz = q w c p σ θ ν t 23 / 38
24 Boundary conditions Wall functions Local friction velocity Turbulence properties ε P = Π k,p = u3,w ε P = u,w = Cµ ku,w l ε z = u 3,w 2zP z 0 [ z 0 κu P ln(z P /z0) + κ(z P z0)/l max ( 1 1 ) dz l m κl ( 2zp z0 1 2κ(z P z0) ln Cµ k u 2 Π k,p,w z0 ) + 1 l max 1 κl ] 24 / 38
25 OpenFOAM implementation Boundary conditions Version 1.6 provided derivedfvpatchfields for turbulence models 1 epsilonwallfunction 2 nutroughwallfunction Version 1.7 provided fixedshearstress and inow for U and ε 1 atmboundarylayerinletvelocity 2 atmboundarylayerinletepsilon What is needed: 1 Modied ε and G k wall functions 2 A wall function for proper shear stress for ABL 3 A way to cleanly handle a z0 distribution 4 Ecient IO between libraries 25 / 38
26 Regarding consistency Handling distributed roughness Two options: 1 Treat it like a wall function parameter 2 Treat it as a separate eld Constructor... zzerofield_ = new volscalarfield ( IOobject ( "z0", db().time().timename(), patch().boundarymesh().mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), patch().boundarymesh().mesh() ); z0_ = (*zzerofield_).boundaryfield()[patch().index()]; / 38
27 Regarding consistency Model parameters Boundary conditions often require same information 1 Directives (#) and macro substitutions ($) 2 IOdictionary (RASProperties and surfacelayerproperties) Bundle boundary conditions together with explicit links to dictionaries Constructor... : epsilonwallfunctionfvpatchscalarfield(p, if), z0_(p.size(), 0.0), kappa_(0.40), beta_(0), L_(1e100), lmax_(1e100) { readsurfacelayerdict(); zzerofield_ = new volscalarfield / 38
28 Regarding consistency Model parameters ε gradient Shear stress ν t wall function ν t gradient Θ gradient surfacelayer- Properties RASProperties ε inlet U inlet ε wall function 28 / 38
29 OpenFOAM implementation buoyantsimplefoam1d 1 Drop pressure equation simplefoam... //p.storepreviter(); // Pressure-velocity SIMPLE corrector { #include "UEqn.H" //#include "peqn.h" #include "thetaeqn.h" } // Update boundary conditions U.correctBoundaryConditions(); theta.correctboundaryconditions(); turbulence->correct(); / 38
30 OpenFOAM implementation buoyantsimplefoam1d 2 Add potential temperature equation ( ) d νt dθ = 0 dz σ θ dz thetaeqn.h fvscalarmatrix thetaeqn ( - (1.0/sigmaTheta) * fvm::laplacian(turbulence->nut(), theta) ); thetaeqn.relax(); 30 / 38
31 1D simulations Stable conditions: QW = 20 W/m 2, u,sl = 0.28 m/s, z 0 = 0.3 m, l max = 8 m 10 1 z/lmax 10 0 Proposed model Apsley and Castro Analytical U/u Figure: Comparison of proposed model, Apsley & Castro and analytical solutions 31 / 38
32 1D simulations Stable conditions: QW = 20 W/m 2, u,sl = 0.28 m/s, z 0 = 0.3 m, l max = 8 m z/lmax Proposed - G k Proposed - ε Proposed - G b Analytical - G k Analytical - ε Analytical - G b Normalized TKE budget Figure: Comparison of proposed model, Apsley & Castro and analytical solutions 32 / 38
33 1D simulations Neutral length-limited conditions: QW = 0 W/m 2, u,sl = 0.65 m/s, z 0 = 0.3 m, l max = 36 m 10 1 z/lmax 10 0 Proposed model Apsley and Castro Analytical U/u Figure: Comparison of proposed model, Apsley & Castro and analytical solutions 33 / 38
34 The Cε3 coecient Calibration The proposed ε equation can be recast in the form proposed by Rodi for modeling stratied ows: This yields Cε3 = (C ε2 Cε1) Cε1 [ Pε = Cε1(1 + Cε3R f )Π kε 1 (1 + R f )2 ( lm l max + 1 k )( 1 l m l max ) 3 ] 1 R f 34 / 38
35 The Cε3 coecient Comparison 10 Cε3Cε1/(Cε2 Cε1) Betts & Haroutunian (1983) Rodi (1987) Kitada (1987) Burchard & Baumert (1995) Apsley & Castro (1997) Freedman & Jackson (2003) Alinot & Masson (2005) Sumner & Masson (2011) z/l Figure: Comparison of proposals for C ε3 35 / 38
36 Conclusions Summary A slightly modied k ε model for neutral and stable, length-limited surface layer ow A corresponding suite of boundary conditions A new Cε3 expression has been derived that ensures similarity theory for a stable atmosphere is a solution of the model equations 36 / 38
37 Conclusions Future work Extend to two and three-dimensional ow cases Prepare tutorial and libraries for release to extend project 37 / 38
38 Thank you Questions? 38 / 38
On the validation study devoted to stratified atmospheric flow over an isolated hill
On the validation study devoted to stratified atmospheric flow over an isolated hill Sládek I. 2/, Kozel K. 1/, Jaňour Z. 2/ 1/ U1211, Faculty of Mechanical Engineering, Czech Technical University in Prague.
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 informationτ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization
τ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization Note that w *3 /z i is used to normalized the TKE equation in case of free
More information1 The Richardson Number 1 1a Flux Richardson Number b Gradient Richardson Number c Bulk Richardson Number The Obukhov Length 3
Contents 1 The Richardson Number 1 1a Flux Richardson Number...................... 1 1b Gradient Richardson Number.................... 2 1c Bulk Richardson Number...................... 3 2 The Obukhov
More informationThe implementation of Two-equation SGS turbulence model and Wall damping function
The implementation of Two-equation SGS turbulence model and Wall damping function Yeru Shang Thermo-fluid Mechanics Research Centre, The University of Sussex, Brighton, UK 2017-11-21 Yeru Shang The implementation
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
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 informationPRELIMINARY STUDY OF COMPUTATIONAL SETUP FOR URBAN STREET CANYONS. by MUHAMMAD NOOR AFIQ WITRI, M.Eng
PRELIMINARY STUDY OF COMPUTATIONAL SETUP FOR URBAN STREET CANYONS by MUHAMMAD NOOR AFIQ WITRI, M.Eng 1 CONTENTS 1.Introduction 2.Building Configuration 3.Boundary Condition 4.Previous Works 5.Summary 2
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 informationPANS turbulence model implementation
CFD with OpenSource Software A course at Chalmers University of Technology Taught by Håkan Nillson PANS turbulence model implementation Developed for OpenFOAM-2.3.x Author: Guglielmo Minelli Peer reviewed
More informationSpring Semester 2011 March 1, 2011
METR 130: Lecture 3 - Atmospheric Surface Layer (SL - Neutral Stratification (Log-law wind profile - Stable/Unstable Stratification (Monin-Obukhov Similarity Theory Spring Semester 011 March 1, 011 Reading
More information(Wind profile) Chapter five. 5.1 The Nature of Airflow over the surface:
Chapter five (Wind profile) 5.1 The Nature of Airflow over the surface: The fluid moving over a level surface exerts a horizontal force on the surface in the direction of motion of the fluid, such a drag
More informationviscousheatingsolver Martin Becker
viscousheatingsolver Adding temperature transport and viscous heating to simplefoam Solver and Test Case Martin Becker martin_becker@dhcae-tools.de 11/09/2012 Abstract The modifications of the simplefoam
More informationStatistical 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 informationA new k-epsilon model consistent with Monin-Obukhov similarity theory
Downloaded from orbit.dtu.dk on: Feb 1, 018 A new k-epsilon model consistent with Monin-Obukhov similarity theory van der Laan, Paul; Kelly, Mark C.; Sørensen, Niels N. Published in: Wind Energy Link to
More informationNumerical 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 informationCFD with OpenSource software A course at Chalmers University of Technology Taught by HÅKAN NILSSON. Project work: intersettlingfoam
Taught by HÅKAN NILSSON Project work: intersettlingfoa Developed for OpenFOAM-2.2.0 Author: Pedra Rain Deceber, 2013 2/21/2014 1 Presentation Outline: The objective of this presentation Introduction to
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 informationSergej S. Zilitinkevich 1,2,3. Division of Atmospheric Sciences, University of Helsinki, Finland
Atmospheric Planetary Boundary Layers (ABLs / PBLs) in stable, neural and unstable stratification: scaling, data, analytical models and surface-flux algorithms Sergej S. Zilitinkevich 1,,3 1 Division of
More informationBoundary 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 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 informationA Thorough Description Of How Wall Functions Are Implemented In OpenFOAM
Cite as: Fangqing Liu.: A Thorough Description Of How Wall Functions Are Implemented In OpenFOAM. In Proceedings of CFD with OpenSource Software, 2016, Edited by Nilsson. H., http://www.tfd.chalmers.se/~hani/kurser/os_cfd_2016
More informationA new source term in the parameterized TKE equation being of relevance for the stable boundary layer - The circulation term
A new source term in the parameteried TKE equation being of releance for the stable boundary layer - The circulation term Matthias Raschendorfer DWD Ref.: Principals of a moist non local roughness layer
More informationCalculations on a heated cylinder case
Calculations on a heated cylinder case J. C. Uribe and D. Laurence 1 Introduction In order to evaluate the wall functions in version 1.3 of Code Saturne, a heated cylinder case has been chosen. The case
More informationHow to implement your own turbulence model (1/3)
How to implement your own turbulence model 1/3) The implementations of the turbulence models are located in $FOAM_SRC/turbulenceModels Copy the source of the turbulence model that is most similar to what
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 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 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 informationInitial and Boundary Conditions
Initial and Boundary Conditions Initial- and boundary conditions are needed For a steady problems correct initial conditions is important to reduce computational time and reach convergence Boundary conditions
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 informationStable Boundary Layer Parameterization
Stable Boundary Layer Parameterization Sergej S. Zilitinkevich Meteorological Research, Finnish Meteorological Institute, Helsinki, Finland Atmospheric Sciences and Geophysics,, Finland Nansen Environmental
More information2.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 informationModelling diabatic atmospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL
Modelling diabatic atospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL Florian Vendel 1, Guillevic Laaison 1, Lionel Soulhac 1, Ludovic Donnat 2, Olivier Duclaux 2,
More informationPart I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows
Physical modeling of atmospheric boundary layer flows Part I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows Outline Evgeni Fedorovich School of
More informationLand/Atmosphere Interface: Importance to Global Change
Land/Atmosphere Interface: Importance to Global Change Chuixiang Yi School of Earth and Environmental Sciences Queens College, City University of New York Outline Land/atmosphere interface Fundamental
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 informationThe Stable Boundary layer
The Stable Boundary layer the statistically stable or stratified regime occurs when surface is cooler than the air The stable BL forms at night over land (Nocturnal Boundary Layer) or when warm air travels
More informationOpenFOAM Course Final Assignment: Tutorial for Natural Convection Boundary Layer
OpenFOAM Course Final Assignment: Tutorial for Natural Convection Boundary Layer Abolfazl Shiri 751028-2051 1 Introduction The aim of this project is to simulate the natural convection boundary layer along
More informationTurbulence 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 informationDirect 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 informationTurbulent transport of the energy in the entrainment interface layer
Turbulent transport of the energy in the entrainment interface layer Marta K. Kopec 1, S. P. Malinowski1, Z. P. Piotrowski2 1 nstitute 2 nstitute of Geophysics, Faculty of Physics, University of Warsaw,
More informationHow to implement your own turbulence model(1/3)
How to implement your own turbulence model(1/3) The implementations of the turbulence models are located in $FOAM_SRC/turbulenceModels Copythesourceoftheturbulencemodelthatismostsimilartowhatyouwanttodo.Inthis
More informationFall Colloquium on the Physics of Weather and Climate: Regional Weather Predictability and Modelling. 29 September - 10 October, 2008
1966-10 Fall Colloquium on the Physics of Weather and Climate: Regional Weather Predictability and Modelling 9 September - 10 October, 008 Physic of stable ABL and PBL? Possible improvements of their parameterizations
More informationMicro-Scale CFD Modeling of Packed-Beds
Micro-Scale CFD Modeling of Packed-Beds Daniel P. Combest and Dr. P.A. Ramachandran and Dr. M.P. Dudukovic Chemical Reaction Engineering Laboratory (CREL) Department of Energy, Environmental, and Chemical
More informationPARAMETERIZATION OF TURBULENCE DAMPING IN SEDIMENT-LADEN FLOW
Civil Engineering Department Katholieke Universiteit Leuven PARAMETERIZATION OF TURBULENCE DAMPING IN SEDIMENT-LADEN FLOW Report No. HYD/ET/00/COSINUS/3 December 2000 Erik A. TOORMAN A contribution to
More informationAdding electric conduction and Joule heating to chtmultiregionfoam Niklas Järvstråt
CFD with open source software, assignment 3: Adding electric conduction and Joule heating to chtmultiregionfoam Niklas Järvstråt Reviewed by Arash Eslamdoost 1 Table of contents 1.Physical background and
More informationReynolds Averaging. We separate the dynamical fields into slowly varying mean fields and rapidly varying turbulent components.
Reynolds Averaging Reynolds Averaging We separate the dynamical fields into sloly varying mean fields and rapidly varying turbulent components. Reynolds Averaging We separate the dynamical fields into
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 informationAtmospheric Boundary Layer Studies with Unified RANS-LES and Dynamic LES Methods
5st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 7 - January 23, Grapevine (Dallas/Ft. Worth Region), Texas AIAA 23-747 Atmospheric Boundary Layer Studies with
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 informationBoundary layer processes. Bjorn Stevens Max Planck Institute for Meteorology, Hamburg
Boundary layer processes Bjorn Stevens Max Planck Institute for Meteorology, Hamburg The Atmospheric Boundary Layer (ABL) An Abstraction (Wippermann 76) The bottom 100-3000 m of the Troposphere (Stull
More informationCHAPTER 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 informationA Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM
206 9th International Conference on Developments in esystems Engineering A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM Hayder Al-Jelawy, Stefan Kaczmarczyk
More informationEffect of modification to tongue and basic circle diameter on vibration in a double-suction centrifugal pump
5th International Conference on Information Engineering for Mechanics and Materials (ICIMM 2015) Effect of modification to tongue and basic circle diameter on vibration in a double-suction centrifugal
More informationLaminar Flow. Chapter ZERO PRESSURE GRADIENT
Chapter 2 Laminar Flow 2.1 ZERO PRESSRE GRADIENT Problem 2.1.1 Consider a uniform flow of velocity over a flat plate of length L of a fluid of kinematic viscosity ν. Assume that the fluid is incompressible
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 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 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 informationADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS
ADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS Radi Sadek 1, Lionel Soulhac 1, Fabien Brocheton 2 and Emmanuel Buisson 2 1 Laboratoire de Mécanique des Fluides et d Acoustique,
More informationModelization of saturated sand flux
Modelization of saturated sand flux O. Durán 1 and H. Herrmann 1, 1 Institute for Computer Physics, University of Stuttgart, 7569 Stuttgart, Germany. Departamento de Física, Universidade Federal do Ceará,
More informationAtmospheric Boundary Layers
Lecture for International Summer School on the Atmospheric Boundary Layer, Les Houches, France, June 17, 2008 Atmospheric Boundary Layers Bert Holtslag Introducing the latest developments in theoretical
More informationRhie-Chow interpolation in OpenFOAM 1
Rhie-Chow interpolation in OpenFOAM 1 Fabian Peng Kärrholm Department of Applied Mechanics Chalmers Univesity of Technology Göteborg, Sweden, 2006 1 Appendix from Numerical Modelling of Diesel Spray Injection
More informationComparison of Numerical Prediction of Pressure Coefficient on Rectangular Small Building
May 26, Volume 3, Issue 5 Comparison of Numerical Prediction of Pressure Coefficient on Rectangular Small Building Neel M. Patel, 2 Satyen D. Ramani PG Student, Department of Civil Engineering, 2 Asst.
More informationBOUNDARY-LAYER METEOROLOGY
BOUNDARY-LAYER METEOROLOGY Han van Dop, September 2008 D 08-01 Contents 1 INTRODUCTION 3 1.1 Atmospheric thermodynamics.................... 3 1.2 Statistical aspects of fluid mechanics................
More informationMODELING THE NEUTRAL ATMOSPHERIC BOUNDARY LAYER BASED ON THE STANDARD k-ε TURBULENT MODEL: MODIFIED WALL FUNCTION
The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 29, Taipei, Taiwan MODELING THE NEUTRAL ATMOSPHERIC BOUNDARY LAYER BASED ON THE STANDARD k-ε TURBULENT MODEL: MODIFIED WALL FUNCTION
More informationPart III: Modeling atmospheric convective boundary layer (CBL) Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, USA
Physical modeling of atmospheric boundary layer flows Part III: Modeling atmospheric convective boundary layer (CBL) Outline Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, USA
More informationAtmospheric stability parameters and sea storm severity
Coastal Engineering 81 Atmospheric stability parameters and sea storm severity G. Benassai & L. Zuzolo Institute of Meteorology & Oceanography, Parthenope University, Naples, Italy Abstract The preliminary
More informationAccommodating LES to high Re numbers: RANS-based, or a new strategy?
Symposium on Methods 14-15 July 2005, FOI, Stockholm, Sweden, Accommodating LES to high Re numbers: RANS-based, or a new strategy? K. Hanjalić Delft University of Technology, The Netherlands Guest Professor,
More informationIntroduction. Haukur Elvar Hafsteinsson. OpenFOAM - Porous Media 1. The Governing Equations. The Class. The solver. Porous media as a flow control
Introduction The Governing Equations The Class The solver Porous media as a flow control Porous media in cylindrical coordinates Available at the course home page: http://www.tfd.chalmers.se/ hani/kurser/os_cfd_2008/
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 informationTutorial for the supercritical pressure pipe with STAR-CCM+
Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
More informationConsistent Two-Equation Closure Modelling for Atmospheric Research: Buoyancy and Vegetation Implementations
Boundary-Layer Meteorol (2012) 145:307 327 DOI 10.1007/s10546-012-9726-5 ARTICLE Consistent Two-Equation Closure Modelling for Atmospheric Research: Buoyancy and Vegetation Implementations Andrey Sogachev
More informationImplementation of the Quasi-Normal Scale Elimination (QNSE) Model of Stably Stratified Turbulence in WRF
Implementation of the Quasi-ormal Scale Elimination (QSE) odel of Stably Stratified Turbulence in WRF Semion Sukoriansky (Ben-Gurion University of the egev Beer-Sheva, Israel) Implementation of the Quasi-ormal
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 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 informationTutorial for the heated pipe with constant fluid properties in STAR-CCM+
Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities
More informationOutline: 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 informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More information2σ e s (r,t) = e s (T)exp( rr v ρ l T ) = exp( ) 2σ R v ρ l Tln(e/e s (T)) e s (f H2 O,r,T) = f H2 O
Formulas/Constants, Physics/Oceanography 4510/5510 B Atmospheric Physics II N A = 6.02 10 23 molecules/mole (Avogadro s number) 1 mb = 100 Pa 1 Pa = 1 N/m 2 Γ d = 9.8 o C/km (dry adiabatic lapse rate)
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 informationCFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer
Accepted for publication in International Journal of Heat and Mass Transfer CFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer Thijs Defraeye a, *,
More informationCOMMENTS ON "FLUX-GRADIENT RELATIONSHIP, SELF-CORRELATION AND INTERMITTENCY IN THE STABLE BOUNDARY LAYER" Zbigniew Sorbjan
COMMENTS ON "FLUX-GRADIENT RELATIONSHIP, SELF-CORRELATION AND INTERMITTENCY IN THE STABLE BOUNDARY LAYER" Zbigniew Sorbjan Department of Physics, Marquette University, Milwaukee, WI 5301, U.S.A. A comment
More informationDEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M.
DEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M. Sahu 1, Kishanjit Kumar Khatua and Kanhu Charan Patra 3, T. Naik 4 1, &3 Department of Civil Engineering, National Institute of technology,
More informationModeling Separation and Reattachment Using the Turbulent Potential Model
Modeling Separation and Reattachment Using the urbulent Potential Model Blair Perot & Hudong Wang Department of Mechanical Engineering University of Massachusetts, Amherst, MA ABSRAC A new type of turbulence
More informationTutorial XiFoam. Applied Mechanics Department Chalmers University of Technology. CFD with OpenSource software, assignment 3
Applied Mechanics Department Chalmers University of Technology CFD with OpenSource software, assignment 3 Tutorial XiFoam Version: OpenFOAM-1.7.x Author: Ehsan Yasari Peer reviewed by: Josef Runsten Jelena
More informationMATHEMATICAL MODELING AND NUMERICAL SOLUTION OF 3D ATMOSPHERIC BOUNDARY LAYER
, Vol, Pt, Special Issue Proceedings of International Conference RDAMM 585 MATHEMATICAL MODELING AND NUMERICAL SOLUTION OF D ATMOSPHERIC BOUNDARY LAYER L. Beneš, K. Kozel Department of Technical Mathematics,
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 informationImplementation of Turbulent Viscosity from EARSM for Two Equation Turbulence Model
CFD with OpenSource software A course at Chalmers University of Technology Taught by Håkan Nilsson Project work: Implementation of Turbulent Viscosity from EARSM for Two Equation Turbulence Model Developed
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 informationTurbulent Convection in Air
SMR.1771-33 Conference and Euromech Colloquium #480 on High Rayleigh Number Convection 4-8 Sept., 2006, ICTP, Trieste, Italy ------------------------------------------------------------------------------------------------------------------------
More informationHigh Reynolds Number Wall Turbulence: Facilities, Measurement Techniques and Challenges
High Reynolds Number Wall Turbulence: Facilities, Measurement Techniques and Challenges Beverley McKeon Graduate Aeronautical Laboratories California Institute of Technology http://mckeon.caltech.edu Chapter
More informationESCI 485 Air/sea Interaction Lesson 3 The Surface Layer
ESCI 485 Air/sea Interaction Lesson 3 he Surface Layer References: Air-sea Interaction: Laws and Mechanisms, Csanady Structure of the Atmospheric Boundary Layer, Sorbjan HE PLANEARY BOUNDARY LAYER he atmospheric
More informationLECTURE 2: FLUID MECHANICS
LECTURE 2: FLUID MECHANICS Introduction Conservation of mass and momentum General types of flow Laminar vs. turbulent flow Shear Stress Reach-average shear stress Bed roughness and reach average flow velocity
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 informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationModelling 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 informationExternal Flows. Dye streak. turbulent. laminar transition
Eternal Flos An internal flo is surrounded by solid boundaries that can restrict the development of its boundary layer, for eample, a pipe flo. An eternal flo, on the other hand, are flos over bodies immersed
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 informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More information