Problems and Solutions
Problems and Solutions
Problems and Solutions
Instantaneous Mean
The instantaneous value of any flow variable can be decomposed into mean + fluctuation.
v
u
is decomposed into
mean + fluctuation
The first is useful in deriving theoretical results but becomes cumbersome in general use. The
notation being used is, hopefully, obvious from the context.
In experimental work and in steady flow the “mean” is usually a time mean, whilst in
theoretical work it is the probabilistic (or “ensemble”) mean. The process of taking the mean
of a turbulent quantity or a product of turbulent quantities is called Reynolds averaging.
Engineers are usually only interested in the mean flow. However, turbulence must still be
considered because, although the averages of individual fluctuations (e.g. u ′ or v ′ ) are zero,
the average of a product (e.g. u ′v ′ ) is not and may lead to a significant net flux.
Consider mass and momentum fluxes in the y direction across surface A. For simplicity,
assume constant density.
7.4.1 Continuity
ρ vA
Mass flux: vA
Average mass flux: vA
The only change is that the instantaneous velocity is replaced by the mean velocity.
The mean velocity satisfies the same continuity equation as the instantaneous velocity.
7.4.2 Momentum
ρ vA
x-momentum flux: ( vA)u = (uv) A
Average x-momentum flux: ( vA)u = (u v + u ′v ′) A u
The average momentum flux has the same form as the instantaneous momentum flux except
for additional fluxes u ′v ′A due to the net effect of turbulent fluctuations. These additional
terms arise because of the averaging of a product of fluctuating quantities.
A net rate of transport of momentum u ′v ′A from lower to upper side of an interface ...
• is equivalent to a net rate of transport of momentum − u ′v ′A from upper to lower;
• has the same dynamic effect (i.e. same rate of transfer of momentum) as a stress (i.e.
force per unit area) of − u ′v ′ .
This apparent stress is called a Reynolds stress. In a fully-turbulent flow it is usually much
larger than the viscous stress.
Other Reynolds stresses ( − u ′u ′ , − v ′v ′ , etc.) emerge when considering the flux of the
different momentum components in different directions.
The mean velocity satisfies the same momentum equation as the instantaneous velocity,
except for additional apparent stresses: the Reynolds stresses − u i′u ′j
The nature of the turbulent stress can be illustrated by considering the motion of particles
whose fluctuating velocities allow them to cross an interface.
If particle A migrates upward (v′ > 0) then it tends to retain its original y
momentum, which is now lower than its surrounds (u′ < 0).
B
If particle B migrates downward (v′ < 0) it tends to retain its original v'
momentum which is now higher than its surrounds (u′ > 0).
A
In both cases, − u ′v ′ is positive and, on average, tends to reduce the
momentum in the upper fluid or increase the momentum in the lower
fluid. Hence there is a net transfer of momentum from upper to lower U
fluid, equivalent to an additional mean stress.
Velocity Fluctuations
Normal stresses: u ′ 2 , v ′ 2 , w′ 2
Shear stresses: v ′w′, w′u ′, u ′v ′
(In slightly careless, but extremely common, usage both − u ′v ′ and u ′v ′ are referred to as
“stresses”.)
v
Most turbulent flows are anisotropic; i.e. u ′ 2 , v ′ 2 , w′ 2 are different. u
Turbulent kinetic energy: k = 12 (u ′ 2 + v ′ 2 + w′ 2 )
root-mean-square fluctuation u′ 2
k
Turbulence intensity: = = rms = 3
mean velocity U U
In general, the advection of any scalar quantity φ gives rise to an additional scalar flux in the
mean-flow equations; e.g.
vφ = v φ + v ′φ′ (2)
123
additional flux
Again, the extra term is the result of averaging a product of fluctuating quantities.
At high Reynolds numbers, turbulent fluctuations cause a much greater net momentum
transfer than viscous forces throughout most of the flow. Thus, accurate modelling of the
Reynolds stresses is vital.
PRODUCTION large
eddies
by mean flow
It is common experience that turbulence can be transported (i.e. carried) with the flow.
(Think of the turbulent wake behind a vehicle or downwind of a large building).
It can be proved mathematically (try the Turbulent Boundary Layer notes if you are really
desperate to know how!) that:
• The source term for an individual Reynolds stress u i u j transport equation has the form:
net source = production + redistribution − dissipation
where:
– production Pij is determined by mean velocity gradients;
– redistribution ij transfers energy between stresses via pressure fluctuations;
– dissipation ij involves viscosity acting on fluctuating velocity gradients.
There are also “advection” terms (turbulence carried with the flow) and “diffusion” terms
(if turbulent stresses vary from point to point).
u2 v2 w2 REDISTRIBUTION
by pressure
fluctuations
DISSIPATION by viscosity
• The production terms for different Reynolds stresses involve different mean velocity
gradients; for example, the rate of production (per unit mass) of u1u1 = u 2 and
u1u 2 = uv are, respectively,
∂U ∂U ∂U
P11 = −2(uu + uv + uw )
∂x ∂y ∂z
(3)
∂V ∂V ∂V ∂U ∂U ∂U
P12 = −(uu + uv + uw ) − (vu + vv + vw )
∂x ∂y ∂z ∂x ∂y ∂z
(Exercise: by “pattern-matching” write production terms for the other stresses).
• Because:
(i) mean velocity gradients are greater in some directions than others,
(ii) motions in certain directions are selectively damped (e.g. by buoyancy forces or
rigid boundaries),
turbulence is usually anisotropic, i.e. u 2 , v 2 , w 2 are all different.
• In practice, most turbulence models do not actually solve transport equations for all
turbulent stresses, but only for the turbulent kinetic energy k = 12 (u 2 + v 2 + w 2 ) ,
relating the other stresses to this by an eddy-viscosity formula (see Section 8).
A flow for which there is only one non-zero mean velocity gradient, ∂U/∂y, is called a simple
shear flow. Because they form a good approximation to many real flows, have been
extensively researched in the laboratory and are amenable to basic theory they are an
important starting point for many turbulence models.
For such a flow, the first of (3) and similar expressions show that P11 > 0 but that
P22 = P33 = 0, and hence u 2 tends to be the largest of the normal stresses because it is the
only one with a non-zero production term. On the other hand, if there is a rigid boundary on
y = 0 then it will selectively damp wall-normal fluctuations; hence v 2 is the smallest of the
normal stresses.
If there are density gradients (for example in atmospheric or oceanic flows, in fires or near
heated surfaces) then buoyancy forces will either damp (stable density gradient) or enhance
(unstable density gradient) vertical fluctuations.
Mixing layer
u2
∆y
v 2
U ui u j
∆U
y y
Wake
(plane or axisymmetric)
∆y
uv u2
U ui u j
∆U
y y
Jet
(plane or axisymmetric)
∆y
u2
v2
uv
U ui uj
∆U
y y
constant-stress layer u2
-uv v2
U ui uj
Even though the overall Reynolds number Re = U ∞ L/ may be large, and hence viscous
transport much smaller than turbulent transport in the majority of the flow, there must be a
thin layer very close to the wall where the local Reynolds number based on distance from the
wall, Re y = u y/ , is small and hence molecular viscosity is important.
Wall Units
An important parameter is the wall shear stress w (drag per unit area). Like any other stress
this has dimensions of [density] × [velocity]2 and hence it is possible to define an important
velocity scale called the friction velocity uτ (also written u*):
u = w / (4)
From u and it is possible to form a viscous length scale l = /u. Hence, we may define
non-dimensional velocity and distance from the wall in so-called wall units:
U yu
U+ = , y+ = (5)
u
Viscous Sublayer
Very close to a smooth wall, turbulence is damped out by the presence of the boundary. In
this region the shear stress is predominantly viscous. Assuming constant shear stress,
∂U
w =
∂y
y
⇒ U= w (6)
i.e. The mean velocity profile in the viscous sublayer is linear. This is generally a good
approximation in the range y+ < 5.
Log-Law Region
At large Reynolds numbers, the turbulent part of the shear stress dominates throughout most
of the boundary layer so that on dimensional grounds, since u and y are the only possible
velocity and length scales,
∂U u
∝
∂y y
Integrating, and putting part of the constant of integration inside the logarithm (to make its
argument dimensionless):
1 u y
U = u ( ln + B) (7)
(von Kármán’s constant) and B are universal constants with experimentally-determined
values of about 0.41 and 5 respectively.
Using the definition of wall units (equation (5)) these velocity profiles are often written in
non-dimensional form:
(viscous sublayer) U + = y +
1 (8)
(log layer) U + = ln y + + B
Experimental measurements indicate that the log law actually holds to a good approximation
over a substantial proportion of the boundary layer. (This is where the logarithm originates in
common friction-factor formulae such as the Colebrook-White formula for pipe flow).
Consistency with the log law is probably the single most important consideration in the
construction of turbulence models.
• A means of specifying the Reynolds stresses (and hence solving the mean flow
equations) is called a turbulence model or turbulence closure.
Q2. The accepted critical Reynolds number in a round pipe (based on bulk velocity and
diameter) is 2300. At what speed is this attained in 5-cm-diameter pipe for (a) air; (b) water?
Q3. Sketch the mean velocity profile in a pipe at Reynolds numbers of (a) 500; (b) 50 000.
What is the shear stress along the pipe axis in either case?
Q4. Explain the process of flow separation. How does deliberately “tripping” a developing
boundary layer help to prevent or delay separation on a convex curved surface?
Q5. The following couplets are measured values of (u,v) in an idealised 2-d turbulent flow.
Calculate u , v , u ′2 , v ′2 , u ′v ′ from this set of numbers.
(3.6,0.2) (4.1,–0.4) (5.2,–0.2) (4.6,–0.4) (3.4,0.0)
(3.8,-0.4) (4.4,0.2) (3.9,0.4) (3.0,0.4) (4.4,–0.3)
(4.0,-0.1) (3.4,0.1) (4.6,-0.2) (3.6,0.4) (4.0,0.3)
Q6. The rate of production (per unit mass of fluid) of u 2 and uv are, respectively,
∂U ∂U ∂U
P11 = −2(uu + uv + uw )
∂x ∂y ∂z
∂V ∂V ∂V ∂U ∂U ∂U
P12 = −(uu + uv + uw ) − (vu + vv + vw )
∂x ∂y ∂z ∂x ∂y ∂z
(a) By inspection, write down similar expressions for P22, P33, P23, P31, the rates of
production of v 2 , w 2 , vw and wu respectively.
(b) Write down expressions for P11, P22, P33 and P12, P23, P31 in simple shear flow (where
∂U/∂y is the only non-zero mean velocity gradient). What does this indicate about the
relative distribution of turbulence energy amongst the various Reynolds-stress
components? Write down also an expression for P(k), the rate of production of
turbulence kinetic energy.