More Related Content
Similar to 8015082 (20)
Recently uploaded (20)
8015082
- 1. 1 Lecture 8 - Turbulence Applied Computational Fluid Dynamics Instructor: Andr辿 Bakker 息 Andr辿 Bakker (2002-2005) 息 Fluent Inc. (2002)
- 2. 2 Turbulence What is turbulence? Effect of turbulence on Navier-Stokes equations. Reynolds averaging. Reynolds stresses.
- 3. 3 Instability All flows become unstable above a certain Reynolds number. At low Reynolds numbers flows are laminar. For high Reynolds numbers flows are turbulent. The transition occurs anywhere between 2000 and 1E6, depending on the flow. For laminar flow problems, flows can be solved using the conservation equations developed previously. For turbulent flows, the computational effort involved in solving those for all time and length scales is prohibitive. An engineering approach to calculate time-averaged flow fields for turbulent flows will be developed.
- 4. 4 Time What is turbulence? Unsteady, aperiodic motion in which all three velocity components fluctuate, mixing matter, momentum, and energy. Decompose velocity into mean and fluctuating parts: Ui(t) Ui + ui(t). Similar fluctuations for pressure, temperature, and species concentration values.
- 5. 5 Examples of simple turbulent flows jet mixing layer wake Some examples of simple turbulent flows are a jet entering a domain with stagnant fluid, a mixing layer, and the wake behind objects such as cylinders. Such flows are often used as test cases to validate the ability of computational fluid dynamics software to accurately predict fluid flows.
- 6. 6 Transition The photographs show the flow in a boundary layer. Below Recrit the flow is laminar and adjacent fluid layers slide past each other in an orderly fashion. The flow is stable. Viscous effects lead to small disturbances being dissipated. Above the transition point Recrit small disturbances in the flow start to grow. A complicated series of events takes place that eventually leads to the flow becoming fully turbulent.
- 7. 7 Transition in boundary layer flow over flat plate
- 8. 8 Transition in boundary layer flow over flat plate Turbulent spots Fully turbulent flowT-S waves
- 9. 9 Turbulent boundary layer Merging of turbulent spots and transition to turbulence in a natural flat plate boundary layer. Top view Side view
- 10. 10 Turbulent boundary layer Close-up view of the turbulent boundary layer.
- 11. 11 Transition in a channel flow Instability and turbulence is also seen in internal flows such as channels and ducts. The Reynolds number is constant throughout the pipe and is a function of flow rate, fluid properties and diameter. Three flow regimes are shown: Re < 2200 with laminar flow. Re = 2200 with a flow that alternates between turbulent and laminar. This is called transitional flow. Re > 2200 with fully turbulent flow.
- 12. 12 Transition in a jet flow
- 13. 13 Large StructureSmall Structure Large-scale vs. small-scale structure
- 14. 14 Weddell Sea off Antarctica
- 15. 15 Alaska's Aleutian Islands As air flows over and around objects in its path, spiraling eddies, known as Von Karman vortices, may form. The vortices in this image were created when prevailing winds sweeping east across the northern Pacific Ocean encountered Alaska's Aleutian Islands
- 16. 16 Alexander Selkirk Island in the southern Pacific Ocean20 km
- 17. 17
- 18. 18 Smoke ring A smoke ring (green) impinges on a plate where it interacts with the slow moving smoke in the boundary layer (pink). The vortex ring stretches and new rings form. The size of the vortex structures decreases over time.
- 19. 19 Simulation species mixing Arthur Shaw and Robert Haehnel, 2004
- 20. 20 Homogeneous, decaying, grid-generated turbulence Turbulence is generated at the grid as a result of high stresses in the immediate vicinity of the grid. The turbulence is made visible by injecting smoke into the flow at the grid. The eddies are visible because they contain the smoke. Beyond this point, there is no source of turbulence as the flow is uniform. The flow is dominated by convection and dissipation. For homogeneous decaying turbulence, the turbulent kinetic energy decreases with distance from grid as x-1 and the turbulent eddies grows in size as x1/2 .
- 21. 21 Re = 9.6 Re = 13.1 Re = 30.2 Re = 2000 Re = 26 Re = 10,000 Flow transitions around a cylinder For flow around a cylinder, the flow starts separating at Re = 5. For Re below 30, the flow is stable. Oscillations appear for higher Re. The separation point moves upstream, increasing drag up to Re = 2000.
- 22. 22 Turbulence: high Reynolds numbers Turbulent flows always occur at high Reynolds numbers. They are caused by the complex interaction between the viscous terms and the inertia terms in the momentum equations. Laminar, low Reynolds number free stream flow Turbulent, high Reynolds number jet
- 23. 23 Turbulent flows are chaotic One characteristic of turbulent flows is their irregularity or randomness. A full deterministic approach is very difficult. Turbulent flows are usually described statistically. Turbulent flows are always chaotic. But not all chaotic flows are turbulent.
- 24. 24 Turbulence: diffusivity The diffusivity of turbulence causes rapid mixing and increased rates of momentum, heat, and mass transfer. A flow that looks random but does not exhibit the spreading of velocity fluctuations through the surrounding fluid is not turbulent. If a flow is chaotic, but not diffusive, it is not turbulent.
- 25. 25 Turbulence: dissipation Turbulent flows are dissipative. Kinetic energy gets converted into heat due to viscous shear stresses. Turbulent flows die out quickly when no energy is supplied. Random motions that have insignificant viscous losses, such as random sound waves, are not turbulent.
- 26. 26 Turbulence: rotation and vorticity Turbulent flows are rotational; that is, they have non-zero vorticity. Mechanisms such as the stretching of three-dimensional vortices play a key role in turbulence. Vortices
- 27. 27 What is turbulence? Turbulent flows have the following characteristics: One characteristic of turbulent flows is their irregularity or randomness. A full deterministic approach is very difficult. Turbulent flows are usually described statistically. Turbulent flows are always chaotic. But not all chaotic flows are turbulent. Waves in the ocean, for example, can be chaotic but are not necessarily turbulent. The diffusivity of turbulence causes rapid mixing and increased rates of momentum, heat, and mass transfer. A flow that looks random but does not exhibit the spreading of velocity fluctuations through the surrounding fluid is not turbulent. If a flow is chaotic, but not diffusive, it is not turbulent. The trail left behind a jet plane that seems chaotic, but does not diffuse for miles is then not turbulent. Turbulent flows always occur at high Reynolds numbers. They are caused by the complex interaction between the viscous terms and the inertia terms in the momentum equations. Turbulent flows are rotational; that is, they have non-zero vorticity. Mechanisms such as the stretching of three-dimensional vortices play a key role in turbulence.
- 28. 28 What is turbulence? - Continued Turbulent flows are dissipative. Kinetic energy gets converted into heat due to viscous shear stresses. Turbulent flows die out quickly when no energy is supplied. Random motions that have insignificant viscous losses, such as random sound waves, are not turbulent. Turbulence is a continuum phenomenon. Even the smallest eddies are significantly larger than the molecular scales. Turbulence is therefore governed by the equations of fluid mechanics. Turbulent flows are flows. Turbulence is a feature of fluid flow, not of the fluid. When the Reynolds number is high enough, most of the dynamics of turbulence are the same whether the fluid is an actual fluid or a gas. Most of the dynamics are then independent of the properties of the fluid.
- 29. 29 Kolmogorov energy spectrum Energy cascade, from large scale to small scale. E is energy contained in eddies of wavelength 了. Length scales: Largest eddies. Integral length scale (k3/2 /竜). Length scales at which turbulence is isotropic. Taylor microscale (15僚u2 /竜)1/2. Smallest eddies. Kolmogorov length scale (僚3 /竜)1/4 . These eddies have a velocity scale (僚.竜)1/4 and a time scale (僚/竜)1/2 . 2 3 2 2 2 is the energy dissipation rate (m /s ) is the turbulent kinetic energy (m /s ) is the kinematic viscosity (m /s) k 竜 僚 Integral scale Taylor scale Kolmogorov scale Wavenumber Log E
- 30. 30 Vorticity and vortex stretching Existence of eddies implies rotation or vorticity. Vorticity concentrated along contorted vortex lines or bundles. As end points of a vortex line move randomly further apart the vortex line increases in length but decreases in diameter. Vorticity increases because angular momentum is nearly conserved. Kinetic energy increases at rate equivalent to the work done by large-scale motion that stretches the bundle. Viscous dissipation in the smallest eddies converts kinetic energy into thermal energy. Vortex-stretching cascade process maintains the turbulence and dissipation is approximately equal to the rate of production of turbulent kinetic energy. Typically energy gets transferred from the large eddies to the smaller eddies. However, sometimes smaller eddies can interact with each other and transfer energy to the (i.e. form) larger eddies, a process known as backscatter.
- 31. 31 t2 t3 t4 t5 t6 t1 (Baldyga and Bourne, 1984) Vortex stretching
- 32. 32 External flows: Internal flows: Natural convection: 5 105xRe along a surface around an obstacle where 袖 UL ReL where Other factors such as free- stream turbulence, surface conditions, and disturbances may cause earlier transition to turbulent flow. L = x, D, Dh, etc. 108 1010 Ra 袖留 硫 3 TLg Ra Is the flow turbulent?
- 33. 33 Turbulence modeling objective The objective of turbulence modeling is to develop equations that will predict the time averaged velocity, pressure, and temperature fields without calculating the complete turbulent flow pattern as a function of time. This saves us a lot of work! Most of the time it is all we need to know. We may also calculate other statistical properties, such as RMS values. Important to understand: the time averaged flow pattern is a statistical property of the flow. It is not an existing flow pattern! It does not usually satisfy the steady Navier-Stokes equations! The flow never actually looks that way!!
- 34. 34 Experimental Snapshot Example: flow around a cylinder at Re=1E4 The figures show: An experimental snapshot. Streamlines for time averaged flow field. Note the difference between the time averaged and the instantaneous flow field. Effective viscosity used to predict time averaged flow field. Effective ViscosityStreamlines
- 36. 36 Velocity decomposition Velocity and pressure decomposition: Turbulent kinetic energy k (per unit mass) is defined as: Continuity equation: Next step, time average the momentum equation. This results in the Reynolds equations. ' ' pPp += += :Pressure :Velocity uUu 0 0:;0 = === U Uuu div divdivaverageTimediv :flowmeantheforequationcontinuity ( ) ref i U k T wvuk 2/1 3 2 222 )( ''' 2 1 = ++= :intensityTurbulence
- 37. 37 錚 錚 錚 錚 錚 錚 + ++ =+ z wu y vu x u SUgraddiv x P Udiv t U momentumx Mx )''()''()'( )()( )( : 2 袖 U 錚 錚 錚 錚 錚 錚 + ++ =+ z wv y v x vu SVgraddiv y P Vdiv t V momentumy My )''()'()''( )()( )( : 2 袖 U 錚 錚 錚 錚 錚 錚 + ++ =+ z w y wv x wu SWgraddiv z P Wdiv t W momentumz Mz )'()''()''( )()( )( : 2 袖 U Turbulent flow - Reynolds equations
- 38. 38 These equations contain an additional stress tensor. These are called the Reynolds stresses. In turbulent flow, the Reynolds stresses are usually large compared to the viscous stresses. The normal stresses are always non-zero because they contain squared velocity fluctuations. The shear stresses would be zero if the fluctuations were statistically independent. However, they are correlated (amongst other reasons because of continuity) and the shear stresses are therefore usually also non-zero. Reynolds stresses 2 2 2 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' xx xy xz yx yy yz zx zy zz u u v u w u v v v w u w v w w 錚 錚金 錚 錚 錚 歎 錚 歎 錚 歎= = 錚 歎 錚 歎錚 歎 錚 歎 錚 錚 錚 錚
- 39. 39 Continuity: Scalar transport equation: Notes on density: Here is the mean density. This form of the equations is suitable for flows where changes in the mean density are important, but the effect of density fluctuations on the mean flow is negligible. For flows with Ti<5% this is up to Mach 5 and with Ti<20% this is valid up to around Mach 1. 0)( =+ U div t Turbulent flow - continuity and scalars 錚 錚 錚 錚 錚 錚 + +陸=陸+ 陸 陸陸 z w y v x u Sgraddivdiv t )''()''()''( )()( )( U
- 40. 40 Closure modeling The time averaged equations now contain six additional unknowns in the momentum equations. Additional unknowns have also been introduced in the scalar equation. Turbulent flows are usually quite complex, and there are no simple formulae for these additional terms. The main task of turbulence modeling is to develop computational procedures of sufficient accuracy and generality for engineers to be able to accurately predict the Reynolds stresses and the scalar transport terms. This will then allow for the computation of the time averaged flow and scalar fields without having to calculate the actual flow fields over long time periods.
Editor's Notes
- #15: The image dimensions are ~240km by ~320 km The image reveals two large ocean circulation features, called eddies, at the northernmost edge of the icepak in the Weddell Sear off Antarctica. These eddies play an important role in ocean circulation and the transportation of heat toward the pole. The ocean eddies are 40km to 60km in diameter. The dark areas are new ice and the lighter green areas are small ice floes that are swept along by the ocean currents. First year seasonal ice is shown in the darker green areas in the lower right corner. The open ocean to the north (top) is blue.
- #16: from: http://landsat.gsfc.nasa.gov/earthasart/vonkar.html Von Karman VorticesImage taken 7/4/2002 by Landsat 7 As air flows over and around objects in its path, spiraling eddies, known as Von Karman vortices, may form. The vortices in this image were created when prevailing winds sweeping east across the northern Pacific Ocean encountered Alaska&apos;s Aleutian Islands.
- #17: http://landsat.gsfc.nasa.gov/earthasart/vortices.html Karman VorticesImage taken 9/15/1999 Each of these swirling clouds is a result of a meteorological phenomenon known as a Karman vortex. These vortices appeared over Alexander Selkirk Island in the southern Pacific Ocean. Rising precipitously from the surrounding waters, the island&apos;s highest point is nearly a mile (1.6 km) above sea level. As wind-driven clouds encounter this obstacle, they flow around it to form these large, spinning eddies.
- #18: http://landsat.gsfc.nasa.gov/earthasart/karman.html Karman Vortices near BroutonaImage taken 6/5/2000 These Karman vortices formed over the islands of Broutona, Chirpoy, and Brat Chirpoyev (&quot;Chirpoy&apos;s Brother&quot;), all part of the Kuril Island chain found between Russia&apos;s Kamchatka Peninsula and Japan.