Professional Documents
Culture Documents
Hydraulic Structures - Lecture Notes
Hydraulic Structures - Lecture Notes
LECTURE NOTES
F. O. K. Anyemedu
KUMASI
MARCH - 2007
1
Chapter 1; Hydraulic Structures
Hydraulic Structures
Hydraulic structures are engineering structures constructed for the purposes of
harnessing and using water resources (groundwater, surface water, lakes, sea,
etc) or for the prevention of the negative and destructive actions (floods, shore
erosion, etc) of water on the surrounding environment.
There are a large variety of hydraulic structures to serve the many purposes for
which water resources are put to use. There are several classifications of
hydraulic structures, however, the most important is the classification by function
as given below.
2
Table 1. Classification of Hydraulic Structures by function
There are three common classification schemes for dams. According to the
function performed, dams are classified into:
i) storage dams for impounding water for developmental uses.
ii) Diversion dams for diverting streamflow into canals or other
conveyance system and
iii) Detention dams to hold the water temporary to retard flood flows
3
ii) ii) non-overflow dams, which are not designed to be overtopped.
The most common classification is one based on the materials of which dams are
made. This classification makes further sub-classification by recognizing the
basic type of design, such as concrete gravity or concrete arch dams. Types of
dams include:
i) Earthfill dams
ii) Rockfill dams
iii) Concrete dams
a) concrete gravity dams
b) concrete arch dams
c) concrete buttress dams
iv) Stone masonry
a) stone-masonry gravity dams
b) stone-masonry arch dams
v) Timber dams
vi) Steel coffer dams
4
Chapter 2: EARTHFILL DAMS
Earthfill dams
Earthfill dams are the most common type of dam built to any height. They are
designed as a non-overflow section with separate spillway. The reason for such
wide spread use of earthfill dams are:
• The foundation requirements are not as rigorous as other dams
• Local available soil is the main construction material
• High skill not required
• No special plants are required. Most earth-moving machines can be used
•
5
In this type of dams, successive layers of moistened or damp soils are laid one
over the other. Each layer not exceeding 20 cm in thickness is properly
consolidated at optimum moisture content, only then is the next layer laid.
6
Zoned Earth Dams
These are dams with the central portions called core or hearting made from
materials which are relatively impervious. The thickness of the core wall is made
sufficiently thick to prevent leakage of water through the body of the dam.
7
The criteria for the design of earth dams are:
1. Sufficient spillway capacity and freeboard are provided so that there is no
danger of overtopping of the dam
2. Seepage flow through the embankment is controlled so that the amount lost
does not interfere with the objective of the dam and there is no erosion or
sloughing of soil. In this respect, seepage line should remain well within the
downstream face of the dam and the portion of the dam on downstream side of
the impervious core should be well drained.
3. Uplift pressure due to the seepage underneath is not enough to cause piping
4. The slopes of the embankment are stable under all conditions of reservoir
operation, including rapid drawdown and during steady seepage under full
reservoir.
5. The stresses imposed by the embankment upon the foundation are less than
the strength of material in the foundation with a suitable factor of safety.
6. The upstream face is properly protected ((stone pitching, riprap, revetment)
against erosion caused by wave action, and the downstream face is protected
(counter-booms, turfs) against the action of rain
Embankment Materials
a. Earth-fill materials.
(1) While most soils can be used for earth-fill construction as long as they are
insoluble and substantially inorganic, typical rock flours and clays with liquid
limits above 80 should generally be avoided. The term “soil” as used herein
includes such materials as soft sandstone or other rocks that break down into soil
during handling and compaction.
(2) If a fine-grained soil can be brought readily within the range of water contents
suitable for compaction and for operation of construction equipment, it can be
used for embankment construction. Some slow-drying impervious soils may be
unusable as embankment fill because of excessive moisture, and the reduction
of moisture content would be impracticable in some climatic areas because of
anticipated rainfall during construction. In other cases, soils may require
additional water to approach optimum water content for compaction. Even
ponding or sprinkling in borrow areas may be necessary. The use of fine-grained
soils having high water contents may cause high porewater pressures to develop
in the embankment under its own weight. Moisture penetration into dry hard
borrow material can be aided by ripping or plowing prior to sprinkling or ponding
operations.
8
characteristics) may be difficult to use in an embankment. However, this depends
upon local climatic conditions and the size and layout of the work, and must be
assessed for each project on an individual basis. The cost of using drier material
requiring a longer haul should be compared with the cost of using wetter
materials and flatter slopes. Other factors being equal, and if a choice is possible,
soils having a wide range of grain sizes (well-graded) are preferable to soils
having relatively uniform particle sizes, since the former usually are stronger, less
susceptible to piping, erosion, and liquefaction, and less compressible. Cobbles
and boulders in soils may add to the cost of construction since stone with
maximum dimensions greater than the thickness of the compacted layer must be
removed to permit proper compaction. Embankment soils that undergo
considerable shrinkage upon drying should be protected by adequate
thicknesses of non-shrinking fine-grained soils to reduce evaporation. Clay soils
should not be used as backfill in contact with concrete or masonry structures,
except in the impervious zone of an embankment.
(4) Most earth materials suitable for the impervious zone of an earth dam are
also suitable for the impervious zone of a rock-fill dam. When water loss must be
kept to a minimum (i.e., when the reservoir is used for long-term storage), and
fine-grained material is in short supply, resulting in a thin zone, the material used
in the core should have a low permeability. Where seepage loss is less
important, as in some flood control dams not used for storage, less impervious
material may be used in the impervious zone.
a) Hydraulic failure:
Hydraulic accounts for over 40% of earth dam failure and may be due to one or
more of the following:
9
.
F.W.L
10
Figure 3. Beginning of downstream failure
.
H.F.L
.
Channel formation
11
ii) Piping through foundation: When highly permeable cavities or fissures or
strata of gravel or coarse sand are present in the dam foundation, it may lead to
heavy seepage. The concentrated seepage at high rate will erode soil which will
cause increase flow of water and soil. As a result, the dam will settle or sink
leading to failure.
piping
c) Structural Failure:
About 25% of failure is attributed to structural failure, which is mainly due to
shear failure causing slide along the slopes. The failure may be due to:
i) Slide in embankment: When the slopes of the embankments are too steep, the
embankment may slide causing failure. This might happen when there is a
sudden drawdown, which is critical for the upstream side because of the
development of extremely high pore pressures, which decreases the shearing
strength of the soil. The downstream side can also slide especially when dam is
full. Upstream embankment failure is not as serious as downstream failure.
ii) Foundation slide: When the foundation of an earthfill dam is composed of fine
silt, clay, or similar soft soil, the whole dam may slide due to water thrust. If
seams of fissured rocks, such as soft clay, or shale exist below the foundation,
the side thrust of the water pressure may shear the whole dam and cause its
12
failure. In such failure the top of the dam gets cracked and subsides, the lower
slopes moves outward and forms large mud waves near the dam heel.
iii) Faulty construction and poor maintenance: When during construction, the
compaction of the embankment is not properly done, it may lead to failure.
iv) Earthquake may cause the following types of failure to earthfill dams;
1 cracks may develop in the core wall, causing leakages and piping failure.
2. slow waves may set up due to shaking of reservoir bottom, and dam may
fail due to overtopping
3. settlement of dam which may reduce freeboard causing failure by
overtopping
4. slidding of natural hills causing damage to dam and its appurtenant
structures
5. Fault movement in the dam site reducing reservoir capacity and causing
overtopping.
6. Shear slide of dam
7. The sand below foundation may liquefy
8. Failure of slope pitching.
B = H/5 + 3
B = 0.55(H)1/2 + H/5
13
B = 1.65(H + 1.5)1/3
B = 1.67(H)1/2
Balustrades are provided at the end of the roads to prevent car falling off the
slopes.
2. Side slopes of dam: side slope of dams must satisfy the static stability.
However, since the stability computations can be done only after defining the
profile of the dam and determining the seepage line, it becomes necessary to
give an initial side slopes. Initial slope could be taken from the tables below.
(Taken from “hydraulic structures” –N. P. Rosanova)
Slope Material Side slopes depending on height of dam
of dam Less than 5m From 5m - From 10-15m
10m
Upstream Clayey 2 2.5 3
Sandy 2.5 – 2 3 – 2.5 3
Downstream Clayey 1.5 1.75 1.75
With filter Sandy 2 2 2
Downstream Clayey 1.75 2 2.25
without filter Sandy 2 2.25 2.25
In low head dams, usually one and constant side slope is used, however in
medium and high head dams, different side slopes are usually adopted to reduce
the volume of the dam.
3. Berms: Berms are constructed at both the upstream and downstream side of
the dam for the purposes of observing the conditions of protections at the
slopes and their repairs and also for increasing the width of dam at the base
14
with aim of increasing seepage length. It is also done when constructions
coffer dams are made part of the body of the dam. At the downstream side,
berms are done at an interval of 10 – 15m high. Width of berm is taken
between 1 – 2 m.
4. Free Board
Normal freeboard is the vertical distance between the normal pool level and the
crest of the dam. Minimum freeboard is the vertical distance between the high
flood level and the crest of dam.
The minimum height of freeboard is taken as 1.5 hw where hw is given by:
hw = 0.032 (V.F)1/2 + 0.763 – 0.271(F)1/4 for F, 32 km -----------(X)
and hw = 0.032 (V.F)1/2 for F . 32km -------------------------------------(Y)
where hw = wave height ( height of water from top to trough of waves in meters)
V = velocity of wind in km/hr
F = fetch or straight length of water expanse in km.
Free board values as recommended by U.S.B.R are given in table below.
15
Upstream side protection: For protecting the upstream slope from deterioration
and damage from wave action, the slope is covered with different protective
material.
Rock riprap, either dry dumped stone boulders or hand-packed stone boulders
can be made. Stone pitching provided at slope of 1.5 : 1 to 2 : 1 for ordinary soil
material of dam and 3 : 1 for poor soil material. The revetment stones are fixed
at the toe of the dam to prevent slipping of the embankment. The thickness of the
stone pitching is usually more than 60 cm. In most cases the stone pitching is
placed over gravel then sand cushion. Big sizes stones with their broader face
downwards are packed with each other by means of hammer.
Concrete, reinforced concrete slabs, steel plates, bituminous material pavement,
brick tile pavement could also be used. However, extensive survey done by US
Corps of Engineers in the 1940s o over 100 dams show that dry dumped riprap
stone pitching has the best efficiency in terms of failure rate.
16
Downstream protection
One of the cheapest, simple and effective methods of protecting the downstream
from rainfall and wind action is by planting green grass (turfs) on the slope.
Counter-boom could also be done.
17
Figure 4. Downstream protection
6. Drainage: Drainage in earth dams is meant for lowering the seepage curve;
prevent seepage water from flowing onto the downstream slope, and
conveying seepage water through the body of the dam to the downstream
part of the dam. By its function, dam drainage must have two parts; an intake
structure (drainage trench) that allows seepage water from the body and
foundation of the dam, while at the same time preventing deformation due to
seepage and conveyance structure that transports the seepage water from
the dam. However, in most drainage, it is difficult to see clearly these two
parts.
18
Hydrotechnical construction practice has work out many drainage systems
depending on the type of dam, materials of the foundation and body of dam.
Among some of the most commonly used drainage systems are:
i) Drainage prism: with many positive sides (advantages) but requires the use of
large quantity of stones
ii) A type of drainage prism in which the filter material of the drainage system is
laid to extend to a certain height on the downstream side. Such drainage
system is used when there could be rise of the tail water above the crest of
the prism.
iii) Flat horizontal drainage: It requires much smaller quantity of stones and
simplifies construction. It has the advantage of draining both the foundation
and body of dam and it is used mainly when the foundation is made up of
saturated material
iv) Combination of horizontal drain with the prism.
v) Horizontal piped drainage: consists of a pipe (tube), laid parallel to base of
the slope of dam.
vi) Horizontal stone drainage: a type of horizontal piped drain in which instead of
the pipe, a stone prism is used.
19
Seepage through dam
Filtration (seepage) computation through dams are carried out with the aim of:
• determining the position of the seepage depression curve
• finding the value of the gradient and velocity of filtration
• determine the filtration flow (discharge)
The seepage pattern through a dam is shown above. The pattern is the same
irrespective of the material (sand, clay, loam) of the dam, though the rate of
seepage will depend on soil type. The emergence of seepage lines on the
downstream slope tends to make the downstream slope unstable. Either the
downstream slope has to be made very flat or the seepage must be diverted
away from the downstream slope. The second alternative is favoured because it
is economical.
The diagram below shows the seepage lines for a homogeneous earthfill dam
with a horizontal gravel filter on the base of the dam at the downstream side. The
flow lines
20
enter the blanket vertically. Cassagrande (1937) has shown that the phreatic
line, which is the topmost seepage line, quite closely approximates a parabola.
The parabola
Insert drawing
21
intersects the water surface at A such that AB = 0.3 CB. Near the upstream
face, the phreatic line diverges from the parabola and join B perpendicularly.
The focus of the parabola is at F and the directrix passes through D. Every point
on the parabola with origin at the focus F can be given by
y 2 − y 02
x= ----------------1
2 y0
At the point A, x = d and y = H. Substituting in equation 1 yields
y0 = H 2 + d 2 − d ------------(2)
y0 , which is also the distance FD, is determined from eq. (2)
How to get the point D graphically. With A as the centre and AF as radius, draw
an arc to cut the water surface at J. From J drop a perpendicular to the base at
D. D is the directrix. Locate the mid-point of FD at E which is the vertex of the
parabola.
The phreatic line can be drawn by equation (1).
For an approximate estimate of seepage, the Darcy’s law at point G , per unit
length of dam
dy
q = KiA = K y 0 -------------- (3)
dx
dy y0
From eq. (1) dx = --------------------(4)
2 xy 0 + y 02
Since at G, x= 0, and y = y0 then
dy/dx = 1 ------------------------------(5)
Hence q = K (1) (y0) ------------------------------(6)
The flow through the section at G is the same as at any other section.
Consequently, the total seepage through the length L of the dam is given by:
Q = K(y0)(L) ------------------------------(7)
Worked example
A homogeneous earthfill dam has a top width of 30ft and a height of 100ft with a
freeboard of 10ft. The side slopes are 1V: 2.5H. It has a horizontal drainage
blanket at the base that extends from the downstream toe to a distance of 100ft.
The embankment has a permeability of 1.5 x 10-5 ft/sec. Determine the seepage
through the dam.
Solution
H = 90 ft ; CB = 90 (2.5) = 225ft ; AB = 0.3(225) = 67.5ft and CA = 157.5 ft
d = base length – blanket length – CA
22
= 530 -100 – 157.5 = 272.5 ft.
Then from Darcy’s law of flow through porous medium, considering unit
thickness,
∆h
∆q = k . ( b.1) --------------------a)
l
23
equipotential
lines
Flow lines
h b N b
Or Total flow q = ∑ ∆q = k . .N t = kh t . ----------- d)
Nd l Nd l
Where Nt = total number of flow channels in the net.
When the flow field is square, then b = l and the discharge becomes
Nt
q = kh ---------------------f)
Nd
24
The submerged unit weight soil is given by
G −1
γ sub = s γw ----------(Y)
1+ e
where Gs = specific gravity of soil; e = void ratio.
For a seepage line at a gradient, I, the upward seepage force per unit volume is
iγw . When the two forces are in balance,
Gs − 1
i= -----------------(Z)
i+e
Eq. (Z) is known as the critical gradient and equal to unity. A gradient of slightly
higher than unit value will cause piping or sand particles to be in an unstable
condition known as quicksand. The actual gradient at the downstream end of the
dam is evaluated from the flow net by dividing the head difference between the
lst two potential lines by the distance between these potential lines. This should
be less than unity.
An empirical approach, the creep ratio, L/H is computed; here L is the length
along the surface of contact between the soil and the base of the structure. This
ratio is kept at 4 for gravels and 18 for sand and silt.
The following are some common measures for controlling seepage through the
dam and embankments.
1. Prevention of Seepage through Foundation
i) By providing drainage trenches
ii). By providing downstream seepage berms
iii). By providing impervious blanket layer on upstream slope
iv). By providing impervious cutoff.
25
ii). By providing toe filter
iii). By providing filter downstream of toe
iv). By providing downstream coarse section
v). By providing chimney drains extending upwards into the embankment
Chimney drain
Design of Filter
The design of filters should be done in such a way that all the seepage water
through the dam is effectively drained off. The filter consists of several layers.
The first layer of the filter which comes in contact with the seeping water consists
of fine sand material. Subsequent layers of filter are made of sand of increased
fineness. The last layer of the filter is made of gravels. The soil of the earth dam
and the foundation material surrounding the filter are known as the base material.
The filters of filter drains are known as reverse or inverted filter.
The U.S. Bureau of reclamation, Washington 1960 recommended the following
for materials to be used for filters;
i) Filter material should be fine and poorly graded so that the voids in
the filter are small and thus prevent base material from entering the
filter
ii) The filter material should be coarse and pervious in relation to the
base material. This aspect facilitates rapid removal of seeping
water without building up any seepage forces within the filter.
iii) The filter material should be coarser than the perforations of
openings in the drain pipes, so that filter material is not lost in the
drains. The perforations (openings in the pipes drains should be
adequate to admit all seeping water safely.
iv) The thickness of filter material should be sufficient to provide a
good distribution of all particle sizes, also throughout the filter. The
thickness should be adequate to provide safety against piping.
Terzaghi has recommended the following two requirements which should be
fulfilled by the filter:
a). The D15 size of the filter material must not be more than 4 to 5 times the D85 of
the base material. This prevents the foundation material from carrying through
the pores of the filter material.
26
b). The D15 size of the filter material should be at least 4 to 5 times the D16 of the
base material. This keeps seepage forces within the filter to permissible levels.
The above criteria can be expressed as follows:
The above criteria has been modified by USBR in “ Design of small Dams ” as:
a)
a) D15 of filter material = 5 to 40
D15 of base material
d). The grain size curve of the filter material should be about parallel to the curve
of the base material.
Earthfill dams usually fail due to the sliding of large soil mass along a curved
surface.
The most common method used for examining the stability of earthfill dam
embankment slopes is the Swedish slip circle or (the slice) method.
27
P
β
WL
S
α
H
i
R 4.5 H
Q
For a homogeneous soil, the centre of critical slip centre lies on the line PQ. The
coordinates of point Q is H downwards from toe and 4.5 H horizontal away from
toe as shown above. The location of point P is done with the help of directional
angles i, α and β given in the table below.
First extend the surface of the fluid WL. Draw angle i from the toe of dam to cut
the surface of fluid at S. From S draw the angle β and from the toe draw angle α
to intercept angle β at the point P. Now by joining QP, we obtain the line on which
the centre of critical slip centre lies.
Now on line QP we assume point O and an arc AB is drawn with radius OA. The
earth mass above the arc is divided into a number of vertical segments known as
slices. The number of slices may be from 6 to 12. Neglecting the forces between
slices, each slice is assumed to act independently as a column of soil of unit
thickness and width. The weight of each slice is assumed to be acting at its
centre. The weight is resolved into two components N and T (ie. normal
28
component passing through the centre of the arc O, and will not cause any
driving moment on the slice. But the tangential component T causes a driving
moment of T x R). Near the toe of the dam, some of the slices may cause
resisting moment and in such cases, T is negative.
O
R O
A
R
N
α
B
W
Let thew slices be numbered 1, 2, …6 and their weights be numbered W1, W2…
W6
Now N = W cos α and T = W sin α
From Coulomb’s equation, the resisting force
FR = c x ΔL + N tan φ ……….1)
where c = unit cohesion ; ΔL = curved length of slice ;
φ = angle of internal friction of soil.
F .S =
[
M R R c ∑ ∆L + tan φ ∑ N
=
]
MD R∑ T
F .S =
[c∑ ∆L + tan φ ∑ N ] .................(3)
∑T
The values of ∑N and ∑T are generally calculated in a tabular form.
29
Now if w1, w2, …….wn are the weight of slices, then
∑N = cos α (∑w) and ∑T = sin α∑(w)
The factor of safety is computed for the point O with equation (3). Different
points on the line QP are chosen and the procedure gone through to deyermine
the factor of safety. Various factor of safety so obtained are plotted as ordinates
on the corresponding centres and a smooth curve of F.S drawn. The centre
corresponding to the lowest factor of safety is the required centre of critical slip
circle.
The critical condition for downstream slope occurs when the reservoir is full and
seepage is at its maximum rate. The seeping water below the phreatic surface
exerts a pore pressure on the soil mass. Therefore, the shearing strength of the
slices of the critical arc within the range of pore-pressure is reduced.
The net shear strength of the affected slice is
If dam is suddenly emptied, it may cause critical conditions for the stability of the
upstream slope. In this condition, the soil pores remain filled with water causing
30
the water level to remain the same as before the sudden drawdown. The water
weight within the soil tends to slide the upstream slope along a circular arc. The
hydrostatic force acting along the upstream slope when reservoir is full is also
removed with the sudden drawdown. Therefore, the tangential component of the
saturated weight is the main force causing disturbance to the upstream slope.
The shear resistance is considerably reduced due to development of pore
pressure on the likely slip surface.
The factor of safety can be calculated by the formula (6) above.
At dam site, the silt and clay or fine, loose cohesionless material that form the
foundation may have good imperviousness, but are weak in shear and always
require checking.
For increasing the shear area and keeping the shear stress within permissible
limits, the embankment slopes are flattened or berms are added on both sides of
the dam.
0.4B
Toe
h1 h
h2 J L Foundation material
Firm stratum B
The method for determining the factor of safety, which is approximate is based
on the assumption that earthfill materials have an equivalent liquid weight, which
would produce the same shear stress as the material will develop itself.
The horizontal shear on left is given by:
31
h12 − h22 φ'
P= .γ m tan 2 45 o − ...................(7)
2 2
φ '
where γ m . tan 2 45 o − .........................equivalent.liquid .unit.weight
2
γ d ( h1 − h2 ) + γ f .h2
γm =
h1
Chapter 3
GRAVITY DAMS
32
Basically, gravity dams are solid concrete structures that maintain their stability
against design loads from the geometric shape and the weight and strength of
the structure. Generally, they are constructed on a straight axis, but may be
slightly curved or angled to accommodate the specific site conditions. Gravity
dams typically consist of a non-overflow section(s) and an overflow section or
spillway. They are constructed with masonry or concrete but of late conventional
concrete or roller-compacted concrete are popular.
The two general concrete construction methods for concrete gravity dams are
conventional placed mass concrete and roller-compacted concrete (RCC).
33
the concrete approaches its annual temperature cycle. Control measures involve
precooling and postcooling techniques to limit the peak temperatures and control
the temperature drop. Reduction in the cement content and cement replacement
with pozzolans have reduced the temperature-rise potential. Crack control is
achieved by constructing the conventional concrete gravity dam in a series of
individually stable monoliths separated by transverse contraction joints. Usually,
monoliths are approximately 50 feet wide.
Site Selection
a. General. During the feasibility studies, the preliminary site selection will be
dependent on the project purposes. Purposes applicable to dam construction
include navigation, flood damage reduction, hydroelectric power generation, fish
and wildlife enhancement, water quality, water supply, and recreation. The
feasibility study will establish the most suitable and economical location and type
of structure.
b. Selection factors.
(1) A concrete dam requires a sound bedrock foundation. It is important that the
bedrock have adequate shear strength and bearing capacity to meet the
necessary stability requirements. The foundation permeability and the extent
and cost of foundation grouting, drainage, or other seepage and uplift control
measures should be investigated.
34
location would minimize the concrete material requirements and the associated
costs.
(3) The criteria set forth for the spillway, powerhouse, and the other project
appurtenances will play an important role in site selection. The relationship and
adaptability of these features to the project alignment will need evaluation along
with associated costs.
(5) The method or scheme of diverting flows around or through the damsite
during construction is an important consideration to the economy of the dam. A
concrete gravity dam offers major advantages and potential cost savings by
providing the option of diversion through alternate construction blocks, and
lowers risk and delay if overtopping should occur.
3-3. Loads
a. General. In the design of concrete gravity dams, it is essential to determine the
loads required in the stability and stress analysis. The following forces may affect
the design:
(1) Dead load. (2) Water Pressure (Headwater and tailwater pressures).
(3) Uplift. (4) Temperature. (5) Earth and silt pressures. (6) Ice pressure.
(7) Earthquake forces. (8) Wind pressure. (9) Subatmospheric pressure.
(10) Wave pressure. (11) Reaction of foundation.
b. Dead load.
35
Water pressure is the major external force acting on the dam. As the water is
stored in the reservoir, and stands against the body of the dam, it exerts
horizontal pressure on the dam.
P=1/2γh
Dam
h/3
Head water
h
Tail water
P1
P2
Ww Ww
Fig. G.D. 2 Water pressure for slanted upstream face and water at tailrace
When the upstream face of the dam is vertical, the water pressure P = ½ γh2 and
acting at h/3 from base.
When the upstream face is slanted, the water pressure is resolved in two
components:
Horizontal oressure P = ½ γh2 and acting at h/3 from base.
Vertical pressure Ww = weight of water on slanted side and acting at centre of
gravity of volume of water.
Uplift Pressure.
It is the second major external force acting upwards on the dam.
36
Uplift pressure resulting from headwater and tailwater exists through cross
sections within the dam, at the interface between the dam and the foundation,
and within the foundation below the base. This pressure is present within the
cracks, pores, joints, and seams in the concrete and foundation material. Uplift
pressure is an active force that must be included in the stability and stress
analysis to ensure structural adequacy. These pressures vary with time and are
related to boundary conditions and the permeability of the material.
(a) General. The uplift pressure will be considered as acting over 100 percent of
the base. A hydraulic gradient between the upper and lower pool is developed
between the heel and toe of the dam. The pressure distribution along the base
and in the foundation is dependent on the effectiveness of drains and grout
curtain, where applicable, and geologic features such as rock permeability,
seams, jointing, and faulting. The uplift pressure at any point under the structure
will be tailwater pressure plus the pressure measured as an ordinate from
tailwater to the hydraulic gradient between upper and lower pool.
(b) Without drains. Where there have not been any provisions provided for uplift
reduction, the hydraulic gradient will be assumed to vary, as a straight line, from
headwater at the heel to zero or tailwater at the toe. Determination of uplift, at
any point on or below the foundation, is demonstrated in Figure G.D.3 above.
(c) With drains. Uplift pressures at the base or below the foundation can be
reduced by installing foundation drains. The effectiveness of the drainage system
will depend on depth, size, and spacing of the drains; the character of the
37
foundation; and the facility with which the drains can be maintained. This
effectiveness will be assumed to vary from 25 to 50 percent, and the design
memoranda should contain supporting data for the assumption used.
Along the base, the uplift pressure will vary linearly from the undrained pressure
head at the heel, to the reduced pressure head at the line of drains, to the
undrained pressure head at the toe, as shown in Figure G.D.4
Where the line of drains intersects the foundation within a distance of 5 percent
of the reservoir depth from the upstream face, the uplift may be assumed to vary
as a single straight line, which would be the case if the drains were exactly at the
heel. This condition is illustrated in Figure G.D.5
If the drainage gallery is above tailwater elevation, the pressure of the line of
drains should be determined as though the tailwater level is equal to the gallery
elevation.
38
Figure G.D.5. Uplift distribution with foundation drains near upstream face
(e) Zero compression zones. Uplift on any portion of any foundation plane not in
compression shall be 100 percent of the hydrostatic head of the adjacent face,
except where tension is the result of instantaneous loading resulting from
earthquake forces. When the zero compression zone does not extend beyond
the location of the drains, the uplift will be as shown in Figure G.D.6. For the
condition where the zero compression zone extends beyond the drains, drain
effectiveness shall not be considered. This uplift condition is shown in Figure
G.D.7
39
Figure G.D.7. Uplift distribution cracked base with drainage, zero
compression zone extending beyond drains(3-5)
.
When an existing dam is being investigated, the design office should submit a
request to CECW-ED for a deviation if expensive remedial measures are
required to satisfy this loading assumption.
(b) RCC concrete. The determination of the percent uplift will depend on the mix
permeability, lift joint treatment, the placements, techniques specified for
minimizing segregation within the mixture, compaction methods, and the
treatment for watertightness at the upstream and downstream faces. A porous
upstream face and lift joints in conjunction with an impermeable downstream face
may result in a pressure gradient through a cross section of the dam
considerably greater than that outlined above for conventional concrete.
Wave Pressure
40
The portions of the dam is subjected to the impact of waves, which is produced in
water surface during wind
Water surface
Pw
hw 5/3hw
2.4γh
Figure G.D. 8
Wind Load
Wind load in stability analysis is usually ignored
Earth pressures against the dam may occur where backfill is deposited in the
foundation excavation and where embankment fills abut and wrap around
concrete monoliths. The fill material may or may not be submerged. Silt
pressures are considered in the design if suspended sediment measurements
indicate that such pressures are expected. Whether the lateral earth pressures
will be in an active or an at-rest state is determined by the resulting structure
lateral deformation.
41
Earthquake Forces
(1) General.
(a) The earthquake loadings used in the design of concrete gravity dams are
based on design earthquakes and site-specific motions determined from
seismological evaluation. As a minimum, a seismological evaluation should be
performed on all projects located in seismic zones.
(b) The seismic coefficient method of analysis should be used in determining the
resultant location and sliding stability of dams. In strong seismicity areas, a
dynamic seismic analysis is required for the internal stress analysis.
(c) Earthquake loadings should be checked for horizontal earthquake
acceleration and, if included in the stress analysis, vertical acceleration. While an
earthquake acceleration might take place in any direction, the analysis should be
performed for the most unfavorable direction.
42
Figure G.D. 9. Seismically loaded gravity dam, nonoverflow monolith
(a) Inertia of concrete for horizontal earthquake acceleration. The force required
to accelerate the concrete mass of the dam is determined from the equation:
Pe = Max =( W/g)αg = W α
(b) Inertia of reservoir for horizontal earthquake acceleration. The inertia of the
reservoir water induces an increased or decreased pressure on the dam
concurrently with concrete inertia forces. FigureG.D.9 shows the pressures and
forces due to earthquake by the seismic coefficient method. This force may be
computed by means of the Westergaard formula using the parabolic
approximation:
where
Pew = additional total water load down to depth y (kips)
Ce = factor depending principally on depth of water and the earthquake vibration
period, t , in seconds e
h = total height of reservoir (feet)
43
Westergaard's approximate equation for Ce, which is sufficiently accurate for all
usual conditions, in pound-second feet units is:
51
Ce =
2
h
1 − 0.72
1,000t e
Analysis of the stability and calculation of the stresses are generally conducted
at the dam base and at selected planes within the structure. If weak seams or
planes exist in the foundation, they should also be analyzed.
Dams are designed for the most adverse combination of load conditions as have
reasonable probability of simultaneous occurrence. The following basic loading
conditions are generally used in concrete gravity dam designs (see Figure G.D
10).
44
(1) Load Condition No. 1 - unusual loading condition - construction.
(a) Dam structure completed. (b) No headwater or tailwater.
(a) Pool elevation at top of closed spillway gates where spillway is gated, and at
spillway crest where spillway is ungated. (b) Minimum tailwater. (c) Uplift. (d) Ice
and silt pressure, if applicable.
(a) Pool at standard project flood (SPF). (b) Gates at appropriate flood-control
openings and tailwater at flood elevation. (c) Tailwater pressure. (d) Uplift. (e)
Silt, if applicable. (f) No ice pressure.
(4) Load Condition No. 4 - extreme loading condition - construction with operating
basis earthquake (OBE).
(5) Load Condition No. 5 - unusual loading condition - normal operating with
operating basis earthquake.
(6) Load Condition No. 6 - extreme loading condition - normal operating with
maximum credible earthquake.
(7) Load Condition No. 7 - extreme loading condition - probable maximum flood.
(a) Pool at probable maximum flood (PMF). (b) All gates open and tailwater at
flood elevation. (c) Uplift. (d) Tailwater pressure. (e) Silt, if applicable. (f) No ice
pressure.
b. In Load Condition Nos. 5 and 6, the selected pool elevation should be the one
judged likely to exist coincident with the selected design earthquake event. This
means that the pool level occurs, on the average, relatively frequently during the
course of the year.
45
Stability Considerations
a. General requirements. The basic stability requirements for a gravity dam for all
conditions of loading are:
(1) That it be safe against overturning at any horizontal plane within the structure,
at the base, or at a plane below the base.
(3) That the allowable unit stresses in the concrete or in the foundation material
shall not be exceeded.
Characteristic locations within the dam in which a stability criteria check should
be considered include planes where there are dam section changes and high
concentrated loads. Large galleries and openings within the structure and
upstream and downstream slope transitions are specific areas for consideration.
b. Stability criteria. The stability criteria for concrete gravity dams for each load
condition are listed in Table G.D -1.
Overturning Stability
a. Resultant location.
The overturning stability is calculated by applying all the vertical forces (SV) and
lateral forces for each loading condition to the dam and, then, summing moments
(SM) caused by the consequent forces about the downstream toe. The resultant
location along the base is:
Resultant location =
∑M
∑V
b. Criteria. When the resultant of all forces acting above any horizontal plane
through a dam intersects that plane outside the middle third, a noncompression
zone will result.
For usual loading conditions, it is generally required that the resultant along the
plane of study remain within the middle third to maintain compressive stresses in
the concrete. For unusual loading conditions, the resultant must remain within the
middle half of the base. For the extreme load conditions, the resultant must
remain sufficiently within the base to assure that base pressures are within
prescribed limits.
46
Sliding Stability
τ F ( σ tan φ + c ) 1 ( w − u ) tan φ CA
F .S . = = = +
τ τ P Fφ Fc
Table G.D.2: Partial factors of safety against sliding (taken from Textbook of
Water and Power Engineering, R.K. Sharma & T.K. Sharma)
47
x
h P
P e
W
h/3 R
B B
2 2
Bending moment MB
W 6.W .e W 6.e
+= = 1 +
A A.B A B
where W = weight of dam, A = area of dam section, e = eccentricity
The elementary profile of the gravity dam to bear only the water pressure, will be
triangular in section as shown below. The width of the profile will be zero at the
water surface, where the pressure is zero and maximum at the base, where the
pressure is maximum. Thus the shape of the elementary profile is the same as
that of the hydrostatic pressure distribution. When the reservoir is empty, the
only force acting is self weight (W) of the dam acting at a distance B/3 from the
heel. It is the maximum possible inner-most position of the resultant so that no
tension develops and provides the maximum possible stabilizing force against
overturning without causing tension at toe under empty dam condition. If any
triangular profile other than the right-angled is provided, its weight will act closer
to the upstream face to provide a higher stabilizing force but will cause tension to
develop at the toe.
Vertical stresses developed when dam is empty will be:
W 6.e
Pmax = 1 + ...at..toe
A B
W 6.e
Pmin = 1 − ....at..heel
A B
Now when reservoir is full and downstream empty, forces acting on elementary
profile will be:
48
h
P
P
h
3 W R K
B B
J 3 3
γh
u Uplift pressure
a) Stress criteria. When the reservoir is empty, there is no tension in the dam,
the resultant is acting at the inner 1/3rd point J. When the reservoir is full, for no
tension, the resultan must pass at outer 1/3rd point K.
Now taking moment of all forces about point K,
h2 + c.B2 – B2 sg = 0 or
49
By considering the force triangle, using similar triangles, we have:
(W –u )/P = (h/3)/(B/3)
1 1
Bhsgγ w − γ .c.B.h h
2 2 =
1 2 B
γh
or 2
B ( sg − c ) = h
2 2
h
B=
( sg − c )
When the uplift force is not considered, c = 0 or B = h/ (sg)1/2
P = μ(W – u)
or 1/2γw.h2 = μ(1/2.B.h.sg. γw – 1/2. c.B. γw.h
B = h/μ(sg –c)
B = h/μsg ----------------------------------------------- ©
The base width ‘B” of the elementary profile should be greater of the widths
obtained in equation (B) or (C)
W − u 2(W − u )
p= (1 + 1) =
B B
2
p = [1 / 2 B.h.sgγ w − 1 / 2c.B.γ w .h] = γ w .h.( sg − c ) − − − − − −( E )
B
The corresponding stress at the heel is:
50
W − u
p= (1 − 1) = 0 − − − − − − − − − − − − − − − −( F )
B
When the reservoir is empty, the only force acting on the elementary profile is its
weight, acting through J. In this case, the maximum compressive stress at the
heel = W/B(1+1) = 2W/B and the corresponding normal stress at toe is
W/B(1-1) =0
The elementary profile of a gravity dam is a triangle with maximum water surface
at its apex. This profile is only theoretical one. For meeting the practoical
requirements certain changes have to be made namely: i) for communication,
road has to be provided and therefore a top width; ii) for wave action, free board
above the high flood level must be provided.
The addition of the above will cause the resultant force to shift towards the heel.
Earlier, when the reservoir was empty, the resultant was passing through the
inner middle third. The above changes will shift it towards the heel, crossing the
inner middle third point; this will create tension in the toe. To prevent this tension,
some concrete is added in the dam body towards upstream side.
Free board
h
Dam
a
16
Figure G.D. : Practical dam profile
51
Principal and Shear Stress in Dam
Need to redraw
Let dr, ds and db be the lengths of AB, AC and BC respectively, and let
p = be intensity of water pressure
σ1 = principal stress on plane AB
τ = shear stress
Now considering unit length of the element ABC of the dam, the normal forces on
the planes AB, AC and BC are σ1.dr σ1.dr, pv.db and p.ds sinα respectively.
Resolving all the forces in the vertical direction,
Or pv = p.sin 2 α + σ1.cos 2 α
52
σ1 = p v. sec 2 α –p.tan 2 α ---------------------- (H)
σ1 = p v. sec 2 α -----------------------------(I)
In the same way,considering the hydrodynamic pressure for the upstream side,
in the horizontal direction;
Τ.db = σ1.dr.sin α + p.ds cos α
τ = (σ1-p).sin α.cos α.
or τ = (p v - p) tan α
The shear stress for the upstream side has the same value but with reversed
direction
For the upstream side, τ = - (p v - p) tan α
53
ii) Analytical method.
We shall consider only the analytical method.
For the analytical method, the following steps are carried out:
a) Considering unit length of the dam, all vertical loads are determined the
algebraic sum of all vertical forces ∑V is calculated
b) Considering unit length of the dam, all horizontal forces are determined
and their algebraic sum ∑H calculated.
c) The sum of the righting moments (∑MR) and the sum of overturning
moments ∑MO) at the toe of the dam are calculated.
The difference between the algebraic sum of the overturning and righting
moments is determined i.e ∑M = ∑MR - ∑MO
d) The location of the resultant force FR from the toe of the dam is also
calculated by the following:
xe =
∑ M − − − − − − − − − − − − − − − (J )
∑V
e) Now the eccentricity of the resultant force is determined by:
B
e = − x e − − − − − − − − − − − − − − − −( K )
2
N σT =
∑ V 1 + 6.e − − − − − − − − − − − − − ( L)
B B
N σH =
∑V 1 − 6.e − − − − − − − − − − − − − − − (M )
B B
h) The principal and shear stresses at the toe and heel are determined from
σ1 = p v. sec 2 α -----------------------------(I)
τ = - (p v - p) tan α
54
( F .S .) =
∑M R
∑M O
µ ∑V
Sliding . factor − ( F .S ) =
∑H
µ ∑ V + B.q
Shear. friction. factor − ( F .S ) =
∑H
LOW AND HIGH GRAVITY DAMS
Low Dam: is one of limiting height such that the resultant of all forces passes
through the middle third and the maximum compressive stress at the toe does
not exceed the permissible limit i.e
σ1 =γwH(s-c+1) = fc
f
H= − − − or
γ w ( s.g − c + 1)
or limiting height
f
H=
γ w ( s.g + 1)
The limiting height, for the usual stress of dam material i. e. γw = 9.81 x1000
kg/m3, ρ = 2.4 and f = 30 kg/cm2 (or 300tonnes /m2) gives:
300 x1000
H= = 8.99m
9.81x1000( 2.4 + 1)
s = Specific gravity of material; c = coefficient of uplift pressure, γw = specific
weight of water; f = permissible compressive stress of material.
High Dam: A dam the height of which exceeds the limiting height of low dams is
termed as high. In a high dam, the allowable stress are often exceeded if the
resultant of all the forces were to pass through the middle third; to avoid
excessive stresses the resultant is maintained still near the centre of the base for
which purpose the downstream slope is flattened and the upstream slope is also
provided with a batter.
55
Reservoir level
Limiting height
of low dam
H1
Height of
high dam
Upstream B
batter
Before starting, one must establish whether it is low or high dam using the above
relationship.
.
O
h1 = 2a S s − c
h K
Q
J S
h
B1 =
S −c
s
B
Top width is chosen according to Creager must be about 14%of the dam’s
height.
56
Free Board: is provided on the basis of height of waves and other practical
considerations. In practice a free board of between (3 – 5)% of dam height is
used
h
B1 = − − − or
Ss − c
Base width B1 I is given by
h
B1 =
µ ( S s − c)
WORKED EXAMPLES
The diagram below shows the cross section of a masonry dam. Determine the
stability of the dam. Also determine the principal stress at the toe and heel of the
dam. Take unit weight of dam material as 2250 kg/m3, density of water 1000
kg/m3 and the permissible shear stress of joint = 15 kg/cm2. Assume value of
coefficient of friction μ = 0.75
SOLUTION
A. Vertical forces
i) Self weight of dam = (12 + 2.5)/2 x 15 x 1 x 2250 = 244,687.5 kg
A B
D
β
15m α
D
C
1 2.5 8.5
57
iv) Therefore ∑V = 244687.5 + 7500 – 90,000 = 162.187.5 kg
xi) Factor of safety against sliding = μ∑V = 0.75 x 162,187.5 = 1.08 >1.0 safe
∑H 112,500
xii) Shear friction factor =μ∑V + b.q = 0.75x 162,187.5 + 12 x 15 x 104 = 17.08
∑H 112,500
Stress calculation
58
B 12
Example 2: From the data given below, design a stone masonry gravity dam of
practical profile.
59
4.5
1157.0 M N
1155.5 1.5
L
14m O
25m K
22m
Q DAM
R
1130.5
J S
P
B1 = 16m
0.3
The design of the dam can be done with respect to the details.
Depth of water = 1155.5 – 1130.5 = 25m
Top width of dam = 14% height of dam = 14?100 x 26.5 = 3.71m
Assume a roadway width of 4.5 m
Therefore provide top width MN of dam = 4.5m
h 25
Base width of dam PS = B1 = = = 15.81m ; assume 16m
( Ss − c ) 2.5 − 0
Extra width JS = MN/16 = 4.5/16 = 0.281m (assume 0.3m)
Vertical distance LK = 2a Ss = 2 x 4.5 2.5 = 14.23 assume 14m
Vertical distance LR = 3.1a Ss = 3.1x 4.5 2.5 = 22m
Example 3
What should be the maximum height of elementary profile of a dam, it the safe
limit os stress on the masonry should not exceed 350 tonnes/m2. Assume weight
of masonry 2.4 tonnes/m3. Determine the base width also. Determine H and B if
uplift intensity factor is 0.67 and factor of safety is 2
Solution: The limiting height of elementary profile of a masonry dam
60
fc 350 x1000 xg
H= = = 103m
γ w ( sg + 1) 1000 xg (2.4 + 1)
H 103
Base width B = = = 66m
Ss − c 2.4
fc 350 x1000 xg
ii) F.S =2, c = 0.67: H = = = 64m
F .Sxγ w ( Ss − c + 1) 2.0 x1000 xg ( 2.4 − 0.67 + 1)
H 64
B= = = 49m
Ss − c 2.4 − 0.67
Example 4.
A concrete gravity dam has maximum water level 305.0m, bed level 225.0m, top
required level of dam 309.0m, downstream face slope starts at required level
300.0m, downstream slope 2:3, tail water is nil, upstream face of dam is vertical,
centre line of drainage gallery is 8m downstream of upstream face, uplift
pressure is 100% at heel, 50% at line of gallery and zero at toe, specific gravity
of concrete is 2.4. Considering only weight, water pressure and uplift, determine
i) Maximum vertical stresses at toe and heel of dam
ii) major principal stresses at toe of dam and
iii) Intensity of shear stress on a horizontal plane near the toe
Solution:
Height of dam H = 309.0 – 225.0 = 84m
Depth of water h = 305.0 – 225.0 = 80m
1. Top width of dam = 14% of height = 0.14 x 84 = 12m
2. Bottom width of dam = 12 + (300 – 225) x2/3 = 62m
Calculating weight and moments by considering unit length of dam
61
309
305 8 54
300
W1
U1
α u2
40
80
W2
80
u1
225
62
62
Principal stress at toe σ = p n sec 2 α – p’ tan 2α p’ = zero because tail wate is
zero; tan α = 2/3 sec 2 α = 1 + tan 2 α = 1 +(2/3)2 =13/9
Therefore σ = p n sec 2 α= 138 γw x13/9=199 γw N/m2
Example 5. A concrete gravity dam has maximum reservoir level 150.0m, base
level of dam = 100.0m, tail water elevation 110.0m, base width of dam 40m,
location of drainage gallery 10m from upstream face which may be assumed as
vertical. Compute the hydrostatic thrust and the uplift force per metre length of
dam at its base level. Assume 50 % reduction in net seepage head at the
location of the drainage gallery.
Solution:
Free board = 5% of dam height = 0.05x50= 2.5 adopt 3m
Dam height h = 50 + 3 =53m
Top width a = 14%xh = 0.14 x 53 = 7.5m adopt 7.5m
1 1 30
1
143
W1
α U3
u2
30
U4
50
W2 50
u1
1
40
Designation x γw Force x γw Moment arm Moment about toe
X γw
63
Weight of dam
W1=7.5x53x1x2.4 =954 32.5+7.5/2=36.25 34583
W2 = 32.5 x43/2 x1 x2.4 =1677 32.5x2/3=21.67 36341
Weight of tail water
W3 = 8x10/2x1x1 40 8 x 1/3 = 2.67 107
∑V1 =2671 ∑M1 (+) 71031
Uplift pressure
U1=20 x 10/2 100 30 + 10 x2/3=36.67 3667
U2=10 x30 x1 300 (30 +10/2) = 35 10500
U3=10 x30 x 1 300 (30/2) = 15 4500
U4=20 x 30/2 x 1 300 (30 x 2/3) = 20 6000
∑V2= 1000 ∑M2 =24,667
Water pressure
U/S=γwh2/2=γwx1.10.10/2 1250 50/3 (-) 20,833
D/S=γwh2/2=γw.1.10.10/2 50 10/3 (+)167
∑H3= 1200 ∑M3 (-) 20,666
64
Chapter 4
Spillway is a passageway to convey past the dam flood flows that cannot be contained in
the allotted storage space or which are in excess of those turned into the diversion
systems. Spillways function infrequently, at times of flood or sustained high runoff, when
other facilities are inadequate. However its ample capacity is of prime importance for the
safety of the dam and other hydraulic structures. Hydraulic aspects of spillway design
relates to design of the three spillway components: control structure, discharge channel
and terminal structure. The control structure regulates outflows from the reservoir and
may consist of a sill, weir section, orifice, tube or pipe. Design problems here relate to
determining the shape of the section and computing discharge through the section. The
flow released from the control structure is conveyed to the streambed below the dam in a
discharge channel. This can be the downstream face of the overflow section, a tunnel
excavated through an abutment, or an open channel along the ground surface. The
channel dimensions are fixed by the hydraulics of channel flow. Terminal structures are
energy dissipating devices that are provided to return the flow into the river channel
without serious scour or erosion at the toe of the dam.
Spillways are usually located at the normal flood level and the water just overflows the
crest when it gets to that level. However, some spillways are equipped with gates, which
are temporary barriers installed over the permanent crest of the spillway for storing
additional water during period of low water seasons. All small flows exceeding the
barrier top level are allowed to pass over the barrier, but during large flood flows, the
barrier is removed and full spillway capacity is used to discharge the flood water.
Types of Spillways
Spillways may be classified, depending upon the type of structure as:
i) Side channel spillway
ii) Straight drop spillway
iii) Overflow or Ogee spillway
iv) Chute or Trough spillway
v) Shaft (Morning glory) spillway
vi) Siphon spillway
Straight Drop (free overfall) Spillway. In this type, water is allowed to fall freely from a
low weir and vertical fall structure. In some cases the crest of the spillway is extended in
65
the shape of overhanging lip which keeps the discharge away from the straight drop
section. Ogee (Overflow) Spillway: is a special form of a weir whose shape is made to
conform to the profile of lower nappe of a ventilated sheet of water falling from a sharp-
Upstream
Earthfill dam
downstream
crested weir. The profile is so shaped that the discharging water always remains in touch
with the spillway surface.
66
Spillway
spillway
x 1.85 = 2 y ( h )
0.85
--------------------------- (S-1)
where x and y are the co-ordinates of the crest profile measured from the apex of the
crest, and h is the design water head.
Design head h
x
O
a = 0.175h
r2 b = 0.292h
r1
r1 = 0.50h
r2 = 0.20h
67
0.724( x + 0.27h )
1.85
y= ------------------------ (S-2)
( h ) 0.85
For larger discharges, (eg flows beyond H.F.L), the nappe may leave the ogee profile and
will cause negative pressure resulting in cavitation and increase of discharge.
Q = C.Le.(he)3/2 ------------------------------(S-3)
68
Example 9.10
Design an overflow spillway section for a design discharge of 50,000cfs. The upstream
water surface is at Elv. 800 and the channel floor is 680. The spillway, having a vertical
face, is 180ft long.
Solution
1. Assume C = 3.95
Q 50,000
he3 / 2 = = = 70.32
2. from the discharge equation CL 3.95(180 )
he = 17.1
3. Depth of water upstream = 800 – 680 = 120 ft.
Velocity of approach vo = 50,000/120(180) = 2.31ft/sec
Velocity head = v2/2g = (2.31)2/(32.2) = 0.08 ft
4. Maximum water head = 17.1 – 0.08 = 17.0ft
5. Height of crest, P = 120 – 17.0 = 103ft.
6. Since he <30 ft, design head, hd = 17.1/1.42 = 12.0ft
7. P/hd = 103/12 = 8.58 > 1.33, high overflow section
8. Downstream quadrant of the crest shape
X (select)ft Y (computed) ft
5 1.18
10 4.25
15 9.00
20 15.30
30 32.40
X (selected) ft Y (computed) ft
1.0 0.09
2.0 0.39
3.0 1.10
3,36 2.00
69
CHUTE OR TROUGH SPILLWAY
Chute spillway is a type of spillway in which the discharge is conveyed fro a reservoir to
the downstream river level through a steep open channel placed either along the dam
abutment or through a saddle. In a chute spillway, the velocity of flow ia always greater
than the critical. The name chute applies, regardless of the control device used at the
head of the chute, which can be an overflow crest, a gated orifice or a side channel crest.
The chute spillway consists of four parts: an entrance channel, a control structure or crest,
the sloping chute and a terminal structure. The entrance structure is usually an open
channel of sub-critical flow. The critical velocity occurs when the water passes over the
control.
Flow in the chute is maintained at supercritical stage until the terminal structure DE is
reached. It is desirable that from B to C, where a heavy cut is involved, the chute is
placed on a light slope. From C to D it follows the steep slope and ends with an energy
dissipating device placed at the bottom of the valley D. The axis of the chute is kept as
straight as possible; otherwise, the floor has to be super-elevated to avoid the piling up
high-velocity flow around the curvature.
It is preferable that the width of the control section, the chute, and the stilling basin are
the same.
To prevent hydrostatic uplift under the chute, a cutoff wall is provided under the control
structure and a drainage system of filters and pipes is provided. When the stilling basin is
operating, there is substantial uplift under the lower part of the chute and upstream part of
the stilling basin floor. The floor must be made sufficiently heavy or anchored to the
foundation.
Chute section
A
Terminal Structure
B
C
70
It is important that the slope of the chute in the upstream section BC chute, should be
sufficiently steep to maintain a supercritical flow to avoid the formation of a hydraulic
jump. Therefore slope should always be above the critical slope.
Fro the Manning’s formula;
1
Q= AR 2 / 3 I 1.2 -------------------------- ( S-5 )
n
h3 = q2/g ----------------------------------(S-6 )
By putting this into Manning’s formula, the critical slope could be determined. The slope
to be given to the chute must be greater than this critical slope.
A review of existing spillways indicate that the actual slopes of the upstream section of
the chute are 1 or 2 % or more.
Chute Sidewalls
Except for the diverging or converging sections, chute channels are designed with parallel
vertical sidewalls, commonly of reinforced concrete 30 to 45cm thick; designed as
retaining wall. The height of the walls is designed to contain the depth of flow for the
spillway design flood by gradually flows hydraulics equation. Free board or allowance is
made for pier end waves, roll waves and air entrainment. In view of uncertainties
involved in the evaluation of surface roughness, pier end waves, roll waves and air
entrainment, the following empirical equation is added to the computed depth of water
surface profile
In side-channel spillway, the overflow weir is placed along the side of the discharging
channel, so that the flow over the crest falls into a narrow channel section (trough)
opposite the weir, turns through a right-angle, and continues I the direction approximately
parallel to the weir crest
71
Side channel trough
A
Side channel crest
A
Dam
Control section
Trough
Chute
Side-channel spillway: plan and section
This type of weir is adoptable to certain special conditions, such as when a long overflow
crest is desire but the valley is narrow, or where the overflow are most economically
passed through a deep narrow channel or a tunnel. The crest is similar to an overflow or
ordinary weir section. A control section downstream from the trough is achieved by
constricting the channel or elevating the channel bottom to induce critical flow.
Downstream from the control section functions as a chute spillway. Thus the side-
channel design is concerned only with the hydraulic action in the trough upstream of the
control section, where varied flow takes place. The water surface profile in the trough is
determined from the momentum equation applied to gradual varied flow. If flow rate and
velocity in left is Q1 and V1 while those at the right is Q2 and V2, then by applying the
momentum equation, the rate of change of momentum in the reach dx is equal to the
external forces acting in the reach.
Rate of momentum at left(upstream) = ρQ1V1
Rate of momentum at downstream = ρQ2V2
Change in momentum = ρ(Q2V2 – Q1V1) ------------------------------------------(S -8)
Due to smallness of friction forces and weight component in the direction of flow, the
only force acting is the hydrostatic pressure force.
The resultant hydrostatic force Fp1 –Fp2 = γ.(A1. h1 – A2h2) = γ(Q1/V1 – Q2V2)dy ---(S-9)
72
Crest
dy
P1
P2
dx
Q2 V1 V2 V ( Q − Q1 )
dy = + (V1 − V2 ) + 1 2 ---------------------------(S-10)
g Q1 Q2 Q2
Worked example:
Design a side-channel trough for a spillway of 100ft length for a maximum discharge of
2500 cfs. The side-channel trough has a length of 100ft and a bottom slope of 1ft in 100
ft. A control section of 10 ft width is placed downstream from the trough with the bottom
of the control at the same elevation as the bottom of the trough floor at the downstream
end.
Solution
1/. Critical depth at the control section, yc = (q12/g)
2. q1 = 2500/10 = 250cfs/ft.
3. yc = (2502/32.2)1/3 = 12.44ft
6. For a side-channeltrough, assume a trapezoidal section with 2V:1H slope and a 10ft
bottom width. Also, assume that the transition loss from the end of side-channel trough to
the control section is equal to 0.2 of the difference in velocity heads between the end of
the transition.
73
7. The following energy equation may be written between the trough end and the control
section;
y100 + h100 = yc + hc + 0.2(hc – h100) or
y100 + 1.2h100 = 12.44 + 1.2(6.27) = 19.96 --------- (S-11)
where the subscript c refers to “critical” and “100” refers to the distance of the trough
from the upstream end of the spillway. Equation (S-11) is solved by trial and error.
Assume y100 = 19.2; then A = 376.3 ft2, v = Q/A = 2 500/376.3 = 6.64 ft/sec, and h100 =
6.64)2/2(32.2)0.68. Thus Eq. (S-11) is satisfied.
8. The known values above relate to section 1 at the downstream end of the trough.
Section 2 is taken at the upper end of a selected increment dx, 25 ft in this case. A value
of the change in water level, dy, is assumed for each reach, all terms of equation (S-10)
are evaluated, and dy is computed, as in table below. The assumed and computed values
of dy should match; otherwise, a new value is assumed for dy.
9. The process is repeated until the upstream end of the channel is reached.
1 2 3 4 5 6 7 8
dx Bottom dy Water y A Q= v
(selected Level (Assume) Level (col 4-col 2) q(L-dx) = Q/A
)
D/S end 100.0 119.2 19.20 376.3 2500 6.64
25 100.25 1.0 120.2 19.95 398.5 1875 4.71
0.62 119.82 19.57 387.2 4.84
25 100.50 0.50 120.32 19.82 394.6 1250 3.17
0.40 120.22 19.72 391.6 3.19
25 100.75 0.25 120.47 19.72 391.6 625 1.60
15 100.90 0.10 120.57 19.67 390.2 250 0.64
0.07 120.54 19.64 389.3 0.64
9 10 11 12 13 14
Q1 + Q2 Q2 – Q1 v1 + v2 v2 – v1 dy Remarks on
computed assumed dy
4375 -625 11.35 -1.93 0.63 High
4325 -625 11.48 -1.80 0.61 Ok
3125 -625 8.01 -1.67 0.41 High
3125 -625 8.03 -1.65 0.41 Ok
1875 -625 4.79 -1.59 0.24 Ok
875 -375 2.24 -0.96 0.07 High
875 -375 2.24 -0.96 0.07 ok
74
6. Q1 + Q 2 (col. 9) = Col. 7 + value in col. 7 at the preceding station (section)
MORNING GLORY OR SHAFT SPILLWAYS
In this type of spillway, water enters over a horizontal circular crest and then drops
through a circular shaft, and then through a horizontal tunnel or conduit. This type of
spillway is suitable under the following situations:
1. For damsites with steeply rising abutments particularly where a diversion tunnel
can be utilized as discharge carrier.
2. For sites in narrow canyons.
3. For sites where there is inadequate space for locating other type of spillway.
Advantages: The main advantage of such type of spillways are as follows:
1. The nearly maximum capacity may be attained at relatively low heads.
2. It is ideal for sites where maximum spillway overflow is to be limited.
This type of spillway consists of four parts: 1) a circular weir at the entry, 2) a flared
transition conforming to the shaft of the lower nappe of a sharp-crested weir, 3) a vertical
drop shaft, and 4) a horizontal or near-horizontal outlet conduit or tunnel. As the head
increases, the control shifts from weir crest, to drop shaft, and to outlet conduit.
Condition 1: a free-discharging weir prevail as long as the nappe forms to converge into
the shape of a solid jet.
Condition 2: weir crest is drowned out. The US Bureau of Reclamation (1977) indicated
that this condition is approached when Hd/Rs >1, where Hd is the design head and Rs is
the radius of the crest.
Condition 3: Spillway is flooded out, showing only a slight depression and eddy at the
surface. Under condition 3, the head rises rapidly for a small increase in discharge. Thus
the design is not recommended under this condition (i.e., under the design head, the outlet
conduit should not flow more than 75% full).
The U.S. Bureau of Reclamation suggested that the following weir formula may be used
for flow through the shaft spillway entrance regardless of the submergence, by adjusting
the coefficient to reflect the flow conditions.
75
Where C = discharge coefficient related to Hd/Rs and P/Rs , where Hd = design head and
P = crest height from the outlet pipe; Rs = radius to the circular crest; H = head over the
weir ( Gupta page 495)
Equation (S-15) can be used to determine the crest size (radius), Rs for a given discharge
under the maximum head.
Equation (S-13) can be used to determine the shape of the transition (drop shaft) that is
required to pass the design discharge with the maximum head over the crest.
Worked Example: A shaft spillway is to discharge 2000 cfs under a design head of 10 ft.
Determine the minimum size of the overflow crest. Also determine the shape of the
transition if the control section is 4 ft below the crest level.
Solution
1. Since the coefficient C is related to P and Rs, assume that P/Rs > 2 and determine
Rs by trial and error.
2. Try Rs = 7 ft.
Hd/Rs = 10/7 = 1.43
From the curve of C as a function of Hd/Rs, we choose C = 1.44
3. From equation (S-15)
or R = 9.14/(Ha)1/4
Ha (select), ft R (ft)
14 4.73
16 4.57
18 4.44
20 4.32 SIPHON SPILLWAYS
A siphon spillway is a short enclosed duct whose longitudinal section is curved. When
flowing full, the highest point in the spillway lies above the liquid level in the upstream
reservoir, and the pressure at that point must therefore be sub-atmospheric. This is the
essential characteristics of a siphon. Siphon spillways can be saddle type or volute type
but the volute is not very common. The siphon has usually three parts, 1) the inlet or
mouth; 2) throat, and 3) lower limb.
The siphon spillway functions as follows:
When the water level exceeds the crest level, the water commences to spill and flows
over the downstream slope in much the same way as a simple Ogee spillway. As the
water rises further, the entrance is sealed off fro the atmosphere. Air is initially trapped
76
within the spillway, but the velocity of flow of water tends to entrain the air (giving rise
to aeration of the water) and draws it out through the exit. When all the air has been
expelled, the siphon is primed and is therefore acting as a simple pipe. \There are thus
three possible operating conditions depending on upstream depth.
1. Gravity spillway flow
2. Aerated flow (during self priming)
3. Pipe flow (after priming).
Operational problems with siphon spillways. The aerated condition is unstable and is
maintained while the siphon begins to prime. In a simple siphon, a small change in H
produces a sharp increase or decrease in the discharge through the spillway. Depending
on the discharge entering the reservoir, the siphon could go through the following cycle:
a) if the spillway is initially operating with gravity flow, then the upstream
(reservoir) level must rise,
b) when the upstream level has risen sufficiently, the siphon primes and the
spillway discharge increases substantially,
c) the upstream level falls until the siphon de-primes and its discharge drops.
The cycle (a) to (c) is then repeated. This obviously can give rise to radical surges and
stoppages in the downstream flow.
Other potential problems encountered with siphon spillway are:
i) blockage of spillway entrance by debris ( this problem could be overcome by
submerging the inlet of the hood into the water or installing a trash-
intercepting grid in front of the intake).
ii) Substantial foundations required to resist vibrations during operation of
siphon.
iii) Waves arriving in the reservoir during storms may alternately cover and
uncover the entry, thus interrupting smooth siphon action
The discharge of a saddle siphon can be calculated from the following formula;
Q = CA 2 gH --------------------------- (S-16)
where A = area of cross section at crown; A = L x b
L = length of hood (going into the paper)
b = height of throat
H = operating head
H = Reservoir level – downstream tail water level if outlet is submerged
H = Reservoir level - downstream centre of outlet if the outlet is discharging
freely.
C = Coefficient of discharge. Its average value may be taken as 0.65.
77
Air vent
Crown
De-primer
Throat
Hood
Crest
Mouth
H1
H2
Spillway gates are the temporary barrier installed over the permanent crest of the
spillway, for storing additional water during dry weather season. The small flows in
excess above the spillway gates is allowed to pass over the gates, but in case of large
flood, the spillway gates are opened and the full capacity is used to remove excessive
flood water.
Spillway gates can be provided on all types of spillways except siphon spillway.
Following types of spillway gates are commonly used:
a) Flash board gates
b) Stop logs or needle gates
c) Radial gates
d) Drum gates
e) vertical lift gates
78
Flash Board Piers
Spillway crest
Stop log
79
Chapter 5
RESERVOIRS
Reservoirs are the most important man-made storage elements in water systems
because their capacity and operational schedules determine the rates and
volumes of flow in streams. They have enabled humans to make the desert
bloom and flourish and provide water supplies for large and concentrated
populations. A reservoir is usually created by the construction of a dam across a
flowing stream. Storing water is necessary because when water occurs naturally
in streams and is not stored, it is sometimes not available when needed.
Reservoirs solve this problem by capturing water when it occurs and making it
available at latter times. However, with today’s environmental awareness,
reservoirs are sometimes negative symbols of “man taking dominion over nature”
and work against sustainable development as they interfere with natural
ecosystems.
While the concept of a reservoir may bring to mind a large body of water, many
small reservoirs are also in service. These include urban water tanks, farm
ponds, regulating lakes, and small industrial or recreational facilities. These
small reservoirs can have important cumulative effects especially in rural areas
A reservoir is created with the impounding of runoff from the catchment upstream
by the construction of a dam across a river or stream. Storage is done during the
Max. Elevation
Dead storage
period the flow is in excess of the demand for release during the lean period.
Storage reservoirs are constructed for the purposes of 1) flood control or 2)
conservation of water. Reservoir construction may be contemplated under one of
the following conditions:
80
i) Running water is available in sufficient quantity, but running to waste or causing
damage. In this case, an attempt is made to utilize the available water for
developmental purposes of the area.
ii) Water is needed all the time either for the generation of hydro-power to meet
the requirement of power for development or for irrigation , however the flow in
the stream is unable to meet the demand all the time.
a) Availability of water
81
Investigations for Reservoir Planning
The contour interval is 2.5 m but at the site of dam an interval of 1.5 m is used
and a map on a scale of 1/250 to 1/500.
From the contour map of the reservoir area, the water spread of the reservoir at
any elevation may be directly determined by measuring the area at that contour
with a planimeter. From these measurements, an area – elevation curve can be
drawn. The capacity of the reservoir may be determined by taking contour areas
at equal interval and summing up these areas by any of the following methods:
Elevation, H(m)
V = f(H)
A = f(H)
Area (Km2)
82
dam, reservoir size, capacity of power house and other installations and 2) to
determine the hydrograph of the worst flood to estimate the spillway size.
The storage capacity of site should be sufficient to the meet the requirement for
which it is designed.
i) Site where the width of river is narrow but rapidly widening upstream
ii) Site where sufficient quantity of water is available. This will depend on the
intensity of rainfall, run-off and catchment area.
iii) As much as possible, site should be such that water can flow under gravity
from reservoir to demand points
iv) Site must be close to demand points
v) The geological conditions at the site should permit minimum percolation
losses, with maximum run-off.
vi)Site must be close to useful construction materials.
vii) The topography should be favourable for site for spillway.
viii) Site should be such that the run-off water has minimum percentage of
sediment
ix) Site must be free from such minerals and salts, which may make the water
unfit for the purpose, for which the dam built.
Classification of Reservoirs
a) Storage reservoir
b) Flood Control (flood mitigation) Reservoir
c) Distribution Reservoir
d) Single Purpose reservoir,
e) Multipurpose Reservoir.
Storage Reservoir: are small storage capacity reservoir, which are constructed
to store water for meeting the demand for water supply, irrigation, power
generation etc.
Flood Control (Flood Mitigation) Reservoir: The main function of such
reservoirs is to temporary store flood waters and release it at such a safe rate
that it may not flood downstream side.
Distribution Reservoir: These are small capacity reservoirs used for water
supply of towns. These reservoirs may store raw water for treatment. As the
water is required at varying rate during the day, they permit the water treatment
plant to function at a constant rate
83
Single Purpose Reservoir: It is a reservoir to serve only one purpose which
may be: a) Municipal Water Supply, b) Power generation; c) Flood Control;
d) irrigation; e) Navigation; f) Recreation.
Multipurpose Reservoir: In this reservoir, the storage and release cater for a
combination of two or more purposes.
STORAGE ZONES AND LEVELS OF A RESERVOIR
1. Dead Storage: us about 10 – 25 % of the gross storage and it is provided to
cater for sediment deposition by the impounded sediment –laden waters. Usually
this volume is below sluiceway, and therefore below it, the reservoir is not
susceptible to release water by the built in outlet means. The dead storage is
equivalent to the volume of sediment expected to be deposited in the reservoir
during the design life of the reservoir.
2. Live (Active or Useful) Storage: It is the storage capacity above the inactive
storage, which constitutes useable portion of the total storage.It is thus the
difference of gross storage capacity and the sum of dead storage capacity and
inactive capacity Live storage has to be sufficient so that the project is successful
for i) 75% of its life period in an irrigation project; ii) 90% for hydro- power and iii)
100% for water supply.
3. Flood Storage: is the storage contained between the normal reservoir level
and the full reservoir level. It is the storage space provided in a reservoir for
storing flood water temporary to moderate the releases downstream.
4. Valley Storage: is the storage in the river in floods after it gas overflowed its
banks. It is important element in the design of large size flood control reservoirs
where it could be a significant proportion in respect to the reservoir storage
volume.
5. Surcharge Storage: It is the storage between the Normal full reservoir level
and the maximum possible level in the reservoir. This storage is usually difficult to
control since it depends on maximum floods, rains and the resulting run-offs.
6. High Flood Level (Maximum Pool Level): The level to which the water will rise
during the design flood
7. Full Reservoir Level The level to which water will rise during ordinary
conditions of operation of reservoir.
Minimum Pool Level: The level to which water from the reservoir in ordinary
conditions may be drawn.
Storage Capacity and Yield
Yield: It is the volume of water that can be supplied from the reservoir I a specific
interval of time. The time interval may vary from a day to a year.
Dependable or Firm Yield: The maximum guaranteed supply of water during the
worst dry period.
84
Secondary Yield: The quantity of water available in excess of dependable yield
during the flood period.
Designed Yield: It is the quantity of water for which the project is designed after
ascertaining the availability of water. It is usually kept lower than firm yield.
Average Yield: It is the arithmetic average of safe yield and secondary yield over
a long period.
Reservoir Storage Capacity Determination
The storage capacity of a reservoir to meet the demand of continuous supply is
determined with the help of observed discharge data of a stream on which the
dam is to be constructed. The flow values for the driest years in as long a period
as is available, eg, 25 to 30 years are used.
There are two approaches for the determination of the size of reservoir storage.
The simplified method, which are commonly used in planning stage studies,
comprise mass curve analysis.
The detail method, used at the time of developing reservoir operating plans,
performs a sequential reservoir routing of the historical flows, which will not
be treated here.
A third, by Bar Graph Method. In this method, the inflows of the driest year are
plotted as ordinates against time as abscissa to get a stepped graph as bar
graph. The area under the bar graph represents total volume of inflow into
reservoir. The average demand is likewise plotted. The area between the two
plots indicates either surplus or deficit volume. The area of maximum deficit
between the demand and inflow represents the minimum storage required.
Surplus
Deficit
Surplus
Max. deficit
STORAGE
85
Simplified Procedure for Reservoir Storage Capacity.
There are two methods of analyses: i) the sequential mass-curve method and ii)
the non-sequential mass-curve method.
Mass curves are useful in reservoir studies since they provide a ready means of
determining storage capacity necessary for a particular average drawoff.
difference between the inflow and outflow to the reservoir during successive
periods are evaluated, the maximum value of which is the required storage.
86
demand volume is determined. The maximum value of cumulative excess of
demand over inflow represents the minimum storage necessary to meet the
demand. The cumulative excess inflow volume starting from each demand
withdrawal from the storage is also determined which indicates the filling up of
the reservoir and the volume in excess of storage to be spilled over.
Example: In the table below, (col. 2) gices the yield of a river, (col. 3 & 4) gives
total losses including evaporation, percolation, etc., and total estimated
consumption of water. Determine the storage capacity of the proposed reservoir.
1 2 3 4 5 6 7
Month Yield of Total Total Total Surplus Total
River Losses consumptio surplus
(3 + 4) (2 – 5)
n
June 0.0385 0.0863 0.1883 0.2746 - -
July 0.4013 0.0784 0.1851 0.2633 0.1378 -
Aug 0.4327 0.0678 0.1782 0.2460 0.1863 -
Sept. 0.4832 0.0615 0.1513 0.2128 0.2704 -
Oct 0.3384 0.0482 0.1482 0.1964 0.2902 -
Nov. 0.3156 0.0257 0.1215 0.1472 0.1684 -
Dec. 0.2983 0.0173 0.0987 0.1160 0.1823 -
Jan 0.2357 0.0195 0.1142 0.1337 0.1020 -
Feb. 0.2013 0.0202 0.1348 0.1550 0.0463 1.3837
Mar. 0.1605 0.0823 0.1562 0.2085 - -
April 0.1114 0.0845 0.1732 0.2577 - -
May 0.0713 0.1032 0.1784 0.2816 - -
June 0.0401 0.0975 0.1895 0.2870 - -
July 0.3247 0.0767 0.1883 0.2650 0.0597 -
Aug 0.3652 0.0681 0.1790 0.2471 0.1181 -
Sep 0.2830 0.0600 0.1551 0.2151 0.0679 -
Oct 0.2218 0.0515 0.1503 0.2018 0.0200 -
Nov 0.1943 0.0301 0.1284 0.1585 0.0358 -
Dec 0.1814 0.0165 0.1010 0.1175 0.0639 -
Jan 0.1703 0.0183 0.1186 0.1369 0.0334 -
Feb 0.1415 0.0197 0.1207 0.1404 0.0011 03997
87
Mar 0.1137 0.0515 0.1531 02046 - -
April 0.0945 0.0726 0.1802 0.2528 - -
May 0.0687 0.0985 0.1854 0.2839 - -
June 0.0312 0.1013 0.1903 0.2916 - -
July 0.4215 0.0695 0.1900 0.2595 0.1620 -
Aug 0.4183 0.0603 0.1805 0.2408 0.1775 -
Sept. 0.3967 0.0587 0.1570 0.2157 0.1810 -
Oct 0.3818 0.0494 0.1525 0.2019 0.1799 -
Nov 0.2323 0.0286 0.1312 0.1598 0.0725 -
Dec 0.1738 0.0108 0.1135 0.1243 0.0495 -
Jan 0.1621 0.0143 0.1218 0.1361 0.0260 0.8484
Feb 0.1358 0.0195 0.1256 0.1451 - -
Mar 0.1148 0.0613 0.1583 0.2196 - -
April 0.0887 0.0688 0.1831 0.2519 - -
May 0.0595 0976 0.1867 0.2843 - -
8 9 10 11 12 13
Deficit Total deficit Cumulative Cumulative Difference Carry over
(5-2) yield Draft (10-11)
0.2361 - 0.0385 0.2746 0.2361 Emptying
- - 0.4398 0.5381 0.0983 of reserv.
- - 0.8725 0.7481 0.0884 Filling of
- = 1.3557 1.9969 0.3588 Reservoir
- - 1.6941 1.1933 0.5008 Overflow
- - 2.0097 1.3405 0.6692 Overflow
- - 2.3080 1.4565 0.8515 Overflow
- - 2.5437 1.5902 0.9535 Overflow
- - 2.7450 1.7452 0.9998 Emptying
0.048 - 2.9055 1.9537 0.9518 Emptying
0.1463 - 3.0169 2.2114 0.8055 Emptying
0.2103 - 3.0882 2.4930 0.5952 Emptying
88
0.2469 0.8876 3.1283 2.7800 0.3483 Emptying
- - 3.4530 3.0450 0.4080 Emptying
- - 3.8182 3.2921 0.5261 Emptying
- - 4.1012 3.5072 0.5940 Emptying
- - 4.3230 3.7090 0.6140 Emptying
- - 4.5173 3.8675 0.6498 Emptying
- - 4.6987 3.9850 0.7137 Emptying
- - 4.8690 4.1219 0.7471 Emptying
- - 5.0105 4.2623 0.7482 Emptying
0.0909 - 5.1242 4.4669 0.6573 Emptying
0.1583 - 5.2187 4.7197 0.4990 Emptying
0.2125 - 5.2874 5.0036 0.2838 Emptying
0.2604 0.7221 5.3185 5.2952 0.0234 Emptying
- - 6.7401 5.5547 0.1854 Filling
- - 6.1584 5.7955 0.3629 Filling
- - 6.5551 6.0112 0.5439 Filling
- - 6.9369 6.3121 0.7238 Filling
- - 7.1692 6.3729 0.7963 Overflow
- - 7.3430 6.4972 0.8458 Overflow
- - 7.5051 6.6333 0.8718 Overflow
0.0093 - 7.6409 6.7784 0.8625 Overflow
0.1048 - 7.7557 6.9980 0.7577 Emptying
0.1632 - 7.8444 7.2499 0.5945 Emptying
0.2248 0.5021 7,9039 7.5342 0.3697 Emptying
Graphical Method: there are two graphical procedures. In the first method,
(Mass Curve method), It is accumulated separately as a mass inflow curve and
Ot as a mass yield curve. For a constant draft rate, the yield curve is a straight
line having a slope equal to the draft rate. At each high point on the mass inflow
curve, a line is drawn parallel to the yield curve and extended until it meets the
inflow curve. The maximum vertical distance between the parallel yield line and
the mass inflow curve represents the required storage. Assuming the reservoir is
full at A; and going from A to F along the inflow curve. From A to B, the draft is
more than inflow resulting in a lowering of reservoir; from B to C, the inflow is
89
higher than the draft, but not enough to refill the reservoir; from C to D, the draft
is more, once again causing a drop in reservoir level; from D to E, however, the
inflow is very high, thus filling the reservoir at E.
Cumulative volume F
Inflow curve
G
C
Storage
D Yield/demand
Curve
A
B
Time
∑(I –O)
Storage volume
Time
90
RESERVOIR SEDIMENTATION
Every stream (river) carries some sediment load along with the flow. The larger
solid particles roll along the bed as bed load. Smaller particles are kept in
suspension by the upwards component of the turbulent forces and will only settle
by gravitational force when velocity reduces. Such particles are called
suspended load.
When sediment laden water reaches a reservoir, in the vicinity of a dam, the
velocity and the turbulence are greatly reduced due to enlarged area of the
channel. The larger suspended particles and most of the bed load gets deposited
in the head reaches of the reservoir. The smaller particles remain in suspension
for longer periods and are deposited further down in the reservoir. Very fine
particles may remain in suspension for still longer period and some may pass
over the dam along with the water discharged through spillway or sluiceways.
If the water stored in the reservoir is clear and the inflow is charged with
sediment or is muddy, the heavier water with sediment will flow along the channel
bottom towards the dam under the influence of gravity and the clear light water
will flow on the upper surface of the turbid water. This condition is known as
stratified flow and the under flow of sediment laden water is known as density
current. The process whereby sediments are deposited at the bed of the
reservoir due to reduced turbulence and velocity is known as reservoir
sedimentation. The deposition of sediments will reduce the water storing capacity
of the reservoir and if this process continues for a long time, a stage will be
reached when the whole reservoir will be silted and would be rendered useless.
91
Estimation of the Rate of Sedimentation
The amount of sediment load carried by a stream is determined by taking
samples of water carrying silt at various depths. The samples are filtered,
sediment removed and dried. The sediment load per day may be computed by:
CxQ
Wd = 86400 x tones per day
1000
where Wd = dry weight of sediment
Q = mean daily discharge in cumecs
C = dry sediment weight in grams per litre of water.
Khosla’s Studies.
In the absence of any information on sediments, one could use the studies of Dr.
A. N. Khosla. On the basis of studies on some catchment areas, Khosla
suggested the rate of silting per 100 square kilometer of catchment to be 0.036
million cubic metres.
Recent studies have shown that Khosla’s value is on the low side. The average
sediment rate of reservoirs in USA is given below:
92
Trap efficiency (η) = Sediment deposited in reservoir .
Total sediment brought into reservoir by river
It has been observed that most reservoirs trap about 95 – 100% sediment
brought to it. It is not possible to reduce it below 90 % irrespective of sediment
control measures taken.
Capacity-Inflow Ratio
Capacity inflow ratio is the ratio of the capacity (volume) of a reservoir to the total
inflow (volume) of water in a given time (usually one year).It is an important factor
as detailed observations of reservoir sedimentation has shown that trap
efficiency is a function of this factor.
Thus: η
trap efficiency(η) = f (Reservoir capacity)
Total inflow
The figure below shows a graph of trap efficiency and capacity inflow ratio based
on observations of existing reservoirs. From the figure, it will be noted that for a
given inflow rate, the trap efficiency decreases with age of reservoir as the
capacity of the reservoir reduces due to sediment deposition. Hence, the rate of
silting is higher in the initial stages and it decreases as silting takes place.
Complete filling of the reservoir may take a long time, but actually the useful life
of the reservoir is terminated when its capacity is reduced to 20 % of the
designed capacity or sediment is so much as to prevent the from serving
intended purpose.
For small reservoirs on large rivers having large inflow rates, the trap efficiency
trap efficiency is low, because the capacity inflow ratio is very small. The silting in
such reservoirs is low, because most of the sediments are passed to the
downstream.
93
The procedure for determination of the life of a reservoir is as follows:
i) The capacity of the reservoir is determined. Next, the capacity inflow ratio and
the trap efficiency are determined for the full capacity of the reservoir with the
help of the curve.
ii) Divide the total capacity into 10 parts. Assuming that 10% capacity has been
reduced due to deposits of sediments, determine the trap efficiency for the
reduced capacity (90% of the total capacity) and the inflow ratio.
iii) For the above 10% interval of capacity, determine the average trap efficiency
by taking into account the trap efficiency (η) obtained in the above steps.
iv) By collecting water samples and drying the sediments, determine the
sediment inflow rate. The total annual sediment collected in the reservoir per
year is multiplied with the trap efficiency determined in step (iii) above. The
sediment quantity so obtained is converted into hectare-meter units deposit per
year.
v) Now the reservoir volume interval ( i.e. 10% of the capacity is divided by the
sediment deposited per year, obtained in step (iv) above. It will give the number
of years required to fill this volume interval of 10% capacity.
vi) Now repeat the procedure (i) to (v) above for further capacity intervals i.e.
80%, 70%, 60%etc of the total reservoir capacity. The total life of the reservoir
will be the total number of years required to fill each of the volume intervals.
Examples
The following data is available from a reservoir regarding its trap efficiency and
capacity inflow ratio;
Cap 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Inflow
Ratio
Trap 87 92 93 95.5 96 96 96.5 97 97.5 98
efficiency
Find the probable life of reservoir with an initial reservoir capacity of 40 million
cubic meters and with an annual average flood inflow of 80 million cubic meters.
The annual sediment inflow is found to be 22 x 104 tonnes. The specific weight of
sediment may be taken as 1.15 tones per m3. The useful life of reservoir may be
assumed when its capacity is reduced to 20% of its designed capacity.
Solution.
Let’s assume that 20%of the initial capacity is filled in the first interval, i.e. 20/100
x 40 x 106 = 8 x 106m3 capacity is filled at a time.
94
For 0.5 C-I-R the trap efficiency is 96%.
At the end of the first interval, the capacity of reservoir = (40 – 8) x 106 = 32 x
106m3
At the end of the first interval the C-I-R = 32/80 = 0.4
Similarly, the silting period for the other intervals can be calculated. The rest of
the calculations are done in the table below
Example 2.
Below are data regarding trap efficiency and capacity inflow ratio of a reservoir
Capacity 0. 0. 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Inflow ratio 8 1
Trap 85 86 92 94 95 9 96 9 97 9 9
efficiency 5.5 6.5 7.5 7.7
Determine the probable life of the reservoir with an initial reservoir capacity of 25
million m3, if the annual inflow is 50 million m3 and the average annual sediment
inflow is 300,000 tonnes. Assume the density of the sediment as 1250kg/m3. The
95
useful life of the reservoir will terminate when 85 % of the initial capacity is filled
with sediment.
Solution:
Average annual sediment inflow = 300,000 x 1000 kg
Volume of sediment inflow = 300,000,000/1250 = 0.24 x 106
Initial reservoir capacity = 25 x 106m3
Annual inflow = 55 x 106m3
Therefore Initial capacity = 25 x 106 = 0.5
Capacity 50 x 106
The capacity inflow ratio goes on decreasing as the sediment takes place and
trap efficiency also decreases. The volume of interval chosen is 5 x 106 i.e. 20%
of initial capacity. To obtain the 85% filling of initial capacity the volume interval in
last two readings is taken as 5%.
Similarly, for the calculation of the years to fill the last 5% capacity is determined
by taking 1 = . 1 . = 4.87 years
St 0.2.52
Example 3. A reservoir has a drainage area of 49,000 acres. The reservoir has a
capacity of 10,000 acre-ft. Streamflow runoff averages 15.5 in. per year. The
annual sediment production is 9.5 x 106 ft3. Determine the life of the reservoir,
assuming that the life is over when 80% of the original capacity is lost
Solution
1. Annual inflow I = (49,000) (15.5/12) = 63,290 acre-ft
96
2. Annual sediment production = 9.5 x 106 ft3 or 218.5 acre-ft
3. Unfilled capacity (20%) = (20/100) (10,000) = 2,000 acre-ft
4.
ii) The ratio of reservoir capacity and drainage size: If the storage capacity of the
reservoir is much less than the annual inflow into reservoir, large quantity of
water will flow out of reservoir and with it large volume of sediment.
iii) Design of Reservoir: For the release of density currents, adequate outlets
such as sluices should be provided in the design of the dam. I this way, sufficient
quantity of sediment may be flashed fro the reservoir.
97
Post-Construction Methods:
i) Mechanical Stirring of Sediment: The sediments which are deposited in the
reservoir bed are scrapped and stirred by mechanical equipments and kept in
moving state. The sediments mixed water is pushed towards the sluices and is
flashed out.
ii) Excavation: It is a very expensive method. Iy is expensive to excavate and
also expensive to dispose of excavated material.
iii) Dredging: It is also expensive as it require the purchase and maintenance of
costly machines as dredger.
OPERATION OF RESERVOIRS
After a dam is built, it must be operated correctly. The key person is the operator,
who makes decisions about when to release or store water. The operator may be
a part-time worker who lives near a lake and who occasionally operates a gate,
or a highly trained engineer working at a remote location who makes system
decisions based on computer forecast. In the past, reservoir operating decisions
were made by RULES CURVES, which provided the operator with simple
guidelines about how much water to release and what lake levels to maintain. As
the science of forecasting and the use of computers has become more complex,
however, reservoir operation has become more sophisticated. It is not
uncommon to have a reservoir control center where operators use computers to
monitor weather forecasts furnished from satellite data and simulate future
demands for water to make decisions about water releases. They may also be
bound by legal requirements to release water for downstream users, including
fish and wildlife.
98
Qi – Q0 = DS
99
Snohomish County PUD No. 1 (SnoPUD) operates Spada Reservoir using two rule
curves based on storage (Figure 8.3). In addition,”poweroff" storage may be used during
extreme conditions (Spada reservoir has yet to reach the poweroff storage). Comparing
the reservoir storage to the Full Storage condition (FS), Upper normal rule curve, and
Lower normal rule curve determine the state of the reservoir. If reservoir storage is
between the rule curves (State 3), then a variable amount of flow (based on a ratio of
storage to lower rule curve and time of year) is conveyed through the Powerhouse
pipeline. If storage in Spada is below the Lower normal rule curve (State 4) then
minimum flow to meet instream flow requirements and M&I demand is conveyed. The
maximum flow of 1300 cfs is conveyed if the reservoir is above the Upper normal rule
curve (States 1&2). No flow is conveyed if storage is below the "Power off" level.
FLOOD ROUTING
Flood routing is the process whereby the shape of a flood hydrograph is modified
as it passes through a hydraulic system (such as a river channel, reservoir, lake,
etc). The hydrograph resulting after passing through the hydraulic system is one
with considerably reduced peak and enlarged time base. The process can be
represented as:
100
Flood routing in fact is a technique of determining the extent of change in the
inflow hydrograph as it passes through a hydraulic system and it is done for the
purposes of:
In flood routing, all the components of the routing process (inflow, storage and
outflow) are connected by the storage (continuity) equation in the form:
[ Inflow volume in time Δt] + [Change in volume of water stored by the hydraulic
system during time Δt] = [Outflow volume in time Δt]
where I and O are the average inflow and outflow rates within the time Δt.
The hydraulic system can be a reservoir (lake) or a section of a river channel,
and so the routing process is, accordingly classified into two namely: reservoir
routing and channel (streamflow) routing.
In the above equation (1) I- the inflows are known quantities but the outflows, O
and the storage, S are unknown parameters. To solve the equation, either both O
and S have to be related to a common unknown parameter or S has to be
defined in terms of d to a common unknown parameter or S has to be defined in
terms of O. The former approach is applicable to reservoir routing whiles the
latter to channel routing.
Equation (1) when expressed in the differential form can be integrated, however,
the terms in the equation have a form that is not amenable to direct solution.
Therefore, the numerical solution expressed as numerical approximation can be
written as :
I 1 + I 2 S1 − S 2 O1 + O2
+ = ------------------ (2)
2 ∆t 2
101
where subscript 1 and 2 denote beginning and end of routing period Δt
( I 1 + I 2 ) + 2S1 + O1 − 2O1 = 2S 2
+ O2 -------------- (3)
∆t ∆t
Since the outflow and the storage are both functions of water surface elevation,
the storage equation becomes a relation between the known inflow and the
unknown elevation, from which elevation can be computed as a function of time.
Procedure
1. What is available at the beginning.
i) An inflow hydrograph I =f(t)
ii) An elevation –capacity curve H =f(S) or elevation – area curve of the
reservoir
H =f(A)
iii) An elevation – outflow curve H = f(O)
4. Select the time interval Δt for the routing. This is usually taken as the interval
of time given for the inflow hydrograph. Otherwise, the time interval is so chosen
so that one does not miss the peak values.
5. At the beginning of the routing, I1 = I2 =O1. With O1 known, read off H from H
=f(O) curve. With H now known, read off 2S/ Δt +O from the H = f(2S/ Δt + O)
curve. So at time t = 0, the following quantities are known; H, O2 and (2S2 + O2)
and they will become the quantities at the end of a fictitious previous period.
102
Time Inflow (I1 + I2) (2S/ Δt + -2O1 (2S2/ Δt+O2) H O
O)
0 I1 - - - X* H* O*
6 I2 I1+I2 X* -2O* X H O
16 I3
18 I4
26 I5
30 I6
6. The values in row corresponding to time t = 0, being the values at the end of a
routing period is transferred to become the initial values on the left hand-side of
the routing equation (3) for the next time step of the routing period.
7. The left hand side of the routing equation now has known values that yield a
value for (2S2 + O2). With this known, H is read from the plot and O is also
determined.
8. The procedure is repeated until all the inflow values are used up.
Storage data
Elevation(ft) 50 50.5 51.0 51.5 52.0 52.5 53.0 53.5
Storage(acre-ft) 231 247 277 313 353 400 452 509
Solution:
1. The discharge data, which is the outflow, is computed from the given equation
Q = 3 x (35) x H2/3 and the storage data are listed in tabular form.
Water surface Head on crest Storage Outflow (cfs) (2S/ Δt + O)
Elevation (ft) (ft) (acre-ft) Q=3 x (35) x Δt = 0.5 day
3/2
H
50 0 231 0 465.6
50.5 0.5 247 37.1 536.4
51.0 1.0 277 105.0 664.4
51.5 1.5 313 192.9 825.3
52.0 2.0 353 297.0 1009.4
52.5 2.5 499 415.0 1224.3
53.0 3.0 452 545.6 1459.5
53.5 3.5 509 687.5 1715.7
103
The routing is done in the table below
1 2 3 4 5 6 7 8
Time Inflow I 1 + I2 (2S1+O1) -2O1 (2S2+O2) H O
0 0 - - - 465.6 50 0
0.5 70 70 465.6 0 535.6 50.45 35
1.0 185 255.0 535.8 -70 720.6
1.5 360 545.0 720.6
2.0 480
2.5 300
3.0 165
3.5 80
4.0 0
4.5 0
5. After the computation, a plot of the inflow and the outflow hydrographs are
made from which the extent of peak reduction can be seen and when the peak
occurs.
Elevation (m) 54 56 58 60 62 64 66 68 70 72
Surface area 4.12 4.19 4.28 4.36 4.45 4.58 4.71 4.88 5.05 5.22
(km2)
104
Chapter 6
The construction work done at the river or canal for diversion of water to the off-
taking canal is known as head-work. Depending on the purpose, head-works
may be:
1. storage head-work or
2. diversion head-work
Storage head-work is constructed to store water for some purpose
The main purpose of diversion head-work is to divert the required quantity of
water into the off-taking canals for irrigation purpose. A diversion head-works
may also serve the following purposes:
105
i) to raise the water level in the river for increasing its command
ii) to regulate the intake of water into the canal
iii) to control silt entry into the canal
iv) to store water for meeting emergency requirements
v) to prevent fluctuations in the level of supply of the river.
i). Diversion head-work provides an obstruction across a river, so that the water
level is raised and water is diverted to the channel at required level. The increase
water level helps the flow of water by gravity and results in increasing the
commanded area and reducing the water fluctuations in the river.
ii). Diversion head-work may serve as silt regulator into the channel. Due to the
obstruction, the velocity of the river decreases and silt settles at the bed. Clear
water with permissible percentage of silt is allowed to flow through the regulator
into the channel.
iii).To prevent the direct transfer of flood water into the channel.
In rivers with a stable base flow and a high enough water level throughout the
year in relation to the bed level of the intake canal, one can resort to run-off-river
water supply (Figure 38 and Example 1 in Figure 37).
106
A simple offtake structure to control the water diversion is sufficient. The offtake
should preferably be built in a straight reach of the river (Figure 39). When the
water is free from silt, the centre line of the offtake canal could be at an angle to
the centre line of the parent canal. When there is a lot of silt in the system, the
offtake should have a scour sluice to discharge sediments or should be put at a
90° angle from the parent canal. If it is not possible to build the offtake in a
straight reach of the river, one should select a place on the outside of a bend, as
silt tends to settle on the inside of bends. However, erosion usually takes place
on the outside of the bend and therefore protection of the bank with, for example,
concrete or gabions might be needed.
107
The functions of the offtake structures are:
To pass the design discharge into the canal or pipeline
To prevent excessive water from entering during flood
108
Considering these functions, the most important aspect of the structure is the
control arrangement, which can be a gate, stop logs, or other structures. When
the gate is fully opened, the intake behaves like a submerged weir (Figure41)
and its discharge is given by equation 23
109
A weir should be located in a stable part of the river where the river is unlikely to
change its course. The weir has to be built high enough to fulfill command
requirements. During high floods, the river could overtop its embankments and
change its course. Therefore, a location with firm, well defined
banks should be selected for the construction of the weir. Where possible, the
site should have good bed conditions, such as rock outcrops. Alternatively, the
weir should be kept as low as possible. Since weirs are the most common
diversion structures, their design aspects will be discussed below.
Divide wall
Fish ladder
Head Regulator
Weir
CANAL
Under sluices
Marginal Bunds
d/s GUIDE BUND
1. Weir (Barrage) may be masonry; rockfill or cement concrete. All weirs normally
have the following components: i) body wall usually with shutters; ii) upstream
rough stone or boulder pitching; iii) upstream curtain wall; iv) upstream
impervious apron; v) crest shutters; vi) downstream impervious apron; vii)
downstream curtain wall; vii) downstream apron for channel bed.
110
Shutters
Rockfill weir
111
General Equation of Weirs.
The flow over any type of weir is given by
112
Where: Q = Discharge (m3/sec) C1 = Coefficient related to condition of
submergence and crest shape (Figure 43) C2 = Coefficient related to crest
shape (Figure 44)
B = Weir length, i.e. the weir dimension across the river or stream (m) H = Head
of water over the weir crest (m)
113
114
Design of weir
a) Design of Weir Walls. Weir walls are design as a retaining wall or masonry
dam of solid gravity type. The two main loads to be considered are the horizontal
water thrust, P, and the self – weight of the weir body wall, W. The resultant of
forces P and W should pass through the middle third of the base of the weir, for
the weir to be safe.
Bligh has recommended the following:
115
H +d
b= −−−−−−−−−−−−−
S (D.H -1)
and ..a = 0.552( H + d ).. − − − − − − − − −
From the above recommendations, stability are checked under the following
conditions:
I) The water in upstream is at crest level or the crest of shutter and there
is no flow over the weir
II) The weir is submerged and water is passing over it.
III) Water is passing over weir crest and the weir is discharging with a
clear overfall.
Shutter
Body of weir
u/s apron
d/s apron
Lu / s = 2.208 CH -----------------------( )
where C = creep coefficient
H = head of water stored against the body wall
The thickness of this apron is about 30 cm placed over 30 to 50 cm thick
concrete foundation.
116
4 H −h
Thickness of apron: tD / s = -------------------(D.H - 3)
3 SG − 1
Where tD/s = is in metres; H = total head causing seepage.
h = head lost by the creeping water u to the point where thickness is to be
calculated.
SG = specific gravity of material of the apron.
The floor of the weir may burst due to inadequate weight of the weir and the
heavy uplift pressure. This bursting of the floor reduces the effective length of
the impervious floor, which will resulting increasing exit gradient, and can cause
failure of the weir.
Rapture of floor could be prevented by:
i) providing adequate length of impervious floor
ii) providing sufficiently thick impervious apron
iii) by providing a pipe at the upstream end to reduce the uplift pressure at the
downstream end.
117
If H is the total loss of head, then the loss of head per unit length of the creep ©
shall be
H H
C' = =
L [ B + 2( d1 + d 2 + d 3 ) ]
Bligh called the loss of head per unit length of creep as Percolation coefficient.
The reciprocal, (L/H) of the percolation coefficient is known as the coefficient of
creep C.
d1 d2
B=L d3
Design Criteria.
For design, Bligh proposes the following two important criteria:
118
1. Safety Against Piping: For safety against, the creep length L should be
sufficient to provide a safe hydraulic gradient, depending on the type of soil.
Where L =safe creep length; H = H1 – H2 Total loss of head for the length of
creep;
H1 = upstream head; H2 = downstream head; C = creep coefficient
2. Safety against Uplift: The uplift pressure = γwH! ; where H! = the uplift pressure
head at any point of the apron.
The downward resisting force per unit length of apron FR = t. γw . S
where t = thickness of the apron floor; S = specific gravity of floor material.
Equating the two forces, we determine the limiting floor thickness
γw H! = t. γw .S
H! = t . S
Or H! – t = t . S – t = t(S – 1)
H! −t h
Or t= ≅ - the limiting thickness
S −1 S −1
where h is the ordinate of the hydraulic gradient line measured above the top of
the floor.
Now applying a safety factor of 4/3
4 h
t=
3 S −1
Divide Wall
This is a long solid wall of stone masonry built at right angles to the axis of the
body of the weir. Its function is to divide the river channel into two portions and to
create a still pond near the canal head regulator end. It is usually built in stone
masonry but it could be done in earth with stone pitching or in concrete. Because
of the divide wall a still pond is created. This helps the deposition of silt before
the water is allowed to pass through the head regulator. It also protects the
regulator from the direct turbulence of river flow.
119
Concrete slab
1.5
Stone Masonry
Stone pitching
Concrete foundation
Divide Wall
Many rivers carry substantial sediment loads, especially during the rainy season,
in the form of sand, silt, weeds, moss and tree leaves. Approximately 70% of all
suspended and bed load sediments travel in the lower 25% of the flow profile.
While suspended silt can be beneficial to the scheme by adding nutrients to the
farmland, coarse sediments usually cause problems once they are blocked by a
weir or other diversion structure. Headworks have to be adapted to these
sediment loads to avoid silting of canals and structures. A properly-designed
intake should divert only the relatively clean upper part of the water flow into the
canal and dispose of the lower part down the river. A sluice should
therefore be incorporated into the diversion structure design. It should be placed
in line with the weir near the canal intake (Figure 51). Its seal level is generally
placed at the river bed level while the floor to the intake gate should be located
higher (Figure 52).
The control arrangement in the scour sluice generally consists of a series of stop
logs (timber, concrete) or a sluice gate. This arrangement allows the water to be
raised when there are very few or no sediments in the water. During the flood
season, the sluice is permanently open or opened at
regular intervals so that depositions of sediments can be flushed away. The guide
wall prevents lateral movement of sediments deposited in front of the weir and
separates the flow through the sluice and the flow over the weir.
Under Sluices (Scour Sluices): These are openings provided in the body wall of
the weir. Their main function is to prevent the obstructions to the flow of water
through the main sluice. They also transport deposited silt in front of the head
regulator at the upstream side, to the downstream side, thus preventing the entry
of bed silt in the canal. They create a clear, unobstructed river channel at the
approach portion of the head regulator. They also help to reduce maximum flood
level. They are provided near the wing walls of the weir. They are operated by
120
means of gates provided for this purpose. The gates are operated by levers
provided at the top of the weir.
In designing the scour sluices, the following are to be noted:
i) the design flow should be at least twice the discharge in the off-take
channel
ii) the crest of the under sluice should be one metre lower than the crest
of the head regulator.
iii) The downstream portion of the under sluice must be well protected to
prevent downstream erosion.
The head regulator consists of a raised crest with abutments on both sides. Piers
may be constructed on the crest to create a number of spans and to support a
roadway and a platform at its top to operate the gates.
Marginal Bunds: are earthen embankments constructed parallel to the river along
the banks. They are located near the approaches of the diversion head works.
Their main function is to control the floods and prevent submergence of the areas
behind the hydraulic structures. Their design and construction are similar to earth
dams.
Guide Banks (or Bell’s bunds): are constructed almost parallel to the flow of the
river, near the hydraulic structure. They are bunds with a straight length and a
curved end at their flanks. Their function is to guide the flow into and protect the
hydraulic structure by controlling and regulating the flow of the river.
121
Marginal Bunds(embankment)
Stone pitching
Marginal Embankment
122
Chapter 7
DISCHARGE MEASUREMENT
Almost any kind of obstacle that partially restricts the flow of water in an irrigation
canal and provides a free fall, to ensure that upstream and downstream flow are
independent, can be used as a measuring device, provided that it can be
calibrated. Standard structures, which have already been accurately described
and calibrated, exist. Weirs, flumes and orifices are the devices that are normally
used for discharge measurement.
Conservation of mass
Q = A x V = Constant
Conservation of energy
Equation 38
123
Equation 38 sums up the pressure head, velocity head and gravitational head to
give the total head. For an open canal, the pressure head equals the water depth
h. When there is frictional loss along the flow path, an expression for frictional
head loss must be included. Thus applying the Bernoulli Equation to two
successive cross sections along a flow path results in:
Equation 39
The numbers 1 and 2 refer to the first and second cross section. HL is the
frictional head loss.
Specific energy
The concept of specific energy is used in the analysis of critical flow. At any
cross-section of a canal, the energy with respect to the canal bed is referred to
as specific energy. It is derived from the Bernoulli Equation according to the
following equation:
124
Where: E = Specific energy (m), h = Depth of flow (m) g = Gravitational
force (9.81 m/sec2)
The cross-sectional area varies with the depth of flow only if the geometry of the
canal is constant. Therefore, for a given discharge the specific energy is a
function of depth alone
Rectangular canal
125
Plotting E against h for different values of (Q/b) gives curves as shown in
diagram above.
The curves show that, for a given discharge and specific energy, there are two
alternate depths of flow, which coincide at a point where the specific energy is a
minimum for a given discharge. Below this point, flow is physically not possible.
At this point flow is critical and it occurs at critical depth and velocity. At a greater
depth, the velocity is low and flow is sub-critical. At the lesser depth, the velocity
is high and flow is super-critical. For sub-critical flow, the mean velocity is less
than the
velocity of propagation of stream disturbances such as waves. Thus, stream
effects can be propagated both upstream and downstream. This means that
downstream conditions affect the behaviour of flow. When flow is supercritical,
the velocity of flow exceeds the velocity of propagation. Consequently, stream
effects (for example, waves) cannot be transmitted upstream, and downstream
conditions do not affect the behaviour of the flow. For critical flow, the specific
energy is a minimum for a given discharge. In this case, a relationship exists
between the minimum specific energy and the critical depth. This relationship is
found by differentiating Equation 41 with respect to h, while Q remains constant.
This gives:
Froude Number
Where: Fr = 1 for critical flow, Fr = > 1 for super-critical flow, Fr = < 1 for sub-
critical flow
If a structure is built in a canal which has sub-critical flow, it may cause the flow
to pass through the critical to the super-critical state. This means that the state
126
upstream of the structure becomes independent of the state downstream. This
can either be achieved if the structure narrows the canal, which means
increasing the (Q/b)- ratio without altering the specific energy, or if it raises the
canal bed, which means reducing the specific energy without altering the
discharge per unit width. That is how critical flow is obtained with a measuring
device. A control section in a canal is a section that produces a definitive
relationship between water depth and discharge.
Hydraulic jump
Equation 43
6.6.2. Weirs
127
A weir in general is any regular obstruction over which flow occurs. Weirs, especially the
sharp-crested ones are commonly used for measuring large and small open flows in the
field or laboratory.
Head race
H Ho
vo
z
δ
P
Tail race
LB
B
The zone of fluid before the wall of the weir is called the head race while the zone of
water after the weir is called the tail race
The distance LB from the wall of the weir to the point where a fall in the level of stream is
noticeable, known as the drawdown section is given by:
LB = (3- 5)H
H - measured at the section B-B is called the geometric head which is the rise of fluid
above the crest of the weir
z – geometric drop in the level of fluid at the weir. It is the difference between the head
race and the tail race.
vo – approach velocity
Ho – The total head which is given by the sum of the geometric head and the velocity
head.
zo – total drop in head over the weir which is the sum of the geometric drop plus the
velocity head.
128
Classification of weirs
a) Rectangular
b) Triangular
c) Trapezoidal
d) Circular
e) Parabolic
b) Broad – crested weir 2H ≤ δ < 8H: In broad-crested weir, the falling nappe is
supported on the crest and does not allow the nappe to fall free.
a) Suppressed weir: - When the length of the weir is the same as the width of the
channel (stream) It is so called because in this case, the effect of sides or ends of the weir
on contraction of the nappe is eliminated or suppressed. The width of the nappe is the
same as the width of the stream.
b) Weirs with end contractions: - In this case, the length of the weir is less than
the width of the channel. Therefore there is a lateral contraction of the nappe so that the
length of the nappe is less than the length of the weir crest.
4. By the level of downstream water level relative to the crest of the weir.
a) Un-submerged weir: - When the water level on the downstream side of the weir
is below the crest of the weir, it is said to be un-submerged.
b) Submerged weir: When the water level on the downstream side of the weir
rises above the level of the crest, the weir is said to be submerged.
The flow of fluid over weirs is extremely complex phenomenon and virtually impossible
of exact analytical solution. As result, all weir formula are derived by making lots of
129
simplification of the problem which leads to an approximate solution which is somehow
rectified by the use of experimental coefficients derived from experimental set ups.
To derive a simplified relationship for a weir, the following assumptions are made:
ii) all fluid particles move horizontally as they trave pass the weir crest
iv) the influence of viscosity, surface tension, turbulence are all neglected
v12/2g 2
1 1 h
H H
dh
2
Weir crest
v1
P
B
Fig. Simplified flow over weir with the above assumptions
u = cv 2gh..
dQ = Bdhcv 2 gh
The theoretical flow passing over the weir is obtained by integrating the above equation
over the geometric head.
130
H
2
Q = ∫ Bcv 2 gh dh = Bcv 2 g H 3 / 2
0
3
To obtained the actual flow rate over the weir due to the above assumptions, we introduce
an experimental coefficient of discharge
2
Q= CD B 2 g H 3 / 2
3
If the approach velocity vo is appreciable (>1.0), then instead of geometric head one has
to use the total head Ho
As can be seen, the formula for flow over weir seems simple and easy to use however, the
difficulty is with the correct determination of the discharge coefficient. In the laboratory,
this could be determined quite easily and accurately, however, it becomes extremely
difficult to determine this accurately in the field. A few empirical formulae exist for the
determination of the discharge coefficient for a rectangular weir.
Contracted Weirs
For weir with end contractions, Francis observed that the effect of the contraction is to
reduce the length of nappe by 0.1H. Thus if there are n end contractions, then the length
of the nappe will be reduced by 0.1nH and will be equal to (B – 0.1nH) and the discharge
equation will then be given as:
2
Q= CD ( B − 0.1nH ) 2 g H 3 / 2
3
A triangular weir (notch) is particularly useful where the discharge is to vary over a large
range and the same accuracy is desired for both small and large discharges because of the
following advantages:
131
i) The coefficient of discharge for a triangular weir is practically independent of the head.
This is because, for all heads the ratio of the head to the wetted length of crest is constant.
But in a rectangular weir, the ratio of head to the wetted length of crest is not constant.
ii) When the discharge is small, a triangular notch provides a greater head than the
rectangular notch. Hence head measurement can be done more accurately over the
triangular notch than over the rectangular notch.
iii) When the flow rate is small, there is the tendency of a clinging nappe in a rectangular
notch. However, for the same flow rate over a triangular notch, the head will be greater
and clinging nappe will be avoided.
b h
dh
H
θ
The equation for the discharge over the vee-notch can be obtained just like we did in the
rectangular weir by considering an elementary strip of width b and height dh.
Considering the notch as a large orifice, the velocity of flow is given by:
u = cv 2gh..
dQ = bdhcv 2 gh
θ
Putting b into the discharge equation, dQ = 2( H − h) tan
2
cv 2 gh dh
132
H
θ θ
Q = ∫ 2 tan
2
( )
cv 2 g Hh1 / 2 − h3 / 2 dh =
8
15
cv 2 g tan H 5 / 2
2
0
The above equation represents the theoretical discharge without consideration of the
contraction of the flow area. The actual discharge is obtained by applying the coefficient
of discharge CD. to obtain:
θ
5
8
Q= CD 2 g tan .H 2
15 2
Q = 1.42 H5/2
Q = 1.4 H5/2
2. By King
Q = 1.343 H2.47
The coefficient of discharge is made up mainly of the contraction coefficient and the
velocity coefficient. The contraction coefficient of a notch depends on the wetted
perimeter. For a rectangular notch, the wetted perimeter does not vary with head because
the base, which forms part of the wetted perimeter remains constant for all heads.
Consequently, the coefficient of contraction and hence the coefficient of discharge will
not be constant for all heads. In a triangular notch, there is no base and therefore the
contraction is due to the sides only. The wetted perimeter depends on the length of the
sides which in turn depend on the head. Consequently, the coefficient of discharge is
fairly constant in a triangular notch.
Trapezoidal Weir
133
H
θ/2
θ/2
Trapezoidal Weir
The trapezoidal notch is equivalent to two notches; one rectangular with length b and the
other a triangular with an apex angle of θ. The total discharge over the trapezoidal notch
would be equal to the sum of discharges over the equivalent rectangular and triangular
notches.; thus
Qtrap = Q1 + Q2 =
θ θ
5
2 8 2 8
C Db 2 g H 3 / 2 + C D 2 g tan H 2 = C D* 2 g H 3 / 2 b + H tan
3 15 2 3 15 2
Cippoletti Weir
The principle behind the Cippoletti weir is based on the fact that in a rectangular weir,
due to end contractions, the effective length of the weir is reduced and consequently the
discharge. On the other hand, by the addition of side slopes to a rectangular weir, gives an
increase in discharge. Therefore if the side slopes of a trapezoidal weir are so adjusted
such that the reduction in discharge due to end contractions is just equal to the increase in
discharge due to the addition of side slopes, the net effect I discharge is zero and
therefore the standard weir formula without end contraction can be used.
The Cippoletti Weir is a trapezoidal weir with side slopes 1 horizontal to 4 vertical such
that the discharge in effect equal to that of a rectangular weir without end contraction.
3
2
For Cippoletti weir, Q= C Db 2 g H 2
3
Cipolletti gave the equation for the discharge through such a weir as:
Q = 1.86 b H3/2
Submerged Weir
134
A submerged weir (notch) is one in which the level of fluid in the downstream side of the
weir is above the crest of the weir.
1 H
A
H1 B
2 H2
The flow over such a weir can be obtained by dividing the flow into two portions:
i) Flow over the upper part of the section (above line AB) of depth (H1 – H2)
may be considered as a free from a weir into the atmosphere and
ii) Flow through the lower part of depth (below line AB) may be considered as a
discharge through a submerged orifice.
The discharge over the upper portion behaving as a free weir is given by
2 2
C D b 2 g ( H1 − H 2 )
3/ 2
Q1 = = CD b 2 g H 3 / 2
3 3
Q2 =C D bH 2 2 g ( H1 −H 2 ) =C DbH 2 2 gH
2
Q= C Db 2 g H 3 / 2 + C DbH 2 2 gH
3
135
0.385
H n
Q = Q1 1 − 2
H1
Broad-Crested Weir
Broad-crested weirs (long-based weirs) are weirs that have crest lengths that are
sufficiently longer (δ > 0.67H) to prevent the nappe from springing free. These weirs are
usually constructed of concrete, have rounded edges, and are capable of handling much
larger discharges than sharp-crested weirs.
Some advantages of broad-crested weirs.
i) It is cheaply and easily fabricated (can be done in-situ or off-site)
ii) It is easy to install
iii) It possesses a wide modular range;
iv) It produces minimum afflux(increases in upstream depth due to weir
installation)
v) It requires a minimum of maintenance.
Flow over a broad-crested weir must pass through the critical depth somewhere near the
downstream corner. Broad-crested weirs are characterized by the presence of two drops
in the water level over the weir.
1 2
ZB(o) ZB
Ho H
Zc
v h
136
By taking Bernoulli’s equation in sections 1 and 2 , the velocity of flow on the weir is
obtained as
v =cv 2 gZ B0 =cv 2 g ( H o −h )
dq [
d c v h 2 g ( H o − h) ]
= =0
dh dh
1
(
d h Ho − h )= Ho − h − 2
h
=0
dh Ho − h
2 2 4
Q = b. H o cv 2 g H o − H o = bcv 2 g H o3 / 2 = 0.385cv b 2 g H o3 / 2
3 3 27
Q = CDb 2 g H o3 / 2
q = CD 2 g H o3 / 2
137
Ventilation of Weirs
When flow takes place over a suppressed weir with wing walls on the downstream side,
the free access of air below the nappe is prevented. The air entrapped between the nappe
and the weir face is gradually evacuated by the dynamic action of the flowing fluid. This
eventually leads to subatmospheric pressure in the underside of the nappe bringing about
increases in the discharge.
Investigations by Bazin have shown that if the flow is sufficient to prevent the air having
access to the underside of the nappe, it may assume one of the forms shown below
Sub-atmospheric
pressure
Atmospheric
Pressure
Fully aerated (free, ventilated) nappe: In this type, the weir discharges free as the
pressure below the nappe is maintained at atmospheric by proper ventilation.
Drowned (depressed) nappe: When the underside of the nappe is poorly aerated, the air
enclosed between the nappe and the downstream face of the weir gradually mixes up with
the flowing fluid and is carried away. This causes the pressure under the napper to fall
below atmospheric. The pressure difference on the two sides of the nappe-(the upper side
exposed to atmospheric while the bottom side has sub-atmospheric pressure) –causes the
nappe to deflect towards the weir.
Clinging nappe: This type occurs when there is no ventilation and all the air under the
napper has been carried away by the dynamic action of the flowing fluid. The nappe is
forced to adhere to the downstream face of the weir by the atmospheric pressure on the
upper side. In this type, the discharge is 20 to 30 per cent more than in free nappe.
Knowledge of these types of nappes is necessary especially when the weir is used for floe
(discharge) measurement. All the discharge formulae have been derived for free nappe
138
and therefore some adjustment in the measured flow must be made when other forms of
nappe occurs during discharge measurement. For discharge measurement it is always
advisable to ensure that the nappe is fully aerated.
FLUMES
Although weirs are the cheapest and simplest structures for flow measurement in open
channels, however the relatively high losses caused by weirs and the tendency of
sedimentation of suspended particles due to reduced velocities in the vicinity of the weir,
could in certain cases pose important limitations Flumes provide a convenient alternative
to weirs for flow measurement in open channels where high head losses and
sedimentation are of concern. Such cases include flow measurement in wastewater
treatment plants and irrigation channels with suspended particles
Flumes are devices in which the flow is locally accelerated by a streamlined lateral
contraction in the channel sides. A flume has: a convergent section, in which the
flow accelerates, a throat section and a divergent section, in which the
flow returns to normal
Flumes are of two types: non-modular or the venturi flume and the modular or the
standing wave flume.
In the non-modular flume, the velocity at the throat is maintained below the critical value
so that no standing wave is produced. However, in the modular type, the flume is
designed such that the velocity of flow at the throat is greater than the critical velocity
thereby resulting in a standing wave within the flume.
1
2
b1 b2
PLAN
1 2
v12/2g
v22/2g
h1
h2
Critical depth
ELEVATION
139
When the width of a channel is reduced while the bed remains flat, the discharge per unit
width increases. If the losses are neglected, the specific energy remains constant. If the
conditions are made such that the free surface does not pass through the critical depth, the
arrangement forms a Venturi flume.
Referring to the diagram above, for continuity of flow,
h1 + v12/2g = h2+v22/2g
Substituting v1 for continuity equation,
h1-h2 = v22/2g[1-(h2b2/h1b2)2]
2 g ( h1 − h2 )
v2 =
b h 2
1 − 2 2
h1b
2 g ( h1 − h2 )
Q = v2 A2 = b2 h2 .. − − − −Theoretical..disch arg e
b h 2
1 − 2
b h1
2 g ( h1 − h2 ) A1 A2
Q = CDb2 h2 = CD 2 g ( h1 − h2 )
b h 2 A12 − A22
1 − 2 2
b1h1
CD = 0.95 to o.99
Therefore to measure the flowrate, water depth are measured at two locations, one at
section 1-1 and the other at section 2-2.
h1 + v12/2g = h2 + v22 = H
140
v2 = (2g(H-h2)1/2
Q =b2 (
2 g Hh22 −h 3 ) =b 2 (
2 g Hh2 −h 3 )
1/ 2
(
d Hh22 − h3 )
=0
dh
Therefore 2 Hh2 − 3h22 = 0
2
h2 = H
3
Thus for the flow to be maximum, the depth at the throat section should be two-thirds the
total energy.
3h
For this condition, the velocity v2 = 2 g ( H − h2 ) = 2 g 2 − h2 = gh2
2
This implies that for maximum discharge, the depth at the throat must be equal to the
critical depth and the velocity critical velocity.
Substituting h2 = 2/3H in the discharger expression, we shall obtain
1/ 2
2 2 2
Q = b2 2 g H H − H = 2 gb2 H 3 / 2 = 1.705b2 H 3 / 2 = 1.705b2 ( h1 + v12 / 2 g )3 / 2
3 3 3 3
Although flat bed modular flume is easier to construct, it is sometimes necessary to raise
the invert in the throat to attain critical conditions. For the flat or raised floor, the throat
length should ideally be sufficient to ensure that the curvature of water surface is small to
make the water surface parallel to the invert. Such ideal flumes become relatively long
and expensive. A more compact, short and less expensive groups of flumes have been
constructed. However, the surface water profile varies rapidly and as such the theoretical
analysis for ordinary flumes is not wholly applicable. Therefore empirical relationship are
usually developed for such flumes (e.g; Parshall Flume), which make each such flume
standard with its unique calibration curve.
141
The Parshall flume consists of a converging section with a level floor, a throat section
with a downward sloping floor and a diverging section with an upward sloping floor.
Flume sizes are known by their throat width and each size has its own characteristics,
which is unique (see Table 31).
The flow through the Parshall flume can occur either under free flow or under submerged
flow conditions. Under free flow the rate of discharge is solely dependent on the throat
width and the measured water depth, ha. The water depth is measured at a fixed point in
the converging section.
The upstream water depth-discharge relationship, according to empirical calibrations, has
the following general form:
Q = K x (ha)u
Table 32 gives the values for K and u for each flume size.
When the ratio of gauge reading hb to ha exceeds 60% for flumes up to 9 inches, 70% for
flumes between 9 inches and 8 feet and 80% for larger flume sizes, the discharge is
reduced due to submergence. The upper limit of submergence is 95%, after which the
flume ceases to be an effective measuring device because the head difference between ha
and hb becomes too small, such that a slight inaccuracy in either head reading results in a
large discharge measurement error
142
143
The discharge under submerged conditions is:
Qs = Q - Q c
Figure 77 gives the corrections Qc for submergence for flumes with 6 inch,
9 inch and 1 foot throat width. The correction for the 1 foot flume is made
applicable to other sizes by multiplying the correction Qc for the 1 foot by
the factors given in Figure 77 (1 foot flume).
144
145
146
147
148
The weir is the most practical and economical device for water measurement.
Weirs are simple to construct, easy to inspect, robust and reliable. Discharge
measurement weirs can either be sharp-crested (Figure 69, 70, 71) or
broadcrested (Figure 72).
149
Sharp-crested weirs
Sharp-crested weirs, also called thin plate weirs, consist of a smooth, vertical, flat
plate installed across the channel and perpendicular to the flow (Figure 69). The
plate obstructs flow, causing water to back up behind the weir plate and to flow
over the weir crest. The distance from the bottom of
the canal to the weir crest, p, is the crest height. The depth of flow over the weir
crest, measured at a specified distance upstream of the weir plate (about four
times the maximum h1), is called the head h1. The overflowing sheet of water is
known as the nappe.
Thin plate weirs are most accurate when the nappe springs completely free of
the upstream edge of the weir crest and air is able to pass freely around the
nappe. The crest of a sharp-crested weir can extend across the full width of
channel or it can be notched. The most commonly used notched ones are:
Rectangular contracted weir
Trapezoidal (Cipoletti) weir (Figure 70)
Sharp sided 90° V-notch weir (Figure 71)
The type and dimensions of the weir chosen are based on the expected
discharge and the limits of its fluctuation. For example, a V-notch weir gives the
most accurate results when measuring small discharges and is particularly
adapted to the measurement of fluctuating discharges.
Calibration curves and tables have been developed for standard weir types. The
conditions and settings for standard weirs are as follows:
i. The height of the crest from the bottom of the approach canal (p) should
preferably be at least twice the depth of water above the crest and should in no
150
case be less than 30 cm. This will allow the water to fall freely, leaving an
airspace under and around the jets.
ii. At a distance upstream of about four times the maximum head a staff gauge is
installed on the crest with the zero placed at the crest elevation, to measure the
head h1.
iii. For the expected discharge, the head (h1) should not be less than 6 cm and
should not exceed 60 cm.
iv. For rectangular and trapezoidal weirs, the head (h1) should not exceed 1/3 of
the weir length.
v. The weir length should be selected so that the head for the design discharge
will be near the maximum, subject to the limitations given in (ii) and (iii).
vi. The thickness of the crest for sharp-crest weirs should be between 1-2 mm.
In sediment-laden canals, a main disadvantage of using weirs is that silt is
deposited against the upstream face of the weir, altering the discharge
characteristics. Weirs also cannot be used in canals with almost no longitudinal
slopes, since the required difference in elevation between the water
levels upstream and downstream side of the weir is not available. Discharge
equations for weirs are derived by the application of the Continuity and Bernoulli
Equations (Equation 12 and 38 respectively). In each case, a discharge
coefficient is used in order to adjust the theoretical discharge found by laboratory
measurement
Where:
151
Q = Design discharge over weir (m3/sec) b = Length of weir crest (m) h =
Design water depth measured from the top of the weir crest (m) Table 28 gives
discharge data related to length of crest, b, and water head, h, over a weir.
Where:
Q = Design discharge over weir (m3/sec) b = Length of weir crest (m) h =
Design water depth measured from the top of the weir crest (m)
152
V-notch weir
A V-notch weir has two edges that are symmetrically inclined to the vertical to
form a notch in the plane perpendicular to the direction of flow. The most
commonly used V-notch weir is the one with a 90° angle. Other common V-
notches are the ones where the top width is equal to the vertical depth (1/2 x 90°
V-notch) and the one where the top width is half of the vertical depth (1/4 x 90° V-
notch) (Figure 71). The V-notch weir is an accurate discharge-measuring device,
particularly for discharges less than 30 l/sec, and it is as accurate as other types
of sharp-crested weirs for discharges from 30 to 300 l/sec (U.S. Department of
Interior, 1975). To operate properly, the weir should be installed so that the
minimum distance from the canal bank to the weir edge is at least twice the head
on the weir. In addition, the distance from the bottom of the approach canal to the
point of the weir notch should also be at least twice the head on the weir (U.S.
Department of Interior, 1975).
The general and simple discharge equation for a V-notch
153
Where:
Q = Design discharge over the weir (m3/sec) θ = Angle included between the
sides of the notch (degrees)
h = Design water depth (m)
Table 30 gives discharge data for the three common V-notches related to water
depth (head) and angle°.
154
Broad-crested weir
A broad-crested weir is a broad wall set across the canal bed. The way it
functions is to lower the specific energy and thus induce a critical flow (Figure 72
Flumes
Discharge measurement flumes are extensively used in irrigation schemes
mainly because they:
Can be used under almost any flow condition
155
Have smaller head-losses than weirs, thus are more
accurate over a large flow range
Are insensitive to the velocity of approach
Are relatively less susceptible to sediment and debris transport
However, major disadvantages of flumes include the relative large sizes and the
accurate manufacturing/construction workmanship required for optimum
performance (James, 1988).
A canal section that causes flow to pass from sub-critical through critical to the
super-critical state forms a control and the discharge is a single valued function
of the upstream water level. Critical flow can be achieved by raising the canal
bed, thereby reducing the specific energy, or by decreasing the canal width,
thereby increasing the discharge per unit width (see Section 6.6.1). This latter
technique is the one used by flumes. A flume has:
Super-critical flow passing from the throat will return to sub-critical flow
downstream of the flume. This occurs due to the development of a hydraulic
jump, which is induced within the divergent section by a sill or other barrier.
Where there is sufficient head available, the divergent section of the
flume could be avoided as the flow could fall freely in a stilling basin. In this case,
weirs could also be used. However, if canals are expected to carry a lot of
sediment, the flume should be the better choice. Flumes are most commonly
rectangular or trapezoidal in cross-section. The former type is the most simple to
construct, but if the canal cross-section is not rectangular there is a risk that
unpredictable flow patterns will result from an abrupt change of cross-section.
Trapezoidal flume
Whenever the canal section is not rectangular, trapezoidal flumes such as those
shown in Figure 79, are often preferred, especially for measuring smaller
discharges. A typical trapezoidal flume has an approach, a converging section, a
throat, a diverging and an exit section. A minimum transition will be required. An
additional advantage is the flat bottom, which allows sediment to pass through
fairly easily. Furthermore, the loss in head may be less for comparable
discharges. Trapezoidal flumes are particularly suited for installation in concrete-
lined canals. The flume should normally be put on top of the lining, thus
constricting the flow section to the extent required for free flow conditions over a
156
whole range of discharges up to the canal design discharge. As a rule of thumb,
one can say that the lower the canal gradient the higher the elevation of the
flume above the canal bed level. The flow characteristics of the flume can be
determined experimentally. This allows for the calibrations of the flume.
Cut-throat flume
The cut-throat flume has a converging inlet section, throat and diverging outlet
section. The flume has a flat bottom and vertical walls (Figure 80). It is preferable
157
to have the cut-throat flume operating under free flow conditions. This facilitates
measurements and ensures a high degree of accuracy.
Free flow conditions through the cutthroat flume are described by the following
Where:
Q = Discharge (m3/sec) C = Free flow coefficient ha = Upstream water depth
(m),
K = Flume length coefficient W = Throat width (m)
158
159
Orifices
Orifices, such as gates and short pipes, are also used as water measuring
devices (Figure 82). However, they do not offer any advantage over the use of
weirs or flumes. Furthermore, their calibrations are not as accurate nor as stable
as other types of measuring devices. For weirs the discharge is proportional to
the head above the crest raised to the power 3/2 (Equations 44, 45, 46,
48). Therefore, they are sensitive to the fluctuations in the upstream water level.
For orifices, including gates and short pipes, the discharge is proportional to the
head of water above the crest raised to the power 1/2, as shown by Equation 34
(see Section 6.1.3). Therefore, they are less sensitive to small fluctuations of the
upstream water level. Under submerged conditions both the upstream and
downstream sides of the structure need water level recordings. For free flow
conditions, the discharge is a function of the upstream water depth alone.
160
The general discharge equation for a free flow orifice is (Equation 34):
Where:
161
Q = Design discharge through orifice (m3/sec) C = Design coefficient
(approximately 0.60)
A = Cross-sectional area of the orifice (m2) g = Gravitational force (9.81 m/sec2)
h1 = Water depth upstream of orifice over reference level (m) (Figure 83)
Where:
Q = Design discharge through orifice (m3/sec), C = Discharge coefficient, which
is 0.63 for sluice gates and submerged orifices and 0.85 for short pipes, A =
Cross-sectional area of the orifice (m2) g = Gravitational force (9.81 m/sec2), h1
= Water depth upstream of orifice over reference level (m), h2 = Water depth
downstream of the structure (m)
Current meters
Current meters are used to measure the velocity in a canal, from where the
discharge can be calculated using the Continuity Equation 12 (see Section 5.1).
Most current meters have a propeller axis in the direction of the current. The
flowing water sets the propeller turning. On a meter, forming part of the
equipment, the number of revolutions per time unit can be read and, by means of
162
a calibrated graph or table, the velocity can be determined. A well known type of
current meter is the Ott instrument C31 for velocities up to 10 m (Figure 85).
Propeller meters are reliable and accurate, but rather expensive. In measuring
the velocities, the number of points per vertical and the number of verticals per
cross-section should be determined. For this purpose, the quantity of work and
the time required should be weighed against the degree of accuracy (Euro-
consult, 1989). For example, measurements can be taken at 10 cm horizontal
distance over the cross-section and at 0.2h and 0.8h depth at each 10 cm (h is
the water depth). The velocity is the average of the velocity at 0.2h and 0.8h
depth. If the water depth is less than 0.5-0.6 m, one reading can be done at 0.6h.
Then, for each vertical the flow per unit width can be calculated according to q =
v av x h (Figure 86a). These qs are distributed over the total width (Figure 86b)
and the area between the q-line and the water surface gives the total discharge.
It is also possible to establish the discharge per section and to consider the sum
of the discharges in the sections as the total discharge.
163
Discharge measurement in pipelines
Several types of devices can be used to measure the discharge in pipelines. This
section will discuss differential pressure and rotating mechanical meters, as they
are the ones commonly used.
Venturi tube
The pressure drop between the inlet and throat is created as water passes
through the throat. In the section downstream of the throat, the gradual increase
in cross-sectional area causes the velocity to decrease and the pressure to
increase. The pressure drop between the Venturi inlet and the throat is related to
the discharge, as follows
164
Where:
Q = Discharge (l/min), C = Flow coefficient, D = Diameter of upstream section
(cm), d = Diameter of contraction (cm)
P1 = Pressure in upstream section (kPa), P2 = Pressure in contraction (kPa), K
= Unit constant (K is 6.66 for Q in l/min, d and D in cm, and P1 and P2 in kPa)
There are many types of rotating mechanical flow meters used in pipelines.
These flow meters normally have a rotor that revolves at a speed roughly
proportional to the discharge and a device for recording and displaying the
discharge and total volume. The rotor may be a propeller or axial flow turbine, or
a vane-wheel with the flow impinging tangentially at one or more points.
Calibration tests are usually needed to accurately relate rotor revolutions to the
flow. The lowest discharge that can be accurately measured by a rotating
mechanical flow meter depends on the amount of bearing friction that can be
tolerated while the occurrence of cavitation often establishes the largest flow rate
that can be measured (see Module 5). Head loss through most rotating
mechanical discharge meters is moderate.
165
Chapter 8
Weir or Anicut
It is a solid wall of masonry or concrete constructed along the width of the river
from one bank to the other to raise the water level on the upstream side so that
proper supplies of water may be assured to the canals. Usually water flows over
the crest of the weir but sometimes shutters 1m or more high are provided on the
crest of the weir. During the dry period the shutters are raised to raise the water
level and during rainy season, the shutters are dropped to allow flood waters to
pass without causing damage on the upstream side.
Barrage
The function of a barrage is the same as that of a weir. In a barrage, a low height
weir is constructed first and then according to design piers are constructed at
regular intervals along the width of the river. The gaps between the piers are
closed by means of gates, which can be lifted or lowered down mechanically or
manually or both. Barrages are preferred now-a-days as they offer greater
flexibility in operation and better control. In this case, water level can be raised
many times more than weirs and at will, more discharge can be passed by lifting
the gates, but in case of weirs, the crest is fixed and only a fixed amount of
discharge can be passed.
When the foundation of a hydraulic structure is pervious, seepage will take place
as long as differential head exists across the structure.
Seepage of water below the structure may result in failure due to piping or
pressure uplift.
Referring to the diagram below, assume that there are three cut-offs down
beneath the foundation of a weir with depths y1, y2, and y3 at A B and C. Then the
total seepage length (creep length)
Lc = L1 + L2 + 2(y1 + y2 + y3)
If H is the total head difference downstream and upstream, then the hydraulic
gradient is given by:
Hydraulic gradient = H / Lc = head lost per unit length travelled by seepage
water.
Head loss across cut-off A = (H/Lc) x 2y1
166
B = (H/Lc) x 2y2
C = (H/Lc) x 2y3
Based on the above, Bligh suggested the following to prevent failure by seepage.
Weir or barrage
Hydraulic gradient
H
A h Q
B C
y1 L1 y2 L2
y3
The ordinate at any point on the bottom of the floor to the hydraulic gradient
represent the uplift at that point. This uplift must be balanced by the weight of
the floor if uplift is to be prevented. Therefore the floor must be of sufficient
thickness to provide the necessary resisting force to this uplift force.
If we let t be the required floor thickness to balance uplift force , then
167
Weight of floor = γ . SG. t
Uplift force = γ.h
where γ = weight density of water, SG specific gravity of material of floor.
In the limiting condition,
.γ. SG. t = γ.h
or t = h/SG.
H1
H2
H2
H1
Uplift pressure on foundation
If h’ is the ordinate to the hydraulic gradient line considered from the top of the
floor, then: h = h’ + t
And t = h/SG = (h’ + t)/SG
Or t = h’/(SG-1)
This is the limiting thickness of the floor to withstand uplift.
A. N. Khosla’s Theory
After studying a lot of dam failure constructed based on Bligh’s theory, Khosla
came out with the following;
Seeping water below a hydraulic structure does not follow the bottom profile of
the impervious floor as stated by Bligh but each particle traces its path along a
series of streamlines.
168
For steady flow, seepage in a homogeneous soil obeys the Laplacian equation :
d 2Q d 2Q
+ 2 =0 ---------------------- ( )
dx 2 dz
where Q = K . h - flow potential
h = residual head at any point within the soil and K coefficient of
permeability as defined by Darcy.
The above equation represents two sets of curves STREAMLINES and
EQUIPOTENTIAL (VELOCITY POTENTIAL) LINES which intercept each other
orthogonally.
Streamlines
h
A B C D
Equipotential
lines
Flow network
Streamline
A streamline is the path traced out by particles of water seeping through the
subsoil. Every particle traces out its own path which represents the streamline for
that particle.
The first streamline, immediately below the hydraulic structure follows the bottom
profile of the structure and will be almost the same as Bligh’s creep path. The
others streamlines will be curves as they proceed through the pervious soil. If an
impervious boundary intervenes, then the last streamline near the boundary
follows the boundary.
Equipotential Lines.
Let us consider that there is no water on the downstream side. If we assume the
downstream bed as the datum, then the on the upstream side is h. All the floor at
the upstream side AB and every particle of water entering the subsoil on the
upstream side will be having a head h. This head will be gradually utilize as the
particle traces out the streamline and when it emerges at the downstream end,
the entire head h will have been utilized so that the head at the exit end will be
169
zero. At any intermediate position along the streamline, the particle will be
having a head h1 which is the residual head; which will be dissipated by the
particle while traveling the remaining path on the streamline. Thus every
streamline moving from a head of h to a head of zero will pass through a residual
head of h1. If all points on all streamline with a residual head of h1 are joined, we
shall obtain a curve called an equipotential line. Thus equipotential line is a line
joining points of equal residual head.
Flow Nets
The streamlines and equipotential lines always intercept at right angles. In a flow
field, the combination of the two groups of lines form a network in the flow field
called the flow net. The space enclosed between any two adjacent streamlines
and equipotential lines is known as a field. Though it is possible to draw an
infinite streamlines and equipotential lines, for the sake of practical convenience
and easy interpretation, only a limited number is usually drawn in such a way that
every field becomes an elementary square.
By drawing the flow net, all characteristics such as flow rate, velocity etc can be
obtained.
Every particle of water while seeping through the sub-soil, at any position will
exert a force f, which will be tangential to the streamline at any point. As the
streamlines bend upward, the tangential force f will be having a vertical
component f1. Also at that point, there will be a downward force W due to the
submerged weight of the soil particle. Thus at that point there will be two forces
on the particle; one upward vertical component of f, and the other, the
submerged weight. It is evident that if the soil particle is not to be dislodged, then
the submerged weight must be greater than the upward vertical component of f.
170
The upward vertical component force at any point is proportional to the water
pressure gradient dp/dx.
Hence for stability of the soil and for the prevention of erosion and piping, the
seeping water when it emerges at the downstream side, at the exit position, the
force f1 should be less than the submerged weight W. In other words the exit
gradient at the downstream end must be safe.
If at the exit point at the downstream side, the exit gradient is such that the force
f1is just equal to the submerged weight of the soil particle, then that gradient is
called critical gradient. Safe exit gradients = 0.2 to 0.25 of the critical exit
gradient.
Values of safe exit gradient may be taken as:
0.14 to 0.17 for fine sand
0.17 to 0.20 for coarse sand
0.20 to 0.25 for shingle
H H
C E E C
b1
b
d b d
ii) A straight horizontal floor depressed below the bed, but without any
vertical
cut-off.
171
ii) A straight horizontal floor of negligible thickness with a sheet pile at
some
intermediate point
The mathematical solution of the flow-nets of the above profiles have been given
in the form of curves. From the curves, percentage pressures at various key
points E, C, E1, C1 etc) be determined. The important points to note are:
i) Junctions of pile with the floor on either side{E, C (bottom), E1, C1
(top) }
ii) Bottom point of the pile (D), and
iii) Junction of the bottom corners (D, D’) in case of depressed floor
The percentage pressures at the key points of a simple forms will become valid
for any complex profile, provided the following corrections are effected:
i) correction for mutual interference of piles
ii) correction for the thickness of floor
iii) correction for slope of the floor.
H H
E C
d
1
b1
D d
D b
b
Depressed floor
Sheet pile at intermediate position
172
173
Correction for Mutual Interference of Piles
174
D d + D
C = 19 x --------------------------- ( )
b1 b
This correction is positive when the point is considered to be at the rear of the
interfering pile and negative for points considered in the forward or flow direction
with the interfering pile.
For example, correction for pressure at C1 for pile line (1) by the interference of
pile line (2) is positive as pile line (1) is to the rear of the interfering pile line (2).
Similarly, correction for pressure at E2 for pile line (2) due to the interference of
pile line (1) is negative, because E2 is in the forward or flow direction of
interfering pile line (1).
The interference effect will not be present on the intermediate pile if:
i) the outer pile is equal to or longer than the intermediate pile and
175
ii) if the distance between the intermediate pile and outer interfering
pile is less than twice the length of the outer pile
1 2 3
E1 C1 bs C3
D3
D
d
b1
D1
D2
b D3
E1* C1*
E1 C1
Depressed floor
D1
176
The corrected pressures at E1 should be less than the computed pressure t E1*.
Therefore the correction for the pressure at E1 will be negative. And so also is for
pressure at C1.
Correction for Slope of Floor.
A correction for a sloping impervious floor is positive for the down slope in the
flow direction and negative for the up slope in the direction of flow.
The correction factor must be multiplied by the horizontal length of the slope and
divided by the distance between the two poles between which the sloping floor
exists.
In the diagram above, correction for slope can be applied only to point E2. As the
point E2 is terminating at the descending slope in the direction of flow, the
correction will be positive. The value of correction will be:
C.F. x bs/b1
Where C.F. =correction factor
bs = horizontal length of sloping floor
b1 = horizontal distance between the pile lines
Exit Gradient
b d
177
For the standard form consisting of a floor of a length b, and a vertical cut-off of
depth d, the exit gradient at its downstream end is given by:
1
exit gradient GE = (H/d) x
π λ
By referring to plate (17.3 ), for any value of α = b/d, the corresponding value of
1
can be read off. When H and d are given, GE should be easily calculated.
π λ
The value of exit gradient GE should be within safe limits as given below.
First consider plate 17.1. Values of 1/α = d/b are plotted as abscissa and
percentage pressure φ = P/H x 100 are plotted as ordinates. There are three
curves, one for QD and one for QE for sheet piles fixed at the ends of the floor and
the last one for φD for depressed floor
To find the percentage pressure at the points C1 and D1 of the upstream pile,
calculate
For this value of 1/α determine the value of φD from the curve
Then subtract this value of φD from 100 to get φD1
ΦD1 = 100 – φD
Similarly, determine the value of φE for the value of 1/α. Subtract this φE from
100 to get φC1 φC1 = 100 – φE
Example: Let b = total length of floor = 60.0 m and d = depth of u/s pile = 6.0
Then 1/α = d/b = 6/60 = 0.1
For this value of 1/α, we read the values of φD = 20 and φE = 28.
Then φD1 = 100 – 20 = 80
φC1= 100 – 28 = 72
These percentages pressures must be corrected for mutual interference.
Next, let us consider plate 17.2. This gives pressure curves for φC, φD for sheet
pile not at end condition. The ratios b1/b are plotted as abscissa and φC = Pe/H x
178
100 plotted as ordinates on the left portion of graph and φD = Pe/H x 100 plotted
on the right side of graph for different values of α.
To find φE for any value of α and base ratio b1/b, read φC for base ratio (1-b1/b) for
that value and subtract this value from 100.
Thus for finding φE for base ratio b1/b = 0.3 and α = 3.0
(1 – b1/b) = 1 – 0.3 = 0.7 for α = 3.0
Qc = 20
ΦE = 100 – Qc = 100 – 20 = 80%
ΦD can be calculated for values of b1/b less than 0.5.
Similarly finding φD for b1/b = 0.13 and α = 2.0
(1 – b1/b) = 0.87
Corresponding to 0.87 and α = 2.0 value of φD’ = 40
ΦD = 100 – 40 = 60%
Example. Determine the percentage pressures for various key points shown in
diagram below. Calculate the exit gradient and plot the hydraulic gradient line for
the pond level on the u/s side with no water on downstream side.
▼110
▼107
▼105.0
▼104
▼102.0
0.5m
6 2 25
(1) ▼100
(2)
0.5m
D1 14 6 6
▼99 D2
(3)
45m
60m ▼94.0
D3
Solution:
1. Upstream pile No. 1
Total length of floor b = 60 m
Depth of u/s pile below the top of floor = 105 – 99 = 6.0 m
α = b/d = 60/6 = 10.0
1/α = 0.1
179
From plate 17.1, φC1 = 100 – φE = 100 – 28 = 72 %
ΦD1 = 100-20 = 80 %
For these values, three corrections must be applied .
i) Correction for mutual interference of piles, for pressure at C1
φC1 is affected by the presence of intermediate pile No. (2)
Let D = depth of pile No. (2) = 104.00 – 99.00 = 5.0 m
d = depth of pile No. (1). = 104.00 – 99.00 = 5.0 m
b1 = distance between the two piles = 14.0 m
b = total length of floor = 60.0 m
Then the correction =
D d +D 5 5+ 5
19 . = 19 . = 19 x0.6 x0.1 = 1.9%(+ve)
b1 b 14 60
Point C1 is to the rear of the next pile No. (2). Hence correction is positive.
1.0
E1 C1 E2 C2
▼99.0
D1 D2
φC1 calculated from curves corresponds to point C1! At the top of the floor. But we
want the pressure at C1 at the bottom of the floor. Pressure atC1 will be more
than that at C1!, therefore the correction will be positive.
Pressure at C1! = φC1 = 72%
Pressure at D1 = φD1 = 80%
Distance from C1! To D1 = 1050. – 99.0 = 6.0 m
Hence pressure correction at C1 = (80 – 72) x (105.0-104.0)
(105.0 – 99.0)
= 8/6 x 1 = 1.67 % (+ve)
180
ΦD1 = 80%
Pressure is observed at E2 ! and the direction of flow is from E2! To E2. Hence,
pressure at E2 will be less than that at E2!. Correction is therefore negative.
iii) Correction at E2 due to slope is zero as E2 is neither the beginning nor the end
of a slope.
Hence corrected pressure at E2 = (73% - 1.9% - 1.17%) = 69.93%
181
i) Correction at C2 due to pile interference. Pressure at C2 is affected by
pile line (3) and as C2 is to the rear of pile line (3), correction will be
positive. The amount of this correction is given by:
D D+d
Correction = x
b1 b
Where D = depth of pile No. 3, the effect of which is considered below
the level at which the interference is desired = 104.0 – 99.0 = 10 March
2007
d = depth of pile No. (2), which is affected = 104.0 – 99.0 =
5.0m
b1 = distance between pile (2) and pile (3) = 45m
b = total floor length = 60.0m
10 10 + 5
= 19
45 60
Correction
2 15
= 19 = 2.28%(+ve)
9 60
182
Note that for downstream pile at end, the values of φD and φE are taken directly
from the curves for the required value of 1/α.
a) Correction for φE3
i) Correction due to mutual interference of pile. The point E3 gets affected by
pile No.2. E3 is in the forward direction of flow with reference to pile No. 2,
D D+d
correction is negative. Value of correction = x
b1 b
where D = depth of pile No. 2 which influences point E3 of pile No. 3, the depth
being considered below the level at which interference is desired = 100.0 –
99.0 = 1.0m
d = depth of pile No. 3 which is affected = 100.0 – 94.0 = 6.0m
b1 = distance between the two piles = 45.0m
b = total floor length = 60 m
1.0 1 + 6
therefore correction = = x = 0.262%(−ve)
45 60
iii) Correction due to floor thickness. From the direction of flow of
seepage water, it is clear that pressure at E3 at the bottom of the
floor is less than that at E3!, top of the floor which is the value got by
Khosla’s curves. Hence, correction to be applied for φE3 should be
negative.
Value of this correction = 32% - 22% x 2.00
(102.0-94.0)
= (10 x 2)/8 = 2.5% (-ve)
iv) Correction due to slope. Point E3 is neither starting point nor the
terminus of a slope. Hence, slope correction = 0
Hence, corrected pressure at E3 = φE3=(32%-0.262%-2.5%)=
29.238%
The corrected pressures at all the key points can be shown in the table below
Exit Gradient
183
1
= 0.16
π λ
1
Hence, exit gradient GE = (H/d) x = (8/8) x 0.16 = 0.16
π λ
This exit gradient is safe even for fine sand.
184
Chapter 9
Drain
Drain
Canal Canal
Super passage
Aqueduct
If the bed level of the canal is not very much higher than that of the drainage,
then the cross-drainage work is called siphon aqueduct
Similarly, if the difference between the bed level of the canal and natural drainage
in not very much, the natural drainage being at a higher level, then the cross-
drainage work is called siphon super passage.
The selection of any particular type of cross-drainage work depends on the
following:
1. Bed level of canal and bed level of natural drainage;
2. Discharge of the canal and the drainage
3. The foundation conditions of the site
4. Availability of materials and labour
5. availability of existing modes of communication.
Aqueduct
185
An aqueduct is a structure constructed at the position, where a canal has to
cross a drainage at a sufficiently high level. It is very similar to a bridge but
instead of a roadway or railway, canal water flows above a natural drainage.
Irrigation water may be taken through a pipe over the drainage water if the
section of the irrigation canal is very small. The pipe is supported in position at
the ends by masonry walls and in between by supports
Tail wall
Pipe Head wall
Support
Type I. In this type of aqueduct, the sides of the aqueduct are completely
earthen embankment with full earthen slopes. The length of the culvert through
which the discharge is passing should be adequate to accommodate the water
way of the canal and also the bottom width of the side embankments. In this
case the original canal section is retained and no fluming of canal section is
done.
Type II.
Type II is very much like type one except that the outer section is reduced by
supporting it by masonry or concrete wall.
Type III.
In this type of aqueduct, the canal banks are discontinued at the aqueduct
portion and the canal water is transported in a masonry or concrete trough
constructed above the drainage culvert. In this case the length of the culvert is
reduced because fluming is done. This type becomes necessary where big
canals have to cross big natural drainage with large catchment areas.
Selection of a suitable type. In order to select a suitable type, it is necessary to
understand the following terms:
186
1. Culvert length: is the width of the aqueduct, measured along the drain.
The culvert depends upon the shape and size of the canal section
2. Length of aqueduct: is the length measured perpendicular to the drain. It
is evidently equal to the width of the drain
3. Bank connection: consist of masonry (concrete) wings, etc required to
connect the regular section of the canal to its modified section over the
aqueduct.
In the type I, since the canal is fully in earthen section, the culvert length is
maximum. Hence the cost per unit length of aqueduct will be maximum.
However, bank connections are not require. Hence, the selection of this type
depends on the relative cost of bank connections and that of the aqueduct
proper. In all cases the cost of bank connection is independent of the length of
the aqueduct. \Therefore, the type I is suitable when the length of the aqueduct is
small and the cost of bank connections would be large in comparison to the
saving from the reduction in the width of work if type III were adopted. On the
contrary, in type III, the culvert length is minimum. Hence the cost per unit length
of the aqueduct is minimum. However, the cost of bank connections will be
additional. Therefore, type III is suitable where the length of culvert is large. Type
two is suitable for intermediate conditions.
187
permissible i.e the length between abutments may be increased by 20% to
compensate for contraction of waterway due to the width of the piers.
iii) Velocity of flow through the barrel. This may be taken between 2.0 -3.0 m/s.
Lower velocities may result in silting in the barrel and higher velocities will cause
abrasions of
barrel surface by rolling grit.
Iv0. Height of Opening: The depth of flow in the barrel can be calculated as the
waterway discharge and velocity are fixed. Sufficient margin (free board) should
be provided between the H>F>L and the bottom of canal bed.
Total height of opening = (Depth of flow in barrel + clearance (freeboard)
v). Number of Spans: The total width of the culvert or total length of aqueduct
already been fixed. The number of spans for the culvert should be fixed
considering structural strength and economy of design. More spans may be used
in case of arched culverts.
vi). Contraction of canal waterway: Canal waterway in the aqueduct portion is
reduced in type III by fluming. Fluming ratio is the ration between the restricted
width of canal in the aqueduct portion and the normal width of the canal. This is
taken as ½. However, this reduction in width should not results in the velocity in
the trough exceeding 3m/s or exceeding the critical value bringing the flow to
supercritical flow, which could result in the formation of hydraulic jump in trough –
not desirable.
vii). Length of Contraction: This is the horizontal distance in which the normal
width of the canal is gradually reduced to the contracted width in the aqueduct
portion. A convergence of 2 horizontal to 1 lateral could be assumed.
Viii). Length of expansion: the expansion ration is assumed as 3 horizontal to 1
lateral. This decides the length of the expansion downstream side of the
aqueduct.
ix). Bank connections : these consist of two sets of wing walls for the canal called
canal wings and two sets of wing walls for the drainage called the drainage
wings.
Canal wing walls. These are to be provided on the upstream and downstream
sides of the aqueduct. They protect and retain the earth in the canal banks. The
wings should be constructed on sound foundation.
The Drainage wing are provided on the upstream and downstream sides of the
barrel. They protect and retain the natural slopes of the drain and increase the
seepage path and control the exit gradient. The drainage wings should be taken
deep into the foundation, well below the maximum scour depth. They must be
taken back well into the top of the guide banks.
188
Full supply discharge = 20cumecs; full supply depth = 1.5m; Bed width = 10.0m;
Side slopes = 1.5:1. The canal has to be provided with inspection paths of 3.0 m
width on both the banks. The channel is unlined. Bed level of the channel =
103.00m. Good foundation is available at site. Missing data may be suitably
assumed.
Solution:
The bed level of the canal = 103.0m
The bed level of the drainage = 100.0m
The High Flood Level in the drainage = 1.0m so that level of H.F.L = 101.0
The difference between canal bed level and H.F.L = 103.0 – 101.0 = 2.0m
Hence an aqueduct would have to be designed.
Determination of flood discharge: by Dicken’s formula
Q = C . M3/4
Where C = is a constant = 20; M = catchment area = 20.0km2
TYPE I AQUEDUCT
Road
Canal
H.F.L
DRAIN
CULVERT LENGTH
189
TYPE II AQUEDUCT Masonry wall
Road
Canal
H.F.L
DRAIN
CULVERT LENGTH
190
The pier must be taken 1.5 times the scour depth below the bed of drain,
which is 1.5 x 2.4 = 3.6 m.
therefore total height of pier hp = 3.6 + 1.0 = 4.6 m
s/2
r
B
A
The canal over the drain is flumed with a fluming ratio of 1:2
Hence, bed width of the canal over the drain bc/d = 10 x ½ = 5.0 m
191
The section over the drain will be rectangular and is to be built in R.C.C. The
height of side walls is kept 0.5 m higher than the full supply discharge F.S.D
x) Transitions
Upstream transition is done in 2:1 splay in contraction
Therefore length of contraction transition = (10 – 5)/2 x 2 = 5 m
Upstream transition is done in 3:1 splay expansion
Length of expansion = (10 – 5)/2 x 3 = 7.5 m
xi) An inspection path of 3.0 m width is provided on both sides of the aqueduct.
x = 3/ 2 1 − --------------------- ( )
Bo − B 3f / 2 B x
where Bx = channel width at any section X-X at a distance x from the flumed
section.
L = length of transition
Bo = Bed width of the normal channel section
Bf = bed width of the flumed section.
192
A
B C D
Expansion
Contraction Transition
Transition
Bo Bx Bf
PLAN
A
L B C
v42/2g
v12/2g 2
v2 /2g 2
v3 /2g D
y1 y2 y4
SECTION
Step II. The energy loss due to expansion of section between section CC and
DD may be taken as equal to 0.3(V32-V42)/2g
Step III. Total energy line at section CC = {T.E.L. at section DD + 0.3(V32 – V42) /
2g}
Water surface elevation at section CC, = {T.E.L. at section CC – V32/2g}
Bed level at section BB = {T.E.L at section BB – V22/2g}
Step IV The channel section in the trough from section CC to BB remains
constant. The only loss of head in this section is due to pipe friction, which can
be computed by the Manning’s formula
Q = i/n . A R2/3(I)1/2
Therefore T.E.L at section BB = (T.E.L at section CC + loss of head)
W.S.E at section BB = (T.E.L at section BB – V22/2g)
Bed level at section BB = (W.S.E. at section BB – y2)
193
Step V. The loss of energy between section AA and BB due to contraction
= 0.2(V22 – V12) / 2g
T.E.L at section AA = [T.E.L. at section BB + 0.2(V22 – V12) / 2g]
W.S.E at section AA = [T.E.L at section AA – V12/2g]
Bed level at section AA = W.S.E. – y1]
Step VI After calculating the T.E.L and W.S.E and bed levels of all sections, the
total energy line may be drawn and the bed line also.
Step VII. The drop in the W.S.L. between the two adjacent sections is due to:
i)drop in the energy line between the two sections’
ii) increase velocity head for contraction and decreased velocity head for
expansion.
Step VIII) When the water surface profile has been plotted over the whole length,
the velocity head can be determined by measuring the vertical distance between
the T.E.L and the water surface line at any point.
The velocity head can be converted into equivalent velocity by V = (2gh)1/2
Step IX. Discharge at any point,
Q = A.V.
With the velocity known, and the flow also available, the dimensions of the
section can easily be calculated.
If the section is rectangular, then A = b x h
If the section is trapezoidal, A = (b + mh)h, where m = side slope =cot θ; θ =
angle of drain.
SIPHON AQUEDUCT
In siphon aqueduct, the difference between the bed level of the drainage and the
bed level of the canal is not much so the bed level of the drainage is depressed
at the site of the crossing so that there is sufficient clearance between H.F.L of
the drainage and the bed of the canal. The drainage water passes below the
canal through the depressed portion which makes it works like siphon.
L V 2 Va2
h = 1 + f 1 + f 2 . . −
R 2g 2g
194
Retaining walls
Canal banks u/s H.F.L
d/s H.F.L.
F.S.L
Arch thickness
F.S.L
d/s H.F.L u/s
After fixing the velocity through the barrel of the siphon the head ‘h’ required to
generate that much velocity can be determined by equation above. As the
downstream H.F.L of the drain remains unchanged, the upstream H.F.L. can be
obtained by adding ‘h’ to the downstream H.F.L. ‘h’ is known as the afflux.
195
Since the barrel is depressed below the bed level of the drainage, it will be below
the normal water table in the surrounding area. Due to this, a static uplift
pressure will be exerted on the barrel.
Solution:
The flow is large and therefore the drainage size is assumed large and therefore
we choose type III aqueduct which might be economical. The canal at the
aqueduct shall be flumed and taken in concrete trough.
196
= (22 – 12)/2 x 2 = 10 m
Providing a splay of 3 : 1 in the expansion the length of the expansion transition
= (22 – 12) /2 x 3 = 15 m
Length of the flumed portion from abutment to abutment
= total drainage water-way = 94.5 m
Design of Levels at Various Sections {refer to diagram on page (8)}
a) Section DD. Cross-section area A = B + 1.5h)h = (22+1.5x1.5)x1.5
=36.75 m2
Canal velocity = Q/A = 30/36.75 = 0.82m/s
Velocity head = V2/2g = (0.82)2/2x9.81 = 0.034m
Therefore Relative level (R.L) of water surface
= R.L of bed + water depth = 94.5 +1.5 = 96.0 m
Relative level of the Total Energy Line (T.E.L) = 96.0 + 0.034 = 96.034 m
b) Section CC. Cross section area of canal = 12 x 1.5 = 18 m2
Velocity = Q/A 30/18 = 1.67 m/s
Velocity head = V2/2g = (1.67)2/2 x 9.81 = 0.142 m
Loss of head in expansion from section CC to DD
= 0.3(V32 – V42)/2g = 0.3(0.142 – 0.034) = 0.032
Therefore Level of T.E.L at section CC = T.E.L at section DD + Head loss
= 96.034 + 0.032 = 96.066 m
R.L of water surface at CC = T.E.L – velocity head = (96.066 –
0.142)=95.924 m
Therefore R.L. of canal bed for constant depth,
= R.L. of water surface – water depth = (95.924 – 1.5) = 94.424 m
197
= R.L of T.E.L at BB + head loss = (96.11 + 0.0215) = 96.325m
R.L. of water surface
= R. L. of T.E.L – V12/2g = 96.325 – (0.82)2/2 x 9.81= (96.325 – 0.034) = 96.291 m
Design of Contraction Transition: The design of the transition will be done on the
basis of Chaturvedi’s formula
L.Bo3 / 2 B f
3/ 2
x = 3/ 2 1 −
Bo − B 3f / 2 B x
here Bo = 22 m; L = 10m; Bf = 12 m
10.(22) 3 / 2 12 3 / 2 12 3 / 2
or x= 1 − = 16.741 −
(22) 3 / 2 − (12) 3 / 2 B x B x
The values of x for various values of Bx are calculated from the above equation,
and are tabulated below. The distance x is measured from flumed section BB as
shown in page 8.
Bx 12 13 14 15 16 17 18 19 20 21 22
x 0 1.89 3.38 4.76 5.86 6.81 7.62 8.34 8.96 9.5 10
L.Bo3 / 2 B f
3/ 2
x = 3/ 2 1 −
Bo − B 3f / 2 B x
here, L = 15m; Bf =12; B0 =22
15 x(22) 3 / 2 12
3/ 2
x= 1 −
(22) 3 / 2 − (12) 3 / 2 B x
12 3/ 2
x = 25.111 −
B x
The values of x, for various values of Bx are calculated in the table below.
Bx 12 13 14 15 16 17 18 19 20 21 22
x 0 2.84 5.18 7.14 8.8 10.21 11.44 12.5 13.44 14.26 15
198
The trough shall be divided into two equal compartments by 30 cm thick concrete
wall. The road way 6.0 m wide shall be carried over one of the compartments
The height of trough = 1.5 + 0.5 m free board = 2.0 m.
L V 2
h = 1 + f 1 + f 2 .
R 2g
V = velocity through barrel ( = 1.99m/sec)
f1 = coefficient of loss of head at entry ( 0.505)
f2 = a(1 + b/R) a and b are picked from table given above.
a = 0.00316 and b = 0.030
199
Assuming the subsoil water level is up to the bed level of the drain, the static
uplift on the floor, will be = (92.5 -90.5) = 2.0 m
To be continued on 475
200
Chapter 10
DESIGN OF CHANNELS
Aim of channel design: is to arrive at a cross-section, which can carry the design flow
without either scouring or silting problems.
Design Parameters
For the design of channels, the following data shall be available:
i) Design discharge Q; ii) Surface and soil properties (roughness coefficient); silt factor f
The design: consists of determining the following four factors:
i) Area of cross-section, A of the channel
ii) Hydraulic mean depth or hydraulic radius R
iii) Velocity of flow, v
iv) Longitudinal slope of the bed, S.
To begin the design, the following two equations are considered (available)
a) Q = Av - continuity equation
b) V = f(n, R, S) – Flow equation, Manning, Chezy, Kutter.
In the above, there are four unknowns and therefore to solve them, we require two
additional equations.
To get the additional two equations, the following procedures may be adopted
i) By providing a channel of best discharging section, an equation between area
A and hydraulic mean depth may be obtained.
ii) Based on consideration of scouring and silting, a limiting equation for velocity
of flow may be obtained.
iii) Fixing the longitudinal slope based on available ground slope.
iv) Based on experience, fixing a suitable width-depth (B/D) ratio not necessary
that which gives best section.
1) To keep silt in suspension, the silt supporting power is directly proportional to the
bed width of the stream and does not depend on the wetted perimeter.
2) The limiting velocity which does not produce silting or scouring may be called
critical velocity. The relation between critical velocity and depth is
vc = 0.55xD0.64
And after modification we will arrive
v = 0.55 x m x D0.65 where m = v/vc
201
v - velocity of flow in channel, m = critical velocity ratio (C.V.R.), m varies from 0.7
for fine sand to 1.3 for very coarse sand
To determine the mean velocity of the flow in the channel, Kennefdy used the Kutter’s
formula: v = C (RS)1/2 where
1 0.00155
23 + +
C= n S
0.00155 n
1 + 23 + S
R
The difficulty with Kennedy’s formula is that he did not give any formula to fix the slope
S but has to be fixed based on the natural ground slope which makes it difficult to obtain
the best section for the given discharge.
Worked example
202
v The hydraulic radius R = A/P = 65/35.97 = 1.806m
1 0.00155
23 + 0.0225 + 1
5000
The Kutter’s constant C = = 49.63
0.00155 0.0225
1
+ 23 + 1.806
1
5000
1 0.00155
23 + 0.0225 + 1
5000
C = = 51.05
0.00155 0.0225
1 + 23 + 1 2.14
5000
The computation that were done previously could be done by using Garret’s
diagrams, which give a graphical solution for Kennedy’s and Kutter’s equation.
Garret’s diagram has discharge plotted as abscissa. The ordinates on the left side
indicate the slope and those on the right side, the depths in the channel and critical
velocity vc. Bed width lines are shown as dotted.
The procedure for using the Garret’s diagrams is as follows:
A) Follow the discharge curve for a given discharge and note its point of intersection
with the given slope line (which are horizontal).
203
B) From this point of intersection, draw a vertical line. This will intersect many bed
width curves.
C) For any bed width, consider the point of intersection of the vertical line and read
out the correspoiding water depth D and critical velocity vc on the right end of
diagram.
D) Calculate the velocity of flow v for the assumed bed width and its corresponding
water depth.
E) Determine the critical velocity ratio (C.V.R = v/vc)
F) The value of C.V.R should be the same as one given. Otherwise, the procedure
must be repeated with other values of bed widths got by the intersection of the
vertical on the bed width curves for the given slope, till the value of C.V.R. is very
near the given value
Garret’s diagram are drawn for side slopes of ½:1. The curves can be used for any
other value of roughness n by making use of the small nomogram provided at the top
of curves. When any other value of n is to be made use of, see whether the value is to
the left or right of the arrow shown on the nomogram and note the shift. When the
vertical line is drawn through the intersection of the slope line and discharge curve, it
is with respect to n = 0.0225 which corresponds to the arrow on the nomogram. The
vertical line has to be shifted either to the left or to the right by the same amount as
the required value of n is shifted on the nomogram, and then the intersection points on
the bed width curves considered.
Solution:
Refer to the Garret’s diagram. For a slope 1/5000, it corresponds to 0.2 as the slope is
marked in 1/1000 units.
The horizontal from slope 0.2 cuts the discharge curve for 2.5 m3/s at some point.
From that intersection point, draw a vertical line to cut many bed widths curves. Let
us consider the bed widths curves of b = 6.0; b = 7.0m; b = 8.0m; b = 9.0m;
b = 10m. From these points of intersection, horizontal lines to the right may be
drawn and the corresponding values “D” and vc read out. The rest of the computation
are done in the table below.
204
From the table, we see that No. 3 gives very near value to the one given. Hence
channel dimensions should br: B = 8.0m; D = 0.67m
205
LACEY’S REGIME THEORY
Since all the above conditions cannot usually be satisfied by artificially constructed
channel, it follows that it cannot be in “true regime”. Either in initial or final regime
206
and finer, the shape of the section becomes narrower. Therefore when a channel is
constructed with small cross-section, and steeper slope than necessary, scour starts and
continues till final regime is established. Similarly, in a channel with a bigger cross-
section and flatter slope, silting starts and continues till final regime is established.
207
By using equation (L-1), (L-2) and (l-7) ie the fundamental Lacey’s equations, we shall
obtain
S = (f3/2)/(4980R1/2) ---------------(L-8)
S = (f5/3)/(3340Q1/6) ---------------(L-10)
v = (2/5 xf x R)1/2
v2 = 2/5 x f x R
R = 5/2 (v2/f) -----------------------(L-11)
208
Design of Channels by Lacey’s Theory
Example: Design a regime channel for a discharge of 50m3/s, with silt factor f = 1.0
using Lacey’s theory.
209
Design of channel by use of Lacey’s diagram
Example: determine the channel dimensions and slope for a channel with a discharge of
53m3/s and silt factor of 1.0 using Lacey’s diagram
Referring to plate 16.5, there are discharge curves for 50 m3/s and 55m3/s. By
interpolation, the discharge curve for 53 m3/s can be drawn. This curve intersects the
curve for silt factor f = 1.0 at some point. For that point, the abscissa will be 30.0m. This
gives the bed width b = 30.0m
The ordinate for the intersection point = 2.0m. This gives the depth.
Now to determine the regime slope, refer to plate 16.6. Get the intersection point of the
curve for silt factor f = 1.0 and the discharge curve for Q = 53 m3/s. The ordinate of the
intersection point represents the slope.
The ordinate in this case = 0.15 representing a slope of 0.15 in 1000. t
Therefore S = 0.15/1000 = 1/ 6667
210