Week 8
Week 8
Hydraulics
Prof. Mohammad Saud Afzal
Department of Civil Engineering
Gradually Varied Flow (GVF)
• The flow in a channel is termed GRADUALLY VARIED, if the
flow depth changes gradually over a large length of the
channel.
The cross- sectional shape, size and bed
slope are constant
• Assumptions
The channel is prismatic.
• Assuming 𝜶 = 𝟏, we get
𝑽𝟐
𝑯=𝒛+𝒚 +
𝟐𝒈
• Differentiating both the sides of the
above equation w.r.t x
𝒅𝑯 𝒅𝒛 𝒅𝒚 𝒅 𝑽𝟐
= + + (Eq. 1)
𝒅𝒙 𝒅𝒙 𝒅𝒙 𝒅𝒙 𝟐𝒈
Represents the water surface
slope w.r.t the channel bed
Represents energy slope
Represents bottom slope
𝒅𝑯
= −𝑺𝒇 𝒅𝒛
𝒅𝒙 = −𝑺𝟎
𝒅𝒙
• Further,
𝒅 𝑽𝟐 𝒅 𝑸𝟐 𝒅𝒚
=
𝒅𝒙 𝟐𝒈 𝒅𝒚 𝟐𝒈𝑨𝟐 𝒅𝒙
or
𝒅 𝑽𝟐 −𝑸𝟐 𝒅𝑨 𝒅𝒚
=
𝒅𝒙 𝟐𝒈 𝒈𝑨𝟑 𝒅𝒚 𝒅𝒙
𝒅𝑨
= 𝑻, where T is the top-
𝒅𝒚
width of the channel
• So we can rewrite Eq. 1 as
𝒅𝒚 𝑸𝟐 𝑻 𝒅𝒚
−𝑺𝒇 = −𝑺𝟎 + −
𝒅𝒙 𝒈𝑨𝟑 𝒅𝒙
𝑸𝟐 𝑻
or NOTE: = 𝑭𝒓 𝟐 , where 𝑭𝒓 is Froude Number
𝒈𝑨𝟑
𝒅𝒚 𝑺𝟎 −𝑺𝒇
= 𝑸𝟐 𝑻 Differential Equation of GVF
𝒅𝒙 𝟏− 𝟑
𝒈𝑨
Classification of Flow Profiles
• If 𝑸, 𝒏 and 𝑺𝟎 are fixed, then the normal depth 𝒚𝟎 and the
critical depth 𝒚𝒄 are fixed.
Depth obtained from uniform flow
equations
• Three possible relationships that may exist between 𝒚𝟎 and 𝒚𝒄
are:
𝒚 𝟎 > 𝒚𝒄
𝒚 𝟎 < 𝒚𝒄
𝒚 𝟎 = 𝒚𝒄
• Further, 𝒚𝟎 does not exist when:
Region 2 – Space between the top line and the next lower
line.
Solution: 𝒃 = 𝟏𝟎 𝒎 𝒚 = 𝟏. 𝟓 𝒎 𝑽 = 𝟏 𝒎/𝒔
𝑺𝟎 = 𝟏/𝟒𝟎𝟎𝟎 𝑺𝒇 = 𝟎. 𝟎𝟎𝟎𝟎𝟒
𝑨 = 𝒃 × 𝒚 = 𝟏𝟎 × 𝟏. 𝟓 = 𝟏𝟓 𝒎𝟐 𝑻 = 𝒃 = 𝟏𝟎 𝒎 𝑸 = 𝑨𝑽 = 𝟏𝟓 𝒎𝟑 /𝒔
𝟏
𝒅𝒚 𝑺𝟎 −𝑺𝒇 𝒅𝒚 𝟒𝟎𝟎𝟎
− 𝟎.𝟎𝟎𝟎𝟎𝟒 𝒅𝒚
= 𝑸𝟐 𝑻
= 𝟏𝟓𝟐 ×𝟏𝟎
= 𝟐. 𝟐𝟓 × 𝟏𝟎−𝟒
𝒅𝒙 𝟏− 𝟑 𝒅𝒙 𝟏−
𝒅𝒙
𝒈𝑨 𝟗.𝟖𝟏×𝟏𝟓𝟑
Problem- 2
• A rectangular channel with a bottom width of 4 m and a bottom slope of 0.0008 has a
discharge of 1.5 m3/s. In a gradually varied flow in this channel, the depth at a certain
location is found to be 0.30 m. Assuming Manning’s n = 0.016, determine the type of
GVF profile.
Solution: 𝒃=𝟒𝒎 𝒚 = 𝟎. 𝟎𝟑 𝒎 𝐐 = 𝟏. 𝟓 𝒎𝟑 /𝒔
𝑺𝟎 = 𝟎. 𝟎𝟎𝟎𝟖 𝐧 = 𝟎. 𝟎𝟏𝟔
𝟏/𝟑 𝟏/𝟑
𝑸 𝟏.𝟓 𝒎𝟑 𝒒𝟐 𝟎.𝟑𝟕𝟓𝟐
Now, = = 𝟎. 𝟎𝟑𝟕𝟓 /𝒎 𝒚𝒄 = =
𝒃 𝟒 𝒔 𝒈 𝟗.𝟖𝟏
= 𝟎. 𝟐𝟒𝟑𝒎
𝟏 𝟏/𝟐 𝟏 𝟒𝒚𝟎 𝟐/𝟑
Now, 𝑸 = 𝑨𝑹𝟐/𝟑 𝑺𝟎 𝟏. 𝟓 = 𝟒 × 𝒚𝟎 (𝟎. 𝟎𝟎𝟎𝟖)𝟏/𝟑
𝒏 𝟎.𝟎𝟏𝟔 𝟒+𝟐𝒚𝟎
𝟐/𝟑
𝟒 𝒚𝟎 𝒚𝟎
𝟏. 𝟓 = × 𝟒𝟐/𝟑 × (𝟎. 𝟎𝟎𝟎𝟖)𝟏/𝟐
𝟎.𝟎𝟏𝟔 (𝟒+𝟐𝒚𝟎 )𝟐/𝟑
𝟐/𝟑
𝒚𝟎
= 𝟎. 𝟎𝟖𝟒𝟐
(𝟒+𝟐𝒚𝟎 )𝟐/𝟑
From trial and error
𝒚𝟎 = 𝟎. 𝟒𝟐𝟔 𝒎
𝒚𝟎 > 𝒚𝒄 (Mild slope)
Also 𝒚 𝟎 > 𝒚 > 𝒚𝒄
𝑴𝟐
Class Question
A wide rectangular channel has a Manning’s coefficient of 0.018. For a discharge
𝒎𝟑
intensity of𝟏. 𝟓 /𝒎, identify the possible types of gradually varied flow profiles
𝒔
produced in the following break in the grade of the channel.
𝑺𝟎𝟏 = 𝟎. 𝟎𝟎𝟎𝟒 𝒂𝒏𝒅 𝑺𝟎𝟐 = 𝟎. 𝟎𝟏𝟔
𝒎𝟑
Solution: Discharge intensity q = 1.5 /𝒎
𝒔
𝟏/𝟑 𝟏/𝟑
Critical depth 𝒚𝒄 = 𝒒𝟐 Τ𝒈 = 𝟐
𝟏. 𝟓 Τ𝟗. 𝟖𝟏 = 𝟎. 𝟔𝟏𝟐 𝒎
• Neglect any wall shear stress τw, within relatively short segment between the
sections
Where
y 21b y1
F1 pc1 A1 pc1
2 2
y 2 2b y2
F2 pc 2 A2 pc 2
2 2
b is the channel width
Hydraulic Jump Derivation
• Momentum equation can be written as
y 21 y 2 2 V1 y1
(V2 V1 ) Eq. 19
2 2 g
• Conservation of mass ( continuity ) gives
y 21 y 2 2 V1 y1 V1 y1 V12 y1
( V1 ) ( y1 y2 ) Eq. 21b
2 2 g y2 gy2
Hydraulic Jump Derivation
y2 2 y2
( ) ( ) 2 Fr1 0
2
Eq. 21c
y1 y1
Where Fr1 is upstream Froude number
y2 1
(1 (1 8 Fr12 )
y1 2
Hydraulic Jump Derivation
• Solution with minus sign is neglected ??, Thus
y2 1
(1 (1 8 Fr12 ) Eq. 22
y1 2
• We can also obtain 𝒉𝑳 / 𝒚𝟏 by using Eq. 21
• The result is
hL y2 Fr12 y1 2 Eq. 23
1 [1 ( ) ]
y1 y1 2 y2
Hydraulic Jump Derivation
• This means that y2/y1 cannot be less than 1 and Froude number upstream Fr1 is
always greater than 1 for hydraulic jump to take place.
V12 V22
y1 y2
2g 2g
V12 V22
y1 1 y 2 1
2 gy1 2 gy 2
Since V2
F 2 Froude number y1 1 F22 2 2 F22
gy
y 2 1 F1 2 2 F12
2
Also Q2 Q2
F 2
1
2
and F 2
2
2
B g y13 B g y 23
y2 1 1
(1 (1 8 Fr1 ) (1 (1 8 * 3.92 2 ) 5.07
2
y1 2 2
hL 0.671 m
Class Question
1) Prove that energy loss in a hydraulic jump occurring in a rectangular channel is
( y2 y1 ) 3
hL Eq. 24
4 y1 y2
The loss of mechanical energy that takes place in a hydraulic jump is calculated
by the application of energy equation (Bernoulli’s equation). If loss of total
head in the pump is 𝒉𝑳 , then we can write by Bernoulli’s equation neglecting
the slope of the channel.
𝒚𝟏 + 𝑽𝟐𝟏 Τ𝟐𝒈 = 𝒚𝟐 + 𝑽𝟐𝟐 Τ𝟐𝒈 + 𝒉𝑳
𝒒𝟐 𝟏 𝟏
𝒉𝑳 = 𝒚𝟏 − 𝒚𝟐 + 𝟐 − 𝟐 𝒒 = 𝑽𝟏 𝒚 𝟏 = 𝑽𝟐 𝒚 𝟐
𝟐𝒈 𝒚𝟏 𝒚𝟐
𝒒 𝑽𝟏
From Eq 21.c we are putting 𝑽𝟏 = (𝑭𝒓𝟏 = )
𝒚𝟏 𝒈𝒚𝟏
𝟐
𝒚𝟏 𝒚𝟐𝟐 + 𝒚𝟐𝟏 𝒚𝟐 𝒒𝟐 𝟐𝒒
= 𝒚𝟏 𝒚𝟐𝟐 + 𝒚𝟐𝟏 𝒚𝟐 − =𝟎
𝟒 𝟐𝒈 𝒈
𝒚𝟏 𝒚𝟐𝟐 + 𝒚𝟐𝟏 𝒚𝟐 𝟏 𝟏
𝒉𝑳 = 𝒚𝟏 − 𝒚𝟐 + 𝟐− 𝟐
𝟒 𝒚𝟏 𝒚𝟐
Solution:
F1 10 .0 and E L 3.20 m
y2
y1
3
Energy loss EL
4 y1 y 2
E L y 2 y1 1 3.20 13.651 1
3 3
37.08
y1 4 y 2 y1 y1 4 13.651
3.20
(i) y1 =depth before the jump = 0.0863 m
37.08
y 2 =depth after the jump =13.651×0.0863= 1.178 m
V1 V1
(ii) F1 10.0 V1 9.201 m / s
gy1 9.81 0.0863
Discharge intensity q V y
1 1 9.201 0.0863 0.7941 m 3
/s/m
V2 q 0.7941
(iii) Froude number after the jump F2 0.1983
gy2 y 2 gy2 1.178 9.81 1.178
Class Question
A rectangular channel has a width of 1.8 m and carries a discharge of 1.8 at a depth of
0.2 m. Calculate (a) the specific energy, (b) depth alternate to the existing depth and
(c) Froude numbers at the alternate depths.
Solution:
Let y1 0.20 m Existing depth
Area A1 By1 1.8 0.20 0.36 m 2
1.80
Velocity V1 Q A1 5.0 m 2 / s
0.36
E1 y1
V12
0.20
5.0
2
1.4742 m
(a) Specific energy 2g 2 9.81
(b) Let y 2 =depth alternate to y1
2 1.8
2
Then E2 E1 y2
V 2
1.4742 y2 1.4742, as V1 A1 V2 A2
2g 2 9.81 1.82 y 22
By trial and error, y 2 1.45
Solution:
V2
q 2.5 Y1
V1 10.0 m / s V1
y1 0.25
V1 10.0
Initial Froude number F1 6.386
gy1 9.81 0.25
(i) The sequent depth ratio y 2 y1 is given by y2 1
1 1 8F12
y1 2
1 1 8 6.386
y2 1 2
0.25 2
y 2 y1 3 2.136 0.2503
EL 3.141 m
4 y1 y 2 4 2.136 0.250
Class Question
A hydraulic jump occur in a horizontal triangular channel. If the sequent depths in
this channel are 0.60 m and 1.20 m respectively, estimate (i) the flow rate, (ii)
Froude number at the beginning and end of the jump and (iii) energy loss in the
jump.
Solution:
(i) Consider a triangular channel of side slope
m=1
m horizontal: 1 vertical in fig (in the present y
case m=1)
P=pressure force= A y my m y 3 3
2 y
3 Q 2 Q 2
M=Momentum flux =
A my 2
For a hydraulic jump in horizontal, frictionless channel P1 M 1 P2 M 2
m y13 Q2 m y 23 Q2
2
3 my 1 3 m y 22
Q2 1 1
m y12
gm 3
y 2 y1
3
y 22 3
On simplifying Q 2
m 2
y13 3 1 2 y14
where
y2
g
3 2 1 y12
y1
y2
In the present problem m=1, 1.2 0.6 2.0
y1
Q 2 1 0.6 2 3 1 2 2
5
0.24192 Q 1.541 m 3 / s
g 3 2 1
2
Q
II. For triangular channel F as such
A gAT
21.541
2
6.222
2
2 2 2 F
9.81 1 0.6
Q T Q 2 my 2Q 1
F2 3 2 5
5
2 6
gA gm y gm y F1 2.494
A1 1 0.6 0.36 m
2 2
E L 0.6
4.281
2
1.2
1.070
2
1.534 1.258 0.276 m
2 9.81 2 9.81
Class Question
Water flows in a wide channel at q =10 𝒎𝟑 /(s.m) and 𝒚𝟏 = 1.25 m. If the flow
undergoes a hydraulic jump, compute (a) 𝒚𝟐 , (b) 𝑽𝟐 , (c) 𝑭𝒓𝟐 , (d) 𝒉𝒇 , (e) the percentage
dissipation, (f) the power dissipated per unit width, and (g) the temperature rise due
to dissipation if 𝑪𝒑 4200 J/(kg. K).
Solution: 𝒒 𝟏𝟎𝒎𝟑 /(s.m)
𝑽𝟏 = = = 𝟖. 𝟎 𝒎/𝒔
𝒚𝟏 𝟏. 𝟐𝟓 𝒎
The upstream Froude number is therefore
𝑽𝟏 𝟖. 𝟎
𝑭𝒓𝟏 = 𝟏/𝟐
= 𝟏/𝟐
= 𝟐. 𝟐𝟖𝟓
𝒈𝒚𝟏 𝟗. 𝟖𝟏(𝟏. 𝟐𝟓)
𝟐𝒚𝟐
This is a weak jump. The depth 𝒚𝟐 is obtained from = −𝟏 + (𝟏 + 𝟖𝑭𝒓𝟐𝟏 )𝟏/𝟐
𝒚𝟏
𝟐𝒚𝟐 𝟏
𝟐 𝟐
= −𝟏 + (𝟏 + 𝟖(𝟐. 𝟐𝟖𝟓) ) = 𝟓. 𝟓𝟒
𝒚𝟏
𝟏
𝒚𝟐 = 𝒚𝟏 𝟓. 𝟓𝟒 = 𝟑. 𝟒𝟔 𝒎
𝟐
𝑽𝟏 𝒚𝟏 𝟖. 𝟎(𝟏. 𝟐𝟓)
The downstream velocity is 𝑽𝟐 = = = 𝟐. 𝟖𝟗 𝒎/𝒔
𝒚𝟐 𝟑. 𝟒𝟔
𝑽𝟐 𝟐. 𝟖𝟗
The downstream Froude number is 𝑭𝒓𝟐 = 𝟏/𝟐
= 𝟏/𝟐
= 𝟎. 𝟒𝟗𝟔
𝒈𝒚𝟐 𝟗. 𝟖𝟏(𝟑. 𝟒𝟔)
(𝟑. 𝟒𝟔 − 𝟏. 𝟐𝟓)𝟑
As expected, 𝑭𝒓𝟐 is subcritical, the dissipation loss is 𝒉𝒇 = = 𝟎. 𝟔𝟐𝟓 𝒎
𝟒(𝟑. 𝟒𝟔)(𝟏. 𝟐𝟓)
𝑵 𝒎𝟑
Power= 𝝆𝒈𝒒𝒉𝒇 = 𝟗𝟖𝟎𝟎 𝟑 𝟏𝟎 𝟎. 𝟔𝟐𝟓 𝒎 = 𝟔𝟏. 𝟑 𝒌𝑾/𝒎
𝒎 𝒔.𝒎
Finally, the mass flow rate is 𝒎ሶ = 𝝆𝒒 (1000 kg/𝒎𝟑 )[10 𝒎𝟑 /(s.m)] 10,000 kg/(s m), and the
temperature rise from the steady flow energy equation is
ሶ 𝒑 ∆𝑻
Power dissipated = 𝒎𝒄
𝑾 𝑱
𝟔𝟏𝟑𝟎𝟎 = 10,000 kg/(s m) 𝟒𝟐𝟎𝟎 . 𝑲 ∆𝑻
𝒎 𝒌𝒈