Department of Civil Engineering: Term Paper
Department of Civil Engineering: Term Paper
ROLL: 6413
COURSE NO : CE 801
SUBMITTED BY : SUBMITTED TO :
MD. TANVIR ANJUM ENGR.RABINDRA
RANJAN SAHA
ID NO. :WUB 10/18/75/6413
PEng. Associate Professor and Head
BATCH :75/A ;ROLL :6413 Department of Civil Engineering
TERM PAPER
OPEN CHANNEL FLOW
CE801
1.0. INTRODUCTION:
Open-channel flow, a branch of hydraulics and fluid mechanics, is a type of
liquid flow within a conduit or in channel with a free surface, known as a
channel. The other type of flow within a conduit is pipe flow. These two types
of flow are similar in many ways but differ in one important respect: the free
surface. Open-channel flow has a free surface, whereas pipe flow does not.
1.1. OBJECTIVES:
The objective of this course is to introduce Open Channel Flow to students,
explaining the types of open channel and their behaviours, the causes and
principles of such behaviours, and applications open channels, enabling the
students to identify the open channels, and to analyse, design and manage
some of the types.
2.2. Application:
Common topics of design for hydraulic engineers include hydraulic
structures such as dams, levees, water distribution networks, water
collection networks, sewage collection networks, storm water management,
sediment transport, and various other topics related to transportation
engineering and geotechnical engineering. Equations developed from the
principles of fluid dynamics and fluid mechanics are widely utilized by other
engineering disciplines such as mechanical, aeronautical and even traffic
engineers.
Related branches include hydrology and rheology while related applications
include hydraulic modeling, flood mapping, catchment flood management
plans, shoreline management plans, estuarine strategies, coastal protection,
and flood alleviation.
2.4. APPLICATION:
Page 3 of 29
TANVIR ANJUM
ROLL: 6413
Rivers.
Streams.
Canals.
water flowing through a culvert running underneath a street.
flows in sewers.
tunnels water treatment plants.
storm sanitary sewer systems.
industrial waste applications.
sewage treatment plants.
and irrigation systems.
The depth of flow does not change over time, or if it can be assumed to be
constant during the time interval under consideration.
Page 4 of 29
TANVIR ANJUM
ROLL: 6413
Rapidly-varied flow
Gradually-varied flow
o Continuous flow
o Spatially-varied flow
Fr =
Page 5 of 29
TANVIR ANJUM
ROLL: 6413
where U is the mean velocity, D is the characteristic length scale for a
channel's depth, and g is the gravitational acceleration. Depending on the
effect of viscosity relative to inertia, as represented by the Reynolds number,
the flow can be either laminar, turbulent, or transitional. However, it is
generally acceptable to assume that the Reynolds number is sufficiently
large so that viscous forces may be neglected.
The expanded Bernoulli equation the mechanical energy per unit volume of
fluid moving along a streamline, v2/2 +p + g h is constant. This can be
written a little more conveniently for our purposes as energy per unit weight
of fluid Ew. Because weight equals volume multiplied by ρg
E = v2/2g + / +h
2.8. CRITICAL FLOW WITH CRITERIA:
For a given specific energy and discharge per unit width q, there are two
possible (real) depths of flow, and that transition from one depth to the other
can be accomplished under certain situations. These two depths
represented on the two different limbs of the E-y curve separated by the
crest c, are characteristic of two different kinds of flow; a rational way to
understand the nature of the difference between them is to consider first the
flow represented by the point c. Here the flow is in a critical condition,
Page 6 of 29
TANVIR ANJUM
ROLL: 6413
poised between two alternative flow regimes, and indeed the word “critical "
is used to describe this state of flow; it may be defined as the state at which
the specific energy E is a minimum for a given q.
(i) The specific energy and specific force are minimum for the given
discharge.
(iii) For a given specific energy the discharge is maximum at the critical
flow.
(iv) The velocity head is equal to half the hydraulic depth in a channel of
small slope.
(v) The velocity of flow in a channel of small slope with uniform velocity
distribution, is equal to the celerity of small gravity waves (C gh = ) is
shallow water caused by local disturbance.
Critical flow may occur at a particular section or in the entire channel, then
the flow in the channel is called "Critical flow". y f A,D c = ( ) for a given
discharge.
For a prismatic channel for a given discharge the critical depth is constant
at all sections of a channel. The bed slope which sustains a given discharge
at a uniform and critical depth is called "Critical slope Sc". A channel slope
causing slower flow in sub critical state for a given discharge is called "sub
critical slope or mild slope". A slope greater than the critical slope is called
steep slope or super critical slope.
Page 7 of 29
TANVIR ANJUM
ROLL: 6413
In fluid dynamics, the Chézy formula describes the mean flow
velocity of turbulent open channel flow. The formula is
V= C (RS)
Where V=mean velocity
R=hydraulic radius
S= slope of the energy line or bed slope and
C= a factor called Chezy’s
In 1889 the Irish Engineer’ Robert Manning presented a formula for design
of channel in FPS which is latter modified as well known form:
V=(1/n)(R2/3 S1/2)in SI unit
Where, v=mean velocity in mps
R=Hydraulic radius in meter
V=(1.49/n)(R2/3 S1/2) in FPS
As per Chezy,V=C
Where,n=Manning’s n
C=[1/n{(R2/3 S1/2)]/
C=[1/n]R(-)1/
Page 8 of 29
TANVIR ANJUM
ROLL: 6413
2.11.NORMAL DEPTH, NORMAL VELOCITY AND NORMAL
DISCHARGE:
Manning’s formula in FPS unit [A R2/3= nQ/(1.49 S)] is used for computation
of normal depth & normal velocity by using algebraic method.
2.12.Example :5-3
A channel having a bottom width 20 m, side slopes 2:1 , a depth of water 6
m and channel slope, S = 0.005 but a discharge of 3400 m3/s is observed
under given conditions (Figure below). Compute the values of Kutter’s n and
Manning’s n .
Page 9 of 29
TANVIR ANJUM
ROLL: 6413
Y=6
m
B=20
Solution :
Given,
b = 20 m ,
y=6m ;
side slope, z = 2
discharge, Q = 3400 m3/s
S = 0.005
Page 10 of 29
TANVIR ANJUM
ROLL: 6413
C = [ 23 + 0.00155/S + 1/ n] / [1+ {23 + 0.00155 / S}× {n/ √ R)}]
Again , C = V/√(RS) = 17.7/(√(4.10 × 0.005) = 123.62,
say C = 124
Putting the respective values of S, R and C in the above equation
= [ 23 + 0.00155/0.005 + 1/ n] / [1+ {23 + 0.00155 / 0.005}× {n/ √ 4.1)}]
124 = ( 23 + 0.31 + 1/n)/ {1+(23 + 0.31)(0.49n)
124 = ( 23.31+ 1/n) / (1 + 11.42 n
124 (1+ 11.42n) = (23.31n +1)/n
124 (n + 11.42n2) = 23.31 n+ 1
100.69 n + 1416.08 n2 - 1 = 0
1416.08 n2 +100.69 n - 1 = 0
The equation is a quadratic equation and has two values of n), like
ax2 + bx +c =0
x = [- b ± √ ( b 2– 4ac) ]/2a
n = [- 100.69 ± √ {100.692 – 4×1416.08 × (- 1)}]/
[2 x 1416.08]
n1 = [- 100.69 ± 125.71] /[ 2 x 1416.08]
n1 = [ -100.69 + 125.71]/[2832.16] = 0.009
n2 = [ - 100.69 – 125.71 ]/ [2832.16] = - 0.08
The value of n2 can not be negative
Hence the value of Kutter’s n = 0.009
2.13.Example 5-6
A trapezoidal channel has base width b = 20 m and side slopes 1 H :1 V.
The channel Manning’s roughness coefficient n = 0.014 and channel bottom
slope S = 0.0002 . Compute
Solution
Given, n = 0.014 S = 0.0002, and Q = 430 m3/s and side slopes; 1: z = 1:1
yn
B=20
Q = (1/n) A R⅔ √S
Page 12 of 29
TANVIR ANJUM
ROLL: 6413
A R⅔ = nQ/√S
= 0.014 ×430 /(√0.0002)
= 425.68
A = (b+zy)y = (20+y)y
P = b +2y√(1+z2) = 20 + 2y√2
Now LHS = A R⅔
= (b+zy)y [{(b+zy)y} / {b +2y√(1+z2)}]2/3
Using algebraic method, by trial and error
y A P R=A/P(m) R⅔ LHS= RHS Remarks
(m) (m2) (m) A R⅔
5.0 125.00 34.1 3.66 2.37 296.25 425.68 Far
6.15 160.82 37.4 4.30 2.66 427.78 y -
closest
2.14.Example 5- 4(A):
Page 13 of 29
TANVIR ANJUM
ROLL: 6413
A trapezoidal channel with b = 20 ft, z = 2, S0 = 0.0016 and n = 0.025,
carries a discharge of 400 cfs. Compute the normal depth and velocity and
critical depth. Hence state the flow properties.
Solution
To be calculated (a)normal depth (b) normal velocity and (c)critical depth (d)
state of flow
yn
B=20
Given
A R⅔ = n Q /(1.49√S) .............................Eq-1
From Eq-1
RHS = n Q /(1.49√S)
LHS = A R⅔
Page 14 of 29
TANVIR ANJUM
ROLL: 6413
Using Algebraic Method: Trial and Error method:
LHS =A R⅔
LHS =A R⅔
LHS =A R⅔
= [(20 + 2*3.3) *3.3] * [{(20+ 2* 3.3) 3.3} /{20+ 2 * 3.3 √( 1 +22 )}] ⅔
=162.6
LHS =A R⅔
Page 15 of 29
TANVIR ANJUM
ROLL: 6413
Q = 400 cfs (given) , We know,
Normal velocity, Vn = Qn /A
Q=AV
1. Normal discharge Q
2. Mean velocity of flow V
3. Normal depth y
4. Coefficient of roughness, n
5. Channel slope, S
6. Geometric elements that depend on the shape of the channel section,
such A and R When any four of the above variables are
given,remaining two can be computed by using the above two
equations are as cited below:
Page 16 of 29
TANVIR ANJUM
ROLL: 6413
1 To compute the normal discharge: In practical applications, the
capacity of a given channel or form the construction of a synthetic
rating curve of the channel
2 To determine the velocity of flow. It is often required the study of
scouring and silting effects in a given channel.
3 . To compute the confident of roughness: - The computation is used
to ascertain the roughness coefficient in a given channel, the
coefficient thus determined may be used in another similar channel.
4. To compute the channel slope – This is required to adjust the slope
of the given channel
5. To determine the dimensions of the channel This computation is
mainly required for design purposes.
Page 17 of 29
TANVIR ANJUM
ROLL: 6413
V2/2g = Q2 / 2gA y2
dA = T (y) dy
[ ]=- =
=- =-
dy/dx =
2.16.Example 7-1
A trapezoidal channel having b=20 m, side slope, 2:1, S0 = 0.0016 and n =
0.025 carries a discharge of 400 cumecs. Compute the backwater profile
created by a dam which backs up the water to a depth of 5 m immediately
behind the dam. The upstream end of the profile is assumed at a depth
equal to 1% greater than the normal depth (4.20 m). The energy coefficient
⍺ = 1.10
Solution:
Page 18 of 29
TANVIR ANJUM
ROLL: 6413
n = 0.025
Q = 400 cumecs
⍺ = 1.10
Z
B=20
Kn = Q / √S0
Z = √(A3 / T)
Zc = Q/ √(g/⍺)
Area A = (b + zy)
Putting the respective value in the above two equations for D and A, we get
Page 19 of 29
TANVIR ANJUM
ROLL: 6413
D = y(10+y)/ (10 + 2y) and A = y(20 + 2y)
V2 = Q2 / A2 = [400 / ( 20+2y)y] 2
V2 /2g = D /2
we get
K = 1×A R⅔ /n
Z = √(A3 / T)
R⅔ = 3.54⅔ = 2.323
dx/dy = 1008
Cross section 1
dy = 5 - 4.9 = 0.10m
Page 21 of 29
TANVIR ANJUM
ROLL: 6413
Average of dx/dy between the depth 5m and 4.9 m
y = 4.80 m
Average of
dx = 112 m
Similarly,
Profile draw
Length of slope =
= 927/625.0008 = 1.48 m
Page 22 of 29
TANVIR ANJUM
ROLL: 6413
The rapid changes in depth of flow from high stage to low stage causes
steep depression in the water surface and the phenomenon that causes an
abrupt change in the channel slope or cross section is known as hydraulic
drop.
2. Hydraulic jump
The local phenomenon at which the rapid change in the depth of flow is from
a low stage to a high stage there is abrupt rise of water surface is called the
hydraulic jump
Initial depth: The depth of flow in The special case of the hydraulic drop
that occurs where the bottom of flat channel is discontinued. As the free over
fall enters the air in the form of a nappe, there will be no reverse curve in the
water surface until it strikes some object at a lower elevation.
channel before jump is called the initial depth (y1).
Sequent depth: The depth of flow in channel after the jump is called the
sequent depth(y2 )
Practical application of hydraulic jump
1. To dissipate energy in water flowing over dams , weirs, and other
hydraulic structures.
Page 24 of 29
TANVIR ANJUM
ROLL: 6413
To indicate special flow conditions, such as the existence of
supercritical flow or the presence of a control section so that a gauging
station may be located
To mix chemicals used for water purification
To aerate water for city water supplies and
To remove air pocket from water-supply lines and thus prevent air
locking.
Hydraulic jump in rives are often used as recreation as for examples -
fun/sport by Kayaker (Kayak is a kind of light boat –Kayaker means
kayak race) and canoers i.e. who travels in canoe (a small, light,
narrow boat, pointed at both ends and moved using a paddle)
Two types
1.Undular jump
2.Direct jump
Undular jump:
If the jump is low, that is the change in depth is small, the water will not
rise obviously and abruptly but will pass from the low stage to the high
stage through a series of undulations gradually diminishing in size.
This low jump is called undular jump.
Direct jump:
When the jump is high that is, the change in depth is great, the jump is
called a direct jump. The energy content in the flow after the jump is
appreciably less than that before the jump.
Page 25 of 29
TANVIR ANJUM
ROLL: 6413
(2) F1 = 1 to 1.7, the water surface shows undulations and the jump is called
an undular jump
(3) F1 = 1.7 to 2.5 a series of small rollers develop on the surface of the jump
but downstream water surface remains smooth. The velocity throughout is
fairly uniform and the energy loss is low. This jump may be called weak jump
(4) F1 = 2.5 to 4.5, an oscillating jet entering the jump bottom to surface to
and back again with no periodically. Each oscillation produces a large wave
of irregular period which very commonly in canals, can travel for miles doing
unlimited damage to earth banks and ripraps. This jump may be called an
oscillating jump
(5) For F1 = 4.5 to 9.0 , the downstream extremity of the surface roller and
the point at which the high-velocity jet tends to leave the flow occur at
practically the same vertical section. The action and position of this jump are
least sensitive to variation in tail water depth. This jump is well-balanced and
the performance is at its best. The energy dissipation ranges from 45 to 70
%. This jump is called a steady jump
(6) F1 ≥ 9.0, the high velocity jet grabs of water rolling down the front face of
the jump, generating waves downstream and a rough surface can prevail.
The jump action is rough but effective since the energy dissipation may
reach 85%. This jump may be called a strong jump
2.18.Example 12-1
Page 26 of 29
TANVIR ANJUM
ROLL: 6413
Water flows in a rectangular channel that contains a sluice gate. The
upstream depth of flow is 1.50 m, the sluice gate opening is 0.62 m and the
unit discharge is 1.0 m2/s. Determine the flow depth downstream of the
sluice gate and draw specific energy curve.
Solution
Given; The upstream depth, y1 = 1.50 m ; The sluice gate opening , = 0.62
m
E = y + (q2/2y2 g )
Page 27 of 29
TANVIR ANJUM
ROLL: 6413
E2 = y2 + q2/(2y22 g ) = 0. 20 + 12/ (2×0.202 g ) = 1.52 m
2.19.Example 12-3
Channel section rectangular, initial depth (y1) 1.20 m. Calculate height of
hydraulic jump and length of jump. Also determine sequent depth. Velocity is
9 m/s.
Solution
or
Page 28 of 29
TANVIR ANJUM
ROLL: 6413
Length of jump
3.0. CONCLUSION’S:
After completing this course we have learn about open channel flow and its
classification,Velocity and pressure distributions, energy equations and
specific energy,Critical flow and control, principles of flow measurement and
devices,Concepts of uniform flow, Chezy's and Manning's equations and
their application,Computation of uniform flow, momentum equation and
specific momentum,Hydraulic jump and its,mechanism,Theory and analysis
of gradually varied flow,Computation of flow profile,Design of channel and
hydraulic structures. Now in real life we make full use of this course. After
the completion of this course we will be able
1.to understand flow of water and different flow patterns, flow mechanism
and channel geometry
2. to develop the open channel flow equations from the basic conservation
equation
3. to solve open channel flow problems through the selection of appropriate
equations and apply the physical mechanism of hydraulic jump
4 .to determine gradually varied flow and design channel with hydraulic
structures.
Page 29 of 29