Concept of Uniform Flow: Lecture Note For Open Channel Hydraulics
Concept of Uniform Flow: Lecture Note For Open Channel Hydraulics
CAPTER Four
UNIFORM FLOW
Concept of Uniform Flow
Uniform flow is referring the steady uniform flow
Steady flow is characterized by no changes in time.
Uniform flow is characterized by the water cross section and depth remaining constant
over a certain reach of the channel.
For any channel of given roughness, cross section and slope, there exists one and only one
water depth, called the normal depth , at which the flow will be uniform.
Uniform equilibrium flow can occur only in a straight channel with a constant channel
slope and cross-sectional shape, and a constant discharge.
The energy grade line Sf, water surface slope Sw and channel bed slope S0 are all parallel,
i.e. Sf=Sw=So
For steady uniform channel flow, channel slope, depth and velocity all remain constant
along the channel
V2/2g
P1
V Wsin
Y0
Y0 P2
W
Y0
A
Y0
L
A
ALS 0 0 PL 0 S 0 RS 0 ------------------ (4.2)
P
Where R= A/P = hydraulic radius , which is a length parameter accounting for the shape
of the channel. And it plays a very important role in developing flow equations which are
common to all shapes of channels.
Expressing the average shear stress 0 = kV2, where k=a coefficient which depends
on the nature of the surface and flow parameters. Equation (4.2) can be written as
1
Where C = a coefficient which depends on the nature of the surface
k
Equation (4.3) is known as Chezy formula and the coefficient C is known as the chezy
coefficient.
We remember that for pipe flow, the Darcy –Weisbach equation is
L V2
hf f
D 2g
Where hf = head loss due to friction in a pipe of diameter D and length L
f = Darcy-Weisbach friction factor
By Belete B. AAiT Department of Civil Engineering October 2010 Page 3 of 26
Lecture Note for Open Channel Hydraulics
VD
For smooth pipes, f is found to be a function of the Reynolds number ( Re ) only.
v
For rough turbulent flows, f is a function of the relative roughness (/D) and types of
roughness, which independent of the Reynolds number.
So for the case of Open Channel, we can be considered to be a conduit cut into two. The
hydraulics radius would then be appropriate length parameter and prediction of friction
factor f. So equation of head loss due to friction is written as
L V2 8g
hf f V . R . h f / L -------------------- (4.4)
4R 2g f
Note that for uniform flow in an open channel hf/L = slope of the energy line = Sf = S0, it
may be seen that equation (4.4) is the same as Chezy formula, equation (4.3). with
C 8g / f --------------------------------------------------------------(4.5)
Equation (4.5) be use to develop different charts and empirical formulas that relate C
with Re.
4 RV 4R
Re and
v s
1
1.80 log Re 1.5146 ----------------------------------------------- (4.6)
f
1 21.25
And 1.14 2.0 log s 0.9 -------------------------------------------- (4.7)
f 4 R Re
Equation (4.7) is valid for 5000 ≤ Re ≤ 108 and 10-6 < s/4R < 10-2
Generally, the open channels that are encountered in the field are very large in size and
also in the magnitude of roughness elements. Due to scarcity of reliable experimental or
field data on channels covering a wide range of parameters, values of s are not available
to the same degree of confidence as for pipe materials. However, the following table will
give the estimated values of s for some common open channel surfaces.
Example 4.1: A 2.0m wide rectangular channel carries water at 20oc at a depth of 0.5m. The
channel is laid on a slope of 0.0004. Find the hydrodynamic nature of the surface if
the channel is made of
A. Very smooth concrete
B. Rough concrete
C. Estimate the discharges in the channel in both case using chezy formula with
Dancy-Weisbach f.
Solution:
B=2.0m R= A/P = (0.5*2)/(2+2*0.5)=0.33m
Y=0.5m o= RSo = (9.81x103)x0.333x0.0004 = 1.308N/m2
S0 = 0.0004 v* = shear velocity = o/ = (1.308/103) = 0.03617m/sec
T= 20oc Kinematic viscosity () @ 20oc = 10-6m2/s
a). For a smooth concrete surface
from the Equivalent roughness table, s = 0.25mm = 0.00025m
s v*/ = (0.00025*0.03617)/(10-6) = 9.04 > 4 but less than 60 Transition
b) for a rough concrete
From the Equivalent roughness table, s = 3.5mm = 0.0035m
s v*/ = (0.0035x0.03617)/(10-6) =126.6 > 60 Rough
c). case (i) = smooth concrete channel
s= 0.25mm and s/4R = (0.25)/(4x0.33x103) = 1.894x 10-4
1 21.25
1.14 2.0 log s 0.9
f 4 R Re
Since Re is known we use the moody chart or iterative trial and error methods, we get
By Belete B. AAiT Department of Civil Engineering October 2010 Page 5 of 26
Lecture Note for Open Channel Hydraulics
f= 0.0145
C= (8g/f) = (8x9.81/0.0145) = 73.6
V = C(RS0) = 73.6x(0.333x0.0004) = 0.850m/s
Q= AV= (2x0.5)x0.850 = 0.85m3/se
Case (ii) Rough concrete Channel
s= 3.5mm and s/4R = (3.5)/(4x0.33x103) = 2.625x 10-3
f = 0.025
C = (8g/f) = (8x9.81/0.025) = 56.0
V = C(RS0) = 56x(0.333x0.0004) = 0.647m/s
Q= AV= (2x0.5)x0.647 = 0.647m3/se
1. Pavlovsik formula: C
R x
n
in which x 2 .5 n 0 .13 0 .75 R
n 0 .10 and
n= manning’s coefficient. This formula appears to be in use in Russia
2. Bazin’s formula C 87.0 , in which M = a coefficient dependent on the surface roughness
1 M / R
1 0 . 00155
3. Ganguillet and Kutter Formula 23
, in which n= manning’s
n So
C
0 . 00155 n
1 23
S o R
coefficient
, ---------------------(4.14)
Where n = Equivalent roughness.
From Equ. (4.14),
------------------------------------------------------- (4.15)
This equation gives a means of estimating the equivalent roughness of a channel having multiple
roughness types in its perimeters
Table 4.3. Equations for equivalent Roughness Coefficient (adopt from K. Subramanya 2010)
No Investigators ne Concept
1 Horton (1933); Einstein (1934) 2/3 Mean Velocity is
n
1
3/2
i Pi constant in all subareas
P
2 Pavloskii (1931), Muhlhofer(1933) 1/ 2 Total resistance force F
Einstein and Banks (1950)
1
n i2 Pi is sum of subarea
P
resistance force, Fi
3 Lotter (1932)
PR 5 / 3 Total discharge is sum
P i R i5 / 3
of subarea discharge
ni
4 Yen (1991)
n P i i
Total shear velocity is
P weighted sum of subarea
shear velocity
Example 4.2: An earthen trapezoidal channel (n = 0.025) has a bottom width of 5.0 m, side slopes
of 1.5 horizontal: 1 vertical and a uniform flow depth of 1.10 m. In an economic study
to remedy excessive seepage from the canal two proposals, a) to line the sides only
and, b) to line the bed only are considered. If the lining is of smooth concrete (n =
0.012), calculate the equivalent roughness in the above two cases.
1
Where K AR 2 / 3 is called the conveyance of the channel and express the discharge capacity of
n
the channel per unit longitudinal slope. The term nK = AR2/3 is also called the section factor for
uniform flow computations.
The basic variables in uniform flow problems can be the discharge Q, velocity of flow V, normal
depth y0, roughness coefficient n, channel slope S0 and the geometric elements (e.g. B and side
slope m for a trapezoidal channel). There can be many other derived variables accompanied by
corresponding relationships. From among the above, the following five types of basic problems are
recognized.
Table 4.5 Problem Types and the given and required variables in uniform flow
Hydraulic Radius
Hydraulic radius plays a prominent role in the equations of open-channel flow and therefore, the
variation of hydraulic radius with depth and width of the channel becomes an important
consideration. This is mainly a problem of section geometry.
Rectangular Channel
Area A = Byo
Wetted Perimeter P = B + 2yo
Y0
Hydraulic Radius R= A/P
y ---------- (4.18)
R o
B y o
1 2
B
a). wide Rectangular Channel
As yo/B, the aspect ratio of the channel decrease, R yo. Such channels with large bed-widths as
compared to their respective depths are known as wide rectangular channels. In these channels, the
hydraulics radius approximates to the depth of flow.
Considering a unit width of a wide rectangular channel, A= yo, R = yo and B = 1.0
3/5
1 5 / 3 1/ 2 qn
Discharge per unit width Q/B = q yo S o yo --------------------------- (4.19)
n S o
This approximation of a wide rectangular channel is found applicable to rectangular channels
yo/B < 0.02.
b). Rectangular Channels with yo/B ≥ 0.02
Qn
AR 2 / 3 , and 2/3
By o
5/3
yo / B
5/3
B8/ 3
For these channels AR
So B 2 y o 2/3
1 2 y o / B 2/3
Qn AR 2 / 3
8/3
o 5/3
o ------------------------------------------------------ (4.20)
So B8 / 3 B 1 2 o 2 / 3
y0
Where 0 , Tables of (o) Vs o will provide a non-dimensional graphical aid for general
B
Qn
application. Since , one can easily find yo/B from this table for any combination of Q,
So B8/ 3
n, So, and B in a rectangular channel.
1 Area = A = (B + myo)yo
Wetted Perimeter = P= (B +2yom2+1)
m
Hydraulic Radius
yo
A B my o y o
B
R
P
B 2 m 2 1y o
Qn B myo 5 / 3 y o5 / 3
AR 2/3
, Non-dimensionalising the variables,
So B 2 m 2 1yo 2/3
AR 2 / 3 Qn 1 m o 5 / 3 o5 / 3
( o , m) ----------------------------------------- (4.21)
B8/ 3 So B8/ 3 1 2 m 1y o
2
2/3
Where o = yo/B.
Equation (4.21) represent as curves or tables of vs o with m as the third parameter to provide a
general normal depth solution aid. It may be noted m=0 is the case of rectangular.
1 r
m r yo m=cot
r =yo r =yo
B
By Belete B. AAiT Department of Civil Engineering October 2010 Page 18 of 26
Lecture Note for Open Channel Hydraulics
Referring the above figure the full supply depth = normal depth at design discharge = yo at normal
depth
For standard trapezoidal case
1
Area = A = Byo +myo2 + yo2 = (B +yo)yo , where m m tan 1
m
Wetted perimeter =P = B +2myo +2yo = B +2yo
( B y0 ) y o
Hydraulics radius =R= A/P =
B 2 yo
1 ( B y o ) 5 / 3 y o5 / 3 1 / 2
By manning’s formula Q S o ---------------------------------------- (4.22)
n ( B 2 yo ) 2 / 3
Qn 5 / 3 (1 o ) 5 / 3 o5 / 3
( o ) ------------------------------------------------------ (4.23)
S o1 / 2 B 8 / 3 (1 2 o ) 2 / 3
For standard triangular case similarly we can drive
1
Q (y o2 )( y o / 2) 2 / 3 S o1 / 2 ------------------------------------------------------------------------- (4.24)
n
Example 4.3 A standard lined trapezoidal cannel section is to be designed to convey 100m3/sec of
flow. The side slope is to be 1.5H:1V and the manning’s coefficient n= 0.016. The longitudinal
slope of the bed is 1 in 5000m. If a bed width of 10m is preferred what would be the normal depth?
Solution = m+ tan-1(1/m)= 1.5 +tan-1(1/1.5) yo
Then 0 , which imply that
Q= 100m3/sec = 2.088 B
m= 1.5 Qn 5 / 3 100 X 0.016 X ( 2.088) 5 / 3 o B 0.74 X 10.0
1 0.8314 y0
n= 0.016 S o1 / 2 B 8 / 3 (0.0002 )1 / 2 (10 .0) 8 / 3 2.088
y 0 3.554
So = 0.0002
(1 o ) 5 / 3 o5 / 3
B= 10m 1 0 . 8314
(1 2 o ) 2 / 3
yo =?
By solving trial and error
o = 0.74
my my
a). For a given side slope m, what will be the water depth y for best hydraulic trapezoidal cross-
section?
For a given A , P = Pmin
If the depth of flow is confined to the deep channel only (y < h), calculation of discharge by using
Manning’s equation is very simple. However, when the flow spills over the flood plain (y > h), the
problem of discharge calculation is complicated as the calculation may give a smaller hydraulic
radius for the whole stream section and hence the discharge may be underestimated. The following
method of discharge estimation can be used. In this method, while calculating the wetted perimeter
for the sub-areas, the imaginary divisions (FJ and CK in the Figure) are considered as boundaries for
the deeper portion only and neglected completely in the calculation relating to the shallower portion.
1. The discharge is calculated as the sum of the partial discharges in the sub-areas; for e.g. units 1,
2 and 3 in Figuer
2. The discharge is also calculated by considering the whole section as one unit, (ABCDEFGH
area in Figure), say Qw.
3. The larger of the above discharges, Qp and Qw, is adopted as the discharge at the depth y.
Where A and R are in general, functions of the geometric elements of the canal. If the canal is of
trapezoidal cross-section,
-------------------------------------- (4.25)
Equ. (4.25) has six variables out of which one is a dependent variable and the rest five are
independent ones. Similarly, for other channel shapes, the number of variables depends upon the
channel geometry. In a channel design problem, the independent variables are known either explicitly
or implicitly, or as inequalities, mostly in terms of empirical relationships. The canal-design practice
given below is meant only for rigid-boundary channels, i.e. for lined and unlined non-erodible
channels.
Canal Section
Normally a trapezoidal section is adopted. Rectangular cross-sections are also used in special
situations, such as in rock cuts; steep chutes and in cross-drainage works. The side slope, expressed as
m horizontal: 1 vertical, depends on the type of canal, (i.e. lined or unlined, nature and type of soil
through which the canal is laid). The slopes are designed to withstand seepage forces under critical
conditions, such as;
1. A canal running full with banks saturated due to rainfall,
2. The sudden drawdown of canal supply.
Usually the slopes are steeper in cutting than in filling. For lined canals, the slopes roughly
correspond to the angle of repose of the natural soil and the values of m range from 1.0 to 1.5 and
rarely up to 2.0. The slopes recommended for unlined canals in cutting are given in Table (4.7).
In addition to the maximum velocities, a minimum velocity in the channel is also an important
constraint in the canal design. Too low velocity would cause deposition of suspended material, like
silt, which cannot only impair the carrying capacity but also increase the maintenance costs. Also, in
unlined canals, too low a velocity may encourage weed growth. The minimum velocity in irrigation
channels is of the order of 0.30 m/sec.
By Belete B. AAiT Department of Civil Engineering October 2010 Page 25 of 26
Lecture Note for Open Channel Hydraulics
Free Board
Free board for lined canals is the vertical distance between the full supply levels to the top of lining
(Fig. 4.17). For unlined canals, it is the vertical distance from the full supply level to the top of the
bank.
This distance should be sufficient to prevent overtopping of the canal lining or banks due to waves.
The amount of free board provided depends on the canal size, location, velocity and depth of flow.
If any other value of m is use, the corresponding value of B/y0 for the efficient section would be
In large channels it is necessary to limit the depth to avoid dangers of bank failure. Usually depths
higher than about 4.0 m are applied only when it is absolutely necessary. For selection of width and
depth, the usual procedure is to adopt a recommended value.
Example 4.4.