Nezu 1977博士论文
Nezu 1977博士论文
by Iehisa NEZU
July 1977
Department of C i v i l Engineering
Kyoto University
Kyoto 606, Japan
_ i
n
CONTENTS
Page
ABSTRACT
K
CHAPTER 1 INTRODUCTION
-ii-
CHAPTER 5 EXPERIMENTAL C O N S I D E R A T I O N ON
MECHANISM OF T U R B U L E N C E PRODUCTION
5.1 Introduction
5 2 Probability distributions of the velocity fluctuations
and t h e i n s t a n t a n e o u s R e y n o l d s s t r e s s
5.2.1 P r o b a b i l i t y d e n s i t y f u n c t i o n o f v e l o c i t y and
i t s h i g h e r moments
5.2.2 P r o b a b i l i t y d e n s i t y f u n c t i o n o f instantaneous
Reynolds s t r e s s
5.2.3 Conditional p r o b a b i l i t y d i s t r i b u t i o n o f
instantaneous Reynolds s t r e s s
5 3 I n t e r n a l s t r u c t u r e o f R e y n o l d s s t r e s s and
mechanism o f t u r b u l e n c e - p r o d u c t i o n
5 3.1 I n t e r n a l s t r u c t u r e o f Reynolds s t r e s s
5 3.2 R e l a t i v e i n t e n s i t y o f t h e d i f f e r e n t events
5.3.3 Mechanism o f t u r b u l e n c e - p r o d u c t i o n i n t h e w a l l
r e g i o n and i t s p r e d i c t i o n
5.3.1* E f f e c t o f w a l l roughness upon b u r s t i n g events
5'.3.5 Behaviour o f v e l o c i t y f l u c t u a t i o n s
when t u r b u l e n c e occurs
5.4 Periodic c h a r a c t e r i s t i c s of the bursting phenomenon
5.4.1 B r i e f r e c a p i t u l a t i o n of the bursting period
obtained p r e v i o u s l y
5.U.2 D i s c r i m i n a t i o n c r i t e r i o n and d e f i n i t i o n o f
bursting period .
5.4.3 B u r s t i n g p e r i o d and i t s p r o b a b i l i t y c h a r a c t e r i s t i c s
5.5 Internal structure of bursting phenomenon
5.5.1 Dependence o f t h e b u r s t i n g phenomenon upon t h e i n n e r
and outer parameters
5.5.2 F u r t h e r suggestions on b u r s t i n g phenomenon
5.6 Visual observation of bursting phenomenon
56 l Instantaneous v e l o c i t y p r o f i l e s i n v e r t i c a l plane
5'.6.2 Instantaneous v e l o c i t y p r o f i l e s i n h o r i z o n t a l plane
CHAPTER 6 PHYSICAL MODEL OF WALL TURBULENCE
6.1 Introduction
6.2 IT-eddy model
6.2.1 C o n s t i t u t i o n o f a II-eddy model and i t s f o r m u l a t i o n
6.2.2 Turbulence i n t e n s i t i e s and Reynolds stresses
evaluated by a II-eddy
6.2.3 S p e c t r a l d e n s i t y f u n c t i o n s o f energy
6.2.4 Comparison o f t u r b u l e n c e i n t e n s i t i e s w i t h
the experimental data
6.3 Renewal model
6.3.1 Formulation o f a renewal model
6.3-2 Further c o n s i d e r a t i o n o f renewal model by t a k i n g
account o f t h e p r o b a b i l i t y d i s t r i b u t i o n o f
bursting period
6.3.3 Turbulent c h a r a c t e r i s t i c s evaluated by t h e present
renewal model
6.3.4 Explanation o f t h e b u r s t i n g process by t h e renewal
model
6.4 Combined model
-iv-
PART 2 TURBULENT STRUCTURE OVER PERMEABLE BED
1 0 9
CHAPTER 7 INTRODUCTION
1 5 1
CONCLUDING REMARKS
-V -
ACKNOWLEDGEMENTS IS6
— over
-vi-
TuR BULENT STRUCTURE I N OPEN-CHANNEL EIOWS
Abstract
flows are still very insufficient, compared with those of air flows
- 1-
(5) Modelling of turbulent shear flow
Part 1 deals with the turbulent structure over smooth and rough
solid beds, and Part 2 deals with that over permeable porous bed .
average . Chapter 5 deals with (3) and (4). That is, the i
turbulent structure in terms of short- time average or the bursting
- 2-
we propose some new physical models of wall turbulence ·on a basis
s uction or injection through the bed , and intend to grope our way
by such a transpiration.
were given.
- 3-
PART 1 T u r b u l e n t s t r u c t u r e o v e r s o l i d bed
Chapter 1 Introduction
i n v e s t i g a t i o n s on s o l i d - w a l l t u r b u l e n c e i n a i r f l o w ( boundary l a y e r
R o t t a ( l 9 7 2 ) , Hinze(l975) and o t h e r s .
-k-
References (Chapter 1)
l ) Rouse, H. & Ince, S. : H i s t o r y o f h y d r a u l i c s , State Univ. o f Iowa
1957¬
2) Monin, A.S. & Yaglom, A.M. : Statistical fluid mechanics,
MIT Press, v o l . 1 , 1971.
3) Monin, A.S. & Yaglom, A.M. : Statistical fluid mechanics,
MIT Press, v o l . 2 , 1975-
-5-
C h a p t e r 2 T h e o r e t i c a l c o n s i d e r a t i o n on t u r b u l e n t s t r u c t u r e
2.1 Introduction
x
dx dy p dx dx dy dy
d
± { - l s + v JL)=- l (2.) 3
Oy dy
'
l
P= P j r U - y J c o s f l -fiC7 -7J) (2«4)
2
Since t h e w a l l shear s t r e s s x 0 i s d e f i n e d as r //>={-CITH- i> dul d )
0 y | _ C = ut
(U, i s t h e f r i c t i o n v e l o c i t y ) , t h e Reynolds s t r e s s d i s t r i b u t i o n
i s g i v e n by
-üv 1 dll*
where, ( = y/ , u\=u/u,
A and R,= wju>\.
-6-
where, i s t h e h u l k mean v e l o c i t y .
Defining E= *( au/dy) 1
( d i r e c t - d i s s i p a t i p n o f t h e mean f l o w ) , ps-S au/d y
or
(2 9>
r>+«.7**-£ -
+
where, y = U,y/y , and t h e ( + ) - s u f f i x denotes t h e dimensionless quan
t i t y by U # and v i n t h e f o l l o w i n g .
Now, by assuming t h e P r a n d t l ' s m i x i n g - l e n g t h model, t h a t i s :
-ms l
t (du/d y
y (2«io) i
, and i f R w i s very l a r g e ,
When y +
i s very s m a l l , o r , (2.12) y i e l d s a v e l o c i t y distri-
+ +
b u t i o n i n t h e viscous sublayer, i . e . t/ = y +
U <D (2-15)
-7-
F i g . 2.3 shows t h e f o l l o w i n g r e s u l t s , t h a t i s :
(2-16)
2.2.2 Budget o f t u r b u l e n t e n e r g y
The basic equations o f t u r b u l e n t components can be w r i t t e n :
9i dU dU
„ vv . — UkU
> , +
j T7oV i aT~\ / au,\ / a . \
u
* 3 j r
) T~ Ö
a _
T + — 1
J +2u( — : 1 f — i )
. *t " So:,/ % j J U j
p
= 0 (2-20)
djc.
I n t h e f l o w f i e l d o f F i g . 2 . 1 , (2.20) can be s i m p l i f i e d as f o l l o w s
dU_ p du _ 1 d l d#\
dy +
pl^~ t
^ U y V - " - V
T y ) (2-21)
^ a
y 2 dy \ p Q ) (2-22)
!
P dm _ l a/— g \ B
fi dz 2 dy\ y
dy
d •(2-23)
2
du.s
+( (2-24)
Oy I \ Q z
+ —
~pü düv "I
d (2-25)
" yJ
-8-
Denoting 9** " 2+
+ •* • • - « . + «. + ». ( 2
" 2 6 )
5, ±iH
5,1 2
-rZ* + - *a f \(2ïy)Il ( 2 - 2 7 )
y
J 0 U + t)dy=U *U t M = (T /fi)'U
0 m (2*31)
Now, denoting
r
' n '
=
2 '
T -**'*
2 '
T s
y
- "* '
2 J
W.=-4^>
fi ox fi dy
^=-~-"
fi dz
(2-32)
J
( 2 . 2 l ) - ( 2 . 2 5 ) become n e a r l y t o
P + PV, = (, + dT T1 I dy
PV = e, + 9 ( 7 ,
+R)/dy (2-33)
Ï rI
py> = «s + ar^/a.
+ + lu+J
' Tx j[Ty Tx) =Ty{" j) ^ .
-9-
s m a l l e r - s c a l e t u r b u l e n c e , and so t h e t u r b u l e n t dissi-
p a t i o n i s n e a r l y given by
2.3 S p e c t r a l d e n s i t y f u n c t i o n o f t u r b u l e n t e n e r g y
2.3.1 R e l a t i o n between one- a n d t h r e e - d i m e n s i o n a l s p e c t r a
Define t h r e e - d i m e n s i o n a l spectrum E(k) and one-dimensional
spectrum S(k) i n t h e f o l l o w i n g : f f- f-
j" ?-EU)<*=£
0 (2-39) ^ (2M0)
Then, i n t h e i s o t r o p i c t u r b u l e n c e we can o b t a i n
-10-
the f i n a l stage where n i s dominant (-7 power l a w ) ,
(c) Inertial subrange ( i n t e r m e d i a t e - s c a l e eddy)
When t h e Reynolds number R = *'Lj » i s very l a r g e , t h e r e i s
L
LJ ,=
1 aV(2/KC)**R * L . LjX= a V * / / Ï 5 • ( 2 / *C)** « ^ L (2-52)
-11-
By comparison w i t h (2.47) and (2.55),
K„= ( 8 / 9 )( 5 5 C / 1 8 ) -V* (2'57)
where, B' = ( 5 6 7 / 5 5 ) - y a . 4
t . s 10(1,/?) .
f2J Denoting m
J' ' j ) rj k= S k d fi
-12-
i n t o t h e s m a l l e r - s c a l e eddies via a cascade process, and e v e n t u a l l y
dissipated i n t o heat energy. This phenomenological change
process of energy may he analogous t o t h a t o f 'water s e r v i c e ' , as
shown i n F i g . 2.9. According t o t h i s model, t h e r a t e e o f
t u r h u l e n t d i s s i p a t i o n can be evaluated from t h r e e d i f f e r e n t methods:
(A) 'Reservior'-method (corresponds to productive subrange)
Use e=K-u"/L x (Eq.(2.50)).
(B) 'Conduit '-method (corresponds to inertial subrange)
Use « = sjk)/c)** (Eq. (2.47)).
(C) 'Terminal tube'-method (corresponds to viscous subrange)
l
Use e=isvu' j~k'sj.k)dk (Eq. (2.43)).
3 !
K
~ = ™.K= D
U
2 ( 3 r / 2 ) / ( 5 5 C / 1 8 ) "V* , ( r = 0.5) (2-63)
-13-
2.4 S u b d i v i s i o n o f f l o w f i e l d a n d i t s t u r b u l e n c e c h a r a c t e r i s t i c s '
2.4.1 S u b d i v i s i o n o f o p e n - c h a n n e l f l o w f i e l d
Since t h e r e may e x i s t an analogy between s p e c t r a l and w a l l -
t u r b u l e n c e spaces, as described i n F i g . 2.13, an open-channel f l o w
f i e l d can be a l s o d i v i d e d i n t o t h r e e subregions (see F i g . 2.14),
i n t h e same manner as s p e c t r a l space.
(1) mil region (y+ <^ 100)
This region corresponds to an inner layer in boundary layer,
and is dominated by and v (inner parameters). As men-
tioned later, the bursting phenomenon occurs most remarkably
in this region.
(2) Free-surface region (0.6 < £ <, 1.0)
This region corresponds to an outer layer, where the turbul-
ent characteristics are dominated by U m a x and h (outer
parameters).
(3) Equilibrium region (lOOPff <, E, <, 0.6)
This is an intermediate region between (1) and (2) subregions,
where a dynamically equilibrium state for turbulent energy
budget, i.e. P - e, is nearly realized.
+
2.4.2 Wall r e g i o n (, «s 100)
We have confirmed e x p e r i m e n t a l l y t h a t
tv/O* =A { •( y V . (Ai i s c o n s t a n t ) (2-64)
L s
l L u*/<< = * x
+
• (y )'/ 2
. ( A i s constant)
2 (2-65)
From (2.50),
2 66
»'/(/, = /(, •(/)"'/, A, = UtAjKVA (- >
u'/U has a -1/6 power law since A
# 3 becomes n e a r l y constant when R L
-14-
Next, assuming phenomenologically that
JL _ , + <wf_ ( A14 i s a p r o p o r t i o n a l c o n s t a n t )
A (2-67)
u
* <Ly
we can also o b t a i n from (2.66)
+ +
U = A lA -{y ) •/*) ,
t e ( y > 30)+
(2-68) .
2.4.3 F r e e - s u r f a c e r e g i o n (0.6£{£-i-0)
We have confirmed e x p e r i m e n t a l l y t h a t
( 2 7 l )
'
-15-
From (2.11) and (2.75), we can approximately o b t a i n
L
2 ' f'+l J (2-76)
0 r
' ^ ^ s i n g v^T^t/, as a c h a r a c t e r i s t i c v e l o c i t y scale
nearly - '
3
th/U t = C {->
t ( 2 . 7 9 )
-16-
Next, since t h e energy "budget i n t h i s r e g i o n i s
u e >
u{- dt u,'
by s u b s t i t u t i n g (2.78) we can o b t a i n
+ 7
,/={/+ - t/ = - C . C 2 ( ln( 1 - / P ? ) + V ' l -f) - In f 3 (2-82)
+
Ur=U+
-tf =- * - ' {ln(l-vT=7) + V l - f ) * r (2-83)
A war
And a l s o , (2.14) o f P r a n d t l can be extended i n t o t h i s r e g i o n ,
which r e s u l t s i n Vp^u*^-u +
=- inf (2-84)
Consequently, {/ = « c . .—£ £_ (2-85)
" (2-1)
Since Ci i s about 3 (see Chapter 4 ) and thus KC\ -1, (2.82) can be
deduced by averaging (2.83) o f Karman and ( 2 . 8 4 ) o f P r a n d t l weighted
with a r a t i o of 2 t o - 1 . This f a c t suggests t h a t t h i s r e g i o n
overlaps t o some extent w i t h t h e w a l l r e g i o n . F i g . 2.18 shows
U, U
N K and Up obtained from t h e above. And, t h e i r agreement
i s also good.
2.4.5 O v e r a l l d i s t r i b u t i o n o f mean v e l o c i t y
Since we w i l l deal i n d e t a i l w i t h t h e o v e r a l l d i s t r i b u t i o n s o f
turbulence c h a r a c t e r i s t i c s i n Chapter 4, we here consider only
those o f mean v e l o c i t y . Because a f u l l y developed open-channel
f l o w may be dominated by an a c t i v e motion o r l o c a l s i m i l a r i t y , as
w i l l be shown i n Chapter 4, t h e previous c o n s i d e r a t i o n on a basis
of t h e m i x i n g - l e n g t h t h e o r y i s s t i l l v a l i d . So, i t i s an essen-
t i a l work t o determine a c c u r a t e l y t h e d i s t r i b u t i o n o f m i x i n g - l e n g t h .
Fig. 2.19 shows t h e numerical s o l u t i o n s o f ( 2 . 1 l ) obtained by using
(2.77) m u l t i p l i e d by t h e damping f a c t o r o f ( 2 . 1 3 ) , t o g a t h e r w i t h
the experimental values (see Chapter 4 ) . A good agreement b e t -
ween t h e both i s noticeable. And, i t i s concluded t h a t even
(2.l4) o f Prandtl i s applicable f a i r l y accurately i n y +
> 30.
2.5 T u r b u l e n c e c h a r a c t e r i s t i c s o v e r rough bed
I n t h i s s e c t i o n we have e s t i m a t e d t h e e f f e c t o f wall-roughness
upon t h e t u r b u l e n c e c h a r a c t e r i s t i c s , by comparison w i t h those over
the smooth bed o b t a i n e d i n previous s e c t i o n s . However, some
estimations o b t a i n e d i n t h i s s e c t i o n a r e o m i t t e d here since t h e
e f f e c t o f wall-roughness w i l l be i n d e t a i l i n v e s t i g a t e d i n Chapter k.
-IB-
ü and v, i n order t o p r e d i c t t h e magnitude o f t h e c o n t r i b u t i o n
t o t h e Reynolds s t r e s s from each event.
9 1 1 <2 91)
^-hwït ^ " -" '
Expanding <D(£, 7j) in a Taylor series about £ = rj = 0 yields
, 7 ) = + n ( I , ( ê , 7 ) l £ 0 ( 2 9 3 )
^ n l o ^ ( ^ ^') ' ' ''''" ' *
Thus i f and Q correspond to the coefficients in Taylor expansions of 0 ( £ , 17)
Jfc jk
and In C>(£, 17) respectively. The relations between the moments and the cumu-
lants are successively obtained by making use of (2.91)-(2^93) •
Now, since M = ¥ = 0, M = I = 0, M = & = 1, M = S = 1 and J l f =
10 01 20 02
2
u
Ql0 = M -2,
lQ Q =M +3R,
3l 31 Q 22 = M -2R*-l 2l ƒ
Q Jk for j < k can be obtained by merely exchanging j and k in the terms of Q jk for
\n®(£, ) V = - m 2
- M Z v + y 2
) + . Z J i + k
£ f P v[
k
- (2.95)
Through an inverse transformation of (2.90) in which the terms of <!>(£;, 17) of less
than fourth order are taken into account, p{ü, Ö) can be written as
-kZiwfafiSSP (3.96)
-19-
-*<M)[i+J^2feM)]. (.
2 97)
G M
= 2n(l-R>)i V\- eX
2(1-JP) )' (2.98)
and i7 . (Ö, 0) is a Herraite polynomial in two variables.
u
where R (ü, v) is an adjoint Hermite polynomial in two variables (see Frenkiel &
jk
Klebanoff 1973).
The probability distribution of one variable is much simpler and is derived in
the same manner as (2.96):
4
0 „ 3' 1 / H\ 2
the Gram-Charlier distribution becomes the same as the Gaussian one, and thus it
may be said that Q gives a measure of the skewness or intermittency of the
jk
distribution.
Now consider the probability distribution p {w) of the normalized Reynolds w
R R 2 w
c : : r i ( ) r- / v + w m
e X P C X P
7r(l-^)i \l-^jJ 0 \ 2(1-*») )
dü ,
1+ S ^\{H {%-Rwiü)+B {-%RwlÜ)}
jk ik
-£•(2.102)
Since H jk is an odd function for when j + k is odd and vice versa, the third-order
Cumulants in the correction term of p (w), i.e. the second term in (2.102), vanish. w
The third-order cumulants, which are closely connected with turbulent diffusion,
as will be mentioned later, are much more important quantities than the fourth-
order cumulants. Because they disregarded this third-order cumulant, some
previous studies such as those by Lu & Willmarth (1973) and Antonia & Atkinson
(1973) could find little obvious relation between the bursting process and the
probability distribution of the Reynolds stress.
-20-
From the above description, a conditional probability distribution should be
introduced in order to evaluate the effect of the third-order cumulants, while, for
simplicity, the fourth-order cumulants, which are less important, may be omitted.
We shall denote the probability distributions of each event shown in F i g . 2.21
by p ) (outward interaction), p (u>) (ejection), p (w) (inward interaction)
l { w t 3
2 2
where K is the i>th-order modified Bessol function of the second kind. Hence
v
+ Z>+
W -( ^^
+D+
f + { w ) m §.c"K \t\)
i{ TJ^{<*+*> {T ) ))' (2.113)
-21-
3
S and S are the skewness factors of u and v respectively, as S =ft and S = t) .
u v u
3
v
2 2
as D = ü.v and D = v.ft, ,and here we shall call them the diffusion factors.
u v
for the sequence of the bursting process and that they are closely connected
with the turbulent diffusion.
Next, taking into account the partition level H in the diagram of Reynolds
stress w = uvjuv according to Lu & Willmarth (1973), as shown in F i g . 2 . 2 1 ( b ) ,
the contributions to the Reynolds stress can be associated with one of five events
including a hole event when |u>| < H. The hole event is labelled event 5. Then the
time fraction T^H) and the contribution to the Reynolds stress RS^H) corres-
ponding to each event can be represented by
I Pi(w)dw (» = 2,4),
T (H) =
4 (2.116)
r-u
p (w)dw
t (t = l,3),
J -00 /
TJH) = f p (w) dw = 1 - S
w (t h e h o l e e v e n f c
) (2.117)
j-n i=i
ƒ
H
( ƒ
"-H
wpi{w)dw > 0 (i = 2,4), j
J -ii f=i
When H = 0, the above equations describe the contribution of each event given
by F i g . 2.21 (a ) . When H > 0, it may be expected that the characteristics of each
event such as skewness and intermittency can be made clear. Also, some relation-
ships between the coherent vortex motion with turbulent production which was
observed by Corino & Brodkey (1969) or Kim etal. (1971) through flow visualiza-
tion and the data obtained in this study through point measurements may be
discussed by varying H as a parameter.
-22-
The average value Z,(B) o f u which produces t h e sweeps w i t h w=H
i s a l s o given i n t h e same manner, S u b s t i t u t i o n o f (2.97) and
(2.109) i n t o (2.120) y i e l d s
/ z r , ( z ) - «,(z)
(2-121) = 1 - (2-122)
X o ( Z ) / r - «*(Z)
+ <?<») / 2
, fljisCOx + a o . O / e + B 130,, +
fa.U) = ~ { ( B . , Z + B o i Z ) X o U ) +
2
and ( 3 - 2B X>os) / 6
ft,BltU* + 3)(<?>o+30,i)/6 +
«,(Z) = 2*#"(B) • « p {-RZ\/R
2
(2-123) (3B + l)(<? 03 + 30 )/6 1 1
l
fl„ = ( 3 B + l ) ( 0 » + 3 0 I l )/6 +
B(B'+ 3)(0 0 S + 30„)/6
(2-124)
and Ji]w^B)
P o f u which produces t h e e j e c t i o n and t h e sweeps w i t h
w > H, r e s p e c t i v e l y , are g i v e n i n t h e f o l l o w i n g s by i n t e g r a t i n g
(2.97) i n respect o f ¥ > H :
2 1 2 6 )
p,(£UèB) = i « p (- f ) ^ U U > " *,<»)> exp ( -
-23-
References (Chapter 2)
1) Rotta, J. C. : Turbulente Strömungen, B. G. 18) Kline, S. J. et al. (edi.) : Computation of turbulent
Teubner, 1972. boundary layers - - - 1968 AFOSR-IFP-Stanford
2) van Driest, E . R. : On turbulent flow near a wall, Conf., vol.1, 1968.
J. Aeron. Sci., vol.23, pp.1007-1011, 1956. 19) Reichardt, H . : Vollstandige Darstellung der
3) Schlichting, H . : Boundary layer theory (6-th edi.), turbulenten Geschwindigkeitsversteilung in glatten
McGraw-Hill, pp.560-625, 1968. Leitungen. ZAMM, Bd.31, pp.208-219, 1951.
4) Schubauer, G. B. : Turbulent processes as ob- 20) Davies, P. O. A. L & Yule, A. J. : Coherent
• structures in turbulence, J. Fluid Mech., vol.69.,
served in boundary layer and pipe, J. Applied
PP.513-537, 1975.
Phys., vol.25, pp.188-196, 1954.
21) Kline, S. J . , Reynolds, W. C , Schraub, F. A. &
5) Tennekes, H. & Lumley, J. L : A first course in
Runstadler, P. W. : The structure of turbulent
turbulence, MIT Press, 1972.
boundary layers, J . Fluid Mech., vol.30, pp.741¬
6) Bradshaw, P. : 'Inactive' motion and pressure
773,1967.
fluctuation in turbulent boundary layers, J. Fluid
Mech., vol.30, pp.241-258, 1967. 22) Kim, H. T., Kline, S. J. & Reynolds, W. C : The
production of turbulence near a smooth wall in
7) Hussain, A. K. M. F. & Reynolds, W. C. : The
a turbulent boundary layer, J. Fluid Mech.,
mechanics of an organized wave in turbulent
vol.50, pp.133-160, 1971.
shear flow, J. Fluid Mech., vol.41, pp.241-258,
23) Corino, E . R. & Brodkey, R. S. : A visual investi-
1970.
gation of the wall region in turbulent flow, J .
8) Batchelor, G. K. : Homogeneous turbulence,
Fluid Mech., vol.37, pp. 1-30, 1969.
Camb. Univ. Press, 1953,
24) Grass, A. J . : Structural features of turbulent flow
9) Inoue, E . : On the structure of wind near
over smooth and rough boundaries, J. Fluid
the ground, Bull, of the National Institute
Mech., vol.50, pp.233-255, 1971.
of Agricultural Sciences, Ser. A , No. 2, 1952.
25) Brodkey, R. S., Wallace, J. M. & Eckelmann, H. :
10) Nakagawa, H., Nezu, I. & Ueda, H. : Turbulence Some properties of truncated turbulence signals
of open channel flow over smooth and rough in bounded shear flows, J. Fluid Mech., vol.63,
beds, Proc. of J S C E , No.241, pp.155-168, pp.209-224, 1974.
1975. 26) Nakagawa, H. & Nezu, I. : Prediction of the
11) Tchen, C. M. : On the spectrum of energy in contributions to the Reynolds stress from the
turbulent shear flow. Jour, of R.N.B.S., vol.50, bursting events in open channel flows, J. Fluid
pp.51-62, 1953. Mech., vol.80, pp.99- 128, 1977.
12) Hinzc, J. O. : Turbulence, McGraw-Hill, pp.451¬ 27) Frenkiel, F . N. & Klebanoff, P. S. : Probability
566, 1959. distributions and correlations in a turbulent
13) von Karman, T. : Progress in the statistical theory boundary layers, Phys. of Fluids, vol.16, pp.
of turbulence, Proc. of N.A.S., vol.34, pp.530¬ 725-737, 1973.
539, 1948. 28) Lu, S. S. & WUlmarth, W. W. : Measurements of
14) Monin, A. S. & Yaglom, A. M. : Statistical fluid the structure of the Reynolds stress in a turbulent
mechanics, MIT Press, 1971. boundary layer, J. Fluid Mech., vol.60, pp.481-
15) Bradshaw, P. : Conditions for the existence of 511, 1973.
an inertial subrange in turbulent flow, A.R.C., 29) Antonia, R. A. & Atkinson, J. D. : High-order
R. & M., No.3603, 1967. moments of Reynolds shear stress fluctuations
16) Zagustin, A. & Zagustin, K. : Analytical solution in a turbulent boundary layer, J. Fluid Mech.,
for turbulent flow in pipes, La HouiHe Blanche, vol.58, pp.581-593. 1973.
No.2, pp.1 13-1 18, 1969.
17) Townsend, A. A . : Equilibrium layers and wall
turbulence, J. Fluid Mech., vol.11, pp.97-120,
1961.
-2k-
C h a p t e r 3 Method o f t u r b u l e n c e m e a s u r e m e n t and i t s a n a l y s i s
3.1 Introduction
The method o f t u r b u l e n c e measurements, techniques o r instruments
can be d i v i d e d b r o a d l y i n t o two groups : one i s point-measurement
t h a t a d e t e c t i n g element i s i n t r o d u c e d i n t o t h e f l o w i n g f l u i d , and
another i s f l o w v i s u a l i z a t i o n t h a t a t r a c e r or other i n d i c a t o r i s
i n t r o d u c e d i n t o t h e f l u i d t o make t h e f l o w p a t t e r n v i s i b l e .
I n t h i s study, i . e . turbulence-measurements i n water f l o w , as t h e
former method a h o t - f i l m anemometer and as t h e l a t t e r method a
hydrogen-bubble t r a c e r are adopted here , r e s p e c t i v e l y .
3.2 P r o p e r t i e s o f d u a l - s e n s o r h o t - f i l m anemometers
3.2.1 P r i n c i p l e o f h o t - f i l m anemometer
As shown i n F i g . 3.1, i n a h o t - w i r e or h o t - f i l m anemometer
two b a s i c a l l y d i f f e r e n t methods can be a p p l i e d : (a) c o n s t a n t - c u r r e n t
method and ( b ) constant-temperature method. Since we adopted a
constant-temperature anemometer manufactured by DISA, o n l y (b)-method
i s considered here.
W e l l , an equation o f heat t r a n s f e r from c y l i n d r i c a l h o t - f i l m
sensor (see F i g . 3.2) was g i v e n by Kramers:
Since ., _ EL/R* /„ 2L
2 ai
we can o b t a i n E - E\ = A'u (3-3)
0
A' = Const. X W ( 0 / O - ' P -*R,{T,- T„)
r
a
(3-4)
E s Const. Xkl P ° *R,( T ~ T )
2
T € (3-5) K
-25-
3.2.2 Directional sensitivity o f hot-film sensor
We adopted a dual-sensor h o t - f i l m prohe(DISA 55A89, as shown
i n F i g . 3.U) i n order t o measure t h e t u r b u l e n c e i n t e n s i t i e s ,
Reynolds s t r e s s and o t h e r s .
W e l l , according t o Hinze, t h e e f f e c t i v e v a l u e o f U can be
obtained a p p r o x i m a t e l y from t h e r e l a t i o n : W=uKco*'*+khm**) (.
3 7)
!
(£ -^) =o fl \ "TTT J
Denoting . „, _ r(4'-£.) f*>
I /i 2 ( <=1, 2 )
F F s (3-10)
, I
2(cos *i + i Y s i n * )
r < ( a ) =
r m o r -0-12)
(3-11) ^
An experiment o f t h e d i r e c t i o n a l s e n s i t i v i t y o f t h i s probe has
been performed i n a homogeneous water f l o w (U=3^.3 cm/s, u'/U=0.0U).
Fig. 3.6 shows t h e e x p e r i m e n t a l v a l u e s , t o g a t h e r w i t h t h e curves o f
(3.11). And, t h e values o f k 2
c a l c u l a t e d by (3.12) are a l s o
shown i n F i g . 3-7. From these r e s u l t s , t h e average values o f k
became equal t o ki=0.12 and k = 0 . 2 1 , which were o f t h e same order t o
2
-26-
3.2.3 C a l i b r a t i o n o fdual-sensor h o t - f i l m anemometer
F i r s t l y , we examined t h e temperature dependence o f t h e e l e c t r i c
r e s i s t a n c e o f t h e h o t - f i l m , which r e s u l t s are shown i n F i g . 3.8.
Fig. 3-9 shows t h e o u t p u t v o l t a g e s o f anemometers i n t h e s t i l l water
(i.e. U=0). These r e s u l t s agree well with ( 3 - 5 ) , and g i v e
El =0.212(33.40- 7 „ ) . £ ^ = 0 . 1 7 8 ( 3 7 . 4 2 - 7 » ) (3-18)
3.3 E f f e c t o f v a r i o u s f a c t o r s on t h e p r o p e r t i e s o f h o t - f i l m anemometer
3.3.1 E f f e c t o f t h e n o i s e
F i g . 3.13 shows t h e data processing system c o n s t r u c t e d i n t h i s
study. Consequently, i t i s apprehended t h a t some noises may be
i n v o l v e d i n t o t h e analogue data o f h o t - f i l m anemometers u n t i l they
are converted i n t o t h e d i g i t a l data.
Wow, i f U and Tw vary i n f i n i t e s i m a l l y , t h e output v o l t a g e E
v a r i e s by AE= (OE/OU) AU+( °E/dT ) a ÓT W (3-21)
-27-
CU i s t h e s e n s i t i v i t y f a c t o r o f t h e anemometers, and F i g . 3-lU
shows an example o f C^-curves c a l c u l a t e d from ( 3 . l 8 ) - ( 3 . 2 0 ) and
(3.22). I f the noise i s expressed by AE, AE/E describes t h e
n o i s e - s i g n a l r a t i o of t h e analogue data. From t h i s f i g u r e ,
AU/U becomes n e a r l y equal t o (5-7)AE/E. That i s , N-S r a t i o o f
even o n l y 1% i n t h e instruments(see F i g . 3-13) causes t h a t o f
(5-7)$ i n t h e e v a l u a t i o n o f v e l o c i t y .
Thus, we made some devices t o minimize t h e noise
i n v o l v e d i n t o t h e analogue data.
B = ^ L _ ( ^ + { , « ^ L ) (3-23) .
T X 1
V C
^^2\ '-E- ^ifj (3-25) .
-28-
Since t h e v a r i a t i o n o f t h e s e n s i t i v i t y f a c t o r Cj_ a g a i n s t t h e water
temperature changes i s very small as seen i n F i g . 3-lU, AU/U o r AV/U
(ex. t h e r e l a t i v e t u r b u l e n c e i n t e n s i t y u'/U o r v'/U) can be e v a l u a t e d
f a i r l y a c c u r a t e l y i n s p i t e o f ATw.
In c o n s i d e r a t i o n o f t h e above, a s t a b l e water temperature was
here maintained d u r i n g t h e o p e r a t i n g p e r i o d by c i r c u l a t i n g water
throughout t h e day before t h e t e s t began i n order t o d i m i n i s h t h e
e f f e c t o f i t s changes on t h e h o t - f i l m as much as p o s i b l e , and i n
• s u i t AT became w i t h i n only 0.5°C a t t h e end o f a l l measurement t e s t s
W
3.3.3 E f f e c t o f t h e nonhomogeneity o f v e l o c i t y a l o n g t h e s e n s o r
When t h e l e n g t h I o f t h e h o t - f i l m sensor i s l a r g e r than t h e
microscale X o f t u r b u l e n c e , t h e e f f e c t o f t h e nonhomogeneity o f
v e l o c i t y along i t s sensor should n o t be neglected. I n t h e same
manner as t h e d e s c r i p t i o n by Hinze, t h e f o l l o w i n g r e l a t i o n s as t o
V-type sensors(cf. F i g . 3-5) are obtained approximately:
(cf. 2.2) and our sensor has 1=1 mm, t h e above e f f e c t may be
w i t h i n only 1% when h=8 cm.
3.3.4 Other e f f e c t s
To d i m i n i s h t h e e f f e c t o f i m p u r i t i e s i n t h e water upon t h e
characteristics o f t h e h o t - f i l m , t h e suspended m a t e r i a l s i n t h e
f l o w were f i l t e r e d by gauzes. And, v a r i o u s devices have been
done as t o i n s u l a t i o n badness o f probe connector by t h e submerged
water, s e l f - v i b r a t i o n o f t h e experimental channel and o t h e r s .
-29-
3.4 H y d r o g e n - b u b b l e t e c h n i q u e
3.4.1 P r i n c i p l e o f h y d r o g e n - b u b b l e t r a c e r a n d i t s a n a l y s i s
As shown i n F i g . 3 . l 6 , a f i n e p l a t i n u m w i r e o f d=50 um diameter
i s s t r e t c h e d i n the water i n t h e d i r e c t i o n o f y- or z - a x i s .
This w i r e forms t h e negative e l e c t r o d e o f a pulsed D.C. c u r c u i t o f
(100-500) v o l t s , where i n r e s u l t hydrogen bubbles w i t h a diameter
of ( 0 . 5 - 1 . 0 ) d are produced and t h e y form t h e markers f o r f l o w v i s u a -
lization. I n order t o o b t a i n a marking i n space (combined t i m e -
s t r e a k markers), t h e w i r e has s h o r t s e c t i o n s a t n e a r l y r e g u l a r i n t e r -
v a l s (about 2 mm) coated f o r i n s u l a t i o n .
Fig. 3-17 shows t h e data processing system f o r hydrogen-bubble
method. Any frame i n l6mm f i l m s i s enlarged by t h e f i l m - p r o j e c t o r
and t h e coordinates o f some hydrogen-bubbles can be read t o t h e l i m i t
of 0.1mm by t h e instrument o f d i g i t a l c o o r d i n a t o r . F i g . 3.18 shows
an a n a l y s i s method o f instantaneous v e l o c i t y p r o f i l e ( s i n g l e frame
method). Denoting t h r e e b u b b l e - p o i n t s near t h e r e f e r e n c e a x i s x r
References (Chapter 3 )
1) Hinze, J. O. : Turbulence, McGraw-Hill, pp.73- 3 ) Schraub, F. A., Kline, S. J . , Henry, J., Runstadler,
1 4 1 1 9 5 9
• - P. W. & Littcll, A. : Use oNiydrogcn bubbles for
2) Bradshaw, P. : An introduction to turbulence and quantitative determination of time-dependent
its measurement, Pergamon Press, pp.85-133, velocity fields in low-speed water flows, Trans, of
1 9 7 1
- ASME, BE, pp.429-444, 1965.
-30-
C h a p t e r 4 E x p e r i m e n t a l c o n s i d e r a t i o n on t u r b u l e n t s t r u c t u r e
and e n e r g y budget
4.1 Introduction
I n the recent t u r b u l e n c e research a ' s e l f - c o n s i s t e n t ' characte-
r i s t i c of t u r b u l e n c e has been n o t i c e d t h a t a t u r b u l e n c e - c h a r a c t e r i s t i c
may be more c l o s e l y r e l a t e d t o another turbulence-characteristic
than mean-velocity one, t h a t i s , t h e r e may be a stronger i n t e r r e l a -
t i o n among the t u r b u l e n c e c h a r a c t e r i s t i c s .
I n t h i s chapter, we i n v e s t i g a t e i n d e t a i l t h e t u r b u l e n c e
i n t e n s i t i e s , Reynolds s t r e s s , s p e c t r a l d i s t r i b u t i o n s , t h e r a t e s o f
p r o d u c t i o n , d i f f u s i o n and d i s s i p a t i o n o f turbulence,and others
i n open-channel f l o w s from t h e v i e w p o i n t o f t h e s e l f - c o n s i s t e n c y
of t u r b u l e n c e o r s i m i l a r i t y law mentioned above, and then we make
c l e a r s y s t e m a t i c a l l y the e f f e c t s of t h e Reynolds number, the Froude
number and the w a l l roughness upon these t u r b u l e n t s t r u c t u r e and
the energy budget.
4.2 Experimental equipment and data a n a l y s i s
Three groups o f experiments on two-dimensional fully developed
t u r b u l e n t f l o w i n an open channel were conducted i n a t i l t i n g flume
15m l o n g , 50cm wide and 30cm deep, as shown i n F i g . k.l. The chan-
nel slope could be changed by a d j u s t i n g two j a c k s so t h a t normal
f l o w could be obtained. A few b a f f l e s and screens t o prevent t h e
occurrence o f l a r g e - s c a l e disturbances were set up a t t h e entrance o f
the channel, and consequently a f u l l y developed t u r b u l e n t f l o w was
obtained a t the t e s t s e c t i o n 9-5m downstream o f t h e entrance.
As described i n Table U . l , f i r s t group (a) c o n s i s t e d o f t h r e e
d i f f e r e n t kinds o f the Reynolds number ReEU h/v, i . e . R e = ( l . 2 ,
m 2.9,
1+
l+.6)xl0 , where the Froude number Fr^U /Vgh m and the w a l l e q u i v a l e n t
sand roughness k+Ek U*/v were kept t o be n e a r l y c o n s t a n t , i . e . Fr=0.T
s
-32-
2
q i n t h e e q u i l i b r i u m r e g i o n , as s i m i l a r t o t h e p r o p e r t i e s o f t h e
other kinds o f energy such as t h e r m a l energy. Then,
1 d
Jll = li!2Üj =-;<o-(4.3)
2
9 dC d(
The above conception may be a l s o s i m i l a r t o t h a t o f Townsend(196l)
or Zagustin e t a l . ( l 9 6 9 ) - By assuming t h a t X i s constant and
a p p l y i n g t h e r e s u l t s o f Il-eddy model proposed i n Chapter 6, t h e
f o l l o w i n g u n i v e r s a l functions of turbulence i n t e n s i t i e s are obtained
X( w xe
from (i+.3): u' / <tf \ ** n f \ '-n ( ^ <* n\
I n t h i s s e c t i o n , we discuss e x p e r i m e n t a l l y t h e v a l i d i t y o f
(h.7) by s y s t e m a t i c a l l y v a r y i n g t h e Reynolds, t h e Froude numbers
and t h e w a l l roughness ( c f . Table k.l).
-33-
(3) Reynolds-stress method
The friction velocity can be evaluated from (2.5) by using
the measured Reynolds stress distribution.
(4) Velocity gradient method
By definition, v* = * du/a | y y=0
where,
"*=ƒƒ u(x. y, t ) exp— i ( kx + at) dxdt, "*=jj v(x,y. t ) exp—i ( kx+wt) dxdt (4-9)
4.3.3 E f f e c t o f t h e F r o u d e number on t u r b u l e n c e i n t e n s i t i e s
F i g . h.3 shows t h e e f f e c t o f t h e Froude number on t h e t u r b u l e n c e
intensities. Except f o r Case G-2, t h e experimental values o f u'/U #
J 2
u ' / V = 0.55, i.'7ir =0.17. L/79^0.28 (4-13)
4.3.4 E f f e c t o f t h e w a l l r o u g h n e s s on t u r b u l e n c e i n t e n s i t i e s
Fig. U.5 shows t h e e f f e c t o f t h e w a l l roughness on t h e t u r b u -
lence i n t e n s i t i e s . These experimental values agree w e l l w i t h
those o f Grass(1971) o b t a i n e d by u s i n g t h e hydrogen-bubble technique
I n t h e r e g i o n o f £ >0.3, t h i s roughness e f f e c t may become f a i r l y
weak, and consequently t h e experimental values s t i l l c o i n c i d e w e l l
w i t h (U.7) independently o f t h e w a l l roughness. I t should be,
however, n o t i c e d t h a t t h i s e f f e c t appears remarkably near t h e w a l l .
So, t h e experimental values near t h e w a l l are i n d e t a i l shown
+
against y i n F i g . U.6. The values o f u'/U* i n t h e w a l l r e g i o n
g r a d u a l l y decrease w i t h increase o f k , b u t v'/U* and w'/U are
s #
-36-
+
l a y e r y =10-30 by t h e e j e c t i o n s and sweeps due t o t h e f l o w i n s t a b i -
lity (see Chapter 5 ) , a b u f f e r - l a y e r over a rough bed disappears
p e r f e c t l y or p a r t l y i n t o t h e roughness elements, and consequently
the turbulence would be produced i n another way. Actually, the
p o s i t i o n y^~(15-20) where u'/U* a t t a i n s maximum becomes under t h e
a n d t n u s
t o p o f roughness element when +/b, > +, i . . kg >T0 >
k y e
+
i s examined against y i n F i g . h.l. I t should be n o t i c e d t h a t
U7£/,)(t/)"'/ has r e a l l y a u n i v e r s a l c h a r a c t e r i s t i c i r r e s p e c t i v e o f
+ s
-37-
4.4 R e l a t i v e t u r b u l e n c e i n t e n s i t y and R e y n o l d s s t r e s s
4.4.1 R e l a t i v e t u r b u l e n c e i n t e n s i t y u'/U
The c h a r a c t e r i s t i c s o f only u'/U are i n v e s t i g a t e d here since
v'/U and w'/U are e a s i l y evaluated from (h.12).
Fig. k.8 shows t h e values o f u'/TJ o b t a i n e d by v a r y i n g o n l y
t h e Reynolds number Re. I n % < 0 . 1 , u'/U decreases remarkably
w i t h increase o f Re, but i n £ >0.1 t h i s e f f e c t o f Re becomes weaker.
And, such a c h a r a c t e r i s t i c i s a l s o seen i n the experimental curves
o b t a i n e d by L a u f e r ( l 9 5 l ) and B l i n c o e t a l . ( l 9 7 l ) . Wow, t h e f o l l o w i n g
r e l a t i o n can be obtained from (2.lU) and (h.l):
u! _ 2.3 e x p ( - f ) j
U £-'lnf+^(V(,) ' A{RJ=— l n * , + 5.5 (4-18)
1 1
U c In f - l t " ln(i,/A)-C (i/))r
-38-
Fig. k.H shows the values o f Reynolds s t r e s s o b t a i n e d at each
group i n Table h.l. They show a universal distribution indepen-
d e n t l y o f t h e Reynolds, t h e Froude numbers and t h e w a l l roughness.
Now, from (2.5) and ( 2 . l i t ) we o b t a i n -m _ i _
4.4.3 C o r r e l a t i o n c o e f f i c i e n t s o f R e y n o l d s s t r e s s
A c o r r e l a t i o n c o e f f i c i e n t o f Reynolds s t r e s s i s d e f i n e d as
Rs-üv/u'v', which i s a measure i n d i c a t i n g i t s s e l f - c o n s i s t e n c y o f
t u r b u l e n c e or s i m i l a r i t y .
Fig. h.12 shows the values o f R o b t a i n e d a t each group, i . e .
(a), (b) and ( c ) i n Table h.l. I n t h e w a l l r e g i o n and t h e f r e e -
surface r e g i o n r e s p e c t i v e l y , R increases and decreases monotonously
w i t h y/h, w h i l e i n t h e e q u i l i b r i u m r e g i o n i t remains n e a r l y constant
, i . e . R =0.*+. The comparison w i t h previous data i n boundary l a y e r
and pipe f l o w s shown i n F i g . h.12 i n d i c a t e s t h a t R d i s p l a y s the
u n i v e r s a l c h a r a c t e r i s t i c s , i r r e s p e c t i v e of t h e f l o w c o n d i t i o n s and
the roughness. Now, such a curve o f R can be p r e d i c t e d from
v
(I+.7) and (U.20), as f o l l o w s : R- ~~ - [izl±zll*i*ï? u «1
uV 2.92 exp ( - 2 f) •
-39-
c a l c u l a t i o n developed by Bradshaw et a l . ( l 9 6 7 ) .
More, from (h.l) and (h.20), we also obtain
MmzSL m (i-f)-U*.f)- ( j
» 7 9.56 e x ( - 2 f )
P
-ho-
depth, i . e . F r = l , and then decrease w i t h Fr i n t h e s u p e r c r i t i c a l f l o w .
Such a s o - c a l l e d resonant hehaviour near t h e c r i t i c a l - d e p t h may be
due t o an i n a c t i v e motion, as mentioned l a t e r .
By the way, t h e equation o f f r e e surface p a t t e r n can be l i n e a r l y
approximated and i n r e s u l t i t becomes
2
a%' F SH' IF *RS NS dk' 2s dk'
i i _ gpt) — 2B —- — ( —- ) - =0 (4-23)
V P T 2 2
dx U H dxdt K
U j dt h d x Uh m dt >
y u
where, J »~ * .
So, t h e existence c o n d i t i o n s o f t h e small amplitude wave o f (h.2h)
are given by <6 = yv/ 2. 4=2/)-N/2 (4-26) 4 - fi± J^p- 1) + j F * (4-27)
Consequently, denoting /? =e
r 2)* +fi-fin)-*- (4-29) ,
-1*1-
to the t u r b u l e n t energy, i . e . Q= — (4.30)
ƒ * ( / . / 2 ) ? dy
31)
Applying ( I t . 7 ) , (*t.30) y i e l d s a
= {-^) ' JJ\—) U
*
Next, t h e s p e c t r a l a n a l y s i s o f t h e surface-wave f l u c t u a t i o n
was done by making use o f t h e F.F.T. method (the sample s i z e N was
chosen N=lt096), and i t s raw r e s u l t s were smoothed by a d i g i t a l
2
t r i a n g u l a r f i l t e r , and then they were normalized by h' .
Consequently, these a n a l y s i s o f t h e frequency spectrum F ( f ) s a t i s f i e d
w i t h i n e r r o r s o f 2-3% t h e normalized c o n d i t i o n t h a t f
[ F(f)df=l •
J 0
Now, these r e s u l t s o f F ( f ) are shown i n F i g . h.16 as a parameter
o f Fr. When Fr i s so small as(T) or .2'' shown i n F i g . U.l6, F(f)
has an e v i d e n t peak a t f = 2 - 3 (Hz), and f o r f > f , i t seems t o
m m
-1+2-
measured by using a gauge-type pressure t r a n s d u c e r . The h y d r a u l i c
conditions o f t h i s experiments (Group P) are i n d i c a t e d i n Table h.2.
F i g . U.1T shows t h e i n t e n s i t y p'^vp" " o f w a l l - p r e s s u r e f l u c t u a -
2
where, V i s t h e dimensionless L a p l a c i a n by h.
I f t h e t u r b u l e n t s t r u c t u r e i s dominated by t h e a c t i v e p a r t , t h e u n i -
+
v e r s a l expressions o f u =fin.( f . * , ) and v*fV = D -expO f ) , become v a l i d , m t
and consequently t h e f o l l o w i n g u n i v e r s a l f u n c t i o n i s a l s o o b t a i n e d
from ( U . 3 3 ) : P'/ u>
fi = '/x
P 0 = c(c, R,) (4-34)
J
Or, A = / ' ( _ i ) . ( 0 .fl.)= fun. ( 7 i , ) - f
r G r (4.35)
m
(k.38) i s s a t i s f a c t o r y w i t h e r r o r s o f 5% i n t h e r e g i o n o f £ >0.1
f o r t h e smooth bed since u'/U <0.1 (see F i g s . U.8-U.10).
Next, a p p l y i n g t h a t u/u =C' mai
1
and Ljk=ST (see 4 . 6 . 4 ) , (U.39)
becomes ky>(l/7) o r u >(7/f)^
X (4 . 40)
i
from F ( f ) o b t a i n e d i n the same manner (F.F.T. method) as F i g . k.l6.
Judging from t h i s f i g u r e , t h e sampling number N was reasonably chosen
N=4096. F i r s t l y , the d i s t r i b u t i o n of S (k) satisfys f a i r l y w e l l
u
4.6.2 S p e c t r a l d i s t r i b u t i o n n o r m a l i z e d by L x
I n order t o i n v e s t i g a t e s y s t e m a t i c a l l y the s p e c t r a l d i s t r i b u -
tions on the d i f f e r e n t f l o w c o n d i t i o n s , i t i s necessary t o normal-
ize them by a s u i t a b l e c h a r a c t e r i s t i c scale. As t h i s s c a l e ,
L (macro-scale) or ri(micro-scale) i s u s u a l l y chosen.
x Here, we
adopted Lx because of an easy comparison w i t h t h e t h e o r e t i c a l r e s u -
l t s obtained i n 2.3.
Now, assuming t h a t (2.48) i s v a l i d up t o k=0, we can o b t a i n
L =(x/2)S (0)
x U (4-42)
However, f o r t h e p r e s e n t , (k.k2) i s considered t o he a rough esti-
mate because o f (4.40). And, a l l t h e analyzed r e s u l t s were p l o t t e d
by making use o f X-Y p l o t t e r i n t h e computer c e n t e r , i n t h e manner
as described i n F i g . 4 . 2 1 .
F i g . 4.22 shows t h e s p e c t r a l d i s t r i b u t i o n s a t t h e t y p i c a l mea-
sured p o i n t s o f each r e g i o n f o r t h e smooth bed. Since R# ^Ej_ i s
n e a r l y equal t o 500, t h e t h e o r e t i c a l curve c a l c u l a t e d w i t h RL=500
in 2.3.2 i s d e s c r i b e d by a t h i c k l i n e i n F i g . 4 . 2 2 . Although t h e
e f f e c t o f t h e n o i s e i n v o l v e d may appear i n t h e h i g h e r wave-number
p a r t o f i - i ^ 7 o J t h i s observed s p e c t r a l d i s t r i b u t i o n s show a u n i v e r -
x
I f t h e value o f L x i s a c c u r a t e , Ym must c o i n c i d e w i t h Y t -
So, i f t h e value Ym measured from F i g . 4 . 2 4 i s d i f f e r e n t from t h e
s
t h e o r e t i c a l value y t > ^x i - c o r r e c t e d i n t h e f o l l o w i n g so t h a t Ym
:
coincides w i t h Y t L &l -
x m lr /rj**
t (4-45)
-46-
Now, t h e s p e c t r a l d i s t r i b u t i o n s normalized by u s i n g L x corrected
i n t h i s manner (see F i g . 4.21) a r e r e p l o t t e d i n F i g . 4.26. Actually,
such c o r r e c t e d d i s t r i b u t i o n s show a good agreement w i t h t h e t h e o r e -
t i c a l curves w i t h RL=1000 or 2000 i n t h e i n e r t i a l subrange whose
observed values are most r e l i a b l e . And, (a)-group i n Table h.l
a l s o showed t h e s i m i l a r r e s u l t s .
To sum up, a l l these data have confirmed t h a t t h e i n e r t i a l sub-
range where t h e -5/3 power l a w was v a l i d had t h e u n i v e r s a l c h a r a c t e r -
i s t i c s i r r e s p e c t i v e o f t h e Reynolds, t h e Froude numbers and t h e w a l l
roughness.
This f a c t suggests t h a t t h e i s o t r o p i c r e l a t i o n , i . e . a l o c a l i s o t r o p y
may be f a i r l y v a l i d i n t h e viscous subrange and t h a t t h e tendency
toward i s o t r o p y may be a l i t t l e s t r o n g e r f o r t h e rough bed t h a n t h e
smooth bed. This suggestion may not be c o n t r a d i c t o r y t o t h e r e s u l t s
o b t a i n e d i n 4.3 (e.g. F i g . 4 . 6 ) .
Next, a l t h o u g h t h e spectrum S ( k ) o f t h e Reynolds s t r e s s
u v i s zero
f o r any wave-number k i n t h e s t r i c t l y i s o t r o p i c t u r b u l e n c e , i t i s
e s s e n t i a l t o s (t)*o
uv i n t h e shear f l o w . F i g . 4.28 shows t h e co-
spectra S ( k ) o f -üv o b t a i n e d by F.F.T. method.
u v Though t h e r e a r e
f a i r l y l a r g e s c a t t e r s i n t h e observed v a l u e s , i t can be recognized
t h a t S ( k ) becomes smaller t h a n S ( k ) i n t h e range o f L -k>30 .
u v u x
, t o g a t h e r w i t h t h e t h e o r e t i c a l curves. Up t o t h e i n e r t i a l subrange
(the i n c r e a s i n g p a r t o f D (*)-*'/!>), t h e observed values agree v e r y w e l l
u
-hi-
an agreement "between t h e b o t h i s comparatively good although t h e data
i n t h i s range may be l i t t l e accurate.
4.6.4 C h a r a c t e r i s t i c e d d y - s c a l e s and t h e i r r e l a t i o n s
F i g . 4 . 3 0 shows t h e macro-scale L x evaluated from (k.k2) and i t s
c o r r e c t e d value L by u s i n g ( 4 . 1 + 5 ) , which were measured by a s i n g l e -
x
From t h i s f i g u r e , t h e f o l l o w i n g r e l a t i o n s are o b t a i n e d i r r e s p e c t i v e
k=B
o f F r :
Ljh = B-t+ , (f<0.6) (4.48) , U > - (f^0.6) (4-49)
-1+8-
B = 1 . 1 f o r t h e lower v e l o c i t y group. These r e s u l t e x p l a i n s w e l l
the v a r i a t i o n s i n F i g s . 4.30 and 4.31 ( o r see F i g . 4.34).
c a n t e
Next, t h e micro-scale A x d e f i n e d as x s u' / J(du/dxV
x evaluated
from t h e f o l l o w i n g t h r e e methods.
(1) Spectrum method : Evaluation from the dissipation spectrum D (k), u
c
• f - ( f ) - ' - y s s
4
As f o r t h e higher v e l o c i t y group (Re - 3 x l 0 ) , approximating t h a t A =A X
-49-
r e g i o n i s o b t a i n e d i f (4.49) i s reasonably used i n s t e a d o f (U.48)
i n (U.51). To sum up, i t should be n o t i c e d t h a t t h e m i c r o - s c a l e X x
4.7 T u r b u l e n t e n e r g y budget i n o p e n - c h a n n e l f l o w
4.7.1 T u r b u l e n t e n e r g y d i s s i p a t i o n
We (1975) have already examined the s u p e r i o r i t y among t h r e e k i n d s
o f the e v a l u a t i o n method o f d i s s i p a t i o n e proposed i n 2.3.3.
To be b r i e f , t h e accurate values o f L x and X x cannot be e a s i l y o b t a i n e d
, especially i n very l a r g e Reynolds number, as have been p o i n t e d out
i n F i g s . 4.30 and 4.32. Consequently, both (A) and (C) methods seem
-50-
t o "be rough e s t i m a t i o n . On the o t h e r hand, because t h e -5/3 power
law was recognized c l e a r l y i n every s p e c t r a l d i s t r i b u t i o n (see 4 . 6 ) ,
i t can be concluded t h a t t h e (B)-method, i . e . (2.1+7), i s most r e l i a -
ble. This c o n c l u s i o n i s a l s o g i v e n by Grant e t a l . ( l 9 6 2 ) , Bradshaw
(1967), Lawn(1971) and o t h e r s .
Now, F i g . 1+.35 (a) and (b) show t h e dimensionless d i s t r i b u t i o n s
o f t u r b u l e n t energy d i s s i p a t i o n eh/U^ e v a l u a t e d from (2.1+7) f o r t h e
higher and lower v e l o c i t y groups, r e s p e c t i v e l y . The observed values
of eh/uj i n both groups i n d i c a t e a s i m i l a r v a r i a t i o n , t h a t i s a u n i -
v e r s a l d i s t r i b u t i o n i r r e s p e c t i v e o f t h e Reynolds, t h e Froude numbers
and t h e w a l l roughness. Above a l l , i t should be n o t i c e d t h a t the
e f f e c t of roughness upon eh/lkj i s much smaller than upon u'/U o r L /h. # x
-51-
because o f the small v a r i a t i o n o f K vs. R L (see F i g . 2.11), and t h i s
3
shows the independence o f the d i s s i p a t i o n eh/U # on the Reynolds number
, as mentioned above.
Next, F i g . U.36 shows the dimensionless d i s t r i b u t i o n s o f d i s s i p a -
+
t i o n e Eev/u£ i n the w a l l r e g i o n . Although we here could not obtain
the data i n the viscous sublayer, Laufer(l954) succeeded i n the measu-
+
rement very near the w a l l i n a i r p i p e f l o w , i . e . y <10, and hence h i s
data are r e p l o t t e d i n F i g . 4.36. The observed values o f e +
also
show w e l l a universal characteristic. That i s , they a t t a i n the
maximum near the sublayer edge Ó+-10, and then decrease monotonously
w i t h increase o f y+. The t h e o r e t i c a l curves o f (4.56) and (2.78)
are i n d i c a t e d i n F i g . 4.36 by s o l i d l i n e s ( l ) and J ) , r e s p e c t i v e l y .
The observed values show a b e t t e r agreement w i t h (2.78) than (it. 56)
+
because (It. 7) o r (It.56) becomes i n a p p l i c a b l e near the w a l l o f y <50
as mentioned i n 4.3.3. On the o t h e r hand, when R* i s very l a r g e /
(2.78) becomes n e a r l y equal t o (2.6k) i n the w a l l r e g i o n , and thus
+
i t i s confirmed t h a t the d i s s i p a t i o n e over the sublayer edge can be
given by (2.6k) o r (2.78), namely the - 1 power law o f y . +
Then,
A! becomes n e a r l y equal t o C^.O.
4.7.2 P r o d u c t i o n and d i f f u s i o n o f t u r b u l e n t e n e r g y
The t u r b u l e n t energy p r o d u c t i o n PEE -üv3u/8y i s e a s i l y o b t a i n e d
from the measured Reynolds s t r e s s and mean v e l o c i t y d i s t r i b u t i o n s .
Over the e q u i l i b r i u m r e g i o n , the t h e o r e t i c a l d i s t r i b u t i o n i s g i v e n
from (2.1k) and (4.20), as f o l l o w s : W /, 1\
= 1 f U f ) ( 4 5 7 )
77?» l - - - I F 7 " ) " - -
The observed values o f P and i t s t h e o r e t i c a l curve o f (It.57) are
both shown i n Figs. It.38 and It.ho, as mentioned l a t e r .
I t i s necessary f o r the e x p l a n a t i o n o f the t u r b u l e n c e p r o d u c t i o n
+
mechanism t o i n v e s t i g a t e the d i s t r i b u t i o n o f i t s p r o d u c t i o n P EPv/U^
+
i n the w a l l r e g i o n . However, since P i n t h i s r e g i o n could not b * e
-52-
In d e t a i l , t h e t u r b u l e n t p r o d u c t i o n P, t h e t u r b u l e n t d i s s i p a t i o n e
and t h e d i r e c t - d i s s i p a t i o n E are n e a r l y equal t o each o t h e r a t t h e
+
sublayer edge y ~ 10 ( c f . 2.2.1 or F i g . 2.2). I n t h e sublayer
+
( y < _ 1 0 ) , they show E> e> P, i . e . t h e d e f i c i e n c y o f t u r h u l e n t energy.
+
In t h e b u f f e r l a y e r (10<_ y <_ 3 0 ) , they show P> e, i . e . t h e s u f f i c i e n c y
+
of t u r b u l e n t energy. And, i n t h e l o g r e g i o n ( y > 30, see F i g . 2.15),
they show P- E , i . e . t h e e q u i l i b r i u m s t a t e o f t u r b u l e n t energy.
z
Next, we consider t h e t u r b u l e n t energy d i f f u s i o n T = q v/2. r
r S
r i - l = U / 2 + „*) v
l
(4-59)
T /U*
T = f'(5f + l)/2 (4-61)
4.7.3 T u r b u l e n t e n e r g y budget
Since t h e t u r b u l e n t p r o d u c t i o n P, t h e d i s s i p a t i o n e and t h e d i f f u -
s i o n T= 9T /9y were o b t a i n e d i n t h e previous s e c t i o n , t h e r e l a t i o n s h i p
r
-54-
e x i s t s near t h e f r e e s u r f a c e . Between two ranges t h e r e e x i s t s
t h e i n t e r m e d i a t e range where t h e energy f l o w s under a dynamically
e q u i l i b r i u m s t a t e so t h a t P - e, as already shown i n F i g . 2.13.
By t h e way, from ( 2 . 2 8 ) , t h e t u r b u l e n t energy equation i n t h e
w a l l r e g i o n i s given by Pv _ tv Tv i?v a 1
( 7 \
+ ( 4 6 5 )
u; ~ u; u; 7JJ" ö«
y [ ï U J J *
+
Expanding t h e v e l o c i t y i n a Taylor s e r i e s about y = 0 y i e l d s
n
£-WWAW.\> +«>"> .
S l n c e + i
J(d* /a -y = o.3 and
y AdwVdyV =o.i a c c o r d i n g t o Hinze (1975),
+5
(h.6l) becomes J / 2 U * = 0.05 y + 2
2
+ o( y ) (4-68)
2
Consequently, t h e r a t e v = - 9 ( , ~ / 2 ^ ) / ö
7 y
+i
o f work by t h e viscous stress
becomes n e a r l y equal t o V =-0.1 v e r y close t o t h e w a l l .
T Thus,
Vm i n t h e w a l l r e g i o n must not be n e g l e c t e d because i t i s o f t h e same
order as t h e p r o d u c t i o n P +
or t h e d i s s i p a t i o n e , as seen i n F i g . +
U.36.
F i g . U.39 r e p l o t s t h e t u r b u l e n t energy budget i n t h e w a l l r e g i o n ,
where t h e values o f P, e, T and V T were o b t a i n e d by L a u f e r ( l 9 5 U ) .
P i s roughly balanced by e. The sum o f T and V T i s r o u g h l y balanced
by R'. I n d e t a i l , t h e sublayer ( y <10) +
i s t h e d e f i c i e n c y range
(P < e) o f t h e t u r b u l e n t energy. The b u f f e r l a y e r (10< y < +
30)
i s i t s excess range (P > e ) , and the l o g r e g i o n ( y >30) +
i s i t s equi-
l i b r i u m range (P - e ) , as mentioned p r e v i o u s l y . The excess t u r b u -
l e n t energy i n t h e b u f f e r l a y e r i s t r a n s p o r t e d i n t o t h e sublayer by
t h e pressure d i f f u s i o n (R' < 0), and a l s o i n t o t h e f r e e - s u r f a c e r e g i o n
by t h e k i n e t i c energy d i f f u s i o n (T < 0). Consequently, t h e t u r b u l e n t
budget among them balances i t s e l f .
Next, F i g . h.kO shows t h e t u r b u l e n t energy budget over t h e rough
bed. Although t h i s budget r e l a t i o n over t h e rough bed i s n e a r l y
s i m i l a r t o t h a t over t h e smooth bed, t h e roughness e f f e c t appears
i n t h e d i f f u s i o n terms T and R' near t h e w a l l . Though t h e t u r b u -
lence p r o d u c t i o n mechanism over t h e rough bed i s , a t p r e s e n t , s t i l l
unknown, i t s mechanism may be q u i t e d i f f e r e n t from t h a t over t h e
smooth bed, since t h e b u f f e r l a y e r ( t h e excess range of t h e t u r b u l e n t
energy) e x i s t s no l o n g e r over t h e rough bed. A clue t o make c l e a r
-55-
i t s mechanism w i l l be given i n Chapter 5.
L a s t l y , F i g . h.kl shows t h e r e l a t i o n s h i p between t h e p r o d u c t i o n P
and t h e d i s s i p a t i o n e. Although t h e r e are a l i t t l e s c a t t e r s i n t h e
observed values, t h e values o f |P-e|/e become n e a r l y equal t o 0.2 i n
the e q u i l i b r i u m r e g i o n , and then i n t h e f r e e - s u r f a c e region they
increase up t o 1 . 0 , i . e . P = 0. The curve (T) i n F i g . h.kl which
was c a l c u l a t e d from (U.56) and ( 4 . 5 7 ) , agrees w e l l w i t h t h e observed
values, and thus i t explains w e l l t h e r e l a t i o n s h i p between P and e.
References (Chapter 4)
14) Chen, C K. & Roberson, J. A. : Turbulence in
1) Raichlen, F. : Some turbulence measurements
in water, Proc. of ASCE, EM-2, pp. 73-97, 1967. wakes of roughness elements, Proc. of ASCE,
2) McQuivey, R. S. & Richardson, E. V. : Some H Y - l . p p . 5 3 - 6 7 , 1974.
turbulence measurements in open-channel flow, 15) Antonia, R. A. & Luxton, R. E . : The response
Proc. of ASCE, HY-1, pp. 209-223, 1969. of a turbulent boundary layer to a step change in
surface roughness, J. Fluid Mech., vol.48, pp. 721
3) Imamoto, H . : Hydraulic studies on turbu-
- 761, 1971.
lence in open-channel flow, Doctral thesis
presented to Kyoto University, 1971 (in Japa- 16) Engelund, F. : Dispersion of floating particles
nese). in uniform channel flow, Proc. of ASCE, HY-4,
4) Blinco, P. H. & Partheniades, E : Turbulence pp. 1149-1162,1969.
characteristics in free surface flows over smooth 17) Laufer, J. : Investigation of turbulent flow in a
and rough boundaries, J. Hydraulic Research, two-dimensional channel, NACA, TR-1053, 1951.
vol.9, pp.43-69, 1971. 18) Bremhorst, K.. & Walker, T. B. : Spectral measure-
ments of turbulent momentum transfer in fully
5) Nakagawa, H., Nezu, I. & Ueda, H. : Turbulence of
developed pipe flow, J. Ruid Mech., vol.61,
open channel flow over smooth and rough beds,
pp. 173-186, 1973.
Proc. of JSCE, No.241, pp. 155-168, 1975.
19) Hinze, J. O. : Turbulence (2-nd edi.), McGraw-
6) Townsend, A. A. : Equflibrium layers and wall
Hill, pp.586-770, 1975.
turbulence, J. Fluid Mech., vol.11, pp.97-120,
1961. 20) Corcos, G. M. : The structure of turbulent pressure
field in boundary layer flows, J. Fluid Mech.,
7) Luufer, J. : The structure of turbulence in fully'
vol.18, pp.353-378, 1964.
developed pipe flow, NACA, TR-1174, 1954.
21) Monin, A. S. & Yaglom, A. M. : Statistical fluid
8) Bradshaw, P. ; 'Inactive' motion and pressure
Mechanics, MIT Press, vol.2, pp.337-652, 1975.
fluctuations in turbulent boundary layer, J. Fluid
Mech., vol.30, pp.241-258, 1967. 22) Lin, C. C. : On Taylor's hypothesis and the
9) Nakagawa, H. & Nezu, I. : On a new eddy model acceleration terms in the Navier-stokes equation,
in turbulent shear flow, Proc. of JSCE, No.23I, Quar. of Applied Math., vol.x, pp.295-306,
pp.61-70, 1974. 1953.
0) Ljatkher, V. M. : Calculation of spectra of tur- 23) Bradshaw, P. : Conditions for the existence of
bulent pulsations in uniform flows, Proc. of 12-th an inertial subrange in turbulent flow, A.R.C.,
Cong, of IAHR, vol.2, B-3, 1967. R. & M., No.3603, 1967.
24) Grant, H. L . , Stewart, R. W. & MoiUiet, A. :
1) Smuiek, R. : Discussion on 'Measurement of
Turbulence spectra from a tidal channel, J. Fluid
turbulence in water' by Richardson et al., Proc
Mech., vol. 12, pp. 241-268, 1962.
of ASCE, HY-1, pp.519-523, 1969.
25) Townsend, A. A. : The structure of turbulent
2) Monin, A S. & Yaglom, A. M. : Statistical fluid
shear flow, Camb. Univ. Press, 1956.
mechanics, MIT Press, vol.1, pp.257-416, 1971.
26) Lawn, C. J. : The determination of the rate of
3) Grass, A. J. : Structural features of turbulent flow
dissipation in turbulent pipe flow, J. Fluid Mech.,
over smooth and rough boundaries, J. Fluid
vol.48, pp.477-505, 1971.
Mech., vol.50, pp. 233-255, 1971.
-56-
Chapter 5 Experimental consideration on mechanism o f
turbulence production
5.1 Introduction
The physics o f -wall-turbulence have now been e s t a b l i s h e d i n a t
l e a s t a q u a l i t a t i v e sense through experiments and r e f i n e d data analy-
sis, as p o i n t e d out by M o l l o - C h r i s t e n s e n ( l 9 7 1 ) • K l i n e e t al.(1967)
and Kim e t a l . ( l 9 7 l ) discovered t h e b u r s t i n g phenomenon, i n which
e j e c t i o n s played a predominant r o l e , by making v i s u a l s t u d i e s o f t h e
mechanism o f t u r b u l e n c e p r o d u c t i o n by t h e hydrogen-bubble technique.
A l s o , Corino & Brodkey(1969) and Grass(1971) found by f l o w v i s u a l i z a -
t i o n t h a t sweeps are prominent near t h e w a l l . This suggests t h a t
sweeps may be as important as e j e c t i o n s f o r t u r b u l e n c e p r o d u c t i o n .
This was v e r i f i e d by Nychas e t a l . ( l 9 7 3 ) and Offen & K l i n e ( 1 9 7 4 , - 7 5 ) •
On t h e basis o f these q u a l i t a t i v e r e s u l t s , some researchers
(Lu & W i l l m a r t h 1973, Brodkey e t a l . 1974,etc.) have attempted t o
o b t a i n more q u a n t i t a t i v e knowledge about t h e s t r u c t u r e o f t h e Reynolds
s t r e s s , i . e . turbulence p r o d u c t i o n , by u s i n g - p o i n t measurements.
In t h i s chapter, we a l s o e x p e r i m e n t a l l y i n v e s t i g a t e t h e mechanism
of turbulence p r o d u c t i o n i n open-channel flows by means o f b o t h c o n d i -
t i o n a l point-measurements and f l o w v i s u a l i z a t i o n , on t h e basis o f t h e
t h e o r e t i c a l p r e d i c t i o n s o f 2.6 and t h e previous r e s u l t s mentioned above
5.2 P r o b a b i l i t y d i s t r i b u t i o n s of t h e v e l o c i t y f l u c t u a t i o n s and
the instantaneous Reynolds s t r e s s * )
5.2.1 P r o b a b i l i t y d e n s i t y f u n c t i o n o f v e l o c i t y and
i t s h i g h e r moments
The data p l o t t e d i n F i g . 5-1 are t h e p r o b a b i l i t y d e n s i t i e s measu-
r e d i n each r e g i o n o f t u r b u l e n t f l o w over a smooth bed. Data were
sampled by d i v i d i n g t h e r e g i o n -3 < ü, v < 3 i n t o 40 c e l l s , where
ü=u/u' and v=v/v'. When t h e skewness f a c t o r S and t h e f l a t n e s s
f a c t o r F are g i v e n , t h e Gram-Charlier d i s t r i b u t i o n (2.96) i s d e t e r -
mined by using t h e f o l l o w i n g r e l a t i o n s from (2.9k):
-58-
value on t h e p o s i t i v e s i d e and i t s longer t a i l on t h e n e g a t i v e s i d e ,
and v i c e versa f o r p ( v ) .
The above d e s c r i p t i o n o f t h e f l o w over t h e smooth bed a p p l i e s
i n t h e case o f t h e rough beds t o o . The values o f t h e f o u r t h - and
f i f t h - o r d e r moments are shown i n F i g . 5-3, as w e l l as t h e data ob-
t a i n e d by Lawn(l97l). I n t h e w a l l and e q u i l i b r i u m r e g i o n s , t h e
f o u r t h - o r d e r moments h a r d l y d e v i a t e from t h e Gaussian d i s t r i b u t i o n
and thus t h e t h i r d - o r d e r d i s t r i b u t i o n may be v a l i d i n these regions.
I n t h e f r e e - s u r f a c e r e g i o n , however, t h e t h i r d - o r d e r approximation
may y i e l d s i g n i f i c a n t e r r o r s because o f a considerable d e v i a t i o n
from t h e Gaussian d i s t r i b u t i o n .
Judging from F i g . 5-3, t h e f i f t h - o r d e r moment w i d e l y deviates
=
from Gaussian v a l u e , Mso 0. The f o l l o w i n g can be deduced i n t h e
same way as (2.9k): <? = M -10M X (5-2)
The values o f M 50 c a l c u l a t e d from (5-2) by s e t t i n g Q 5 0 equal t o zero
are shown i n F i g . 5.3, and they agree f a i r l y w e l l w i t h t h e observed
values. This i s t h e reason why t h e h i g h e r - o r d e r cumulants can be
neglected as described i n 2.6. Since a l l t h e odd-order moments o f
a Gaussian d i s t r i b u t i o n are equal t o zero, i t i s suggested t h a t t h e
t h i r d - o r d e r Gram-Charlier d i s t r i b u t i o n should replace t h e Gaussian
one. A l s o , i t seems t h a t i n t h e w a l l r e g i o n t h e f i f t h - o r d e r moment
( i n g e n e r a l , odd-order moments) i s i n f l u e n c e d by t h e roughness,
which w i l l be discussed l a t e r i n d e t a i l .
-59-
between t h e t h e o r e t i c a l and experimental values was o b t a i n e d ,
e s p e c i a l l y i n t h e e q u i l i b r i u m r e g i o n , and any small discrepancy b e t -
ween t h e two was due t o a l a r g e d e v i a t i o n o f t h e f o u r t h - o r d e r moments
from Gaussian as shown i n F i g . 5.3- Consequently, t h e u n c o n d i t i o n a l
p r o b a b i l i t y d i s t r i b u t i o n p ^ w ) o f Reynolds s t r e s s can be represented
by one d i r e c t l y d e r i v e d from a Gaussian d i s t r i b u t i o n w i t h h i g h accu-
racy, as v e r i f i e d by Lu & ¥illmarth(l9T3) and Antonia & Atkinson(1973).
I t may be noted from F i g . 5.4 t h a t t h e p r o b a b i l i t y d i s t r i b u t i o n o f
values o f uv has a very sharp peak a t v=0 and a v e r y l o n g t a i l exten-
ding t o l a r g e values o f |w| . I n a t h e o r e t i c a l equation p ^ w ) becomes
i n f i n i t e a t w=0 because K (0)= 0
00
and becomes l a r g e r on t h e p o s i t i v e
side o f t h e curve than on t h e negative side because y j _ J ™ ^ (tcjdtossl
F M 16if
uv = f ( i 2 2 ZW» <* + 3i-^13 + 2é(M M 21 03 +M M )
12 30
2 2
+ 18(M + 2M^ + 2M\ ) - 30B M
2 2 22
The theoretical values of S uv and F uv were calculated from (5.3) and (5. h) by
using t h e experimental values o f t h e moments as i n d i c a t e d i n F i g . 5-5.
The c a l c u l a t e d values are i n comparatively good agreement w i t h t h e
observed ones. Some o f t h e discrepancy i n t h e F u v (8-th order mom-
ent o f u or v ) diagram may be due t o t h e e r r o r i n t r o d u c e d by neglec-
t i n g t h e h i g h e r - o r d e r cumulants.
As p r e v i o u s l y mentioned, t h e d i f f e r e n c e between t h e t h i r d - o r d e r
d i s t r i b u t i o n and t h e f o u r t h - o r d e r one i s n o t l a r g e . Both S u v and
-60-
F u v show a gradual v a r i a t i o n w i t h y/h up t o t h e e q u i l i b r i u m r e g i o n
(that i s , S - -2, Fuv u v = 10-20 i n t h i s r e g i o n ) , s i m i l a r t o t h e d i s -
t r i b u t i o n o f t h e h i g h e r - o r d e r moments shown i n F i g . 5-3, b u t i n t h e
f r e e - s u r f a c e r e g i o n they increase a b r u p t l y and then decrease towards
the surface. Although good agreement between t h e behaviour o f Suv
and Fuv i n open-channel f l o w and i n boundary-layer f l o w can be seen
up t o t h e e q u i l i b r i u m r e g i o n , some discrepancies s i m i l a r t o those i n
Fig. 5'2 are evident i n t h e f r e e - s u r f a c e r e g i o n corresponding t o t h e
outer l a y e r . Since t h e absolute values o f Suv and Fuv are compara-
t i t e l y l a r g e i n t h e f r e e - s u r f a c e r e g i o n , t h e s t r o n g asymmetry and
i n t e r m i t t e n c y i n t h e Reynolds-stress f l u c t u a t i o n s may appear, as
i n f e r r e d from F i g . 5-4.
The e f f e c t o f roughness on Suv and F v seems t o appear o n l y i n t h e
U
-61-
axis. Also, S u < 0, S v > 0 and D u < 0, D v > 0 except i n t h e imme-
diate v i c i n i t y of the w a l l . Another c h a r a c t e r i s t i c seen i n F i g . 5.6
i s t h e remarkable e f f e c t o f roughness i n t h e v i c i n i t y o f t h e w a l l .
I n case o f t h e smooth bed, S and D b o t h vary so g r a d u a l l y t h a t they
remain n e a r l y constant when y/h i s below t h e e q u i l i b r i u m r e g i o n .
For a rough bed, t h e absolute values o f S and D decrease towards t h e
w a l l so r a p i d l y t h a t they become zero a t y/h~0.1 and t h e n , changing
s i g n , increase upon approaching t h e w a l l . These trends are shown
by t h e dashed l i n e s i n F i g . 5.6. Since i n t h e middle o f t h e e q u i l i -
b r i u m r e g i o n t h e values o f S and D f o r a rough bed always c o i n c i d e
w i t h those f o r a smooth bed, i t can be concluded t h a t roughness has
an e f f e c t on t h e values o f S and D i n t h i s r e g i o n a t l e a s t .
Although t h e r e i s a c l o s e resemblance between our curves and Lawn's
f o r a smooth bed, t h e d i f f e r e n c e i n magnitude over t h e whole depth
between t h e two may be caused by t h e d i f f e r e n c e i n data p r o c e s s i n g .
As mentioned by F r e n k i e l & Klebanoff(1967), t h e analog method used by
Lawn cannot be expected t o o b t a i n such accurate values o f t h e h i g h e r -
order moments as t h e d i g i t a l procedure adopted here.
S , D , S~ and D~ can be e v a l u a t e d from (2.115) f o r t h e values o f
+ +
-62-
+ +
(4) Both S and D beoome nearly equal to zero 3 though the approvi-
mation of =0 is worse than that of Pd~.
+ +
Next, by s u b s t i t u t i n g t h e values o f S , D , S" and D~ i n t o (2.109)
f
-(2.112) v i t h t h e Bessel f u n c t i o n K ] / ( t ) equal t o (7r/2/)J e~ , t h e con-
2
-63-
5.3 I n t e r n a l s t r u c t u r e o f R e y n o l d s s t r e s s a n d
mechanism o f t u r b u l e n c e - p r o d u c t i o n * )
5.3.1 I n t e r n a l s t r u c t u r e o f Reynolds s t r e s s
As d e f i n e d i n 2 . 6 , t h e f r a c t i o n o f time Tj_(H) and t h e c o n t r i b u t i o n
to t h e Reynolds s t r e s s RSj_(H) corresponding t o each event i n t h e flow-
over a smooth bed are shown i n F i g s . 5.10 ( a ) , (b) and (c) f o r a t y -
p i c a l p o i n t i n t h e w a l l , e q u i l i b r i u m and f r e e - s u r f a c e regions r e s -
pectively. T h e o r e t i c a l curves were obtained from (2.1l6)-(2.119),
w h i l e our experimental data were analysed by almost t h e same method
of c o n d i t i o n a l sampling as t h a t used by Lu & W i l l m a r t h ( l 9 T 3 ) :
5 5)
T (n) = UmAjr^l (t.H)dt
t T t =-£-.2 {/,('.//)}> <'
T
RSXH)=±r lim ^rf "(') •(»)/,('.//)«" («-(') / , ( ' , / / ) ) , (5-6)
Lo : otherwise
I n t h e w a l l r e g i o n ( F i g . 5.10 a) and t h e e q u i l i b r i u m r e g i o n (
Fig. 5.10 b ) , t h e agreement between t h e experimental data and t h e
p r e d i c t e d values i s f a i r l y good over a wide range o f h o l e s i z e H,
so t h a t i t may be expected t h a t t h e t h i r d - o r d e r p r o b a b i l i t y distri-
b u t i o n represents t h e c o r r e c t p i c t u r e f o r a sequence o f b u r s t i n g
processes. I n t h e f r e e - s u r f a c e r e g i o n ( F i g . 5.10 c ) , however, a d i s -
crepancy between t h e two appears f o r l a r g e values o f H, owing t o t h e
neglect o f h i g h e r - o r d e r terms.
Although t h e time occupied f o r H < 1 amounts t o about a h a l f o f
the t o t a l t i m e , t h e corresponding c o n t r i b u t i o n t o t h e Reynolds s t r e s s
i s o n l y a few p e r c e n t , which suggests t h a t w ( t ) has a l a r g e i n t e r -
mittency. Though t h e inward i n t e r a c t i o n shows s l i g h t l y l a r g e r values
than t h e outward i n t e r a c t i o n , they both become n e g l i g i b l y small a t
H -5, and consequently t h e negative c o n t r i b u t i o n t o t h e Reynolds s t r e -
ss disappears. The sweep event decreases r a p i d l y w i t h H and beyond
-65-
t i c a l c o n s i d e r a t i o n s i n 2.6 may be able t o e x p l a i n t h e b u r s t i n g phe-
nomenon i n some d e t a i l f o r b o t h rough and smooth beds.
I n case o f t h e smooth bed, t h e r e e x i s t s t h e r e l a t i o n e j e c t i o n >
sveep > i n v a r d i n t e r a c t i o n > o u t v a r d i n t e r a c t i o n v i t h i n t h e observed
range o f y/h. The d i f f e r e n c e betveen t h e t v o i n t e r a c t i o n s , however,
+
i s almost n e g l i g i b l y s m a l l , because S ^ 0 on t h e basis o f t h e sym-
metry betveen S u and S v and between and D .v I n the wall region
the i n t e n s i t i e s o f a l l events decrease w i t h i n c r e a s i n g y/h, and good
agreement was o b t a i n e d between our r e s u l t s and those by Brodkey e t a l .
(19TM. Although t h e values o f t h e e j e c t i o n given by Lu & W i l l m a r t h
(1973) agree w e l l w i t h our data, t h e values o f t h e sweep show some
differences. This discrepancy w i l l be discussed l a t e r , w i t h t h e
+
characteristic q u a n t i t i e s i n t h e w a l l r e g i o n represented by a y
parameter.
Next, i n t h e e q u i l i b r i u m r e g i o n , t h e i n t e n s i t y o f each event i s
n e a r l y constant i r r e s p e c t i v e o f y/h. Since i n t h i s r e g i o n a dynamic
e q u i l i b r i u m e x i s t s between t h e t u r b u l e n c e p r o d u c t i o n and d i s s i p a t i o n ,
and a s i m i l a r i t y i n t u r b u l e n t s t r u c t u r e e x i s t s independently o f e x t e r -
n a l boundary c o n d i t i o n s ( c f . Chapter U), i t i s expected t h a t t h e b u r -
s t i n g process and t h e accompanying t u r b u l e n c e p r o d u c t i o n may a t t a i n a
s t a b l e e q u i l i b r i u m s t a t e i n t h i s r e g i o n , r e s u l t i n g i n almost constant
i n t e n s i t y f o r each event. The r a t e s o f i n t e n s i t y c o n t r i b u t e d by e j e -
c t i o n s and sweeps are about 75% and 60% r e s p e c t i v e l y , and t h e excess
Reynolds s t r e s s balances t h e sum o f t h e negative r a t e s c o n t r i b u t e d by
the inward and outward i n t e r a c t i o n s . S i m i l a r c h a r a c t e r i s t i c s were
v e r i f i e d by p o i n t measurements by Lu & W i l l m a r t h o r Brodkey e t a l . and
by flow v i s u a l i z a t i o n by Corino & Brodkey(1969) or Kim e t a l . ( l 9 7 l ) .
In the free-surface r e g i o n t h e r e l a t i v e i n t e n s i t y o f each event
r a p i d l y increases w i t h y/h. Near t h e f r e e surface b o t h ejections
and sweeps show a p o s i t i v e s t r e s s r a t e o f over 100*. w h i l e t h e n e g a t i -
ve stress brought about by t h e i n t e r a c t i o n s increases t o such an e x t e -
nt t h a t t h e d i f f e r e n c e s between t h e p o s i t i v e s t r e s s and t h e negative
s t r e s s are equal t o t h e n e t Reynolds s t r e s s , whose absolute v a l u e ,
-66-
however, becomes v e r y s m a l l , i . e . -üv^O.
The r e l a t i o n s between the i n t e n s i t i e s o f each event i n case o f the
rough beds and the smooth bed are almost the same, as shown i n F i g .
5.12 ( b ) , so t h a t we can confirm the o b s e r v a t i o n by Grass(1971) t h a t
both e j e c t i o n s and seeps e x i s t i r r e s p e c t i v e o f the roughness c o n d i t i -
ons. But, i n the range from the w a l l t o the middle of the e q u i l i b -
rium r e g i o n , where the roughness e f f e c t on the t u r b u l e n t structure
appears ( c f . Chapter k), some d i f f e r e n c e s i n the i n t e n s i t y p r o f i l e s
f o r the rough and smooth beds can be observed. I t i s noteworthy
t h a t , c o n t r a r y t o the case o f a smooth bed, the i n t e n s i t y of ejections
decreases towards the w a l l t o become n e a r l y equal t o t h a t o f the sweeps
at y / h - 0 . 1 , and sweeps may become more i n t e n s e than e j e c t i o n s i n the
v i c i n i t y of the w a l l as was observed by Grass(1971) u s i n g the hydrogen
bubble technique. This f a c t t h a t b o t h e j e c t i o n s and sweeps, which
are the predominant events i n the b u r s t i n g phenomenon, may be greatly
a f f e c t e d by the roughness c o n d i t i o n i s v e r y important and w i l l be
discussed i n d e t a i l l a t e r .
I n the same manner, F i g . 5-13 shows the c o n t r i b u t i o n s t o Reynolds
s t r e s s from the d i f f e r e n t events f o r the h i g h e r - v e l o c i t y group ex-
periments. I n s p i t e o f the wide v a r i a t i o n o f Fr, the values o f each
RS^ coincide f a i r l y v e i l w i t h each o t h e r , and also they agree w e l l
w i t h the t h e o r e t i c a l curves and the observed values of F i g . 5.12 (a).
Next, the f r a c t i o n s o f time occupied by each event on the smooth
and rough beds are shown i n Figs. 5.14 (a) and (b), respectively,
t o g e t h e r w i t h the t h e o r e t i c a l curves c a l c u l a t e d from (2.116).
U n l i k e RSj_, the magnitude o f the f r a c t i o n of t i m e T± s a t i s f i e s the
r e l a t i o n sweep > e j e c t i o n > o u t w a r d - i n t e r a c t i o n >_ i n w a r d - i n t e r a c t i o n .
Also, these r e s u l t s agree very w e l l w i t h those o b t a i n e d by Brodkey et
al.(197*0. Each f r a c t i o n of time i s n e a r l y constant i r r e s p e c t i v e o f
y/h except i n the f r e e - s u r f a c e region, that i s , T 2 or T 4 and Ti or T 3
-67-
as t h a t o f RS-^ t h a t t h e b u r s t i n g process near t h e f r e e surface may
consist o f smoother and more i s o t r o p i c events.
L a s t l y , t h e t u r b u l e n t energy p r o d u c t i o n P=-üv 8U/3y i s p l o t t e d vs.
the flow depth i n F i g . 5«15- The p o s i t i v e and negative time-average
production Pp and P , r e s p e c t i v e l y , are given by
n
Of course, P=Pp-P i s t h e n e t p r o d u c t i o n .
n The c o n t r i b u t i o n o f nega-
t i v e production t o t h e net p r o d u c t i o n i s comparatively s m a l l up t o
the e q u i l i b r i u m r e g i o n , b u t i n t h e f r e e - s u r f a c e r e g i o n i t becomes o f
the same order as t h a t o f t h e n e t p r o d u c t i o n though i t s absolute va-
lue i s very small. I t can be concluded, t h e r e f o r e , t h a t t h e energy
interchange from turbulence t o t h e mean f l o w should n o t be n e g l e c t e d
i n t h e free-surface region. This p r o p e r t y has been suggested by
Hino e t al.(1975) i n an o s c i l l a t i n g pipe f l o w , t o o .
5.3.3 Mechanism o f t u r b u l e n c e - p r o d u c t i o n i n t h e w a l l r e g i o n
and i t s p r e d i c t i o n
Since t h e t u r b u l e n c e occurs almost i n t h e e j e c t i o n and sweep events
and moreover most remarkably near t h e w a l l , we i n v e s t i g a t e t h e c o n t r i -
b u t i o n s o f t h e e j e c t i o n and sweep i n t h i s region.
Fig. 5-16 shows t h e d i s t r i b u t i o n s o f RSj_ and T^ f o r t h e e j e c t i o n
and sweep i n t h e w a l l r e g i o n , whose data were obtained i n t h e e x p e r i -
ment Case A - l (R*=600). A l s o , t h e r e s u l t s o b t a i n e d by Brodkey e t a l .
(R*=195) and Lu & Willmarth(R*=l800) are shown i n F i g . 5.16. Except
+
f o r t h e outer l a y e r ( y ~ R#), a good agreement among these data i s
recognized, and thus t h e f o l l o w i n g u n i v e r s a l c h a r a c t e r i s t i c s are
obtained. F i r s t l y , f o r t h e w a l l and e q u i l i b r i u m r e g i o n s , i . e .
30 <_ y +
<_ 0.6R#, RS^ and T^_ become n e a r l y constant. That i s , acco-
r d i n g t o our experimental data,
TÏ= 0.29 ± 0.012 , r„= 0.34 ± 0.008 1
f (5-9)
ftSj = 0.77 ± 0.029 , R$ = 0.57 ± 0.029 J
-68-
consequently t h e i r c o n t r i b u t i o n s may become constant i r r e s p e c t i v e o f
+
y. This e s s e n t i a l c h a r a c t e r i s t i c s should be taken i n t o account i n
m o d e l l i n g t h e t u r b u l e n c e , which w i l l be considered i n Chapter 6.
Next, i t should be n o t i c e d t h a t t h e r e l a t i o n i n magnitude between
+
the e j e c t i o n and sweep may reverse a t y =10, which corresponds t o
the edge o f t h e viscous sublayer. That i s , i n t h e sublayer t h e
sweeps may become more dominant than t h e e j e c t i o n s . Actually,
because t h e sublayer i s t h e range o f t h e d e f i c i e n c y of turbulent ener-
gy, as mentioned i n 4.7, t h e sweep motion ( u >0, v <0) may p l a y an
important r o l e i n t h i s energy t r a n s f e r .
W e l l , such an i m p o r t a n t f a c t t h a t t h e mechanism o f t u r b u l e n c e -
+
p r o d u c t i o n becomes d i f f e r e n t a t t h e edge o f sublayer, i . e . y =10,
can be w e l l e x p l a i n e d by t h e t h e o r y o f 2.6, as f o l l o w s : F i r s t l y , we
+
assume t h a t even i n t h e r e g i o n o f y <20,
(1) R becomes n e a r l y constant. A c t u a l l y , R =0.36 according t o
Eckelmann(1974).
(2) S +
=0 and S" =S V =-S u , j u d g i n g from F i g . 5.2.
(3) D +
=0 and D"/S~ -0.7, as i n f e r r e d from 5.2.3.
On t h e b a s i s o f t h e above assumption, RSi /RS and T4/T2 can be t 2
5.3.4 E f f e c t o f w a l l r o u g h n e s s upon b u r s t i n g e v e n t s
The r a t i o RS4/RS2 o f t h e Reynolds stress o f a sweep t o t h a t o f an
e j e c t i o n i s p l o t t e d i n F i g . 5.19, f o r each degree o f roughness, w i t h
the r e s u l t s o b t a i n e d by Wallace e t a l . ( l 9 7 2 ) , Lu & W i l l m a r t h ( 1 9 7 3 )
-69-
and Brodkey e t a l . ( l 9 7 U ) . I n our experiments, f o r y +
>100 and case
of a smooth bed RS4/RS2 remains n e a r l y constant, the Reynolds s t r e s s
of sweeps being about 10% o f t h a t o f the e j e c t i o n s , showing a good
agreement w i t h data by o t h e r i n v e s t i g a t o r s . For y +
<100, our values
increase towards the w a l l l i k e those o f Wallace e t a l . ( l 9 7 2 ) o r Brodkey
et a l . (197*0, w h i l e the values by Lu & Willmarth(1973) show the reverse
tendency, decreasing towards the w a l l . Lu & W i l l m a r t h suggested
t h a t t h i s discrepancy might be due t o the d i f f e r e n c e i n the Reynolds
number. However, from F i g . 5.17 i t can be seen t h a t RS /RS may 4 2
-70-
W < 7 (5-io)
5.3.5 B e h a v i o u r o f v e l o c i t y f l u c t u a t i o n s when t u r b u l e n c e o c c u r s
I t i s an e s s e n t i a l work how t o connect t h e knowledge o f b u r s t i n g
events obtained by point-measurements w i t h t h a t o b t a i n e d by v i s u a l
methods. Since t h e l e v e l H i s expected t o f u l f i l l t h e f u n c t i o n o f
a f i l t e r by which t h e e j e c t i o n motion can be d e t e c t e d from t h e e j e c -
t i o n event s i g n a l s , as p o i n t e d o u t i n 2 . 6 , we here examine t h i s sugg-
estion i n the following.
An example o f t h e c o n d i t i o n a l p r o b a b i l i t y d i s t r i b u t i o n s o f u, i . e .
p (u
e |w>H) and p ( u |w>H) c a l c u l a t e d from (2.125) and (2.126) as w e l l
s
-71-
t i v e l y , j u s t a t t a i n s t h e l e v e l H, t o g e t h e r v i t h t h e t h e o r e t i c a l cur-
ves c a l c u l a t e d from (2.121) and (2.122). Although t h e r e i s a l i t t l e
d i f f e r e n c e between t h e observed and t h e o r e t i c a l values o v i n g t o using
the d i g i t a l s i g n a l s o f v ( t ) , agreement betveen t h e t v o i s comparative-
l y good. Of course, u <0 and u > 0 , and t h e l a r g e r H i s , t h e more
e s
i n t e n s i v e t h e v e l o c i t y f l u c t u a t i o n s accompanied v i t h t u r b u l e n c e - p r o d u -
;
c t i o n become. At H= 5 vhere t h e i n t e r a c t i o n events almost disappear,
Fig. 5.22 i n d i c a t e s t h a t u =-1.0~ -1.5
e and u = 1.0- 1.5•
s
-72-
5.4 P e r i o d i c c h a r a c t e r i s t i c s o f t h e b u r s t i n g phenomenon*)
5.4.1 B r i e f r e c a p i t u l a t i o n of the bursting period obtained previously
K l i n e e t a l . ( l 9 6 7 ) , Corino & Brodkey(1969), Kim e t a l . ( l 9 7 l ) , Nychas
et a l . ( l 9 7 3 ) , Offen & K l i n e ( l 9 7 4 , -75) and others found by f l o w v i s u a l -
i z a t i o n t h a t a sequence o f t h e b u r s t i n g events had a q u a s i - c y c l i c p r o -
cess. That i s t o say, i t shows a p e r i o d i c motion on t h e average i n
space and t i m e , b u t n o t p e r f e c t l y p e r i o d i c a t one place i n time nor a t
one time i n space. For example, Corino & Brodkey p o i n t e d out i n F i g .
2.20 t h a t there were v a r i a t i o n s o f t h e sequence, and t h a t a l l o f t h e
steps d i d not appear a l l t h e time or i n t h e exact fashions described,
but on t h e average i t proceeded.
Thus, i n order t o r e v e a l t h e governing parameters o f t h e b u r s t i n g
phenomenon, i t i s necessary t o i n v e s t i g a t e i t s p e r i o d i c c h a r a c t e r i s t i -
cs. The v i s u a l method whereby t h e b u r s t i n g p e r i o d i s determined by
counting t h e frequency o f occurrences o f t h e coherent motions near t h e
w a l l v i s u a l i z e d on high-speed movie f i l m s , i s simple and p l a i n .
However, i t needs l a b o r i o u s work and more or l e s s i n v o l v e s s u b j e c t i v e
judgement. On t h e o t h e r hand, some researchers as w e l l as the author
have shown t h a t t h e e x i s t e n c e o f a sequence o f b u r s t i n g events such as
e j e c t i o n s , sweeps and i n t e r a c t i o n s can be also d e t e c t a b l e even i n t h e
f l u c t u a t i n g v e l o c i t y s i g n a l s obtained by h o t - w i r e or h o t - h i l m anemo-
meters when a c o n d i t i o n a l sampling technique i s reasonably used.
I f i t i s , t h e r e f o r e , p o s s i b l e t o e s t a b l i s h a reasonable c r i t e r i o n f o r
d i s c r i m i n a t i n g t h e e j e c t i o n s or t h e sweeps from these v e l o c i t y s i g n a l s ,
the b u r s t i n g p e r i o d can be a l s o evaluated e a s i l y from t h e data a n a l y s i s
of point-measurement s i g n a l s by using a high-speed d i g i t a l computer.
Of course, these r e s u l t s should be compared w i t h t h e v i s u a l data, since
i t may be f a i r l y difficult t o detect a c c u r a t e l y t h e coherent motions
extending i n space only by one or a few h o t - f i l m s .
F i r s t l y , Rao e t a l . ( l 9 7 l ) estimated t h e mean b u r s t i n g p e r i o d T B
from t h e s i n g l e - h o t - w i r e s i g n a l s o f u ( t ) i n a boundary l a y e r by u s i n g
-73-
a special c r i t e r i o n f o r discrimination. They proposed t h e f o l l o w i n g
experimental r e s u l t s .
T t / , V = 0.65 Rf
B (5-12) a n d T U /i =Z2
0 max t (5-13)
where 6 i s t h e boundary l a y e r t h i c k n e s s .
Consequently, (5-12), (5-13) and (5.1*0 are almost t h e same. (
T
B max/ "becomes n e a r l y equal t o k from (5-13) since 5/6 f o r t h e
u 5
#
-75-
and the w a l l roughness upon the b u r s t i n g p e r i o d s y s t e m a t i c a l l y .
Now, i t may be n o t i c e d i n F i g s . 5.27 and 5.28 that T 2 shows a
n e a r l y l i n e a r increase w i t h H when H<5, and a more remarkable increase
when H>5 although the slope o f T (H) 2 increases more or l e s s c o n t i n u -
ously. This tendency o f T (H) 2 may be r e l a t e d t o the f a c t t h a t the
i n t e r a c t i o n events scarcely c o n t r i b u t e t o the p r o d u c t i o n of Reynolds
stress when l e v e l H reaches about 5- This suggests t h a t the e j e c t i o n
s i g n a l s w ( t ) or the sweep s i g n a l s w i j ( t ) w i t h H<_5
2 contain a part , of
the i n t e r a c t i o n motions corresponding t o the events (l)-(4) i n F i g .
2.20. This suggestion might be also i n f e r r e d from t h e r e s u l t s o f a
c o n d i t i o n a l sampling technique, as already shown i n 5.3.5.
From the above phenomenological c o n s i d e r a t i o n s , t h e c o n t r i b u t i o n s
of the i n t e r a c t i o n motions should be removed from t h e sequence o f t h e
b u r s t i n g process i n order t o evaluate the p e r i o d o f only t h e e j e c t i -
ons or sweeps. However, t h e r e i s a t present q u i t e a l a c k o f know-
ledge about the c o n t r i b u t i o n s of the i n t e r a c t i o n - l i k e motions which
may be contained i n the s i g n a l s o f w ( t ) or wi»(t).
2 Hence, we now
propose a t e n t a t i v e assumption t h a t the i n t e r a c t i o n - l i k e c o n t r i b u t i o n s
may be o f the same order as those o f the i n t e r a c t i o n event s i g n a l s
w j ( t ) or w ( t ) , since the c o n t r i b u t i o n s o f |w^(t) |<_5
3 might be roughly
equal t o each o t h e r , owing t o a d e t e c t i o n of old-born or new-born
small b u r s t i n g motions. By assuming t h a t the number Nj_ o f n occurren-
ces of these i n t e r a c t i o n - l i k e motions w i t h a l e v e l o f H i s roughly
given by an average of those o f two i n t e r a c t i o n events: t h a t i s ,
Nj_ =(N +N )/2, the r e v i s e d b u r s t i n g p e r i o d i s d e f i n e d as f o l l o w s :
n 1 3
T=
2 T/(N -N ) 2 in for ejections, T^E T/(N -Ni )
4 n f o r sweeps.
T U
2 max / h i s shown against H i n F i g . 5.29 (.(T -Tj_ ) denotes T )
e n 2 for
the case corresponding t o F i g . 5.27. When H i s s m a l l , T2U /h i s max
-76-
as t h e minimum value o f T ( H ) and Ti+(H), r e s p e c t i v e l y , because i t
2
VC/mox , ^ m n x __ ( T _ . 5 _ 3 . f j ) (5-15)
h n
(5.15) shows t h e same order as ( 5 . l U ) or Lu & W i l l m a r t h ' s results i n
a boundary l a y e r , though a q u a n t i t a t i v e comparison among these data
cannot be done reasonably because o f t h e d i f f e r e n c e s o f t h e f l o w con-
d i t i o n s and t h e e v a l u a t i o n methods. I t should be n o t i c e d t h a t t h e
ejection period T e becomes n e a r l y equal t o t h e sweep p e r i o d T , as s
-TT-
f i e d v i t h t h e e j e c t i o n or sveep one, T ~ T ~ T-g . e s
the f o l l o v i n g e x p l a n a t i o n s . That i s , T= T - T = T . e s B
-78-
p r o b a b i l i t y o f t h e b u r s t i n g p e r i o d can be described approximately
by a log-normal d i s t r i b u t i o n , i . e . (5.16), i r r e s p e c t i v e o f Re, Fr
and t h e w a l l roughness.
Then, t h e f o l l o w i n g equation can be o b t a i n e d e a s i l y from (5.16):
-'fi
1>/2
Denoting T =T,
B C-pE( (T-T) 2
) and I|>EO"B/TT
(5-18)
5.5 I n t e r n a l s t r u c t u r e o f b u r s t i n g phenomenon
5.5.1 Dependence o f t h e b u r s t i n g phenomenon upon t h e i n n e r and o u t e r
parameters
For t h e p r e s e n t , from the previous and present experimental data,
we can approve t h e o p i n i o n t h a t t h e b u r s t i n g p e r i o d may be c o n t r o l l e d
by the outer r a t h e r than t h e i n n e r parameters. Consequently, t h e
t u r b u l e n t s t r u c t u r e i n t h e w a l l r e g i o n o f an open-channel f l o w may be
c h a r a c t e r i z e d not o n l y by the i n n e r parameter (see Chapters 2 and 4 ) ,
but a l s o by t h e outer parameter t h a t c h a r a c t e r i z e s t h e f r e e - s u r f a c e
r e g i o n s , which have been a l r e a d y recognized i n a boundary l a y e r f l o w .
-79-
Nov, t h e streamvise, the v e r t i c a l and t h e t r a n s v e r s e s p a t i a l
scales o f t h i s coherent motion are denoted by X , \ l 2 and X , respec-
3
-80-
the boundary l a y e r edge, b u t we cannot f i n d any bulge motions i n an
open-channel f l o w . We i n f e r here t h a t t h e r e may e x i s t a close r e l a -
t i o n between t h e b u r s t i n g mechanism and t h e breakdown mechanism o f
the l a r g e s t - s c a l e eddy w i t h scale Lo. I f the Strouhal s i m i l a r i t y
might be v a l i d f o r b o t h mechanisms, t h e f o l l o w i n g r e l a t i o n c o u l d be
obtained. T u /h^{Si Y
B max K
l
U J h ) (5-20)
5.5.2 F u r t h e r s u g g e s t i o n s on b u r s t i n g phenomenon
I t may be f a i r l y w e l l a s s e r t e d from t h e v i s u a l and h o t - f i l m measure-
ments t h a t t h e breakdown o f t h e l a r g e s t - s c a l e eddy i s i n t e r m i t t e n t ,
and t h e n generates most o f t h e t u r b u l e n t energy or t h e Reynolds stress
(cf. F i g s . 2.5 & 2.9). On t h e o t h e r hand, as i n d i c a t e d by Sandborn
(1959)5 t h e breakdown o f t h e s m a l l e s t - s c a l e eddy i s also i n t e r m i t t e n t ,
and t h e n d i s s i p a t e s most o f t h e t u r b u l e n t energy.
Among t h e e v a l u a t i o n methods o f t h e b u r s t i n g p e r i o d mentioned
p r e v i o u s l y , t h e methods o f Rao et a l . ( l 9 7 l ) and Ueda & H i n z e ( l 9 7 5 )
are based on a m i c r o s c a l e i n t e r m i t t e n c y , w h i l e t h e methods o f Kim e t
al.(l97l), Lu & W i l l m a r t h ( 1 9 7 3 ) , t h e present author as w e l l as t h e
v i s u a l method are based on a macroscale i n t e r m i t t e n c y . Since almost
the same c h a r a c t e r i s t i c s o f t h e b u r s t i n g p e r i o d have been e v a l u a t e d
from these two d i f f e r e n t methods, i t i s i n f e r r e d t h a t t h e breakdown
of a macroscale eddy may have a close r e l a t i o n and i n t e r a c t i o n w i t h
t h a t o f a m i c r o s c a l e eddy, and t h a t t h e t u r b u l e n t s t r u c t u r e may be
essentially c h a r a c t e r i z e d by b o t h t h e macro- and m i c r o - s c a l e s .
Consequently, an energy cascade process i n which t h e t u r b u l e n t energy
of macroscale eddies i s g r a d u a l l y t r a n s f e r e d t o t h a t o f m i c r o s c a l e
eddies, would be more c o m p l i c a t e d than we have p r e v i o u s l y understood.
Now, assuming t h a t t h e r e e x i s t s a s e l f - s i m i l a r i t y i n t h e breakdown
of eddies o f t h i s cascade process, log-normal c h a r a c t e r i s t i c s o f t h e
bursting process or t h e coherent motion r e p o r t e d here c o u l d be
explained by t h e eddy model o f Gurvich & Yaglom(l96T)• Thus, i t i s
suggested t h a t t h e b u r s t i n g phenomenon may be c l o s e l y r e l a t e d t o
the breakdown mechanism o f both t h e macro- and m i c r o - s c a l e eddies.
By t h e way, some researchers have a l r e a d y found t h a t t h e r e e x i s t e d
also the h i g h and low speed s t r e s k s , whose spacing i s n e a r l y equal t o
2h, on t h e f r e e surface of an a c t u a l r i v e r , t h a t i s t o say, t h e b o i l
phenomenon. Therefore, although t h e scale between t h e b u r s t i n g and
the b o i l phenomena i s q u i t e d i f f e r e n t when t h e Reynolds number i s very
l a r g e , t h e both phenomena are very s i m i l a r t o one another. Hence, we
may be able t o expect t h e p o s i b i l i t y o f an i n t e r r e l a t i o n between t h e
b o t h phenomena, and furthermore we are now i n v e s t i g a t i n g t h e s p a t i a l
s t r u c t u r e o f b u r s t i n g phenomenon by o b t a i n i n g t h e c o n d i t i o n a l s p a t i a l
c o r r e l a t i o n s i n order t o make c l e a r t h e above suggestions.
5.6 V i s u a l o b s e r v a t i o n s o f b u r s t i n g phenomenon * )
5.6.1 I n s t a n t a n e o u s v e l o c i t y p r o f i l e s i n v e r t i c a l p l a n e
The v e l o c i t y p r o f i l e s o f t h e l o w e r - v e l o c i t y group i n Table U . l ( c )
were v i s u a l i z e d by t h e hydrogen-bubble method i n order t o o b t a i n t h e
distinct pictures. The t i m e - l i n e s o f hydrogen bubbles were photo-
graphed by t h e l6mm high-speed cine-camera. F i g . 5-39 shows some
t y p i c a l photographic i l l u s t r a t i o n s o f t h e e j e c t i o n and sweep motions
f o r t h e rough bed (Case D - l ) . The p i c t u r e s o f No.l - k i n d i c a t e t h e
e j e c t i o n motions, w h i l e t h e p i c t u r e s o f No.5 - 9 i n d i c a t e t h e sweep
motions. For example, t h e p i c t u r e o f No.l probably shows t h e f i r s t
stage o f the b u r s t i n g , where a low-speed s t r e a k i s l i f t e d from t h e
wall. The p i c t u r e o f No.2 probably shows t h e o s c i l l a t o r y motion i n
the streamwise d i r e c t i o n , and then t h e v i o l e n t b u r s t o r breakup occurs
i n No.3 p i c t u r e . These p i c t u r e s o f t h r e e stages o f t h e b u r s t i n g
motion agree very w e l l w i t h Kim e t a l ' s o b s e r v a t i o n s . Also, the
e j e c t i o n motion may occur more v i o l e n t l y than sweep motion does.
By means o f t h e method described i n 3.4.1, t h e instantaneous velo-
-82-
c i t y p r o f i l e s (ü, v) vere analyzed, and some examples o f these data
f o r t h e smooth bed are shovn i n F i g . 5.1*0, t o g e t h e r v i t h the mean
v e l o c i t y d i s t r i b u t i o n U(y) described by a t h i c k s o l i d l i n e . Since
according t o (5.15) the b u r s t i n g p e r i o d i n t h i s case i s n e a r l y equal
to T = ( l . 5~3.0)h/U = ( l ~ 1 . 5 ) s e c , t h e e i g h t successive p i c t u r e s i n
B max
F i g . 5.40 (a) & (b) are expected t o show one sequence o f the b u r s t i n g
process. The lov-speed p a r t (u=u-U <0) tends t o be s u r e l y l i f t e d up
and e j e c t e d from t h e w a l l .
I n order t o examine t h e behaviours o f t h e low-speed and t h e h i g h -
speed p a r t s , we have performed the c o n d i t i o n a l l y ensemble average o f
the instantaneous v e l o c i t y ( u , v ) , i n t h e same manner as Grass(1971)
did. The t o t a l successive p i c t u r e s o f l 6 o frames were d i v i d e d i n t o
8 groups; t h a t i s , one group c o n s i s t e d o f 20 successive frames, whose
elapsed time T G was 3-7 seconds ( c f . F i g . 5.40). Thus, i t i s expec-
t e d t h a t each group contains at l e a s t one b u r s t i n g process or t h e
coherent motion, since t h e value o f T G i s l a r g e r than t h e mean b u r s -
t i n g p e r i o d T-g. The frames o f p i c t u r e where t h e v e l o c i t y u i n d i c a t e d
the maximum or t h e minimum at any measured p o i n t Pj_ ( i = l - 12) were
chosen from each group, and then these instantaneous v e l o c i t y p r o f i l e s
were ensemble averaged. Figs. 5.41 and 5.42 show the c o n d i t i o n a l l y
averaged instantaneous v e l o c i t y p r o f i l e s a t i n s t a n t s when u became
minimum and maximum a t any p o i n t Pj_, r e s p e c t i v e l y . As seen i n F i g .
5.41, t h e low-speed p a r t drew b l a c k has a s t r o n g c o r r e l a t i o n w i t h
v >0, and s u r e l y i t shows the e j e c t i o n motion. I n t h e same manner,
F i g . 5-42 shows t h e sweep motion. The both f i g u r e s i n d i c a t e t h a t t h e
l i f t u p - e j e c t i o n motions may e x i s t up t o t h e f u r t h e r r e g i o n from t h e
w a l l , compared w i t h t h e sweep motions. This suggests a p o s i b i l i t y o f
the i n t e r r e l a t i o n between the e j e c t i o n ( b u r s t i n g ) and the b o i l pheno-
mena, as mentioned previously.
+
=R /3=100, t o o .
W However, as y increases, the streaks i n t e n d t o be-
come weaker and also i t s spacing may become l a r g e r (see our paper des-
c r i b e d i n the margin o f p. 82). At present, i t i s q u i t e unknown
whether t h i s i s due t o the f a c t t h a t the b u r s t i n g phenomenon may become
weaker near the f r e e surface o r t h i s may i n d i c a t e an i n t e r r e l a t i o n
w i t h the b o i l phenomenon. Thus, we are now performing the f u r t h e r
d e t a i l e d researches about these t h i n g s .
References (Chapter 5 )
4
1) Davies, P. O. A. L. & Yule, A. J. : Coherent ) Willmarth, W. W. & Lu.S. S. : Structure of the
structures in turbulence, J. Fluid Mech., vol.69, Reynolds stress near the wall, J. Fluid Mech.,
vol
pp.513-539, 1975. - 5 5 , pp.65-92, 1972.
5
2) MoUo-Christensen, E. : Physics of turbulent now, ) Wallace, J. M., Eckelmann, H. & Brodkey, R. S. :
w
A1AAJ., vol.9, pp. 1217-1228, 1971. ^ - l region in the turbulent shear flow, J.
3) Kline, S. J., Reynolds, W. C , Schraub, F . A. & Fluid Mech., vol.54, pp.39-48, 1972.
6 )
Runstadler, P. W. : The structure of turbulent Nychas, S. G., Hershey, H. C. & Brodkey, R. S. :
A
boundary layers, J. Fluid Mech., vol.30, pp.741- visual study of turbulent shear flow, J . Fluid
7 7 3 , 1 9 6 7
' Mech., vol.61, pp.513-540, 1973.
-84-
7) Offen, G. R. & Kline, S. J. : Combined dye- 19) Zaric, Z. : Analyse statistique de la turbulence
streak and hydrogen-bubble visual observation of pre" d'une paroi par e'chantillonnage .-onditionnel,
a turbulent boundary layer, J. Fluid Mech., Comptes Rendus, t-269, Se'ric-A, pp.513-515,
vol.62, pp. 223-239, 1974. 1972.
8) Lu, S. S. & Willmarth, W. W. : Measurements of 20) Echelmann, H. : The structure of the viscous
the structure of the Reynolds stress in a turbulent sublayer and the adjacent wall region in a turbulent
boundary layer, J. Fluid Mech., vol.60, pp.481 - channel flow, J. Fluid Mech., vol.65, pp.439¬
511, 1973. 459, 1974.
9) Brodkey, R. S., Wallace, J. M. & Eckelmann, H. : 21) Rao, K. N., Narasimha, R. & Narayanan, M. A. B. :
Some properties of truncated turbulence signals The 'bursting' phenomenon in a turbulent bound-
in bounded shear flows, J. Fluid Mech., vol.63, ary layers, J. Fluid Mech., vol.48, pp.339-352,
pp. 209-224, 1974. 1971.
10) Nakagawa, H. & Nezu, I. : Prediction of the 22) Lu.S. S. & Willmarth, W. W. : Measurement of
contributions to the Reynolds stress from the the mean period between bursts, Phys. of Fluids,
bursting events in open-channel flows, J. Fluid vol.16, pp.2012-2013, 1973.
Mech., vol.80, pp.99-1 28, 1977. 23) Laufer, J. & Narayanan, M. A. B. : Mean period
11) Frenkiel, F. N. & Klebanoff, P. S. : Probability of the turbulent production mechanism in a
distributions and correlations in a turbulent boundary layer, Phys. of Fluids, vol. 14, pp. 182¬
boundary layer, Phys. of Fluids, vol. 16, pp. 726¬ 183, 1971.
737, 1973. 24) Hinze, J. 0. : Turbulence (2-nd edi.), McGraw-
12) Bayazit, M. : Free surface flow in a channel of Hill, pp.586-770, 1975.
large relative roughness, J. Hydraulic Research, 25) Clark, J. A. & Markland, E. : Flow visualization
vol. 14, pp. 115-126, 1976. in turbulent boundary layers, Proc. of ASCE,
13) Gupta, A. K. & Kaplan, R. E . : Statistical charac- HY-10, pp.1653-1664, 1971.
teristics of Reynolds stress in a turbulent boundary 26) Gupta, A. K., Laufer, J. & Kaplan, R. E. : Spatial
layer, Phys. of Fluids, vol. 15, pp. 981-985, 1972. structure in the viscous sublayer, J. Fluid Mech.,
14) Antonia, R. A. & Atkinson, J. D. : Higher-order vol.50, pp.493-512, 1971.
moments of Reynolds shear stress fluctuations in 27) Sandborn, V. A. : Measurements of intermittency
a turbulent boundary layer, J. Fluid Mech., of turbulent motion in a boundary layer, J. Fluid
vol.58, pp.581-593, 1973. Mech., vol.6, pp. 221-240, 1959.
15) Lawn, C. J. : The determination of the rate of 28) Ueda, H. & Hinze, J. 0. : Fine-structure turbu-
dissipation in turbulent pipe flow, J. Fluid Mech., lence in the wall region of a turbulent boundary
vol.48, pp.477-505, 1971. layer, J. Fluid Mech.,. vol.67, pp. 125-143, 1975.
16) Corino, E. R. & Brodkey, R. S. : A visual investi- 29) Mollo-Christensen, E. : Intermittency in large-
gation of the wall region in turbulent flow, J. scale turbulent flows, Annu. Review of Fluid
Fluid Mech., vol.37, pp. 1-30, 1969. Mech., vol.5, pp. 101-1 18, 1973.
17) Kim, H. T , Kline, S. J. & Reynolds, W. C. : The 30) Gurvich, A. S. & Yaglom, A. M. : Breakdown of
production of turbulence near a smooth wall in eddies and probability distribution for small-
a turbulent boundary layer, J. Fluid Mech., scale turbulence, Phys. of Fluids, vol.10, pp.
vol.50, pp. 133-160, 1971. S59-S65, 1967.
18) Grass, A. J. : Structural features of turbulent flow 31) Einstein, H. A. & Li, H. : The viscous sublayer
over smooth and rough boundaries, J. Fluid along a smooth boundary, Proc. of ASCE, EM-2,
Mech., vol.50, pp.233-255, 1971. pp. 1-27, 1956.
-85-
Chapter 6 P h y s i c a l model o f w a l l turbulence
6.1 Introduction
Even before t h e existence c f t h e b u r s t i n g phenomenon vas 'Ï.J**.-?*r*-»"»*'r«3 r,
*) Refer t o our paper {On a new eddy model in turbulent shear flow)
p u b l i s h e d i n Proc. of JSCE, No. 231, pp.61-70, 1974.
-87-
flow. I t s legs AB and CD are i n a plane p a r a l l e l t o t h e x-y p l a n e ,
and i t s t o p BC i s p a r a l l e l t o t h e z-axis. The c o n d i t i o n 9=0 may
represent an i n c i p i e n t stage o f t h e Il-eddy which coincides w i t h a
r
liïirji;; eddy or an attached eddy on t h s ".. ali. However, since t h e
i n c i p i e n t stage o f t h e Il-eddy or t h e horseshoe v e r t e x i s very com-
p l i c a t e d because o f t h e n o n - e q u i l i b r i u m s t a t e ( c f . F i g . 6.7), an
e q u i l i b r i u m c o n d i t i o n under which t h e H-eddy i s l i f t e d up from t h e
bottom and f u l l y developed i s considered here. That i s t o say, some
of t h e horseshoe v o r t i c e s i n t h e w a l l r e g i o n s u r v i s e w i t h o u t t h e v o r -
+
t e x breakdown during t h e i r development, even i n t h e region o f y > 1 0 0 ,
and then they a t t a i n t h e e q u i l i b r i u m s t a t e i n t h e t u r b u l e n t structure
, as mentioned i n Chapters h and 5. Since AB »BC because o f BC* -
At - 100 ( c f . Fig. 6.6), t h e c o n t r i b u t i o n o f t h e v o r t e x BC t o t h e
t u r b u l e n t s t r u c t u r e can be ignored except t h a t t h e v o r t e x BC w i l l
suffer the l i f t force. The v o r t e x tubes AB and CD may be reasonably
assumed t o have an e l l i p t i c a l cross s e c t i o n w i t h a long radius a i n
the x - d i r e c t i o n and a short radius b i n t h e z - d i r e c t i o n , since t h e y
are d i s t o r t e d by t h e v o r t e x - s t r e t c h i n g e f f e c t . Consequently, i t can
be assumed t h a t t h e primary motion o f t h e Il-eddy describes an e l l i -
p t i c steady c i r c u l a t i o n w i t h angular v e l o c i t y w i n t h e v o r t e x axis.
From t h e above assumptions, t h e motion o f any Il-eddy i s f o r m u l a t e d
z )
0 i s probably much smaller than t h a t o f t h e e l l i p t i c motion,
since t h e Il-eddy a t t a i n s t h e e q u i l i b r i u m s t a t e . That i s t o say,
the Il-eddy c o n s i s t s o f a primary motion w i t h an e l l i p t i c steady
c i r c u l a t i o n i n t h e v o r t e x a x i s , and a secondary motion, i . e . t h e
o v e r a l l deformation motion o f t h e eddy, w i t h a p e r t u r b a t i o n o f i t s
vortex tube due t o t h e v o r t e x - s t r e t c h i n g e f f e c t .
By t h e way, since t h i s secondary motion d i s t o r t s t h e l e n g t h I (=
AB=CD) o f t h e legs and t h e i n c l i n a t i o n angle 9, i t i s very difficult
to determine t h e v e l o c i t y U 0 o f t h e secondary motion a c c u r a t e l y .
So, t h e f o l l o w i n g s i m p l i f i e d f o r m u l a t i o n i s considered here. The
Il-eddy i s s u r e l y i n f l u e n c e d by t h e v o r t e x - s t r e t c h i n g since t h e r e i s
the v e l o c i t y shear 3U/3y=3U/3y even though t h e r e l a t i v e v e l o c i t y
U EU-U i s very s m a l l .
C Then, t h e e q u i l i b r i u m Il-eddy ABCD is infini-
t e s i m a l l y perturbed t o AB'C'D by t h i s v o r t e x - s t r e t c h i n g e f f e c t , as
shown i n F i g . 6.9 (consequently, z =0).
0 The head BC s u f f e r s t h e
l i f t f o r c e L and the drag f o r c e D. So long as t h e IT-eddy under t h e
e q u i l i b r i u m c o n d i t i o n does not instantaneously disappear by t h e
s t r e t c h i n g , some apparent resistances F must be thought t o work upon
AB and CD o f t h e eddy against t h e forces o f L and D, which probably
r e s u l t from complex i n t e r a c t i o n s between t h e mean f l o w and t u r b u l e n c e .
Now, d e s c r i b i n g t h e v o r t e x tube BC as a c i r c u l a r c y l i n d e r w i t h a
radius R and a l e n g t h li , t h e f o l l o w i n g equation can be o b t a i n e d
from F i g . 6.9. - Re
FI7!
2
] AFCOB6 + Fsine. jO) +ÓD=O ]
ó'ê) ~ \c DKÓB)
dy
Consequently, t h e secondary motion Uo=(Ax,Ay,0) can be p r e d i c t e d
from ( 6 . 5 ) . Considering t h a t t h e apparent r e s i s t i n g f o r c e i s i n c l u -
ded o n l y i n t h e term A, i t i s suggested t h a t a p e r t u r b a t i o n along a
v o r t e x l i n e may be represented by At. So, 6 i s assumed t o be i n -
dependent o f time t f o r s i m p l i f y i n g t h e f o l l o w i n g a n a l y s i s .
Then, from ( 6 . l ) - ( 6 . 6 ) , „(,)= a<Usin& ,i s i n 9 + ^ < u cos (*> t + d)<™6
0 0 0
where, j~ k
dU
(6,8)
'T^Hrt'^^'iy-
and, Ao and 6 are constants.
(6.7) are t h e basic equations o f t h e Il-eddy model, whose f i r s t
and second terms describe t h e primary and secondary motions, respe-
ctively. I t i s n o t i c e d t h a t f o r l a r g e r v e l o c i t y g r a d i e n t dU/dy o r
l a r g e r angular v e l o c i t y w t h e v o r t e x l i n e may be more e a s i l y r a i s e d
w i t h a longer p e r i o d o f t h e p e r t u r b a t i o n and i n an extreme case t h e
Il-eddy may be broken down by quick s t r e t c h i n g w i t h o u t any v i b r a t i o n .
However, i t i s d i f f i c u l t t o evaluate t h e a c t u a l v a l u e o f con a t p r e -
sent.
6.2.2 T u r b u l e n c e i n t e n s i t i e s and R e y n o l d s s t r e s s e s e v a l u a t e d by
a n-eddy
Assuming t h a t t h e d i r e c t c o r r e l a t i o n between t h e primary and
secondary motions i s n e g l i g i b l e , i . e . co>>w 05 the turbulence i n t e n s i -
t i e s and Reynolds stresses c o n t r i b u t e d by any n-eddy can be e v a l u a -
t e d from ( 6 . 7 ) , as f o l l o w s :
(6.10) o b v i o u s l y i n d i c a t e s t h a t t h e r e i s no c o r r e l a t i o n betveen u
and v and betveen v and v , v h i c h coincides v i t h t h e c h a r a c t e r i s t i c s
of tvo-dimensional shear f l o v . Since 0 <_ 6 <_ TT/2 and -uv" > 0,
(6.10) gives (A a> )
0 0
2
(6-11)
y^—
( 6 - 1 2 )
/ { M W l H - C W W f l ) {(aft^«+U,«jW»,
6.2.3 S p e c t r a l d e n s i t y f u n c t i o n s o f energy
The above i n v e s t i g a t i o n has been l i m i t e d o n l y t o a s p e c i f i e d eddy
element, b u t t h e c h a r a c t e r i s t i c s o f t u r b u l e n c e as a v h o l e c o n t r i b u t e d
by a l l o f t h e eddies should be made c l e a r . I n order t o a t t a i n t h i s
purpose, t h e conception o f energy spectrum about eddy s c a l e X must
be i n t r o d u c e d i n t h e same manner as Chapter 2, t o g e t h e r v i t h adopt-
ion o f space vave number i n E u l e r i a n expression i n s t e a d o f t h e above
Lagrangian form.
Although a p a r t o f t h e primary motions may produce t u r b u l e n c e
-91-
due t o t h e i n t e r a c t i o n w i t h t h e mean f l o w , i t i s considered t h a t
almost o f t h e primary motions i n t h e e q u i l i b r i u m r e g i o n have a
t u r b u l e n c e cascade process where a l a r g e r scale eddy s u c c e s s i v e l y
transports i t s t u r b u l e n t energy i n t o a smaller scale eddy, because
i n t h i s r e g i o n t h e r e e x i s t s an e q u i l i b r i u m s t a t e o f t h e t u r b u l e n t
energy and any eddy element s a t i s f i e s t h e s i m i l a r i t y law, as shown i n
Figs. 2.9 & 2.13.
Here, we adopt Heisenberg's concept o f t h e e f f e c t i v e v i s c o s i t y i n
the cascade process, as f o l l o w s . Since t h e e f f e c t i v e v i s c o s i t y v T
W~Ax • (6-15)
2 * 3
(3./X)
,2 = ••(6.17)
J
(1+. )
By t h e way, i n t h e range ( k 0 k <_ k ) where t h e cascade process
e x i s t s , t h e Reynolds s i m i l a r i t y i s v a l i d , and furthermore t h e t u r -
b u l e n t s t r u c t u r e i t s e l f has a s i m i l a r i t y i n t h e e q u i l i b r i u m r e g i o n ,
which i s t h e same conception as ' s e l f - s i m i l a r i t y 1
mentioned i n 5 . 5 .
Consequently, t h e r a t i o e may be regarded as constant because t h e
shape o f any Il-eddy i s s i m i l a r t o each o t h e r .
2
Therefore, a spectrum d e n s i t y f u n c t i o n G(k) o f to w i l l be o b t a i n e d
for k 0 < k <k do, 2
A ( 3 t / K ) *" 3
(irsintf) 'X h H
t t (6-i8)
d l 8
(i+..>*
Now, on t h e assumption t h a t u 0 i s so i n a p p r e c i a b l e compared w i t h
w i n t h e cascade process, t h e spectra o f t u r b u l e n c e i n t e n s i t i e s can
be w r i t t e n , as f o l l o w s : ,
r*(t) = -^(j.in«) G(i) £ ,(t) S I
(6-19)
-y- (7 • tal)" J
That i s .
3
/
3(1 + .
_ _£_W£>*
w ( B _ j ) X («M >.. K r X f
3 +,J)>S
(1
2__W^ ( B . ,X.X K
f r .
-93-
n
£„(*)«(«/"'to"*)*"', £ ' „ ( * ) - ( < r ~ V V 0 , («/•"'O* ( -2i)
6
vhere, *• = [ ^ J*
?i y .
e S k m
' ' «'* = \ E u ( k U k = l c t ^ i ^ K - k ^ - * ) ( 6 . 2 3 )
-97-
3
a
;i=
1
.r Ji -idy= /-ï{i-exp(-e )} !
(6-31)
J dx a 4 "t
(6-33)
y=v = o
11 dt-U*
(6-34)
'T
Then,
(6-35)'
-98-
Since u*= f(3Z/d ) . y y 0 , we can also o b t a i n from (6.30)
• 'Mf^-^jï
6.3.2 F u r t h e r c o n s i d e r a t i o n o f renewal model by t a k i n g a c c o u n t o f
the p r o b a b i l i t y d i s t r i b u t i o n o f b u r s t i n g period
The above r e s u l t s represent t h e t u r b u l e n t c h a r a c t e r i s t i c s during
-
( 6 . 3 3 ) becomes j« n
0
+ +
U^U a>(y )=U - j^GU)dsJ t(9)dr 0
i
oer (6-38)
l t 1
^J^j" Gls)dsj (eril0)) dT-U*
t t (6-39)
(^) = T
^—r-l\~ — " f ' ^ ' ^ + l w / ) ] •• (6-40)
K
U T
J
{\+/3){2aU'a
*Y 0
L
J - - io°o' Jo * jf
I I
+
+ ^(y ) - { l - e r f ( ^ r ' i ^ l ) } ] (6-41)
4
iV
where, f " g^jg expf-a'))^
v
L a s t l y , t h e b u r s t i n g p e r i o d T-g i s obtained from (6.37),
as f o l l o w s : T „ 2 . -J-s
1 giM _ 4r u +1
0
+ J
(6-42)
-99-
PART 2 Turbulent s t r u c t u r e over permeable bed
Chapter 7 Introduction
In Part 1 we have i n v e s t i g a t e d t h e o r e t i c a l l y and experimen-
t a l l y the problems ( l ) - ( 5 ) mentioned p r e v i o u s l y (see page 1 ) , which
are very basic and important i n t u r b u l e n c e research o f two-dimen-
s i o n a l s o l i d open-channel f l o w s . However, v a r i o u s t u r b u l e n t
phenomena which o f t e n appear i n h y d r a u l i c engineering and environ-
mental sciences are very complicated, and then even i n two-dimensi-
onal open-channel flows the f o l l o w i n g boundary c o n d i t i o n s are
furthermore added t o the f l o w d e a l t w i t h i n Part 1.
(1) A seepage f l o w e x i s t s under t h e porous bed. Consequently, an
i n t e r a c t i o n w i t h main and seepage flows cannot be n e g l e c t e d
i n open-channel f l o w .
(2) Suction or i n j e c t i o n e x i s t s s t e a d i l y through t h e porous bed.
Consequently,the t u r b u l e n t s t r u c t u r e may be v a r i e d e s s e n t i a l l y .
( 3 ) The bed surface i s d i s c o n t i n o u s , uneven or wavy.
(4) The bed i s moving sand surface. Consequently, an i n t e r a c t i o n
between f l o w and moving bed should be considered ( Sediment
t r a n s p o r t by t u r b u l e n t f l o w ) .
(5) Suspended m a t e r i a l s e x i s t i n open-channel f l o w . Consequently,
the t u r b u l e n t s t r u c t u r e may be v a r i e d e s s e n t i a l l y .
In t h i s p a r t , we i n v e s t i g a t e the t u r b u l e n t s t r u c t u r e s w i t h
boundary c o n d i t i o n s o f (1) and (2) i n Chapter 8 and Chapter 9,
r e s p e c t i v e l y , and we consider the fundamental c h a r a c t e r i s t i c s
o f turbulence over a permeable bed, by comparing w i t h t h e r e s u l t s
obtained over a s o l i d bed (Part 1 ) .
-109-
Chapter 8 I n t e r a c t i o n between main a n d s e e p a g e f l o w s * )
8.1 Introduction
When t h e bed p e r m e a b i l i t y i s small enough t o n e g l e c t t h e
e f f e c t o f seepage f l o w under the bed, t h e main f l o w shows t h e chara-
c t e r i s t i c s over a s o l i d rough bed, which have been already discussed
i n Part 1. When the e f f e c t o f seepage f l o w on the main f l o w
cannot be n e g l e c t e d , t h e r e e x i s t s a hydrodynamic i n t e r a c t i o n b e t -
ween them through t h e porous bed. Compared w i t h s o l i d rough bed,
the momentum exchange becomes more a c t i v e , and consequently an
a d d i t i o n a l shear s t r e s s w i l l be produced by t h i s i n t e r a c t i o n .
For example, Lovera e t al.(1969) r e p o r t e d t h a t t h e f r i c t i o n - f a c t o r
o f f l a t permeable bed i n a n a t u r a l r i v e r was much l a r g e r than t h a t
o f rough bed given by Nikuradse. A l s o , Monin e t al.(1971)
p o i n t e d out t h a t s i m i l a r phenomenon was observed i n the wind f l o w
over a vegetated ground. These a d d i t i o n a l s t r e s s may be caused
by an i n t e r a c t i o n between main and seepage f l o w s .
The seepage f l o w near t h e pore surface i s e x c i t e d t o be
t u r b u l e n t by t h e pressure f l u c t u a t i o n s o f main f l o w . Otherwise,
the seepage f l o w becomes laminar and Darcy law i s a p p l i c a b l e .
Then, there i s an analogy between smooth and permeable f l o w f i e l d s ,
as shown i n Fig.8.1. That i s , t h e viscous sublayer i s also
d i s t u r b e d by the pressure f l u c t u a t i o n s o f main f l o w , as have been
mentioned i n 4 . 7 . 3 . Consequently, an a d d i t i o n a l shear s t r e s s i s
probably induced i n the main f l o w by t h e d i s t u r b e d seepage f l o w .
To sum up, t h e r e e x i s t s a hydrodynamic i n t e r a c t i o n mechanism
between main and seepage f l o w s , as shown i n F i g . 8.2. This
chapter deals w i t h these feed-back systems, i n order.
-110-
8.2 Theoretical consideration
8.2.1 Basic equations o f seepage flow
The l o c a l v e l o c i t y Vp i n the porous media must obey the f o l l o -
wing Navier-Stokes equation.
The i n e r t i a l f r i c t i o n i s a n o n l i n e a r term, as f o l l o w s :
((»,'D» ) 5 - | , | , -(8-4)
vA
where, K i s the i n t r i n s i c p e r m e a b i l i t y o f the medium, and C i s a
experimental constant.
D e f i n i n g the Reynolds number « = T l l / x o f seepage f l o w as the A v 9
F r x CR
kl\-- - ^-jt + ^ (8-5)
-111-
( 8 6 )
y-component : f = "9 i* ~ y > '
x-component : (1 +cit )q K H =(K /»
9 ) /, (8-7)
*JL- _ = *(£-)-„„,
v
dt dx fi
(8-11)
8
lr=- l(l)-a„,
v
dt dy fi '
„„ . __£ï-?ï*
v (8-12)
(8-13)
a y ox ,
-112-
8.2.2 A n a l y s i s o f t u r b u l e n t f l u c t u a t i o n s o f seepage flow
We consider the t u r b u l e n t f l u c t u a t i o n s o f seepage f l o w e x c i t e d
by the pressure f l u c t u a t i o n p 0 on t h e boundary (y=0) o f main f l o w
(see Fig. 8.1 ( b ) ) . That i s , the boundary c o n d i t i o n s o f seepage
flow are assumed as f o l l o w s :
h
, at y— p (8-i5)
!
where, Re{ } denotes t h e r e a l p a r t o f { }. ?o'(=v^ö ) i s t h e r.m.s.
value o f pressure f l u c t u a t i o n s i n t h e main f l o w . h i s the
P
thickness o f porous media, a t which t h e t u r b u l e n t f l o w i s attenuated
and becomes laminar, as mentioned p r e v i o u s l y . Then, on these
boundary c o n d i t i o n s , (8.11)-(8.12) are e a s i l y solved, as f o l l o w s :
C-h < y < 0)
p
P {
' «h(^/ )J (8-16)
9 { X m i ( i j r + < u < ) } ( 8 , 1 7 )
'^' (7^717) c^vv
. nk* /^(po'/^sinh^U.+y) .
?
>-*'t-(^V' coshC^) * exp.(^)} (8-18)
where,
1
By using t h e r e l a t i o n o f ,» = 1 [?,?,•+«,{ J, }] =-J- ö^V
(q x and q* are complex conjugates), t h e f o l l o w i n g approximations
are obtained i n t h e r e g i o n very close t o the porous boundary, i . e .
0 ~ |y| « h . — n'k*
P d = f - T T T \ ( p ^ y (8-19)
1
n k*
l
ll = 7 = = = f = = T = = = . ( ? ^ ) - ( 8 - 2 0 )
-113-
( 8 2
That i s , ?0 * . « ) = ^ V * . « ) and
(4 J'J^V*' » ' °
y :
" 9 , V ( a , W ) ( a ^ + of )
E (M+O
0 = -\ /»U)COS <(/i+t )dt (8-27)
-114-
Now, assuming t h a t t h e a u t o - c o r r e l a t i o n has an exponential
t y p e , we o b t a i n j r - 1tI/i)< -,„ d r = ? ( 8 . 2 9 )
T h e n > { 8 < 3 1 )
' ' ^ ^ T i . ^ . ) '
However, t h e c a l c u l a t e d values o f q« and q' u s i n g (8.22), (8.23) and
(8.31) are d i v e r g e n t , because t h e higher wave number components o f
(8.31) are t o o l a r g e . (Same example appears i n (2.43) by using
(4.47).) Hence, i t may be b e t t e r t o use t h e f o l l o w i n g f i t t i n g
f u n c t i o n which s a t i s f i e s t h e c o n d i t i o n o f A'(0)=0, since (8.28) o f
the f r o z e n turbulence i s probably v a l i d when t i s very s m a l l .
)= i s a n
^' T+7J777 (^2 experimental constant) -(8-32)
Then, 1 , 1 v
£ o U + D = — e x p ( - — l / i + l l ) - (8'33)
( i ,„ ) = + -(8.34)
I ( , , ) ] (8 36)
' - - ^ ^ ) ^ T 7 i 7 V r T T ^ ^ ' - ^ '
0 J
where, =«/VC«,"-D> • * =i-W°.>°
-115-
If (8.30) i s used i n s t e a d o f (8.32), 2b 2 i n t h e denominator o f (8.45)
i s replaced o n l y by TTbj .
Where, _ <n
w n i c n
~Jo V(.a "l + t ) { a 7 " + t ) U + t ) ' i s easily integrated
as follows:
f
( i ) «/S«! >l ( When \ > a ( > a { , t h e lower s i g n i s used)
/ = J = j , a n - ( ; a
- L ^ ) A
7 v
Ian (8'46)
where, t h e c o e f f i c i e n t a i s a f u n c t i o n o f t h e w a l l roughness k*
and t h e Reynolds number Re, i. e. ; a= G(Re, kg ) from (4.34).
-118-
The i n t e g r a l time s c a l e Z can be approximately given by
ï~o.4h/u 0 , since v Z e i s n e a r l y equal to the m a c r o s c a l e L x i n the
main flow * J .
From the above c o n s i d e r a t i o n s , we can o b t a i n the f o l l o w i n g
formula which e s t i m a t e s the induced s t r e s s T j .
r
- ^ ^ ^ . - ( (8-50)
r 0 »(2A,) *Co Wo
where, .
F i g . 8.5 g i v e s t h a t 2>i=0.12 or b =0.16, and thus Trbi=0.38 or
2
-119-
8.3 Experimental consideration
8.3.1 Experimental set-up and procedure
As shown i n F i g . 8.11, t h e glass beads o f diameter d =1.25 cm m
i n F i g . 8.12.
Now, d e f i n i n g the f r i c t i o n f a c t o r f K o f seepage flow as
,
follows: /. /A: (2/VT)-(^/2 ,)
S
j ( Darcy-Weisbach type, c f . ( 2 . 7 ) ) ,
is applicable.
W e l l , the w a l l s h e a r s t r e s s x 0 or the f r i c t i o n v e l o c i t y U #
should be e v a l u a t e d as a c c u r a t e l y as p o s s i b l e , s i n c e t h e induced
s t r e s s i s probably not so l a r g e . Although t h e e v a l u a t i o n
method of U # has been a l r e a d y d i s c u s s e d i n 4.3.1, (2) log-law
method should not be adopted h e r e . I t ' s because t h e Karman
constant K i s not n e c e s s a r i l y equal to K=0.4. A c t u a l l y , Chu
et a l . ( 1 9 7 3 ) have o b t a i n e d K=0.27 i n a permeable a i r p i p e flow.
T h e r e f o r e , we have e v a l u a t e d the v a l u e s o f from (3)
R e y n o l d s - s t r e s s method. F i g . 8.14 shows an example o f the
d i s t r i b u t i o n o f Reynolds s t r e s s which was measured over the
permeable porous bed by the d u a l - s e n s o r h o t - f i l m anemometer.
Although t h e r e i s some s c a t t e r i n g o f t h e measured data n e a r t h e
bed because of t h e l o c a l e f f e c t o f t h e roughness elements, these
2
values o f -üv/U . agree f a i r l y w e l l w i t h (4.20). Consequently,
we can conclude t h a t t h e Reynolds-stress method adopted here i s
reasonably v a l i d .
Next, F i g . 8.15 shows t h e v e l o c i t y d e f e c t law o f the main
flow. A good agreement between the observed values and (8.47)
i s s u r e l y recognized. Consequently, the values o f U 0 and (8.48)
are reasonably e v a l u a t e d from (8.47).
I n Table 8.2 which summarizes the above r e s u l t s , i t should be
n o t i c e d t h a t the Karman constant K tends t o decrease as the porous
bed becomes l o o s e r . W e l l , i t i s well-known t h a t the Karman
constant can be v a r i e d by t h e o u t e r c o n d i t i o n s . For example,
these phenomena have been observed i n f l o w s w i t h suspended m a t e r i -
als, w i t h s u c t i o n o r i n j e c t i o n (see next c h a p t e r ) or w i t h loose
boundaries. Ippen(1971) p o i n t e d out t h a t t h e Karman constant
i s v a r i e d by the boundary c o n d i t i o n s near t h e bed. That i s t o
say, the Karman constant may decrease as t h e momentum exchange o f
t u r b u l e n c e becomes more a c t i v e near t h e bed, and v i c e versa.
By t h e way, P h i l l i p s ( 1 9 6 7 ) h a s e v a l u a t e d A =0.12 w i t h d e v i -
m
References ( Chapter 8 )
6) Arbhabhirama, A. & Dinoy, A. A. : Friction factors
1) Lovera, F. & Kennedy, J. F. : Friction-factors for
and Reynolds number in porous media flow,
flat-bed flows in sand channels, Proc. of ASCE, Proc. of ASCLVHY-6,pp.901-911, 1973.
pp. 1227-1234, 1969. 7) Corcos, G. M. : The structure of turbulent pres-
v u
2) Monin, A. S. & Yaglom, A. M. : Statistical fluid " ' 1.' . . . , n r l l l i H M p r h
w,T n ., ->C7 i n 10-71 sure field m boundary layer flows, J. Fluid Mech.,
mechanics, MIT Press, vol.1, pp.257-327, 1971.
•>v w . . «, #\ i_ r r vol. 18, pp.353-378, 1964.
3) Mi es, J. W. : On the generation of surface waves n x " .,, „, , ,„ , . . , « . . . ,,„„,„.
' ' „ , ,. ,, , o t 8) Willmarth, W. W. & Wooldridge, C. t. : Measure-
by shear flows, J. Fluid Mech., vol.3, pp.185-
r r
~ , , „,„ .„„ M
' ments of the fluctuating pressure at the wall
204, 1957.
beneath a thick turbulent boundary layer, J.
4) Phillips, O. M. : The maintenance of Reynolds
Fluid Mech., vol. 14, pp. 187-210, 1962.
stress in turbulent shear flow, J. Fluid Mech.,
9) Lighthill, M. J. : Physical interpretation of the
vol.27, pp. 131-144, 1967.
mathematics theory of wave generation by wind,
5) Chu, Y. & Gelhar, L. W. : Turbulent pipe now
J. Fluid Mech., vol.14, pp.385-398, 1962.
with granular permeable boundaries, Ralph M.
nParsons Labo.,
i u MIT, n . vNo.
U I T Report i 148,
u o 1972.
ion ' ° ) Ippen, A. T. : A new look at sedimentation in
r r
-124-
transpiration. However, the r e c e n t r e s e a r c h o r new intention
is probably concerned w i t h the t u r b u l e n c e c h a r a c t e r i s t i c s w i t h
t r a n s p i r a t i o n and the i n f l u e n c e of t r a n s p i r a t i o n on the w a l l - t u r b u -
lence mechanism or b u r s t i n g phenomenon. For example, r e f e r to
the papers by Andersen e t a l . ( 1 9 7 5 ) ( b o u n d a r y l a y e r ) , Schildknecht
et a l . ( 1 9 7 5 ) ( p i p e flow) and p r e s e n t a u t h o r s ( 1 9 7 5 ) ( o p e n - c h a n n e l ) .
In order to examine the e f f e c t of s u c t i o n o r i n j e c t i o n on the
t u r b u l e n t s t r u c t u r e i n open-channel flow, t h e p r e s e n t c h a p t e r
f i r s t l y e s t a b l i s h e s the e v a l u a t i o n method o f the f r i c t i o n velocity
and makes c l e a r s y s t e m a t i c a l l y the mean v e l o c i t y distribution,
i n c l u d i n g even the c a s e o f s t r o n g s u c t i o n or i n j e c t i o n a t which
r e l a m i n a r i z a t i o n or s e p a r a t i o n phenomenon may occur, respectively.
Next, the t u r b u l e n c e c h a r a c t e r i s t i c s such as t u r b u l e n c e intensities
, s p e c t r a l d i s t r i b u t i o n , t u r b u l e n t energy budget and so on a r e
investigated.
"-*-£(•£)* <»•»
Since (2.4) i s s t i l l v a l i d because of ( 2 ) , the f o l l o w i n g e q u a t i o n
for the shear s t r e s s x ( y ) i s obtained from (2.1) by u s i n g (2.4)
and (9.1).
-125-
j = -j--gI.y+v,U+r(y) (9.2)
US = TJ =gl,h—!L-^u>dy
P (9.5)
and t = y/h .
The second and t h i r d terms o f (9.4) are t h e a d d i t i o n a l ones caused
by t r a n s p i r a t i o n v . 0 Of course, when v = 0 ( s o l i d bed), these o
Ummx aX I (9.7)
e,ao=j%'d(
(9.3) becomes v( )= { Q . e . c o - O . ^ o e i C O J f / m . x ( 9 . 8 ) .
y
2
becomes U+ =gh , sin 9
t t t since t h e normal f l o w depth h 0 i s formed
(where, U^o i s t h e f r i c t i o n v e l o c i t y when v = 0 . ) . o
-126-
Substituting (9.9) and (9.10) f o r ( 9 . 8 ) , we can obtain
+ +
r + = _ L _ = ( l - 0 + V ( ^ - ^ m . x ) +*(« (9.11).
+
The c o r r e c t i o n term • t f ) « r * « W - r a ) i s c a l c u l a t e d as f o l l o w s : *)
C 9 1 2 )
;-e,m -
flss
where, ([/*/t7*.)*=r/r( .=0)| .
w (= and the ( + ) - s u f f i x denote the
dimensionless q u a n t i t y by U * and v, as mentioned p r e v i o u s l y .
Consequently, when t h e mean v e l o c i t y d i s t r i b u t i o n nCO i s known,
the c o r r e c t i o n term * ( Q can be e a s i l y e v a l u a t e d . For example,
when t h e 1/7-power v e l o c i t y law i s a p p l i c a b l e , (9.12) becomes
+
for Fr« 1, as f o l l o w s : (7/9)|t>. £/ . +f(f"-i) 1^0.072v,*U +
m x mmx y
which i s f a i r l y small * J .
(9.15)
Then, r/ = P (l - f ) + v ( u- f l/oux )
0
increases.
0 I f
shear
w e notice
s t r e s s a t the f r e e s u r t a c e \y n j
f„ÏÏ^A
moderate r a t e o f t r a n s p i r a t i o n , i . e . (9.13) ,
</tl«-l"
t h e n we can obtain
( f ) e c o ( 9 1 4 )
i_e,ci)) -iv{ ( - i V W h - ^ ' - ' '
+
y (f)= (
ro i 9 i a n d (9 14*1 i s n e g l i g i b l e i n the
The d i f f e r e n c e between (9.12) and ( . y . i ^ j
s u b c r i t i c a l flow ( F r « l ) •
-127-
fi becomes n e a r l y equal t o zero a t 3- 0.3, where a s e p a r a t i o n o f t h e
f l o w may appear. On t h e o t h e r hand, fi increases f o r s u c t i o n , and
thus t h e v e l o c i t y g r a d i e n t dU/ dy\, .=US/v = a i s o increases. The
curves o f vt=&/vfï vs.$ show t h e c h a r a c t e r i s t i c s t h a t V o i n c r e a s e s
r a p i d l y f o r i n j e c t i o n , w h i l e i t approaches about -0.1 f o r suction,
where a r e l a m i n a r i z a t i o n may occur.
F i g . 9.5 shows t h e r e l a t i o n s h i p between f = f f i and Re, as a 0
We c a l l i t t h e 3 - f i method.
9.2.6 Some t u r b u l e n c e c h a r a c t e r i s t i c s
(1) Reynolds s t r e s s distribution
When t h e c o r r e c t i o n term 0(£) i n (9.11) i s n e g l i g i b l e , t h e
Reynolds s t r e s s i s given by
— ^ +
= ( i - f ) + ^ ( i / - f i / ™ x )-~r (9.39).
-134-
Hence, as a parameter o f $=v /U* o 0 f o r any Re, fi i s c a l c u l a t e d from
( 9 . 3 7 ) , and then -ïïv/U^o i s c a l c u l a t e d from (9.39) w i t h R^=/f72 Re
and Vo = 3/>^- An example o f t h e s e r e s u l t s i s p l o t t e d f o r Re=5xl0 l1
In c o n t r a s t w i t h the r a t i o fi of t h e w a l l shear s t r e s s , y i n c r e a s e s
for i n j e c t i o n , while y decreases f o r suction. At B = -0.1 (v - 0
( 2 ) T u r b u l e n t e n e r g y budget
The t u r b u l e n t energy p r o d u c t i o n P i s given by
t 9 4 1 )
i-*(=§)(£) - -
Fig. 9.9 shows the a b s o l u t e v a r i a t i o n s of P which were c a l c u l a t e d
from (9.41) by u s i n g (9.18) and (9.39). Except f o r near the w a l l ,
the t u r b u l e n t p r o d u c t i o n P a l s o i n c r e a s e s w i t h an enlargement of
i n j e c t i o n , w h i l e i t d e c r e a s e s with an enlargement of s u c t i o n .
-135-
On t h e o t h e r hand, the opposite v a r i a t i o n c h a r a c t e r i s t i c s appear
near t h e w a l l . I t ' s because t h e Reynolds s t r e s s increases f o r
i n j e c t i o n , w h i l e t h e viscous s t r e s s increases f o r s u c t i o n , as
compared w i t h n o - t r a n s p i r a t i o n ( v = 0 ) . o
By newly d e f i n i n g t h e t u r b u l e n t d i f f u s i o n T r as T r H ( „ + «.)O72 0 ,
(9.43) coincides w i t h (2.30). Although (9.42) o r (9.43) i s n o t
d i r e c t l y examined i n t h i s chapter, i t i s suggested t h a t T /U* r 0
+ !
E =
0 £«-/fJ = E* • fl
0 (9.44),
!/o = y f t > yVvfi"
-136-
(9.45) corresponds t o (2.8) or ( 2 . 9 ) . Defining % = & / £ „ (where
to i s the value a t v = 0 ) ,o (9.45) becomes
< = ^ / ¥ ( (9.46)
9.3 E x p e r i m e n t a l c o n s i d e r a t i o n on mean v e l o c i t y d i s t r i b u t i o n
9.3.1 E x p e r i m e n t a l s e t - u p and p r o c e d u r e
Fig. 9.11 ( a ) shows t h e e x p e r i m e n t a l s e t - u p i n which t h e
r e c i r c u l a t i n g t i l t i n g flume i s 15 m long, 50 cm wide and 30 cm deep.
About 9 m downstream from t h e channel e n t r a n c e , an apparatus o f
t r a n s p i r a t i o n flow was s e t up ( s e e F i g . 9.11 ( b ) ) , whose porous
p l a t e made o f s i n t e r e d p l a s t i c s o f 0.65 mm diameter was 104 cm long,
49 cm wide and 1.5 cm t h i c k . I n order t o u n i f o r m l y g i v e any
t r a n s p i r a t i o n v e l o c i t y v , a b u f f e r zone made o f g l a s s beads and
0
-137-
The mean v e l o c i t y d i s t r i b u t i o n s were measured i n d e t a i l by a
s i n g l e - s e n s o r h o t - f i l m anemometer (DISA 55A83, s e e F i g . 3 . 3 ) .
Four kinds o f experiments were performed (Case H-1, H-2, H-3
and H-4), a s d e s c r i b e d i n T a b l e 9.1. Case H-1 and H-4 were t h e
most important experiments, which correspond t o Case A-1 and G-1
in Table 4.1, r e s p e c t i v e l y . I n the case o f H-1, H-2 and H-3,
two i n j e c t i o n s and two s u c t i o n s were g i v e n . I n t h e c a s e o f H-4,
t h r e e i n j e c t i o n s and t h r e e s u c t i o n s were g i v e n . The maximum
quantity Q t of these t r a n s p i r a t i o n s was l e s s than t h e 10% o f t h e
main flOW quantity Q. xkfctoM Hydraulic data for experiments with suction or injection
r a t e s were g i v e n , P •0.243 -0.202 -0.137 -0.094 -0.050 0.0 0.049 0.089 0.148 0.198 0.249
v: •0.111 -0.102 -0.084 -0.066 -0.041 0.0 0.062 0.136 0.319 0.589 1.036
as d e s c r i b e d i n •12.6Ï -10.5 -7.1 -4.9 -2.6 0.0 2.5 4.6 7.7 10.2 12.9
Q /Q
t
Case H-2 Case H-3
Table S I . This 9 8 1 4 5 8 7 1 4 5
Run
2.066 1.684 1.369 1.062 0.799 2.307 2.166 1.907 1.509 1.147
chapter i s s u p p l e - B
-0.111 -0.052 0.0 0.056 0.111 -0.047 -0.031 0.0 0.051 0.104
+
-0.074 -0.042 0.0 0.072 0.189 -0.039 -0.027 0.0 0.065 0.172
mented by t h e s e q ° / Q
-9.9" -4.6 0.0 5.0 9.9 -5.BZ -3.8 0.0 6.3 12.8
r e s u l t s of Table S I ,
i f necessary, i n the following sections,
9.3.2 Wall s h e a r s t r e s s o r f r i c t i o n v e l o c i t y
1.0
v e l o c i t y , U/U
3
vs. y/h. For i n j e -
max
c t i o n , t h e lower v e l o c i t y r e g i o n i s
l i f t e d up from t h e w a l l , w h i l e f o r » • ...
suction the higher v e l o c i t y region i s
attracted toward t h e w a l l . Surely,
1.0
«/>Wx
the e f f e c t o f t r a n s p i r a t i o n on t h e Fig. S3 Mean velocity distribution with
suction or injection
-138-
t u r b u l e n t s t r u c t u r e i s recognized remarkably near t h e w a l l , even
when 3= v /U\ 0 o i s below ±0.1.
Well, i n our previous research( t h i s d o c t o r a l d i s s e r t a t i o n ,
1977), the f r i c t i o n v e l o c i t y U + could be evaluated by t h e extended
methods o f (2) and (3) i n 4.3.1. I n t h e (2) log-law method,
the friction velocity was determined by (9.22), where t h e Karman
constant < was given by (9.24). On t h e o t h e r hand, i n t h e (3)
Reynolds-stress method, t h e f r i c t i o n v e l o c i t y was determined by
comparing (9.39) w i t h t h e measured Reynolds s t r e s s d i s t r i b u t i o n s .
Fig. 9.13 shows t h e r e s u l t s evaluated by both these methods.
In t h i s experiments, they agree w i t h each o t h e r w i t h i n t h e e r r o r o f
1
— 1 5 % . F i g . 9.14 shows t h e v a r i a t i o n o f U* vs. 3,
together w i t h t h e t h e o r e t i c a l curves o f (9.37). Also, a good
agreement between t h e experimental and t h e o r e t i c a l values i s
recognized i n F i g . 9.14.
In t h e c o n s i d e r a t i o n o f t h e above, we could conclude i n Paper
( I ) t h a t the 3-ft method described i n 9.2.5 i s t h e most reasonable
and the e a s i e s t f o r e v a l u a t i o n o f U^, because (9.37) i s r e a l l y v a l i d .
The values o f U A described i n Table S I were determined by t h i s
3-ft method.
9.3.3 Mean v e l o c i t y d i s t r i b u t i o n
Fig. 9.15 shows t h e dimensionless mean v e l o c i t y d i s t r i b u t i o n
+
U =U/U i n the case o f H-1 and H-4 ( o l d d a t a ) .
+ These experimen-
tal values agree f a i r l y w e l l w i t h t h e t h e o r e t i c a l curves o f (9.18)
or the b i l o g - l a w o f (9.21).
Now, F i g . S5 shows t h e dimensionless mean v e l o c i t y d i s t r i b u -
t i o n w i t h i n j e c t i o n i n t h e case o f H-1 (new d a t a ) . Run 1 i s a
case w i t h o u t t r a n s p i r a t i o n , i . e . v = 0 , whose data c o i n c i d e very
o
-139-
throughout t h e whole r e g i o n from the bed up t o t h e f r e e s u r f a c e .
That i s t o say, i t i s confirmed t h a t t h e m i x i n g - l e n g t h model o f
(9.17) i s v a l i d even i n t h e case o f an open-channel f l o w w i t h
injection. Above a l l , even i n t h e case o f t h e l a r g e i n j e c t i o n
r a t e v£ =1 where s e p a r a t i o n c h a r a c t e r i s t i c s may appear, a good
agreement i s recognized. A l s o , an examination o f v a l i d i t y o f
(9.32) i s shown i n F i g . 9.16. S u r e l y , (9.32) agrees f a i r l y w e l l
w i t h the observed values. Then, i t i s o b t a i n e d t h a t K =0.33,
0
+ + + + +
U „ t / l = C, In U „ y l + C , U 0 l • (-0.K ^<~0.04 ) .. (9 . 47)
-141-
He obtained Ci=0.06 and C = l l i n a t r a n s p i r e d boundary l a y e r .
2
-142-
from (9.10) by using evaluated by t h e 3-fi method. These c a l -
c u l a t e d values are shown by s o l i d l i n e i n F i g . S9.
Although there i s some s c a t t e r i n g i n t h e observed data, they
agree f a i r l y w e l l w i t h t h e c a l c u l a t e d curves o f (9.10). Consequ-
e n t l y , we can s u r e l y conclude t h e v a l i d i t y o f t h e 8-fi method and
the b i l o g - l a w .
9.3.5 R e l a m i n a r i z a t i o n by t h e s t r o n g s u c t i o n
From t h e above c o n s i d e r a t i o n , i t has been c l a r i f i e d t h a t t h e
mean v e l o c i t y d i s t r i b u t i o n i n an open channel w i t h t r a n s p i r a t i o n
except f o r t h e s t r o n g s u c t i o n s a t i s f i e s t h e b i l o g - l a w o r t h e mixing
l e n g t h model f a i r l y w e l l throughout t h e f l o w depth. On t h e contra-
ry, as t h e s u c t i o n r a t e becomes s t r o n g e r , t h e main f l o w cannot keep
the f u l l y developed t u r b u l e n t s t a t e any longer and consequently
a r e l a m i n a r i z a t i o n may occur, as mentioned p r e v i o u s l y . I t i s well
known t h a t such a r e l a m i n a r i z a t i o n also occurs i n a boundary l a y e r
w i t h a s t r o n g l y favourable pressure g r a d i e n t (dp/dx <0, F i g . 9.3).
Hence, i t i s d e s i r a b l e t o e x p l a i n these r e l a m i n a r i z a t i o n phenomena
systematically.
By t h e way, Huffman e t al.(1972) p o i n t e d o u t t h a t t h e o u t e r
boundary c o n d i t i o n s such as t r a n s p i r a t i o n and pressure g r a d i e n t ,
a f f e c t e d t h e i n n e r l a y e r p r i m a r i l y v i a t h e shear s t r e s s gradient
3 /9y i n t h e inner l a y e r .
T I n t h e same manner as they d i d , t h e
+
average value o f Sr/Sy over t h e w a l l r e g i o n o f 0 < y < 6 =100
i s now d e f i n e d by Ox/3y>. Then, t h e f o l l o w i n g r e l a t i o n i s
obtained from (9.11). nt-MbUV 1
[ c m H-1)
" KtlulMrlutton
3+
'\»W
<Z>(3) +
» : (1»»il4, )/R.
(9.49) 1 : Fr-0.17
x
2 : Ft-0.60
3 : Fr'O.eo,
« : Fr-1.20
Fig. S10 shows t h e curves o f (9.49) 5
O
: Fr-1.50
: Q /Q-0.1
t
f o r v a r i o u s Froude numbers.
Although the d e t a i l e d d i s c u s s i o n i s Fijr.SIO Average gradient of the shear stress
S10
o m i t t e d here (see P a p e r ( I ) ) , t h e
-143-
c r i t i c a l value vt =-0.085 f o r r e l a m i n a r i z a t i o n corresponds to
+ +
<3i /3y > =-0.009 i n the case o f small Froude number.
On the o t h e r hand, Patel e t al.(1968) also discovered t h a t the
r e l a m i n a r i z a t i o n o f a c c e l e r a t e d t u r b u l e n t f l o w occurred a t a c r i t i c a l
value o f shear s t r e s s g r a d i e n t , independently o f o v e r a l l Reynolds
number. That i s , when the shear s t r e s s could be described by
r = W y ••• (9-50), Two-d1»enslonal Flows
tuil Jet
C l o s e d Channel
A P a t e l 1 Head ( 1 9 6 9 )
by a parameter o f d^va/pUS =-0.009 A Laufer (1950)
Boundary L a y e r
• J u l l e n et *1.(1970)
+ +
This value i s equal t o < 3 x / 8 y > c
• Narayanan et al.(1969)
Open C h a n n e l
— • Authors
-0.009, which i s very remarkable. 1 : Re-SCOO
2 : Re'5000
3 : Re-11600(Case H-1)
A*-27
+
t i o n o f the damping f a c t o r A again- Fr-0.2
Fr-0.8
-144-
w h i l e i t i n c r e a s e s r a p i d l y w i t h an enlargement ofsuction.
The e x p e r i m e n t a l curve o f ( T ) f o r an open channel i s n e a r l y g i v e n by
+
f=27/(l+4.1 V ) Cf. < 0 ) | (9.51).
+
1=27/(1+2.71-.+) ( v , > 0 ) J
9.4 E x p e r i m e n t a l c o n s i d e r a t i o n on t u r b u l e n c e c h a r a c t e r i s t i c s
-145-
Acknowledgements
Academic c a r e e r
Kyoto, JAPAN
-156-
CONTENTS
( FIGURES AND T A B L E S )
Page
PART 1 TURBULENT STRUCTURE OVER SOLID BED
1
1 2
^ " c l r c o n s . a e r a t l o n on t u r b u i e n t s t r u c t u r e
Chapter 2 8
Method o f t u r b u l e n c e measurement and i t s a n a l y s e
Chapter 3
F x D e r i m e n t a l c o n s i d e r a t i o n on 13
Chapter 4 t u r b u l e n t s t r u c t u r e and e n e r g y b u d g e t
F x n e r i m e n t a l c o n s i d e r a t i o n on 30
Chapter 5 S a n i s m of t u r b u l e n c e - p r o d u c o n t l
43
P h y s i c a l model o f w a l l turbulence
Chapter 6
rhavter 7 Introduction 50
cZL s i n t e r a c t i o n between main and s e e p a g e f l o w s 54
cwter 9 Turbulent s t r u c t u r e with transpiratmn
8m£
Ch.1
Applled Physic fluid Experimental
M e t h o t
& Hl thematlct Dynamics Hydraulics
Ch.2 Theoretical Approach Ch 3 ' °' Turbulence
Measurements
Boundary L a y e r Theory O
tr DJ
GOttlnqen-school
DJ h
•school
rt.
„
t n
i
* »
Measurements o f
Basic Quantities
Turbulence
of
Ch ^ S t r u c t u r e of
Reynolds-stress
Fluctuating
art
H
* Turbulence I n t e n s i t y 0 Probability Density Function Mot-wire Anemoaetcr Amerlkan-school rn
9 Reynolds S t r e s s 0 C o n t r i b u t i o n s of Etch Event z:
% A c t i v e and I n a c t i v e M o t i o n s • F r i c t i o n o f Time by E t c h E v e n t
Theory of Isotropic H
Turbulence
cr
r . / Turbulent E n t i t y
\ C h g Hech«n1sa of
Cawbrldoe-school _
3
r+
~S (TJ
oo
H
X)
C h
- *• Budget
Tu rbuIence-produc11on O n
Locally
Turbulence
Isotropic
Q- R+ n
P Spectral Function
•> T u r b u l e n t P r o d u c t i o n ,
0 B e h a v i o r of Con*rent Notions
• P e r i o d o f B u r s t i n g Phenomenon O H
D i s s i p a t i o n and D i f f u s i o n 1 T M O - f a c e d n « s s of T u r b u l e n t Run lan-school
r+ r+ C
Structure —I. ~J
Experiments of o C
Free-turbulence 3 O
Experiments of r+
C Flow V l s u a l l i a t l o n of
Ch. 6 Hodel of T u r b u l e n t Flow Turbulence Wall-turbulence -j
Boundary-layer
t> n-eddy Model % Instantaneous Velocity P r o f i l e t Pipe flow o
? Renewal Model ft l o w - and H i g h - s p e e d S t r e a k Transpiration
< if
* Combined Model
HACA-school
Flov
fD
1 z
Hot-wire. Pressure-transducer etc. T u r b u l e n t Open- i/i
I
c h a n n e l Flow
Nean-
o r
veloctty >
Ch. 7 Introduction Profile
z
cr m
Q_
Interaction between Turbulent Structure -n
ch.8 Main and S e e p a g e Flows
cK.g with Trintpfrattoft
r
o
ft Model of Seepage Flow « Effect of Suction Flow
p> E v a l u a t i o n o f Induced Stress « E f f e c t of Injection Flov
M.
SERVO AMPLIFW
G
Outjjxrt to incfiesbng n
DISA 55A83
8
Fig. 3.9 Output voltages of anemometers in s t i l l water.
(b)
o
F i g 3 J 0 C a l i b r a
Fig. 3.11 Temperature c h a r a c t e r i s t i c s - tion curves of dual-sensor.
of c a l i b r a t i o n c o e f f i c i e n t A (TW).
6
1- in-, f-3 m
CO
Sensibility Factor ol 1-st Sensor
ro
O O
-h 01
I C
t/> —«
r-l- cv
r+
l/i rt>
(Tj O .
l/l O
O tJJ
-! —'
60 80 cr
-S
U icrry )
s
CD
O
ro
on
DISA 55005
Hot-Film Anemometer Mufti-Channel
Analog 5k.
Data 5 F.H. Data Recorder
DISA 55D05
H o t - F i l m Anemometer / •
Low-Pass Filter
Digital Data
\ <
FAC0M
A-0
U-200
Converter
1 '
FAC0H 230-75
Large Computer
60
CI
U( TM>i
F i g . 3.13 Data processing system
of point-measurements.
Fig. 3.15 Temperature s e n s i t i v i t y of hot-film sensors.
Fig. 3.16 Flow v i s u a l i z a t i o n system (orientation of H -bubble wires
2
(b)
Pulse Generator Reference
Ajtis
Sequence of
Kot i o n Picture
IS™ Film
Digital Projector
Coordinator
y(U-a
Data Card FACOM 230-75
Urge Computer
12
chapter 4 E x p e r i m e n t a l c o n s i d e r a t i o n on
t u r b u l e n t s t r u c t u r e and e n e r g y budget
Head tank
\f
15 m -
•9.5 I B -
ValvSI Test section
[Screen
1
Ground pool Rectangular
Pump weir
( a ) H y d r a u l i c c o n d i t i o n s t h a t t h e Froude number i s n e a r l y c o n s t a n t ( S m o o t h b e d ) .
Slope
h U Urn Re Fr Q Tw
Case max
(cm) (cm/s) (cm/s) (cm/s) (cm/s) (cm/s) xlO 3
xlO 4
xlO" 4
(1/s) PC)
A-3 3.14 46.1 40.7 2.206 2.152 2.132 0.643 1.18 0.733 14.8 6.4 17.18
fl-S 5 . 9 9 60.0 53.8 2.485 2.485 2.507 1.33 2.89 0.703 10.7 16.1 15.75
A-6 8.00 67.0 60.1 2.638 2.638 2.732 2.03 4.62 0.689 9.5 24.1 18.60
G-1 8.01 48.3 40.1 2.152 2.468 1.896 1.63 3.03 0.455 4.58 16.1 17.78
G-2 5.49 73.2 58.5 3.138 3.808 3.177 1.76 3.23 0.7S8 18.8 16.1 20.51
G-3 4.14 101.8 77.2 3.748 4.954 3.856 1.59 3.27 1.21 36.6 16.0 21.27
C-4 3.15 118.6 101.2 4.711 5.462 5.180 1.49 3.19 1.82 86.9 15.9 20.35
G-5 2.20 169.4 144.8 6.162 9.066 7.950 1.39 3.27 3.12 293.0 15.9 • 21.43
B-l 7.94 17.7 15.5 0.895 0.895 0.938 5.66 0.98 0.175 1.13 6.1 11.35 9
C-l 7.83 15.7 13.2 0.989 0.989 0.999 7.37 0.98 0.150 1.23 5.2 17.95 48
B-l 7.63 17.2 13.9 1.267 1.336 1.416 8.32 0.86 0.160 2.77 5.3 12.30 136
13
3.0
Case Re U'/U, v ' / U . l a u f e r ( 19541 Pioe flow
3
A-5 2.9x10, 05 O Q — R e ^ . ' x l O
14
Case Re Fr R. Laufe. (1954) Pipe flow
4
• A-1 l . l x l O 0.17 600 — R.=10O0
O G-1 3 . Ox 0.45 1600
— R,=8700
A G-2 3.2x 0.80 1800
-9- G-3 3.3x 21 1600 © =2.3exo(-y*/R,)
$ G-4 3.2x 1.82 1500 4
© • A . f y ^ *
0 G-S 3.3x 3.12 1400
Grass(1971)
Re-6700 Ff=0.20
- 0
21
85
T) u'/U.=2.30exp(-5)
5 w7U."1.63exp(-l)
v7U.=l.Z7exp(-5)
©
i©
m
0.0
0.0 0.5 .0
=y/h
15
Author Grass (1971)
0.0 L
0
50 y 100
16
Case Re
B l i n c o e t a 1 . < 1 9 7 1 ) Open c n d n n e l
Re-0.9xlrf
Re-2.4xl0
T h e o r e t i c a l curve
R ; - 6 0 O frtxn u ' / U . - 2 . 3 e x o ( - § | and
R.-2OO0 U/U.=2.Sln(|).A(R.)
u'/U
0.0
0.01 0.1 | 1.0 o.oi |=y/h
= y / h
17
(a)
Fr=0.7i
Case U.vcm/s)
(a) frta.7i
Case Re
4
O A-3 1.18xl0
0> A-5 2 . 8 9 x ••
0.0 J L
0.0 0.1 0.2 0.5 1.0
y/h
5
^A-3 1.0
°- -«V/U?
(b)
10'
•V. 3-2
Case Fr Calculated
1.0 (b) G-1 0.46
R =1600
G-2 0.80
•So
4
.2 » 10
G-3 1.21
Case U.(cm/s) Fr G-4 1.82
G-1 2 152 0.46 O G-5 3.12
G-2 3. 138 0.80
y/h -Bremhorst et a l . (1973)
G-3 3. 743 1.21 0 4
R = (2.3 - 9.0) x 10 (Pipe)
G-4
G-S
4 . 711
£. 1 6 2
1.82 0 e
y/hs 3.12
'O - — Theore t i c a l curve 0.5
0.5
0.0
0.0
y/h, y/h.
o.o
0.0 0.5 .0
-uv/Uj;
•
0
§
C-l
0-1
48
136
0.989
1.267
•
0
O
B-l
C-l
3-1
B r e m h o r s t e t a 1 . I. 1 9 7 3 )
Re"I 2 . 3 - 9 . 0 1 x 1 0 (Pipe)
Lu et al.(1973)
Theoretical curves Ca1culated :
X s
Re .4.0xl0 '8oundary 4
laver)
R '600
0.5
0.0 0.0
0.0 0.0 0.1 0.2 0.5 ! .0
0.5 1.0 y/h
-uv/uf
18
(c)
Case 0 Re Fr Case h u; 1
Re Fr
"max Urn iC "max m U
P-l 8 .11 49.20 39 .64 1 .908 3.20 0. 445 H-1 7.92 51.82 40.95 1 .885 3.13 0 465
P-2 7,02 57.04 45 .76 2 ,303 3.18 0. 552 H-2 7.22 55.47 45.12 2 .336 3.13 0. 536
P-3 6 .16 63.16 52 .09 2,.558 3.17 0. 670 H-3 6.37 68.80 50.20 2 .601 3.11 0, 635
P-4 5 .35 66.78 54 .72 2 ,913 3.14 0, 723 H-4 5.57 72.68 57.57 2 .842 3.08 0, 779
P-5 5,. 73 76.81 55 96 3 ,246 3.09 0. 747 H-5 5.29 69.15 59.71 3 .119 3.04 0. 829
P-6 5 .21 81.69 61 .33 3. 444 3.12 0, 858 H-6 5.28 80.48 60.89 3..786 3.19 0, 846
P-7 4 .58 85.79 69. 89 3. 525 3.16 1. 043 W-7 4.24 97.80 75.16 3 .902 3.16 1. 166
P-8 4, 02 96.34 79. 91 3. 560 3.11 1. 273 H-8 3.43 107.49 91.57 4 .590 3.24 1. 568
P-9 3. 82 100.91 83. 94 3..704 3.17 1. 372 W-9 2.98 117.05 107.32 5 .038 3.28 1. 986
i 1.99;
i it"!
O : vV« IP/h
2
• .fl'33.3 (f.'fr)
F i g . 4.14
9/ Intensity of
surface wave
CI 8, —4*3.0
f1uctuations.
0.5 2.0
19
Spectra of Surface Wave
Data F i l e
Spectra" of Wall Pressure fluctuations u(I) and v(I) ,1 = 1 ,s|
B**J.LIIO*
N=4Q96
0.W5Ï
Lx= — Su(0)
if
Lx (Lx)
Normalized by
Lx (Lx)old
S(k) Printing,
for k-Lx
Lx X-Y plotting
by X-Y p l o t t e r or f i 1 i n g on MT
Read
Vt-T^-Xk-Lx)-
Lx modified by
H
Lx/Lx=(Y / t Y l n )
21
for smooth bed (Case A-1).
Fig. 4.24 Spectral distributions normalized by Lx
(higher-velocity group).
23
Fig. 4.27 Normalized spectra s (k)/Lx of v.
v
24
•I II L 1 11
Observed Data by Dual-hot-film Probe •
(a) Case Fr
Observed Data by Single-hot-film Probe 1.0 G-1 0.46
C a s e F r X Raichlen(1967), Fr=0.25 G-2 0.80
G-3 1.21
O G-1 0.46 A HcOuivey et al.(1969), Fr=0.36
G-4 1.82
-e-G-2 0.80
G-5 3.12
• G-3 1.21
0) G-4 1.82
-ft- G-5
1.0
1.0
26
Fig. 4.35 Distributions of turbulent energy dissipation eh/ul
27
Fig. 4.37 Turbulent energy diffusion Tr=q v/2i)\.
2
28
30.0
3
pn/u;
Rough Bed (Case 0 - 1 )
Smooth Bed (Case A - 1 )
£h/u2 O : f (Production)
P (Production)
• .' 6 (Dissipation)
£ (dissipation)
X : T (Diffusion)
T (Diffusion)
: R'-P-e-T
R'-P-e-T
Q. a,
•r—
- C c
3
A5 ° * ' CO O
C •!-
O +->
^ ••' - •r- u
ƒ* o +-> 3
03 ~ d
\ R— O
lb* = CU 5 ¬
0 £ CL
i 5 \
0 0 -* tr * -
J y °*
O
II
a / L«E- —
+ 8> °
1 * ! ! • | O
'
29
chapter 5 E x p e r i m e n t a l c o n s i d e r a t i o n on
mechanism o f t u r b u l e n c e - p r o d u c t i o n
1 ] 1 I 11 It 1 1 I I 1 1 III 1 ! M U M
PROÓAT-HTTJ OMIT- FUNCTIONS of U AND * FLUCTUÜ'ONS AUTNORS GUOT* 1 XJOLIN(L972}
» o A
• A
OPEN CRIINNEL SOU^DA-F LAYER
a
RE.l.TÏLU' R**'l. ilo'
A
i.
A
A «SI * A , . • M • \*
(GAUSS•IN)
^~^—ia. 0 0% 0 A
O ° 0°°A 0%
A
A
A
A
KREP]INF197]] OI1 CMRTNTL
: £ÏP<RI[R«NT| 1 CURVES OF $^ ANRJ F —
A RTV-t.800 m<3 7,100
A
A
A
A *
(GAUSSIAN) A ^RV(
1 1 1 11 11 1 1 1 1 1 ! 1 II 1 1 I I I ! ! !
+ 2 3
'1 10 y 10 10
1
0=u/u and v=v/v 1
AUTHORS LI^I(1971}
SROCT* {A, BOUAF- [01 S-MTN PICK
1
0
(-1.0 ! * *
1.0
+
(GAUSS IAN)
0.0
C.o or r A , 2
_ 1
ou 0.6 0.
V \ AJ A y/h
-2.0
v
-4.QH t v v: . ,
5-:» TO**)
-6.0
OOSER.-C A
CALCSTUF:
-3.0
Fig. 5.4
P r o b a b i l i t y density d i s t r i b u t i o n s
Fig. 5.3 D i s t r i b u t i o n s of the 4 - t h of Reynolds-stress f l u c t u a t i o n s .
and 5-th order moments.
30
""">"" " " toom
31
PROFCIHL I MY DUTRLHULTON OF REYNOLD. STRAIT
Smóotfi Ilril ( A - 1 ) WÏOIH LAD
0.1 |
WIL IFLFRI.ON/• JL. : CALEULATFD CURVAI (Y/TT>0.LT))
Qj) ILNIR F R.IT T («LI IL» T T«II.'(L BY . ,11 F. | VW.T KvlPtvvj •WNUL CURVTI
© ©
HRRNELAS SLR*>\\ OR(U|H«'D LIY UIL. (VRNT )-/N-0 JJO (*•.»,•,)
(J*<T1ON IMERECT ION
L-RNL
IVANT
0.05 ©LNTR<*((( ION ®
(IN-.RJI)
LAAOL
0 --0 01S
- 0 3 - 1 l l .. 1 I I I | L _ l 1 I 1 I I l lI Q
1 N • ».?9
' 5 10 H 15 20 1
" N /I,
•0.5 l—_1 1 L _ L _
- I — I — I — J 1 1 I L_L L_ I L__J _ J I
10 H 15 "lo
Fig. 5.10 Fractional contributions to Reynolds stress and
fraction of time occupied by each event (smooth bed)
an am
! il L I'
F i g . 5.11
F r a c t i o n a l c o n t r i b u t i o n s to Reynolds s t r e s s and
f r a c t i o n of time occupied by each event (rough bed).
AUTHOR,(<X*N CHANNEL) ORODKO, 1| «I.(1»H.)
Obitrvcd Calculaud ( O i l Chanavel)
(B) U4 (CALL D)
.0
©.
©
i ,i i ii
0.0 0.1 0.2 1.0 0.0 0.1 0.2 0.5 y/h 1.0
pi/si (I9'J)
RJFCTITM(?L O O A
S«EED (•! • •
r CALCULATED CURVES
\~ 0.2
0.0
J I Mil I I I II I u J I1 111II
1
10 IO- 10 J
CO
C
0=0000 O o O
•R-
+J
//// !! | If
1 -O CD
/'/// 1 •R—
////
s-
+J
1 C
j O 4-
fi I
O
1—
RO
A
O
C
O
t ƒ
1 O
*R—
L/V
CO
i 4->
<J
CO
O)
s_
-t->
CO
S- TO
<4-
"O V) CD
-a
r—
+->O
RO SZ
>> O
tu
u
F—- cc
RA
RA O
c_>
CO
LD
LN
CD
•I—
U_
34
OI
N s-
OPEN CHANNEL (BE.L.ULO*) SMOOTH OH L.HANNE L( HE* - LOW) P(v|w>H) CONDT TTONAL PROBDBK 11
P(Qlw>H)
CASE t[ CALCULATED ^7:
ET «1. (19/?) VOOU. UCD ( C*V» A )
1.6 O 0 /\ : B.ODKEY ET AL.[19M)
O
0 SMOOTH BOUNDARY LAYER(RE-I4.fi 10*)
(ü 1)6.2 ^: TU 1 WLLTMERTH( 1973)
RS4/RS2 . / 1'
\ \ 1 JF?CTLOC.
[VENT /
1.2 - \H-0 \
H^ — \
1.0
ip a a
1
1 '1
o.a 1
v/v'
0.6
Fig. 5.21 Conditional probability distributions of u and v .
I I Uliffi I 1 I I I Mil I 1 I I II
10 i<y 10 J
to
en
0.5 1.0
\ <
c n>
- -s
DJ
D J ua
r + ro
H at
^ —•
• c 1 n 11 1 1 111111
fD
t/i o-°sQ;fCt. y/h
0.5 1.0
0.5 1.0
O
•C^T' '1*
52 c •
i!
—
- o c
•r- O
+J •«-
> O en
2
Ii -• c
— - > O u
•r- O
+-> rO
(_> M—
OJ t-
•r-J 3
O) to
I
4_ OJ
O OJ
t-
S-
O
• r - aj
> -c
(O +->
jCNJ
=
CD C
CO •<-
i i;
cv
f
\
ji
1
•
1i
1 1 1 < '
36
Fig. 5.29
Fig. 5.27 Mean period f of
e Mean period (fe-fi ) of
n
37
AUTOCORRELATION OF
ID
O SMOOTH BED (C.IE A->)
• HOUUH BED (CASE 0-1)
He Ru(T)
° »* S • O o. o • •
o, I I I I I I I I I I oI I I I I I I I I I
0.0 0. 1 0.2 0.5 1.0 0.0 0.1 0.2 0.5 ,. 1.0
!• n
y/h ?.S 1.0
oo
EN
C-L 0.4& O
~i—I—i—i—I—r i—i—i—r ~ i — i — i — r
J>
cP
1 u2 I
DJ —1
=S
CL
-s
—I.
o" ƒ t> O p. h
IN % IÏ
% 8
In fB
-s
£ —1
ro
ro<
ro
ro
' — 1
r+
—<-
O Fig. 5.32 Mean period of bursting phenomenon
13 (higher-velocity group).
—-
n 3 J nn
^^(Noniialued)
0.0 O? (.o 2 0
10.0
J I I I I I I ' I
1.0
5.0 2
- -0 -1.0 0.0 1.0 2 0
L T / T )
°9< - (NorroaUzed)
1.0
°-° LogfT/ï) °- 5
Fig. 5.36 Normalized probability distributions of period T
Fig. 5.35 Probability distributions of 99.99 1 1 1 1 1 1 99.99
period of bursting phenomenon I I I'I I I I 1 1
I I I I I I I J 1 1 1 1 1
f
4
1.0
101 10
KIM ET EL. (1971) (VI,UIL METHOD)
R»-
1
o /
J I I I I I I I I 0.1 0.1
0.0
0.0 0.1 0.2
I I I I ?
, k 1 0 0.0 0.1 0.2 0.5 1.0 0.01
I i!i Ii I i I .2.0
IiI Q0 1
I I ' I I II I I I I I I I I
y/ti •3.0 -2.0 -,JD 0 1.0 3.0 -30 -2 0 1.0 O 1.0 2.0 3.0
L o g f j i / X«)
Fig. 5.37 Coefficient of variation o»/T of
period of bursting phenomenon. Fig. 5.38 Probability distributions of X 3
of bursting phenomenon.
1
(il) Instantaneous Velocity P r o f i l e ü and v Case A-1 R e ^ l . O x l O ' . F r = 0 . 1 5
t 23 4 8 Ss e c t - 23 670 t 23 853 , t 24 03 7,
41
I hu. Conditions, Ki»bl00. li,*300
y*-?l(l, y/h-0.7
42
i i : Li-
Vortex line,
a v e r s e mode.
3
DJ
•a
rt
ON
Pressure
measured at 0
O!
Flow R a p i d l y
A c c t l a n t f d to
I n i t i a l Condition
du« t o t n a rrjy and
HotMntuM T r a m f a r
fro* VorUx Syltti
-J'
tr
e
ro
n
Subitantlal Portion of
(nargy T r a n i f t r r c d to
Slow H o l t c u l a r T r a m f a r o f Cnaroy and
Homantin (and M a i l o r H a a t . I f A p p l l c a b f i ) ro
V o r t u S v i t a * D*g<narat«i f r o a 'Mowing' W i l l t o S y b t a y t r , MHARA I t
to T u r b g i a n c t dut to Accuvulatai u n t i l lr*akdown
Vortaa I r a a k - u p naar
WIL), and l l U H t a w i U l y
01sslooted t h a n by
VUeotHy
x-ut
Fig. 6.6 Conceptual model of the turbulence near the wall Fig. 6.8 A n-eddy model in the moving-axes.
during a ' c y c l i c ' process (after J.O. Hinze).
— : r m
L+AL n-eddy Model ( e=0.6 )
t
uJ + AtO
UJLIU „ O r ,
70.0 80.0
~fir y //•//////////////> X
Fig. 6.9 Secondary motion of a IT-eddy Fig. 6.11 Relative turbulence intensity predicted
(vortex-stretching). by n-eddy model.
I 0.45 O A-1 0 O
O • 34 MIMI
?
J
0.80
l.?l
0>
-O-
5 0- 1 9
C I 48
25 i. fl? (a b 1- 1 IJ6
1.0 3.1? 0 (j
w'/u'
Sneepluck u u i e j i f >
d|iproachini| itirtau'
v'/u' . :„ ^ . B a B
y llliinl ~V~ 0" 0
0.5 is
_ B
oo
n - E d d y Hodel ( e-0.6 J
» : 45*s6s78"
[dge ol sublayer ~» /.'LLLRIJLLLL t LLJ
/ 1 : 45-S8190-
0.0 1 1 1 1 1 1 1 1 1
30 0.0 0.5 1.0 0.0 1.0
J WETJL y/h 0.5 y/h
Fig. 6.10 Angle of eddy inclination Fig. 6.12 Comparison between n-eddy model and
to the wall (after J. Sternberg)
experimental r e s u l t s .
0=0.7, 3=0.01
% R.
E c k e t , T i a n n
21o)
A2A G r a s s
1.000 Laufer
30 y AO 50
Fig. 6 . 1 3 Mean velocity distributions when u ; is varied. F i g > 6 1 5 Turbulence i n t e n s i t i e s when [it is varied.
30 y AO
~n r
II I
20 -uv/. 30
y*=A0
403
(20),
y*=30
201
Sweep Ejection^
y =20 Interaction
40.5
(10)
10 r=i5
26.5
(5)
y*=io
r=5
0.0 0.5 1.0
r='/ T
3.0
'a
2.0
0.0' 10 20 30 40 50
Fig. 6.24 Bursting period estimated Fig. 6.26 Turbulence intensity estimated by combined model
by renewal model.
Combined model
A = 10 B=0.3
17.7-1
Height roughly y +
~ (10 - 15) y . ~ &
Quasi-laminar — Turbulent —
Interface Turbulent Nonturbulent
u u
Main Flow Main Flow Loa Profile
Log Profile
* Velocity
I n t e r a c t i o n F e e d b a c k S y s t e m >
1
I n d u c e d F l o w
Seepage
M a i n Flow
S t r e s s
E x c i t e D i s t u r b a n c e s
Fig. 8.2
Hydrodynamic interaction mechanism
between main and seepage flows.
Fig. 8.4 Turbulent permeable
open-channel flow f i e l d .
Willmarth et al.(196?)
«<SyU (Boundary Layer)
O 5.00
• 0.63
Corcos (1964)
J/0 (Pipe Flow)
e o.is
<D 0.20
O 0.75
0 11.2
50
n
51
'w
>> p
+J Ü
•1- (H
on 0
c o
01
+J <
C -5
•R— •
ö
0)
S- M
3 Q)
00 -u
on
0
S -) «s
s
-
O-
1 Q)
I — Si
1—
ro M-
3 O
01 C
on O
0) - r -
I — +->
C CJ
o e
•I- 3
on f-
cOJ ro
E _
•r- on
O ro
O
IT- ^3
.U
CO B
X
,
cn
•I—
U-
Fig. 8.12
Relationship between seepage velocity g m Fig. 8.13
and hydraulic gradient i .
e Relationship between RK and
52
3
O
ro
E
3
ro
O
CU
CU
a LO
CO
C A S E A
R U N
' • A . a
f r o m b\ from b t
poroul
C A S E C n (cm) strata E - 1 0.0672 5.5 0.184 0.156
F - 1 0.1228 5.5 0.096
E 1.08x10-' 0.172 0.39 0.083
1.2S 3
Table 8.3
Evaluated values of the induced s t r e s s .
C A S E h U'
hplh t.-ilh £ UJV. RK /
R U N (cm) V
(crn/sec)
E - 3 8.01 0.42 4.4x10* 0.63 3.90x10-' 5.152 0.334 5.768 0.942 0.567 0.497 1.761
F - 3 8.33 0.41 4.4x10' 0.58 7.50x10-' 5.421 0.270 4.479 4.900 0.541 0.389 1.958
53
chapter 9 T u r b u l e n t s t r u c t u r e w i t h transpiration
© ffi
® ®
Fig. 9.4 E f f e c t s o f t r a n s p i r a t i o n
on the wall shear s t r e s s .
54
55
Fig. 9.9 Fig. 9.10
Absolute production of turbulent Absolute variations of d i r e c t - d i s s i p a t i o n
energy with transpiration. and production by transpiration 3.
56
CO <— CTL CO
—• US CNJ R-.
L—O N CO UI CO
en ts
CO i — CO CO
O R- CN
Ull
• CNJ
O UD
ex er
CO UJ Cfi CN
•— ^ IN
o
a.
CM J- c
O
o o
M •4- -R-
4->
ro ro
4 - > S-
CN ra •!-
VO
i CD CO "O Q.
1/1
O O c
CN • i - ro
"1 t — S-
th t
drau
O ^ EG
O >vr-3
3=
•— C\J N
«5 i/l
CD CC CN -CR
u R** 00 -ER C~3
-Q
(O
<si
rO
i—
I
CNJ
I
m
I
Fig. 9.12 Mean velocity p r o f i l e
with t r a n s p i r a t i o n .
• (-0.092)
57
(a) (b)
U only
CASE H-1
RUN
1 o 0
4
4 0> 0. 0598
5 & 0. 1417
8 0 -0. 0410
9 •© -0. 0676
yu./y
CASE HR-1
u 8
RUN * . ' _ . .
30 ^.-i(yi.o.v.v-io)
MX .
10 0 0.616x10.148
11 • 0.802 • 0.198
12 A 0.973 • 0.249 *Jl
1 O 0
r.'
I • O.CHW
S « 0.1*17
a * -o.wio
9 » -0.0*76
1.0 10 100
58
Fig. 9.20 Distributions of Reynolds s t r e s s with t r a n s p i r a t i o n .
59
CASE
CASE H-1 RUN » .
vi ! O 0.0
3 0 0.033
0.0 4 <j> 0.061
0.073
5 - O 0.158
-G- 0 . 1 7 6
7 • -0.022
E «JJ - 0 . 0 3 7
8 C> - 0 . 0 3 2
-Q- - 0 . 0 6 0
9 -«--0.050
Andersen(1975)
Andersen(1975)
• 0.049
• 0.049
A -0.083
A -0.033
u7U,'CVe>p(-*)
5 J
^ 0.0
O 0.0
CP 0.073 0.073
-e- 0.177
0,= 3.0 0.177
V -0.036 è -0.036
-0.060 -O- - 0 . 0 6 0
Andersenf1975) Andersen(1975)
A -0.083 A -0.083
60
2.0
k (cm") 10'
1
y/h=0.05
S
•# 5
I
*
0.0
0.060
0-1 «2
-V
<t> 8 -0.041 V
- » 9 -0.068
1
: TD RMANS
. fomjla
«
10
10
10' •0° k(cm') 10' 10 2
Lxk
62
1.0 U (a) Macro-scale
V u
(b') Micro-scale
* 9
<Ü1
i 111 0.01J 1—I I I I 111 J : l l M I I
0.1
« 0.1 1.0
y/h
4 o
CASE H-1
RUN v+
• 1 0.0
0 4 0.060
-9- 5 0.142
0> 8 -0.041
-©•9 -0.068
0.01 r - j — i — i i i i 111 — i — i — i |
hO.OI 0.1
63