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! Wall bounded turbulence! Second order closure! Direct Numerical Simulations! Large-ed simulations! Summary! Most engineering problems involve turbulent flows. Such flows involve are highly unstea and contain a large range of scales. However, in most cases the mean or average motion is well defined.! Above: A jet in a cross flow! Left: cross section of a jet! Flow over a sphere! Instantaneous flow past a sphere at R = 15,000.! The drag depends on the separation point! Dye in water shows a laminar boundary layer separating ahead of the equator and remains laminar for almost one radius. It then becomes unstable and quickly turns turbulent. ONERA photograph, Werle 1980. From "An Album of Fluid Motion," by Van Dyke, Parabolic Press.!
Instantaneous flow past a sphere at R = 30,000 with a trip wire The classical experiment of Prandtl and Wieselsberger is repeated here, using air bubbles in water. A wire hoop ahead of the equator trips the boundary layer. It becomes turbulent, so that it separates farther rearward than if it were laminar (compare with top photograph). The overall drag is thereby dramatically reduced, in a way that occurs naturally on a smooth sphere only at a Reynolds numbers ten times as great. ONERA photograph, Werle 1980. Examples of Reynolds numbers:! Flow around a 3 m long car at 100 km/hr:! Re = LU 3 27.78 = = 5.5 106 5 v 1.5 10 Flow around a 100 m long submarine at 10 km/hr:! Re = LU 100 2.78 = = 2.78 10 8 v 10 6 Kinematic viscosity! (~20 C)! Water ν = 10-6 m 2 /s! Air ν = 1.5 10-5 m 2 /s! 1km/hr = 0.27778 m/s! Water flowing though a 0.01 m diameter pipe with a velocity of 1 m/s! Re = LU v = 0.01 1 10 6 =10 4 It can be shown that for turbulent flow the ratio of the size of the smallest ed to the length scale of the problem! δ L O(Re 3 / δ 4 ) L O(Re 1/ 2 ) In 3D! In 2D! If about 10 grid points are needed for Re=10 (the driven cavity problem)! Re! 3d!! 2d! 10 3!~ 300 3! ~ 100 2! 10 4! ~ 2000 3! ~ 300 2! 10 5!~ 10000 3! ~ 1000 2! Reynolds Averaged Navier-Stokes (RANS): Only the averaged motion is computed. The effect of fluctuations is modeled! Large Ed Simulations (LES): Large scale motion is fully resolved but small scale motion is modeled! Direct Numerical Simulations (DNS): Every length and time scale is fully resolved! Largest computations today use about 4000 3 points! Reynolds Averaged Navier-Stokes Equations! To solve for the mean motion, we derive equations for the mean motion by averaging the Navier-Stokes equations. The velocities and other quantities are decomposed into the average and the fluctuation part! a = A + a'
a = A + a' Defining an averaging procedure that satisfies the following rules:! < a > = A < a' > = 0 < a + b > = A + B < ca > = ca < a > = A This will hold for spatial averaging, temporal averaging, and ensamble averaging! There are several ways to define the proper averages! For homogeneous turbulence we can use the space average! < a > = 1 L L adx 0 For stea turbulence flow we can use the time average! < a > = 1 T adt T 0 For the general case we use the ensemble average! < a > = a r (x,t) ensambles a = A + a' Start with the Navier-Stokes equations! t u + uu = 1 p + ν 2 u Decompose the pressure and velocity into mean and fluctuations:! u = U + u' p = P + p' Or, in general, for any dependant variable:! a = A + a' < a >=A < a' >=0 < ca >=ca < a >= A Applying the averaging to the Navier-Stokes equations results in:! t U + UU = - 1 P + ν 2 U + < u'u'> < u'u' > < u'v' > < u'w' > < u'u'>= < u'v' > < v'v' > < v'w' > < u'w' > < v'w' > < w'w' > Reynoldʼs stress tensor! Physical interpretation! < uv > Fast moving fluid particle! Slow moving fluid particle! Net momentum transfer due to velocity fluctuations! Closure:! Since we only have an equation for the mean flow, the Reynolds stresses must be related to the mean flow.! No rigorous process exists for doing this!! THE TURBULENCE PROBLEM!
Outline! Turbulence Modeling II! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions! Wall bounded turbulence! Second order closure! Direct Numerical Simulations! Large-ed simulations! Summary! Introduce the turbulent ed viscosity! Zero and One equation models! U < u'u'> ij = ν i T + U i j j where! = l 2 0 t 0 Zero equation models! One equation models! Prandtlʼ mixing length! = l 0 2 du Smagorinsky model! = l 2 0 2S ij S ij Baldvin-Lomaz model! ( ) 1/ 2 = l 2 0 ( ω i ω i ) 1/ 2 l 0 = κy S ij = 1 U i + U j 2 j i ω i = U i U j j i = k 1/ 2 t 0 Where k is obtained by an equation describing its temporal-spatial evolution! However, the problem with zero and one equation models is that t 0 and l 0 are not universal. Generally, it is found that a two equation model is the minimum needed for a proper description!
Two equation models! To characterize the turbulence it seems reasonable to start with a measure of the magnitude of the velocity fluctuations. If the turbulence is isotropic, the turbulent kinetic energy can be used:! ( ) k = 1 < u'u' > + < v'v' > + < w'w' > 2 The turbulent kinetic energy does, however, not distinguish between large and small eddies.! To distinguish between large and small eddies we need to introduce a new quantity that describe! Usually, the turbulent dissipation rate is used! ε ν u' i u' i j j Smaller eddies dissipate faster! Solve for the average velocity! t U + UU = 1 P + ( ν + ) 2 U Where the turbulent kinematic ed viscosity is given by! = C µ k 2 ε The exact k-equation is:! k t + U k U j = τ i ij ε + ν k 1 j j j j 2 u ' iu ' ' i u j where! τ ij = u ' ' i u j 1 p'u ' j The exact epsilon-equation is considerably more complex and we will not write it down here.! Both equations contain transport, dissipation and production terms that must be modeled! The general for for the equations for k and epsilon is:! k t + U k = D k k + production dissipation ε t + U ε = D ε ε + production dissipation These terms must be modeled! Closure involves proposing a form for the missing terms and optimizing free coefficients to fit experimental data!
The k-epsilon model! Dk Dt = + (ν + C U 2 ) k - τ i ij ε j Dε Dt = (ν + C ε 3 ) ε + C 4 k τ U i ε 2 ij C 5 j k Turbulent! transport! Production! Dissipation! Here! = C k 2 ε and! τ ij =< u ' ' i u j >= 2 3 kδ U ij ν i T + U j j i C 1 = 0.09; C 2 =1.0; C 3 = 0.769; C 4 =1.44; C 5 =1.92 Two major numerical difficulties! The equations may be stiff in some regions of the flow requiring very small time step. This can be overcome by an implicit scheme.! In reality k goes to zero at the walls. In simulations this usually takes place so close to the wall that it is not resolved by the grid. To overcome this we usually use a wall function or a damping function! Other two equation turbulence models:!!rng k-epsilon!!nonlinear k-epsilon!!k-enstrophy!!k-l o!!k-reciprocal time!!etc! Turbulent transport of energy and species concentrations is modeled in similar ways.! For temperature we have:! T t + ut = α 2 T u = U + u' T =< T > +T' < T > + U < T >= α 2 < T > < UT > t Gradient Transport Hypothesis:! < UT > α T < T > Spreading rates:! Model Predictions!!! exp! k-e!!cmott! Plane jet!0.10-0.11! 0.108!!0.102! Round jet!0.085-0.095! 0.116!!0.095! Mixing layer!0.13-0.17! 0.152!!0.154!
From: C.G. Speziale: Analytical Methods for the Development of Reynolds-stress closure in Turbulence. Ann Rev. Fluid Mech. 1991. 23: 107-157! From: C.G. Speziale: Analytical Methods for the Development of Reynolds-stress closure in Turbulence. Ann Rev. Fluid Mech. 1991. 23: 107-157! Results! Wall bounded turbulence! From: C.G. Speziale: Analytical Methods for the Development of Reynolds-stress closure in Turbulence. Ann Rev. Fluid Mech. 1991. 23: 107-157! Wall bounded turbulence! Fundamental assumption: determined by local variables only! Mean flow! Define a shear velocity:! v * = = µ du,, ν [ kg /ms 2 ], [ kg /m 3 ], m 2 /s [ ] Only the mean shear rate and the properties of the fluid are important! = du,, ν Normalize the length and velocity near the wall! u + = u v * y + = y v* v Called wall variables!
For parallel flow! Near the wall v * = 0! d dp < u'v' >= dx + d µ du the fluid knows nothing about what drives it. = µ du Thus we ignore Integrate from the wall to y:! the pressure u + = u v * y d < u'v' >= d µ du gradient! y + = y v* 0 v Resulting in:! < u'v' >= µ du < u'v' >= µ du Very close to the wall:! < u'v' > 0 so approximately! µ du = Integrating! u(y) = µ y Using the nondimensional values! ( ) 2 u v = y * µ v = v * * µ y v = v* y * ν or:! u + = y + Very close to the wall! v * = = µ du u + = u v * y + = y v* v Further away from the wall:! µ du v * = 0 so approximately! < u'v' >= = µ du du Taking:! < u'v' >= u + = u 2 du where! ν l v * T = l o and! 0 = κy yields! 2 du du < u'v' >= l o = du κy Giving:! κy du = 2 Or simply assume:! u' y du y + = y v* v We have! κy du = Using the nondimensional values! integrating! giving! du + κy + = v * =1 + du + = 1 κ + y + u + = 1 κ ln y + + C v * = = µ du u + = u v * y + = y v* v Thus, the velocity near the wall is! u + Viscous sub-layer! Velocity versus distance from wall! κ = 0.4 C = 5.5 1 u + = κ ln y + + C u + = y + 10! Buffer layer! Outer layer! y + v * = = µ du u + = u v * y + = y v* v For a practical engineering problem! L = 1m; U = 1m/s; ν = 10-6 (water)! The Reynolds number is therefore:! Re = LU v =106 For a flat plate, the average drag coefficient is! C D = 0.592Re 1/ 5 where! C D = F D 1 U 2 2 LW
Or! The average shear stress is therefore! And we find! C D = 0.0037 F D = LW = C 1 U 2 D 2 = 3.74 v * = 0.06 The average thickness of the viscous sub-layer is 10 in units of y + :! Thickness of the viscous sub-layer! y = 10ν 10 10 6 = =1.667 10 4 m = 0.1667 mm ν * 0.06 Find the thickness of the boundary layer! δ L = 0.37Re 1/ 5 δ L = 0.0233m = 23.3mm To resolve the viscous sublayer at the same time as the turbulent boundary layer would require a large number of grid points! To deal with this problem it is common to use wall functions where the mean velocity is matched with an analytical approximation to the viscosus sublayer.! For a reference, see: Patel, Rodi, and Scheuerer, Turbulence Models for Near-Wall and Low Reynolds Number Flows: A Review. AIAA Journal, 23, page 1308! Second order closure! Derive equations for the Reynolds stresses:! The k-epsilon and other two equation models have several serious limitations, including the inability to predict anisotropic Reynolds stress tensors, relaxation effects, and nonlocal effects due to turbulent diffusion.! For these problems it is necessary to model the evolution of the full Reynolds stress tensor! The Navier-Stokes equations in component form:! u i t + u u i j = - 1 p + ν 2 u i Multiply the equation by the velocity! u u i i t + u iu j = - 1 p + ν 2 u i and averaging leds to equations for! u i u j t
The new equations contain terms like! u i u i u j which are not known. These terms are therefore modeled! The Reynolds stress model introduces 6 new equations (instead of 2 for the k-e model. Although the models have considerably more physics build in and allow, for example, anisotrophy in the Reynolds stress tensor, these model have yet to be optimized to the point that they consistently give superior results.! Turbulence Modeling III! Grétar Tryggvason! Spring 2010! For practical problems, the k-e model or more recent improvements such as RNG are therefore most commonly used! Outline! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions! Wall bounded turbulence! Second order closure! Direct Numerical Simulations! Large-ed simulations! Summary! Direct Numerical Simulations! In direct numerical simulations the full unstea Navier-Stokes equations are solved on a sufficiently fine grid so that all length and time scales are fully resolved. The size of the problem is therefore very limited. The goal of such simulations is to provide both insight and quantitative data for turbulence modeling! Channel Flow! Wall! Flow direction! Streamwise velocity! Periodic streamwise and spanwise boundaries!
Streamwise vorticity! Channel Flow! Turbulent shear stress! Turbulent eddies generate a nearly uniform velocity profile! Streamwise vorticity! Turbulence are intrinsically linked to vorticity, yet laminar flows can also be vortical so looking at the vorticity is not sufficient to understand what is going on in a turbulent flows. Several attempts have been made to define properties of the turbulent flows that identifies vortices (as opposed to simply vortical flows.! One of the most successful method is the lambda-2 method of Hussain.! Visualizing turbulence! u v u = w u u y z v v y z w w y z ( ) = 1 2 S = 1 2 u + T u ( ) = 1 2 Ω = 1 2 u - T u 2 u v + u y w + u z 0 v u y w u z u y + v 2 v y w y + v z u y v 0 w y v z u z + w v z + w y 2 w z u z w v z w y 0 It can be shown that the second eigenvalue of! S 2 + Ω 2 define vortex structures! Referece: J. Jeong and F. Hussain, "On the identification of a vortex," Journal of Fluid Mechanics, Vol. 285, 69-94, 1995.! λ 2 λ 2 = 0.2 Other quantities have also been used, such as the second invariant of the velocity gradient:! Q = u i j u j i λ 2 = 0.3
Large Ed Simulations! Unstea simulations where the large scale motion is resolved but the small scale motion is modeled. Frequently simple models are used for the small scale motion. Most recently some success has been achieved by intrinsic large ed simulations where no modeling is used but monotonicity is enforced by the methods described in the lectures on hyperbolic methods! In the simplest case, the Smagorinsky ed viscosity is used in simulation of unstea flow, thus resulting in a viscosity that depends of the flow.! = l 2 0 ( 2S ij S ij ) 1/ 2 Summary! S ij = 1 U i + U j 2 j i Since the viscosity increases, the size of the smallest flow scales increases and lower resolution is needed! Turbulence models are used to allow us to simulate only the averaged motion, not the unstea small scale motion.! Turbulence modeling rest on the assumption that the small scale motion is universal and can be described in terms of the large scale motion.! Although considerable progress has been made, much is still not known and results from calculations using such models have to be interpreted by care!! For more information:! D. C. Wilcox, Turbulence Modeling for CFD (2 nd ed. 1998).! The author is one of the inventors of the k-ω model and the book promotes it use. The discussion is, however, general and very accessible, as well as focused on the use of turbulence modeling for practical applications in CFD!