Chapter 1
Chapter 1
Chapter One
When the flow takes place in a channel or pipe such that the water has a free surface exposed to
the atmosphere, we spoke of open channels, culverts, spillways, and similar human made
structures are designed & analyzed by the method of open channel hydraulics.
The primary differences b/n the confined flow in pipes & open channel flow is that the pipe flow
is closed channel, which is the top surface is covered by solid boundary, it is not exposed to
atmospheric pressure but open channel flow is exposed to atmospheric pressure. In open
channels the cross-sectional area of the flow is variable that depends on many parameters of the
flow. For this reason hydraulic computations related to open channel flow are more complicated.
The prime motivating force (the force causing motion) for open channel flow is gravity or the
slope provided at the bottom (bed).
EL Hf
Hf
Y1 HGL V 2
EL
2g
HGL
Y1
Y2
Y2
Z1
Z2 Z1
Z2
Despite the similarity between these two flows it is much more difficult and complex to solve
problems of the open channel case. This is due to the fact that the flow condition in open channel
flow varies as per time and place. When we say the flow condition it includes depth of flow,
cross-sectional area and slope of the channel. In turn the depth of flow, discharge and slope of
the channel and water surface are related to each other.
In addition the bed roughness varies greatly leading the selection of friction coefficient to
uncertainty. The cause of flow in open channel the gravitational forces and viscous shear forces
along the channel wetted perimeter resists flow.
Flow cross section Known(fixed) Unknown in advance because the flow depth is unknown
Characteristics Velocity deduced Flow depth deduced simultaneously from solving both
flow parameters from continuity continuity & momentum equations
V p y
0 , 0 and 0
t t t
B. Unsteady flow: - Here the fluid characteristics vary with time such that
V p y
0 0 0
t , t and t
C. Uniform flow: - A space as a criterion is used. Open channel flow is said to be uniform if
the depth of flow, velocity remains constant or the same at every section of the channel.
Uniform flow may be steady or unsteady, depending on whether or not the depth changes
with time.
V y
0 0
s , and s
D. Non uniform flow: - In case when the velocity, depth of flow in a channel changes with
space:
V y
0 0
s , and s
E. Steady uniform flow: - The depth of flow does not change during time interval and space
under consideration.
F. Steady non-uniform flow:- Depth varies with distance but not with time. This type of flow
may be either (a) gradually varied or (b) rapidly varied. Type (a) requires the application of the
energy and frictional resistance equations while type (b) requires the energy and momentum
equations.
G. Unsteady uniform flow: - This is a flow in which the depth is varying time but not with
space.
H. Unsteady non uniform flow: - Is the flow in which the depth is varying with space and time.
Natural channels: These channels naturally exist without the influence of human beings.
E.g. Rivers, streams, tidal estuaries, aqueducts.
Aqueducts are underground conduits which carry water with free surface.
Artificial channels: Such channels are formed by man’s activity for various purposes.
E.g. irrigation channel, navigation channel, sewerage channel, culverts, power canal……
etc.
The above two channels can have either of the following features:
Generally the natural channels fall into the non-prismatic group. That is why intensive study of
the behavior of flow in natural channels requires other fields of studies like, sediment transport,
geomorphology, hydrology, river engineering.
The commonly needed geometric properties are shown in the figure below and defined as:
Depth (y) – the vertical distance from the lowest point of the channel section to the free surface.
Stage (z) – the vertical distance from the free surface to an arbitrary datum
Area (A) – the cross-sectional area of flow, normal to the direction of flow
Wetted perimeter (P) – the length of the wetted surface measured normal to the direction of
flow
Surface width (B) – width of the channel section at the free surface
Hydraulic radius (R) – the ratio of area to wetted perimeter (A/P)
Hydraulic mean depth (D) – the ratio of area to surface width (A/B)
The uniform flow formulas given above show that for a given slope and roughness, the velocity
increases with the hydraulic radius. Therefore, for a given area of water cross-section, the rate of
discharge will be a maximum when R is a maximum, which is to say, when the wetted perimeter
and so the frictional resistance is a minimum. Such a section is called the most efficient cross-
section.
Equation (b) indicates that for the given area of cross section and slope the discharge Q will be
maximum when R is maximum.
We can conclude that for most efficient and economical channel section the wetted perimeter
should be minimum & also frictional resistance, o is minimum.
A = BY…………………………. (i)
P = B + 2Y………………………. (ii)
dp d
0 ( A / Y 2Y ) 0
dY dY
A 20
Y
A 2Y 2 B * Y
Thus, the rectangular channel is most efficient and economical when the depth of water is one-
half of the width of the channel and the discharge will be maximum.
Accordingly, the most efficient channel shape is the semi-circle. The usual shape for new
channel & canal is the rectangular or trapezoidal such that the inscribed semi-circle is tangential
to the bed & sides.
Although first developed for solid bodies they are equally applicable to fluids. A brief
description of the concepts are given below.
Conservation of matter
This says that matter cannot be created nor destroyed, but it may be converted (e.g. by a
chemical process.) In fluid mechanics we do not consider chemical activity so the law reduces to
one of conservation of mass.
Conservation of energy
This says that energy cannot be created nor destroyed, but may be converted form one type to
another (e.g. potential may be converted to kinetic energy). When engineers talk about energy
"losses" they are referring to energy converted from mechanical (potential or kinetic) to some
other form such as heat. A friction loss, for example, is a conversion of mechanical energy to
heat. The basic equations can be obtained from the First Law of Thermodynamics but a
simplified derivation will be given below.
Conservation of momentum
The law of conservation of momentum says that a moving body cannot gain or lose momentum
unless acted upon by an external force. This is a statement of Newton's Second Law of Motion:
Force = rate of change of momentum
In solid mechanics these laws may be applied to an object which is has a fixed shape and is
clearly defined. In fluid mechanics the object is not clearly defined and as it may change shape
constantly. To get over this we use the idea of control volumes. These are imaginary volumes of
fluid within the body of the fluid. To derive the basic equation the above conservation laws are
applied by considering the forces applied to the edges of a control volume within the fluid.
1.4.1 The Continuity Equation (conservation of mass)
For any control volume during the small time interval t the principle of conservation of mass
implies that the mass of flow entering the control volume minus the mass of flow leaving the
control volume equals the change of mass within the control volume.
If the flow is steady and the fluid incompressible the mass entering is equal to the mass leaving,
so there is no change of mass within the control volume.
Considering the control volume above which is a short length of open channel of arbitrary cross
section then, if is the fluid density and Q is the volume flow rate then mass flow rate is
Q and the continuity equation for steady incompressible flow can be written
Q entering = Q leaving
As, Q, the volume flow rate is the product of the area and the mean velocity then at the upstream
face (face 1) where the mean velocity is u1 and the cross-sectional area is A
Q entering = A1U1
Similarly at the downstream face, face 2, where mean velocity is U2 and the cross-sectional area
is A2
Then: Q leaving = A2U2
Consider the forms of energy available for the above control volume. If the fluid moves from the
upstream face 1, to the downstream face 2 in time t over the length L.
The work done in moving the fluid through face 1 during this time is
Work done = P1A1L
Where
p is pressure at face 1
The mass entering through face 1 is
Mass entering = 1A1L
If z1 is the height of the centroid of face 1, then the potential energy of the fluid entering the
control volume is:
PE = mgz1 = 1A1 Lg z1
The total energy entering the control volume is the sum of the work done, the potential and the
kinetic energy:
We can write this in terms of energy per unit weight. As the weight of water entering the control
volume is 1 A1L g then just divide by this to get the total energy per unit weight:
2
p1 u1
Total energy per unit weight = Z1
g 2 g
At the exit to the control volume, face 2, similar considerations deduce
2
p u
Total energy per unit weight = 2 2 Z 2
g 2 g
If no energy is supplied to the control volume from between the inlet and the outlet then energy
leaving = energy entering and if the fluid is incompressible
1 = 2 =
2 2
p1 u1 p u
So Z 1 = 2 2 Z 2 = H= constant---------------------------------eqn. 1.2
g 2 g g 2 g
This is the Bernoulli equation.
Note:
1. In the derivation of the Bernoulli equation it was assumed that no energy is lost in the control
volume - i.e. the fluid is frictionless. To apply to non frictionless situations some energy loss
term must be included
2. The dimensions of each term in equation 1.2 has the dimensions of length (units of meters).
For this reason each term is often regarded as a "head" and given the names
3. Although above we derived the Bernoulli equation between two sections it should strictly
speaking be applied along a stream line as the velocity will differ from the top to the bottom of
the section. However in engineering practice it is possible to apply the Bernoulli equation
without reference to the particular streamline
In deriving the above momentum and energy (Bernoulli) equations it was noted that the velocity
must be constant (equal to V) over the whole cross-section or constant along a stream-line.
Clearly this will not occur in practice. Fortunately both these equation may still be used even for
situations of quite non-uniform velocity distribution over a section. This is possible by the
introduction of coefficients of energy and momentum, and β respectively.
The values of and β must be derived from the velocity distributions across a cross-section.
They will always be greater than 1, but only by a small amount consequently they can often be
confidently omitted – but not always and their existence should always be remembered. For
turbulent flow in regular channel does not usually go above 1.15 and will normally be below
1.05. We will see an example below where their inclusion is necessary to obtain accurate results.
To determine the values of and the velocity distribution must have been measured (or be
known in some way). In irregular channels where the flow may be divided into distinct regions
may exceed 2 and should be included in the Bernoulli equation.
The figure below is a typical example of this situation. The channel may be of this shape when a
river is in flood – this is known as a compound channel.
When uniform flow occurs gravitational forces exactly balance the frictional resistance forces
which apply as a shear force along the boundary (channel bed and walls).
Uniform flow is the result of exact balance between the gravity and friction force
A L sin = o .P.L
Solving for o ,
A
o = .S R.S ………………………………… eqn. 1.4
P
The shear stress is assumed proportional to the square of the mean velocity,
Or o= kV2
V2 = RS , Let C 2 -constant (b/c &k- are constant)
k k
A very many studies have been made of the evaluation of C for different natural and manmade
channels. These have resulted in today most practicing engineers use some form of this
relationship to give C:
1
R6
C
n
This is known as Manning’s formula, and the n as Manning’s n.
Substituting equation 1.5 in to the above equation gives velocity of uniform flow:
1 2 1
V= R 3 S0 2 ……………………………………………… eqn. 1.6
n
1 2 1
Q= AR 3 S 0 2
n
1.5.3 Conveyance
Channel conveyance, K, is a measure of the carrying capacity of a channel. The K is really an
agglomeration of several terms in the Chezy or Manning's equation:
So,
Use of conveyance may be made when calculating discharge and stage in compound channels
and also calculating the energy and momentum coefficients in this situation.
Specific Energy
It is total available energy in a given open channel flow, taking the bed of the flow channel as
datum line.
For any cross section, shape, the specific energy (E) at a particular section is defined as the
energy head to the channel bed as datum. Thus,
V2
E Y …………………………………………….. (1)
2g
For a rectangular channel, the value of flow per unit width is Q/B = q, and average velocity
qB q
V Q
A BY Y
2
q y
E y y q
2
…………………………………… (2)
2g 2 gy 2
q2
( E y)Y
2
(For the case of constant q)………………………… (3)
2g
A plot of E vs y is a hyperbola like with asymptotes (E-y) = 0 i.e. E = y and y = 0. Such a curve
is known as specific energy diagram.
The two alternative depths represent two totally different flow regimes slow & deep on the upper
limp of the curve (sub critical flow) & fast and shallow on the lower limb of the curve.(super
critical flow)
Flow over a raised hump -Application of the Specific energy equation.
The specific energy equation may be used to solve the raised hump problem. The figure below
shows the hump and stage drawn alongside a graph of Specific energy E against y.
Apply the Bernoulli equation between sections 1 and 2. (Assume a horizontal rectangular
channel
These points are marked on the figure. Point A on the curve corresponds to the specific energy at
Point 1 in the channel, but Point B or Point B' on the graph may correspond to the specific
energy at point 2 in the channel.
All point in the channel between point 1 and 2 must lie on the specific energy curve between
Point A and B or B'. To reach point B' then this implies that Es1 –Es2 >Z which is not physically
possible. So point B on the curve corresponds to the specific energy and the flow depth at section
2.
From these graphs we can identify several important features of rapidly varied flow.
Q2
E y ………………………………………….. (1) [V=Q/A]
2gA2
dE Q 2 dA
1 3 ……………………………………….. (2)
dy gA dy
2
Qmax Tc
3
1 …………………………………………….. (3)
gAc
3
gAc
From eqn. (3) Q 2 and substitute in eqn. (11) then
Tc
Ac
Ec y c ………………………………………….. (4)
2Tc
Q 2T
eqn. (3) can be solved by trial & error for irregular section by plotting f ( y ) and critical
gA3
depth occurs for the value of y which makes f(y) =1
Its physical significance is the ratio of inertial forces to gravitational forces squared
Its value determines the regime of flow – sub, super or critical, and the direction in which
disturbances travel
A transition in general form may have a change of channel shape, provision of a hump or a
depression, contraction or expansion of channel width, or any combination.
The concepts of specific energy and critical energy are useful in the analysis of open channel
transition problems.
Bed raising===hump
Width narrowing===constriction
Width widening===expansion
b/ Bed raising
i.e. E2=minimum
c/ Width narrowing
The specific energy @1 is equal to the specific energy @2, since no loss in energy.
Choking condition:-
DZ<=DZmax B2<=Bmini
The upstream water surface elevation is not affected by the conditions at
section 2 till a critical stage is first achieved.
DZ>DZmax B2>Bmini
the upstream water depth is different from Y1 and this condition is called
choking
The concept of this principle is applied when it’s not possible to apply the energy principle. One of
the major situations where momentum principle is applied is during analysis of rapidly varied flow
(RVF). In rapidly varied flow there is an excess amount of energy loss. Therefore, it is not possible
to analyze such type of flow using energy equation. Good example for RVF is Hydraulic Jump.
Hydraulic Jump
By far the most important of the local non-uniform flow phenomena is that which occurs when
supercritical flow has its velocity reduced to sub critical. There is sudden rise in water level at
the point where hydraulic jump occurs (Rapidly varied flow). This is an excellent example of the
jump serving a useful purpose, for it dissipates much of the destructive energy of the high –
velocity water, thereby reducing downstream erosion. The turbulence with in hydraulic jumps
has also been found to be very useful & effective for mixing fluids, & jumps have been used for
this purpose in water treatment plant & sewage treatment plants.
Y2
V2
V1
Y1
Lj
The momentum of water passing through section (1) per unit time is given as:
p1 rQV1
QV1 ………………………………………. (i)
t g
p2 rQV2
QV2 …………………………………………. (ii)
t g
P
Q (V2 V1 ) ………………………………………. (iii)
t
Q
A1Y 1 A2 Y 2 (V2 V1 ) …………………………………… (v)
g
Q Q Q
A1Y 1 A2 Y 2
g A2 A1
2
Q 1
A1Y 1 A2 Y2 A 1 A .......... .......... .......... .......... .........( iv)
g 2 1
Q2 Q2
A1Y 1 A2 Y 2 = Constant. …………… (vii)
gA1 gA2
M1 M2
M1and M2 are the specific forces at section (1) & (2) indicates that these forces are equal
before & after the jump.
Y2 = sequent depth
A2= By2
Y1 Y
Y 1 ,Y 2 2
2 2
Q2 y
By2 * 2
Q y
By1 1
gBy1
2 gBy2 2
Q2 By 2 Q2 By 2
1 2 .......... .......... .......... .......... .......... .......... ..(viii)
gBy1 2 Bgy2 2
Flow per unit width of q= Q/B Q=qB, then eqn. (viii) becomes
q 2 B 2 By12 q 2 B 2 By22
Bgy1 2 Bgy2 2
q2 1 1 y22 y12
………………………………… (.ix)
g 1
y y 2 2
2q 2
y1 y2
y y12
2
2
g y2 y1
2q 2
y1 y2 ( y1 y2 )......... .......... .......... .......... .......... .......... ........( x)
g
2q 2
y2 y12 y1 y22 0.......... .......... .......... .......... .......... .........( xi)
g
y2
2
y 2q
2
y1 2 .......... .......... .......... .......... .......... .( xii)(a )
2 2 gy2
y1
2
y1 2q
2
y2 .......... .......... .......... .......... ........( b)
2 2 gy2
y2 8q 2
y1 (1 1 )......... .......... .......... .......... .......... ..( c)
2 gy23
y1 8q 2
y2 (1 1 .......... .......... .......... .......... .......... ...( xii)(d )
2 gy13
2
y1 8q
1 (1 1 3 .......... .......... .......... .......... .......( xii)(e)
y2 2 gy 2
2
y2 8q
1 (1 1 3 .......... .......... .......... .......... .......... ( f )
y1 2 gy 1
q
V1 V2 y2 q
F1 , F2
gy1 gy2 gy2 gy23
y1 1
Therefore (1 1 8F22 )......... .......... .....( g )
y2 2
y2 1
(1 1 8F12 .......... .......... .......... .......... ..( h)
y1 2
q2 q2
y1 y .......... .......... .......... ......( b)
2 gy12 2 gy22
2
q 2 y22 y12
( y2 y1 )......... .......... ....( c)
2 g y12 y22
2q 2
From eqn. (x) substituting: y1 y2 ( y1 y2 ) in to this eqn. & by rearranging:
g
hlf E
y2 y1 3 .......... .......... .......... .......... ......( 2)
4 y1 y2
Therefore power lost = Q hlf (kw)………………… (3)
Types of Hydraulic jump
Hydraulic jumps are classified according to the upstream Froude number and depth ratio
F1 Y2/y1 Classification