Gradually-Varied Flow (GVF) AUTUMN 2020 3.1 Normal Flow Vs Gradually-Varied Flow
Gradually-Varied Flow (GVF) AUTUMN 2020 3.1 Normal Flow Vs Gradually-Varied Flow
Gradually-Varied Flow (GVF) AUTUMN 2020 3.1 Normal Flow Vs Gradually-Varied Flow
EGL (en
V 2 /2g
ergy gra
de line)
Friction
slope Sf
h
Geome
tric slop
e S0
In normal flow the downslope component of weight balances bed friction. As a result, the water
depth ℎ and velocity 𝑉 are constant and the total-head line (or energy grade line) is parallel to
both the water surface and the channel bed; i.e. the friction slope 𝑆𝑓 is the same as the geometric
slope 𝑆0 .
sluice hydraulic
gate weir change
jump
of slope
As a result of disturbances due to hydraulic structures (weir, sluice, etc.) or changes to channel
width, slope or roughness the downslope component of weight may not locally balance bed
friction. As a result, the friction slope 𝑆𝑓 and bed slope 𝑆0 will be different and the water depth
ℎ and velocity 𝑉 will change along the channel.
The gradually-varied-flow equation gives an expression for dℎ⁄d𝑥 and allows one to predict
the variation of water depth along the channel.
𝑉2 𝑉2
𝐻 = 𝑧𝑠 + = 𝑧𝑏 + ℎ + (1)
2𝑔 2𝑔
where 𝑧𝑠 is the level of the free surface and 𝑧𝑏 is the level of the bed.
Although not crucial, we make the small-slope assumption and make no distinction between
the vertical depth ℎ (which forms part of the total head) and that perpendicular to the bed,
ℎ cos θ (which is used to get the flow rate).
Differentiate (2):
d𝐻 d𝑧𝑏 d𝐸
= + (4)
d𝑥 d𝑥 d𝑥
Define:
d𝐻
= −𝑆𝑓 (5)
d𝑥
d𝑧𝑏
= −𝑆0 (6)
d𝑥
𝑆𝑓 is the downward slope of the energy grade line, or friction slope; (more about how this is
calculated later). 𝑆0 is the actual geometric slope. Then
d𝐸
= 𝑆0 − 𝑆𝑓 (7)
d𝑥
Thus, the specific energy only changes if there is a difference between the geometric and
friction slopes, i.e. between the rates at which gravity drives the flow and friction retards it.
Otherwise we would have normal flow, in which the depth and specific energy are constant.
Equations (5) and (7) are two forms of the gradually-varied-flow equation. The third, and most
common, form rewrites d𝐸 ⁄d𝑥 in terms of the rate of change of depth, dℎ⁄d𝑥 .
dℎ 𝑆0 − 𝑆𝑓
= (9)
d𝑥 1 − Fr 2
Since the flow (and hence the velocity profile and bed stress) vary only gradually with distance,
friction is primarily determined by the local bulk velocity 𝑉. The local friction slope 𝑆𝑓 can
then be evaluated on the “quasi-uniform-flow” assumption that there is the same rate of energy
loss as in normal flow of the same depth; e.g. using Manning’s equation:
1 2/3 1/2
𝑉= 𝑅 𝑆
𝑛 ℎ 𝑓
Inverting for the friction slope:
𝑛2 𝑉 2 𝑛2 𝑄 2
𝑆𝑓 = 4/3
= 4/3 (10)
𝑅ℎ 𝑅ℎ 𝐴2
Both 𝐴 and 𝑅ℎ (which depend on the channel shape) should be written in terms of depth ℎ.
In general, the deeper the flow then the smaller the velocity and friction losses. Qualitatively,
greater depth lower velocity smaller 𝑆𝑓 ;
smaller depth higher velocity greater 𝑆𝑓 .
For any given discharge (but independent of slope) there is a particular critical depth ℎ𝑐 , at
which Fr = 1. For example, in a wide or rectangular channel, ℎ𝑐 = (𝑞 2 /𝑔)1/3.
For any given discharge and slope there is a particular normal depth ℎ𝑛 , associated with
uniform flow. For example, in a wide channel, ℎ𝑛 = (𝑛𝑞/√𝑆0 )3/5. The normal depth is the
depth to which the flow would tend given a long enough undisturbed fetch.
A downward slope is classed as steep if the normal depth is less than the critical depth (i.e. the
normal flow is supercritical) and mild if the normal depth is greater than the critical depth (i.e.
the normal flow is subcritical).1 Note that, in principle, if rarely in practice, a particular slope
could be described as either steep or mild depending on the flow rate.
In general, given the actual depth ℎ, critical depth ℎ𝑐 and normal depth ℎ𝑛 , simple inspection
of signs of numerator and denominator on the RHS of the GVF equation
dℎ 𝑆0 − 𝑆𝑓
=
d𝑥 1 − Fr 2
indicates whether depth is increasing or decreasing with distance. In the special case 𝑆𝑓 = 𝑆0
we have dℎ⁄d𝑥 = 0; i.e. normal flow. Otherwise:
• 𝑆0 − 𝑆𝑓 > 0 if and only if ℎ is greater than normal depth; (slower, hence less friction)
• 1 − Fr 2 > 0 if and only if ℎ is greater than critical depth
Hence,
dℎ
<0 (depth decreasing) if and only if ℎ lies between normal and critical depths.
d𝑥
Considering the whole range of possibilities allows a two-character classification of slopes (S1,
M3 etc) where:
• the first character is S, C, M, H, A for Steep, Critical, Mild, Horizontal, Adverse;
• the second character is 1, 2, 3 depending on where ℎ lies with respect to ℎ𝑐 and ℎ𝑛 .
Typical profiles and examples of where they occur are given in the table overleaf.
Note that:
• the normal depth (where 𝑆𝑓 = 𝑆0) is approached asymptotically (dℎ⁄d𝑥 → 0);
• the critical depth (where Fr = 1) is approached at right angles (dℎ⁄d𝑥 → ∞).
• above normal depth a backwater asymptotes a horizontal surface (M1, S1).
1
An alternative statement is that a downward slope is steep if it exceeds the critical slope (that slope at which the
normal flow is critical).
CRITICAL hc= hn C1
C1 ℎ > ℎ𝑐 = ℎ𝑛
(undesirable;
undular unsteady C3
flow) C3 ℎ𝑐 = ℎ𝑛 > ℎ
A control point is a location where there is a known relationship between water depth and
discharge (aka “stage-discharge relation”). Examples include critical-flow points (weirs,
venturi flumes, sudden changes in slope, free overfall), sluice gates, entry or discharge to a
reservoir. A hydraulic jump can also be classed as a control point. Control points often provide
a location where one can start a GVF calculation; i.e. a boundary condition.
(ii) Given a long-enough undisturbed fetch the flow will revert to normal flow.
(iii) A hydraulic jump occurs between regions of supercritical and subcritical gradually-
varied flow at the point where the jump condition for the sequent depths is correct.
(iv) Where the slope is mild (i.e. the normal flow is subcritical), and any downstream control
is a long way away, a hydraulic jump can be assumed to jump directly to the normal
depth.
normal M1
normal
hc
hn h1
h2 M3 hydraulic hn
WEIR
jump
CP CP
normal M1
normal
hn h1
h2 M3 hydraulic hn
jump
CP CP
S1
normal
hn
h1
h2 S3 normal
CP hn
CP
normal
hn
RESERVOIR
CP
hc S2
normal
RESERVOIR
normal M1 CP
hn
RESERVOIR
normal
M2
hn CP
hc
critical
Analytical solutions of the GVF equation are very rare and it is usual to solve it numerically.
The process yields a series of discrete pairs of distance 𝑥𝑖 and depth ℎ𝑖 . Intermediate points can
be determined, if required, by interpolation.
All methods employ a discrete approximation to one of the following forms of GVF equation:
d𝐻 𝑉2
= −𝑆𝑓 (total head 𝐻 = 𝑧𝑠 + 2𝑔) (11)
d𝑥
d𝐸 𝑉2
= 𝑆0 − 𝑆𝑓 (specific energy 𝐸 = ℎ + ) (12)
d𝑥 2𝑔
dℎ 𝑆0 − 𝑆𝑓
= (depth h) (13)
d𝑥 1 − Fr 2
In any of these the friction slope can be obtained by inverting the friction law; e.g. for
Manning’s equation:
𝑛2 𝑉 2 𝑛2 𝑄 2
𝑆𝑓 = 4/3
= 4/3 (14)
𝑅ℎ 𝑅ℎ 𝐴2
and the Froude number is
𝑉
Fr = (15)
√𝑔ℎ̅
where ℎ̅ = 𝐴/𝑏𝑠 is the mean depth (= actual depth for a rectangular or wide channel).
(upstream control)
CP
• backward in 𝑥 if the flow is subcritical
flow
(downstream control).
h4
There are two main classes of method: h2
h3
h0 h1
d𝐻
= −𝑆𝑓
d𝑥
This is solved as a standard-step method (find depth ℎ at specified distance intervals Δ𝑥). The
equation is discretised as
𝐻𝑖+1 − 𝐻𝑖 𝑆𝑓,𝑖 + 𝑆𝑓,𝑖+1
= −( ) (16)
Δ𝑥 2
solving sequentially for ℎ1 , ℎ2 , ℎ3 , … starting with the depth at the control point ℎ0 .
Since both 𝐻 and 𝑆𝑓 are functions of ℎ, the method operates by adjusting ℎ𝑖+1 iteratively at
each step so that the LHS and RHS of (16) are equal.
This is a good method, but since it requires iterative solution at each step it is better suited to a
computer program than hand or spreadsheet calculation.
d𝐸 𝑉2
= 𝑆0 − 𝑆𝑓 , where 𝐸 =ℎ+
d𝑥 2𝑔
The equation is then discretised (d𝑥/d𝐸 ≈ Δ𝑥/Δ𝐸 ) and rearranged for distance increments:
Δ𝐸
Δ𝑥 = (17)
(𝑆0 − 𝑆𝑓 )av
where
Δ𝐸 = 𝐸𝑖+1 − 𝐸𝑖 (18)
There are various ways of estimating the average slope difference. The example to follow uses
the average of values at depths ℎ𝑖 and ℎ𝑖+1.
dℎ 𝑆0 − 𝑆𝑓
=
d𝑥 1 − Fr 2
Here we shall solve this by a direct-step method (find displacement 𝑥 at specified depth
intervals Δℎ).
The function on the RHS is first written as a function of ℎ. The equation is then discretised
(d𝑥/dℎ ≈ Δ𝑥/Δℎ ) and rearranged for distance increments as:
d𝑥
Δ𝑥 = ( ) Δℎ (19)
dℎ av
As before, the bracketed term on the RHS can be taken as the average of values at the start and
end of an interval or (my preference) by evaluation at the interval mid-depth:
1
ℎmid = (ℎ + ℎ𝑖+1 )
2 𝑖
Comment. Different authors adopt different ways of solving the GVF equation numerically,
particularly in choosing whether to use the specific-energy or depth form, and how to form the
average derivative (e.g., average of values at the ends of the interval or simply the single value
at the midpoint). All should give the same answer when the step size Δℎ becomes very small,
but may differ for the larger step sizes typical of hand calculations. The specific-energy form
seems to be slightly more common in the literature, but my own tests suggest that the depth
form, with derivative evaluated just once at the mid-point of the interval, gives slightly better
results for large step sizes. Note that surface profiles become highly curved near critical points
and more steps, with a smaller Δℎ, should be used there.
(a) Determine the normal depth at this flow rate in the 5 m-wide channel.
(b) Show that critical conditions occur at the narrow section.
(c) Determine the depth just upstream of the narrowed section, where the width is 5 m.
(d) Determine the distance upstream to where the depth is 5% greater than the normal
depth, using two steps in the gradually-varied-flow equation.
(a)
Given:
𝑏 =5m (𝑏min = 1.2 m)
𝑆0 = 0.0002
𝑛 = 0.02 m−1⁄3 s
𝑄 = 10 m3 s −1
𝑏√𝑆 ℎ5/3
𝑄=
𝑛 (1 + 2ℎ/𝑏)2/3
Rearranging as an iterative formula for ℎ to find the normal depth at the channel slope 𝑆0 :
3/5
𝑛𝑄
ℎ=( ) (1 + 2ℎ/𝑏)2/5
𝑏√𝑆0
Here, with lengths in metres:
ℎ = 1.866(1 + 0.4ℎ)2/5
Iteration (from, e.g., ℎ = 1.866) gives normal depth:
ℎ𝑛 = 2.453 m
Answer: 2.45 m.
(b) To determine whether critical conditions occur, compare the total head in the approach flow
with that assuming critical conditions at the throat.
The total head, assuming normal flow and measuring heights from the bed of the channel is
𝑉𝑛2 𝑄2
𝐻𝑎 = 𝐸𝑎 = ℎ𝑛 + = ℎ𝑛 + = 2.487 m
2𝑔 2𝑔𝑏 2 ℎ𝑛2
At the throat the discharge per unit width is
Since the approach-flow head 𝐻𝑎 is less than the critical head 𝐻𝑐 (the minimum head required
to pass this flow rate through the venturi), the flow must back up and increase in depth just
upstream to supply this minimum head. It will then undergo a subcritical to supercritical
transition through the throat. The total head throughout the venturi is 𝐻 = 𝐻𝑐 = 2.880 m.
(c) In the vicinity of the venturi the total head is 𝐻 = 2.880 m. Just upstream (where width
𝑏 = 5 m), we seek the subcritical solution of
𝑉2 𝑄2
𝐻 = 𝑧𝑠 + =ℎ+
2𝑔 2𝑔𝑏 2 ℎ2
Rearrange for the deeper solution:
𝑄2
ℎ=𝐻−
2𝑔𝑏 2 ℎ2
Here, with lengths in metres:
0.2039
ℎ = 2.880 −
ℎ2
Iterate (from, e.g., ℎ = 2.880) to get the depth just upstream of the venturi:
ℎ = 2.855 m
Answer: 2.85 m
h0 = 2.855
(d) Do a GVF calculation (subcritical, so physically it
should start at the fixed downstream control and work h =12.716
h 2= 2.576
upstream, although mathematically it can be done the
other way) from the pre-venturi depth (ℎ = 2.855 m) Step 2 Step 1
to where ℎ = 2.576 m (i.e. 1.05 × ℎ𝑛 ). Using two
steps the depth increment per step is
x2 x1 x0 raised bed
2.576 − 2.855
Δℎ = = −0.1395 m
2
Both depth and specific-energy methods are shown on the following pages.
Δ𝑥 d𝑥
≈( )
Δℎ dℎ mid
dx
Δx ≈ ( ) Δh
dh mid
(‘mid’ means mid-point of interval, half way between ℎ𝑖 and ℎ𝑖+1 ; sometimes written ℎ𝑖+1⁄2 )
2
𝑉2 𝑄 2 /𝑏 2 0.4077
Fr = = =
𝑔ℎ 𝑔ℎ3 ℎ3
Manning’s equation (see earlier) gave
𝑏√𝑆 ℎ5/3
𝑄=
𝑛 (1 + 2ℎ/𝑏)2/3
Assuming that the rate of loss of energy (𝑆𝑓 ) at a general depth ℎ is the same as the channel
slope that would give normal flow at that depth, rearrangement for the slope gives
𝑛𝑄 2 (1 + 2ℎ/𝑏)4/3 −4
(1 + 0.4ℎ)4/3
𝑆𝑓 = ( ) = 16 × 10
𝑏 ℎ10/3 ℎ10/3
Hence,
0.4077
d𝑥 1 − Fr 2 1−
= = ℎ3
dℎ 𝑆0 − 𝑆𝑓 (1 + 0.4ℎ)4/3
(2 − 16 × ) × 10−4
ℎ10/3
With
d𝑥
Δ𝑥 = ( ) Δℎ and Δℎ = – 0.1395 m
dℎ mid
working may then be set out in tabular form. (All depths assumed to be in metres.)
d𝐸
= 𝑆0 − 𝑆𝑓
d𝑥
d𝑥 1
=
d𝐸 𝑆0 − 𝑆𝑓
Δ𝑥 1
≈
Δ𝐸 (𝑆0 − 𝑆𝑓 )av
Δ𝐸
Δ𝑥 ≈
(𝑆0 − 𝑆𝑓 )av
(‘av’ is taken as the average of values calculated at start and end of each interval in ℎ.)
Here:
𝑉2 𝑄 2 /𝑏 2 0.2039
𝐸 =ℎ+ =ℎ+ =ℎ+
2𝑔 2𝑔ℎ2 ℎ2
and the same expression as before may be used for 𝑆𝑓 , so that:
(1 + 0.4ℎ)4/3
𝑆0 − 𝑆𝑓 = 2 × 10−4 − 16 × 10−4
ℎ10/3
With
Δ𝐸
Δ𝑥 =
(𝑆0 − 𝑆𝑓 )av
working may then be set out in tabular form. (All depths assumed to be in metres.)
𝑖 ℎ𝑖 𝑥𝑖 𝐸𝑖 𝑆0 − 𝑆𝑓 Δ𝐸 (𝑆0 − 𝑆𝑓 )av Δ𝑥
0 2.855 0 2.880 6.61810–5
–0.136 5.67910–5 –2395
1 2.716 –2395 2.744 –5
4.74010
–0.137 3.59310–5 –3813
2 2.576 –6208 2.607 –5
2.44510
(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) Assuming no losses at the sluice, find the volume flow rate 𝑄.
(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.
(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.