STRESS TRANSPORT MODELLING 2
|
|
- Lucy Turner
- 6 years ago
- Views:
Transcription
1 STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1
2 Introduction In the previous lecture we introduced the stress transport equations, and some relatively simple models for closing them. We noted that the generation terms are represented exactly, and that in many flows the pressure-correlations are then the most important terms that require modelling. The (linear) models developed for the pressure correlations were based on simple return-to-isotropy ideas. They required additional corrections (involving wall-distances) to account for wall-proximity effects. In this lecture we look at more advanced approaches. In particular, how models can be built to ensure they satisfy the 2-component turbulence limit. STRESS TRANSPORT MODELLING 2 p.2
3 2-Component Turbulence 2-component turbulence can be found at fluid-fluid interfaces. Velocity fluctuations normal to surface vanish, but those parallel to it need not. Consequently, v 2 = 0, but u 2 and w 2 are non-zero. Fluid 2 Fluid 1 y x Near-wall turbulence is close to 2-component. Buoyancy/rotation effects can drive turbulence towards 2- component limit. Channel flow DNS of Kim et al (1987) Although many flows not close to 2-component limit, a TCL model should have better physics built into it, and might be hoped to perform better than simpler models even away from this limit. STRESS TRANSPORT MODELLING 2 p.3
4 Realizability Realizability (Schumann 1977, Lumley 1978): ensures model cannot predict unphysical values eg. negative normal stresses. Applied by ensuring model behaves correctly when one stress vanishes ie. in 2-component limit. At RSM level one can control the rate of change of stresses in the 2-component limit. A necessary condition for realizability is that Dv 2 /Dt = 0 when v 2 = 0: to prevent v 2 being driven to negative values. Shih & Lumley (1985) developed realizable models for φ ij, φ iθ in the second-moment equations. Modelling at UMIST has followed similar lines, but used slightly different constraints for reasons given later. This lecture describes methods of including 2-component limit in stress transport modelling, illustrated by model development done at UMIST. STRESS TRANSPORT MODELLING 2 p.4
5 Stress Transport Equations Transport equation for u i u j can be written: Du i u j Dt = P ij + φ ij + d ij ε ij Production P ij = u i u k U j x k u j u k U i x k Diffusion typically modelled by simple GGDH: d ij = x k [( ) ] k νδ lk + c s ε u ui u j ku l x l Modelling of pressure-strain (φ ij ) and dissipation (ε ij ) will be considered in this lecture. STRESS TRANSPORT MODELLING 2 p.5
6 Wall/Surface Asymptotes One can examine the immediate near-wall behaviour by writing Taylor series expansions: u 1 = a 1 + b 1 y + c 1 y 2 + u 2 = a 2 + b 2 y + c 2 y 2 + u 3 = a 3 + b 3 y + c 3 y 2 + At solid wall, a 1 = a 2 = a 3 = 0. Furthermore, continuity leads to b 2 = 0. This leads to u 2 1, u2 3 y2, but u 2 2 y4. The near free-surface behaviour and limiting forms of all terms in the transport equations can be examined in a similar fashion. STRESS TRANSPORT MODELLING 2 p.6
7 Wall Asymptotes At a solid wall, for example, one can show that: d ν ij φ ij + d p ij ε ij u 2 2νb νb 1c 1 y + 4νb 1 c 1 y + 2νb νb 1c 1 y + v 2 12νc 2 2 y2 + 4νc 2 2 y2 + 8νc 2 2 y2 + w 2 2νb νb 3c 3 y + 4νb 3 c 3 y + 2νb νb 3c 3 y + uv 6νb 1 c 2 y 2 + 2νb 1 c 2 y + 4νb 1 c 2 y + STRESS TRANSPORT MODELLING 2 p.7
8 Anisotropy Invariants Near the 2-component state, turbulence has different characteristics than in regions close to isotropy. Anisotropy invariants can be used to give some information on turbulence structure: A 2 = a ij a ij A 3 = a ij a jk a ki a ij = u i u j /k (2/3)δ ij Realizable states lie within a triangle in A 2 -A 3 space. Lumley s flatness parameter A = 1 9/8(A 2 A 3 ). In isotropic turbulence A = 1 In 2-component turbulence A = 0 STRESS TRANSPORT MODELLING 2 p.8
9 Constructing TCL Models One approach might be to use A (or other invariants) to make coefficients in the linear models vanish in the 2-component limit. But one cannot get different limiting behaviour for different components this way. Alternatively, we can construct higher-order tensor forms of models, and choose the resulting model coefficients in order to give desired behaviour. The pressure-strain is often the most influential modelled term: φ ij p ρ ( ui + u ) j x j x i STRESS TRANSPORT MODELLING 2 p.9
10 Modelling φ ij As outlined earlier, one can obtain a Poisson equation for fluctuating pressure p: 1 ρ p x 2 i = 2 x l x k (u l u k u l u k ) 2 U k x l u l x k Solving this in terms of a Greens function gives an expression for φ ij : φ ij = 1 ( ) 3 u k u l u i + 3 u k u l u j dv 4π V ol r l r k r j r l r k r i r }{{} φ ij1 1 ( ) 2 u l u i + 2 u l u j U k dv 2π V ol r k r j r k r i r l r }{{} φ ij2 STRESS TRANSPORT MODELLING 2 p.10
11 Modelling φ ij2 Assuming the velocity gradient varies reasonably slowly, one can model φ ij2 as φ ij2 = ( ) Xkj li + X lj Uk ki where Xkj li = 1 x l 2π V ol 2 u l u i r k r j dv r The tensor X li kj can then be modelled in terms of the Reynolds stresses. A general linear tensor form (with required symmetries) is: X li kj = αu l u i δ kj + β (u i u j δ lk + u l u j δ ik + u l u k δ ij + u i u k δ lj ) +γu k u j δ il + ξk (δ lj δ ik + δ lk δ ij ) + ηkδ il δ kj From the exact integral, we want to ensure that: Continuity: Xki li = 0, li Normalization: Xkk = 2u lu i STRESS TRANSPORT MODELLING 2 p.11
12 A Linear φ ij2 Model Applying these constraints leads to (Launder et al 1975) φ ij2 = γ (P ij (1/3)δ ij P kk ) 30γ 2 ( Ui k + U ) j 55 x j x i 8γ 2 11 (D ij (1/3)δ ij D kk ) where D ij = u i u k U k / x j u j u k U k / x i. However, one cannot satisfy 2-component limit for any value of the free coefficient γ. STRESS TRANSPORT MODELLING 2 p.12
13 Non-Linear φ ij2 Modelling 1 The two-point correlations in integral for Xkj li functions of the stresses. One might thus include higher order terms in X li kj : are not simple linear X li kj k = λ 1δ li δ kj + λ 2 (δ lj δ ki + δ lk δ ij ) + λ 3 a li δ kj + + λ 10 (a lm a mj δ ki + a lm a mk δ ij + a im a mj δ lk + a km a mi δ lj ) + + λ 13 a li a km a mj + λ 14 a kj a lm a mi + + λ 20 a mn a mn (a lj δ ki + a lk δ ij + a ij δ lk + a ki δ lj ) Inclusion of cubic terms allows the 2-component limit to be satisfied: φ 222 = 2Xk2 l2 U k = 0 when u 2 2 x = 0 l Application of this, together with continuity and normalization, yields a system of equations for the λ s, which when solved gives a model with 2 free coefficients. STRESS TRANSPORT MODELLING 2 p.13
14 Non-Linear φ ij2 Modelling 2 The resulting model can be written: φ ij2 = 0.6 (P ij (1/3)δ ij P kk ) + 0.3a ij P kk [ [ uk u j u l u i Uk U ] l u lu k k x l x k k [ ]] U j U i u i u k + u j u k x l x l c 2 [A 2 (P ij D ij ) + 3a mi a nj (P mn D mn )] {( 7 + c 2 15 A ) 2 (P ij δ ij P kk /3) [a ij 0.5 (a ik a kj δ ij A 2 /3)] P kk 4 [( ui u m 0.05a ij a lk P kl k P mj + u ) ] ju m k P u l u m mi (2/3)δ ij k P ml [ ] [ [ ul u i u k u j u l u m u k u m Ul k 2 (1/3)δ ij k 2 6D lk + 13k + U ]] k x k x l u } lu i u k u j k 2 (D lk P lk ) Fu (1988) set c 2 = 0. Later work at UMIST has included non-zero values for both c 2 and c 2, to improve wall-bounded flows. STRESS TRANSPORT MODELLING 2 p.14
15 Modelling φ ij1 Including non-linear terms in Rotta s (1951) return to isotropy model: φ ij1 = c 1 εa ij c 1ε(a ik a kj A 2 δ ij /3) c 1εA 2 a ij Inclusion of a cubic term does not increase generality since the Cayley Hamilton theorem gives (a ik a kl a lj A 3 δ ij /3) = A 2 a ij /2. One can apply the 2-component limit behaviour in a similar manner to φ ij2 : finding coefficients that give φ 221 = 0 when u 2 2 = 0. However, the resultant model does not perform well in shear flows. An alternative is to make the coefficients functions of A, to ensure they vanish in the 2-component limit. The current UMIST form is: φ ij1 = c 1 ε [a ij + c 1(a ij a jk A 2 δ ij /2)] f Aεa ij with c 1 = 3.1(A 2 A) 1/2, c 1 = 1.1 and f A = A1/2. STRESS TRANSPORT MODELLING 2 p.15
16 Dissipation Modelling 1 In high-re-number flow, ε ij is normally assumed isotropic: ε ij = (2/3)εδ ij However, this form does not satisfy the 2-component limit (at a wall ε 22 = 0). A form ε ij = εu i u j /k does vanish if u 2 2 = 0. One can combine these two forms to give with f ε a suitable function of A. This does give ε 22 = 0 when u 2 2 = 0. ε ij = (2/3)εδ ij f ε + u iu j k ε(1 f ε) But note that it still does not satisfy the exact wall or surface-limiting behaviour for all components. STRESS TRANSPORT MODELLING 2 p.16
17 Dissipation Rate Modelling 2 The dissipation rate ε is typically obtained from its own separate transport equation: Dε Dt = c εp kk ε1 2k c ε2 ε 2 k + x k [( ) ] k ε νδ lk + c ε ε u ku l x l Anisotropy invariants can be introduced into the coefficients to make the model respond to different types of turbulence structure. The current UMIST form is: c ε2 = 1.92/( A 1/2 2 A) c ε1 = 1.0 STRESS TRANSPORT MODELLING 2 p.17
18 Scalar Flux Modelling In some flows it is sufficient to employ a simple gradient-diffusion type of model for the scalar fluxes. In others one might need to employ a full second-moment closure. A similar treatment can be applied to the scalar flux equations to ensure realizability. The transport equation for u i θ can be written: Du i θ Dt = P iθ + φ iθ + d iθ ε iθ Production is now P iθ = u i u k T x k u k θ U i x k Diffusion can often conveniently be modelled via the GGDH. The isotropic dissipation assumption gives ε iθ = 0. The only term remaining is the pressure-scalar gradient correlation φ iθ. STRESS TRANSPORT MODELLING 2 p.18
19 Pressure-Scalar Gradient 1 In similar fashion to φ ij, we can obtain φ iθ = 1 3 u k u l θ dv 1 2 u l θ U k dv 4π V ol r l r k r i r 2π V ol r k r i r l r }{{}}{{} φ iθ1 φ iθ2 The rapid part, φ iθ2, can be modelled as φ iθ2 = 2b l ki U k x l where b l ki = 1 4π V ol 2 u l θ r k r i dv r Linearity principle: φ iθ is linear in the scalar, so b l ki stresses, but only linear in scalar fluxes. can be non-linear in We want the model for b l ki to satisfy: Continuity: b k ki = 0, Normalization: bl kk = u lθ. STRESS TRANSPORT MODELLING 2 p.19
20 A Linear φ iθ Model A linear form for b l ki could be written b l ki = α 1 u l θδ ik + α 2 ( uk θδ li + u i θδ lk ) Application of the above constraints leads to the linear QI model (Launder, 1973): φ iθ2 = 0.8u k θ U i x k 0.2u k θ U k x i But this does not satisfy the 2-component limit. STRESS TRANSPORT MODELLING 2 p.20
21 Non-Linear φ iθ2 Model 1 Including all cubic terms in b l ki leads to: b l ki = α 1 u l θδ ik + α 2 ( uk θδ li + u i θδ lk ) + α3 u l θa ik + α 4 ( uk θa li + u i θa lk ) + + α 7 u m θa ml a ik + α 8 u m θ (a mk a il + a mi a kl ) + + u n θa mn (α 13 a ml δ ik + α 14 (a mk δ li + a mi δ lk )) To apply the 2-component limit: Shih & Lumley (1985) employed the Schwarz inequality constraint: ( uα θ ) 2 u 2 α θ 2. This leads to both free coefficients in φ ij being determined as zero. Instead, at UMIST, have ensured that nett contribution P 2θ + φ 2θ = 0 when u 2 = 0. This results in the condition U k x l b l k2 = 1 2 u lθ U 2 x l when u 2 = 0 STRESS TRANSPORT MODELLING 2 p.21
22 Non-Linear φ iθ2 Model 2 Application of continuity, normalization, and 2-component limit results in all coefficients in b l ki being determined: φ iθ2 = 0.8u k θ U i 0.2u k θ U k + (1/3) ɛ x k x i k u iθ P ɛ ) 0.4u k θa il ( Uk x l + U l x k + 0.1u k θa ik a ml ( Um x l + U l x m 0.1u k θ (a im P mk + 2a mk P im ) /k ( Uk a ml + U ) l (amku i θ a miu k θ ) x l x k 0.05a ml [ 7a mk ( u i θ U k x l ) + u k θ U i x l ) u k θ ( )] U i U i a ml + a mk x k x l STRESS TRANSPORT MODELLING 2 p.22
23 Non-Linear φ iθ1 Model A non-linear extension of Monin s (1965) return to isotropy for φ iθ1 : ε φ iθ1 = c θ1 k u iθ c ε θ1 k a iju j θ c ε θ1 k a ika kj u j θ c ε θ1 k A 2u i θ This satisfies realizability for any coefficient values. We found it beneficial to include some mean scalar gradient term for shear flows with strong scalar gradients: φ iθ1 = c θ1 ε k c θ1t Rka ij T x j [ ui θ(1 + c θ1a 2 ) + c θ1a ik u k θ + c θ1a ik a kj u j θ ] where c θ1 = 1.7R 1/2 [ (A 2 A) 1/2] c θ1 = 0.8 c θ1 = 1.1 c θ1 = 0.6 c θ1t = 0.2A 1/2 STRESS TRANSPORT MODELLING 2 p.23
24 Homogeneous Flows 1 Stress anisotropy development in simple shear flow: Note that the basic RSM (from the previous lecture) would return u 2 2 = u2 3. STRESS TRANSPORT MODELLING 2 p.24
25 Scalar fluxes in simple shear flow: Weakly Strained Homogeneous Flows 2 Strongly Strained : Basic RSM : TCL STRESS TRANSPORT MODELLING 2 p.25
26 Free Jets Experiments show that an axisymmetric free jet spreads more slowly than a plane one. Standard models (k-ε and the basic RSM) predict the opposite. Dynamic Field Thermal Field Plane Jet Round Jet Plane Jet Round Jet Experiments Basic RSM TCL Model The TCL model does give the correct relative spreading rates. STRESS TRANSPORT MODELLING 2 p.26
27 Summary In this lecture we have seen how the 2-component limit can be built into second-moment models by ensuring correct behaviour of transport equation when one stress vanishes. Anisotropy invariants can be used to enforce the limit and to make model respond to different types of turbulence. Examples of the development of TCL models for stresses and scalar fluxes have been outlined. It has been shown that TCL models can bring predictive improvements, even in free flows. Further applications will be shown in seminars next week. STRESS TRANSPORT MODELLING 2 p.27
28 References Fu, S., 1988 Computational Modelling of Turbulent Swirling Flows with Second-Moment Closures, PhD Thesis, Faculty of Technology, University of Manchester. Kim, J., Moin, P., Moser, R., 1987 Turbulence statistics in fully developed channel flow at low Reynolds number, J. Fluid Mech., vol 177, pp Launder, B.E., 1973 Scalar property transport by turbulence, Report HTS/73/26, Mech. Eng. Dept., Imperial College, London. Launder, B.E., Reece, G.J., Rodi, W., 1975 Progress in the development of a Reynolds stress turbulence closure, J. Fluid Mech., vol 68, p 537. Lumley, J., 1978 Computational modelling of turbulent flows, Adv. Appl. Mech., vol 18, pp Monin, A.S., 1965 On the symmetry of turbulence in the surface layer of air, Izv. Atm. Oceanic Phys., vol 1, p 45. Rotta, J., 1951 Statistische Theorie nichthomogener Turbulenz, Zeitschrift für Physik, vol 129, pp Schumann, U., 1977, Realizability of Reynolds stress turbulence models, Phys. Fluids, vol 20, pp Shih, T-S., Lumley, J.L., 1985 Modeling of pressure correlation terms in Reynolds stress and scalar flux equations, Report FD-85-03, Sibley School of Mechanical and Aerospace Eng., Cornell University. STRESS TRANSPORT MODELLING 2 p.28
CHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More 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 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 boos Built by scientists, for scientists 3,900 116,000 120M Open access boos available International authors and editors Downloads Our authors
More informationPerspectives on Reynolds Stress Modeling Part I: General Approach. Umich/NASA Symposium on Advances in Turbulence Modeling
DLR.de Chart 1 Perspectives on Reynolds Stress Modeling Part I: General Approach Umich/NASA Symposium on Advances in Turbulence Modeling Bernhard Eisfeld, Tobias Knopp 11 July 2017 DLR.de Chart 2 Overview
More informationProbability density function and Reynolds-stress modeling of near-wall turbulent flows
Probability density function and Reynolds-stress modeling of near-wall turbulent flows Thomas D. Dreeben and Stephen B. Pope Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca,
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
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 information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationElliptic relaxation for near wall turbulence models
Elliptic relaxation for near wall turbulence models J.C. Uribe University of Manchester School of Mechanical, Aerospace & Civil Engineering Elliptic relaxation for near wall turbulence models p. 1/22 Outline
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 informationOn modeling pressure diusion. in non-homogeneous shear ows. By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3
Center for Turbulence Research Proceedings of the Summer Program 1996 63 On modeling pressure diusion in non-homogeneous shear ows By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3 New models
More informationFUNDAMENTAL AND CONCEPTUAL ASPECTS OF TURBULENT FLOWS
FUNDAMENTAL AND CONCEPTUAL ASPECTS OF TURBULENT FLOWS Arkady Tsinober Professor and Marie Curie Chair in Fundamental and Conceptual Aspects of Turbulent Flows Institute for Mathematical Sciences and Department
More informationReynolds stress budgets in Couette and boundary layer flows
Reynolds stress budgets in Couette and boundary layer flows By Jukka Komminaho and Martin Skote Dept. of Mechanics, KTH, SE-1 44, Stockholm, Sweden Reynolds stress budgets for both Couette and boundary
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 informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More 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 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 informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
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 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 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 informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationA MATHEMATICAL MODEL OF DRIP EMITTER DISCHARGE DEPENDING ON THE GEOMETRIC PARAMETERS OF A LABYRINTH CHANNEL *
11 th National Congress on Theoretical and Applied Mechanics, -5 Sept. 009, Borovets, Bulgaria A MATHEMATICAL MODEL OF DRIP EMITTER DISCHARGE DEPENDING ON THE GEOMETRIC PARAMETERS OF A LABYRINTH CHANNEL
More informationEuclidean invariance and weak-equilibrium condition for the algebraic Reynolds stress model
J. Fluid Mech. 2006), vol. 569, pp. 399 408. c 2006 Cambridge University Press doi:10.1017/s0022112006003041 Printed in the United Kingdom 399 Euclidean invariance and weak-equilibrium condition for the
More informationCentro de Estudos de Fenómenos de Transporte, FEUP & Universidade do Minho, Portugal
DEVELOPING CLOSURES FOR TURBULENT FLOW OF VISCOELASTIC FENE-P FLUIDS F. T. Pinho Centro de Estudos de Fenómenos de Transporte, FEUP & Universidade do Minho, Portugal C. F. Li Dep. Energy, Environmental
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationNumerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model
Numerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model Hitoshi Sugiyama 1),Naoto Kato 1), Masahiro Ouchi 1) 1) Graduate School of Engineering,Utsunomiya University
More informationdissipation (P=") in homogeneous shear ow corresponding to the initial conditions (b 11 ) 0 =
Figure 11. Computed model predictions for the time evolution of the ratio of production to dissipation (P=") in homogeneous shear ow corresponding to the initial conditions (b 11 ) 0 = 0:32, (b 22 ) 0
More informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
More informationMATH45061: SOLUTION SHEET 1 V
1 MATH4561: SOLUTION SHEET 1 V 1.) a.) The faces of the cube remain aligned with the same coordinate planes. We assign Cartesian coordinates aligned with the original cube (x, y, z), where x, y, z 1. The
More informationDynamics of the Coarse-Grained Vorticity
(B) Dynamics of the Coarse-Grained Vorticity See T & L, Section 3.3 Since our interest is mainly in inertial-range turbulence dynamics, we shall consider first the equations of the coarse-grained vorticity.
More informationMODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION
MODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION F. Dehoux Fluid Mechanics, Power generation and Environment Department MFEE Dept.) EDF R&D Chatou, France frederic.dehoux@edf.fr
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 informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationOn Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers
Continuum Mech. Thermodyn. (1996) 8: 247 256 On Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers I-Shih Liu Instituto de Matemática Universidade do rio de Janeiro, Caixa Postal
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More informationIn this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case
Section.. Thermoelasticity In this section, thermoelasticity is considered. By definition, the constitutive relations for F, θ, Gradθ. This general case such a material depend only on the set of field
More informationNear-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Near-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods Axel Probst (now at: C 2 A 2 S 2 E, DLR Göttingen) René Cécora,
More informationHydrodynamic Characteristics of Gradually Expanded Channel Flow
International Journal of Hydraulic Engineering 2013, 2(2): 21-27 DOI: 10.5923/j.ijhe.20130202.01 Hydrodynamic Characteristics of Gradually Expanded Channel Flow Edward Ching-Rue y LUO Department of Civil
More informationInterpretation of the mechanism associated with turbulent drag reduction in terms of anisotropy invariants
J. Fluid Mech. (27), vol. 577, pp. 457 466. c 27 Cambridge University Press doi:1.117/s221127583 Printed in the United Kingdom 457 Interpretation of the mechanism associated with turbulent drag reduction
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationAnalysis and modeling of buoyancy generated turbulence using numerical data
Analysis and modeling of buoyancy generated turbulence using numerical data Sharath S. Girimaji Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, VA
More informationRANS simulations of rotating flows
Center for Turbulence Research Annual Research Briefs 1999 257 RANS simulations of rotating flows By G. Iaccarino, A. Ooi, B. A. Pettersson Reif AND P. Durbin 1. Motivation and objectives Numerous experimental
More informationROTATION OF TRAJECTORIES IN LAGRANGIAN STOCHASTIC MODELS OF TURBULENT DISPERSION
ROTATION OF TRAJECTORIES IN LAGRANGIAN STOCHASTIC MODELS OF TURBULENT DISPERSION B.L. SAWFORD CSIRO Atmospheric Research, Aspendale, Vic 3195, Australia Received in final form 27 July 1999) Abstract. We
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 informationTurbulent Flows. quiescent surroundings W U V. r U. nozzle. fluid supply CHAPTER 5: FREE SHEAR FLOWS
quiescent surroundings x W U V r U J θ d nozzle fluid supply Figure 5.1: Sketch of a round jet experiment, showing the polarcylindrical coordinate system employed. 0.2 x/d = 30 U /U J 0.1 x/d = 60 x/d
More informationThe return to isotropy of homogeneous turbulence
J. Fluid Mech. (21), vol. 436, pp. 59 84. Printed in the United Kingdom c 21 Cambridge University Press 59 The return to isotropy of homogeneous turbulence By KWING-SO CHOI 1 AND JOHN L. LUMLEY 2 1 School
More informationCartesian Tensors. e 2. e 1. General vector (formal definition to follow) denoted by components
Cartesian Tensors Reference: Jeffreys Cartesian Tensors 1 Coordinates and Vectors z x 3 e 3 y x 2 e 2 e 1 x x 1 Coordinates x i, i 123,, Unit vectors: e i, i 123,, General vector (formal definition to
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 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 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 informationFLOW-NORDITA Spring School on Turbulent Boundary Layers1
Jonathan F. Morrison, Ati Sharma Department of Aeronautics Imperial College, London & Beverley J. McKeon Graduate Aeronautical Laboratories, California Institute Technology FLOW-NORDITA Spring School on
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationExam in Fluid Mechanics 5C1214
Eam in Fluid Mechanics 5C1214 Final eam in course 5C1214 13/01 2004 09-13 in Q24 Eaminer: Prof. Dan Henningson The point value of each question is given in parenthesis and you need more than 20 points
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 informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
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 information2.2 The Turbulent Round Jet
Canonical Turbulent Flows 13. The Turbulent Round Jet Jet flows are a subset of the general class of flows known as free shear flows where free indicates that the shear arises in the absence of a boundary
More informationProblem Set 4. Eshelby s Inclusion
ME34B Elasticity of Microscopic Structures Wei Cai Stanford University Winter 24 Problem Set 4. Eshelby s Inclusion Chris Weinberger and Wei Cai c All rights reserved Problem 4. (5 Spherical inclusion.
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 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 informationCOMMUTATION ERRORS IN PITM SIMULATION
COMMUTATION ERRORS IN PITM SIMULATION B. Chaouat ONERA, 93 Châtillon, France Bruno.Chaouat@onera.fr Introduction Large eddy simulation is a promising route. This approach has been largely developed in
More 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 informationA SEAMLESS HYBRID RANS/LES MODEL WITH DYNAMIC REYNOLDS-STRESS CORRECTION FOR HIGH REYNOLDS
A SEAMS HYBRID RANS/ MODEL WITH DYNAMIC REYNOLDS-STRESS CORRECTION FOR HIGH REYNOLDS NUMBER FLOWS ON COARSE GRIDS P. Nguyen 1, J. Uribe 2, I. Afgan 1 and D. Laurence 1 1 School of Mechanical, Aerospace
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 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 informationCFD ANALYSIS OF TURBULENT THERMAL MIXING OF HOT AND COLD AIR IN AUTOMOBILE HVAC UNIT
ISTP-6, 005, PRAGUE 6 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA CFD ANALYSIS OF TURBULENT THERMAL MIING OF HOT AND COLD AIR IN AUTOMOBILE HVAC UNIT Hideo Asano ((, Kazuhiko Suga (3, Masafumi Hirota
More informationNumerical modeling of complex turbulent flows
ISTP-16, 5, PRAGUE 16 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA Numerical modeling of complex turbulent flows Karel Kozel Petr Louda Jaromír Příhoda Dept. of Technical Mathematics CTU Prague, Karlovo
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 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 informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationEvaluation Study of Pressure-Strain Correlation Models in Compressible Flow
Journal of Applied Fluid Mechanics, Vol. 9, No. 6, pp. 2685-2693, 2016. Available online at www.jafmonline.net, SSN 1735-3572, ESSN 1735-3645. Evaluation Study of ressure-strain Correlation Models in Compressible
More informationANALYSIS OF PERIODICALLY-FORCED HOMOGENEOUS TURBULENCE IN THE RAPID DISTORTION LIMIT. A Thesis JOSHUA R. O NEIL
ANALYSIS OF PERIODICALLY-FORCED HOMOGENEOUS TURBULENCE IN THE RAPID DISTORTION LIMIT A Thesis by JOSHUA R. O NEIL Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment
More informationAlgebraic proof and application of Lumley s realizability triangle
Algebraic proof and application of Lumley s realizability triangle G.A. Gerolymos and I. Vallet Sorbonne Universités, Université Pierre-et-Marie-Curie (UPMC), place Jussieu, 755 Paris, France Emails: georges.gerolymos@upmc.fr,
More informationTurbulence Solutions
School of Mechanical, Aerospace & Civil Engineering 3rd Year/MSc Fluids Turbulence Solutions Question 1. Decomposing into mean and fluctuating parts, we write M = M + m and Ũ i = U i + u i a. The transport
More informationO. A Survey of Critical Experiments
O. A Survey of Critical Experiments 1 (A) Visualizations of Turbulent Flow Figure 1: Van Dyke, Album of Fluid Motion #152. Generation of turbulence by a grid. Smoke wires show a uniform laminar stream
More informationEXPLICIT ALGEBRAIC MODELS FOR STRATIFIED FLOWS
EXPLICIT ALGEBRAIC MODELS FOR STRATIFIED FLOWS W.M.J. Lazeroms, G. Brethouwer, S. Wallin,2, and A.V. Johansson Linné FLOW Centre, KTH Mechanics, SE- 44 Stockholm, Sweden 2 Swedish Defence Research Agency
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationInitial Efforts to Improve Reynolds Stress Model Predictions for Separated Flows
Initial Efforts to Improve Reynolds Stress Model Predictions for Separated Flows C. L. Rumsey NASA Langley Research Center P. R. Spalart Boeing Commercial Airplanes University of Michigan/NASA Symposium
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2011 Outline Outline Part I First Lecture Connection between time and ensemble average Ergodicity1 Ergodicity
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
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 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 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 informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
More informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
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 information(U c. t)/b (U t)/b
DYNAMICAL MODELING OF THE LARGE-SCALE MOTION OF A PLANAR TURBULENT JET USING POD MODES. S. Gordeyev 1 and F. O. Thomas 1 University of Notre Dame, Notre Dame, USA University of Notre Dame, Notre Dame,
More informationRica Mae Enriquez*, Robert L. Street, Francis L. Ludwig Stanford University, Stanford, CA. 0 = u x A u i. ij,lass. c 2 ( P ij. = A k. P = A ij.
P1.45 ASSESSMENT OF A COUPLED MOMENTUM AND PASSIVE SCALAR FLUX SUBGRID- SCALE TURBULENCE MODEL FOR LARGE-EDDY SIMULATION OF FLOW IN THE PLANETARY BOUNDARY LAYER Rica Mae Enriquez*, Robert L. Street, Francis
More informationCapturing turbulent density flux effects in variable density flow by an explicit algebraic model
Capturing turbulent density flux effects in variable density flow by an explicit algebraic model I. A. Grigoriev, S. Wallin, G. Brethouwer, and A. V. Johansson Citation: Physics of Fluids 1994-present
More informationUncertainty quantification for RANS simulation of flow over a wavy wall
Uncertainty quantification for RANS simulation of flow over a wavy wall Catherine Gorlé 1,2,3, Riccardo Rossi 1,4, and Gianluca Iaccarino 1 1 Center for Turbulence Research, Stanford University, Stanford,
More informationDivergence free synthetic eddy method for embedded LES inflow boundary condition
R. Poletto*, A. Revell, T. Craft, N. Jarrin for embedded LES inflow boundary condition University TSFP Ottawa 28-31/07/2011 *email: ruggero.poletto@postgrad.manchester.ac.uk 1 / 19 SLIDES OVERVIEW 1 Introduction
More informationDIVERGENCE FREE SYNTHETIC EDDY METHOD FOR EMBEDDED LES INFLOW BOUNDARY CONDITIONS
DIVERGECE FREE SYTHETIC EDDY METHOD FOR EMBEDDED LES IFLOW BOUDARY CODITIOS R. Poletto, A. Revell, T. Craft,. Jarrin 1 School of Mechanical Aerospace and Civil Engineering University of Manchester, Manchester,
More informationJournal of Fluid Science and Technology
Science and Technology Effect of Molecular Diffusivities on Countergradient Scalar Transfer in a Strong Stable Stratified Flow (Study on the Linear and Nonlinear Processes by using RDT) Kouji NAGATA, Takashi
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 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 informationFOUR-WAY COUPLED SIMULATIONS OF TURBULENT
FOUR-WAY COUPLED SIMULATIONS OF TURBULENT FLOWS WITH NON-SPHERICAL PARTICLES Berend van Wachem Thermofluids Division, Department of Mechanical Engineering Imperial College London Exhibition Road, London,
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 informationEVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS
EVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS A Aroussi, S Kucukgokoglan, S.J.Pickering, M.Menacer School of Mechanical, Materials, Manufacturing Engineering and
More information