Open Channel Flow
Open Channel Flow
Open Channel Flow
(670441)
Philadelphia University
Faculty of Engineering
Civil Engineering Department
First Semester, 2011/2012
Lecture Title
Water Flow in
Open Channels
3 3
Introduction
An open channel is a duct in which the liquid flows
with a free surface.
Open channel hydraulics is of great importance in civil
engineers, it deals with flows having a free surface,
for example:
Channels constructed for water supply, irrigation,
drainage, and
Sewers, culverts, and
Tunnels flowing partially full; and
Natural streams and rivers.
4
Pipe Flow and Open Channel Flow
Pipe Flow
The liquid completely fills
the pipe and flow under
pressure.
The flow in a pipe takes
place due to difference of
pressure (pressure
gradient),
The flow in a closed
conduit is not necessarily a
pipe flow.
Open Channel Flow
Flow takes place due to
the slope of the channel
bed (due to gravity).
The flow must be
classified as open
channel flow if the liquid
has a free surface.
5
Pipe Flow Open Channel Flow
6
For Pipe flow (Fig. a):
The hydraulic gradient line (HGL) is the sum of the elevation and the
pressure head (connecting the water surfaces in piezometers).
The energy gradient line (EGL) is the sum of the HGL and velocity
head.
The amount of energy loss when the liquid flows from section 1 to
section 2 is indicated by h
L
.
For open channel flow (Fig. b):
The hydraulic gradient line (HGL) corresponds to the water surface
line (WSL); where it subjected to only atmospheric pressure which is
commonly referred to as the zero pressure reference.
The energy gradient line (EGL) is the sum of the HGL and velocity
head.
The amount of energy loss when the liquid flows from section 1 to
section 2 is indicated by h
L
. For uniform flow in an open channel, this
drop in the EGL is equal to the drop in the channel bed.
7
6.1 Classifications of Open Channel Flow
Classification based on the time criterion:
1. Steady Flow (time independent)
(discharge and water depth do not change with time)
2. Unsteady Flow (time dependent)
(discharge and water depth at any section change with time)
Classification based on the space criterion:
1. Uniform flow (are mostly steady)
(discharge and water depth remains the same at every section in
the channel)
2. Non-uniform Flow
(discharge and water depth change at any section in the channel)
8
Non-uniform flow is also called varied flow ( the flow
in which the water depth and or discharge change
along the length of the channel), it can be further
classified as:
Gradually varied flow (GVF) where the depth of the
flow changes gradually along the length of the channel.
Rapidly varied flow (RVF) where the depth of flow
changes suddenly over a small length of the channel.
9
a) Uniform flow are mostly
steady
b) Unsteady uniform flows
are very rare in nature
c) Steady varied flow
(over a spillway crest)
d) Unsteady varied flow (flood wave)
e) Unsteady varied flow (tidal surge)
10
11
6.2 Uniform Flow in Open Channel
Uniform flow in an open channel must satisfy the following main
features:
1. The water depth y, flow area A, discharge Q, and the velocity
distribution V at all sections throughout the entire channel length must
remain constant.
2. The slope of the energy gradient line (S
e
), the water surface slope
(S
ws
), and the channel bed slope (S
0
) are equal.
S
e
= S
ws
= S
0
12
This is possible when the gravity force (W sin u) component equal
the resistance to the flow (F
f
)
|
.
|
\
|
|
.
|
\
|
= =
0
2
0
. . ) ( S
P
A
K
V PL KV ALS
+ +
=
=
20
6.3 Hydraulic Efficiency of open channel sections
Based on their existence, an open channel can be natural
or artificial:
Natural channels: such as streams, rivers, valleys, etc.
These are generally irregular in shape, alignment and
roughness of the surface.
Artificial channels: built for some specific purpose, such
as irrigation, water supply, wastewater, water power
development, and rain collection channels. These are
regular in shape and alignment with uniform roughness of
the boundary surface.
21
Based on their shape, an open channel can be prismatic or non-
prismatic:
Prismatic channels: the cross section is uniform and the bed slop is
constant.
Non-prismatic channels: when either the cross section or the slope
(or both) change, the channel is referred to as non-prismatic. It is
obvious that only artificial channel can be prismatic.
The most common shapes of prismatic channels are rectangular,
parabolic, triangular, trapezoidal and circular; see the next figure.
22
Most common shapes of prismatic channels
23
Most economical section is called the best hydraulic
section or most efficient section as the discharge,
passing through a given cross-sectional area A, slope of
the bed S
0
and a resistance coefficient, is maximum.
Hence the discharge Q will be maximum when the
wetted perimeter P is minimum.
P
const S
P
A
C A S R C A V A Q
e e h
1
* . = = = =
24
Economical Rectangular Channel
D, B A =
B D 2 P + =
D
A
2D P + =
0
dD
dP
=
2 2 2
2 0 2
D
D B
D
A
D
A
dD
dP
= = =
|
.
|
\
|
=
D
B
2 =
2
B
D =
P should be minimum for a given area;
D
D
D D
D D
D B
D B
P
A
R
h
4
2
2 2
2
2
2
=
+
=
+
= =
So, the rectangular channel will be most economical when either:
the depth of the flow is half the width, or
the hydraulic radius is half the depth of flow.
2
D
R
h
=
25
Economical Trapezoidal Channel
D ) D n (B A + =
2
1 2 n D B P + + =
D n
D
A
B =
2
1 2 n D ) D n
D
A
( P + + =
0
dD
dP
= = + + = 0 1 2
2
2
n n
D
A
dD
dP
n
D
A
n 1 2
2
2
+ = +
D
D n B
n
D
D nD) (B
n
2
1 2
2
2
+
= +
+
= +
2
D n 2 B
n 1 D
2
+
= +
or
P B B nD B nD = + + = + 2 2( )
R
A
P
B n D D
B n D
h
= =
+
+
( )
( ) 2
R
D
h
=
2
26
Other criteria for economic Trapezoidal section
When a semi-circle is drawn with the trapezoidal center,
O, on the water surface and radius equal to the depth of
flow, D, the three sides of the channel are tangential to
the semi-circle.
To prove this condition, using the figure shown, we have:
OF OM B n D
B
n D = = + = + sin ( ) sin ( ) sin o o o
1
2
2
2
sino =
+
OF
B
n D
2
using triangle KMN, we have:
sino = =
+
MK
MN
D
D n 1
2
27
OF
B
n D
n
=
+
+
( )
2
1
2
D n
B n D
1
2
2
2
+ =
+
using equation to replace the numerator , we obtain:
OF
D n
n
=
+
+
1
1
2
2
OF D =
Thus, if a semi-circle is drawn with O as center and radius
equal to the depth of flow D, the three sides of a most
economical trapezoidal section will be tangential to the semi-
circle.
28
The best side slope for Trapezoidal section
when n =
1
3
o
= u 60
B D n n = + 2 1
2
( )
+
= +
2
2
1
2
D n B
n D
+ = + + = ) ( 2 2 D n B D n B B P
B
A
D
n D =
A
D
n D D n n = + 2 1
2
( )
D
A
n n
2
2
2 1
=
+
29
Now, from equations:
P B nD = + 2( )
B
A
D
n D =
P
A
D
= 2
P
A
D
A n n
2 2 2
4 4 2 1 = = + ( ) ( )
squaring both sides
= 0
dn
dP
2 4 1 2 1
2
1
2
P
dP
dn
A n n = +
[( ) *( ) ]
2
1
1
2
n
n +
=
+ =
2 2
1 4 n n
n =
1
3
1 3
1 n
= = tanu
= u
o
60
The best side slope is at 60
o
to the horizontal, i.e.; of all trapezoidal
sections a half hexagon is most economical. However, because of
constructional difficulties, it may not be practical to adopt the most
economical side slopes
30
Circular section
o
o
2 sin
8 4
2 2
d d
A =
d r P o o 2 = =
In the case of circular channels, the area of the flow cannot be
maintained constant. Indeed, the cross-sectional area A and the wetted
perimeter P both do not depend on D but they depend on the angle o.
Referring to the figure shown, we
can determine the wetted
perimeter P and the area of flow
A as follows:
Thus in case of circular channels, for most economical section, two
separate conditions are obtained:
1. Condition for maximum discharge, and
2. Condition for maximum velocity.
31
1. Condition for Maximum Discharge for Circular Section:
o
o
= 154 D d = 095 .
(Using Mannings formula)
o
o
= 151 D d = 094 .
2. Condition for Maximum Velocity for Circular Section:
(Using the Chezy formula)
Q AV AC R S C
A
P
S
h
= = =
3
Q C S
A
P
2 2
3
=
|
\
|
.
|
dQ
do
= 0
V C R S C
A
P
S
h
= =
V C S
A
P
2 2
=
|
\
|
.
|
dV
do
= 0
o
o
= 128 75 .
D d = 081 .
32
Variation of flow and velocity with depth in circular pipes
33
6.4 Energy Principle in Open Channel Flow
g
V
y Z Energy Total
2
2
+ + =
The total energy of a flowing liquid per unit weight is given by:
If the channel bed is taken as the datum, then the total energy per unit
weight will be:
g
V
y E
specific
2
2
+ =
Specific energy (Es) of a flowing liquid in a channel is defined as
energy per unit weight of the liquid measured from the channel bed as
datum. It is a very useful concept in the study of open channel flow.
34
E y
V
g
E E
s p k
= + = +
2
2
E
p
= potential energy of flow = y
E
k
= kinetic energy of flow =
g
V
2
2
E y
Q
g A
s
= +
2
2
2
Valid for any cross section
Specific Energy Curve:
It is defined as the
curve which shows the
variation of specific
energy (Es ) with depth
of flow y.
35
Specific Energy Curve (Rectangular channel)
Consider a rectangular channel in which a constant discharge
q = discharge per unit width = = constant ( since Q and B are constants)
Q
B
=
= =
y
q
y B
Q
A
Q
V
k p s
E E
y g
q
y E + = + =
2
2
2
p
E
E
K
E
P
E
s
y
c
36
Sub-critical, critical, and supercritical flow
The criterion used in this classification is what is known by Froude number, Fr, which
is the measure of the relative effects of inertia forces to gravity force:
h
r
D g
V
F =
T
A
Width Surface Water
Area) (Wetted Flow of Area
= =
h
D
T
T
Flow Fr
Sub-critical Fr < 1
Critical 1 = Fr
Supercritical Fr >1
g A
T Q
F
r
3
2
2
=
V = mean velocity of flow of water,
D
h
= hydraulic depth of the channel
37
c
y E
min
E y and y
s1 1 1
'
Referring to the energy curve, the following features can be observed:
1. The depth of flow at point C is referred to as critical depth, y
c
.
(It is defined as that depth of flow of liquid at which the specific energy is
minimum, The flow that corresponds to this point is called
critical flow (Fr = 1.0).
2. For values of E
s
greater than E
min
, there are two corresponding depths.
One depth is greater than the critical depth and the other is smaller then
the critical depth, for example;
These two depths for a given specific energy are called the alternate depths.
y y
c
>
y y
c
<
3. If the flow depth
In this case E
s
increases as y increases.
In this case E
s
decreases as y increases.
the flow is said to be sub-critical (Fr < 1.0).
4. If the flow depth
the flow is said to be super-critical (Fr > 1.0).
38
39
Critical depth, y
c
for rectangular channel
=0
dy
dE
2
2
2 y g
q
y E
s
+ =
Critical depth, y
c
, is defined as that depth of flow of liquid at which the
specific energy is minimum, E
min
,
The mathematical expression for critical depth is obtained by differentiating
energy equation with respect to y and equating the result to zero;
d
dy
y
q
g y
q
g
y
( ) ( ) + = +
=
2
2
2
3
2
1
2
2
0
1 0
2
3
3
2
= =
q
g y
y
q
g
3
1
2
|
|
.
|
\
|
=
g
q
y
c
40
Critical velocity, V
c
for rectangular channel
,
2
3
g
q
y
c
=
y
q
y B
Q
A
Q
V =
= =
V
q
y
c
c
=
OR
=
g
y V
y
c c
c
2 2
3
V g y
c c
=
r
c
c
F
y g
V
= =
1
41
Minimum Specific Energy in terms of critical depth
E y
q
g y
c
c
min
= +
2
2
2
g
q
y
c
2
3
=
E y
y
c
c
min
= +
2
E
y
c
min
=
3
2
y
E
c
=
2
3
min
OR
42
Critical depth, y
c
, for Non- Rectangular Channels
=0
dy
dE
s
0 ) (
2
2
1 )
2
(
3
2
2
2
= = +
dy
dA
A g
Q
A g
Q
y
dy
d
(constant discharge is assumed)
1 0
2
3
=
Q
g A
dA
dy
( )
dA/dy = the rate of increase of area with respect to y = T (top width).
OR
= 0 1
3
2
A g
T Q
Q
g
A
T
2 3
=
condition must be satisfied for the flow
at the critical depth.
Recalling that
=
T
A
D
h
Q
g
A D
h
2
2
=
The equation may also be written in terms of velocity
V
g
D
h
2
2 2
=
The velocity head is equal to one-half the hydraulic depth for critical flow.
43
+ =
2
2
2 A g
Q
y E
s
E y
A
T
s
= +
2
This equation represents
the critical state
E y
A
T
c c
= +
1
2
( )
OR
The general equation for the specific energy in
critical state applicable to channels of all shapes.
Trapezoidal section
Circular section Triangle section
Rectangular section
E
B n y y
B n y
c
c c
c
=
+
+
( )
( )
3 5
2 2
E
d d
c
= +
2
1
16
2 2
( cos )
( sin )
sin
o
o o
o
E y
c c
=
5
4
2
3
c
c
y
E =
44
Constant Specific Energy
The specific energy was varied and the discharge was assumed to be
constant. Let us now consider the case in which the specific energy is
kept constant and the discharge Q is varied.
=
2
2
2 A g
Q
y E
s
Q A g E y
s
= 2 ( )
Q A g E y gA E gA y
s s
2 2 2 2
2 2 2 = = ( ) ( )
The discharge will maximum if
dQ
dy
= 0
Q
dQ
dy
g E A
dA
dy
g y A
dA
dy
A
s
= +
(
2 2 2 2
2
( ) ( )
dA/dy = T
0 2 ) 2 ( 2 ) 2 ( 2
2
= gA yAT g AT E g
s
45
4 4 2 0 E T yT A
s
=
= A y E T
s
) ( 2
E y
A
T
s
= +
2
but E y
Q
g A
s
= +
2
2
2
+ = +
T
A
y
A g
Q
y
2 2
2
2
Q
g
A
T
2 3
=
Thus for a given specific energy, the discharge in a given channel is a
maximum when the flow is in the critical state. The depth corresponding
to the maximum discharge is the critical depth.
46
6.5 Hydraulic Jump
A hydraulic jump occurs when flow changes from a supercritical flow
(unstable) to a sub-critical flow (stable).
There is a sudden rise in water level at the point where the hydraulic
jump occurs.
Rollers (eddies) of turbulent water form at this point. These rollers cause
dissipation of energy.
A hydraulic jump occurs in practice at the toe of a dam or below a sluice gate
where the velocity is very high.
47
General Expression for Hydraulic Jump:
In the analysis of hydraulic jumps, the following assumptions are made:
(1) The length of hydraulic jump is small. Consequently, the loss of head
due to friction is negligible.
(2) The flow is uniform and pressure distribution is due to hydrostatic
before and after the jump.
(3) The slope of the bed of the channel is very small, so that the
component of the weight of the fluid in the direction of the flow is
neglected.
48
Location of hydraulic jump
Generally, a hydraulic jump occurs when the flow changes from
supercritical to subcritical flow.
The most typical cases for the location of hydraulic jump are:
1. Jump below a sluice gate.
2. Jump at the toe of a spillway.
3. Jump at a glacis.
(glacis is the name given to sloping floors provided in hydraulic structures.)
49
The net force in the direction of flow = the rate of change of moment in that direction
) (
1 2
V V
g
Q
=
The net force in the direction of the flow, neglecting frictional resistance and the
component of weight of water in the direction of flow,
R = F
1
- F
2
.
Therefore, the impulse-moment yields
F F
Q
g
V V
1 2 2 1
=
( )
Where F
1
and F
2
are the pressure forces at section 1 and 2, respectively.
A y A y
Q
g
V V
1 1 2 2 2 1
= ( )
A y A y
Q
g A A
1 1 2 2
2
2 1
1 1
= ( )
Q
gA
A y
Q
gA
A y
2
1
1 1
2
2
2 2
+ = +
y
= the distance from the water surface to the centroid of the flow area
50
Comments:
This is the general equation governing the hydraulic jump for any
shape of channel.
The sum of two terms is called specific force (M). So, the equation can
be written as:
M
1
= M
2
This equation shows that the specific force before the hydraulic jump
is equal to that after the jump.
Q
gA
A y
Q
gA
A y
2
1
1 1
2
2
2 2
+ = +
51
Hydraulic Jump in Rectangular Channels
A By
1 1
=
y
y
1
1
2
= A B y
2 2
=
y
y
2
2
2
=
Q
g B y
By
y Q
g B y
By
y
2
1
1
1
2
2
2
2
2 2
+ = + ( )( ) ( )( )
Q
gA
A y
Q
gA
A y
2
1
1 1
2
2
2 2
+ = +
q
Q
B
=
q
g
y y
y y
y y
2
2 1
1 2
2
2
1
2
2
|
\
|
.
|
=
2
2
1 2 2 1
q
g
y y y y = + ( )
y y y y
q
g
2 1
2
2
2
1
2
2
0 + =
using
, we get
52
y
y y q
g y
2
1 1
2
2
1
2 2
2
= +
|
\
|
.
|
+
|
\
|
.
|
y
y y q
g y
1
2 2
2
2
2
2 2
2
= +
|
\
|
.
|
+
|
\
|
.
|
This is a quadratic equation, the solution of which may be written as:
where y
1
is the initial depth and y
2
is called the conjugate depth. Both are called
conjugate depths.
These equations can be used to get the various characteristics of hydraulic jump.
y
y
q
g y
2
1
2
1
3
1
2
1 1
8
= + +
|
\
|
.
|
|
y
y
q
g y
1
2
2
2
3
1
2
1 1
8
= + +
|
\
|
.
|
|
53
y
q
g
c
3
2
=
|
|
.
|
\
|
|
|
.
|
\
|
+ + =
3
1 1
2
8 1 1
2
1
y
y
y
y
c
y
y
y
y
c 1
2 2
3
1
2
1 1 8 = + +
|
\
|
.
|
|
\
|
.
|
|
|
( )
2
1
1
2
8 1 1
2
1
F
y
y
+ + =
y
y
F
1
2
2
2
1
2
1 1 8 = + +
|
\
|
.
|
1
1
1
y g
V
F=
F
V
g y
2
2
2
=
But for rectangular channels, we have
Therefore,
These equations can also be written in terms of Froudes number as:
54
H E E E
L
= = A
1 2
E y
q
g y
s
= +
2
2
2
H y
q
g y
y
q
g y
L
= + +
|
\
|
.
|
1
2
1
2
2
2
2
3
2 2
Due to the turbulent flow in hydraulic jump, a dissipation (loss) of energy
occurs:
Where, E = specific energy
For rectangular channels:
hence,
Head Loss in a hydraulic jump (H
L
):
After simplifying, we obtain
2 1
3
1 2
4
) (
y y
y y
H E
L
= = A
55
h y y
j
=
2 1
L h
j j
~ 6
Height of hydraulic jump (h
j
):
The difference of depths before and after the jump is known as the
height of the jump,
Length of hydraulic jump (L
j
):
The distance between the front face of the jump to a point on the
downstream where the rollers (eddies) terminate and the flow becomes
uniform is known as the length of the hydraulic jump. The length of the
jump varies from 5 to 7 times its height. An average value is usually
taken:
56
6.6 Gradually Varied Flow
Non-uniform flow is a flow for which the depth of flow is varied.
This varied flow can be either Gradually varied flow (GVF) or
Rapidly varied flow (RVF).
Such situations occur when:
- control structures are used in the channel or,
- when any obstruction is found in the channel,
- when a sharp change in the channel slope takes place.
57
Classification of Channel-Bed Slopes
The slope of the channel bed is very important in determining the
characteristics of the flow.
Let
S
0
: the slope of the channel bed ,
S
c
: the critical slope or the slope of the channel that
sustains a given discharge (Q) as uniform flow at the critical
depth (y
c
).
y
n
is is the normal depth when the discharge Q flows as
uniform flow on slope S
0
.
58
S S or y y
c n c 0
= =
S S or y y
c n c 0
< >
S S or y y
c n c 0
> <
S
0
00 = .
S negative
0
=
The slope of the channel bed can be classified as:
1) Critical Slope C : the bottom slope of the channel is equal to the critical slope.
2) Mild Slope M : the bottom slope of the channel is less than the critical slope.
3) Steep Slope S : the bottom slope of the channel is greater than the critical slope.
4) Horizontal Slope H : the bottom slope of the channel is equal to zero.
5) Adverse Slope A : the bottom slope of the channel rises in the direction of the
flow (slope is opposite to direction of flow).
59
60
Classification of Flow Profiles (water surface profiles)
The surface curves of water are called flow profiles (or water surface
profiles).
The shape of water surface profiles is mainly determined by the slope of
the channel bed S
o
.
For a given discharge, the normal depth y
n
and the critical depth y
c
may be calculated. Then the following steps are followed to classify the
flow profiles:
1- A line parallel to the channel bottom with a height of y
n
is drawn and is
designated as the normal depth line (N.D.L.)
2- A line parallel to the channel bottom with a height of y
c
is drawn and is
designated as the critical depth line (C.D.L.)
3- The vertical space in a longitudinal section is divided into 3 zones
using the two lines drawn in steps 1 & 2 (see the next figure)
61
4- Depending upon the zone and the slope of the bed, the water profiles
are classified into 13 types as follows:
(a) Mild slope curves M
1
, M
2
, M
3
.
(b) Steep slope curves S
1
, S
2
, S
3
.
(c) Critical slope curves C
1
, C
2
, C
3
.
(d) Horizontal slope curves H
2
, H
3
.
(e) Averse slope curves A
2
, A
3
.
In all these curves, the letter indicates the slope type and the subscript
indicates the zone. For example S
2
curve occurs in the zone 2 of the
steep slope.
62
Flow Profiles in Mild slope
Flow Profiles in Steep slope
63
Flow Profiles in Critical slope
Flow Profiles in Horizontal slope
Flow Profiles in Adverse slope
64
Dynamic Equation of Gradually Varied Flow
Objective: get the relationship between the water surface slope and other
characteristics of flow.
The following assumptions are made in the derivation of the equation
1. The flow is steady.
2. The streamlines are practically parallel (true when the variation in
depth along the direction of flow is very gradual). Thus the hydrostatic
distribution of pressure is assumed over the section.
3. The loss of head at any section, due to friction, is equal to that in the
corresponding uniform flow with the same depth and flow
characteristics. (Mannings formula may be used to calculate the slope
of the energy line)
4. The slope of the channel is small.
5. The channel is prismatic.
6. The velocity distribution across the section is fixed.
7. The roughness coefficient is constant in the reach.
.
65
H Z y
V
g
= + +
2
2
dH
dx
dZ
dx
dy
dx
d
dx
V
g
= + +
|
\
|
.
|
2
2
Consider the profile of a gradually varied flow in a small length dx of an open
channel the channel as shown in the figure below.
The total head (H) at any
section is given by:
Taking x-axis along the bed of the channel and differentiating the equation with
respect to x:
66
= + +
|
\
|
.
| S S
dy
dx
d
dx
V
g
f 0
2
2
dy
dx
dy
dy
d
dx
V
g
S S
f
+
|
\
|
.
| =
2
0
2
dy
dx
d
dy
V
g
S S
f
1
2
2
0
+
|
\
|
.
|
(
(
=
dy
dx
S S
d
dy
V
g
f
=
+
|
\
|
.
|
0
2
1
2
dH/dx = the slope of the energy line (S
f
).
dZ/dx = the bed slope (S
0
) .
Therefore,
Multiplying the velocity term by dy/dy and transposing, we get
or
This Equation is known as the dynamic equation of gradually varied flow. It
gives the variation of depth (y) with respect to the distance along the bottom of
the channel (x).
67
dy
dx
S S
Q T
g A
f
=
0
2
3
1
dy
dx
dE dx
Q T
g A
=
/
1
2
3
The dynamic equation can be expressed in terms of the discharge Q:
The dynamic equation also can be expressed in terms of the specific energy E :
68
dy
dx
= 0
dy
dx
positive =
dy
dx
negative =
Depending upon the type of flow, dy/dx may take the values:
The slope of the water surface is equal to the bottom
slope. (the water surface is parallel to the channel bed)
or the flow is uniform.
The slope of the water surface is less than the bottom slope
(S
0
) . (The water surface rises in the direction of flow) or the
profile obtained is called the backwater curve.
The slope of the water surface is greater than the bottom
slope. (The water surface falls in direction of flow) or the
profile obtained is called the draw-down curve.
(a)
(c)
(b)
69
Notice that the slope of water surface with respect to horizontal (S
w
) is different
from the slope of water surface with respect to the bottom of the channel (dy/dx).
A relationship between the two slopes can be obtained:
S
bc
ab
cd bd
ab
w
= = =
sin|
S
cd
ad
cd
ab
0
= = ~ sinu
Consider a small length dx of
the open channel.
The line ab shows the free
surface,
The line ad is drawn parallel
to the bottom at a slope of S
0
with the horizontal.
The line ac is horizontal.
Let u be the angle which the bottom makes with the horizontal. Thus
The water surface slope (Sw) is given by
70
dy
dx
bd
ad
bd
ab
= ~
dx
dy
S S
w
=
0
dy
dx
S S
w
=
0
The slope of the water surface with respect to the channel bottom is given by
This equation can be used to calculate the water
surface slope with respect to horizontal.
71
Water Profile Computations (Gradually Varied Flow)
Engineers often require to know the distance up to which a surface
profile of a gradually varied flow will extend.
To accomplish this we have to integrate the dynamic equation of
gradually varied flow, so to obtain the values of y at different locations
of x along the channel bed.
The figure below gives a sketch of calculating the M1 curve over a
given weir.
72
Direct Step Method
One of the most important method used to compute the water profiles is
the direct step method.
In this method, the channel is divided into short intervals and the
computation of surface profiles is carried out step by step from one section
to another.
For prismatic channels:
Consider a short length of channel, dx , as shown in the figure.
dx
73
S dx y
V
g
y
V
g
S dx
f 0 1
1
2
2
2
2
2 2
+ + = + +
S dx E E S dx
f 0 1 2
+ = +
dx
E E
S S
f
=
2 1
0
Applying Bernoullis equation between section 1 and 2 , we write:
or
or
where E
1
and E
2
are the specific energies at section 1 and,
respectively.
This equation will be used to compute the water profile curves.
74
1 2
y y >
1 2
y y <
V
n
R S
f 1 1
2 3
1
1
=
/
V
n
R S
f 2 2
2 3
2
1
=
/
S
S S
fm
f f
=
+
1 2
2
The following steps summarize the direct step method:
1. Calculate the specific energy at section where depth is known.
For example at section 1-1, find E
1
, where the depth is known (y
1
). This
section is usually a control section.
2. Assume an appropriate value of the depth y
2
at the other end of the small
reach.
Note that:
if the profile is a rising curve and,
3. Calculate the specific energy (E
2
) at section 2-2 for the assumed depth (y
2
).
4. Calculate the slope of the energy line (S
f
) at sections 1-1 and 2-2 using
Mannings formula
and
And the average slope in reach is calculated
if the profile is a falling curve.
75
L dx
E E
S S
fm
1
2 1
0
,2
= =
L
E E
S
S S
f f
1 2
2 1
0
1 2
2
,
=
+
|
\
|
.
|
n n
L L L L
, 1 3 , 2 2 , 1
.......
+ + + =
5. Compute the length of the curve between section 1-1 and 2-2
or
6. Now, we know the depth at section 2-2, assume the depth at the next
section, say 3-3. Then repeat the procedure to find the length L
2,3
.
7. Repeating the procedure, the total length of the curve may be obtained.
Thus
where (n-1) is the number of intervals into which the channel is divided.