SlideShare a Scribd company logo
Prof. Dr. Atıl BULU1
Chapter 6
Gradually-Varied Flow in Open Channels
6.1. Introduction
A steady non-uniform flow in a prismatic channel with gradual changes in its water-
surface elevation is named as gradually-varied flow (GVF). The backwater produced by a
dam or weir across a river and drawdown produced at a sudden drop in a channel are few
typical examples of GVF. In a GVF, the velocity varies along the channel and
consequently the bed slope, water surface slope, and energy line slope will all differ from
each other. Regions of high curvature are excluded in the analysis of this flow.
The two basic assumptions involved in the analysis of GVF are:
1. The pressure distribution at any section is assumed to be hydrostatic. This follows
from the definition of the flow to have a gradually varied water surface. As
gradual changes in the surface curvature give rise to negligible normal
accelerations, the departure from the hydrostatic pressure distribution is
negligible.
2. The resistance to flow at any depth is assumed to be given by the corresponding
uniform flow equation, such as the Manning equation, with the condition that the
slope term to be used in the equation is the energy line slope, Se and not the bed
slope, S0. Thus, if in a GVF the depth of flow at any section is y, the energy line
slope Se is given by,
34
22
R
Vn
Se = (6.1)
where R = hydraulic radius of the section at depth y.
6.2. Basic Differential Equation for the Gradually-Varied Flow Water Surface
Profile
Since,
S0 = Channel bed slope for uniform flow depth y0,
dy = Water depth variation for the dx canal reach,
d (V2
/2g) = Velocity head variation for the dx reach.
Writing the energy equation between the cross-sections 1 and 2 (Fig. 6.1.),
Prof. Dr. Atıl BULU2
( )
⎥
⎦
⎤
⎢
⎣
⎡
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+=−
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+=−
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+=−
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+++=++
g
V
dy
d
dx
dy
SS
g
V
dx
d
dx
dy
SS
g
V
ddydxSdxS
dxS
g
V
d
g
V
dyy
g
V
ydxS
e
e
e
e
2
1
2
2
222
2
0
2
0
2
0
222
0
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+
−
=
g
V
dy
d
SS
dx
dy e
2
1
2
0
(6.2)
Figure 6. 1.
( )
( )dyyTdA
ygA
Q
g
V
=→= 2
22
22
( ) ( )
( )
( )
( )
( )
2
2
3
2
4
2
2
2
2
2
2
rF
yT
yA
g
V
ygA
yTQ
ygA
dy
dA
AQ
ygA
Q
dy
d
−=−=
−=−=⎥
⎦
⎤
⎢
⎣
⎡
Substituting this to Equ. (6.2),
2
0
1 r
e
F
SS
dx
dy
−
−
= (6.3)
S0d
d
y
Sed
V
y+d
V+d
V2
/2g+d(V2
/2g
E.L
S0
Se
Prof. Dr. Atıl BULU3
Equ. (6.3) is the general differential equation of the water surface profile for the
gradually varied flows. dy/dx gives the variation of water depth along the channel in the
flow direction.
6.3. Classification of Flow Surface Profiles
For a given channel with a known Q = Discharge, n = Manning coefficient, and S0 =
Channel bed slope, yc = critical water depth and y0 = Uniform flow depth can be
computed. There are three possible relations between y0 and yc as 1) y0 > yc , 2) y0 < yc ,
3) y0 = yc .
For horizontal (S0 = 0), and adverse slope ( S0 < 0) channels,
21
0
321
SR
n
AQ =
Horizontal channel, S0 = 0→ Q = 0,
Adverse channel , S0 < 0 , Q cannot be computed,
For horizontal and adverse slope channels, uniform flow depth y0 does not exist.
Based on the information given above, the channels are classified into five categories as
indicated in Table (6.1).
Table 6.1. Classification of channels
Number Channel
category
Symbol Characteristic
condition
Remark
1 Mild slope M y0 > yc Subcritical flow at normal depth
2 Steep slope S yc > y0 Supercritical flow at normal
depth
3 Critical slope C yc = y0 Critical flow at normal depth
4 Horizontal
bed
H S0 = 0 Cannot sustain uniform flow
5 Adverse slope A S0 < 0 Cannot sustain uniform flow
For each of the five categories of channels, lines representing the critical depth (yc ) and
normal depth (y0 ) (if it exists) can be drawn in the longitudinal section. These would
divide the whole flow space into three regions as:
Region 1: Space above the topmost line,
Region 2: Space between top line and the next lower line,
Region 3: Space between the second line and the bed.
Prof. Dr. Atıl BULU4
Figure (6.2) shows these regions in the various categories of channels.
Figure 6.2. Regions of flow profiles
Depending upon the channel category and region of flow, the water surface profiles will
have characteristics shapes. Whether a given GVF profile will have an increasing or
decreasing water depth in the direction of flow will depend upon the term dy/dx in Equ.
(6.3) being positive or negative.
2
0
1 r
e
F
SS
dx
dy
−
−
= (6.3)
Prof. Dr. Atıl BULU5
For a given Q, n, and S0 at a channel,
y0 = Uniform flow depth,
yc = Critical flow depth,
y = Non-uniform flow depth.
The depth y is measured vertically from the channel bottom, the slope of the water
surface dy/dx is relative to this channel bottom. Fig. (6.3) is basic to the prediction of
surface profiles from analysis of Equ. (6.3).
Figure 6.3
To assist in the determination of flow profiles in various regions, the behavior of dy/dx at
certain key depths is noted by studying Equ. (6.3) as follows:
1
1
1
>→<
=→=
<→>
rc
rc
rc
Fyy
Fyy
Fyy
And also,
00
00
00
SSyy
SSyy
SSyy
e
e
e
>→<
=→=
<→>
1. As 0yy → , 0VV → , 0SSe =
0
0
1
lim 2
00
0
==
−
−
=
→ consF
SS
dx
dy
r
yy
The water surface approaches the normal depth asymptotically.
Prof. Dr. Atıl BULU6
2. As cyy → , 12
=rF , 01 2
=− rF ,
∞=
−
=
−
−
=
→ 01
lim 0
2
0 e
r
e
yy
SS
F
SS
dx
dy
c
The water surface meets the critical depth line vertically.
3. As ∞→y , 000 →→=→= er SFV
0
0
2
0
11
lim S
S
F
SS
dx
dy
r
e
y
==
−
−
=
∞→
The water surface meets a very large depth as a horizontal asymptote.
Based on this information, the various possible gradually varied flow profiles are grouped
into twelve types (Table 6.2).
Table 6.2. Gradually Varied Flow profiles
Channel Region Condition Type
Mild slope
1
2
3
y > y0 > yc
y0 > y > yc
y0 > yc > y
M1
M2
M3
Steep slope
1
2
3
y> yc > y0
yc > y > y0
yc > y0 > y
S1
S2
S3
Critical slope
1
3
y > y0 = y2
y < y0 = yc
C1
C3
Horizontal bed 2
3
y > yc
y < yc
H2
H3
Adverse slope 2
3
y > yc
y < yc
A2
A3
Prof. Dr. Atıl BULU7
6.4. Water Surface Profiles
6.4.1. M – Curves
Figure 6.4
General shapes of M curves are given in Fig. (6.4). Asymptotic behaviors of each curve
will be examined mathematically.
a) M1 – Curve
Water surface will be in Region 1 for a mild slope channel and the flow is obviously
subcritical.
Se < S0 → Mild slope channel
y0 > yc → Subcritical flow
2
0
1 r
e
F
SS
dx
dy
−
−
=
Fr < 1 → Subcritical flow → (1 – Fr
2
) > 0
y > y0 → Se < S0
0>
+
+
=
dx
dy
(Water depth will increase in the flow direction)
Asymptotic behavior of the water surface is;
Water depth can be between (∞ > y > y0) for Region 1. The asymptotic behaviors of the
water surface for the limit values (∞ , y0) are;
Prof. Dr. Atıl BULU8
a) y → ∞ , V → 0, Fr → 0, (1 – Fr
2
) = 1
y → ∞ , V → 0, Se →0
0
0
1
0
lim S
S
dx
dy
y
=
−
=
∞→
The water surface meets a very large depth as a horizontal asymptote.
b) y → y0 , V → V0 , Se → S0
0
1
lim 2
00
0
=
−
−
=
→
r
yy F
SS
dx
dy
The water surface approaches the normal depth asymptotically.
The most common of all GVF profiles is the M1 type, which is a subcritical flow
condition. Obstructions to flow, such as weirs, dams, control structures and natural
features, such as bends, produce M1 backwater curves (Fig. 6.5). These extend to several
kilometers upstream before merging with the normal depth.
Figure 6.5. M1 Profile
b) M2 – Curve
Water surface will be in Region 2 for a mild slope channel and the flow is obviously
subcritical. (Fig. 6.4).
y0 > y > yc
y < y0 → V > V0 → Se > S0 → (S0 –Se) < 0
y > yc → Fr < 1 → (1 – Fr
2
) > 0
−=
−
+
=
−
−
= 2
0
1 r
e
F
SS
dx
dy
(Water depth decrease in the flow direction)
Prof. Dr. Atıl BULU9
Asymptotic behavior of the water surface is;
y0 > y > yc
y → y0, Se → S0 , (S0 –Se) = 0
0
1
lim 2
00
0
=
−
−
=
→
r
yy F
SS
dx
dy
The water surface approaches the normal depth asymptotically.
y → yc , Fr = 1 , (1 – Fr
2
) = 0
∞=
−
=
→ 0
lim 0 e
yy
SS
dx
dy
c
The water surface meets the critical depth line vertically.
The M2 profiles occur at a sudden drop of the channel, at constriction type of transitions
and at the canal outlet into pools (Fig. 6.6).
Figure 6.6
c) M3 – Curve
Water surface will be in Region 3 for a mild slope channel and the flow is obviously
subcritical. (Fig. 6.4).
y0 > yc >0
y < y0 → V > V0 → Se > S0 → (S0 – Se) < 0
y < yc → Supercritical flow → Fr > 1 → (1 – Fr
2
) < 0
+=
−
−
=
−
−
= 2
0
1 r
e
F
SS
dx
dy
(Water depth will increase in the flow direction)
Asymptotic behavior of the water surface is;
Prof. Dr. Atıl BULU10
y → yc → Fr = 1 → (1 – Fr
2
) = 0
∞=
−
=
→ 0
lim 0 e
yy
SS
dx
dy
c
The water surface meets the critical depth line vertically.
y → 0, Se → ∞, (S0 – Se) = ∞
y → 0, ∞==
gy
V
Fr
∞
∞
=
→ dx
dy
y 0
lim (Unknown)
The angle of the water surface with the channel bed may be taken as
3
0
0 ⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
cy
y
S .
Figure 6.7. M3 Profile
Where a supercritical stream enters a mild slope channel, M3 type of profile occurs. The
flow leading from a spillway or a sluice gate to a mild slope forms a typical example
(Fig. 6.7). The beginning of the M3 curve is usually followed by a small stretch of rapidly
varied flow and the downstream is generally terminated by a hydraulic jump. Compared
to M1 and M2 profiles, M3 curves are of relatively short length.
Prof. Dr. Atıl BULU11
Example 6.1 : A rectangular channel with a bottom width of 4.0 m and a bottom slope of
0.0008 has a discharge of 1.50 m3
/sec. In a gradually varied flow in this channel, the
depth at a certain location is found to be 0.30 m. assuming n = 0.016, determine the type
of GVF profile.
Solution:
a) To find the normal depth y0,
0
0
00
00
24
4
242
4
y
y
P
A
R
yyBP
yByA
+
==
+=+=
==
( )
( )
5.0
32
0
35
0
21
0
32
0008.0
24
4
016.0
1
50.1
1
×
+
×=
⎟
⎠
⎞
⎜
⎝
⎛
==
y
y
S
P
A
n
AVAQ
( ) 32
0
35
0
24
0842.0
y
y
+
=
By trial and error, y0 = 0.43 m.
b) Critical depth yc,
m
g
q
y
m
B
Q
q
c 24.0
81.9
375.0
sec375.0
0.4
50.1
312312
2
=⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
=⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
=
===
c) Type of profile,
y0 = 0.43 m > yc = 0.24 m (Mild slope channel, M profile)
y = 0.30 m
y0 > y > yc (Region 2)
Water surface profile is of the M2 type.
Prof. Dr. Atıl BULU12
6.4.2. S – Curves
Figure 6.8
General shapes of S curves are given in Fig. (6.8). Asymptotic behaviors of each curve
will be examined mathematically.
a) S1 – Curve
Water surface will be in Region 1 for a steep slope channel and the flow is obviously
supercritical.
y > yc > y0
y > y0 → V < V0 → Se < S0 → (S0 – Se) > 0
y > yc → Fr < 1 → (1 – Fr
2
) > 0
+=
+
+
=
−
−
= 2
0
1 r
e
F
SS
dx
dy
(Water depth will increase in the flow direction)
Asymptotic behavior of the water surface is;
∞ > y > yc
y → ∞ , V → 0 , Fr = 0 , (1- Fr
2
) = 1
y → ∞ , V → 0 , Se = 0
0
0
2
0
01
0
1
lim S
S
F
SS
dx
dy
r
e
y
=
−
−
=
−
−
=
∞→
Prof. Dr. Atıl BULU13
The water surface meets a very large depth as a horizontal asymptote.
y → yc , Fr = 1 , (1 – Fr
2
) = 0
∞=
−
=
−
−
=
→ 01
lim 0
2
0 e
r
e
yy
SS
F
SS
dx
dy
c
The water surface meets the critical depth line vertically.
Figure 6.9. S1 Profile
The S1 profile is produced when the flow from a steep channel is terminated by a deep
pool created by an obstruction, such as a weir or dam (Fig. 6.9). At the beginning of the
curve, the flow changes from the normal depth (supercritical flow) to subcritical flow
through a hydraulic jump. The profiles extend downstream with a positive water slope to
reach a horizontal asymptote at the pool elevation.
b) S2 – Curve
Water surface will be in Region 2 for a steep slope channel and the flow is supercritical.
y0 < y < yc
y > y0 → Se < S0 → (S0 – Se) > 0
y < yc → Supercritical flow → Fr > 1 → (1 – Fr
2
) < 0
−=
−
+
=
−
−
= 2
0
1 r
e
F
SS
dx
dy
(Water depth will decrease in the flow direction)
Prof. Dr. Atıl BULU14
Asymptotic behavior of the water surface is;
y0 < y < yc
y → yc , Fr → 1 , (1 –Fr
2
) = 0
∞=
−
=
−
−
=
→ 01
lim 0
2
0 e
r
e
yy
SS
F
SS
dx
dy
c
The water surface meets the critical depth line vertically.
y → y0 , V → V0 , Se → S0 , (S0 – Se) = 0
0
1
0
lim 2
0
=
−
=
→
r
yy Fdx
dy
The water surface approaches the normal depth asymptotically.
Figure 6.10. S2 Profile
Profiles of the S2 type occur at the entrance region of a steep channel leading from a
reservoir and at a brake of grade from mild slopes to steep slope (Fig. 6.10). Generally S2
profiles are short of length.
c) S3 – Curve
Water surface will be in Region 3 for a steep slope channel and the flow is supercritical.
0 < y < y0
y < y0 → Se > S0 → (S0 – Se) < 0
y < yc → Supercritical flow → Fr > 1 → (1 – Fr
2
) < 0
+=
−
−
=
−
−
= 2
0
1 r
e
F
SS
dx
dy
(Water depth will increase in the flow direction)
Prof. Dr. Atıl BULU15
Asymptotic behavior of the water surface is;
0 < y < y0
y → 0, Se → ∞ , (S0 – Se) = -∞
y → 0 , ∞==
gy
V
Fr
∞
∞−
=
→ dx
dy
y 0
lim (Unknown)
The angle of the water surface with the channel bed may be taken as
3
0
0 ⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
cy
y
S .
y → y0 , Se → S0 , (S0 – Se) = 0
0
1
0
lim 2
0
=
−
=
→
r
yy Fdx
dy
The water surface approaches the normal depth y0 asymptotically.
Figure 6.11. S3 Profile
Figure 6.12. S3 Profile
Free flow from a sluice gate with a steep slope on its downstream is of the S3 type (Fig.
6.11). The S3 curve also results when a flow exists from a steeper slope to a less steep
slope (Fig. 6.12).
Prof. Dr. Atıl BULU16
6.4.3. C – Curves
Figure 6.13
General shapes of C curves are given in Fig. (6.13). Asymptotic behaviors of each curve
will be examined mathematically. Since the flow is at critical stage, y0 = yc, there is no
Region 2.
a) C1 – Curve
Water surface will be in Region 1 for a critical slope channel.
y0 = yc < y < ∞
y > yc → Se < S0 → (S0 –Se) > 0
y > yc → Subcritical flow → Fr < 1 → (1 –Fr
2
) > 0
+=
+
+
=
−
−
= 2
0
1 r
e
F
SS
dx
dy
(Water depth will increase in the flow direction)
Asymptotic behavior of the water surface is;
y → ∞ , Se → 0 , ( S0 – Se) = S0
y → ∞, V = 0 → Fr = 0 → (1 – Fr
2
) = 1
0
0
01
0
lim S
S
dx
dy
y
=
−
−
=
∞→
The water surface meets a very large depth as a horizontal asymptote.
y → yc , Fr = 1 → (1 – Fr
2
) = 0
y → yc = y0 , Se = S0 = Sc → (S0 –Se ) =0
0
0
lim
0
=
=→ dx
dy
yyy c
(Unknown)
C1 and C3 profiles are very rare and highly unstable
Prof. Dr. Atıl BULU17
6.4.4. H – Curves
Figure 6.14
General shapes of H curves are given in Fig. (6.14). For horizontal slope channels,
uniform flow depth y0 does not exist. Critical water depth can be computed for a given
discharge Q and therefore critical water depth line can be drawn. Since there is no
uniform water depth y0, Region 1 does not exist.
a) H2 – Curve
Water surface will be in Region 2 for a horizontal slope channel.
∞ > y > yc
y > yc → Se < S0 = 0 → (S0 –Se) < 0
y > yc → subcritical flow → Fr < 1 → (1 – Fr
2
) > 0
−=
+
−
=
dx
dy
(Water depth will decrease in the flow direction)
Asymptotic behavior of the water surface is;
y → ∞ , Se → 0 , ( S0 – Se) = S0
y → ∞, V = 0 → Fr = 0 → (1 – Fr
2
) = 1
0
0
01
0
lim S
S
dx
dy
y
=
−
−
=
∞→
The water surface meets a very large depth as a horizontal asymptote.
y → yc , Fr = 1 → (1 – Fr
2
) = 0
y → yc , Se = Sc → S0 = 0, (S0 –Se ) = -Se
∞=
−
=
−
−
=
→ 0
0
1
lim 2
0 e
r
e
yy
S
F
SS
dx
dy
c
The water surface meets the critical depth line vertically.
Prof. Dr. Atıl BULU18
b) H3 – Curve
0 > y > yc
y < yc → Se > S0 = 0 → (S0 –Se) = -Se
y → yc , supercritical flow , Fr > 1 → (1 – Fr
2
) < 0
+=
−
−
=
−
−
= 2
0
1 r
e
F
SS
dx
dy
(Water depth will increase in the flow direction)
Asymptotic behavior of the water surface is;
y → 0, Se → ∞ , (S0 – Se) = -∞
y → 0 , ∞==
gy
V
Fr
∞
∞−
=
→ dx
dy
y 0
lim (Unknown)
y → yc → Fr = 1 → (1 – Fr
2
) = 0
y → yc , Se = Sc → S0 = 0, (S0 –Se ) = -Se
∞=
−
=
−
−
=
→ 0
0
1
lim 2
0 e
r
e
yy
S
F
SS
dx
dy
c
The water surface meets the critical depth line vertically.
Figure 6.15
A horizontal channel can be considered as the lower limit reached by a mild slope as its
bed slope becomes flatter. The H2 and H3 profiles are similar to M2 and M3 profiles
respectively (Fig. 6.15). However, the H2 curve has a horizontal asymptote.
Prof. Dr. Atıl BULU19
6.4.5. A– Curves
Figure 6.16
General shapes of A curves are given in Fig. (6.16). For adverse slope channels, uniform
flow depth y0 does not exist. Critical water depth can be computed for a given discharge
Q and therefore critical water depth line can be drawn. Since there is no uniform water
depth y0, Region 1 does not exist as well as in A curves. A2 and A3 curves are similar to
H2 and H3 curves respectively.
Figure 6.17
Adverse slopes are rather rare (Fig. 6.17). These profiles are of very short length.
6.5. Control Sections
A control section is defined as a section in which a fixed relationship exists between the
discharge and depth of flow. Weirs, spillways, sluice gates are some typical examples of
structures which give rise to control sections. The critical depth is also a control point.
However, it is effective in a flow profile which changes from subcritical to supercritical
flow. In the reverse case of transition from supercritical flow to subcritical flow, a
hydraulic jump is usually formed bypassing the critical depth as a control point. Any
GVF profile will have at least one control section.
Prof. Dr. Atıl BULU20
(d)
Figure 6.18
In the synthesis of GVF profiles occurring in a serially connected channel elements, the
control sections provide a key to the identification of proper profile shapes. A few typical
control sections are shown in Fig. (6.18, a-e). It may be noted that subcritical flows have
controls in the downstream end, while supercritical flows are governed by control
sections existing at the upstream end of the channel section. In Figs. (6.18 a and b) for the
M1 profile, the control section (indicated by a dark dot in the figures) is just at the
upstream of the spillway and sluice gate respectively. In Figs. (6.18 b and d) for M3 and
S3 profiles respectively, the control point is at the vena contracta of the sluice gate flow.
In subcritical flow reservoir offtakes (Fig. 6.18 c), even though the discharge is governed
Prof. Dr. Atıl BULU21
by the reservoir elevation, the channel entry section is not strictly a control section. The
water surface elevation in the channel, will be lower than the reservoir elevation by a
headloss amount equivalent to (1 + ζ )V2
/2g where ζ is the entrance loss coefficient. The
true control section will be at a downstream location in the channel. For the situation
shown in Fig. (6.18 c) the critical depth at the free overflow at the channel end acts as the
downstream control. For a sudden drop (free overflow) due to curvature of the
streamlines the critical depth usually occurs at distance of about 4yc upstream of the drop.
This distance, being small compared to GVF lengths, is neglected and it is usual to
perform calculations by assuming yc to occur at the drop.
(e)
For a supercritical canal intake (Fig. 6.18 e), the reservoir water surface falls to the
critical depth at the head of the canal and then onwards the water surface follows the S2
curve. The critical depth occurring at the upstream end of the channel is the control for
this flow.
Figure 6.19
Prof. Dr. Atıl BULU22
6.5.1. Feeding a Pool for Subcritical Flow
A mild slope channel discharging into a pool of variable surface elevation is indicated in
Fig. (6.19). Four cases are shown.
1) The pool elevation is higher than the elevation of the normal (uniform) water
depth line at B. This gives rise to a drowning of the channel end. A profile of the
M1 type is produced with the pool level at B as control.
2) The pool elevation is lower than the elevation of the normal depth line but higher
than the critical depth line at B. The pool elevation acts as a control for the M2
curve.
3) The pool elevation has dropped down to that of the critical depth line at B and the
control is still at the pool elevation.
4) The pool elevation has dropped lower than the elevation of the critical depth line
at B. The water surface cannot pass through a critical depth at any location other
than B and hence a sudden drop in the water surface at B is observed. The critical
depth at B is the control of this flow.
6.5.2. Feeding a Pool for Supercritical Flow
When the slope of the channel is steep, that is greater than the critical slope, the flow in
the channel becomes supercritical. In practical applications, steep channels are usually
short, such as raft and log chutes that are used as spillways.
Figure 6.20
As the control section in a channel of supercritical flow is at the upstream end, the
delivery of the channel is fully governed by the critical discharge at section 1.
Prof. Dr. Atıl BULU23
1) When the tailwater level B is less than the outlet depth (normal depth) at section
2, the flow in the canal is unaffected by the tailwater. The flow profile passes
through the critical water depth near C and approaches the normal depth by means
of a smooth drawdown curve of the S2 type.
2) When the tailwater level B is greater than the outlet depth, the tailwater will raise
the water level in the upstream portion of the canal to form an S1 profile between j
and b`, producing a hydraulic jump at the end j of the profile. However, the flow
upstream from the jump will not be affected by the tailwater.
Such a hydraulic jump is objectionable and dangerous, particularly when the canal
is raft chute or some other structure intended to transport a floating raft from the
upstream reservoir to a downstream pool. A neutralizing reach may be suggested
for the solution of this problem. (Chow, 1959). (Fig. 6.21)
Figure 6.21. Elimination of hydraulic jump by neutralizing reach.
In this reach the bottom slope of the channel is made equal to the critical slope.
According to a corresponding case of C1 profile in Fig. (6.21), the tailwater levels
will be approximately horizontal lines which intersect the surface of flow in the
canal without causing any disturbance. At the point of intersection, theoretically,
there is a jump of zero height.
3) As the tailwater rises further, the jump will move upstream, maintaining its height
and form in the uniform flow zone nb, until it reaches point n. The height of the
jump becomes zero when it reaches the critical depth at c. In the meantime, the
flow profile reaches its theoretical limit cb`` of the S1 profile. Beyond this limit
the incoming flow will be directly affected by the tailwater. In practical
applications, the horizontal line cb``` may be taken as the practical limit of the
tailwater stage.
Prof. Dr. Atıl BULU24
6.6. Analysis of Flow Profile
A channel carrying a gradually varied flow can in general contain different prismoidal
channel cross-sections of varying hydraulic properties. There can be a number of control
sections at various locations. To determine the resulting water surface profile in a given
case, one should be in a position to analyze the effects of various channel sections and
controls connected in series. Simple cases are illustrated to provide information and
experience to handle more complex cases.
6.6.1. Break in Grades
Simple situations of a series combination of two channel sections with differing bed
slopes are considered. In Fig. (6.22.a), a break in grade from a mild channel to a milder
channel is shown. It is necessary to first draw the critical-depth line (CDL) and the
normal-depth line (NDL) for both slopes. Since yc does not depend upon the slope for a
taken Q = discharge, the CDL is at a constant height above the channel bed in both
slopes. The normal depth y01 for the mild slope is lower than that of the of the milder
slope (y02). In this case, y02 acts as a control, similar to the weir or spillway case and an
M1 backwater curve is produced in the mild slope channel.
Various combinations of slopes and the resulting GVF profiles are presented in Fig. (
6.22, a-h). It may be noted that in some situations there can be more than one possible
profiles. For example, in Fig. (6.22 e), a jump and S1 profile or an M3 profile and a jump
possible. The particular curve in this case depends on the channel and its flow properties.
In the examples indicated in Fig. (6.22), the section where the grade changes acts a
control section and this can be classified as a natural control. It should be noted that even
though the bed slope is considered as the only variable in the above examples, the same
type of analysis would hold good for channel sections in which there is a marked change
in the roughness characteristics with or without change in the bed slope. A long reach of
unlined canal followed by a line reach serves as a typical example for the same. A change
in the channel geometry (the bed width or side slope) beyond a section while retaining
the prismoidal nature in each reach also leads to a natural control section.
Prof. Dr. Atıl BULU25
Prof. Dr. Atıl BULU26
Prof. Dr. Atıl BULU27
Figure 6.22
6.6.2. Serial Combination of Channel Sections
To analyze a general problem of many channel sections and controls, the following steps
are to be applied.
1. Draw the longitudinal section of the system.
2. Calculate the critical depth and normal depths of various reaches and draw the
CDL and NDL in all reaches.
3. Mark all the controls, both the imposed as well as natural controls.
4. Identify the possible profiles.
Example 6. 2: Identify and sketch the GVF profiles in three mild slopes which could be
described as mild, steeper mild and milder. The three slopes are in series. The last slope
has a sluice gate in the middle of the reach and the downstream end of the channel has a
free overfall.
Solution: The longitudinal profile of the channel, critical depth line and normal depth
lines for the various reaches are shown in the Fig. (6.23). The free overfall at E is
obviously a control. The vena contracta downstream of the sluice gate at D is another
control. Since for subcritical flow the control is at the downstream end of the channel, the
higher of the two normal depths at C acts as a control for the reach CB, giving rise to an
M1 profile over CB. At B, the normal depth of the channel CB acts as a control giving
rise to an M2 profile over AB. The controls are marked distinctly in Fig. (6.23). With
these controls the possible flow profiles are: an M2 profile on channel AB, M1 profile on
channel BC, M3 profile and M2 profile through a jump on the stretch DE.
Prof. Dr. Atıl BULU28
Figure 6.23
Example 6.3: A trapezoidal channel has three reaches A, B, and C connected in series
with the following physical characteristics.
Reach Bed width
(B)
Side slope
(m)
Bed slope
S0
n
1
2
3
4.0 m
4.0 m
4.0 m
1.0
1.0
1.0
0.0004
0.009
0.004
0.015
0.012
0.015
For a discharge Q = 22.5 m3
/sec through this channel, sketch the resulting water surface
profiles. The length of the reaches can be assumed to be sufficiently long for the GVF
profiles to develop fully.
Solution: The normal depths and critical water depths in the various reaches are
calculated as:
( )
myBT
myBP
ymyBA
2
12 2
+=
++=
+=
B=4m
y
m
1
T
Prof. Dr. Atıl BULU29
yT
yyP
yyA
24
241124
)4(
232
+=
+=++=
+=
( )[ ]
( ) 32
0
23
5.0
0
35
00
5.0
0
32
24
4
5.22
y
S
n
yy
S
P
A
n
A
AVQ
×+
×
+
=
⎟
⎠
⎞
⎜
⎝
⎛
==
Uniform flow depths for the given data for every reach are calculated by trial and error
method;
Reach A: S0A = 0.0004 , nA=0.015 → y0A = 2.26 m
Reach B: S0B = 0.009 , nB =0.012 → y0B = 0.81 m
Reach C: S0C = 0.004 , nC = 0.015 → y0C =1.17 m
Since the channel is prismoidal (the geometry does not change in the reaches), there will
be only one critical water depth.
( )
( )[ ]
1
481.9
245.22
1
3
2
3
2
=
+×
+×
=
cc
c
c
c
yy
y
gA
TQ
yc =1.32 m
Reach A is a mild slope channel as y0A = 2.26 m > y0c = 1.32 m and the flow is
subcritical. Reach B and C are steep slope channels as y0B, y0C < yc and the flow is
supercritical on both reaches. Reach B is steeper than reach C. The various reaches are
schematically shown in Fig. (6.24). The CDL is drawn at a height of 1.32 m above the
bed level and NDLs are drawn at the calculated y0 values.
The controls are marked in the figure. Reach A will have an M2 drawdown curve, reach B
an S2 drawdown curve and reach C a rising curve as shown in the figure. It may be noted
that the resulting profile as above is a serial combination of Fig. (6.22 d and f).
Prof. Dr. Atıl BULU30
Figure 6.23

More Related Content

PPTX
Design of gravity dams
PPTX
Energy and momentum principles in open channel flow
PDF
Lecture notes 05
PDF
Chapter 7 gvf
PPTX
guide banks.pptx
PPTX
PDF
RESERVOIR SEDIMENTATION
PDF
design of weirs
Design of gravity dams
Energy and momentum principles in open channel flow
Lecture notes 05
Chapter 7 gvf
guide banks.pptx
RESERVOIR SEDIMENTATION
design of weirs

What's hot (20)

PPTX
Gradually varied flow
PPTX
Earthen dam
PPT
Ce533 ch1
PPT
ogee design
PPTX
Storage and distribution reservoirs
PPTX
Diversion Headworks
PDF
Chapter 6 design of irrigation channel1
PDF
Unit 12 River training work.pdf
PPTX
PPTX
Design of Hydraulic Structures
PPTX
Critical flow through an Open channel
PDF
Chapter 7 spillway and energy dissipators
PPTX
Rapidly varied flow
PPTX
Open Channel Flow
PPTX
Hydraulic structures
PPTX
S curve hydrograph
PPTX
Gravity Dam (numerical problem ) BY SITARAM SAINI
PPTX
Gradually varied flow
Earthen dam
Ce533 ch1
ogee design
Storage and distribution reservoirs
Diversion Headworks
Chapter 6 design of irrigation channel1
Unit 12 River training work.pdf
Design of Hydraulic Structures
Critical flow through an Open channel
Chapter 7 spillway and energy dissipators
Rapidly varied flow
Open Channel Flow
Hydraulic structures
S curve hydrograph
Gravity Dam (numerical problem ) BY SITARAM SAINI
Ad

Viewers also liked (20)

PPT
open channel flow
PPT
Chap1 open channel flow
PDF
(Part ii)- open channels
PDF
(Part iii)- open channels
PDF
Open Channel Flows (Lecture notes 04)
PPT
Chap1 open channel flow
PDF
Open channel design
PDF
Lecture notes 06
DOCX
Introduction chapter `1
PPTX
Water surface profiles
PDF
Open channels
PDF
Answers assignment 5 open channel hydraulics-fluid mechanics
PDF
(Part i)- open channels
PDF
05 gradually varied flow slide
PDF
KINEMATICS OF FLUIDS (Chapter 3)
ODP
Aram english
PPTX
Pengurusan khidmat pelanggan
PPT
Caribbean Marketing Director's Primary Purposes Rich Sadler
PDF
TRANSCRIPT_
PPTX
Anahit
open channel flow
Chap1 open channel flow
(Part ii)- open channels
(Part iii)- open channels
Open Channel Flows (Lecture notes 04)
Chap1 open channel flow
Open channel design
Lecture notes 06
Introduction chapter `1
Water surface profiles
Open channels
Answers assignment 5 open channel hydraulics-fluid mechanics
(Part i)- open channels
05 gradually varied flow slide
KINEMATICS OF FLUIDS (Chapter 3)
Aram english
Pengurusan khidmat pelanggan
Caribbean Marketing Director's Primary Purposes Rich Sadler
TRANSCRIPT_
Anahit
Ad

Similar to Gradually-Varied Flow in Open Channels ( Lecture notes 06) (20)

PDF
Gradually Varied Profile & Rapidly Varied Profile - Hydraulics And Hydraulics...
PPTX
Chapter 7:Non Uniform G.V.F
PPTX
Module 2- Non Uniform Flow in Open Channels (HHM).pptx
PPTX
Module-III Varried Flow.pptx GVF Definition, Water Surface Profile Dynamic Eq...
PDF
Chapter 2 open channel hydraulics
PPTX
chapter GVF.pptx
PPTX
Open channel Flow
PDF
Gradually varied flow and rapidly varied flow
PDF
TWO-DIMENSIONAL FLOW OF THE REAL FLUIDS (Lecture notes 07)
PDF
4.gradually varid flow part 2.pdf sjjbndnnnnjj
DOC
Chapter 1
PDF
Flood risk evaluation for Serio river, Italy
PDF
River hydraulic modelling for river Serio (Northern Italy), 2014
PDF
Local Energy (Head) Losses (Lecture notes 03)
PDF
Flows under Pressure in Pipes (Lecture notes 02)
PDF
3. UJRRA_22_04.pdf
PDF
Chapter 6 energy and momentum principles
DOCX
Paper_Horsley
PDF
Non- Newtonian behavior of blood in very narrow vessels
PDF
I027055062
Gradually Varied Profile & Rapidly Varied Profile - Hydraulics And Hydraulics...
Chapter 7:Non Uniform G.V.F
Module 2- Non Uniform Flow in Open Channels (HHM).pptx
Module-III Varried Flow.pptx GVF Definition, Water Surface Profile Dynamic Eq...
Chapter 2 open channel hydraulics
chapter GVF.pptx
Open channel Flow
Gradually varied flow and rapidly varied flow
TWO-DIMENSIONAL FLOW OF THE REAL FLUIDS (Lecture notes 07)
4.gradually varid flow part 2.pdf sjjbndnnnnjj
Chapter 1
Flood risk evaluation for Serio river, Italy
River hydraulic modelling for river Serio (Northern Italy), 2014
Local Energy (Head) Losses (Lecture notes 03)
Flows under Pressure in Pipes (Lecture notes 02)
3. UJRRA_22_04.pdf
Chapter 6 energy and momentum principles
Paper_Horsley
Non- Newtonian behavior of blood in very narrow vessels
I027055062

More from Shekh Muhsen Uddin Ahmed (20)

PPTX
Building information modeling
PDF
Building Information Modeling
PDF
Earthquake Load Calculation (base shear method)
PDF
Structural engineering iii
PDF
Structural engineering ii
PDF
Structural engineering i
PDF
Theory of elasticity and plasticity (Equations sheet, Part 3) Att 9035
PDF
Theory of elasticity and plasticity (Equations sheet part 01) Att 8676
PDF
2-D formulation Plane theory of elasticity Att 6672
PDF
Energy principles Att 6607
PDF
Elasticity problem formulation Att 6582
PDF
Introduction to Theory of elasticity and plasticity Att 6521
PDF
TWO-DIMENSIONAL IDEAL FLOW (Chapter 6)
PDF
BASIC EQUATIONS FOR ONE-DIMENSIONAL FLOW (Chapter 04)
PDF
FUNDAMENTALS of Fluid Mechanics (chapter 01)
PDF
DIMENSIONAL ANALYSIS (Lecture notes 08)
PDF
Hydroelectric power plants (chapter 5)
PDF
Dimensional Analysis - Model Theory (Lecture notes 01)
PDF
Specific Energy (Lecture notes 05)
PPTX
Reinforced slab bridge design(AASHTO allowable stress design method)
Building information modeling
Building Information Modeling
Earthquake Load Calculation (base shear method)
Structural engineering iii
Structural engineering ii
Structural engineering i
Theory of elasticity and plasticity (Equations sheet, Part 3) Att 9035
Theory of elasticity and plasticity (Equations sheet part 01) Att 8676
2-D formulation Plane theory of elasticity Att 6672
Energy principles Att 6607
Elasticity problem formulation Att 6582
Introduction to Theory of elasticity and plasticity Att 6521
TWO-DIMENSIONAL IDEAL FLOW (Chapter 6)
BASIC EQUATIONS FOR ONE-DIMENSIONAL FLOW (Chapter 04)
FUNDAMENTALS of Fluid Mechanics (chapter 01)
DIMENSIONAL ANALYSIS (Lecture notes 08)
Hydroelectric power plants (chapter 5)
Dimensional Analysis - Model Theory (Lecture notes 01)
Specific Energy (Lecture notes 05)
Reinforced slab bridge design(AASHTO allowable stress design method)

Recently uploaded (20)

PPTX
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
PPT
Project quality management in manufacturing
PPTX
Internet of Things (IOT) - A guide to understanding
PDF
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
PDF
July 2025 - Top 10 Read Articles in International Journal of Software Enginee...
PPT
Mechanical Engineering MATERIALS Selection
PDF
R24 SURVEYING LAB MANUAL for civil enggi
PPTX
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
DOCX
ASol_English-Language-Literature-Set-1-27-02-2023-converted.docx
PDF
composite construction of structures.pdf
PPTX
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
PDF
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
PDF
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
PPTX
Construction Project Organization Group 2.pptx
PDF
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
PPT
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
PPTX
UNIT-1 - COAL BASED THERMAL POWER PLANTS
PDF
Embodied AI: Ushering in the Next Era of Intelligent Systems
PDF
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
PDF
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
Project quality management in manufacturing
Internet of Things (IOT) - A guide to understanding
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
July 2025 - Top 10 Read Articles in International Journal of Software Enginee...
Mechanical Engineering MATERIALS Selection
R24 SURVEYING LAB MANUAL for civil enggi
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
ASol_English-Language-Literature-Set-1-27-02-2023-converted.docx
composite construction of structures.pdf
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
Construction Project Organization Group 2.pptx
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
UNIT-1 - COAL BASED THERMAL POWER PLANTS
Embodied AI: Ushering in the Next Era of Intelligent Systems
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk

Gradually-Varied Flow in Open Channels ( Lecture notes 06)

  • 1. Prof. Dr. Atıl BULU1 Chapter 6 Gradually-Varied Flow in Open Channels 6.1. Introduction A steady non-uniform flow in a prismatic channel with gradual changes in its water- surface elevation is named as gradually-varied flow (GVF). The backwater produced by a dam or weir across a river and drawdown produced at a sudden drop in a channel are few typical examples of GVF. In a GVF, the velocity varies along the channel and consequently the bed slope, water surface slope, and energy line slope will all differ from each other. Regions of high curvature are excluded in the analysis of this flow. The two basic assumptions involved in the analysis of GVF are: 1. The pressure distribution at any section is assumed to be hydrostatic. This follows from the definition of the flow to have a gradually varied water surface. As gradual changes in the surface curvature give rise to negligible normal accelerations, the departure from the hydrostatic pressure distribution is negligible. 2. The resistance to flow at any depth is assumed to be given by the corresponding uniform flow equation, such as the Manning equation, with the condition that the slope term to be used in the equation is the energy line slope, Se and not the bed slope, S0. Thus, if in a GVF the depth of flow at any section is y, the energy line slope Se is given by, 34 22 R Vn Se = (6.1) where R = hydraulic radius of the section at depth y. 6.2. Basic Differential Equation for the Gradually-Varied Flow Water Surface Profile Since, S0 = Channel bed slope for uniform flow depth y0, dy = Water depth variation for the dx canal reach, d (V2 /2g) = Velocity head variation for the dx reach. Writing the energy equation between the cross-sections 1 and 2 (Fig. 6.1.),
  • 2. Prof. Dr. Atıl BULU2 ( ) ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ +=− ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ +=− ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ +=− +⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ +++=++ g V dy d dx dy SS g V dx d dx dy SS g V ddydxSdxS dxS g V d g V dyy g V ydxS e e e e 2 1 2 2 222 2 0 2 0 2 0 222 0 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ + − = g V dy d SS dx dy e 2 1 2 0 (6.2) Figure 6. 1. ( ) ( )dyyTdA ygA Q g V =→= 2 22 22 ( ) ( ) ( ) ( ) ( ) ( ) 2 2 3 2 4 2 2 2 2 2 2 rF yT yA g V ygA yTQ ygA dy dA AQ ygA Q dy d −=−= −=−=⎥ ⎦ ⎤ ⎢ ⎣ ⎡ Substituting this to Equ. (6.2), 2 0 1 r e F SS dx dy − − = (6.3) S0d d y Sed V y+d V+d V2 /2g+d(V2 /2g E.L S0 Se
  • 3. Prof. Dr. Atıl BULU3 Equ. (6.3) is the general differential equation of the water surface profile for the gradually varied flows. dy/dx gives the variation of water depth along the channel in the flow direction. 6.3. Classification of Flow Surface Profiles For a given channel with a known Q = Discharge, n = Manning coefficient, and S0 = Channel bed slope, yc = critical water depth and y0 = Uniform flow depth can be computed. There are three possible relations between y0 and yc as 1) y0 > yc , 2) y0 < yc , 3) y0 = yc . For horizontal (S0 = 0), and adverse slope ( S0 < 0) channels, 21 0 321 SR n AQ = Horizontal channel, S0 = 0→ Q = 0, Adverse channel , S0 < 0 , Q cannot be computed, For horizontal and adverse slope channels, uniform flow depth y0 does not exist. Based on the information given above, the channels are classified into five categories as indicated in Table (6.1). Table 6.1. Classification of channels Number Channel category Symbol Characteristic condition Remark 1 Mild slope M y0 > yc Subcritical flow at normal depth 2 Steep slope S yc > y0 Supercritical flow at normal depth 3 Critical slope C yc = y0 Critical flow at normal depth 4 Horizontal bed H S0 = 0 Cannot sustain uniform flow 5 Adverse slope A S0 < 0 Cannot sustain uniform flow For each of the five categories of channels, lines representing the critical depth (yc ) and normal depth (y0 ) (if it exists) can be drawn in the longitudinal section. These would divide the whole flow space into three regions as: Region 1: Space above the topmost line, Region 2: Space between top line and the next lower line, Region 3: Space between the second line and the bed.
  • 4. Prof. Dr. Atıl BULU4 Figure (6.2) shows these regions in the various categories of channels. Figure 6.2. Regions of flow profiles Depending upon the channel category and region of flow, the water surface profiles will have characteristics shapes. Whether a given GVF profile will have an increasing or decreasing water depth in the direction of flow will depend upon the term dy/dx in Equ. (6.3) being positive or negative. 2 0 1 r e F SS dx dy − − = (6.3)
  • 5. Prof. Dr. Atıl BULU5 For a given Q, n, and S0 at a channel, y0 = Uniform flow depth, yc = Critical flow depth, y = Non-uniform flow depth. The depth y is measured vertically from the channel bottom, the slope of the water surface dy/dx is relative to this channel bottom. Fig. (6.3) is basic to the prediction of surface profiles from analysis of Equ. (6.3). Figure 6.3 To assist in the determination of flow profiles in various regions, the behavior of dy/dx at certain key depths is noted by studying Equ. (6.3) as follows: 1 1 1 >→< =→= <→> rc rc rc Fyy Fyy Fyy And also, 00 00 00 SSyy SSyy SSyy e e e >→< =→= <→> 1. As 0yy → , 0VV → , 0SSe = 0 0 1 lim 2 00 0 == − − = → consF SS dx dy r yy The water surface approaches the normal depth asymptotically.
  • 6. Prof. Dr. Atıl BULU6 2. As cyy → , 12 =rF , 01 2 =− rF , ∞= − = − − = → 01 lim 0 2 0 e r e yy SS F SS dx dy c The water surface meets the critical depth line vertically. 3. As ∞→y , 000 →→=→= er SFV 0 0 2 0 11 lim S S F SS dx dy r e y == − − = ∞→ The water surface meets a very large depth as a horizontal asymptote. Based on this information, the various possible gradually varied flow profiles are grouped into twelve types (Table 6.2). Table 6.2. Gradually Varied Flow profiles Channel Region Condition Type Mild slope 1 2 3 y > y0 > yc y0 > y > yc y0 > yc > y M1 M2 M3 Steep slope 1 2 3 y> yc > y0 yc > y > y0 yc > y0 > y S1 S2 S3 Critical slope 1 3 y > y0 = y2 y < y0 = yc C1 C3 Horizontal bed 2 3 y > yc y < yc H2 H3 Adverse slope 2 3 y > yc y < yc A2 A3
  • 7. Prof. Dr. Atıl BULU7 6.4. Water Surface Profiles 6.4.1. M – Curves Figure 6.4 General shapes of M curves are given in Fig. (6.4). Asymptotic behaviors of each curve will be examined mathematically. a) M1 – Curve Water surface will be in Region 1 for a mild slope channel and the flow is obviously subcritical. Se < S0 → Mild slope channel y0 > yc → Subcritical flow 2 0 1 r e F SS dx dy − − = Fr < 1 → Subcritical flow → (1 – Fr 2 ) > 0 y > y0 → Se < S0 0> + + = dx dy (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is; Water depth can be between (∞ > y > y0) for Region 1. The asymptotic behaviors of the water surface for the limit values (∞ , y0) are;
  • 8. Prof. Dr. Atıl BULU8 a) y → ∞ , V → 0, Fr → 0, (1 – Fr 2 ) = 1 y → ∞ , V → 0, Se →0 0 0 1 0 lim S S dx dy y = − = ∞→ The water surface meets a very large depth as a horizontal asymptote. b) y → y0 , V → V0 , Se → S0 0 1 lim 2 00 0 = − − = → r yy F SS dx dy The water surface approaches the normal depth asymptotically. The most common of all GVF profiles is the M1 type, which is a subcritical flow condition. Obstructions to flow, such as weirs, dams, control structures and natural features, such as bends, produce M1 backwater curves (Fig. 6.5). These extend to several kilometers upstream before merging with the normal depth. Figure 6.5. M1 Profile b) M2 – Curve Water surface will be in Region 2 for a mild slope channel and the flow is obviously subcritical. (Fig. 6.4). y0 > y > yc y < y0 → V > V0 → Se > S0 → (S0 –Se) < 0 y > yc → Fr < 1 → (1 – Fr 2 ) > 0 −= − + = − − = 2 0 1 r e F SS dx dy (Water depth decrease in the flow direction)
  • 9. Prof. Dr. Atıl BULU9 Asymptotic behavior of the water surface is; y0 > y > yc y → y0, Se → S0 , (S0 –Se) = 0 0 1 lim 2 00 0 = − − = → r yy F SS dx dy The water surface approaches the normal depth asymptotically. y → yc , Fr = 1 , (1 – Fr 2 ) = 0 ∞= − = → 0 lim 0 e yy SS dx dy c The water surface meets the critical depth line vertically. The M2 profiles occur at a sudden drop of the channel, at constriction type of transitions and at the canal outlet into pools (Fig. 6.6). Figure 6.6 c) M3 – Curve Water surface will be in Region 3 for a mild slope channel and the flow is obviously subcritical. (Fig. 6.4). y0 > yc >0 y < y0 → V > V0 → Se > S0 → (S0 – Se) < 0 y < yc → Supercritical flow → Fr > 1 → (1 – Fr 2 ) < 0 += − − = − − = 2 0 1 r e F SS dx dy (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is;
  • 10. Prof. Dr. Atıl BULU10 y → yc → Fr = 1 → (1 – Fr 2 ) = 0 ∞= − = → 0 lim 0 e yy SS dx dy c The water surface meets the critical depth line vertically. y → 0, Se → ∞, (S0 – Se) = ∞ y → 0, ∞== gy V Fr ∞ ∞ = → dx dy y 0 lim (Unknown) The angle of the water surface with the channel bed may be taken as 3 0 0 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ cy y S . Figure 6.7. M3 Profile Where a supercritical stream enters a mild slope channel, M3 type of profile occurs. The flow leading from a spillway or a sluice gate to a mild slope forms a typical example (Fig. 6.7). The beginning of the M3 curve is usually followed by a small stretch of rapidly varied flow and the downstream is generally terminated by a hydraulic jump. Compared to M1 and M2 profiles, M3 curves are of relatively short length.
  • 11. Prof. Dr. Atıl BULU11 Example 6.1 : A rectangular channel with a bottom width of 4.0 m and a bottom slope of 0.0008 has a discharge of 1.50 m3 /sec. In a gradually varied flow in this channel, the depth at a certain location is found to be 0.30 m. assuming n = 0.016, determine the type of GVF profile. Solution: a) To find the normal depth y0, 0 0 00 00 24 4 242 4 y y P A R yyBP yByA + == +=+= == ( ) ( ) 5.0 32 0 35 0 21 0 32 0008.0 24 4 016.0 1 50.1 1 × + ×= ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ == y y S P A n AVAQ ( ) 32 0 35 0 24 0842.0 y y + = By trial and error, y0 = 0.43 m. b) Critical depth yc, m g q y m B Q q c 24.0 81.9 375.0 sec375.0 0.4 50.1 312312 2 =⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ =⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ = === c) Type of profile, y0 = 0.43 m > yc = 0.24 m (Mild slope channel, M profile) y = 0.30 m y0 > y > yc (Region 2) Water surface profile is of the M2 type.
  • 12. Prof. Dr. Atıl BULU12 6.4.2. S – Curves Figure 6.8 General shapes of S curves are given in Fig. (6.8). Asymptotic behaviors of each curve will be examined mathematically. a) S1 – Curve Water surface will be in Region 1 for a steep slope channel and the flow is obviously supercritical. y > yc > y0 y > y0 → V < V0 → Se < S0 → (S0 – Se) > 0 y > yc → Fr < 1 → (1 – Fr 2 ) > 0 += + + = − − = 2 0 1 r e F SS dx dy (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is; ∞ > y > yc y → ∞ , V → 0 , Fr = 0 , (1- Fr 2 ) = 1 y → ∞ , V → 0 , Se = 0 0 0 2 0 01 0 1 lim S S F SS dx dy r e y = − − = − − = ∞→
  • 13. Prof. Dr. Atıl BULU13 The water surface meets a very large depth as a horizontal asymptote. y → yc , Fr = 1 , (1 – Fr 2 ) = 0 ∞= − = − − = → 01 lim 0 2 0 e r e yy SS F SS dx dy c The water surface meets the critical depth line vertically. Figure 6.9. S1 Profile The S1 profile is produced when the flow from a steep channel is terminated by a deep pool created by an obstruction, such as a weir or dam (Fig. 6.9). At the beginning of the curve, the flow changes from the normal depth (supercritical flow) to subcritical flow through a hydraulic jump. The profiles extend downstream with a positive water slope to reach a horizontal asymptote at the pool elevation. b) S2 – Curve Water surface will be in Region 2 for a steep slope channel and the flow is supercritical. y0 < y < yc y > y0 → Se < S0 → (S0 – Se) > 0 y < yc → Supercritical flow → Fr > 1 → (1 – Fr 2 ) < 0 −= − + = − − = 2 0 1 r e F SS dx dy (Water depth will decrease in the flow direction)
  • 14. Prof. Dr. Atıl BULU14 Asymptotic behavior of the water surface is; y0 < y < yc y → yc , Fr → 1 , (1 –Fr 2 ) = 0 ∞= − = − − = → 01 lim 0 2 0 e r e yy SS F SS dx dy c The water surface meets the critical depth line vertically. y → y0 , V → V0 , Se → S0 , (S0 – Se) = 0 0 1 0 lim 2 0 = − = → r yy Fdx dy The water surface approaches the normal depth asymptotically. Figure 6.10. S2 Profile Profiles of the S2 type occur at the entrance region of a steep channel leading from a reservoir and at a brake of grade from mild slopes to steep slope (Fig. 6.10). Generally S2 profiles are short of length. c) S3 – Curve Water surface will be in Region 3 for a steep slope channel and the flow is supercritical. 0 < y < y0 y < y0 → Se > S0 → (S0 – Se) < 0 y < yc → Supercritical flow → Fr > 1 → (1 – Fr 2 ) < 0 += − − = − − = 2 0 1 r e F SS dx dy (Water depth will increase in the flow direction)
  • 15. Prof. Dr. Atıl BULU15 Asymptotic behavior of the water surface is; 0 < y < y0 y → 0, Se → ∞ , (S0 – Se) = -∞ y → 0 , ∞== gy V Fr ∞ ∞− = → dx dy y 0 lim (Unknown) The angle of the water surface with the channel bed may be taken as 3 0 0 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ cy y S . y → y0 , Se → S0 , (S0 – Se) = 0 0 1 0 lim 2 0 = − = → r yy Fdx dy The water surface approaches the normal depth y0 asymptotically. Figure 6.11. S3 Profile Figure 6.12. S3 Profile Free flow from a sluice gate with a steep slope on its downstream is of the S3 type (Fig. 6.11). The S3 curve also results when a flow exists from a steeper slope to a less steep slope (Fig. 6.12).
  • 16. Prof. Dr. Atıl BULU16 6.4.3. C – Curves Figure 6.13 General shapes of C curves are given in Fig. (6.13). Asymptotic behaviors of each curve will be examined mathematically. Since the flow is at critical stage, y0 = yc, there is no Region 2. a) C1 – Curve Water surface will be in Region 1 for a critical slope channel. y0 = yc < y < ∞ y > yc → Se < S0 → (S0 –Se) > 0 y > yc → Subcritical flow → Fr < 1 → (1 –Fr 2 ) > 0 += + + = − − = 2 0 1 r e F SS dx dy (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is; y → ∞ , Se → 0 , ( S0 – Se) = S0 y → ∞, V = 0 → Fr = 0 → (1 – Fr 2 ) = 1 0 0 01 0 lim S S dx dy y = − − = ∞→ The water surface meets a very large depth as a horizontal asymptote. y → yc , Fr = 1 → (1 – Fr 2 ) = 0 y → yc = y0 , Se = S0 = Sc → (S0 –Se ) =0 0 0 lim 0 = =→ dx dy yyy c (Unknown) C1 and C3 profiles are very rare and highly unstable
  • 17. Prof. Dr. Atıl BULU17 6.4.4. H – Curves Figure 6.14 General shapes of H curves are given in Fig. (6.14). For horizontal slope channels, uniform flow depth y0 does not exist. Critical water depth can be computed for a given discharge Q and therefore critical water depth line can be drawn. Since there is no uniform water depth y0, Region 1 does not exist. a) H2 – Curve Water surface will be in Region 2 for a horizontal slope channel. ∞ > y > yc y > yc → Se < S0 = 0 → (S0 –Se) < 0 y > yc → subcritical flow → Fr < 1 → (1 – Fr 2 ) > 0 −= + − = dx dy (Water depth will decrease in the flow direction) Asymptotic behavior of the water surface is; y → ∞ , Se → 0 , ( S0 – Se) = S0 y → ∞, V = 0 → Fr = 0 → (1 – Fr 2 ) = 1 0 0 01 0 lim S S dx dy y = − − = ∞→ The water surface meets a very large depth as a horizontal asymptote. y → yc , Fr = 1 → (1 – Fr 2 ) = 0 y → yc , Se = Sc → S0 = 0, (S0 –Se ) = -Se ∞= − = − − = → 0 0 1 lim 2 0 e r e yy S F SS dx dy c The water surface meets the critical depth line vertically.
  • 18. Prof. Dr. Atıl BULU18 b) H3 – Curve 0 > y > yc y < yc → Se > S0 = 0 → (S0 –Se) = -Se y → yc , supercritical flow , Fr > 1 → (1 – Fr 2 ) < 0 += − − = − − = 2 0 1 r e F SS dx dy (Water depth will increase in the flow direction) Asymptotic behavior of the water surface is; y → 0, Se → ∞ , (S0 – Se) = -∞ y → 0 , ∞== gy V Fr ∞ ∞− = → dx dy y 0 lim (Unknown) y → yc → Fr = 1 → (1 – Fr 2 ) = 0 y → yc , Se = Sc → S0 = 0, (S0 –Se ) = -Se ∞= − = − − = → 0 0 1 lim 2 0 e r e yy S F SS dx dy c The water surface meets the critical depth line vertically. Figure 6.15 A horizontal channel can be considered as the lower limit reached by a mild slope as its bed slope becomes flatter. The H2 and H3 profiles are similar to M2 and M3 profiles respectively (Fig. 6.15). However, the H2 curve has a horizontal asymptote.
  • 19. Prof. Dr. Atıl BULU19 6.4.5. A– Curves Figure 6.16 General shapes of A curves are given in Fig. (6.16). For adverse slope channels, uniform flow depth y0 does not exist. Critical water depth can be computed for a given discharge Q and therefore critical water depth line can be drawn. Since there is no uniform water depth y0, Region 1 does not exist as well as in A curves. A2 and A3 curves are similar to H2 and H3 curves respectively. Figure 6.17 Adverse slopes are rather rare (Fig. 6.17). These profiles are of very short length. 6.5. Control Sections A control section is defined as a section in which a fixed relationship exists between the discharge and depth of flow. Weirs, spillways, sluice gates are some typical examples of structures which give rise to control sections. The critical depth is also a control point. However, it is effective in a flow profile which changes from subcritical to supercritical flow. In the reverse case of transition from supercritical flow to subcritical flow, a hydraulic jump is usually formed bypassing the critical depth as a control point. Any GVF profile will have at least one control section.
  • 20. Prof. Dr. Atıl BULU20 (d) Figure 6.18 In the synthesis of GVF profiles occurring in a serially connected channel elements, the control sections provide a key to the identification of proper profile shapes. A few typical control sections are shown in Fig. (6.18, a-e). It may be noted that subcritical flows have controls in the downstream end, while supercritical flows are governed by control sections existing at the upstream end of the channel section. In Figs. (6.18 a and b) for the M1 profile, the control section (indicated by a dark dot in the figures) is just at the upstream of the spillway and sluice gate respectively. In Figs. (6.18 b and d) for M3 and S3 profiles respectively, the control point is at the vena contracta of the sluice gate flow. In subcritical flow reservoir offtakes (Fig. 6.18 c), even though the discharge is governed
  • 21. Prof. Dr. Atıl BULU21 by the reservoir elevation, the channel entry section is not strictly a control section. The water surface elevation in the channel, will be lower than the reservoir elevation by a headloss amount equivalent to (1 + ζ )V2 /2g where ζ is the entrance loss coefficient. The true control section will be at a downstream location in the channel. For the situation shown in Fig. (6.18 c) the critical depth at the free overflow at the channel end acts as the downstream control. For a sudden drop (free overflow) due to curvature of the streamlines the critical depth usually occurs at distance of about 4yc upstream of the drop. This distance, being small compared to GVF lengths, is neglected and it is usual to perform calculations by assuming yc to occur at the drop. (e) For a supercritical canal intake (Fig. 6.18 e), the reservoir water surface falls to the critical depth at the head of the canal and then onwards the water surface follows the S2 curve. The critical depth occurring at the upstream end of the channel is the control for this flow. Figure 6.19
  • 22. Prof. Dr. Atıl BULU22 6.5.1. Feeding a Pool for Subcritical Flow A mild slope channel discharging into a pool of variable surface elevation is indicated in Fig. (6.19). Four cases are shown. 1) The pool elevation is higher than the elevation of the normal (uniform) water depth line at B. This gives rise to a drowning of the channel end. A profile of the M1 type is produced with the pool level at B as control. 2) The pool elevation is lower than the elevation of the normal depth line but higher than the critical depth line at B. The pool elevation acts as a control for the M2 curve. 3) The pool elevation has dropped down to that of the critical depth line at B and the control is still at the pool elevation. 4) The pool elevation has dropped lower than the elevation of the critical depth line at B. The water surface cannot pass through a critical depth at any location other than B and hence a sudden drop in the water surface at B is observed. The critical depth at B is the control of this flow. 6.5.2. Feeding a Pool for Supercritical Flow When the slope of the channel is steep, that is greater than the critical slope, the flow in the channel becomes supercritical. In practical applications, steep channels are usually short, such as raft and log chutes that are used as spillways. Figure 6.20 As the control section in a channel of supercritical flow is at the upstream end, the delivery of the channel is fully governed by the critical discharge at section 1.
  • 23. Prof. Dr. Atıl BULU23 1) When the tailwater level B is less than the outlet depth (normal depth) at section 2, the flow in the canal is unaffected by the tailwater. The flow profile passes through the critical water depth near C and approaches the normal depth by means of a smooth drawdown curve of the S2 type. 2) When the tailwater level B is greater than the outlet depth, the tailwater will raise the water level in the upstream portion of the canal to form an S1 profile between j and b`, producing a hydraulic jump at the end j of the profile. However, the flow upstream from the jump will not be affected by the tailwater. Such a hydraulic jump is objectionable and dangerous, particularly when the canal is raft chute or some other structure intended to transport a floating raft from the upstream reservoir to a downstream pool. A neutralizing reach may be suggested for the solution of this problem. (Chow, 1959). (Fig. 6.21) Figure 6.21. Elimination of hydraulic jump by neutralizing reach. In this reach the bottom slope of the channel is made equal to the critical slope. According to a corresponding case of C1 profile in Fig. (6.21), the tailwater levels will be approximately horizontal lines which intersect the surface of flow in the canal without causing any disturbance. At the point of intersection, theoretically, there is a jump of zero height. 3) As the tailwater rises further, the jump will move upstream, maintaining its height and form in the uniform flow zone nb, until it reaches point n. The height of the jump becomes zero when it reaches the critical depth at c. In the meantime, the flow profile reaches its theoretical limit cb`` of the S1 profile. Beyond this limit the incoming flow will be directly affected by the tailwater. In practical applications, the horizontal line cb``` may be taken as the practical limit of the tailwater stage.
  • 24. Prof. Dr. Atıl BULU24 6.6. Analysis of Flow Profile A channel carrying a gradually varied flow can in general contain different prismoidal channel cross-sections of varying hydraulic properties. There can be a number of control sections at various locations. To determine the resulting water surface profile in a given case, one should be in a position to analyze the effects of various channel sections and controls connected in series. Simple cases are illustrated to provide information and experience to handle more complex cases. 6.6.1. Break in Grades Simple situations of a series combination of two channel sections with differing bed slopes are considered. In Fig. (6.22.a), a break in grade from a mild channel to a milder channel is shown. It is necessary to first draw the critical-depth line (CDL) and the normal-depth line (NDL) for both slopes. Since yc does not depend upon the slope for a taken Q = discharge, the CDL is at a constant height above the channel bed in both slopes. The normal depth y01 for the mild slope is lower than that of the of the milder slope (y02). In this case, y02 acts as a control, similar to the weir or spillway case and an M1 backwater curve is produced in the mild slope channel. Various combinations of slopes and the resulting GVF profiles are presented in Fig. ( 6.22, a-h). It may be noted that in some situations there can be more than one possible profiles. For example, in Fig. (6.22 e), a jump and S1 profile or an M3 profile and a jump possible. The particular curve in this case depends on the channel and its flow properties. In the examples indicated in Fig. (6.22), the section where the grade changes acts a control section and this can be classified as a natural control. It should be noted that even though the bed slope is considered as the only variable in the above examples, the same type of analysis would hold good for channel sections in which there is a marked change in the roughness characteristics with or without change in the bed slope. A long reach of unlined canal followed by a line reach serves as a typical example for the same. A change in the channel geometry (the bed width or side slope) beyond a section while retaining the prismoidal nature in each reach also leads to a natural control section.
  • 25. Prof. Dr. Atıl BULU25
  • 26. Prof. Dr. Atıl BULU26
  • 27. Prof. Dr. Atıl BULU27 Figure 6.22 6.6.2. Serial Combination of Channel Sections To analyze a general problem of many channel sections and controls, the following steps are to be applied. 1. Draw the longitudinal section of the system. 2. Calculate the critical depth and normal depths of various reaches and draw the CDL and NDL in all reaches. 3. Mark all the controls, both the imposed as well as natural controls. 4. Identify the possible profiles. Example 6. 2: Identify and sketch the GVF profiles in three mild slopes which could be described as mild, steeper mild and milder. The three slopes are in series. The last slope has a sluice gate in the middle of the reach and the downstream end of the channel has a free overfall. Solution: The longitudinal profile of the channel, critical depth line and normal depth lines for the various reaches are shown in the Fig. (6.23). The free overfall at E is obviously a control. The vena contracta downstream of the sluice gate at D is another control. Since for subcritical flow the control is at the downstream end of the channel, the higher of the two normal depths at C acts as a control for the reach CB, giving rise to an M1 profile over CB. At B, the normal depth of the channel CB acts as a control giving rise to an M2 profile over AB. The controls are marked distinctly in Fig. (6.23). With these controls the possible flow profiles are: an M2 profile on channel AB, M1 profile on channel BC, M3 profile and M2 profile through a jump on the stretch DE.
  • 28. Prof. Dr. Atıl BULU28 Figure 6.23 Example 6.3: A trapezoidal channel has three reaches A, B, and C connected in series with the following physical characteristics. Reach Bed width (B) Side slope (m) Bed slope S0 n 1 2 3 4.0 m 4.0 m 4.0 m 1.0 1.0 1.0 0.0004 0.009 0.004 0.015 0.012 0.015 For a discharge Q = 22.5 m3 /sec through this channel, sketch the resulting water surface profiles. The length of the reaches can be assumed to be sufficiently long for the GVF profiles to develop fully. Solution: The normal depths and critical water depths in the various reaches are calculated as: ( ) myBT myBP ymyBA 2 12 2 += ++= += B=4m y m 1 T
  • 29. Prof. Dr. Atıl BULU29 yT yyP yyA 24 241124 )4( 232 += +=++= += ( )[ ] ( ) 32 0 23 5.0 0 35 00 5.0 0 32 24 4 5.22 y S n yy S P A n A AVQ ×+ × + = ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ == Uniform flow depths for the given data for every reach are calculated by trial and error method; Reach A: S0A = 0.0004 , nA=0.015 → y0A = 2.26 m Reach B: S0B = 0.009 , nB =0.012 → y0B = 0.81 m Reach C: S0C = 0.004 , nC = 0.015 → y0C =1.17 m Since the channel is prismoidal (the geometry does not change in the reaches), there will be only one critical water depth. ( ) ( )[ ] 1 481.9 245.22 1 3 2 3 2 = +× +× = cc c c c yy y gA TQ yc =1.32 m Reach A is a mild slope channel as y0A = 2.26 m > y0c = 1.32 m and the flow is subcritical. Reach B and C are steep slope channels as y0B, y0C < yc and the flow is supercritical on both reaches. Reach B is steeper than reach C. The various reaches are schematically shown in Fig. (6.24). The CDL is drawn at a height of 1.32 m above the bed level and NDLs are drawn at the calculated y0 values. The controls are marked in the figure. Reach A will have an M2 drawdown curve, reach B an S2 drawdown curve and reach C a rising curve as shown in the figure. It may be noted that the resulting profile as above is a serial combination of Fig. (6.22 d and f).
  • 30. Prof. Dr. Atıl BULU30 Figure 6.23