Mostafa Momen. Project Report Numerical Investigation of Turbulence Models. 2.29: Numerical Fluid Mechanics
|
|
- Merry Charles
- 6 years ago
- Views:
Transcription
1 2.29: Numerical Fluid Mechanics Project Report Numerical Investigation of Turbulence Models Mostafa Momen May 2015 Massachusetts Institute of Technology 1
2 Numerical Investigation of Turbulence Models Term Project: Mostafa Momen 1. Introduction: Most real-world flows occurring in nature and engineering applications are turbulent. The boundary layer in the earth s atmosphere, interstellar gas clouds, boundary layers on aircraft wings, natural gas flow in pipes, and water flow in rivers are all examples of turbulent flows. Turbulence happens when the Reynolds number of a flow is high, typically once Re > O(1000). Reynolds number is the ratio of the inertial to viscous forces, hence when the inertial forces are dominating; the flow regime will become turbulent. Furthermore, when the small and large scales of the motion or dissipating and energy containing scales in a flow are well separated from each other the turbulence occurs. To simulate these highly nonlinear turbulent flows numerically, there are three primary methods. Among three available approaches Direct Numerical Simulation (DNS) is the most accurate one but also most computationally expensive (sometime impossible) way in high Reynolds flows. Reynolds-Averaged Navier-Stokes (RANS) models, which only solve for the mean flow and model all turbulent scales are the least accurate and the cheapest methods at the same time which are good for some engineering applications. The DNS grids are very fine and in the size of the Kolmogorov length scale so this method requires enormous computational resources for highly turbulent flows. Thus in most engineering applications, RANS models are used. Large-eddy simulation (LES) technique, which consists of resolving governing equations for the energy containing large scales and modeling the smaller scales, lies between these two approaches. In the following study, we will examine a finite-volume (FV) DNS code with some RANS turbulence models versus finite-difference LES results and analyze their differences. To do so, we will employ the methods and notions learned in class to develop the numerical anlaysis. 2. RANS Closure Models: In the RANS approaches to turbulence, all of the unsteadiness is averaged out implying that all unsteadiness is regarded as part of the turbulence. The RANS and energy equations cannot be solved without information about the various correlation terms that make up the stress tensor. Thus, to close these equations one must model all the turbulence part and Reynolds stresses. It is well known that these terms, which represent turbulent diffusion, are much larger than those corresponding to laminar diffusion except in the immediate vicinity of a wall, and in turbulent wall boundary layers, wakes, jets and more complex flows, these turbulent diffusion terms are of similar magnitude to the convective terms. Early approaches to turbulence modeling include the mixing length assumptions of Prandtl (1945) and eddy-viscosity assumptions of Boussinesq (1877) for wall boundary 2
3 layers and jets. Kolmogorov (1942) and Rotta (1962) proposed models based on partialdifferential equations but, in the absence of digital computers, could not solve them (Cebeci 2005). In general, the RANS based turbulence models are divided into zero, one and two equation models according to the number of the equations that they solve to obtain the turbulent viscosity. The zero-equation models also called algebraic models are calculated directly from the flow variables. An alternative approach for preliminary turbulence modeling is to construct (differential) transport equations for some of the turbulence quantities and to model higher-order terms, which turn out to be triple correlations (Fletcher 1991). There are several two-equation models, which are based on solving differential equations. Three of the most popular and widely used models are the k-ε model of Jones and Launder (1972), the k-ω model of Wilcox (1998) and the SST model of Menter (1994) which blends the k-ε model in the outer region and k-ω model in the near wall region. We may write the main RANS based turbulence models according to their structure as in the following: a) Zero-equation models (Algebraic models) 1- Cebeci-Smith model 2- Baldwin-Lomax model 3- Johnson-King moel 4- A roughness dependent model 5- Smagorinsky model b) One equation models 1- Prandtl s one-equation model 2- Baldwin-Barth model 3- Spalart-Allmaras model 4- Rahman-Siikonen-Agarwal model c) Two equation models 1- k-ε models 2- k-ω models 3- Realisability issues Here we focus on two approaches and show the results of the simulations for these two RANS models: 2.1 Smagorinsky model This is a simple model that computes the turbulent viscosity based on the local derivatives of the velocity field and the local grid size as the following: (1) This formula was first proposed by Smagorinsky (1964) for the subgrid scale eddy viscosity and is sensitive to the grid sizes. 3
4 2.2 k-ε model The k-ε model is one of the most popular turbulence models and is the most common model used in CFD to simulate mean flow characteristics in turbulent flows, even though it does not perform well in some cases such as when pressure gradients are very large. The Reynolds stresses might be written as: u i u j = ν T ui x j + u j x i where k is the turbulent kinetic energy (TKE) defined as : 2 3 δ k (2) ij k = 1 2 u u i i (3) The exact k-ε equations include many unknown terms. For a much more practical approach, the standard k-ε model that is proposed by Launder and Spalding (1974) is used to minimize the unknowns of the problem and present a set of equations which can be applied to a large number of turbulent applications. Thus, this model contains two model equations for the TKE and the dissipation as: Model equation for k : ( ) k t + ku i = ν + ν T k x i x j σ k + P x k ε j (4) The second model equation may be expressed as below: Model equation for ε : ( ) ε t + εu i = ν + ν T ε x i x j σ ε x j ε + C ε1 k P C ε k ε 2 k ε (5) where P k =2ν T.EijEij is the rate of production of TKE and the final eddy viscosity would be: k 2 ν T = C µ ε (6) These two equations (4) and (5) are nonlinearly coupled equations that we need to solve together. 4
5 3. Numerical Simulations: Three numerical simulations have been performed to investigate the turbulence models. First and third cases are implemented using 2.29 FV code in different boundaries and second case is done with the help of a 3D LES code on a homogenous boundary. a) Case 1: 2D FV on a homogenous terrain By employing the 2.29 FV code a wall-bounded flow is simulated with the following numerical details and the numerical domain could be viewed in Fig. 1: - Grid size: Boundary Conditions: No-slip condition on the ground, Time-dependent BC s elsewhere with Periodic BC in horizontal directions - Size of domain: Y=100m, X=2π 100m - Time of simulation: 5s with Dt =0.01s - ν =1.695E-3 (almost 0.01 of the air and real-world atmospheric boundary layer) - Wind profiles: perturbed logarithmic profiles Fig.1. Numerical domain of case 1 showing velocity profiles A time-dependent boundary condition is used and a random noise with about 10% of the mean value is added to the wind profiles to make the flow regime turbulent. The formation of turbulent eddies are clear in Fig. 1 near the ground. If one increases the perturbations, eddies and thus the turbulence would grow accordingly. All the involving gradients in equations (4) and (5) are calculated in the main loop of fluid solver. And then two k-ε model equations are treated fully explicitly at each timestep and solved simultaneously to obtain the following results. Numerical Stability: Note that with these conditions sometime the explicit scheme becomes unstable and the solutions grow in time (especially when the generated noises are high). This instability spreads to the whole domain through dispersion terms and blows up all the solutions. Hence if the scheme yields unstable solutions at just one single point, it will blow up the whole domain in the next steps due to spatial dependences. One way to tackle this problem is to use implicit schemes considering that the equations are coupled and nonlinear which might make it more challenging. 5
6 Fig.2. TKE results obtained from the FV DNS code (left) and k-ε equation (right) in case 1 Figure 2 depicts two TKE plots resulted directly from the FV code (left panel) and also the TKE computed from the k-ε equations (right panel). As it is clear from above, they are not exactly the same. They should not be necessarily identical since the k-ε equations are not exact and this approach employs a simplistic assumption about dissipation equation. However, they are both in the same order of magnitude and are somehow correlated. The k-ε model is more sensitive to the time-varying boundary conditions and as one could observe there are some local maxima around the left BC. To understand why these two are different, the production terms of the DNS and k-ε equations (P k ) which are the largest terms in the TKE budget equations are computed and plotted in Fig. 3. Fig.3. TKE production term for FV DNS code (left) and k-ε model, P K (right) As one could observe, the TKE productions are different in two cases, indicating why the TKE in Fig.2 must not be the same. That is because the real production term is different with the production term resulted from equation (5). Now, turbulent viscosity resulted from the k-ε model in equation (6) is plotted against the Smagorinsky model in Fig. 4 to demonstrate their differences. This figure shows that these tow models are yielding different eddy viscosity values due to their distinct structure, however they both peak near the boundary where the turbulence is high. This is true, but as one could see the Smagorinsky model is dependent on the local velocity derivatives which are highest near the wall and is sensitive to the grid sizes, while in k-ε, eddy viscosity depends on the TKE and dissipation values of equations (4) and (5). Thus, k-ε result somehow looks similar to the TKE depicted in the right panel of Fig.2. 6
7 Fig.4. Turbulent viscosity from Smagorinsky (left) and k-ε models (right) b) Case 2: 3D LES on a homogenous terrain Using a large-eddy simulation code, a wall-bounded flow similar to case 1 is simulated with the following numerical details: - Grid size: Boundary Conditions: No-slip condition on the ground, Periodic in horizontal directions, Stress free top BC - Size of domain: Y=100m, X=2π 100m, Z=2π 100m - Time of simulation: 2 hour warm-up + 20 hour to reach to quasi-steady state condition with Dt =0.1s - Wind profiles: 3D Ekman profiles - Type of simulation: Coriolis forcing with geostrophic wind velocity of about 8m/s - Finite difference methods: fully-explicit 2 nd order Adams-Bashforth for timeadvancement, and pseudo spectral approach in horizontal direction with 2 nd order centered-difference in vertical direction - Subgrid scale model used: Lagrangian scale-dependent dynamic model - Courant-Friedrichs-Levy condition for stability ~ 0.3 To compare the previous results with a LES code, case 2 is implemented as above. Figure 5 displays the LES results after 22-hour simulation in quasi-equilibrium state. The left panel of this figure shows the TKE results which are analogous to Fig.2 in terms of the order of the magnitude of the results and also the fact that both peak near the ground that confirm previous results. Fig.5. LES results: TKE (left) and Reynolds stresses (right) 7
8 Furthermore, this figure also indicates the plot of the dominating Reynolds stresses, which again peak near the wall. Note that all results in this figure are averaged in the z- direction (other horizontal direction) and that is why they look smoother than the previous results. Figure 6 demonstrates the output eddy-viscosity of LES code, which is determined from the following equation: ν T = du dy R R dw dy 2 (7) Fig.6. Turbulent viscosity of LES code As one could observe from above, the turbulent viscosity peaks almost in 50m above the ground, which is different from Fig. 4 where it is highest about 5m above the boundary. This difference could be due to the inherent differences between the LES, DNS and RANS models used here as well as the distinct numerical setup of the problem in case 1 and 2. The sum of Reynolds stresses in LES is shown in Fig.5 that peaks near the ground as expected. However, the velocity gradients or strain terms in equation (7) are also high near the wall and that is why in Fig.6 eddy viscosity does not peak near the ground. Nevertheless, the value of the eddy-viscosity in both cases is in the same order of magnitude and that qualitatively confirms the previous results. c) Case 3: FV on inhomogeneous terrains One of the recent research problems that could be investigated further using these codes is studying the mean and turbulence characteristics in inhomogeneous surfaces, which is the case in most of the real-world situations. Investigating the wind effects and pollution dispersion on urban geometries and canyons is one of the current research efforts. Moreover, understanding the physical mechanisms behind the katabatic and anabatic winds on the steep slopes or mountains which are usually driven by heat and occur every day in sloping areas is vey important. Figure 7 demonstrates two examples of such flows. 8
9 Fig.7. Flow field around a building (left) and mountaintop (right) courtesy of U of Tokyo and UAF Hence to that end, we have used the FV code in case 1 and added a mask layer with an elliptical equation near the ground in the middle to simulate the wind effects on a hill. Note that this process is implemented conveniently in FV code; However its implementation is a little challenging in finite difference codes and may require additional schemes such as immersed boundary methods. We have defined the following numerical domain in Fig. 8. All other numerical details are the same as case 1 except the additional mask layer in the middle. Fig.8. Numerical domain of case 3 showing velocity profiles First results are the TKE plots depicted in Fig. 9. Similar to case 1, both of the plots are analogues with little differences due to the structure of the models. We also observe that both graphs peak near the boundaries and also over the hill. The second result compares the Smagorisnky model with k-ε eddy-viscoity in such slope flows. Figure 10 displays that they are not the same again, but both are high near the wall. The Smagorinsky model even shows higher peak on top of the inhomogeneity due to high velocity gradients. 9
10 Fig.9. TKE results obtained from the FV DNS code (left) and k-ε equation (right) in case 3 Fig.10. Turbulent viscosity from Smagorinsky (left) and k-ε models (right) in case 3 4. Conclusion: The 2.29 FV code was used to produce turbulent flows. To that end logarithmic wind profiles were imposed and a random noise was added to their values at each time-step. Two turbulence models, including Smagorinsky and k-ε, were implemented and tested in two cases using this FV code. The results were satisfactory showing peaks for the TKE and production terms near the ground. k-ε model is relatively easy to implement in a fully explicit way, however it becomes challenging in implicit schemes since it includes two nonlinearly coupled equations. On the other hand, k-ε model uses simplistic dissipation equation and works poorly for flows with strong separations and rotating flows. Case 2 employed a 3D large-eddy simulation code to demonstrate the differences among the 2D FV and 3D LES results. The TKE graphs looked similar and the eddy-viscosity magnitudes were in the same order though there were some differences in their profiles. The last case was aimed to show further practical applications, which are the focus of some of the recent research efforts. The mean and turbulence dynamics on inhomogeneous terrains could be investigated using these numerical codes. In fact in highly turbulent flows, which most of times is the case particularly in geophysical flows, the implementation of DNS is impossible considering the current computational resources. Thus, the provided idea here in using the turbulence models to convert DNS to RANS code will be very helpful for further research developments. 10
11 5. Future work: In order to change the DNS code into RANS, one needs to convert the scalar viscosity defined in the code into a space-dependent turbulent viscosity resulted from the k-ε model in an efficient way. This could be implemented using appropriate optimized algorithms and then the 2.29 FV code could be used in higher Reynolds number and lower viscosity flows. We might also improve the generated turbulence mechanism by using temporally and spatially correlated noises. This will enhance the quality of the generated perturbations and they should be verified with some other results too. After implementing these two steps, the created RANS code could be used in many research disciplines that include high Reynolds flows such as wind over inhomogeneous terrains such as the one investigated in case 3 of this project. References: Boussinesq, J.: Theorie de l'ecoulement tourbillant, Mem. pres. Acad. Sc. XXIII, 46, Cebeci, T., J. P. Shao, F. Kafyeke, and E. Laurendeau, 2005: Computational Fluid Dynamics for Engineers. 328ff. pp. Ferziger, J. H., and I. M. Peric: Computational Methods for Fluid Dynamics. 3rd ed, Fletcher, C. A. J.: Computational Techniques for Fluid Dynamics. 2nd ed. 439 pp, Jones, W. P. and Launders, B.E.: The Predicition of Laminarization with a Two-Equation Model of Turbulence. Int. J. Heat and Mass Transfer 15, , Kolmogorov, A.N.: Equations of Turbulent Motion of an Incompressible Fluid, Izv. Akad. Nauk. SSR Ser. Phys. 6, 56, (English translation, Imperial College, Mech. Eng. Dept. Report ON/6, 1968). Menter, F.R.: Two-Equation Eddy Viscosity Turbulence Models for Engineering Applications. AIAA J. 32, , Prandtl, L.: Uber ein neues Formelsystem fur die ausgebildete Turbulenz, Nachrichten von der Akad. der Wissenschaft in Gottingen, Rotta, J.C.: Turbulent Boundary Layers in Incompressible Flow, Prog. Aero. Sci. 2, 1, Wilcox, D. C: Turbulence Modeling for CFD. DCW Industries, Inc., 5354 Palm Drive, La Canada, Calif.,
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 informationAER1310: 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 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 informationTurbulence 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 informationTurbulence: 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 informationThere are no simple turbulent flows
Turbulence 1 There are no simple turbulent flows Turbulent boundary layer: Instantaneous velocity field (snapshot) Ref: Prof. M. Gad-el-Hak, University of Notre Dame Prediction of turbulent flows standard
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 informationComputational 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 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 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 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 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 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 informationProbability 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 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 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 informationWhat 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 informationAn 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 informationPhysics of turbulent flow
ECL-MOD 3A & MSc. Physics of turbulent flow Christophe Bailly Université de Lyon, Ecole Centrale de Lyon & LMFA - UMR CNRS 5509 http://acoustique.ec-lyon.fr Outline of the course A short introduction to
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 informationROBUST EDDY VISCOSITY TURBULENCE MODELING WITH ELLIPTIC RELAXATION FOR EXTERNAL BUILDING FLOW ANALYSIS
August 11 13, ROBUST EDDY VISCOSITY TURBULENCE MODELING WITH ELLIPTIC RELAXATION FOR EXTERNAL BUILDING FLOW ANALYSIS Mirza Popovac AIT Austrian Institute of Technology Österreichisches Forschungs- und
More informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
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 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 informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationAn 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 informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationOptimizing calculation costs of tubulent flows with RANS/LES methods
Optimizing calculation costs of tubulent flows with RANS/LES methods Investigation in separated flows C. Friess, R. Manceau Dpt. Fluid Flow, Heat Transfer, Combustion Institute PPrime, CNRS University
More informationOn 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 informationOn 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 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 informationLES of wind turbulence and heat environment around dense tall buildings
EACWE 5 Florence, Italy 19 th 23 rd July 2009 LES of wind turbulence and heat environment around dense tall buildings Flying Sphere image Museo Ideale L. Da Vinci Tsuyoshi Nozu 1, Takeshi Kishida 2, Tetsuro
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 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 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 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 informationLARGE 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 informationTurbulent eddies in the RANS/LES transition region
Turbulent eddies in the RANS/LES transition region Ugo Piomelli Senthil Radhakrishnan Giuseppe De Prisco University of Maryland College Park, MD, USA Research sponsored by the ONR and AFOSR Outline Motivation
More informationEulerian models. 2.1 Basic equations
2 Eulerian models In this chapter we give a short overview of the Eulerian techniques for modelling turbulent flows, transport and chemical reactions. We first present the basic Eulerian equations describing
More informationHybrid 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 informationIntroduction to ANSYS FLUENT
Lecture 6 Turbulence 14. 0 Release Introduction to ANSYS FLUENT 1 2011 ANSYS, Inc. January 19, 2012 Lecture Theme: Introduction The majority of engineering flows are turbulent. Successfully simulating
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 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 informationNumerical Simulation of a Blunt Airfoil Wake
6th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 7 Numerical Simulation of a Blunt Airfoil Wake C.J. Doolan School of Mechanical Engineering University of Adelaide,
More informationWind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting
Wind Flow Modeling The Basis for Resource Assessment and Wind Power Forecasting Detlev Heinemann ForWind Center for Wind Energy Research Energy Meteorology Unit, Oldenburg University Contents Model Physics
More informationSG Turbulence models for CFD
SG2218 2012 Turbulence models for CFD Stefan Wallin Linné FLOW Centre Dept of Mechanics, KTH Dept. of Aeronautics and Systems Integration, FOI There are no simple turbulent flows Turbulent boundary layer:
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
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 informationTurbulence Laboratory
Objective: CE 319F Elementary Mechanics of Fluids Department of Civil, Architectural and Environmental Engineering The University of Texas at Austin Turbulence Laboratory The objective of this laboratory
More informationModeling 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 informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,800 116,000 120M Open access books available International authors and editors Downloads Our
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 informationIntroduction 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 informationENGINEERING MECHANICS 2012 pp Svratka, Czech Republic, May 14 17, 2012 Paper #195
. 18 m 2012 th International Conference ENGINEERING MECHANICS 2012 pp. 309 315 Svratka, Czech Republic, May 14 17, 2012 Paper #195 NUMERICAL SIMULATION OF TRANSITIONAL FLOWS WITH LAMINAR KINETIC ENERGY
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 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 informationComparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
Proceedings of 4 th ICCHMT May 17-0, 005, Paris-Cachan, FRANCE 381 Comparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
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 informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal and Vijaya Rani PI: Sarma L. Rani Department of Mechanical and Aerospace
More informationA 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 informationVALIDATION OF LES FOR LOCAL HEAT ENVIRONMENT IN TOKYO -COMPARISON WITH FIELD MEASUREMENT DATA-
The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 29, Taipei, Taiwan VALIDATION OF LES FOR LOCAL HEAT ENVIRONMENT IN TOKYO -COMPARISON WITH FIELD MEASUREMENT DATA- Tsuyoshi Nozu 1,
More informationAerodynamic Analysis of wind turbine
Aerodynamic Analysis of wind turbine Master s thesis in Fluid Mechanics AYYOOB ZARMEHRI Department of Applied Mechanics Division of Fluid dynamics CHALMERS UNIVERSITY OF TECHNOLOGY Gothenburg, Sweden 2012
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More informationThe effect of geometric parameters on the head loss factor in headers
Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract
More informationChuichi Arakawa Graduate School of Interdisciplinary Information Studies, the University of Tokyo. Chuichi Arakawa
Direct Numerical Simulations of Fundamental Turbulent Flows with the Largest Grid Numbers in the World and its Application of Modeling for Engineering Turbulent Flows Project Representative Chuichi Arakawa
More informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
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 informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationNUMERICAL SIMULATION AND MODELING OF UNSTEADY FLOW AROUND AN AIRFOIL. (AERODYNAMIC FORM)
Journal of Fundamental and Applied Sciences ISSN 1112-9867 Available online at http://www.jfas.info NUMERICAL SIMULATION AND MODELING OF UNSTEADY FLOW AROUND AN AIRFOIL. (AERODYNAMIC FORM) M. Y. Habib
More informationABSTRACT OF ONE-EQUATION NEAR-WALL TURBULENCE MODELS. Ricardo Heinrich Diaz, Doctor of Philosophy, 2003
ABSTRACT Title of dissertation: CRITICAL EVALUATION AND DEVELOPMENT OF ONE-EQUATION NEAR-WALL TURBULENCE MODELS Ricardo Heinrich Diaz, Doctor of Philosophy, 2003 Dissertation directed by: Professor Jewel
More informationDirect Numerical Simulation of the Neutrally Stratified Turbulent Ekman Boundary Layer
Journal of the Earth Simulator, Volume 6, October 26, 3 15 Direct Numerical Simulation of the Neutrally Stratified Turbulent Ekman Boundary Layer Katsuhiro Miyashita 1, Kaoru Iwamoto 1 and Hiroshi Kawamura
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 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 informationLecture 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 informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationUNCLASSIFIED UNCLASSIFIED
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADPO 13706 TITLE: Energy- Containing-Range Modeling of Fully-Developed Channel Flow Using a Hybrid RANS/LES Technique DISTRIBUTION:
More informationIntroduction to CFD modelling of source terms and local-scale atmospheric dispersion (Part 1 of 2)
1 Introduction to CFD modelling of source terms and local-scale atmospheric dispersion (Part 1 of 2) Atmospheric Dispersion Modelling Liaison Committee (ADMLC) meeting 15 February 2018 Simon Gant, Fluid
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More 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 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 informationDNS of Buoyancy Driven Flow Inside a Horizontal Coaxial Cylinder
DNS of Buoyancy Driven Flow Inside a Horizontal Coaxial Cylinder Imama Zaidi 1, Yacine Addad 2, Dominique Laurence 1 1 The University of Manchester, School of Mechanical, Aerospace and Civil Eng., M60
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationFluctuation dynamo amplified by intermittent shear bursts
by intermittent Thanks to my collaborators: A. Busse (U. Glasgow), W.-C. Müller (TU Berlin) Dynamics Days Europe 8-12 September 2014 Mini-symposium on Nonlinear Problems in Plasma Astrophysics Introduction
More informationINTER-COMPARISON AND VALIDATION OF RANS AND LES COMPUTATIONAL APPROACHES FOR ATMOSPHERIC DISPERSION AROUND A CUBIC OBSTACLE. Resources, Kozani, Greece
INTER-COMPARISON AND VALIDATION OF AND LES COMPUTATIONAL APPROACHES FOR ATMOSPHERIC DISPERSION AROUND A CUBIC OBSTACLE S. Andronopoulos 1, D.G.E. Grigoriadis 1, I. Mavroidis 2, R.F. Griffiths 3 and J.G.
More informationDNS, LES, and wall-modeled LES of separating flow over periodic hills
Center for Turbulence Research Proceedings of the Summer Program 4 47 DNS, LES, and wall-modeled LES of separating flow over periodic hills By P. Balakumar, G. I. Park AND B. Pierce Separating flow in
More informationWall treatments and wall functions
Wall treatments and wall functions A wall treatment is the set of near-wall modelling assumptions for each turbulence model. Three types of wall treatment are provided in FLUENT, although all three might
More informationProcess Chemistry Toolbox - Mixing
Process Chemistry Toolbox - Mixing Industrial diffusion flames are turbulent Laminar Turbulent 3 T s of combustion Time Temperature Turbulence Visualization of Laminar and Turbulent flow http://www.youtube.com/watch?v=kqqtob30jws
More informationNovel uses of DNS with turbulent separation for RANS models
Novel uses of DNS with turbulent separation for RANS models Gary Coleman*, Chris Rumsey* and Philippe Spalart** *NASA Langley Research Center **Boeing Commercial Airplanes 1 Outline q Flow: spatial separation
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
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 informationA Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case
A Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case M. Weng, P. V. Nielsen and L. Liu Aalborg University Introduction. The use of CFD in ventilation research has arrived
More information