Open Channel Flow
Open Channel Flow
Open Channel Flow
CHAPTER I
Introduction
Open channel flow occurs where ever the flow proceeds with the
liquid surface exposed to constant pressure.
This is the most complex flow type, requiring the solution of energy, momentum
and friction equations through time.
Types of Open Channel
Flume: a channel built (or supported) above the ground to convey fluid
from one point to another. In the field flumes are made of concrete,
wood, sheet metal or masonry. Laboratory flumes are usually made of
wood, metal, glass or a composite of these materials.
WHY?????
velocity distribution in open channels
Determination of energy and
momentum coefficients
To determine the values of α and β the
velocity distribution must have been
measured (or be known in some way).
dA
V13 A1 V23 A2 V33 A3 Vi Ai
3
u 3
3 3
V A 3
V ( A1 A2 A3 ) V Ai
Q V1 A1 V2 A2 V3 A3
V
V Ai i
A A1 A2 A3 A i
2
2
V A 2
V ( A1 A2 A3 ) V Ai
Uniform flow and the Development of
Friction formulae
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).
0.00281 1.811
41.65
C S n
0.00281 n
1 41.65
S R
where n is known as Kutter’s n. The above formula gives C in British units,
which could be converted into metric units (exercise)
The Manning equation
R 2/3
So 1 A5 / 3
V Q
n P 2 / 3 S o1 / 2
n
a few typical values of Manning's n
Q AC RS o A5 / 3
So, K ACR1/ 2
Q KS o1 / 2 nP 2 / 3
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.
Computations in uniform flow
We can use Manning's formula for discharge to calculate steady
uniform flow.
Two calculations are usually performed to solve uniform flow
problems.
1. Discharge from a given depth
2. Depth for a given discharge
In steady uniform flow the flow depth is know as normal depth.
As we have already mentioned, and by definition, uniform flow
can only occur in channels of constant cross-section (prismatic
channels) so natural channel can be excluded.
However we will need to use Manning's equation for gradually
varied flow in natural channels - so application to natural/irregular
channels will often be required.
Optimal Shape of Cross-Section
The most hydraulically-efficient shape of channel is the one which
can pass the greatest quantity of flow for any given area or,
equivalently, the smallest area for a given quantity of flow.
Rapid change can also occur when there is a change from super-
critical to sub-critical flow (see later) in a channel reach at a
hydraulic jump.
The application of Energy Equation
p V 2
Recalling the Bernoulli equation,
g 2 g
z cons tan t
Bernoulli
And assuming a hydrostatic pressure distribution we can write
the pressure at a point on a streamline, A say, in terms of the
depth d (the depth measured from the water surface in a
direction normal to the bed) and the channel slope.
pA =ρgy2
y2
In terms of the vertical distance d y1 cos
cos So, p A gy1 cos 2
y 2 y1 cos 2
Thus we have a cubic function with the only unknown being the
downstream depth, y2. There are three solutions to this but only one is
correct for this situation. We must find out more about the flow before we
can decide which it is.
Specific Energy
Specific energy, Es, is
defined as the energy of the
flow with reference to the Fig. 1.10 Specific energy
channel bed as the datum. (definition sketch)
The concept of specific
energy was first introduced
by Bakmenteff (1918).
With reference to the figure
the total energy of flow
with respect to the channel
bottom is given by
p1 V12
E s1 ( y y1 ) V12
g 2 g y
2g
Specific Energy…
Thus the specific energy at an open channel section is equal to the sum of
the flow depth and the velocity head.
In the above equation V1 denotes the velocity of flow at the point of
interest, in the figure above at point 1. In practice it is easier to use the
average velocity of flow at the section and speak about the specific energy
of the flow at a section.
How ever the velocity of flow changes from point to point with in the flow
and as a result the specific energy changes from stream line to stream line.
It is common to use the average velocity of flow with a correction factor.
The specific energy computed using the average velocity is taken to apply
for all points in the section, i.e. is taken as the specific energy of the
section. For steady flow this can be written in terms of discharge Q
(Q / A) 2 q 2
Es y Es y
2g 2 gy 2
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.
Critical, Sub-critical and super critical flow
The specific energy change with depth was plotted above for
a constant discharge Q, it is also possible to plot a graph with
the specific energy fixed and see how Q changes with depth.
For a fixed discharge:
Vc2 y
E sc y c yc c
2g 2
2
y c E sc
3
Variation of the Discharge with depth for a given specific
energy value
The Froude number
V
The Froude number is defined FN
gDm
for channels as:
Its physical significance is the Inertial Force
FN2
ratio of inertial forces to Gravitational Force
Q Q
Q Q
y1b y2b
y 2 y1 3
E
4 y1 y 2
These are useful results and which can be used in gradually varied flow
calculations to determine water surface profiles.
The higher the upstream Froude number the higher the jump and the
greater the loss of energy in the jump.
CHAPTER FIVE
5.0 RAPIDLY VARIED FLOW
Characteristics of RVF
STEADY NON-UNIFORM FLOW
Pronounced curvature of streamlines.
Abrupt change of flow profile (virtually broken)
Resulting the state of high turbulance
Example: Hydraulic Jump
RVF vs GVF
In view of contrast with UF & GVF the following characteristics should be noted.
Pronounced curvature hydrostatic pressure distribution can not be assumed
Rapid variation in flow regime takes place in a very short distance
Effect of boundary friction is comparatively small, which would play a primary role in a
GVF
In RVF the velocity-distribution coefficients and are much greater than unity and can
not be accurately determined.
Flow is actually confined by separation zones as well as solid boundaries. (Because profiles
could be broken).
Approach to the problem
The theory that assumes:
Parallel flow,
Hydrostatic distribution of pressure- Does not apply in RVF computation.
For RVF of continuous flow profile a mathematical equation can be established,
Approach to the solution of such equation
Graphical method (e.g. flow-net analysis)
Numerical method (e.g. method of relaxation)
No satisfactory general solution has yet been obtained
Practical approach
No general solution yet been found
Various RVF phenomena are treated as isolate cases
Each with own semi-empirical/empirical treatment
Experimental results are used empirically
Flow interpreted qualitatively using energy principle, momentum principle, geometry
plus sometimes dimensional analysis
Three isolated cases
Flow over spillway
Hydraulic jump
Flow under gate
Flow over spillways
Definition:
Spillway: is a structure over or through a dam for discharging flood flows; overflow
channel; opening built into a dam or the side of a reservoir to release (to spill) excess
floodwater.
SHARP-CRESTED WEIR (SCW) VS BROAD CRESTED WEIR (BCW)
BCW
Overflow structure with horizontal crest
above which the deviation from a hydrostatic pressure distribution
because of centripetal acceleration may be neglected.
stream-lines are parallel and straight
Criteria 0.5 H1/L0.07
If 0.07 H1/L the energy loss above the crest can not be neglected
0.5 H1/L so that the hydrostatic pressure distribution can be assumed
L = length of the weir crest in the direction of flow,
H1 total energy head over the weir crest
SCW
Overflow structure (H1/L > 15)
The crest length in the direction of the flow is short enough not to influence the H-Q
relationship of a weir
In practice, 0.002m L so that even at a minimum head of 0.03m the nappe is
completely free from the weir body after passing the weir no adhered nappe can
occur
An air pocket beneath the nappe form from which a quantity of air is removed
continuously by the over falling jet.
Therefore, Precaution is required not to ensure that the pressure in the air pocket is not
reduced. Otherwise resulting undesirable effects:
Owing to the increase of the under pressure the curvature of the over falling jet will
increase, causing increase of the discharge coefficient
Irregular supply of air to the pocket will cause vibration of the jet resulting an
unsteady flow
SCW is the simplest form of overflow spillway
Motto:
Spillways must discharge the peak flow under smallest possible head.
Negative pressure on the crest must be limited to avoid danger of cavitation on the
crest or vibration of the structure.
Theoretically, there should be atmospheric pressure on the crest
In same time t, the particle will travel1 a vertical distance y (taking y is positive downward)
y Vo t Sin gt 2
(2)2
Move the origin up so that it coincides with the peak.
1
y Vo ts in gt 2 C '
2
x
t
V Cos
Eliminating t from 1 & 2, from 12 o
x 1 x
y Vo Sin g C '
V0 Cos 2 Vo Cos
hv Dividing each term by the total head H above the crest
H
Y X gH C'
tan HX 2
H H 2Vo2 Cos 2 H
gH C'
Let
A , 2 CosB=
2 - tan, & C
2V 0 H
2
Y X X
A B C
Then, H H
H
= General equation for the lower surface of the nappe in dimension less term
Since the horizontal velocity component is constant, the vertical thickness of the nappe T may be
assumed constant and T
D
H
Adding a term to the above equation the
2 general equation for the upper surface of
the nappe is Y X X
A B CD
H H
H
These equations are quadratic hence, the nappe Surfaces are theoretically parabolic.
Several experimental studies on the nappe over a sharp-crested weir have been made.
Reputed works has been done by US Bureau of Reclamation, they developed the following
equations for the constants in
hv the general nappe equations.hv hv
2
hv
H H 1.568 0.892 0.127
A = - 0.425 + 0.25 B = 0.411 - 1.603 - H H
hv
H
C = 0.150 – 0.45 D = 0.57 – 0.02 (10m) 2 exp (10m)
Where:
hv = the velocity
hv head of the approach flow
m= - 0.208 H
For high weirs, the velocity of approach is relatively small and can be ignored (hv 0)
A = 0.425
B = 0.055
C = 0.150
D = 0.559
Experimental data have indicated that these equations are not valid When,
X hv
H < 0.5 and that > 0.2
H
For < 0.5, The pressure with in the nappe in the Vicinity of the weir crest is > Patm
because of the convergence of the streamlines.
Consequently, forces other than gravity are acting on the nappe,
which makes the principle of the projectile invalid.
N.B: The above theory and equations apply only if the approach flow is sub critical.
For Supercritical flow, or Fr >1, the nappe profile becomes essentially a function of the
Froude number rather than a function of the boundary geometry as described above.
Aeration of the Nappe
In the preceding discussion the over falling nappe is considered a crated; i.e., The upper
and lower nappe surfaces are subject to full atmospheric pressure.
In practice,
- Usually insufficient aeration below the nappe occurs due to removal of air by over
falling set.
Effects of reduction of pressure
Increase in pressure difference on the spillway itself
Change in the shape of the nappe for which the spillway crest is designed
Increase in discharge, sometimes accompanied by fluctuation or pulsation of the
nappe, which may be very objectionable if the weir or spillway is used for measuring
purposes.
Unstable performance of the hydraulic model
Crest Shape of Overflow Spillways
Earliest shapes were based on a simple parabola designed to fit the trajectory of the falling nappe (the
equation for the lower surface of the nappe).
Bazin’s made comprehensive laboratory investigation for nappe shapes.
the used of Bazin’s data in design will produce a crest shape that conincides with lower surface of as
aerated nappe over a sharp-crested weir.
Such a profile is known as Bazin profile Advantage
Should cause no negative pressure on the crest (the presence of negative pressure will lead to
danger of cavitation damage).
In selecting a suitable profile avoidance of negative pressure should be considered an objective,
along with such other factors as maximum hydraulic efficiency, practicability, stability &
economy.
Extensive experiments on the shape of the nappe over-sharp crested weir were conducted
by U. S Bureau of Reclamation; including Bazin’s,
The Bureau has developed coordinates of the nappe surface for various slope faced weirs
On the basis of the Bureau data, The U.S Army. Corps of Engineers has developed several
standard shapes at its Waterways Experimental Station.
Such shapes designed as the WES standard spillway shapes, can be expressed by the
following equation:-
n 1
X n K Hd Y
X and Y are Coordinates of the crest profile with the origin at the highest point of the
crest.
Hd is the design head excluding the velocity head of the approach flow
K & n are parameters depending on the slope of the upstream face. values of k & N are
given as flows:
Slope of upstream face k n
Vertical 2.000 1.850
3 :1 (V = H) 1.936 1.836
3:2 1.939 1.810
3:3 1.873 1.776
For intermediate slopes: approximate value of k and n may be obtained by plotting the
above values against the corresponding slopes and interpolating
from the plot the required values for any given slope within the plotted range.
The upstream face of the spillway crest may some times be designed to set back, as shown
by the dashed lines
Discharge of WES Spillway
The discharge over a spillway can be computed by an equation in the form of SCW/BCW
Q = CLHe1.5
He- the total energy head on the crest, including the velocity head in the approach canal.
The effect of the approach velocity is negligible when height h of the spillway is greater
than 1.33Hd (h > 1.33 Hd), where the design head exclude the approach velocity head.
Under this condition, i.e. h/Hd > 1.33, He = Hd can be taken (the approach velocity head is
negligible) and the coefficient of discharge C has been found to be C = 2.21 (if is in ft C He =
4.03)
y1 2
This has been verified with experiments
Types of Jump
Hydraulic Jumps on horizontal floor are of several distinct types. They can be conveniently classified
according to Froud Number Fr1 of the incoming flow as follows.
Fr1 =1 critical flow no jump can form
1< Fr1 < 1.7 the water surface shows undulation (undular jump)
1.7 < Fr1 < 2.5 a series of small rollers develop on the surface of the jump,
but the d/s water surface remains smooth.
The velocity throughout is fairly uniform, and the energy loss low.
The jump is called weak jump.
2.5 < Fr1 < 4.5 there is an oscillating jet entering the jump bottom to surface and back again with no
periodicity.
Each oscillation produces a large wave of irregular period which, very common in canals,
can travel for miles doing unlimited damage to earth banks and ripraps.
This jump is called Oscillating Jump.
4.5 < Fr1 < 9.0 steady Jump:-
The down stream 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.
The jump is well balanced and the performance is at its best.
The energy dissipation ranges from 45 to 70%.
Fr > 9.0 Strong jump:-
The high-velocity jet grabs intermittent slugs of water rolling down the front face of the
jump,
generating waves down-stream and a rough surface can prevail.
The jump action is rough but effective since the energy dissipation may reach 85%.
N.B. It should be noted that the ranges of the Froude Number given above for the various
types of jump are not clear-cut but overlap to a certain extent depending on local
conditions.
4 y1 y 2
E
Re lative loss : the ratio
E1
Efficiency: the ratio of the specific energy after the jump to that before the jump is defined
as the efficiency of the jump.
E2
2
3
8F1 1 2 4 F1 1
2
E1 2
8F1 2 F1
2
The relative loss is equal to , E2
1
E1
this also is a dimensionless function. of Fr1.
Height of Jump:- the difference between the depths after and before the jump.
Hj = y2 – y1
Expressing each term as a ratio with respect to initial specific energy. y2
hj
y1
h1 y2 y
1 E1
E1
E1
E1 E1 E1
Where is the relative height,
is the relative initial depth, and
2
hj 1 8 F1 3
2
is the relative sequence depth. E1 F1 2
Length of Jump:
The length of a jump (also length of stilling basin) is empirically given as,
L= K(y2-y1)
Where, k – is a coefficient derived from laboratory and filed experiment. 4.5 < k < 5.5 where the lower k = 4.5 applies
of Fr2 > 10 and the higher k=5.5 for Fr2 < 3.
Gates
Flow inunder
canals are mainly used as water level regulators.
Gates
Sometimes, gates are used as discharge regulator (measuring device).
They are under-shot or underflow structures. Example slice gate, radial gate roller gate
The design of underflow gate focuses on head-discharge relationship (Q-H).
The objective is to minimize head loss; this means that the gate has to be lifted out off the
water for design discharge.
The other concern of the design is the pressure distribution over the gate as a function of
opening and gate form.
The H-Q relationship for gate depends on the shape and dimension of the control section
and the resulting curvature of the streamlines.
For gated structures the control section is defined by the vena contract, being the smallest
cross section just down steam of the gate.
In the vena contract, streamlines are straight and parallel.
In gate flow 3 flow types can be distinguished.
1) Free flow: the opening is relatively small (h1/a >2), and the contraction of the steam-
lines in vertical direction is strong.
The down stream water level (h2) won’t affect the flow underneath the gate and a
hydraulics jump will occur down stream of the vena contra.
The discharge depends up on the gate opening the upstream water level and the
contraction coefficient.
Submerged flow: the d/s water level influences the flow underneath the gate. The
hydraulic jump is drowned and the jet underneath the gate is submerged. The discharge
depends upon the upstream and downstream water level and the gate opening.
The boundary between free and submerged flow is a sharp one, which can be cleanly found
from the gate opening and the two water levels.
Weir flow: on off gate
The equation for a free flow underneath a sharp edged gate is:
Q Cd Ba 2gh1
Cd = discharge coefficient
B = Width of gate opening
a = height of gate opening
h1 = upstream water depth
Cc
Cd
The discharge coefficient Cd a
1 Cc
h1
Where,
CC = Contraction coefficient of the jet depending on the shape of the gate and
d = diameter of the rounded bottom edge
For d
4.7 Cc 0.99 Rounded edged gates
a
For submerged flow, some equations include the difference between the upstream and
downstream depths and others use the upstream water level only. The general equation is
offer given as. Q C Ba 2gh 2 1
Where,
a = vertical opening of the gate( a< 0.67h1)
h1 = Upstream water depth
B = Effective width of the opening
C2 = discharge coefficient. h h
a 0.67 * h1 or 1 0.67 1 1 1.5
For values of a athe discharge follows from the equation for a broad-
.
SHARP CRESTED WEIR
BROAD CRESTED WEIR
Sketch of
a Venturi
flume
Irrigation Canal
Grand Canal of China
Beijing and Hangzhou
Hydroelectric Power Plant
SPILLWAY