Open Hannel Flow Notes Answers
Open Hannel Flow Notes Answers
Section 1.2
Example.
The discharge in a channel with bottom width 3 m is 12 m3 s–1. If Manning’s 𝑛 is 0.013 m–1/3 s
and the streamwise slope is 1 in 200, find the normal depth if:
(a) the channel has vertical sides (i.e. rectangular channel);
(b) the channel is trapezoidal with side slopes 2H:1V.
𝑏 = 3 m (base width)
𝑄 = 12 m3 s −1
𝑛 = 0.013 m−1⁄3 s
𝑆 = 0.005
(a)
Discharge:
𝑄 = 𝑉𝐴
where, in normal flow,
1 2/3 1/2 𝑏ℎ ℎ
𝑉= 𝑅 𝑆 , 𝐴 = 𝑏ℎ, 𝑅ℎ = =
𝑛 ℎ 𝑏 + 2ℎ 1 + 2ℎ/𝑏
Hence,
1 𝑏ℎ5/3
𝑄= 𝑆 1/2
𝑛 (1 + 2ℎ/𝑏)2/3
Rearranging as an iterative formula for ℎ:
3/5
𝑛𝑄
ℎ=( ) (1 + 2ℎ/𝑏)2/5
𝑏√𝑆
Here, with lengths in metres,
ℎ = 0.8316 (1 + 2ℎ/3)2/5
Iteration (from, e.g., ℎ = 0.8316) gives
ℎ𝑛 = 1.024 m
(b) Geometry: trapezoidal cross-section with base width b, surface width 𝑏 + 2 × (2ℎ) and two
sloping side lengths √ℎ2 + (2ℎ)2 = ℎ√5.
𝑃 = 𝑏 + 2ℎ√5
Discharge:
1 2/3
𝑄 = 𝑉𝐴 = 𝑛 𝑅ℎ 𝑆 1/2 𝐴
Hence,
2/3
1 1 + 2ℎ/𝑏
𝑄 = ℎ2/3 ( ) 𝑆 1/2 ℎ𝑏(1 + 2ℎ/𝑏)
𝑛 1 + 2√5ℎ/𝑏
𝑛𝑄 (1 + 2ℎ/𝑏)5/3
= ℎ5/3
𝑏√𝑆 (1 + 2√5ℎ/𝑏)2/3
3/5
𝑛𝑄 (1 + 2√5ℎ/𝑏)2/5
ℎ=( )
𝑏√𝑆 1 + 2ℎ/𝑏
Example.
The discharge in a rectangular channel of width 6 m with Manning’s 𝑛 = 0.012 m−1⁄3 s is
24 m3 s–1. If the streamwise slope is 1 in 200 find:
(a) the normal depth;
(b) the Froude number at the normal depth;
(c) the critical depth.
𝑏 =6m
𝑛 = 0.012 m−1⁄3 s
𝑄 = 24 m3 s −1
𝑆 = 0.005
(a)
Discharge:
𝑄 = 𝑉𝐴
where, in normal flow,
1 2/3 1/2 𝑏ℎ ℎ
𝑉= 𝑅 𝑆 , 𝐴 = 𝑏ℎ, 𝑅ℎ = =
𝑛 ℎ 𝑏 + 2ℎ 1 + 2ℎ/𝑏
Hence,
1 𝑏ℎ5/3
𝑄= 𝑆 1/2
𝑛 (1 + 2ℎ/𝑏)2/3
or, rearranging as an iterative formula for h:
3/5
𝑛𝑄
ℎ=( ) (1 + 2ℎ/𝑏)2/5
𝑏√𝑆
Here, with lengths in metres,
ℎ = 0.7926 (1 + ℎ/3)2/5
Iteration (from, e.g., ℎ = 0.7926) gives
ℎ𝑛 = 0.8783 m
𝑉
Fr =
√𝑔ℎ
where
𝑄
𝑉= (in general)
𝐴
𝑄 𝑞
or 𝑉= = (for a rectangular channel; 𝑞 is the flow per unit width)
𝑏ℎ ℎ
2
(𝑞/ℎ)2 𝑞2
Fr = = 3
𝑔ℎ 𝑔ℎ
The normal depth is supercritical because, when ℎ = ℎ𝑛 , then Fr > 1 (part (b)).
Alternatively (and often more conveniently), the normal depth here is supercritical because
ℎ𝑛 < ℎ𝑐 ; so speed 𝑉 is larger, and depth ℎ is smaller in normal flow than critical flow, so that
Fr ≡ 𝑉/√𝑔ℎ must be greater than 1.
Example.
A 3-m wide channel carries a total discharge of 12 m3 s–1. Calculate:
(a) the critical depth;
(b) the minimum specific energy;
(c) the alternate depths when 𝐸 = 4 m.
𝑏 =3m
𝑄 = 12 m3 s −1
(a)
Discharge per unit width:
𝑄 12
𝑞= = = 4 m2 s −1
𝑏 3
Then, for a rectangular channel:
1/3 1/3
𝑞2 42
ℎ𝑐 = ( ) =( ) = 1.177 m
𝑔 9.81
(c) As 𝐸 > 𝐸𝑐 , there are two possible depths for a given specific energy.
𝑉2 𝑄 𝑞
𝐸 ≡ℎ+ where 𝑉= = (for a rectangular channel)
2𝑔 𝐴 ℎ
𝑞2
𝐸 ≡ℎ+
2𝑔ℎ2
Substituting values in metre-second units:
0.8155
4≡ℎ+
ℎ2
For the subcritical (slow, deep) solution, the first term, associated with potential energy,
dominates, so rearrange as:
0.8155
ℎ =4−
ℎ2
For the supercritical (fast, shallow) solution, the second term, associated with kinetic energy,
dominates, so rearrange as:
0.8155
ℎ=√
4−ℎ
Iteration (from, e.g., ℎ = 0) gives ℎ = 0.4814 m.
A long, wide channel has a slope of 1:1000, a Manning’s 𝑛 of 0.015 m–1/3 s and a discharge
of 3 m3 s–1 per metre width.
(b) Calculate the normal and critical depths.
(c) In a region of the channel the bed is raised by a height of 0.8 m over a length
sufficient for the flow to be parallel to the bed over this length. Determine the depths
upstream, downstream and over the raised bed, ignoring frictional losses. Sketch the
key features of the flow, indicating all hydraulic transitions caused by the bed rise.
(d) In the same channel, the bed is lowered by 0.8 m from its original level. Determine the
depths upstream, downstream and over the lowered bed, ignoring frictional losses.
Sketch the flow.
(a)
(i) Specific energy is the head relative to the local bed of the channel:
𝑉2
𝐸 =ℎ+
2𝑔
In critical conditions Fr = 1 and the specific energy is the minimum for that discharge.
(b)
𝑆 = 0.001
𝑛 = 0.015 m–1⁄3 s
𝑞 = 3 m2 s–1
Normal depth:
Discharge per unit width:
1 2/3 1/2
𝑞 = 𝑉ℎ, where 𝑉= 𝑅 𝑆 (Manning), 𝑅ℎ = ℎ (“wide” channel)
𝑛 ℎ
1 2/3 1/2
𝑞= ℎ 𝑆 ℎ
𝑛
ℎ5/3 √𝑆
𝑞=
𝑛
Critical depth:
1/3 1/3
𝑞2 32
ℎ𝑐 = ( ) =( ) = 0.9717 m
𝑔 9.81
(c) To determine the type of behaviour over the raised bed, compare the total head under critical
conditions (the minimum energy necessary to get over the weir at this flow rate) with that
available in the approach flow.
Critical
ℎ𝑐 = 0.9717 m
3
𝐸𝑐 = ℎ = 1.458 m
2 𝑐
𝑧𝑏 = 0.8 m
𝐻𝑐 = 𝑧𝑏 + 𝐸𝑐 = 2.258 m
Approach Flow
Because the channel is described as “long” it will have sufficient fetch to develop normal flow;
hence the approach-flow head is that for the normal depth (ℎ = 1.236 m):
𝑉𝑛2
𝐻𝑎 = 𝐸𝑎 = ℎ𝑛 +
2𝑔
𝑞2
= ℎ𝑛 +
2𝑔ℎ𝑛2
32
= 1.236 +
2 × 9.81 × 1.2362
= 1.536 m
At the normal depth the available head (𝐻𝑎 ) is less than the minimum required to get over the
bed rise (𝐻𝑐 ). Hence the water depth must increase upstream (“back up”), to raise the head
immediately upstream. Thus:
• critical conditions do occur;
• the total head in the vicinity is the critical head (𝐻 = 𝐻𝑐 = 2.258 m).
Over the raised bed there is a hydraulic transition, so the depth over this is critical: ℎ = ℎ𝑐 =
0.9717 m.
𝑞2
𝐻 =ℎ+
2𝑔ℎ2
0.4587
2.258 = ℎ +
ℎ2
0.4587
ℎ=√
2.258 − ℎ
Iteration (from, e.g., ℎ = 0) gives ℎ = 0.5127 m.
Since the preferred (i.e. normal) depth is subcritical, there must be a downstream hydraulic
jump. (A quick calculation shows that the upstream depth for this jump is greater than ℎ2 , so
there is indeed a length of GVF between the area of bed rise and the jump.)
Answer: depths upstream, over, downstream of the raised bed: 2.16 m, 0.972 m, 0.513 m.
(d) The flow does not require additional energy to pass a depressed section; hence, the total
head throughout is that supplied by the approach flow (𝐻 = 𝐻𝑎 = 1.536 m) and the flow
remains subcritical. The depths just upstream and downstream of the lowered section are those
in the approach flow; i.e. normal depth.
As bed height 𝑧𝑏 decreases, specific energy 𝐸 must increase to maintain the same total head.
In the lowered section:
𝐻 = 𝑧𝑏 + 𝐸
1.536 = −0.8 + 𝐸
𝐸 = 2.336 m
𝑞2
𝐸 =ℎ+
2𝑔ℎ2
0.4587
2.336 = ℎ +
ℎ2
(Note that this is the depth of the water column. The actual surface level here is
𝑧𝑠 = −0.8 + ℎ = 1.445 m
so the overall water level also rises in this section.)
normal RVF
normal
hn
hn
Answer: depths upstream, within, downstream of the lowered section: 1.24 m, 2.24 m, 1.24 m.
Example.
A long channel of rectangular cross-section with width 3.5 m and streamwise slope 1 in 800
carries a discharge of 15 m3 s–1. Manning’s n may be taken as 0.016 m–1/3 s. A broad-crested
weir of height 0.7 m is constructed at the centre of the channel. Determine:
(a) the depth far upstream of the weir;
(b) the depth just upstream of the weir;
(c) whether or not a region of supercritical gradually-varied flow exists downstream of the
weir.
𝑏 = 3.5 m
𝑆 = 0.00125
𝑄 = 15 m3 s −1
𝑛 = 0.016 m−1⁄3 s
𝑧𝑤𝑒𝑖𝑟 = 0.7 m
(a) The depth far upstream is normal since the channel is described as “long”. For normal flow
in a rectangular channel:
𝑄 = 𝑉𝐴
where:
1 2/3 1/2 𝑏ℎ ℎ
𝑉= 𝑅 𝑆 𝐴 = 𝑏ℎ 𝑅ℎ = =
𝑛 ℎ 𝑏 + 2ℎ 1 + 2ℎ/𝑏
Hence,
1 𝑏ℎ5/3
𝑄= 2/3
𝑆 1/2
𝑛 (1 + 2ℎ/𝑏)
or, rearranging as an iterative formula for ℎ:
3/5
𝑛𝑄
ℎ=( ) (1 + 2ℎ/𝑏)2/5
𝑏√𝑆
(b) To establish depths near the weir we need to know the flow behaviour at the weir. Compare
the energy in the approach flow with that under critical conditions.
Approach flow
Critical conditions
1/3
𝑞2 𝑄 15
ℎ𝑐 = ( ) 𝑞= = = 4.286 m2 s−1
𝑔 𝑏 3.5
1/3
4.2862
ℎ𝑐 = ( ) = 1.233 m
9.81
3 3
𝐸𝑐 = 2 ℎ𝑐 = 2 × 1.233 = 1.850 m
𝐻𝑐 = 𝑧𝑤𝑒𝑖𝑟 + 𝐸𝑐 = 0.7 + 1.850 = 2.550 m
Since the head required to flow over the weir (𝐻𝑐 = 2.550 m) exceeds that in the approach
flow (𝐻𝑎 = 2.252 m), the depth just upstream of the weir must increase and the flow back up.
The total head at any position in the vicinity of the weir is 𝐻 = 𝐻𝑐 = 2.550 m.
Just upstream and downstream of the weir (i.e. at undisturbed bed level):
𝑉2 𝑄
𝐻=𝐸 =ℎ+ 𝑉=
2𝑔 𝑏ℎ
𝑄2
𝐻 =ℎ+
2𝑔𝑏 2 ℎ2
0.9362
2.550 = ℎ + (*)
ℎ2
The depth just upstream is the deep, subcritical solution. Hence, rearrange as
0.9362
ℎ = 2.550 −
ℎ2
Iteration (from, e.g., ℎ = 2.550) gives
ℎ = 2.385 m
If the flow in the vicinity of the weir is unaffected by the hydraulic jump the flow goes smoothly
supercritical on the downstream side, with total head 𝐻 = 2.550 m (equation (*)).
Rearranging to get an iterative formula for the supercritical solution:
0.9362
ℎ=√
2.550 − ℎ
Denote by subscripts A and B respectively the conditions upstream and downstream of the
hydraulic jump. On the downstream side conditions may be assumed normal, since the channel
is “long” and hence there is sufficient fetch to develop the preferred depth:
ℎ𝐵 = 2.023 m
𝑉𝐵 = 2.118 m s −1 (from part (b))
𝑉𝐵 2.118
Fr𝐵 = = = 0.4754
√𝑔ℎ𝐵 √9.81 × 2.023
Hence, from the hydraulic-jump relation for the sequent depths:
ℎ𝐵 2.023
ℎ𝐴 = (−1 + √1 + 8Fr𝐵2 ) = (−1 + √1 + 8 × 0.47542 ) = 0.6835 m
2 2
Any gradually-varied supercritical flow downstream of the weir would increase in depth until
a hydraulic jump occurred (see the lectures on GVF). Since the depth downstream of the weir
is already greater than any sequent depth upstream of the hydraulic jump, no such increasing-
depth GVF is possible and the hydraulic jump must actually occur at (or just before) the
downstream end of the weir.
Example.
A reservoir has a plan area of 50 000 m2. The outflow passes over a broad-crested weir of width
8 m and discharge coefficient 0.9. Calculate:
(a) the discharge when the level in the reservoir is 0.6 m above the top of the weir;
(b) the time taken for the level of water in the reservoir to fall by 0.3 m.
(a)
Total head:
𝑉2
𝐻 = 𝑧𝑠 + = 𝑧𝑏 + 𝐸
2𝑔
where levels z can be measured relative to any convenient datum. Relative to the top of the
weir, assuming constant head, still water in the reservoir and critical conditions over the weir:
head upstream = head over weir
1/3
3 𝑞2
ℎ0 = ( )
2 𝑔
Hence, in ideal flow,
𝑞 = √(2/3)3 𝑔ℎ03
3/2
𝑄𝑖𝑑𝑒𝑎𝑙 = 𝑞𝑏 = 𝑏√(2/3)3 𝑔ℎ03 = 1.705𝑏ℎ0 (in metre-second units)
When ℎ0 = 0.6 m,
𝑄 = 12.28 × 0. 63/2 = 5.707 m3 s−1
Consider the change in volume of a tank, water surface area 𝐴𝑤𝑠 . When the water level changes
by dℎ the change in volume is
d(𝑣𝑜𝑙𝑢𝑚𝑒) = 𝐴𝑤𝑠 dℎ
Hence, by continuity,
Separating variables,
dℎ
50000 = −12.28 d𝑡
ℎ3/2
Apply boundary conditions ℎ = 0.6 when 𝑡 = 0 and ℎ = 0.3 when 𝑡 = 𝑇, and integrate:
0.3 𝑇
50000 ∫ ℎ−3/2 dℎ = −12.28 ∫ d𝑡
0.6 0
0.3
ℎ−1/2
50000 [ ] = −12.28𝑇
(−1/2) 0.6
50000 × 2 1 1
( − )=𝑇
12.28 √0.3 √0.6
𝑇 = 4355 s
A venturi flume is placed near the middle of a long rectangular channel with Manning’s 𝑛 =
0.012 m−1⁄3 s. The channel has a width of 5 m, a discharge of 12.5 m3 s–1 and a slope of
1:2500.
(a) Determine the critical depth and the normal depth in the main channel.
(b) Determine the venturi flume width which will just make the flow critical at the
contraction.
(c) If the contraction width is 2 m find the depths just upstream, downstream and at the
throat of the venturi flume (neglecting friction in this short section).
(d) Sketch the surface profile.
𝑛 = 0.012 m−1⁄3 s
𝑏 = 5 m (main channel)
𝑄 = 12.5 m3 s−1
𝑆 = 4 × 10−4
(a)
In the main channel,
𝑄 12.5
𝑞= = = 2.5 m2 s −1
𝑏 5
Critical Depth
1/3 1/3
𝑞2 2. 52
ℎ𝑐 = ( ) =( ) = 0.8605 m
𝑔 9.81
Normal Depth
𝑄 = 𝑉𝐴
where:
1 2/3 1/2 𝑏ℎ ℎ
𝑉= 𝑅 𝑆 , 𝐴 = 𝑏ℎ, 𝑅ℎ = =
𝑛 ℎ 𝑏 + 2ℎ 1 + 2ℎ/𝑏
Hence,
1 𝑏ℎ5/3
𝑄= 𝑆 1/2
𝑛 (1 + 2ℎ/𝑏)2/3
or, rearranging as an iterative formula for ℎ:
3/5
𝑛𝑄
ℎ=( ) (1 + 2ℎ/𝑏)2/5
𝑏√𝑆
(b) The flow will just go critical if the head in the throat (𝐻𝑐 ) is exactly equal to that in the
approach flow (𝐻𝑎 ). Measure heads relative to the bed of the channel in the vicinity of the
venturi.
Critical Head
1/3
3 2
3 𝑞𝑚 𝑄
𝐻𝑐 = 𝐸𝑐 = ℎ𝑐 = ( ) , where 𝑞𝑚 =
2 2 𝑔 𝑏𝑚
(Note that the critical depth is different at the throat to that in the main channel, due to the
narrower width.)
1/3
3 𝑄2
𝐻𝑐 = ( 2 )
2 𝑔𝑏𝑚
Approach Flow
The approach flow is normal, since the channel is “long”. Hence,
ℎ𝑎 = 1.546 m
𝑄 12.5
𝑉𝑎 = = = 1.617 m s−1
𝑏ℎ𝑎 5 × 1.546
𝑉𝑎2 1.6172
𝐻𝑎 = 𝐸𝑎 = ℎ𝑎 + = 1.546 + = 1.679 m
2𝑔 2 × 9.81
12.5
𝑏𝑚 = = 3.370 m
2 3/2
(3 × 1.679) √9.81
At the throat,
𝑄 12.5
𝑞𝑚 = = = 6.25 m2 s −1
𝑏𝑚 2
2 1/3 1/3
3 𝑞𝑚 6.252
ℎ𝑐 = ( ) =( ) = 1.585 m
2 𝑔 9.81
3 3
𝐸𝑐 = ℎ = × 1.585 = 2.378 m
2 𝑐 2
𝐻𝑐 = 𝑧𝑏 + 𝐸𝑐 = 0 + 2.378 = 2.378 m
The head throughout the venturi will be the critical head (𝐻 = 𝐻𝑐 = 2.378 m).
𝑄2
𝐻 =ℎ+
2𝑔𝑏 2 ℎ2
At the throat the depth will be the critical depth there; i.e. ℎ = ℎ𝑐 = 1.585 m.
Upstream
Rearrange for the deep, subcritical solution:
0.3186
ℎ = 2.378 −
ℎ2
Iteration (from, e.g., ℎ = 2.378) gives ℎ = 2.319 m.
Downstream
Rearrange for the shallow, supercritical solution:
0.3186
ℎ=√
2.378 − ℎ
WATER PROFILE
critical
Example.
The water depth upstream of a sluice gate is 0.8 m and the depth just downstream (at the vena
contracta) is 0.2 m. Calculate:
(a) the discharge per unit width;
(b) the Froude numbers upstream and downstream.
ℎ1 = 0.8 m
ℎ2 = 0.2 m
(a) Assuming total head the same on either side of the gate:
𝑉12 𝑉22
𝑧𝑠1 + = 𝑧𝑠2 +
2𝑔 2𝑔
Substituting 𝑧𝑠 = ℎ and 𝑉 = 𝑞 ⁄ℎ:
𝑞2 𝑞2
ℎ1 + = ℎ2 +
2𝑔ℎ12 2𝑔ℎ22
From the given data, in metre-second units:
0.8 + 0.0796𝑞 2 = 0.2 + 1.2742𝑞 2
0.6 = 1.1946𝑞 2
𝑞 = 0.7087 m2 s −1
Example.
A sluice gate controls the flow in a channel of width 2 m. If the discharge is 0.5 m3 s–1 and the
upstream water depth is 1.5 m, calculate the downstream depth and velocity.
𝑏 =2m
𝑄 = 0.5 m3 s −1
ℎ1 = 1.5 m
Use upstream conditions to get total head. Then, assuming no losses, find the supercritical
flow with the same head.
𝑄2 3.186 × 10−3
𝐻 =ℎ+ =ℎ+
2𝑔𝑏 2 ℎ2 ℎ2
Hence,
3.186 × 10−3
1.501 = ℎ +
ℎ2
3.186 × 10−3
ℎ=√
1.501 − ℎ
Iteration (from, e.g., ℎ = 0) gives
ℎ2 = 0.04681 m
𝑄 0.5
𝑉2 = = = 5.341 m s−1
𝑏ℎ2 2 × 0.04681
(a) Find the normal depths in both sections and show that normal flow is supercritical on
the spillway and subcritical on the apron.
(b) Baffle blocks are placed a short distance downstream of the slope transition to provoke
a hydraulic jump. Assuming that flow is normal on both the spillway and downstream
of the hydraulic jump, calculate the force per metre width of channel that the blocks
must impart.
𝑆1 = 1⁄30; 𝑆2 = 1⁄1000
𝑞 = 0.8 m2 s −1
𝑛 = 0.014 m−1⁄3 s
ℎ5/3 √𝑆
𝑞=
𝑛
𝑛𝑞 3/5
ℎ=( )
√𝑆
For the two slopes this gives
ℎ1 = 0.1874 m
ℎ2 = 0.5365 m
For subcritical/supercritical:
Method 1 (use critical depth)
The critical depth is
1/3 1/3
𝑞2 0. 82
ℎ𝑐 = ( ) =( ) = 0.4026 m
𝑔 9.81
(b)
The corresponding velocities are deduced from 𝑉 = 𝑞/ℎ, whence:
𝑉1 = 4.269 m s −1
𝑉2 = 1.491 m s −1
Let 𝑓 be the magnitude of the force per unit width exerted by the fluid on the blocks and, by
reaction, the blocks on the fluid, which is clearly acts in the upstream direction.
1
or ρ𝑔ℎ × ℎ (per unit width)
2
Example.
A downward step of height 0.5 m causes a hydraulic
jump in a wide channel when the depth and velocity
of the flow upstream are 0.5 m and 10 m s–1,
respectively.
h2
(a) Find the downstream depth.
h1
(b) Find the head lost in the jump.
(a) The downstream depth can be deduced from the momentum principle if the reaction force
from the step is known. The approximation is that this is the same as would occur if it were in
equilibrium with a hydrostatic pressure distribution here.
5
1 − ℎ22 = 1.019( − 10)
ℎ2
5.095
11.19 = ℎ22 +
ℎ2
Rearrange for the deep, subcritical solution:
5.095
ℎ2 = √11.19 −
ℎ2
A long, wide channel has a slope of 1:2747 with a Manning’s n of 0.015 m–1/3 s. It carries a
discharge of 2.5 m3 s–1 per metre width, and there is a free overfall at the downstream end. An
undershot sluice is placed a certain distance upstream of the free overfall which determines the
nature of the flow between sluice and overfall. The depth just downstream of the sluice is 0.5 m.
(b) Sketch, with explanation, the two possible gradually-varied flows between sluice and
overfall.
(c) Calculate the particular distance between sluice and overfall which determines the
boundary between these two flows. Use one step in the gradually-varied-flow equation.
(a)
Critical Depth
1/3 1/3
𝑞2 2. 52
ℎ𝑐 = ( ) =( ) = 0.8605 m
𝑔 9.81
Normal Depth
1 2/3 1/2
𝑞 = 𝑉ℎ, where 𝑉= 𝑅 𝑆 , (Manning), 𝑅ℎ = ℎ (“wide” channel)
𝑛 ℎ 0
1 2/3 1/2
𝑞= ℎ 𝑆0 ℎ
𝑛
ℎ5/3 √𝑆0
𝑞=
𝑛
3/5 3/5
𝑛𝑞 0.015 × 2.5
ℎ=( ) =( ) = 1.500 m
√𝑆0 √1/2747
(b) The depth just downstream of the sluice is supercritical (0.5 m < ℎ𝑐 ). However, the
preferred depth is subcritical (ℎ𝑛 > ℎ𝑐 ). Hence, if the channel is long enough then there will
be a downstream hydraulic jump, with the flow depth then decreasing to pass through critical
again near the overfall.
hydraulic jump
normal
supercritical critical
normal
supercritical critical
(c) As the channel shortens, the depth change across the hydraulic jump diminishes. The
boundary between the two possible flow behaviours occurs when the supercritical GVF just
reaches critical depth at the overfall (i.e. the limiting depth change across the hydraulic jump
is zero).
As the flow is supercritical, integrate the GVF equation forward from the downstream side of
the sluice gate (where ℎ = 0.5 m) to the overfall (where ℎ = ℎ𝑐 = 0.8605 m). Use 1 step.
GVF equation:
dℎ 𝑆0 − 𝑆𝑓
=
d𝑥 1 − Fr 2
For the direct-step method invert the GVF equation:
d𝑥 1 − Fr 2 d𝑥
= and Δ𝑥 ≈ ( ) Δℎ
dℎ 𝑆0 − 𝑆𝑓 dℎ
For the working, write the derivative as a function of h; (all lengths in metres).
𝑉 𝑞 𝑞2 0.6371
Fr = = Fr 2 = =
√𝑔ℎ √𝑔ℎ3 𝑔ℎ3 ℎ3
𝑛𝑞 2 1.406 × 10−3
𝑆𝑓 = ( ) =
ℎ5/3 ℎ10/3
Δℎ = 0.8605 − 0.5 = 0.3605
d𝑥
𝑖 ℎ𝑖 𝑥𝑖 ℎ𝑚𝑖𝑑 ( ) Δ𝑥
dℎ 𝑚𝑖𝑑
0 0.5 0
0.6803 217.1 78.26
1 0.8605 78.26
(a) Assuming no losses at the sluice, find the volume flow rate, Q.
(c) Compute the distance from the sluice gate to the hydraulic jump, assuming normal
depth downstream of the jump. Use two steps in the gradually-varied-flow equation.
𝑏 = 2.5 m
ℎ1 = 1.8 m
ℎ2 = 0.3 m
𝑛 = 0.012 m−1⁄3 s
𝑆0 = 0.002
(a) Assuming the same total head on either side of the gate:
𝑉12 𝑉22
𝑧𝑠1 + = 𝑧𝑠2 +
2𝑔 2𝑔
𝑞2 𝑞2
ℎ1 + = ℎ2 +
2𝑔ℎ12 2𝑔ℎ22
𝑞2 1 1
ℎ1 − ℎ2 = ( 2 − 2)
2𝑔 ℎ2 ℎ1
Substituting values:
1.5 = 0.5506𝑞 2
Hence, the flow per unit width is
1.5
𝑞=√ = 1.651 m2 s −1
0.5506
and the total discharge is
𝑄 = 𝑞𝑏 = 1.651 × 2.5 = 4.128 m3 s−1
(b)
Normal depth
𝑄 = 𝑉𝐴
Critical depth
1/3 1/3
𝑞2 1.6512
ℎ𝑐 = ( ) =( ) = 0.6525 m
𝑔 9.81
(c) The depth just upstream of the jump is the sequent depth to the normal depth:
ℎ𝑛 = 0.7389 m
𝑄 4.128
𝑉𝑛 = = = 2.235 m s−1
𝑏ℎ𝑛 2.5 × 0.7389
𝑉𝑛 2.235
Fr𝑛 = = = 0.8301
√𝑔ℎ𝑛 √9.81 × 0.7389
so that the depth just upstream of the jump (call it ℎ𝐽 ) is
ℎ𝑛 0.7389
ℎ𝐽 = (−1 + √1 + 8Fr𝑛2 ) = (−1 + √1 + 8 × 0.83012 ) = 0.5734 m
2 2
We must therefore do a GVF calculation from just downstream of the sluice (where ℎ = 0.3 m)
to just upstream of the hydraulic jump (where ℎ = 0.5734 m).
GVF equation:
dℎ 𝑆0 − 𝑆𝑓
=
d𝑥 1 − Fr 2
For the direct-step method rewrite the GVF equation “the other way up”:
For the working, write the derivative as a function of h; (all lengths in metres).
𝑉 𝑄/𝑏 (𝑄/𝑏)2 0.2779
Fr = = 2
Fr = =
√𝑔ℎ √𝑔ℎ3 𝑔ℎ3 ℎ3
𝑛𝑄 2 4/3 −4
(1 + 0.8ℎ)4/3
𝑆𝑓 = ( 5/3 ) (1 + 2ℎ/𝑏) = 3.926 × 10
𝑏ℎ ℎ10/3
0.5734 − 0.3
Δℎ = = 0.1367
2
Working formulae:
d𝑥
Δ𝑥 = ( ) Δℎ
dℎ 𝑚𝑖𝑑
where
0.2779
d𝑥 1−
= ℎ3 , Δℎ = 0.1367
dℎ (1 + 0.8ℎ)4/3 −4
[20 − 3.926 × ] × 10
ℎ10/3
𝑖 ℎ𝑖 𝑥𝑖 ℎ𝑚𝑖𝑑 (d𝑥/dℎ)𝑚𝑖𝑑 Δ𝑥
0 0.3 0
0.3684 338.7 46.30
1 0.4367 46.30
0.5051 288.1 39.38
2 0.5734 85.68
(a) Calculate the normal depth and critical depth in the channel and show that the channel
is hydrodynamically “steep” at this flow rate.
(b) The depth of flow just downstream of the sluice is 0.4 m. Assuming no head losses at
the sluice calculate the depth just upstream of the sluice.
(c) Sketch the depth profile along the channel, indicating clearly any flow transitions
brought about by the sluice and indicating where water depth is increasing or
decreasing.
(d) Use 2 steps in the gradually-varied flow equation to determine how far upstream of the
sluice a hydraulic jump will occur.
𝑆0 = 0.003
𝑛 = 0.012 m−1⁄3 s
𝑞 = 2 m2 s −1
(a)
Normal depth
1 2/3 1/2
𝑞 = 𝑉ℎ = 𝑅 𝑆 ℎ where 𝑅ℎ = ℎ (wide channel)
𝑛 ℎ 0
𝑛𝑞
= ℎ5/3 (*)
√𝑆0
3/5
𝑛𝑞 0.012 × 2 3/5
ℎ𝑛 = ( ) =( ) = 0.6095 m
√𝑆0 √0.003
Critical Depth
1/3 1/3
𝑞2 22
ℎ𝑐 = ( ) =( ) = 0.7415 m
𝑔 9.81
The normal depth is smaller than the critical depth. Hence the normal flow is supercritical; i.e.
the channel is steep at this discharge.
(b) Assuming the same total head on either side of the gate:
𝑞2 𝑞2
ℎ1 + = ℎ2 +
2𝑔ℎ12 2𝑔ℎ22
22 22
ℎ1 + = 0.4 +
2𝑔ℎ12 2𝑔 × 0. 42
0.2039
ℎ1 + = 1.674
ℎ12
Answer: ℎ1 = 1.59 m.
(c) The depth increases in each of the GVF regions (S1 and S3) shown below.
S1
normal
hn
h1
h2 S3 normal
CP hn
(d) Upstream of the sluice there is GVF between the hydraulic jump and the sluice. This starts
at the subcritical sequent depth of the hydraulic jump (to be found below) and ends at the depth
ℎ1 found in part (b).
GVF equation:
dℎ 𝑆0 − 𝑆𝑓
=
d𝑥 1 − Fr 2
For the direct-step method rewrite the GVF equation “the other way up”:
d𝑥 1 − Fr 2 d𝑥
= and Δ𝑥 ≈ ( ) Δℎ
dℎ 𝑆0 − 𝑆𝑓 dℎ
For the working, write the derivative as a function of ℎ; (all lengths in metres).
𝑉 𝑞 𝑞2 0.4077
Fr = = Fr 2 = =
√𝑔ℎ √𝑔ℎ3 𝑔ℎ3 ℎ3
𝑛𝑞 2 5.76 × 10−4
𝑆𝑓 = ( ) =
ℎ5/3 ℎ10/3
Working formulae:
d𝑥
Δ𝑥 = ( ) Δℎ
dℎ 𝑚𝑖𝑑
where
0.4077
d𝑥 1−
= ℎ3 Δℎ =– 0.351
dℎ [30 − 5.76 ] × 10−4
ℎ10/3
𝑖 ℎ𝑖 𝑥𝑖 ℎ𝑚𝑖𝑑 (d𝑥/dℎ)𝑚𝑖𝑑 Δ𝑥
0 1.594 0
1.419 303.9 –106.7
1 1.243 –106.7
1.068 262.2 –92.0
2 0.892 –198.7
o
25
α = 25°
𝑛 = 0.014 m−1⁄3 s
ℎ = 0.35 m
(a)
1
sinα =
Fr
1 1
Fr = = = 2.366
sinα sin25°
(b)
𝑉
Fr ≡
√𝑔ℎ
(c) The critical slope is that at which the normal depth is equal to the critical depth.
An undershot sluice is placed in a channel with a horizontal bed covered by gravel with a
median diameter of 5 cm and density 2650 kg m–3. The flow rate is 4 m3 s–1 per metre width
∗
and initially the depth below the sluice is 0.5 m. Assuming a critical Shields parameter 𝜏𝑐𝑟𝑖𝑡 of
0.06 and friction coefficient 𝑐𝑓 of 0.01:
(a) find the depth just upstream of the sluice and show that the bed there is stationary;
(b) show that the bed below the sluice will erode and determine the depth of scour.
𝑑 = 0.05 m
ρ𝑠 = 2650 kg m−3 (𝑠 = 2.65)
𝑞 = 4 m2 s −1
τ∗𝑐𝑟𝑖𝑡 = 0.06
𝑐𝑓 = 0.01
(a) The bed is mobile if and only if the bed shear stress exceeds the critical shear stress for
incipient motion. (Alternatively, one could find the velocity at which incipent motion occurs.)
First find the critical shear stress for incipient motion. The local shear stress can be found by
finding velocities using open-channel flow theory and using the friction coefficient.
Upstream of the sluice one must have the same total head:
𝑉2 𝑞2
𝐻 = 𝑧𝑠 + =ℎ+
2𝑔 2𝑔ℎ2
but now we look for the subcritical (deep) solution:
𝑞2
ℎ=𝐻−
2𝑔ℎ2
Answer: depth upstream of sluice = 3.70 m; as demonstrated, the bed here is stationary.
The bed will continue to erode, the flow depth ℎ increasing and velocity 𝑉 decreasing until
τ𝑏 = τ𝑐𝑟𝑖𝑡
At this point,
1
𝑐𝑓 ( ρ𝑉 2 ) = τ𝑐𝑟𝑖𝑡
2
whence
2 τ𝑐𝑟𝑖𝑡 2 48.56
𝑉=√ =√ × = 3.116 m s −1
𝑐𝑓 ρ 0.01 1000
𝑞 4
ℎ= = = 1.284 m
𝑉 3.116
This is the depth of flow. Since the downstream water level is set (at 0.5 m above the
undisturbed bed) by the gate the depth of the scour hole is
1.284 − 0.5 = 0.784 m