Chapter4 Scaling-Turbulent - Similarity - Flow
Chapter4 Scaling-Turbulent - Similarity - Flow
September 8, 2010
∂U ∂U −1 dP ∂ 2U ∂ 0 0
U +V = +ν 2 − (u v ) (1)
∂x ∂y ρ dx ∂y ∂y
and the continuity equation,
∂U ∂V
+ =0 (2)
∂x ∂y
Boundary conditions are U = V = 0 at y = 0 and U → Ue as y →
∞, where Ue is the freestream velocity at the edge of the boundary layer.
Capital letters correspond to mean quantities and lower case letters with
a prime denote fluctuations. Several analytical attempts were made in the
past to find self-similar solutions for equations (1) and (2) for turbulent wall-
bounded flows with adverse pressure gradient. These solutions should result
in self-similarity, which means that flow equations should be independent
of the streamwise coordinate. A scaling approach searches for coordinate
transformation to collapse the data of mean velocity profiles, Reynolds shear
stress, Reynolds normal stresses, etc. By trying to collapse the data to a
1
2
single curve, we are trying are go from partial differential equations of (1)
and (2) to an ordinary differential equation. A flow usually needs two scales
to be described universally,i.e. velocity and length. In unsteady flows, a time
scale may also be necessary. Apart from length and velocity, scalings have
also been sought for statistics like correlations, spectra etc. Despite several
attempts so far, we still do not have a universal scaling for turbulent flows
which can represent all of the known experimental data.
Ue − U y −u0 v 0 y
= f ( ), 2
= g12 ( )
uτ ∆ uτ ∆
u0 2 y v02 y
2
= g11 ( ), 2
= g22 ( )
uτ ∆ uτ ∆
Here ∆ = (δ ∗ Ue /uτ ) whereas f , g12 , g11 and g22 represent the functional
dependence of velocity defect, Reynolds shear, streamwise and wall-normal
stresses respectively. For velocity and Reynolds stresses, uτ has classically
been used as internal scale and Ue as external scale because all wall-bounded
flows have distinct inner and outer regions. Here uτ is the friction velocity
3
p
defined as τw /ρ and Ue is the local freestream flow velocity. Near sepa-
ration, uτ becomes zero and thus the wall-unit scaling becomes difficult to
implement. In most cases, the inner velocity scaling based on uτ was unable
to collapse the velocity and Reynolds stress profiles in the outer layer and
thus an effect of Reynolds number was evident in the normalised profiles.
For length, the classical inner scaling has been uτ /ν and the bound-
ary layer thickness δ as outer scale. Owing to the difficulty of precisely
measuring boundary layer thickness, especially over curved surfaces, many
researchers have instead preferred to use boundary layer momentum thick-
ness θ or displacement thickness δ ∗ as length scales because these are in-
tegral quantites and relatively easy to compute. Clauser (1954b) originally
proposed ∆ = (δ ∗ Ue /uτ ) as the length scale but it was largely replaced
by boundary layer thickness δ. The usual length scales for fully-developed
channel, pipe flows and bumps/humps/ramps are the half-channel height H,
diameter D and maximum height h respectively. According to Pope (2000),
a flow is fully-developed when the velocity statistics no longer change with
the streamwise coordinate.
Castillo & George (2001) reconsidered the similarity analysis of the tur-
bulent boundary layer and suggested outer velocity scales as Ue , Ue2 and
Ue2 dδ/dx for the velocity defect, Reynolds normal stress and Reynolds shear
stress respectively. Thus their similarity analysis took the following form:
Ue − U y −u0 v 0 y
= f ( ), dδ
= g 12 ( ),
Ue δ Ue2 ( dx ) δ
u0 2 y v02 y
2
= g11 ( ), 2
= g22 ( )
Ue δ Ue δ
The validity of Ue , Ue2 and Ue2 dδ/dx as outer velocity scales was put in
question by Maciel et al. (2006a) for flows with strong streamwise pressure
gradients. They pointed out that Reynolds stresses scaled with Ue2 can grow
by an order of magnitude for a flow experiencing a strong favorable pressure
gradient and then a strong adverse pressure gradient like on the suction side
of an airfoil.
Zagarola & Smits (1998) proposed Uc − Ub as the outer velocity scale for
turbulent pipe flow, where Uc is the centreline velocity and Ub is the bulk
mean velocity. They showed that Uc − Ub is better than uτ for collapsing
the velocity defect profiles. For turbulent boundary layers, they proposed
an equivalent scale Uzs = Ue δ ∗ /δ, which proved to be a valid scale for the
outer layer. Their scale not only removed Reynolds number effects but also
4
the upstream condition effects were noted to disappear. Later Maciel et al.
2
(2006b) showed that Uzs is a valid scale for Reynolds stresses. To complement
the Zagarola & Smits (1998) velocity scale Uzs , Maciel et al. (2006a) replaced
uτ by Uzs to get ∆zs = Ue δ ∗ /Uzs , which is the same as boundary layer
thickness δ.
DeGraaff & Eaton (2000) proposed a mixed scaling uτ Ue for the Reynolds
stresses but scaling was found not to be universal in the outer region of the
boundary layer. They did notice however the disappearance of Reynolds
number effects with the mixed scaling.
Stanislas et al. (2008) carried out stereo PIV experiments of flat plate
turbulent boundary layer over a range of Reynolds number Reθ from 7800
to 15000. They argued that if Kolmogorov scales determines the friction
at the wall and if they are more representative of dissipation than , then
they should replace the viscosity ν as a scale. The viscosity shall then be-
come a kind of ‘turbulence thermodynamic’ quantity by contributing towards
Kolmogorov scale (η). They found that very near to the wall, the ratio of
Kolmogorov scales to the integral scales is of order 1 and thus near the wall
the Kolmogorov scale η should be used as a length scale. They plotted the
mean velocity U/Ue as a function of y/η and obtained a good collapse near
the wall. However, a good collapse was not observed in the outer layer and
to use δ as the outer length scale, they defined a new length scale η ∗ in form:
∗ δ
η = min η, . (3)
A
Here A is of the order of 1000. The function was designed such that δ
shall take over as a length scale if Kolmogorov scale η is greater than δ/1000.
Therefore, as summarised above, the existence of a plethora of velocity,
length and stress scales proves that we are still far from having a single
universally accepted set of scales for turbulent wall-bounded flows.
for collapsing the mean velocity profiles based on a “pressure” velocity (see
??):
ν dP 1/3
up = ( ) . (4)
ρ dx
Mellor & Gibson (1966) proposed another pressure velocity scale up for
flows near separation,
s
δ ∗ dp
up = uτ β 1/2 = . (5)
ρ dx
Perry & Shofield (1973) proposed a different velocity scale us , which is
given as,
3/2 s
us umτ δ ∗
= 13.53 (6)
Ue Ue ym
where umτ = (τmax /ρ)1/2 and ym is y at τmax . They claimed that all
velocity scales based upon the local pressure gradient are not appropriate
and instead proposed a velocity scale us which explicitly depends on the
maximum shear-stress. It was determined by a fit to the velocity profile in
the same manner as uτ was determined by a fit to the Clauser plot. When
u0 v 0 + max was greater than or equal to 1.5, then uτ was supposed to be replaced
by us . They also proposed a new length scale ∆ps = 2.86δ ∗ Ue /us . Angele &
Muhammad-Klingmann (2005) found that up and us scaled their data both
before and after separation.
Elsberry et al. (2000) performed experimental studies on a turbulent
boundary layer on the verge of separation (thus undergoing strong adverse
pressure gradient) and showed that the traditional outer layer scaling Ue
did not collapse the stress profiles, although it did collapse the mean veloc-
ity profile. They proposed Ue Uemax as an outer layer scaling for Reynolds
shear stress and Ue 2max for the Reynolds normal stresses, where Uemax is the
maximum freestream velocity,
u02 y v 02 y u0 v 0 y
= g11 ( 0.2
), = g22 ( 0.2
), = g12 ( ).
Ue2 max θReθ 2
Ue max θReθ Ue Uemax θReθ0.2
This scaling was thought to account for the non-equilibrium nature of the
outer layer. They also introduced a new length scale θRe0.2θ , based on the
momentum thickness θ and Reynolds number Reθ = (Ue θ/ν). The new set
of scalings were found to achieve a good collapse of the Reynolds stresses.
6
Ue − U y −u0 v 0 y
= f ( ), 2
= g12 ( )
Uzs ∆ Uzs ∆
u0 2 y v02 y
2
= g11 ( ), 2
= g22 ( )
Uzs ∆ Uzs ∆
locity fluctuations start to collapse in the outer layer at y/δ > 0.2, although
for some data-sets, the collapse is not as good as the mean velocity profiles.
The flows like Skåre & Krogstad (1994), which are in equilibrium (Clauser-
Townsend sense) show self-similarity with the new scale too. But self-similarity
in mean velocity have also been showed by those flows (Aubertine & Eaton
(2005) and Spalart & Watmuff (1993)) which are not in classical equilibrium
(Clauser-Townsend sense). The collapse of the rms velocity fluctuations is
not as good as traditional equilibrium flows but is better than other tradi-
tional velocity scales like uτ , Ue or mixed scalings uτ Ue . This corroborates
the proposition of Maciel et al Maciel et al. (2006a), “a flow can be in equilib-
rium even if a particular similarity theory would define it as non-equilibrium
flow. Conversely, if a flow is found to be in equilibrium according to a given
similarity analysis then it is in equilibrium regardless of the similarity as-
sumptions”.
10
0.9
x= 2.00’
x= 2.75’
0.8 x= 3.50’
x= 4.00’
x= 4.50’
0.7 x= 4.75’
x= 4.92’
0.6
(Ue-u)/Ueβ0.1
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
0.8
(Ue-u)/Ueβ0.1
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
0.9
x= 2.0’
x= 4.0’
0.8 x= 5.5’
x= 7.0’
0.7
0.6
(Ue-u)/Ueβ0.1
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
Figure 4: Velocity defect profiles scaled by Ue β 0.1 for Bradshaw & Ferriss
(1965).
1.2
x= 2.0’
x= 2.5’
x= 3.0’
x= 3.5’
1 x= 4.0’
x= 5.0’
0.8
(Ue-u)/Ueβ0.1
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
0.8
x= 4.0
x= 4.2
x= 4.4
0.7 x= 4.6
x= 4.8
0.6
0.5
(Ue-u)/Ueβ0.1
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
y/δ
Figure 6: Velocity defect profiles scaled by Ue β 0.1 for Skåre & Krogstad
(1994).
1.6
x=0.60
x=0.65
x=0.70
1.4 x=0.75
x=0.80
x=0.85
1.2
1
(U∞-u)/Ueβ0.1
0.8
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
Figure 7: Velocity defect profiles scaled by Ue β 0.1 for Spalart & Watmuff
(1993).
13
0.7
x’= 0.00
x’= 0.25
x’= 0.33
0.6 x’= 0.50
x’= 0.75
0.5
(U∞-u)/Ueβ0.1
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
y/δ
Figure 8: Velocity defect profiles scaled by Ue β 0.1 for Aubertine & Eaton
(2005).
1
x=523
x=723
0.9 x=925
x=1121
0.8
0.7
0.6
(Ue-u)/Ueβ0.1
0.5
0.4
0.3
0.2
0.1
-0.1
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
y/δ
Figure 9: Velocity defect profiles scaled by Ue β 0.1 for Nagano et al. (1993).
14
0.11
x’= 0.25
x’= 0.33
x’= 0.50
0.1 x’= 0.75
0.09
0.08
u’/Ueβ0.1
0.07
0.06
0.05
0.04
0.03
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
Figure 10: Streamwise rms velocity fluctuations scaled by Ue β 0.1 for Auber-
tine & Eaton (2005).
0.2
x=0.60
x=0.65
0.18 x=0.70
x=0.75
x=0.80
0.16 x=0.85
x=0.90
0.14 x=0.95
x=1.00
0.12
u’/Ueβ0.1
0.1
0.08
0.06
0.04
0.02
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
Figure 11: Streamwise rms velocity fluctuations scaled by Ue β 0.1 for Spalart
& Watmuff (1993).
15
0.1
x= 4.0
x= 4.2
0.09 x= 4.4
x= 4.6
x= 4.8
0.08 x= 5.0
x= 5.2
0.07
0.06
u’/Ueβ0.1
0.05
0.04
0.03
0.02
0.01
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
y/δ
Figure 12: Streamwise rms velocity fluctuations scaled by Ue β 0.1 for Skåre &
Krogstad (1994).
0.12
x=523
x=723
x=925
x=1121
0.1
0.08
u’/Ueβ0.1
0.06
0.04
0.02
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
y/δ
Figure 13: Streamwise rms velocity fluctuations scaled by Ue β 0.1 for Nagano
et al. (1993).
16
The boundary layer flows that were defined as being non-equilibrium ac-
cording to the classical similarity theory show self-similarity with the new
scaling. This implies that other formulations of similarity analysis and hence
equilibrium exist and a flow can be in equilibrium even if a particular simi-
larity theory does not apply. For Spalart & Watmuff (1993) and Aubertine &
Eaton (2005), which are non-equilibium (Clauser-Townsend sense) flows, the
velocity defect profiles show a good collpase and rms velocity fluctuations
show a collapse of profiles in the outer region (y/δ = 0.2). On the other
hand, Skåre & Krogstad (1994) is an equilibrium turbulent boundary layer
and self-similarity in rms velocity fluctuations is achieved along the whole
boundary layer thickness although the collapse is better in the outer region.
The scale seems to be valid for the outer region (y/δ > 0.2) of the turbulent
boundary layers with adverse pressure gradient developing over flat surfaces
(without curvature). As the scale involves the wall friction velocity uτ in its
formulation, it can not be used for the separated turbulent boundary layers.
17
80
x= 6.92’
x= 11.00’
x= 12.75’
70 x= 18.58’
x= 23.83’
x= 26.92’
60 x= 29.75’
x= 32.25’
50
(Ue-u)/ui
40
30
20
10
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
y/δ
Figure 14: Velocity defect profiles scaled by uτ (δ ∗ /δ)1/2 for Clauser (1954b).
90
x= 7.50’
x= 9.00’
80 x= 11.00’
x= 12.67’
x= 16.17’
70 x= 19.17’
60
50
(Ue-u)/ui
40
30
20
10
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
y/δ
Figure 15: Velocity defect profiles scaled by uτ (δ ∗ /δ)1/2 for Clauser (1954a).
18
80
x= 2.0’
x= 4.0’
x= 5.5’
70 x= 7.0’
60
50
(Ue-u)/ui
40
30
20
10
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
y/δ
Figure 16: Velocity defect profiles scaled by uτ (δ ∗ /δ)1/2 for Bradshaw (1966).
80
x= 2.0’
x= 4.0’
x= 5.5’
70 x= 7.0’
60
50
(Ue-u)/ui
40
30
20
10
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
Figure 17: Velocity defect profiles scaled by uτ (δ ∗ /δ)1/2 for Bradshaw &
Ferriss (1965).
19
140
x=0.65
x=0.70
x=0.75
120 x=0.80
x=0.85
100
80
(Ue-u)/ui
60
40
20
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
Figure 18: Velocity defect profiles scaled by uτ (δ ∗ /δ)1/2 for Spalart & Wat-
muff (1993).
45
x’= 0.25
x’= 0.33
40 x’= 0.50
x’= 0.75
35
30
(Ue-u)/ui
25
20
15
10
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
y/δ
Figure 19: Velocity defect profiles scaled by uτ (δ ∗ /δ)1/2 for Aubertine &
Eaton (2005).
20
35
x= 4.0
x= 4.2
x= 4.4
30 x= 4.6
x= 4.8
x= 5.0
25 x= 5.2
20
(Ue-u)/ui
15
10
-5
0 0.2 0.4 0.6 0.8 1 1.2 1.4
y/δ
Figure 20: Velocity defect profiles scaled by uτ (δ ∗ /δ)1/2 for Skåre & Krogstad
(1994).
20
x=0.60
x=0.65
18 x=0.70
x=0.75
x=0.80
16 x=0.85
x=0.90
14 x=0.95
x=1.00
12
u’/ui
10
0
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
8
x’= 0.00
x’= 0.25
x’= 0.33
x’= 0.50
7 x’= 0.75
6
u’/ui
2
0 0.2 0.4 0.6 0.8 1 1.2
y/δ
8
u’/ui
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
y/δ
4 Conclusions
New empirically-determined velocity scales Ue β 0.1 and uτ (δ ∗ /δ)1/2 are pro-
posed for the velocity defect (Ue − u), rms velocity fluctuations u0rms and
Reynolds stresses ui uj of the attached turbulent boundary layers undergo-
ing adverse pressure gradient. The validity and universality of the proposed
scale is tested by applying it to seven (numerical and experimental) data-
sets with different experimental techniques, Reynolds number ranges, flow
geometries and pressure gradients. The proposed scales achieve a better col-
lapse of the velocity defect and streamwise rms velocity fluctuation profiles
as compared to the traditional scalings of uτ and Ue in the outer region of
the attached turbulent boundary layers (y/δ > 0.2) undergoing adverse pres-
sure gradient. As expected, the equilibrium turbulent boundary layers show
self-similarity with the new set of scalings too. The self-similarity shown by
classicaly non-equilibrium turbulent boundary layers implies that a flow may
be in equilibrium even if a particular similarity theory does not apply.
References
Angele, K. & Muhammad-Klingmann, B., 2005 Effect of streamwise
vortices on the turbulence structure of a separating boundary layer. Eur.J.
Mech. B/Fluids 24, 539–554.
Castillo, L. & George, W., 2001 Similarity analysis for turbulent bound-
ary layers with pressure gradient: Outer flow. AIAA J. 39, 41–47.
Maciel, Y., Rossignol, K.-S. & Lemay, J., 2006a Similarity in the
outer region of adverse-pressure-gradient turbulent boundary layers. AIAA
Journal 44 (11), 2450–2464.
Maciel, Y., Rossignol, K.-S. & Lemay, J., 2006b A study of a turbu-
lent boundary layer in stalled-airfoil-type flow conditions. Experiments in
Fluids 41, 573–590.
Nagano, Y., Tagawa, M. & Tsuji, T., 1993 Effects of adverse pressure
gradient on mean flows and turbulence statistics in a boundary layer. In
Turbulent Shear Flows 8. F. Durst et al.(eds.),. Springer-Verlag, Berlin,
pages 7–21.
Newman, B., 1951 Some contributions to the study of the turbulent bound-
ary layer near separation. Technical report, Aust. Departement Supply
Rept.ACA-53.
Perry, A. & Shofield, W., 1973 Mean velocity and shear stress distribu-
tion in turbulent boundary layers. Phys. Fluids 113, 2068–2074.
Song, S. & Eaton, J., 2004 Reynolds number effects on a turbulent bound-
ary layer with separation, reattachment, and recovery. Exp. Fluids. 36,
246–258.