Discharge Relations For Hydraulic Structures and Head Losses From Different Components

QI
Koikman, P.A.
Discharge relations for
hydraulic structures an.
IL-uIh/ iiiiU
Discharge relations for hydraulic

structures and head losses from
different components
Desk study for Rijkswaterstaat and

the International Association for Hydraulic Research,
performed in cooperation with
the Delft University of Technology
C67
BSW-CW (
BOtJ WSPEUR WERK
CONSTRUCTIES EN WATER
December 1989
*
Z2331
Q 953 delft hydraulics

Preface
This report "Discharge Relations of Hydraulic Structures and Head Losses from
Different Components" is conceived as part III of the future IAHR manual about
discharge relations and losses in general. Part 1, which is presented by
Dr. D. Miller will deal with losses in internal flow systems, and part II
(Dr. W.M. Hager) will deal with the design of discharge measuring structures.
The content of this part of the rnanual is purely based on literature. It is
providing designers both background information and easily applicable data.
There are much data on discharge relations available and so the selection of
topics and literature was one of the problems the author had to cope with and
sometimes arbitrary decisions had to be taken. The problem was also the non-
uniformity of symhols, and matters would have been easier when not at some
time the water level and at another time the energy level had been introduced
as a reference for the hydraulic condition. It was considered to be unprac-
tical to redesign all figures, but for all the figures the meaning of the
applied symbols is given.
It was inevitable that once again use was made of the design manuals of the
Vickburg Waterways Experiment Station of the US Army and the US Bureau of
Reclamation, the data of which are used worldwide because of their concise and
user-friendly presentation. When presented here, however, all results were
transormed into a dimensionless form. 1 could also make use of the open-
channel-flow textbook by prof. E. Naudascher which fortunately had just been
published.
This part III is certainly not sufficiently complete, to cover all the needs
for design, but while also being a general introduction, a certain conciseness
was aimed at. Only after a certain period of use and after receiving comments
from users gaps may be filled.
1 got help from a few colleagues of DELFT HYDRAULICS and RIJKSWATERSTAAT who
corrected initial errors and indicated omissions. In particular 1 wish to
thank prof. J. Battjes of Delft University of Technology for his close reading
of the equations, symbols and text of the first draft, and prof. P. Novak of
New Castle University for his corrections and proposed additions.
RIJKSWATERSTAAT financed, within the framework of the Hydraulic Structures
Research (BSW-CW), the final editing of this report. This edition has been
published in a limited number of copies, as a DELFT HYDRAULICS' report Q953.
December 1989
P.A. Kolkman
CONTENTS
List of Symbols
GENERAL ASPECTS OF DISCHARGE RELATIONS AND LOSSES ..................1
1.1 Introduction ..................................................1

1.2 Discharge relations of hydraulic structures introduced in far-
field computations ............................................6
1.3 The use of the Bernoulli-, the momentum- and the Carnot equations 8
1.4 Some general remarks before starting the calculation of losses 10
1.5 Modular flow equations ........................................11
1.6 Fully submerged flow ..........................................15
1.7 The intermediate-flow relations at control sectioris ...........19
1.8 Remarks about entrance losses .................................23
1.9 Remarks about trash-rack and bridge pier losses ...............25
1.10 Required and obtainable accuracies ............................27
NUMERICAL PROCEDtffiE FOR COMPUTATION OF WATER LEVEL/DISCHARGE

RELATIONS ..........................................................29
2.1 Introduction ..................................................29

2.2 The search procedure ..........................................31
Appendix of Chapter 2: Iterative search procedure for h 2 using the
ACORN-BBC computer.
DISCHARGE RELATIONS FOR CONTROL SECTIONS ...........................36
3.1 The sharp-crested weir ........................................36

3.2 The weir with circular crest section ..........................42
3.3 The nappe-shaped overflow weir ................................46
Amodular flow.................................................47
B submerged flow...............................................48
C modular flow with piers .....................................49
3.4 The nappe-shaped spillway with gate ...........................54
3.5 Broadcrested weirs and dikes ..................................56
3.6 The morning glory spillway ....................................62
CONTENTS (continued)
3.7 Labyrinth weir 69

3.8 Gates with underflow ..........................................73
3.9 Gates with overfiow ...........................................86
3.10 Gates with combined under- and overfiow .......................90
3.11 The Howell Bunger valve .......................................92
3.12 Conduit outlet ................................................96
4. LOSSES OF COMPONENTS UP- AND DOWNSTREAM OF THE CONTROL SECTION .....97
4.1 Orifice and intake losses .....................................97

4.2 Friction losses ...............................................105
4.3 Slots .........................................................109
4.4 Exit losses and hydraulic jump relations ......................112
A abrupt expansions downstream of control section.............114
B hydraulic jump relations (two-dimensional) ..................116
C hydraulic jump on a slope ...................................120
Dgradual expansions ..........................................122
E hydraulic jump in a gradual expansion.......................125
4.5 Trash racks ...................................................127
4.6 Resistance of bridge piers ....................................132
REFERENCES
LIST OF SYMBOLS
a (gate) opening height (m)

a and b parts of the nappe at weir flow (m)
A culvert area (m 2 )
A orrifice area (ni 2 )
Cb loss coefficient by bridge piers

C contraction coefficient
c
Cd discharge coefficient
CD drag coefficient of bar or pier at internal flow or
at free flow with low velocities, without blockage effect
Cf friction factor
C0 discharge coefficient at free flow (used also as a reference
discharge coefficient)
C reduction factor by submergence
C loss coefficient of trash rack
tr
d gate thickness, diameter of semi-circular gate edge, waterdepth (m)
jet thickness (m)
upstream waterdepth (m)
waterdepth in or just downstream of control section
d1 waterdepth upstream of hydraulic jump (m)
downstream waterdepth (m)
d critical depth = (q2/g)1'3 (m)
d waterdepth in culvert (m)
d waterdepth at submerged supercritical flow (m)
D pipe diameter, height of roofed culvert (m)
Dh hydraulic diameter (4 x hydraulic radius) (in)

F Froude number v/Ij (F 2 refers to section 2 etc)
F1 Froude number upstream of hydraulic jump
g acceleration of gravity (ms- 2 )
G0 height of gate opening (m)

h upstream water level related to crest level at modular flow (m)
h0 upstream head related to weir crest or culvert bottom (m)
h1 downstream head in or just downstream of control section (m)
h2 downstream head further downstream (m)
H upstream energy head (above weir crest) at modular flow (m)
upstream energy head above weir crest or culvert bottom (m)
H energy head in culvert (m)
crest-design value of upstream energy head (m)
Hd
H (effective) energy head upstream of a weir (m)
LIST OF SYMBOLS (continued)
K coefficient in the bridge-resistance formula of Yarnell

K loss coefficient of slots
S
K abutment coefficient
K pier flow-contraction coefficient
L length of crest in flow direction (m)
L bruto length of labyrinth weir (m)
L lip length of gate (m)
L length of hydraulic jump (m)
Lr length of radial JUMP (m)

m 1 , m 2 discharge coefficients
p weir height (m)
P pressure (Pa)
Pa atmospheric pressure (Pa)
q discharge per unit width (m 2 s' )
Q discharge
(m3s1 )
Qg gate discharge (m 3 s' )

Q l real discharge over labyrinth weir (m 3 s- 1 )
Q discharge at a sharp-crested straight weir (m 3 s- 1 )
Q 0 discharge at modular flow (used as reference discharge) (m 3 s- 1 )

r ratio r /r
o 21
r 1 , r 2 radii up- and downstream of radial jump (m)
R bottom step (m)
Rh hydraulic radius (m)

S submergence factor (h 1 /h0 ) or (h 2 /h0 )
S 0 slope
IJ velocity of flow (m S_1)
(m s1)
average pipe flow velocity
(m s1)
V flow velocity
Vc culvert velocity (m s ')
V0 orrifice velocity (m s ')
w width of opening, weir, sluice, culvert or gate (m)
w repeating distance at labyrinth weir (m)
w bar thickness of trash rack (m)
W canal width (m)
W distance between trash rack bars (m)
W Weber number
e
y waterdepth on a slope, measured vertically (in)
Y0 ratio d2/d1
LIST OF SYMBOLS (continued)
coordinate in vertical direction (m)

z
correction factor in the Bernoulli equation for uneven flow
distribution

a divergence angle at one side of diffusor (degrees)
correction factor in the momentum equation for uneven flow
distribution

& boundary layer thickness (m)
displacement thickness of boundary layer (m)
6d
shape coefficient of bridge piers in Rehbock formula
(m)
AH loss of energy head
t

wall shear stress (kg m- 1 s- 2 )
e divergence angle (= 2 a) (rad)

coefficient for loss of energy head
related to V2/2g or to (V 1 - V 2 )2/2g
E e entrance loss coefficient etc.
(kg m- 3
p fluid density )
rate of submergence of submerged hydraulic jump (= d/d1)

1
1. GENERAL ASPECTS OF DISCHARGE RELATIONS AND LOSSES
1.1 Introduction
Hydraulic losses are among the classic research topics in laboratories for
hydraulic research, and a continuous flow of publications has been the
result.
In contrast to publications about internal flow, only a limited number of
publications about losses and discharge relations of free surface flow
structures are of more general value. Because at free surface flow many
more parameters (compared to internal flow) are involved, the chance is
small that exactly similar conditions occur as compared to data from
1 iterature.
Because internal flow is discussed in Part T of the IAHR manual, in the

following only data and computation methods are presented which are related
to free surface flow. The application of scale models is not further
discussed.
The data of this manual will mainly be used for estimating in the pre-
design stage the dimensions of a structure through which a certain
discharge should pass at a given combination of up- and downstream water
levels (or difference of water levels).
Another use of these data can be the introduction of structures such as
sluices, weirs, and so on in far-field computational models of canal
networks and estuaries. Then it is important to have insight into the
nature of the discharge relations with all sorts of combinations of the
upstream and downstream water levels.
The following type of equation is called a "discharge relation"

Q = f(h, h geometry) (la)
2)
in which Q = discharge, h 0 = upstream water level, and h 2 = downstream
water level (in general at a distance where the flow is redistributed over
the flow section).
Instead of water levels, also the energy heads can be applied.

Q = f(H 0 , H 2 , geometry) (ib)
When in a structure one section controls either in total or mainly the

t,
discharge relation, then this is called a control section . In Chapter 3
the discharge relations of control sections are presented.
FJ
When the structure has a certain length (in flow direction) the "hydraulic
losses from different components" are also involved, such as entrance- and
exit-losses, friction losses, losses at a hydraulic jump etc.
A number of these losses are expressed in terms of losses of energy head.
As and example the entrance loss:
AH E V2/2 g (2a)
e e c
where:
AH = loss of energy head at the entrance,
= entrance loss coefficient,
V = culvert velocity and
g acceleration of gravity.
At a free-surface condition this equation becomes:
(Q/wd)2/2g (2b)
AH = e
(w = culvert width, d = water depth in the culvert).
This relation differs markedly from Equation 1, and in equation 2b the

depth d is coupled to both the up- and the downstream head. The following
procedure is suitable to determine the discharge relation of the structure.
The discharge relation of a control section, mostly expressed in the form
of Equation (1), is transformed into:
h = f(Q, h 1 )
(h 1 is the head just downstream of the control section and h is the water
level just upstream of this section).
Now to h is added the effect of the different losses from elements
upstream and from h 1 is substracted the effect of the downstream losses
which results in h 2 . (All these additional losses are found from
expressions similar to equation 2b but then transformed into steps in water
level).
When these calculations are performed for a sufficient number of
combinations of Q, h 0 and h 1 then a new discharge relation of the type of
Equation 1 can be established.
3
In this manual data are presented about discharge relations of control

sections and about losses of components. Discharge relations of complete
hydraulic structure are not presented here as these can only be considered
as a special case. In certain cases the discharge relations of the control
section might be representative for the whole structure. This applies when
the flow at the control section is strongly throttled (discharge over
weirs, over and under gates etc.).
Concerning the discharge characteristics of a control section with free

surface flow, there is already a great complexity: three flow regimes are
distinguished. Later this wil be illustrated in Figures 3 and 4, in which
configurations with openings and overtopping are introduced.
free flow (or modular flow), where the discharge is related to the
upstream head only; Q = f(h 0 ). Because the upstream energy head H 0
(being (h 0 + V/2g)) is related to h 0 and Q, also Q = f(H 0 ) is a unique
relation. In the design it is of interest to know how low the
downstream water level must be to guarantee modular flow.
fully submerged flow, here "fully" means with small head differences.
Fully submerged flow is comparable with internal flow because the
nearly horizontal water level is a boundary which remains independent
of the flow velocity (as long as this velocity is small), so the loss
and the discharge relation are also AH = V 2 /2g, and Q = C d A 12gAH,
with the loss coefficient E and the discharge coefficient Cd being
independent of AH. AH is the loss of energy head from up- to downstream
of the control section. When the flow is overtopping a sill or a gate,
then the flow section "A" depends on the water level.
C. intermediate flow, in between the two former regimes; this is the most
complicated one; the discharge depends on both energy levels H 0 and
(or on both heads h 0 and h 1 ). It comprises the situation of what is
called submerged weir flow, wherein AH determines the discharge but
where the flow section is mainly determined by the downstream level
instead of the upstream level (like Figure 4b). It comprises also the
situation of semi-submerged gate outflow with an eddy on top of the jet
(with thickness d.) with super-critical flow (Froude number > 1 or
V > Id.). And as will shown in section 1.7 (Figure 5) there can also
be mixed flow, where the upper part of the weir flow behaves like free
flow and the lower part like submerged flow.
rd'
The transition between the free flow and fully submerged regimes and the
intermediate regime cannot precisely be predicted.
In Figure 1, a picture is shown of the generalized discharge/water-level

relation, where the three flow regimes are presented. Such a presentation
is valid for the control section, but can also be drawn up for the complete
structure.
te
ter
t modular flow
H 0 =3
log Q Q = f(H)
H 0 =2
H0 =1
\\ \overflow
Q.
underfiow
10 1ogH
Figure 1 Generalized discharge/water-level relation for one geometry (and

also one gate position).
In case the complete structure is considered, the characteristic property

of the fully submerged flow regime remains that discharge Q increases in
proportion to the square root of the head difference. At overflow the
control opening and the discharge vary also with the waterlevel h 0 .
The characteristic property of the free (or modular) flow is again that (at
a given upstream water level h ) the discharge Q is completely independent
0
of the head difference. The relation with the upstream water level can be
in proportion with H312 (overfiow) but also with H (underflow or a
horizontal slit halfway) or a combination of both.
At intermediate flow the characteristics of the total structure are more
variable than the one of the control section.
When the head djfference increases, a dip is seen in the water level just
downstream of the control section. This is related to the recovery of
potential energy further downstream. So the head difference over the
control section is greater than the global head difference.
5
Therefore, if the flow opening remains constant, the discharge increases

more than proportional with the square root of the head difference. This
for instance is observed in Figure 42.
At overfiow the flow section decreases with a lower downstream water level
and this effect will generally be stronger than the effect of the increased
head difference as was mentioned above.
The intermediate flow regime can also show a mixed flow condition at a
number of openings, where some of them have modular flow and some submerged
flow.
Figure 1 applies to one geometry only: at a structure with a gate such a

set of curves should be set up for each gate position. Generally, this is
not done but the result is that simplified presentations show scattered
data!
Resuits of computation of losses in free-surface flow at hydraulic

structures are less precise than those for internal flow. The number of
parameters requires schematizations and often a scale model is used for
final calibration.
Computations serve mainly for dimensioning the structure in the pre-design
stage. When the structure is meant for discharge metering it has to be
specially designed for that purpose (standard design weir), or the
structure must have a reference section with a good flow distribution,
where the discharge is calibrated with measurements of velocities at
different points. With an accoustic discharge meter the flow distribution
can be uneven and, depending on the type of equipment, streamlines need not
be parralel.
In Section 2 js introduced a computer algorithrn which can be used to handle

the type of equations for free surface flow in the intermediate flow
regime. For free or fully-submerged flow, data can be used which are
presented in part II and part 1 in this manual (measuring weirs and
internal flow respectively). Left out are special situations with, for
instance, mixed air/water flow, and the transition of roofed culvert parts
with subroerged and with free flow.
1.2 Discharge relations of hydraulic structures introduced in far-field
computat ions
In one-dimensional far-field computational models (networks), and in two-

dimensional tidal and flow models, the discharge relation of structures,
dikes, and so on needs to be represented.
Suppose that the type of curves of Figure 1 are determined either from
scale-model measurements or from computations. In both cases there is some
difficulty in incorporating the resuits in network or in tidal computation
schemes. The nature of the curves is for each upper water level:
1f Ah (or AH) > P, with P = f(h0 ) then modular flow: Q = f 1 (h 0 ) (or f(H 0 ))
1f Ah < R, with R = f(h 0 ) then submerged flow: Q = f 2 (h 0 ) 1 2gAh
In dynamic networks or tidal computations one prefers to use discharge

relations for a barrier or a closure gap, which can be differentiated
continuously. Expressions with "if" statements (Boolean variables) are less
suitable. The discharge relation can be replaced by an empirical expression
approaching reality in a limited range of hydraulic conditions.
The following type of expression might be suitable to represent the

complete weir flow-discharge relations (which at their turn are obtained by
empiry or by computation) in a versatile way:
1 1 -11fl
Q = [ + (3)
{f2(h0) I2gAh)n [f1(h0)]
is representative for the fully submerged flow condition (where Ah is

small) and f 1 for the modular flow.
In this expression f 1 , f 2 and n can stili be functions of h 0 or Ah.

It can be seen that, for one h 0 value, the left term dominates when Ah is
small and the right one dominates when Ah is great, so similar lines as in
Figure 1 are obtained.
The power n is introduced to describe the curvature between the modular
flow and the submerged flow line in Figure 1.
At gate underflow the intermediate flow does not show a simple sharp
transition between the submerged and modular flow regimes, for which a not-
too-small n must be chosen,but the relation Q 1K1 which is valid for the
7
fully submerged regime, gets steeper in the intermediate flow regime. So a

function as
Q = a + OAH + y AHb (with b > a > 0.5)
should replace the former relation Q IKi. The power n in Equation 3 must
be higher than a or b.
In two-dimensional far-field flow calculations the skew overfiow over dikes

causes a modification in the direction of momentum. This is discussed by
Schnfeld [35]. When the dike is supposed to be smooth, the momentum in the
direction of the width (say "perpendicular to the flow") does not change.
The momentum in flow direction modifies according to what is found for two-
dimensional flow. Now the initial upstream velocity (to be used for
obtaining the energy head) must only be the velocity component in
longitudinal direction of the crest.
1.3 The use of the Bernoulli-, the momentum- and the Carnot equations
For free-surface flow the Bernoulli- and momentum equations can be applied
in the same way as for internal flow. The Bernoulli equation applies only
when the flow converges (accelerates).
The momentum equation can be applied when the pressure distribution is
known. The momentum equation is then suited to calculate the losses in the
zone of decelerating flow. The application of the momentuin equation is of
great value in all those cases where the flow profile is suddenly widening
and where in zone of the separated jet and the eddy on top or underneath
reigns a nearly hydrostatic pressure.
As for the momentum equation: for internal flow (pipe flow) it is common to
use the Carnot equation at abrupt expansions with cylindrical or
rectangular culverts. This equation, which is derived from the momentum
equation, is easier to handle than the momentum equation itself. With free-
surface flow, however, the Carnot equation only applies at low velocity
conditions where the water level is almost straight and parallel to the
bottom. Therefore, in general, with currents having free water surface the
mornentum equation must be applied directly.
Most computations of pressures and losses are based on one-dimensional

considerations, where discharge is supposed to be evenly distributed over
the flow section.
1.12
Bernoulli in l-d computations
1.10 h a
and 1. 1 2g - h2 + 0 2
1.08 i r v'd A
J
A
1.06 omentuzn eq
t Force pQ( 2 V
- l v)
1.04
1 f vd A
19 B
AJ -
A
1.02
definitions
VA
=average velocity over
1.00 section A
C f (p V) = shear Stress
06 .05 .04 .03 .02 .01 0
Figure 2 Correction coefficient a and 0 for a fully developed velocity

distribution as function of wail-friction factor; from
Naudascher [24]
But when the velocity profile is non-uniform the Bernoulli equation cannot
be applied automatically. The transport of energy, proportional to p V 3 ,
must be expressed correctly in the energy balance.

A (sometimes necessary) refinement can be obtained when in the Bernoulli
equation a factor a is introduced into the /2g term. The now refers to
the average velocity over the Section A. The a factor is always greater
than 1. A similar procedure is used in the momentum equation (a factor
applies).
In Figure 2 a and 0 factors and their definitions are presented (from
Naudascher [24], but referring to the work of Rouse) for those conditions
where the velocity distribution is fully developed and adapted to the wall
friction c f .
In a streamlined flow contraction of short length one should apply
a = 0 = 1 because the boundary layers have not yet developed and the flow
is still evenly distributed.
10
1.4 Some general remarks before starting the calculation of losses
In general an outlet work or a sluice must be approached as if they

were a chain of hydraulic losses. A dominant element is the discharge
relation of the control section. This control section is a sill for
instance, or an opening or a narrow flow section. First the discharge
relation of the control section must be found.
When first the modular flow is considered, then only the extra losses
at the upstream canal section play a role, resulting in an additional
rise of the upstream level.
There is a maximum downstream level (close to the control section)

where the modular flow is not affected. Naturally the corresponding
water level farther away depends on losses, change in flow section, and
50 on. It is important to consider separately the water level just
downstream of the structure and the water level farther away where the
flow is distributed over the whole flow section.
Discharge-regulating structures are relatively short. The losses due to

the inlet, slots, shafts, and so on and, in fact, are therefore not
independent of each other. In computated these interactions are
neglected, and hence only a limited accuracy of computed resuits is
obtairied. The importance of interactiori effects at non-modular flow
increases when modular flow is nearly reached: diving jets and other
phenomena cause deviations from the one-dimensional flow.
For submerged flow or for the intermediate flow regime, first the
discharge/water-level relation should also be established for the
control section. Then the additional losses up- and downstream must be
added. The momentum equation can be applied just downstream of this
section, but only for abrupt expansiQns where sufficient information
about the pressure distribution is available. For other configurations
(scarcely available) empirical data have to be used.
11
Complicating factors in the analysis can be points F and G:
F. There can exist conditions at which it is not dear where with modular
flow the control section is located. A flow c ntraction at a gate
recess can suddenly create a control section, and even situations where
several critical sections occur are possible. The trash-rack can create
such a flow contraction that at that location a control section will be
formed resulting in extremely high (and often unexpected) losses.
G. Near the modular flow condition hysteresis effects can occur; another
discharge is found depending on whether a certain downstream water
level is reached by raising or by lowering. With an overflow coridition
starting at a low downstream level the free-falling jet dives and
causes a local dip in the downstream water level, which in turn results
in sustaining the modular flow condition with a diving jet. However,
starting at a high downstream water level the submerged jet does not
tend to dive and hence no local dip in the downstream water level
occurs.
The size of hydraulic.structures is large, so that the influence of the

water viscosity (Reynolds number) is negligible. This does not apply for
the scale-model results.
12
1.5 Modular flow eguations
Figure 3 shows examples of modular flow conditions and their discharge

equations. Approximate atmospheric pressure will reign in the whole flow
section of the vena contracta when sufficient side- and/or under-aeration
occurs. Reality is somewhat different but in the applied equations an
empirical coefficient will incorporate these deviations. Near-hydrostatic
pressure occurs in the vena contracta with two-dimensional flow above a
horizontal bottom or a long-crested sill.
In the case of the vertical partial slit the conditions of atmospheric
pressure and the (assumed) horizontal outfiow in the vena contracta lead to
the following expression for the discharge relation (see Figure 3 case a)
Q = C w (2/3)(Ij) [H 3 / 2 - ( H - a) 3 / 2 } (4)
c 0 0
and when the slit is elongated till the upstream water level (case b) this
equation transforms into
Q = C w (2/3)(Ij) H 3 / 2 (5)
c 0
Here the contraction coefficient is introduced as C the slit width as w

cy
the slit height as a, while H presents the upstream energy head in
relation to the bottom of the slit.
In reality the contraction coefficient has to be replaced by a discharge

coefficient Cd now also including the disturbing effects of vertical
velocity components of vertical contraction.
For the overflow nappe (Figure 3 case c) which is fully aerated, a similar
expression applies but although at the top and bottom atmospheric
conditions exist, the pressure inside of Section A-A will deviate, the
velocities will vary in direction and the effective flow section is
evidently less than the overflow height.
In this case the discharge coefficient Cd, still related to Equation 5, is
an empirical one and no relation is derived from an estimated contraction
coefficient.
13
At the broad-crested wejr or sill (symbols, see Figure 4b) streamlines are
straight and thus there is hydrostatic pressure. Application of the
Bernoulli equation for the broad-crested sill results in:

q = h 1 '12g (H 0 - h 1 ) (6a)
q gets maximal when h 1 = 2/3 H 0 , and for modular flow one gets the equation
presented in Figure 3d:
q 0 = (2/3) H 0 'I(2gH 0 /3) (6b)
(H 0 is again the upstream energy-head level above the sill level and q 0 the
discharge of modular flow per unit width). The critical waterdepth d above
the sill (being (2/3) H 0 ) is the one where application of the Bernoulli
equation (see also Figure 4b) leads to the maximum discharge, which remains
when the downstream level is further lowered.
For flow under a gate with a contraction coefficient C see Figure 3e) the
Bernoulli equation, also based on a hydrostatic pressure distribution,
results in a velocity proportional to 12g (H upstr _ hdtr) giving:
q = C c aI2g(H -C c a) (7)
0
(here a is the height of the opening under the gate)

14
mo,ph.rk pre,ure
g =Cc wf (H 0 1)dz
ti c
t-1 i 2' '1
3/2 3/2
Cc w \/2gI 1 (He)
Pa
ATMOSPIIERIC PRESSURE OVER 9ECTION 1-2 (VERTIcAL FARTTkL 9LIT)
etmopherle presIre
11
1fCCW
----t
ti c -\ Ii; - view
Q C,w2/3/2g 140
3/2
2 \
--
ATMOSPI1ERIC PRESSURE OVER SECTION 1-2 (vertieel lIt)
P.
nerIy ettn.pre5ure and more or les
horizontel velocity In AA
/-
-Fv 2a
3/2
q Cd
D (imIter to cee b)
,
NEARLY ATMOSPIIERIC PRESSURE IN CONTROL SECTION AA (2 dim.)
i/3H
= 2/311 e q = 2/3 1-1 0 \/2g (1I 0 /3Y
(AA control ecLlon)
IIYDROSTATIC PRESSURE IN CONTROL SECT1ON AA (2 dma.)
I - ntmospherle
. AC q c \/2g(tl0 ce )I
L,,. A
ITYI)ROSTATIC PflESSURE IN AA (2 DIM)
Figure 3 Modular flow conditions with atmospheric pressure and with

hydrostatic pressure in the control section
15
A variant for the discharge relation under the gate can be derived directly
from Eq. 7 when instead of the upstream energy level H the upstream water
level h 0 is introduced as reference for the upstream condition.
q = C a 12gh 0 /11 + Ca/h 0) (8)

c
Remark;
This can be proved as foliows:
Equation 7 can be written as q 2 = C 2 a 2 2g (h + (q 2 /h 2 2g) - Ca)
or: q 2 (1 - C2a2/h) = C 2 a 2 2g h 0 (1 - Ca/h) and after division of
both terms by (1 - Ca/h) and taking the square root Equation 8 is
obta ined.
It is important to establish the maximum water level downstream of the

control section for which no reduction of flow occurs. This downstream
water level is different for the case that at modular flow the control
section has a hydrostatic pressure distribution or that there is
atmospheric pressure over the whole section:
- at atmospheric pressure any pressure higher than atmospheric pressure
reduces the discharge, therefore any water level above the crest or the
bottom has an effect.
- at hydrostatic pressure a downstream water level lower than the one at
the control section, has no effect on the discharge, so in the case of
Eq. 6 a downstream water level just downstream of the control section
lower than 2/3 H 0 , gives no discharge reduction.
16
1.6 Fully submerged flow
II
The fully submerged flow regime as i i
t s ndicated
i n i Figure 1 is related
to the limit case of a (nearly) horizontal water level and very low
velocities (low Froude number). Because the flow boundaries are now fully
defined, the data of internal flow losses can be used when these are
available for similar geometries. The plane water level can be seen as a
flow boundary or as a line of symmetry of a mirrored situation, so still
more comparable situations might be found.
The losses at internal flow are expressed in terms of

AH = V2/2g (9)
and this also applies for the fully submerged flow. The influence of AH on
the losses E is negligible now. The total loss E is built up by a series of
local losses (inlet, gate, outlet, friction). The final result can now be
transformed as follows:
AH = AH i + AH + AH + AH etc. (9a)
tot nlet gate friction outlet
v2
AH = etc.) (9b)
friction outlet
tot 2g inlet + gate + +
v=_ 05
I2gAH (9c)
Q = Cd A I2gAH (9d)
At overflow the flow section A varies with the water level. Due to a
variation in flow geometry E and Cd do not remain fully constant either.
The friction losses are calculated in a similar way as found for internal
flow. The Moody diagram, as is commonly used for internal flow in circular
pipes, can be used but the "hydraulic diameter" Dh in this diagram has to
be replaced by the hydraulic radius Rh which is defined as the sectional
area divided by the wetted perimeter.
For a pipe with circular section this gives:
Rh = D/4 (10)
17
A problem with the calculation of the friction is that the structure is not
always long enough for getting the fully developed flow distribution, hence
friction can be smaller or larger.
At weir flow, starting with fully submerged flow and gradually loweririg of
the downstream water level, the intermediate flow condition starts already
at relatively small head differences. This is observed especially at weir
flow (Figure 4b) where the flow section, being the water depth h 1 over the
crest, varies with the downstream head h 2 . As has been indicated before,
this condition is mentioned "submerged flow" (in contrast to modular flow
and also in contrast to fully submerged flow where the water level is
approximately horizontal).
An example of the loss coefficient of a sill in the fully submerged

condition is shown in Figure 4a. The water level being nearly horizontal
and the downstream face being vertical, the pressure distribution is a
hydrostatic one and thus the momentum equation can be applied. But because
of the water level being nearly horizontal also the Carnot equation can be
applied and this lead to an analytical expresson for the discharge q as
function of h0 and AH. This can also be done for few under a gate,
Figure 4c.
The calculation for the losses is for the submerged condition, although
here the flow section varies with the downstream water level, quite
similarly as for the fully submerged condition.
18
Bernoulli eq Momentum or Carnot eq. Hloss of energy head
h2 h 0 and 1-ih0
- - 12
- 2gAH q/h 1 _q/(h 2 p) (Carnot)
h2
q = h0 (h 0 +p)/ p
F'ULLY-SUBMERGED FLOW REGIME (small H)
"Bernoulli" "Momentum" only

- H
Bernoulli equation
2
H0 h0 2g
h2
+ (l)
Momentum equation
INTERMEDIATE FLOW REGIME
pg ( * p) 2 pq 2 h 1 = g(h 2 p) 2 +p 2 q2 "h 2 p)
1)
h 1 eh 0 and h2Ch0 2
h0 bi1 ,-- h2 2gAH= q/C_q/h 2 (Carnot)
q=ce /(l_c 0 'h 0 )
FULLY-SUDMERGED FLOW REGIME
Bernoulli equation
2 22
h0 + 0 q 2/ h 2g =h 1 -e 1q /c 2g
Th
o h h2
(o 1)
Momentum equaUon -
INTERMEDIATE FLOW REGIME
f pg h 12 1 pq2/c = ---- pg h+ 2 pq2/h2
Figure 4 Long-crested weir and gate flow under fully submerged condition
and in the intermediate flow regime
19
1.7 The intermediate-flow relations at control sections
In the intermediate-flow conditions direct application of the Bernoulli

equation (from h 0 to h 1 ) and the momentum equation (from h 1 to h 2 ) leads to
a feasible result for the flow over a long sill (Figure 4b) and for flow
under a gate (Figure 4d). But a discharge relation like Equation 1 cannot
be found directly. Starting with h0 , q and h 1 , the h 2 can be found through
an iterative search or by solving a third-order equation. With the
intermediate flow regime there are certain cases where mixed-flow can be
assumed. A part of the discharge is related to free flow and a part to
submerged flow. These cases are illustrated in Figure 5. When using these
mixed flow equations with constant Cc coefficients resuits may not be very
accurate but certainly useful for a first estimate of the discharge. The
procedure is fairly evident for the case of a half submerged narrow
vertical slit, where the horizontal flow contraction will not vary strongly
and where one can clearly distinguish a free-flow part and a submerged
part. At a nappe or at flow through a horizontal slit the pressure in the
control section is also nearly atmospheric but the (vertical) flow
contraction is more variable.
A usual way of expressing the weir flow discharge in the intermediate flow
regime is to start with the modular flow discharge Q 0 and then apply a
reduction factor C 5 which is a function of the dimensionless submergence
factor S.
S = h 2 /h 0 or h 2 /H 0 (11)
QCsQo (12)
where Q 0 is the discharge at modular flow.
Q0 = f(H 0 ) (13)
By introducing a number of simplifications the reduction factor C due to

submergence can be expressed analytically for both a short-crested and a
long-crested weir.
20
top view
r 3/2 3/21
q b Ib Q = Q. b
0
al situation
UT2 :ebon
3/2- 3/21
s a) s
a
q= b\J2g(sa)
:7
q = q b
0 2 -dimensional situatlon
0_b) 2
a * ( H
q b = b'\Jg(H Q_b)
q = q+ qb
Figure 5 Vertical slit and horizontal nappe with rnixed fiow conditions
21
A: Short-crested weir (Figure 5c):
q 0 is found from the equation shown in the figure when b = 0 is introduced:

q0 = Cd(2/3) (V'2g) H0 3/2 (14)
Then we introduce for the fully submerged condition b = h 2 (this is only

true when the sili is relatively high) and we assume for the discharge
coefficients Cd = Cd = Cd
a b
The height b in Figure 5c can be expressed as H 0 S, and the discharge

relation of Figure 5c becomes:
q = + O= Cd(2/3) (hj) (1-S) 3 ! 2 H03/2 + CdS H 0 12g H 0 (1-S)
on
C = q/q 0 = (1-S) 3 ! 2 + ( 3/2) sf

Cs = (1 + S) 1(1-S) (15)
This analytical expression was compared by Abou Seida and Quarashi [2] with
experiments, see Figure 6.
1.0
broadcrested weir, eq(19)
Ii-
00
08
0.6
0. 4
0.2
Z4-
- experimental points Abou Seida;Quaraishi
o.x eo-.3o
e =45.60*
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Figure 6 Discharge reduction by submergence

22
B: Broad-crested weir
Equation 6 was:
q 0 = (2/3) H 0 (12g H 0 /3 (16)
When the weir height p is large compared to the overflow depth, then in the
situation of Figure 4b the water level above the crest gets equal to the
downstream water level. "Bernoulli" then resuits in:
q = h 2 12g (H 0 - h 2 )
and, introducing S = h2/H M h 2 /H 0 one obtains:
q = S H 0 /2g H 0 (1-S) (17)
(17) hoids for S > 2/3, otherwise q = q 0 (18)
This results in:
= 1 for S < 2/3 and
= (3/2) Sv'3 (1-S) for 2/3 < S < 1 (19)
Figure 6 shows the verification of Equation 15 with experiments on the

sharp-crested weir in vertical and in 30 inclined position, and with a
triangular sloping filling block at the upstream face on 450 and 600 (from
Abou-Seida and Quaraishi [2]).
Again, the above relations for the control section must be extended to the
water level further upstream and downstream for obtaining the discharge
relation for the whole structure. The extra losses upstream will not cause
problems but for calculation of the water level further downstream where
the flow diverges, one could try to apply the momentum equation. However,
when the downstream canal has diverging walls or when there is a sloping
bottom, not all the forces which need be introduced are known; this in
contrast to the sill with vertical back wall. What might be feasible, but
this is not sufficiently verified, is that the rise in water level (by
recovery of potential energy) with a certain slope (estimated to be some-
where in between 1:5 and 1:7) is used to introduce an extra horizontal
23
force in the momentum equation. Assuming a hydrostatic pressure distribu-

tion, the longitudinal force component acting at the diverging wails or
sloping bottom can then be calculated; this force indeed gives a reduction
in energy head losses at diverging walls. This procedure has originally
been suggested by Kooman [ 2 1].
It can be seen that the calculations in the intermediate flow regime do not
contain new elements, but computations should be done carefully. The steps
or slopes in the water level do change the flow geometry so that the losses
cannot be found so easily from literature data. The proposed computer
program of Chapter 2 is suitable to handle the trial-and-error procedure
which is related to the water-level computations in the intermediate flow
regime.
11
1.8 Remarks about entrance losses
At the entrance of a sluice or a culvert the pressure lowers due to the

velocity head V2/2g (wherin V is the culvert velocity). However, there is
the "normal" flow convergence which depends only on the culvert section,
and there is an extra contraction which depends on the mouth shape. The
extra contraction will disappear further downstream but the result is an
extra loss (similar to Carnot losses at a suddenly widened flow section).
Per definition:
H = H + E (V 2 /2g) (20)
0 c e c
or
h0 + V/2g = h + (1 + E ) V2/2 g (21)

c e c
h0 and H 0 are pressure head, respectively energy head upstream from the
culvert, and h and H are related to the ones in the culvert V = culvert
c c c
velocity.
E e is determined from pressure and discharge measurements.

= {(h 0 + V/2g - h)/(V 2 /2g)} - 1 (22)
In general it is assumed that the velocities are evenly distributed over

the culvert section.
However, it has been discussed in section 1.3 that in reality an a factor
should be introduced to take into account the uneven velocity distribution
(Figure 2). Because inlet losses of nicely shaped mouthes are small (in the
order of 0.15 to 0.3 in prototype and still smaller in model, see section
4.1) a few percent greater (V2/2g) will also result in a in reality smaller
E e value (a few time 0.01 less).
In the total analysis of a structure the neglect of a is generally not
important, because at the outlet the extra kinetic energy which initially
is neglected is lost anyhow.
In par 4.1 the relation between the flow contraction coefficient C and the
entrance loss E is discussed.
25
In a closed culvert or at lower water velocities the condition is

comparable with a condition with a sudden flow expansion and hence "Carnot"
can be applied to the entrance losses.
Using the symbols of Figure 7,
1
AH = (V - V ) 2 /2g = (C_ - 1) 2 V 2 /2g
m c c
and hence
_1
= (C - 1) 2 (23)
e c
Vm = V ICc
entrance losses
Figure 7 Entrance conditions.
C varies from about 0.55 (at a an ali-sided sharp-edged entrance) tili

about 0.85 for an elliptically shaped entrance, which should result in
values from 0.03 to 0.67. The low values, however, are never reached in
reality (see Figure 65).
26
1.9 Remarks about trash-rack and bridge-pier losses
Data about trash-racks found in literature refer to one-dimensional flow

without influence of the free water surface. But these data can also be
applied at lower water velocities with a nearly horizontal water level at
free surface flow. In other conditions these data give losses which are too
small. Even modular flow can occur jn the flow contraction.
This is the reason why bridge-pier data are also related to the Froude
number v//j, where V is related to velocity and d to water depth of the
undisturbed condition.
There is a direct relation between the losses of trash-racks and the

hydrodynamic forces exerted on the bars. This relation foliows from the
application of the momentum equation. When the rack is fully submerged then
the inflow momentum equals the outflow momentum and the momentum equation
results in:
p g AH = rack resistance force per unit area of the culvert (24)
(AH = loss of energy head, p = fluid density)
and when free surface flow is involved one can derive from the momentum
equation per unit width:
Vi pgh 1 2 + pqV 1 = Force + pgh 2 + pqV 2 (25)

27
There are three types of approach to establish the trash-rack losses (see
Figure 8):
w/b = w/W
trashrack loss coefficient
C = flow contr.oUon coefficieit
C F= force coefficient
A Wi 4
vo
)2
A i =(v/zg) (Carnot)
(Wtw) 1
B force (per meter 2 )= C(rw/W)4- pV 2 = pgiH
wlth C f (barhape. and Ew/W)

F =
C: H=, t v 2 / 2 g with= t(bar ihape.and EW/W)
Figure 8 Flow conditions at trash-rack and different types of approach for

calculation of losses
Estimation of the flow contraction and establishment of the losses by

applying the momentum equation. The assumption is that in the whole
flow-contraction section an equal pressure exists. (For small blockage
ratios this does not need to be true!) Without any influence of the
free surface (in a culvert or at a horizontal free water surface) the
Carnot equation can be applied.
Estimating the drag coefficient of the bars and transferring the total
force into a hydraulic loss.
Using direct data of trashrack tests where losses are presented as
function of shape and density of the bars.
Mostly the third method is used, but it is advisable to apply also method A
so as to regard the consequences of extra blockage by trash and for
checking the possible effect of the modular or intermediate flow regimes.
1.10 Reauired and obtainable accuracies
In ISO 1438, included in the ISO handbook [18] are discussed the accuracies
of discharge relations in prototype; the accuracy of hydraulic models are
discussed by Kolkman [20]. One of the conciusions is that accuracy
considerations must be set up separately for each specific case.
The possibility of systematic errors and scatter in model results must be
considered both. The definition of error can be related to the actual
discharge or to the maximum discharge.
Even when a scale-model investigation is considered, in the pre-design

stage it is important that a good estimate can already be made about the
discharge capacity of the structure. For this case it is important to know
how much one can rely upon data from literature.
The necessary accuracy is at the pre-design stage, say 5 to 15%; this
percentage is related to other inaccuracies of the design input (water
level, the once-in-a-hundred-years river discharge, etc.), and to
construction costs when for instance too-large gated openings are designed.
Especially in the United States a large number of investigations have been

carried out to establish systematic design data of free (modular) flow
construct ions.
For these modular flow conditions it is certainly possible to obtain an
accuracy in discharge between 3 and 8% of the actual discharge. For
discharge sluices in tidal areas the intermediate flow regime occurs
frequently and here the accuracy of estimated losses is less (up to 20%
errors is the experience).
There are structures with a diffusor-shaped outlet resulting in a high
discharge capacity, but systematic data for the functioning of such
diffusors at free surface flow are missing.
When using the structure for water-discharge control the required accuracy
will be within 3 to 10% of the actual discharge. 1f a greater accuracy is
required than 10 to 15%, a specific scale-model study should be set up
and/or a seperate discharge meter has to be instailed.
In certain cases the structure can be calibrated with the use of propeller
velocity meters, acoustic flow meters or with tracer methods. In these
circumstances the obtainable accuracy is variable and depends largely on
the measuring section, the accessibility, the number of measuring points,
the flow distribution and flow direction, turbulence and so on. When there
MA
is a straight canal or culvert in the structure an accuracy of 2% to 5% can

be reached. Again, a discussion about accuracy of measurements is found in
the ISO handbook [18].
When a high accuracy is required in discharge metering with the structure,
the design must be modified to that purpose. A scale-model investigation
must include an extended up- and downstream basin, the water levels must be
measured at points corresponding to the location in prototype, the scale
must be adequately chosen in relation to the wanted accuracy and a
correction procedure should be used to compensate for scale effects due to
friction. Structures with trash-racks need special research procedures. In
general a structure with a trash-rack is not weli-suited for discharge
metering.
In the presentation of discharge coefficients following here, the influence

of the Reynolds' number (viscosity effects) is not inciuded. For the
prototype structures this is of minor importance, and hence in the pre-
design of the structure it does not play any role. For the choice of the
scale factor in hydraulic modelling and for the interpretation of scale-
model results specialized knowledge is needed; this is beyond the scope of
this manual.
30
2. NUMERICAL PROCEDURE FOR COMPUTATION OF WATER LEVEL/DISCHARGE RELATIONS
2.1 Introduction
ne Wo! pierW/[
entrance gate exit
1 = 1 a it - 3 4
III B
A Ll
Figure 9 Water levels in a discharge sluice
In Figure 9 a sluice with free-surface flow is schematized throughout as a

chain of discontinuities for each of which an equation can be set up
relating water levels and discharge. In this schematisation no interactions
between the sections are introduced.
A is the inlet of the sluice. Sluice width is introduced as w. The inlet
loss coefficient is E .. The relation between h 0 , h 1 and Q can be written
as:

h0 + a0 (Q/w0 h0 2 /2g = h 1 + a 1 (1 + E ( Q/w h 1 ) 2 /2g (26)
)
l )
(w 0 is the width upstream of the sluice, w is the sluice width and is

the entrance loss coefficient).
When the flow is supposed to be evenly distributed over the flow sections
the gate-discharge relation becomes:
h 1 + (Q/h 1 w) 2 /2g = h 2 + (Q/Caw) 2 /2g (27)
Here C is the flow contraction coefficient in opening a.

31
When the computation starts from upstream, there is a certain Q whjch

cannot pass at the upstream water level which was introduced; one can now
for instance by trial and error find the modular flow discharge (h 2 cannot
be smaller then the contracted jet Ca).
The gate loss is produced in the deceleration zone of the flow further
downstream of the gate. This loss can be computed using the momentum
equation, provided that Al (in Figure 9) is sufficiently large:
pgw (0.5h 2 2 ) + pQ 2 /(Caw) = pgw(O.5h 3 2 ) + pQ 2 /(h 3 w) (28)
The exit loss at B (a sudden widening) can also be determined with the
momentum equation, assuming that the water level in the basin near the
sluice is related to the water level h 3 , as follows:
pgw (0.5h 3 2 ) + pQ 2 /(h 3 w) = pgw (0.5h 4 2 ) + pQ 2 /(h 4w) (29)
At a gradual widening an empirical relation between h3, h4 and Q has to be

int roduced.
When a hydraulic jump downstream of the gate is involved one obtains:
pgw [0.5(Ca) 2 1 + pQ 2 /wC a = pgw (0.5h 3 2 ) + pQ 2 /(wh 3 ) (30)
When refined calculations are needed, one can still introduce the a and
factors of Figure 1.3, the friction of the culvert walls, a varying C
value of the gate opening, depending on gate position, and on the water
level h 1 . All these equations are coupled. It is common to solve each
equation separately, and the "next" water level (more up- or more
downstream) is found, mostly by trial and error, from the calculated energy
level.
IJsing a computer one can now solve these equations in a standardized way.
For each of the formulas we index the water level of the left-hand section
as h 1 and at the right hand as h 2 . In the case of a hydraulic jump
computation downstream of a gate, we define Ca as being h 1 for instance.
Each of the formula which were discussed before can be wrjtten as follows:
f(h 1 , h 2 , Q) = 0 (31)
32
Starting the calculations from upstream with a given h 0 an Q, the given

water level h 0 serves as the h 1 of Eq. 31. When also the Q is introduced
then Equation 31 can be written as

f(h 2 ) = 0 (32)
For the search for this h 2 standard procedures are available. In section
2.2 a procedure will be presented, written for the ACORN-BBC or the
Archimedes computer.
When Eq. 32 is solved for h 2 , this is the h 1 to be introduced in the second

formula, and then the next h 2 can be found.
The iterative procedure for the determination of the modular flow

conditions is relatively easy as well; when the discharge is too high in
combination with the chosen upper-water level, then no value of h 2 is
found. Hence the modular discharge is the maximum discharge which stili can
pass.
But a direct way of computation is to set up the equation for modular flow
at the control section and do the same type of calculation in the upstream
direction. Downstream of the control section one proceeds in the normal
way.
In the search procedure for h 2 one has to define the h 2 range where the
search takes place. One limit of this range can be h 2min near zero (not
exactly zero because then certain terms of the equations tend to infinity).
can just be an arbitrarily large value of h 2 . Not an extremely
The 2max
large value, because then the search steps generally get large as well, and
two solutions (critical flow and supercritical flow) can be surpassed in
one step.
When one wants to find the subcritical flow, then the search for h2 has to
start at h and one goes downward; for the supercritical flow it is the
2max
other way around, that is working upward from h2.
33
2.2 The search procedure
Equation (32) is transformed into
f(h 2 ) = Z (33)
in which Z has to become about zero, with a wanted accuracy of c.

A simple but fool-proof procedure would be stepping through the h 2 range
with steps Ah, thus obtaining Z 1 , Z 2 and so on.
Whe n
(z )z
~O (34)
n-1 n
the search is stopped, and the procedure starts again with h 1 and h as
new limits of the search range. The Ah steps must be small enough for the
first search, because near the critical waterdepth the two solutions (sub-
and supercritial are near to each other). The simplest procedure is to stop
the iteration when:
f((h1 + h)/2) < c (35)
The wanted h 2 value is then
h = + h )/2 (36)
2 (h n-1 n
In the appendix to this chapter a much faster procedure is used, based upon
the Newton-Raphson method . At an arbitrary begin-value of h 2 the
derivative is computed

dZ/dh 2 = (37)
(Ah 2 has been chosen as 1/2000 for the h 2 search interval)
When from this begin-value of h 2 , with its corresponding Z value the

tangent line is considered, the Z = 0 corresponds with a new h 2 value.
* At a suggestion from H.K.T. Kuijper, Delft tJniversity of Technology.

34
It is checked whether
Z = f(h 2 )l < wanted accuracy
When IZI is stili too large, h replaces the value of h and the
2-new 2
procedure is repeated.
This procedure is incorporated in a new-defined function, FNfast-implcit'

with the following variables:
M = the nuinber out of a series of self-defined functions of the type

f(h 1 , h 2 , Q) = 0 (like Eq. 31).
h 1 (defined as in Eq. 31).
Begin-h 2
these values determine the search range of h 2
End-h 2
Q (discharge, see Eq. 31)
This procedure is fully adequate to solve h 2 from the Bernoulli or the

momentum equation, but when using other types of equations, this procedure
can evoke problems.
35
APPEFILIX of CUAPTER 2 page 1
41 th the ACORN NRC computer the foliowing function "FNf act impi icit"
bas been developed uhich produces the waterlevel h2 uhen one introduces
the discharge 0, the waterlevel hi and the self defined relation
F(hl, 1i2, 0)= 0. The h2 value is found when the function is called for in
the following weyr
h2= FNfastimplicit(fl, hi, beginh2, endh2, 0)
The numher M is the numher of a series of .self-defined functions which is

referred to. The begin h2 and end_h2 values define the renge in which h2
bas to he found. With e loo value of heginh2 the supercritical solution
is found and with m high value the subcritical one.
As an example the sequent depth of a hydraulic jump is produced.

The momentum equation is as follows:
.5rho'g(h1'2) rho'(q2)/hl= .5'rhog'(h22)4 rho'(q2)/h2
This equation is transforrned, putting all terms in the left hand:
.5'rho'g'(hl2)4 rho'(q2)/h1- .5rho'g'(h22)- rhn'(q2)/h2= 1

and Z hem to he zero!
The input in the form of a program is as folloos:
REMEIIBER program starts here.

PRINT' "OIJTPRINT"
h11: q=15: g=9.81: rho=1000
begin h2=40: RFJI this seans that the subcritical solution is found.
end h2= .01: REN the small begin_h2 or end h2 should be greater than zero
11=1: REM this means that the first self-defined function is used.
h2= FNfest impiicit(i1, hi heginh2, endh2, q)
Pil 1 NT"hr'gi nh2= " begin Ii?; " end_1)2= " endh2 ; " h2= " h2
END
RF.MARM in the following function the momentum equation is represented.

LEFFIJse ii (lef i ned 1 (bi h2 , q
Z= .5'rho'g(h1'2) rh0'(q2)/h1- .5rho'g(h22)- rho'(92)/h2
OUTPR INT
begin h2= 40 end h2= IF.-2 h2= 6.29128557
Other begin and end values of h2 result in:
OUTPR INT
begin h2= 1E-2 end h2= 40 h2= 0.999999999
OUTPR INT
begin h2= 10 end h2= 40 h2= -99.99
Appendix of Chapter 2: Iterative search procedure for h 2 using the /\CORN-BBC

computer.
Appendix of Chapter 2: Iterative search procedure for h 2 using the ACORN-BBC

computer.
36
APPENDIX of CIIAPTER 2 page 2
REHARK the following contains the essential additional functions
REM h2 has to be found

LOCAL endh2, heginh2, wanted_precision, precision, steph2, Pi, heginZ, endZ
LOCAL tangent alpha, former_tangent_sipha
endh2= end_h2: beginh2= begin_h2
wanted_precision=ARS((end_h2-begin_h2)/5OOO)
Pl=0: outof_range=FALSE: stop'FALSE
REPEAT
5 teph2= (endh2-beginh2)/2000
endh2= beginh24steph2 -
h2=heginh2: heginZrFNadditional(Mi,h1 ,h2,O)
h2= endh2s endZ=FNadditional(M7.,hl,h2,O)
1F P =2 UIEN former_tangent_alpha= tangent_alpha
tangent alpha=(endZ-heginZ)/(endh2- heginh2)
1F tangent alpha <>0 THEN trial h2= beginh2- (beginZ/tangent_alpha)
IV P%=2 AND tangent alphaformer_tangent_alpha>=O THEN out_of_range=TRUE
1F (trialh2-beginh2)/(end_h2-hegin_h2)>0 THEN out_of_range=TRIJE
TV (tria1_h2begin_h2)/(end_h2-begin_h2)1 THEN out_of_range=TRUE
endh2=beginh2: heginh2= trial_h2
h2= trialh2: Z= FNadditional(M,h1,h2,O) precision=ABS(Z)
TV precision >= wanted precision OR out of range TFJEN stop= TRUE
UNTIL stop
1F out_of_range TIIEN resu1t=99.99 ELSE result =tral_h2
= resu1t t ttttt*tttt It*ttttttt* t Ititti tt*ttt*ttt,ittttttt tttt
DEF FNadditional(H7,h1,h2,Q) LOCAL Z

TV 11ir1 TIIEN Z= FUselfdefined1(h1,h2,(j)
IV Mi=2 TIIEN Z= FNselfdefined2(hl,h2,0)
TV F17.=3 UIEN Z= FNselfdefined3(hl,h2,O)
Z: RE1IARK indicates that Z is the output of the seifdefined function.
Appendix of Chapter 2; (continued).

37
3. DISCHARGE RELATIONS FOR CONTROL SECTIONS
Remark:
As noted before h is used for the water level and H for the energy head.
When only the upstream water level is involved then h and H refer to the
upstream conditions. Otherwise index 0 refers to the conditions upstream of
the control section, 1 to the conditions in or just downstream, 2 to the
condition further downstream. In U.S. reference like [38], [40] and [ 4 3],
H (the real or the effective energy head) is used in contrast to the
design head Hd which is a reference size whereupon the shape of a nappe-
shaped crest is based.
3.1 The sharp-crested weir
Discharge relations for modular flow and the submerged flow are discussed
separately.
A. Modular flow
Still in use and best-known is the weir-formula developed, and extensively

tested between 1925 and 1930 by Rehbock [30], at the of the Karlsruhe
University. The sharp-crested full-width weir has since then also been
widely applied as a measuring weir in laboratories. In Figure 10 a section
of the weir is presented. When installed precisely, obtainable accuracies
are within 1% [18]. The weir is sensitive to roughness of the upstrearn face
and at an imprecise installation 5 to 7 inaccuracies can occur (according
to Bos [7]). Therefore its application as a discharge measurement structure
in the field is less suitable.
The formula is based on a two-dimensional flow pattern and presents as:
Q = w (1.78 + 0.24 * h/p) (h + 1.1 10)3/2 (39)
The equation is not dimensionless, and the dimensions of the parameters

are: discharge Q is in me/s, the width w, the weir height p, and h are all
in metres; h is the upstream water level above the crest.
The flume width upstream of the weir equals the weir width w.
The correction of 1.1 mm in h caused some problems when it was tried to
transform Equation 39 into a dimensionless form. It certainly contains a
mixture of viscosity and surface-tension effects.
38
For a measuring weir the crest shape has, according to ISO [18], to be
designed as indicated in Figure 10.
Special care is required for the upstream face to be vertical and smooth,
and to get a complete aeration under the nappe.
Sarginson [34] proposed to replace the value 1.1 mm by a factor

representing the influence of the surface tension (o) by the Weber number.
He defined this number We as:

We = pgh 2 /0 (40)
After analyzing theoretically the possible effect of surface tension, he

proposes the following discharge formula:

Q = w(l.81 +.22 h/p +4.22/We) h 1 5 (41)
Also this expression is not dimensionless.

The expression by Sarginson deviates up to 1.5 % from the Rehbock values.
Accepted by the ISO committee on water-flow measurements [18], are nowadays

the experimental results of Kinsvater and Carter [19]. They investigated
the discharge over a sharp-crested weir at different ratios between crest
width (w) and flume width (W). They proposed the following discharge
formula:
Q = C e (2/3) ) w ehe 3/2

(I (43)
This formula resulted from investigations for a wide range of W/w and h/p
ratjos. Ce is presented in Figure 11. we is defined as the effective breath
which equals (w + Kb); he is defined as the effective water-level height
being (h + Kh). The h correction Kh remains 1 mm for all situations, the Kb
can be read in Figure 12. Again, the authors attribute the small Kh and Kb
corrections to the viscosity and surface tension effects. The 1 mm
correction means that the total weir formula is not completely correct in
dimensions (according to the theory of dimensjonal analysis).
The differences between the resuits of Kinsvater et al. and Rehbock, are
illustrated with the following results for w/W = 1. The conciusion is that,
applied in practice, they both produce the same result.
39
h p b Q Rehb. % difference
QKinsv
(m) (m) (m) (ma/s) (ma/s)
.05 .4 .5 .01045 .0104 + 0.5

.1 .4 .5 .0296 .0296 0
.2 .4 .5 .0857 .0855 + 0.2
.2 .6 .8 .1342 .1339 + 0.2
A correction factor for the discharge, in case of insufficient aeration, is

mentioned by Bos [7]. It relates to air pressure underneath of the nappe.
AQ/Q = -0,2 (p 2 /pgh). 0

. 92 (43)
The author compared this theoretically with the resuits of Sarginson in

relation to surface tension effect (which also leads to a pressure in the
water which deviates from the atmospheric one); however, the latter leads
to a 5 times greater correction than Eq. 43. It seems reasonable to apply
Eq. 43, which is based on experiments.
B. Submer2ed flow
The ratio of submergence is defined as:

S = (44)
Now h 0 is the definition for the height of the upstream water level above
the crest, h 2 is the downstream water level, also above the crest.
Several authors have presented a reduction coefficient for the discharge,
where S is the only relevant parameter and where h 0 /p is left out.
The discharge is always compared to Q 0 , the modular discharge.
Villemonte [41] presents:
Q/Q0 = (1 - S 1 . 5 ). 305 (45)
Varshney and Mohanty [39] propose (for .03 < S < 1.0):
Q/Q0 = f103 - 0.27S - S 2 ) - . 059 (46)

3) As it is presented in Chapter 1, Figure 6 and equation 15, Abou-Seida
and Quarashi [2] propose an analytical approach, resulting in

Q/Q = (1 + s) 1(1-S) (47)
It is remarkable that Eq. 46 does not reach Q = Q 0 when submergence is

zero, but this in on purpose: Varshney and Mohanty present experiments
where a kind of jump in the curve is observed between S = 0 and S = .03.
Theoretically it might be imaginable that at zero submergence the discharge
increases, because the last bit of air underneath is sucked out by the jet
and a little underpressure occurs.
All three equations are presented in Figure 13. Because at the modular flow
regime atmospheric pressure occurs under the nappe, any extra water
pressure will reduce the discharge, so a downstream water level higher than
the crest has a direct effect.
In the discussions presented by Villemonte [41] it is interesting to see

that he also points to to Eq. 47, while referring to the Frenchmen Dubuat
(Principles d'Hydraulique, Vol. 1 p.203, 1816). Villemonte shows that his
own experimental results are better (within + and - 5% of his formula). At
his smallest investigated h/p value (p being the weir height) of 0.06, his
experiments yield a higher discharge, especially SO when the submergence is
under 0.4. His higher h/p range (up to 0.25) results in a lower discharge
compared to Equation 45, especially when S > 0.4.
For the moment one must conclude that the bundle of curves in Figure 13
shows the possible scatter in discharge relations of submerged conditions.
crest width 1 - 2 mm
Figure 10 Section of the Rehbock standard weir.

41
0.80
076
0.72
Ce :0,602+0075 h/p
068
0.64
04
0.60
0,56 -
0
T-770
0.4 0.8 1.2 1.6 2.0 24
o voJues of h/p
Figure 11 Ce factor as function of h/p and w/W, also after ISO 1438 [18]
0.005
N
(Y)
Ln
0
0.004
N
0.003 0 0 0
$ 0.002
1 : 001
0.001
expected tolerance on Kb: 0.0003m
-0.002
0 0.1 0.2 03 0.4 0.5 0.6 0.7 08 0.9 1.0
ratio w/W
Figure 12 Values of Kb as function of w/W, from Kinsvater and Carter,

after ISO [18]
42
Ce
after fig 6
(eq.47)
WA
ViiIemonte
Varshney/Mohanty/
\ \
\ :
0.4
0.2
l
0,2 0.4 06 0.8 S

ire
> ratio of submergence
Figure 13 Flow reduction by subrnergence at a sharp-crested weir; different

express ions.
43
3.2 The weir with circular crest section
In contrast to the sharp-crested weir, the semi-circular weir does not have
a fixed point of flow separation. From a certain overfall height onward
(H/R = 1.5 to 1.8) the pressure under the nappe becomes lower than
atmospheric, and from this point on when H/R is increased, aeration can
occur. When aerated the flow separation point shifts backwards, and comes
below the crest level. Just upstream of the separation point the pressure
is still lower than atmospheric (flow separation is retarded, particularly
at the smaller diameters), and related to this retardation and to the low
pressures at the crest, the nappe becomes more curved and the velocity
becomes higher. For these reasons the discharge coefficient is greater than
for the sharp-crested weir. In the region 1.5 <H/R< 2 it is still difficult
to obtain a good aeration in a natural way.
Thus, as a measuring device the semi-circular weir is only used from H/R>
1.5 ori.
Rouv/Indlekofer [32] performed tests on a weir with semi-circular crest
(Figure 14) where, with some of the tests the air is sucked out, so as to
obtain comparative tests with and without atmospheric pressure underneath.
They also refer to a number of published data.
Figure 15 and 16 show the discharge coefficient Cd in the following

discharge relation
Q = Cw(2/3) I H/ 2 (48)
The energy level H is introduced as the relevant parameter. Although a wide

range of the radius (1-15cm) and the weir height (20-95cm) is investigated,
the authors state that only the parameter H/R appeared to be relevant. The
weir height had only an effect within 2%. The test values in Figure 15 are
the average ones, obtained with a radius between 10 and 140 min.
For small H/R values, where the nappe remains attached, the discharge
coefficient is presented in Figure 16 in more detail. At a small radius Cd
decreases and this must attributed to viscosity effects.
Similar results are found by Sarginson [34] for the full-circle cylindrical
crest. (In his definition of the discharge coefficient Sarginson uses the
upper water level h and not the energy level F1). From Figure 17 it can be
observed that for diameters smaller than 86mm the discharge coefficient
44
decreases with the crest radius. This must be a viscosity effect related to
the boundary layer which probably means that also the roughness of the
cylinder surface can have a considerable influence. From h/R=2 (or greater)
the flow separates from the highest crest point. The Cd does not grow and
even tends to become lower for h/R, exceeding a value of about 4. The
protrusion of the cylinder upstream from the weir plate, is causing the
weir height p not having so much influence on Cd, because independent of
the weir height the water is almost stagnant underneath.
Figure 14 Weir section with semi-circular crest and with circular crest, a
as investigated by Rouv/Indlekofer [32] and b by Sarginson [34]
4.0 -
3.5 --
/
3.0 -1-'
1 II
2.5
L
naturally 1 /
H/R J
uttacfled
t
::
naturolly
ottachei
1.0
aerated
0.5
RV 2
0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90
- Cd
Figure 15 Discharge coefficients of semi-circular crest with and without

sucking-out air (from Rouv and Indlekofer [32] and cited authors)
YW
1.
R: 549cm smaH ( H/p)
R = 7.62cm
_._R:10.Olcm
-- R: 14.80cm
after Matthew -'-
0.5 for R:2.54cm

after K napp for
eat Radius
:7
- -,
Q Cd w2,3 H3/2
0.50 0.55 0.60 0.65 0.70

10 Cd
Figure 16 Discharge coefficient with attached nappe at varied crest radius,

from Rouve and Indlekofer [32].
Q:Cd W2/3\[ h 312

0.80
fr CA
v,

0.75
D 9 11
xV00
8 47 _____ _____
070
1/2 R
0.65 rl
0 2R mm symbol
c
75 6 8
0
0
13 o
-.
,x
- 0
/ 25 0
v X
86
SV
114 v
136 u
0.55 0 11 1 ______1 _____
0.1 0.2 0.3 0.5 07 1 2 3 S 10
hfR
Figure 17 Discharge coefficient for the circular cylinder crest, after

Sarginson [34]
3.3 The nappe-shaped overfiow weir
A. Modular flow
The crest shape of the nappe-shaped overflow weir was developed by experiments
in scale models where the location of the lower boundary of an overflow nappe
over a sharp-crested weir was exactly measured. All the radii of the curvature
are related to a reference length for which is chosen the "design head" Hd (in
fact the design energy head). When the energy head of free overflow H 0 equals
the design energy head Hd, then one would expect that on the whole crest
pressure is exactly atmospheric. In reality small deviations occur. The weir
height p is not introduced directly, but in the form of the influence of the
velocity head ha of the approach velocity in relation to Hd. The crest shape
is found in Figures 18 and 19. All data are taken from USBR "Design of small
dams", Ref. [38].
The nappe-shaped weir type has the property that when the real overflow height
(H or H ) equals Hd, the pressure at the crest and the downstream face equals
0
the atmospheric pressure, and hence only small risk of unwanted aeration or
cavitation occurs. When the real H value is smaller than Hd, then pressure is
higher than atmospheric, and when H is greater than Hd the pressure is lower
than atmospheric; the natural jet of the sharp-crested weir is thrown out
farther when H increases. The advantage of this weir type is its high
discharge coefficient. For instance, cornpared to the sharp-crested weir at the
same discharge, the crest is now at the top level of the underside of the
sharp-crested overflow nappe. This means that the equivalent crest height of
the sharp crest lays 12.7 of Hd lower than the nappe-shaped profile. An
additional safety exists for conditions at an extremely high flood, because
then the discharge coefficient becomes still higher as a result of lower crest
pressures (see also middle figure of Figure 20). This advantage, however, is
paid for in the form of an additional risk of cavitation.
The very systematic investigations of this weir type (crest shape and
discharge relations) have been performed by the US Bureau of Reclamation USBR
and the Vicksburg Waterways Experiment Station of the US Army (WES) [43].
47
The upper graph of Figure 20 is a revised version from [33], applying the
definition of Cd, similar as for the sharp-crested weir:
Q = CdW (2/3) (Ij) H 3 / 2 (49)
A few excercises and checks:

It is observed in Figure 20, the upper graph, that with a small H the
discharge coefficient reduces to 0.575. This case is comparable with the long-
crested weir for which was derived the theoretical discharge relation of Eq.
6. This equation equals Eq. 49 when Cd = l/v" = 0.577 is applied, which is
well in agreement with the value of Figure 20.
Also one finds for a great weir height a discharge coefficient of 0.739. When
this is combined with the reduction factor of 0.78 at small HO/Hd (again a
situation similar to the long-crested weir) then 0.739 times 0.78 is read in
Figure 20 and this becomes 0.575 which is also in agreement with the
aforementioned Cd of 0.577.
Also one can see that for H0 = Hd the discharge coefficient for a high weir is
0.739. Reading in Figure 19 that Y/Hd = 0.127 one can derive that the
equivalent sharp-crested weir is 1.127 Hd below the water level. The discharge
coefficient of the sharp-crested weir will be 0.739/(1.127 3 / 2 ) = 0.618 while
from the Rehbock formula (eq. 39) combined with Eq. 49 one can derive C 0 =
l.78/[(2/3Ij)1 = 0.603. The difference of 24% might be explained by the
correction in h of lmm which is proposed by Rehbock and which for small scale
tests results in a little higher value of Cd.
B. Submerged flow
Figure 21 shows the different types of flow which might occur at submerged
flow and the decrease of discharge coefficient. This graph, to be found for
instance in the USBR [38] was modified by Naudascher [24], who introduced P/He
on the horizontal axis. He is the upstream energy head, which generally will
be different from the design head (see in Figure 20 the central figure for the
influence of H/Hd).

48
C. Modular flow with piers
When the weir has modular flow conditions also the effect of piers and
abutments is found in the 't Hydraulic Design Criteria" of WES [ 4 3].
The width w which has to be introduced in Eq. 49:
w = [L -2(NK + K ) H1 (51)
p a e
L = total crest length between abutments and piers

N = number of piers
K p = pier contraction coefficient
K = abutment contraction coefficient
H: = upstream energy level (above crest level)
K and K can be read from Figure 22 and 23.

p a
To check the order of magnitude of such an correction:
when K = .025
p
and K = .1
a
and N = 3 and H /L = .1
e
the correction in crest length becornes:
20 * 0.025 + 0.1) 0.1 = 0.035 (or 3.5 %)

H0 Water surface upstreom from weir drawdown
,7
h = __________
0
2q( P + h0)2
d
,Origin and apen of crest
-- xc --+-, --------x --------
YC
Figure 18 Crest shape: symbols and definitions, from [38]

49
h 0 /Hd
0 004 008 0.12 0.16 0.20
D
0.26
0,24
S.
0.22
020
012
1
0.10
>' 008 u- =11

006
0.04
---- ---
055
0.50
045 ;..-------
0.40
Up
0,35
0.30
- --------
0.25
-- ----
020
0.15
0 0.04 0.08 0.12 0.16 0.20

io h/H
0 004 008 012 016 020
4.
(B) oIues of K
044
h 0 /H
Vert ca

.10
i
St n'
fo~ e - 2 hzontol to 3 ~CticGl
1 ::
(C) vaIue of
1
- 11, / Hd
Figure 19 Crest shape at different slopes of the upstream face, from "Design
of Small Dams USBR" [38].
50
0.75
0
0
0.70
c
q-Cow(2/3)\/ -2- 312
4) 9 Hd
1.)
(when Ho:Hd)
4,
0
0
0.65
4)
0
c
4)
MM ( - t -1h0
ho HO
0.575
0.55
0 0.5 1.0 1.5 2.0 2.5 3.0
0
values of p/Hd
0
0 1.1
in
4..
c
4)
0
4) 1.0
0
0
.4-
0
HO
0
Vj
0.8
D 45
0.20.40.60.8 iii.'
Io ratio of "head -on- crest" to"design - head: H 0 / Hd
1.04 h0 onge with
t
4..
c Hd slope the verticQl
4)
0
.4-
- v 0 1.02
o
u -
0.)
00 1 .00
0
t;
L
t 0.5
b. values of p/H rl
1.0 1.5
Figure 20 Discharge coefficients of nappe-shaped crests at different

conditions and different upstream slopes, (USBR [38]).

51
Per cent reduction

/ of coefficenfs
1.2

. o0 25r:
/
/ J Af T8
BL
Per cent reduc- Percentreduc-
o
7i
1.1
60 cents vs.hd/Ht
:T 1: 88
1.0
cl
0 0.2 0.4
4"
0.9-6+I--- i
_4_(i
e/He pHe 0.6 0.8 0 0.5 1.0 1.5
0.7 7777
1 ~______
l

0 -I--1-J-
1 /
/ t
2011
lJ
R 1 ______
0.5 14 jt1r - - ______ _____
i1 /
_3
04
1 - 2- -
_-Droined jump with dving
-
jet - ____J
c \ N
03
0.2
9 \ 0 --- -
15 - -----
20-
0.1
~
0
"
_ ._____ton5Urf 0 40
-- ojy02l___4 5Q
0
0 0.5 1.0 1.5 2.0 2.5 30 3.5 4.0
Figure 21 Flow reduction by submergence-corrected version from Naudascher

[24], and originally taken from Ven Te Chow [40].
52
1.Sr
00331455A0
Hj
type3a
n ~
14'
67 c res t axis 1 0267
type 1 type 2
c
type2
:rk\
design head 1
. to
0
HL\' type 3
I
cy
L fr.267H1
A0
k
0267 Ha crest OXIS 0267
0
L 20
CI 0.6
types 3 and 3a type 4
0
Et 4 PIER NIOSE SHAPES
_\?'\t____
20 as
riote: pier nose Iccated in same plone
L
as upstrearn face of spi$Iwoy
\\
damensions in parentheses are
for type 3a
0.0
-0.10 -0.05 0.0 0.05 0.10 0.15
Im coeficrent of pier co.itrocti K o
-
0 4.LJ
4
c
__ A
1 type 2c pier
1.4
1
I type 2bpier
type 2a pier
type 2c
0
J
design head
. 1.0 ype 2 pier
N
c '.3
.21
5)
0
3
1 X 5 : 24 08
2 0.8 0 0.5- \ e
67
'.3 0 "A .
5) L
5, 1.0
L type 2b type 2a >
0
c 0.6 0 .ixis of type 2 pier
1.5-
0 type 2 spillway crest
0
0 0
0
5)
- 0.' 2.5L
0
'
0 0 rOtio of horizontal
0 distonce to design head
L
0.0
0.10 0.05 0.0 0.05 01 - 0. 15
- coefficient of pier contraction - K
Figure 22 Pier-contraction coefficient, from WES [43]

ij
basIc eqr.rotorr 12.0

1
(D 0 112
o Cd[L -2 jNK 9 ) (21)1
0
wri ere
O F
loj
- 04
o dschar9e suggested 4 0
Cd lrscnarge coeffrclent
SUgCeSteC desIgn cur4e.1
01' L' liet lengtrr of crest - m
OeSIOn curve rill 4.
N number of prers
1
IL
rt 50
K p pIer contractroo coeffrcrent tO
El
(D K o abutmerit coIltractlon coefflclerit Cr lo
- 0
CE H total heod 041 crest - m
c;t
.180
(3 60
0 leld 0
0 00
legend
rt syrnbol project R WIL. W/H
Pl
9) 1 8
0 o dor000 t 2 1.55 0.96 401
1
symtiol broiect t I WIL WIf-r
/
(3
p doreno 2.1 210 3.77 t2
rt
0j 0
o CW 801 j t 55 0 96
1-" o r rock 1.5 2 67 1 42 1015051 22 2 10 3.77
0 carlyle 2 212 177 IA. 00

0
0 pliulOOtt 1.5 267 1 42 1
Fter f ge t 5 1 2.83 948 20
/0 2.12 1 77
1 pirre flat
1, t. 2
ID
() 02 O certe' hill S 383 948
l
0
(9 gated spllIway wItt pers
goted sillwcy witt piers
0.01 Ln
I-. OC 02 00 02 LO
(2 02 00 02 04 late
note
t radius of dbutmertt - 51
(9 t radius of atiutmerrt - 51
W widtfi of approach 1- eprocuced in
W width of approach reproduced in
model - m
model - -r
gross widtrr of sprllway- m
Lgrass wrath of sprilwot - m
ri debEt of approacll in model - m
4-f design teOd 0(1 clest, M.
H dePth Of approach III model - m
0
S
1 ,
1
1 '1
t-
54
3.4 The napoe-shaned soillwav with 2at
For the weir profile of the former section, the discharge coefficient of a
tainter gate can be read in Figure 24 and the one for a flat gate in Figure
25. The discharge coefficients are based on a pure two-dimensional situation.
The WES [43] presented the discharge coefficients for a sharp-edged tainter
gate as C. The discharge equation reads:

Q = C G w 12g H (51)
G is the gate net opening as indicated in Figure 24, w is the width of the
gate opening and H is the energy head compared to the level of the centre of
the gate opening.
For the flat gate on a spillway crest the equation of discharge starts with
the free flow discharge Q 0 (see for the determination of Q section 3.3) and
then a reduction factor is applied. In Figure 25 it is indicated that the
reduction factor can be written as

Q/Q0 = ( H2" 2 - H1' 2)/H3/ 2 (52)
This equation is based upon a consideration presented in section 1.5, and on

Equations 4 and 5. However, it seems illogical that the gate seat does not
correspond to the spillway crest level, whereas in Eq. 4 the bottom level is
the same as in Eq. 5.
55
110
suggested design curve

for X / Hd = 0
100
suggested design curve

for X / d 0 1 to 0 3
pool
70
for rn u 0
60 0 CG0 wV -
where
G. net gate opening
w z gote wdth
F-4 head to center of gate opening
50 l_ I _
060 070 080 0.90 100
- discharge coefficier/t (C)
Figure 24 Discharge coefficient tainter gate on spillway crest,
from WES [43]
Note 0 = free - flow discharge

0
at head H 0
0.9
QG=discharge at heod H 0
and gate opening G
08 H 1 =H 2 -G 0
0.7
0G
o 0.6
0.5
0.4
lip
0.3 , JH2
02
01
0 DEFINITION SKETCH
0 0.1 0.2 0.3 04 0.5 0.6 0.7 08 0.9 1.0
H 312 - H 312
_ ( eq. 4 )
2 H0312
Figure 25 Discharge coefficient of a flat gate on spillway crest,
from WES [43]

56
3.5 Broad-crested weirs and dikes
At broadcrested weirs the strearrilines above the crest are approximately

parallel and there is a hydrostatic pressure distribution. From the upstream
water level to the level on the top of the crest the Bernoulli equation can be
applied. Two elements cause that a correction coefficient need be introduced:
there is an energy loss by friction and the boundary layer displacement
thickness 8 is a quasi-heigtening of the crest.
The discharge equation for a broad-crested weir with a horizontal crest
(sometimes called a long-crested weir) at modular (free) flow is found from
section 1.5, Eq. a and 6b.
q = h 1 12g (H-h 1 ) (53a)
For the condition of maximum discharge (equal to the conditions of free flow)
hoids also:

= (2/3) (53b)
So the discharge equation becomes (introducing (H 0 -h 1 ) equals the velocity

head)
q = (2/3) H 12g H /3
0 0
(At the crest the waterdepth h 1 equals also the critical depth h 1 = d= q2/g.)
In practice also a correction factor C is introduced so that the discharge

equation becomes:
(54)
Q = Cw (2/3) H 12g H /3
0 0
(w = width)
The discharge equation 54 is a theoretical one; a small correction is needed

due to the boundary layer effect.
Naudascher [24] gives the following theoretical corrections for a long-crested
weir with adequate flow convergence, without extra flow contraction or other
losses. It is based on a reduction of the net flow section by the displacement
thickness of the boundary layer:
6d
c = ( 1-2 ) 0- ) 3/2 (55)

w H
0
57
The displacement thickness of the boundary layer depends on the crest length L
and the Reynolds number (related to crest velocity and crest length).
(56)
for ReL < 3 10 5 the 8 d/L = 1.73 .

for ReL > 3 10 5 the &d/L = 0.037 Re (57)
To give an order of magnitude for the necessary correction for the last
condition he presented the following table for C values, for the ratio W/L =
0.2.
H/L
L(m) 0.05 0.075 0.10 0.125 0.15
0.61 0.883 0.923 0.941 0.953 0.960
1.83 0.923 0.446 0.958 0.966 0.971
3.50 0.932 0.953 0.964 0.970 0.975
15.25 0.953 0.967 0.975 0.979 0.982
Table 3.5.1 Correction coefficient C in Eq. 54.
The downstream waterlevel has only effect when the induced pressure at the
downstrearn crest edge is higher than the hydrostatic pressure belonging to
the modular flow condition. This occurs when h2>hd (see Figure 26 for
notat ions).
Only when the back face of the sill is vertical, the submerged-flow
pressure conditions downstream of this face is assumed to be hydrostatic.
The waterdepth h 2 can then be calculated with the momentum equation (see
section 1). The effect of the downstream water level is expressed as a
reduction which is a function of the submergence factor:
S = h 2 / h 0 (58)
Here h 0 and h 2 are the up- and downstream water levels above the crest
level. The downstream water level at which submergence has influence, is
higher than the critical depth hc, being 66% of H 0 , because of the recovery
58
of potential energy. This recovery is stili more when the downstream sili
face is sloped. In Figures 27 to 29 some data about discharge coefficients
of broad-crested weirs and dikes are shown; they have been taken from DELFT
HYDRAULICS [10], [11] and [12]. At the dike model where also skew
conditions are tested (Figure 28) the downstream flume walls do not
correspond with the natural direction of outfiow, so the discharge
reduction factor CRS (by submergence and by the skew approach flow) in
prototype might be somewhat greater than the measured one. The effect of
the skewness in Figure 28 is als greater than what is expected from the
considerations of Schnfeld [35] mentioned in section 1.2.
Figure 26 Notations: flow conditions at a broadcrested weir.

w
0
0
0
0
L
free flow rIations

rduction by submergence 0
submargnc rotio S ()
DISCHARGE COEFFICIENT
CD 2 1.5 0.5 15
and Q:CS 00
.g. BH0
11
=NMI
10
09
1k%UIll1UUL%tt1 1
07
h 0 /p
06
0.5
0.4
03
iIhii__
02
0,1
0 T
0 02 0.4 06 0.8 10 12 14 1.6 18
, 0 0 /L
110w puttrn unstabI
bgInning of 110w sporoton
Figure 27 Discharge relations of broad-crested sills with variable sili

length, from DELFT HYDRAULICS [11]
Q r CR S d W
V' ' H 0 312
1.0 I.0
tf)
u
0.8 0.9
L)
0
0
0 0
0.6 0.8
Gil-
0
0
0)
0.4 0.7
0.2 0.6
0 0.5 L -
0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 0.7 0.5 0.9 1.0
) 0 31 h 2 /H0
p/b w 0.75 -1.3

H 0 /b = 0.3-1
h 0 "b0.2 -0.5
Figure 28 Discharge coefficients of dikes, (model on scale 1:25) from

DELFT HYDRATJLICS [10]
61
he gi ven re la tionievalid for:
120
-vpntrerlepieg face :2
- 1 treon (opleg foe 1:3

- 2 0.10 m o::l 0.10
- 11 1 /P 1 0 (.0 112 0 1.0
- ( ( /0 4 II) .*ed 3 ((.30e

3.10 Ch oeenended oreerleegth L ? 0.10
2.0 2.5 3.0

Cd
h,j ]HO h2
O .Sh'J___________________________________
1 o o vwc--)g 5
H0
1.5
0 0.5 1.0 55
H0 /L
DEC
0
0
Iirn:t of ee(lnb.Ihty
0
0.7(
0
/
/
rncornrr,qndat rrnae vCIuQ / 1
0.81 -- - 7-7'
t /
o.
1.0(
0 X W 70 00 93 96
) s mrgflcu rotk2 S 3 142,4f.)
SUB1IERGED F1.OW
modular linit
70
60
S = 100 ( 11 2 /Ilo) 2
FREE FT.OW
50
2
-1
7 H c /L
Figure 29 Discharge coefficients of different types of silis with sloping

front and back face, from DELFT HYDRAULICS [121
62
3.6 The morning glory spillway
The morning glory spillway is a spillway which is circular (or a part of a

circie) in plan. It is an option when a tunnel is applied in whatever way,
for instance a diversion tunnel left in place after the dam has been built.
The compactuess of its shape is favourable for its application. To prevent
siphoning of the whole tunnel, which would cause a dangerous situation at
extremely high discharges, a control section is created underneath the
spillway with, for instance, a sudden widening of the tube section, or
where air is let in to guarantee that atmospheric pressure occurs at that
level. The level of this air inlet cuts off the discharge, and even when
the inlet is fully submerged, velocities at the cut-off section are
restricted to v"-
2gh where h is the difference in level between the upstream
energy head and the cut-off section. Figure 31 shows an example of how a
cut-off section is introduced.
It is not advised to apply the morning glory spillway for conditions where
prediction of a maximum required discharge capacity is still very
uncertain. At higher water levels Q increases with IIi, whereas for a
straight spillway the discharge increases with H312. Therefore, the risk of
overtopping the dam is smaller in the latter case.
In Figure 30 is shown the different types of flow which can occur when no
aeration is applied.
At condition 1, the discharge is determined by the overflow capacity of the

crest.
At condition 3, at full siphoning, the total head difference H is
determinative, Q V'2g HT.
At condition 2 the level difference H determines the discharge. H is
a a
between the upstream level and the highest point where air at atmospheric
pressure occurs. But when air pressure is lower than atmospheric, the
effective H is greater.
a
At condition 2 the Ha moves downwards at greater discharge, because the air
is sucked out gradually and the situation is very unstable.
Without the introduction of a control section by widening and/or aeration
it is in fact the downstream tube-end which, in condition 3, is the control
section, while with aeration at the point "orrifice control" of condition 2
the situation is stabilized at that level. See further the discharge curve
of Figure 31.
63
Like the nappe-shaped spillway of section 3.3 the design of the crest of
the morning glory spillway is based upon the atmospheric wall pressure
condition related to a design discharge. And as in the straight spillway,
this condition is established by measuring the lower nappe profile of the
free overfali nappe of a sharp-crested weir (in plan circular) with radius
R, see [42].
In the design of the morning glory spillway the requirements lead to a

compromise. This design process is an iterative one. In the final design
stage model investigations are often applied. The following elements should
be taken into account.
a At the condition of submerged crest and discharge control at the

control section, the discharge has to pass any section without a wall
pressure lower then atmospheric, so the distance z below the upstream
water level follows from the conditions that z equals also the velocity
head; so
Q= 1r R2I (59)
hence at any level
R . = (Q/ (2gz) (60)

min
b As a first estimate in the iteration, the design head can be taken

equal to the demand of the design water level departing from free flow
and maximum discharge. With Figure 32 the outer radius R 5 can be
determined. As a next step H can be read in Figure 33. R and H are
5 5 S
the values which were tested in the sharp-crested weir test of Ref.
[42]. The discharge is free flow when H/R ~ 0.45, and is submerged
and R the weir shape can be designed
when H /R > 0.9 to 1. From H 5
5 S 5
from Figure 33.
c The location where the curves, discussed in a and b, meet can be a

first choice for the level of the discharge cut-off (control) section,
but due to a discharge coefficient smaller than one (0.9 is advised in
[38]) the control section is a little wider than the corresponding

section of the curve obtained in a. Then the transition between the
curves a and b is smoothed. Condition a can lead to a local or to a
total widening of the upper part of the bell mouth as found in b, but
this is not so when H /R is smaller than 0.225 to 0.3.
d s
64
The procedure does not lead straight forwardly to the final design, so
other values of H and R have to be tried also.
d s
The correspondance of the area of the control section with the tube sectjon
further downstream, is also a point of consideration.
However, when other Hd values are tried, the real energy head H is not
equal anymore to the design head. Figure 34 gives the correction for the
discharge (but only for Hd/R = 0.3). In [38] it is advised to use this
curve for other Hd/R values also. But one should consider that at small H
values the discharge relation approaches the one of the broad-crested weir.
Therefore, in Figure 34 could also be completed with reduction factors (at
a small head) for other Hd/R s values than 0.3; they are calculated from
Figure 32 in combination with the theoretical discharge equation for a
broad-crested weir of Equation 54.
Another design element in the last stage of the design is related to

cavitation. In contrast to a straight spillway the discharges smaller than
the ones at the design head, cause wall pressures lower then atmospheric.
Wagner [42] also performed tests with lower pressures than atmospheric
underneath the nappe and these results could serve, by comparison with the
designed shape, for estimation of the related pressures.
65
H0 0 'es H.
1Li
T,- otston tube
TOroot of tr'oflstcrr tube
ltCI,r,ed shaft
Outlef leg of Condtat
Accelerott9
f 10w
Oeceleratng flow
CONOITION L CREST CONTROL
- - -' t ,
-T - asto tube
- OrfCe Cottr'Ol Ot thrOOt of trots,t,on
CONOITION 2 TUBE OR ORIFICE CONTROL
, HydO 1 c 9 ad eet
h
f1,
Ta,lwoter
CONOITION 3. FULL PIPE FLOW
Ppe control, 0 0 (II,- hL ) Condrtron 3-
Pont of change from oriflce

control to pipe flow--
- h '-Errotc tlo
- Tube or orifice Control, 0 F(H 0 'e), cotd,ton 2 ,' rarrge
It
51
1
frlead af wh,cpt tunnel fiows- 0 75 full
1 0 f downstream end -
Po,nt of chotge from wmc to orifice 0ntr0l-
b
- --
C
- Cr'est control, 0 f(H, ),onditrot t
DISC HARGE - SECONO-FEET
Figure 30 Nature of fiow and discharge characteristics of a morning glory

spiiiway, from USBR [381
- R8OO -
--
-00
81:
81
811
>
41
811
80
80
spulwoy 1
0 spilwoy 2
80
0 100 200 300 400 500 600 700 800 900

-* dlschar9e per spulwoy in rn 3 /s
Figure 31 Control section introduced by aeration and enlarged pipe profile

(design Coyne and Bellier, from [4])
67
0.8
0.7
0.6
D )15 NOTE: dotted lines are
1 sedon extrapo1aton
0.5
tr 0.4
r'J
NI
11
0.3
0.2
0.1
0 0.5 1 1.5 2
HdIR52
Figure 32 Discharge relation of the Morning Giory spiliway at the design

head, dimensionless redesign after USBR [38]

W.
L=
1.1c
1.08 4 III IIIIIIIIIIIIIIIILI
1.08
-4 5
NOTE Dat ed Ilr,CS are
Hd - -03 bosed 0fl eetraataton
O data
1.0
-t- .-- -- 1 - ----
1 .0
0.1 0.4 0.8

Rd
1.2 1.6 2.0
TI
H 8
Relationship of - to - for circular sharp-crested weirs.
d s
0 05 1 5 20 05 30 33
- --
05
0b0.
10
\
55
't \
0
-20
25
-30-- -
H
15
Mff
Upp, ,ppe
,..qIigibte opp.00d, I0c1Iy).
OPP 05d
Figure 33 Lower nappe profiles of the morning glory spiliway,

from TJSBR [38]
M
.
1.2
(Y)
curve for Hd/Rs 0.3

T
t' 1.0
I
0
L
6.922
g
0.824
-
0.795
0.6
0 01 0.2 0.3 0.4 0.5
He
R5
Figure 34 Reduction of crest-discharge coefficient for heads different

from the design head; the curve is for Hd/R = 0.3 (from USBR
{38J)
ce
3.7 Labyrinth weir
The labyrinth weir is used for increasing the weir length within a limited
width of the opening. In the following the performance of the sharp-crested
labyrinth weir is presented, following the recommendations of Hay and
Taylor [ 1 6].
Figure 35 shows the Labyrinth weir in plan, and the applied symbols.
The following parameters are importance:

L/w = Length magnification
QN = Normal discharge of straight sharp-crested weir with width w
QL = real discharge over labyrinth weir
h0 = upstream waterdepth + lmm (like eq. 43)
p = weir height above (upstream) bottom
a = angle of the flanks; a equals a at a triangular shape
max
S = submergence factor = h 2 /h 0
R = bottom step (downstream bottom lower)
When h 0 /p is small, the extra weir length tends to be fully effective and
globally Q L/Qn = L/w within certain design limits.
The following recommendations are taken from Hay and Taylor:

- Crest-length magnification, L/w: 1f operation of higher head related to
crest height is envisaged, then length magnification greater than 6
gives little return when designing for h/p ratios exceeding .25.
- Vertical aspect ratios, w/p: It is recommended that the adopted value of
the vertical aspect ratio should not be less than 2 in the case of
trapesoidal weirs (in plan) and not less than 2.5 in the case of
triangular weirs. Nappe interference effects will be negligible on a
weir designed in accordance with these recommandations. Where p is very
high, the value of w/p ratio may be less than 2 providing the ratio on
the head to the cycle pitch h 0 /w does not exceed 0.25.
- Side-wall angle a: This factor is of primary importance in determining
the performance of the labyrinth weir. Furthermore, the correct choice
of the side-wall angle a does not affect the structural costs. For
maximum performance, the greatest value of the side-wall angle should be
adopted in the design, i.e. triangular plan-form weirs should be used
wherever possible.
71
1f the triangular plan-form is unacceptable for other reasons, a

trapezoidal plan-form having a side-wall angle not less than .75 times
the maximum value (given by the triangular plan form), may be used
without incurring a large loss of performance.
- Channel-bed elevation differences: Downstream interference which has
small detrimental effects on performance can be reduced by decreasing
the elevation of the downstream channel bed, and can be completely
eliminated for the operating range recommended here, if the difference
in channel-bed elevation is equal to or greater than the maximum
operating head.
- Aprons: Both upstream and downstream aprons are detrimental to
performance. However, they may be necessary for structural reasons. In
situations where a fall in the channel-bed elevation occurs on the
downstream side of the weir, the size of the downstream aprons will not
or nearly not affect performance, provided the channel depth at all
points of the downstream channels is greater than the expected depth of
flow in the upstream channel.
- Crest sections: Under high flow rates labyrinth weir discharge tends to
be, to some extent, independent of the weir-crest coefficient; the use
of complex, expensive crest sections is unnecessary from the hydraulic
point of view.
- Submergence: Because labyrinth weirs operate under smaller head than a
corresponding linear weir discharging the same quantity of water, in a
situation which usually involves operations under drowned conditions the
use of a labyrinth weir will increase the degree of submergence. It is
not recomrnended that labyrinth weirs be used where they would usually be
subject to operatlon under heavily drowned flow conditions
It has been found by Hay and Taylor that the Villemonte equation for the
effect of submergence (Eq. 45) can be applied here too:
Q = Q 0 [1 - (h2/h0)'.5].385
Figure 36 and 37 give the discharge relations for the triangular-plan and
the trapezoidal-plan labyrinth weir for horizontal bottom and for deepened
downstream bottom.
* Authors comment: This is not literaly true because at submerged flow a

labyrinth weir placed over the whole canal width, downstream of a narrow
sluice, can serve as discharge spreader.
72
w C.
c ho nnel
yj mmy
channel
L:4c.2b 0
developed crest length
PLAN
Figure 35 Labyrinth weir in plan and symbols.
Lfw=6
QL
N
- - - - no downstreom hterference -
2
downstreom Interference
(zero chan9s in bed &evatlon)
0.1 0.2
IN 0.3 0.4 0.5 0.6
h0
Figure 36 Discharge amplification for the triangular plan-form weirs

without aprons with w/p ~ 2.5 and a =
a max
From Hay and Taylor [16].
73
LI _
0L
UN

-- - - no downstreom i- terference
downstreom Inter ference
(zero chonge in bed elevatlon)

0 0.1 0.2 0.3 0.4 0.5 0.6
ho
Figure 37 Discharge relation for the trapezoidal planform weirs without

aprons with w/p ~ 2 and a = 75a. From Hay and Taylor [16].
74
3.8 Gates with underflow
In this section only two-dimensional situations will be considered, without

shafts and slots. In this section results from different authors duplicate,
but the reason of referring to these (nearly equal) results is that they
served for each of the authors as the basis for further analysis;
therefore, it was considered only fair to present also the results from
where they started.
When the upstream energy head is H and the pressure head just downstream
of the gate is h 1 , then the Bernoulli equation leads to the following
discharge equation.
q = C a 'I(2g (H 0 -h 1 )) (61)
c
C is the contraction coefficient of the jet, a is the gate opening.

In the case of the free jet the downstream pressure head h 1 equals Ca.
The flow contraction depends on the geometry: this results in an influence
of D/a in the case of a culvert with height D or, at free surface flow, of
h/a (where h is the upstream waterdepth). But the flow contraction also
depends on the pressure distribution at the top of the jet contraction.
When the jet is subrnerged the piezometric head there is constant; at a free
jet the pressure itself is constant (atmospheric). The influence of the
gravity can then be expressed in terms of the Froude number.
Fr = V /I(g Cca) (62)

c
(Vc is the velocity at the jet contraction)

When Fr is great enough the gravity influence disappears and the submerged
and the free flow have the same contraction. This is found both in the
theoretical study of Rouv and Khader [33] and the experiments of Nago [23]
presented hereafter.
When the upstrearn water has a free surface, the Froude number and the
upstream waterdepth are coupled and Cc can then be expressed in terms of
h/a only. At free outflow from a culvert Cc depends on both D/a and the
Froude number, as shown in Figure 38.
The Cc can be derived theoretically, as for instance published by Rouv and

Khader [33]. Friction effects are neglected in the derivation. The interes-
ting presentation of Figure 38 is not directly suitable for practical use
75
because the C value and the discharge are both needed to calculate the
Froude number, and then C depends on this Froude number. In a design graph
the relation a/D, a/d and a/H 0 should have been presented (d = critical
depth, H0 = energy upstream). A global check consists of recalculating at
different points in Figure 39a the Froude number and the contraction
coefficient C and comparing these resuits with Figure 38: no deviations
were found.
When the contraction coefficient is known, the free flow discharge can
be expressed theoretically by application of the Bernoulli equation as in
Eq. 61.
From dimensional consideration it can be seen that Eq. 61 and the influence
of the Froude number can be combined in a single presentation.
q = Cd av"
2 gH (63)
Cd is a discharge coefficient.
But in general the following expression is used:
Q=CwaI(2gh Q ) (64)
It can be theoretically proved, see section 1.5 Eq. 8, that the discharge
coefficient Cd of Eq. 64 equals:

Cd = c/'/(1 + Ca/h 0 ) (65)
The advantage of the use of Eq. 64 and 65 over Eq. 61 is that no iterative
procedure is needed; from C and a/h 0 the q can be computed.
Cozzo [9] investigated the contraction for sector gates and flat gates with
free surface flow and found systematically that the angle near the edge was
determinative for the discharge coefficient (Figure 39a). The flat gate and
tainter gate show the same results. Figure 39b, derived from the curves of
Figure 39a, shows the discharge related to the upstream energy head.
In Figure 40 the discharge relation for free flow is presented for a tam-
ter gate located in a roofed culvert with free outflow, from WES [43]. The
upstream pressure is expressed in terms of an energy head in the culvert.
In fact, it is the contraction coefficient which is presented in Figure 40.
ffel
From the rnodular flow relations of Cozzo the limiting conditions for which
the downstream water level does not affect the modular flow discharge is
presented in Figure 41. This relation is found by calculating the
contraction coefficient C , from the discharge coefficient Cd (the inverse
of Eq. 65), and then calculating the conjugate depth of the hydraulic jump.
The C is from h /a = 2 on nearly constant.
c 0
In Figure 42 are presented calculations and measurements of the submerged

discharge relations at a vertical sharp-edged gate plate. The measurements
of Henry are found in Rouse [31], and are also discussed in Naudascher
[24].
The resuits of the computations presented by Rouse/Naudascher [24, 31] were
based upon the theoretical values of the contraction coefficient similar to
the ones presented in Figure 38.
Descrepancies were found which were attributed to disturbances like uneven
flow distribution, and so on. The author performed new computations,
presented in Figure 42, based upon slightly smaller C values (see Figure
42 also), and he found a fairly good agreement with the measurements. The
use of adapted Cc values is useful when the discharge relations of complete
structures are performed with the computer program of Chapter 2.
There exist also measurements which are in agreement with these lower
values. For instance, Naudascher [24] presents similar measurements as the
ones of Cozzo, but then from Gentilini. These results have for great h 0 /a
(about 12) values a Cd of 0.585, corresponding with a Cc of about 0.6.
Nago [23] performed a very systematic research on the free flow

characteristics of different types of gate edge. He distinguishes also, as
it was presented by Rhouv and Khader in Figure 38, a flow contraction plus
an effect induced by gravity. Because Nago only considered free surface
flow the results are simpler. The contraction coefficient he obtained
theoretically. He presents the effect of gravity as a correction which can
be added to the contraction coefficient at submerged flow. The correction
is a function of a/h 0 (Figure 44). The calculated discharge coefficients
are compared with measurements, Figure 46.
The results for the three gate types of Figure 44 are presented in
Figure 45.
Boiten [6] investigated vertical plates with round edges at underflow

conditions. The shapes are presented in Figure 47. The discharge relations
presented in Figure 48 show that all the circular shapes have the same
77
discharge coefficient. It is observed that from h 0 /a being about 2.5 that

this parameter has nearly no influence on the contraction coefficient of
the investigated circular shapes. At smaller values this parameter has an
important effect but measurements are then also scattered, due to air-
entrainment.
The discharge coefficient can be much higher for the rounded edge than for
the sharp-edged gate. It is remarkable that Boiten found a contraction
coefficient for the sharp-edged gate which is higher when compared to the
resuits presented above.
The results of a series of tests of reversed Tainter gates in completely

submerged culverts (Pickering [27]) with and without divergence in the
downstream culvert part, are presented in Figure 49 and 50. This is in this
part of the manual one of the few cases which represents internal flow. But
the reversed tainter valve is typically in use in a hydraulic structure
(high head navigation lock) and nowhere else.
Re-*oo 1
Fr= \
0.7
Vc cc
/7 / / / /
cc c
----
0.5
0 0.2 0.4 0.6 0.8 1.0
alD
Figure 38 Contraction as function of Y 1 /a and Froude number, from

Naudascher [24] and after Rouve and Khader [33]

T \
Q Cd a V12 g h 0
- 30T o
0.80
PLE
do.70
.
0.60
.4
3
77 7171 :
4 6 8 10 12 14 16 SECCR GAYE
>
Figure 39a Discharge coefficient for free outfiow as function of plate

angie at the edge, from Cozzo [9]
09
CdH
0.7
uji;.i' 1!
0.6
TJ
05
2 4 6 8 10 12 14 16
H0 /o
Figure 39b Discharge coefficient related to the upstream energy head,

derived form Figure 39a
79
bosc Qguatlon 10
o = CdGO B\/
whr
O = dischargQ
Cd = discharge cofficiQnt 0.8 DEFINIT ION
i
GD= 90t12 Opfliflg
p,
B = width of gate Oponing
II
H = snrgy grad. QII2v.
- (inwrt QIv. CdGO) 06
Remork. in this dehnition
Cd equals the flow contraction
coefficient in eq 7 and 65
-
(0
FA
0.
o 0.4
00
1L'iii
0.2
LEGEND
von Miss (computed
Garrison n-)jdql
0
0.60 0.65 0.70 0.75 0.80 0.85
-* Cd
Figure 40 Discharge coefficient for free outflow of a tainter gate in a

roofed culvert, after WES [43]
09
08
07
0. r 30
os
2 0
05
0 9O
04
0.3
02
2 3 4 5 5 8 10 15 20 30 40 50
0. 11 0/0
Figure 41 Limit condition of modular flow derived from the discharge

relations of Cozzo (Figure 39a)
Cc used in computotion
To;
0.6
0.5
> a/h 0
ho
1
q = C d a/2gh0 computed
meosured
0.7-
0.6
0.5
0.4
Cd
T:: oilF111Ff4 __
0.2
ME _IwlilIliffi__
0
0 2 4 6 S 10 12 14 16 18 20
) h/o
Figure 42 Discharge relation for submerged flow of a sharp-edged vertical

plane gate, calculated with a best-fit relation for the
contraction coefficient compared to measurements by Henry, as
indicated in Rouse [31]
81
hJ
j2
U_ cc o
TYPE -1
q = Cda
ho
Cd C 114 Ca/h 0 (free flow)

CC Q
TYPE - 2
c c fFT
c
h 1h /11
0
- ( C c a/h 0 )' (sribmerged flow)
d
C
c
- (C
c)0
/ C C (free flow)
h' lip h C = (C) (submerged flow)

Lj
TYPE -3
Figure 43 Type of gate edges investigated by Nago and notations [23]
correction for the gravity

x Cc -0.05 (at submerged condition
i CC: C - (Cc)0 .Cc:O and C:(C) 0 )
-0.101
0 0.1 Q2 0.3 0.4 0.5 0.6 0.7
a/h
Figure 44 Correction coefficient for gravity effects on the contraction,

from Nago [21]
EM
r
HEW
(C). (Cc)o
0.8
0.7
0 0.2 0.4 0.6

a/h 0
0.6
TYPE- 1
) o/h 0
TYPE - 2
1.0
45 o
( ).9 Cce
(Cc)o
-
,// /////
0.01
=
0.5
/7/
Q2 0.4 a6 0.8 1.0

) 0/ho
TYPE - 3
Figure 45 Contraction coefficients without gravity correction,

from Nago [23]

83
\ -f- 0
T
hjI
N . ,.;T-i
0.80
30
910
Cd
CL7 0 ho
C CM
77 Cd 065 9
0.60 inchned
gote
[J d1a,
Go 0 0.1
-
1.0
o 10
1 4.0
0.5C ' ' 1
0.55
0 0.1 02 0.3 0.4 05 0.6 0.7 0 01 0.2 0.3 0.4 0.5 0.6
oI h 0
Q /h 0
TYPE - 1 TYPE- 2
L5
0.7 h0
0.11
0 0. 20
Cd o 0.30
o / h0 0 . 11
0.6
0.20
0
G- O~0
*
033 0
9
9
0
lip length
r.11
0 0.5 1.0
LIG
TYPE - 3
Figure 46 Discharge coefficient compared with measurements, from Nago [23]

84
0.3 d
type A

EI type B

type C
Figure 47 Shape of round edges as investigated by Boiten [6]

85
10
2 8 (
E
vvv d/ar2.1 (c)
0.9 4(a)
64 ( b)
18 a.b)
x1.2 (c)
0.8 = - = 0.84 (b-3--
Cc
0.7
_______._ 0.25(c)
' sharpe edge (d/ar0)

0 a
06
( )
air entr&nment
0 THH I I
.5 4 13 b /1U
h/a
m
oLe
cc
~ 07

0 1 2 3 4
Limit values of C c from Fig. A
Figure 48 Discharge relations for round-edged gates [6], A measurements

and B the proposed relation
TYPE t
J
'
TYPE 2
TYPE 3
TYPE 4
86'
TYPE 5
1 /0
36 14 / 68
TYPE 6
J - -t ,
/0
56
TYPE 7
to
56 224
TYPE 8
Figure 49 Culvert configurations with expansions and reversed tainter

valve investigated by Pickering [27]
87
12
11
I-
Numbers of lines
,
refer to lig. 49
ir:
>
4:9 2-
4)
0
4)
0
0
7
0
1,7,8
t,
>
\ V02
VO
o
30 40 50 60 70 80 90 100
, percent valve opening
effect of culvert exponsions on valve

loss coefficients
Figure 50 Energy loss of reversed tainter valve in culverts with and

without divergence, with submerged flow, from Pickering [27]
3.9 Gates with overflow
Two examples of discharge relations of overflowing gates are presented, the

drum gate in Figure 51 and the flap gate in Figure 52. The discharge
coefficients refer to the following equation:

q = Cd (2/3) V(29) H 0 3/2 (66)
Resuits of a model investigation of a flap gate on top of a vertical

(Stoney) gate are presented in Figure 53. The angle a was varied between 27
and 68 degrees. As reference the results of Sarginson are also presented
(see section 3.2).
As for other gate types it might be that, depending on the crest shape,
some discharge data presented in section 3.1, 3.2 or 3.3 can also be used
for overflowing gates.

E;1,J
50
iii
4)
1)
L
20
/
straight inclined
- weirs
IC
Sue
/
S ' /
0.1 02 0,3 04 0.5 06 0.8 1.0
rutio
-20 r
0.55 0.60 0.65 0.70 0.75
) Coq/VHo3/2
Figure 51 Discharge coefficient of a sector-shaped raising drum gate, from

Ven Te Chow [40] referring to original work by Breadly [8]
o 2040 60 0
800 070
-20 0 0 20 60 800
L 125
0,65 0.70 1.00 0.65 0.70
60 0 80
c
a60Lr4f 0 .650 .70
0.55 0.55 0.60 0,65

'-"-
d
-20 1 0
,
60'
1 PLO
: 0.0
.
055 _r! 02 055
060
_1__ 0,50
-200 0.25.1J 40 60 80
0.10
0.55
-' 2i]
4I0!
2 ':::
OL
QzwCd (2/3 )\ H'2
Figure 52 Discharge coefficient of a fully aerated flat-flap gate, from

Naudascher [24] referring to original work by Bss.
flow locally free
Q from crest
3
4 c - 27 49'
- 1 iJ-
: 37 0 30 1
-
V
/
2 r,9
X :47Q
/
r573 OVX
v 1
o :62.30'
/
o :6730'
x 01
/
1x /
NSARGINSON
0.5

meosured ot
0.2 R80mm -
0.2
0.1
0.50 0.55 0.60 0.65 0.70 0.75 0.80
) Cd
Figure 53 Discharge coefficient of a fiat-fiap gate with circular crest,
from Boiten [5]

91
3.10 Gates with combined under- and overfiow
The discharge relation of gates where simultaneously over- and underfiow

occurs, is not fundamentally different from the situation where only one of
the two occurs. Whether the situation is suitable for an estimate of the
discharge depends on how exactly the conditions just downstream of the gate
are known. The overfali nappe should be fully aerated, but that is also the
case for nappe flow without underflow.
Naudascher [24] discusses an example of a calculation where the momentum
equation is applied downstream of the gate under the assumption that just
downstream of the gate hydrostatic pressure occurs, even at the sloping
bottom. This is only justified when no recovery of potential is to be
expected.
Figure 54 shows the conditions, the symbols and the momentum equation as it
was applied to compute the water-level difference (h 2 - h 1 ). The coeffi-
cients Ca and Cc related to the nappe and the jet respectively, follow from
the nappe thickness (in vertical sense) at the end of the gate crest and
the contracted jet thickness (contraction is about equal to the discharge
coefficient).
Figure 55 shows the computed resuits compared to measurements.
H
T
bx L LY
/,/ I /////
(q 0 + q 0 ) 2 q 0 2
(h1+Ay)2_(h2+y)2p
h2+y
Figure 54 Application of momentum equation at combined overflow and

underflow, from Naudascher [24]

92
measured for
0.1.
- V --0,571
y3 - y2 o 0.857 V .=0.77
y3 o 0.2
0.4 0
'&0.15 ,
h 0.3
- 0 0.6
)0.32 -
t 06
B
/ ):054
-
FO 30 -
146 0.2
0
0
o computed for
0.1
0 2 .0.457
0 0.32
n
1.0 1.2 1.4 1.6 1.8
yo I y3
Figure 55 Difference in water level near the gate and further downstream;
calculated and measured, from Naudascher [24]
93
3.11 The Howeli Bunger valve
A Howeli Bunger valve (also called hollow-jet value or cone valve) is

applied at the downstream end of a tunnel. It has the function of spreading
the jet. There is also a very good aeration and hence the cavitation risk
is reduced. The spreading of the jet is favourable for the stilling basin.
The central cone is fixed to the pipe through supports which are internal
vanes. The cylinder which moves around the pipe is the valve, see further
Figure 56.
Gieseke [14] performed the theoretical study on the discharge relations of
the Howell Bunger valves and compared the results with measurements. He
applied potential flow theory and assumed radial symmetrical flow.
The gate opening (S/D), the cone angle (2a) and the cone diameter (D-2b)
was varied (see Figure 56).
For a = 450 and b = 0 the computations were verified with different
experiments, which showed a good agreement (Figure 58). Figures 57 shows
the results of the computations with systematically varied parameters. The
values are discharge coefficients related to the tube section and not to
the gate opening.
a
el
closed position
ote cosed
movable
cylinder
-
cone
VO
D
T -
Lribs
"gote opened
_________
OC
b>O
I1 -
bO
OC
_
ko
b<O
Figure 56 The Howeli Bunger valve and the schematisations investigated by

Gieseke [14]
95
Q 0
iCt Ct 0
4
1.00 1.001
a9c 0.90
080 UUPIPAili!4 0.80
0.7C 0.70
Q6C diUdUj 0.60 NW(iI!iUU

0. SC IUl4P1Ii 0.50
040 liLIiUUIU 0.40 UVZIPiIIRI

0.30 Ip,i'I.I 0.30
0.20 IiAPiUUIIUI 0.20
4IIIII 4AIIIIiI
uumimumi
0.10 0.10
0 IlRIRUulP
0 0.20 0.40 0.60 (180
-J-
1.00 a
0
0 020 0.40 0.80 080 5a
0 D
0
Ctr
0.90
CIBO
0.70- I ?
e l I-
/W
o6oj// I
05C-- -
040 - -
c
0.30 - - - - - - -
0.20
b/D:-O.2
0.10-
0L_ 0 0.20 0.40 0.60 0.80 1.0)a

0
Figure 57 Discharge relations of the Howell Bunger valve for b/D=O,

for b/D = -.2 and for b/D = +.2, from Gieseke [14].
W.
0
- cCcuated
Ctr experirnentul Jonschin
- - -
Vermiton
- 6Sditebern
0.90
0. SC
0.70 --- - :45

b/D:O --
0.60
050
0.40
formula
Knapp [13] -----
0.30
0.20
-/--
0.10
/-
0 - 020 0.40 0.60 0.80 tor.
0
Figure 58 Comparison of calculated values with experiments,

from Gieseke [14]
97
3.12 Conduit outlet
The free flow loss of the conduit outlet depends on the Froude number
similarly to a gate with underfiow at free flow conditions (section 3.8).
1f the water in the fully filled conduit fiows into the air, the discharge
capacity of the whole tube depends on the level at which the atmospheric
pressure can be assumed to act. At high discharge (great Froude number) the
influence of the gravity tends to disappear and the centre of gravity of
the pipe section should be taken as pressure head. Figure 59 shows resuits
for a circular conduit, published by the US army Waterways Experiment
Station [43].
0.
0.
0.
0. Nam =11111111 logo
Yp
UliuIuuII
__ RUIIIIIIIII
0 0.
NEI ....ii:1uHuu
iuiii__iuuiiiiiiuuu
iI __
Emuuu__uiuuiiuuiuuin
IUIIIUIIIIihiIIi
M1111011111111llIllIl
1 2 3 4 5 6 78910
v 4Q
Fo7 1IVD5I2
LEGEND
PP(S5iR( AI(NT SYMBOL DATA SOURCE BOTTOM SUPPORT
EX/T POP T4/ -0---- STATE UWIVERSITY OF IOWA WONE

o DENO,OW MODEL LEVEL
~0
DENISOW PROTOTYPE LEVEL
CAPRSON MODEL PARABOL?C
S Y000HIOGMENY MODEL t 0W 20
T ENID PROTOTYPE PARABOLIC
-Z:Z A FOPT RANOALL MODEL LEVEL
FORT RANDALL PROTOTYPE LEVEL
+ OAHE PROTOTYPE PARABOLIC
Figure 59 Level at which the pressure head is zero, at a circular conduit

exit, from WES [43].
4. LOSSES OF COMPONENTS UP- AND DOWNSTREAN OF THE CONTROL SECTION
4.1 Orifice and intake losses
When the flow contraction is known, both orifice and intake (or entrance)
losses can be calculated with the momenturn equation in the region
downstream of the section with maximum flow contraction (assuming
hydrostatic pressure at this section). Or, at submerged flow, the Carnot
equation can be applied. The flow-contraction coefficient in the orifice or
in the entrance is C with an average velocity in the contraction V and a
culvert velocity V downstream (Figure 60). Now introducing for the average
flow velocity in the orifice V the Carnot equation is:
AH (V V)2/2g (V0/C0 V) 2 /2g (67)

= =
- -
= culvert area
A0 = orrifice area
V. = Q/A 0 V =
'///////////////////t
A]
~
A V
Vr V0/t
orrifice losses
Vm r V IC
entrance losses
Figure 60 Orifice and entrance losses
When the culvert is wide, compared to the orifice, the influence of V can
be neglected. When the culvert section is not widening (intake condition)
then V = V and the entrance loss becomes:
c 0
99
_1
AH = 1< V2/2g = (C - 1)2 V2/2 g (68)
e e c c
The flow contraction coefficient C is derived from the orifice tests

c
presented in USBR [38] and shown in Figure 61, series 1. The conversion
from the discharge coefficient C from Figure 61 to the contraction
coefficient C is done as foliows:
c
At low velocities (low Froude number) it can be derived that:
Ah AH
The discharge relation of Figure 61 can be read as:
AH = Q 2 /C 2 A 2 2g
From (67) it follows that:
_1
AH = (V2/2g)(C - A/s) 2 = (Q 2 /A 2 2g)(C-' - Als)2
0
Combining the two expressions for AH:
c = C S/(S + CA) (69)
From the contraction coefficient the entrance-loss coefficient of a long

culvert can be derived by using Eq. 68:
= (C' - 1)2 (70)
From the C values presented in Figure 61, the Ke values of a culvert

become:
C 0.56 0.58 0.61 0.69 0.83

c
K 0.62 0.52 0.41 0.20 0.04

e
The discharge coefficients of Figure 61 clearly show that the discharge

coefficient C will increase when the culvert length increases. This is
caused by the spreading of the jet into the culvert (two-step flow
expansion) resulting in outlet loss reduction. The small value of C for a
sharp-edged orifice is in agreement with other experiences.
100
Figure 62 is an example showing an inlet of a sloped culvert which can

remain partially filled with air or become fully submered, only depending
on the shape of the entrance and H/D. In the case of the sharp-edged
entrance there is a control section formed at the inlet and due to the
sharp edge the flow section is largely reduced. These elements affect the
discharge capacity much. The line for the full culvert is only an example,
the real discharge curve depends on the culvert slope, the culvert length
and maybe the presence of air caused by vortices.
Figure 63 shows the sharp-edged box culvert inlet with a discharge

behaviour similar to the sharp-edged coridition of the circular culvert of
Figure 62.
The discharge formula which is used in Figure 63 is:

Q = Cd w D I(2gh) (71)
where h is the height of the upstream water level above the bottom, D the
culvert height and w the culvert width. The expression (71) is similar to
the one that is often applied for gates (see Equation 64). The results are
comparable as well. (Figure 39 and 46).
Figures 64 to 66 show intake losses of concrete culverts. The section where

the gates are located is used as reference culvert section.
It can be seen that the nicely-shaped inlet of Figure 64, upper part
without slots, piers etc. has the lowest intake-loss coefficient, namely
0.16. For other conditions the loss coefficients vary between 0.2 and 0.8.
The losses in hydraulic models all seem to be lower than in prototype; an
explanation was not found in the reference [43]. Piers, slots, support
beams, etc. they all are elements that influence the hydraulic losses.

101
E9TRACE
SERIES 1 SERIES 2 SERIES 3 SERIES 4 SERIES 5 SERIES 6 SERIES 1
CONOITIONS
2.9
7-29 x
r[j
-
In

2
C=062
(''-058
-064 C '-065 C-072 C -'077 C'-060 c=061
EIfijtIcotenfrooce' -
2 =084 070 ' 8

A122 o122
NOT ES
11 r !5
C - 066 VOILICSof C gioco ore C - 072 C 0 131 C 023
- 0 61 avero es for the formuk,
VC 29
Loss coeftlderif
p 078
5 - CA'f2ROh
-
d(IoIIIoo to thi, lIg)

pItiIon of Corooi
073 C085
C =0 8RQ1J
CAIS
7 1,
_0 96 _ 0 92 C 0 go =088
C000CIISED 6Oi
- 0 Ii
UNIVE RsiTv or WIScONSIN
eUtLCTIN NO 26
gor. 309. Flow thogh .vbn,e,g.d tvb,. 2*2-0-2331.
Figure 61 Discharge coefficients of submerged rectangular orifices without

and with tubes, from USBR [38], and their derived contraction
coefficients C
c
item
haedge
5 1
Control of inlet-
part f culvert how
for entire range
square-edge inlet.3=0
4
h ,/D
3
Curved A - Square-
edged [niet
rouned
d Ln1et2Qj5
2 Ranie of k,m.d. wtara miz.d ori
itug flop tll occur
,ntrol of [niet - Part full culvert flow when--, 012
Curve B - Rounded inlet

0 LI
Q/g112 D 512
Figure 62 Head-discharge curves for square-edged and rounded inlets for

circular culverts on steep slopes, form USBR [38]
0.60
- U
(3
0 w=width of culvert
-c
0
1fl
0.50
c
Ii
0
9-
0
0 fl44
1.0 2.0 3.0 4.0
hID
Figure 63 Discharge coefficient for submerged box culvert spiliway with

square-edged top opening, dimensionless redesign after USBR [38]
103
AV ER AGE
CONDUIT PROPER INTAKE
LEPIGT Pl REYNOLDS VELOCITY COEFFICIENT
SNAPE PROJECT (i) DIAM (2) NUMBER (2) HEAD (1)
SWIGLE INTAKE (CONCRETE DAM CONDUITS)
PINE FLAT 54 29-3.6X10' 65-61

A
t ,8 (PROTOTYP[) (PROTOTP)
A90 830
C-90. E-O
- -c-- -
F -5LO, G 7
PROFE
L! $ ES 80 83 67X10' 97 007(3)
CI 20M00(L) (MOO(L)
-4: - - - a. 75,8.25
1
C100. E5 7
F.43.c-I4
P AN
DOUBLE INTAKE (EARTH DAM TUNNEL)
- 8 t)ENISON 40 2 X 10 66 019
(PROTOTYP[) (PROTOTVP)

A.250. 8-390
C. 190. 0.200
(- 90.T-530
C I'!0 0 SO SO 47 62-96X10' 61-62 012
Cl 25 MO0L) (M00(L)
PROFILE sci ABOVt)
FT RANOALL(5) 39 07-ISXIO' 6-72 0.25
(0T0T' (PROTOTYPE)
jD A-240,8- 160
C - 230. 0-220
PLAN E - ii o. 7-490
FT RANOALL 39 09-1.0 XI0' 4686 016
(1.25 1400(L) (1.400(L)

(st( A8Ov()
TRIPLE INTAKE (EARTH DAM TUNNEL)
TIONESTA 96 15-41X10' 7-50 033
0400(L)

- 7.5 T -660
D
PROFILE
IPITAKE HEAD LOSS
Va
VV(LQCITV IN C0PUT
PROPER
PLAN
D$4(NSIONS IN PROTOTYPE FE.(T

(OUIVALENT DIAMETER FOR NONCIRCULAR
SECTCNS BASED ON NYDRAULIC RADIUS
DOES NOT INCLIM C,AT(-SLOT LOSSES
L(NGTH Of TRANSITION
() ROOf ClvE MAJOR AXIS KORIZONTAL
Figure 64 Entrance losses of concrete conduits, form WES [43]

104
IFi1'11 PLAN
WAPPAPELL0
025 MODEL)
A31', 8=25'
C=201 , D20'
E= 10', T50'
T
K e (Tfr4PEE GATES)0 50
LOW K (GATES t &3)= 0.37
PROF ILE
11 2 11
T 10 N E STA
PLAN (136 MODEL)
A30, 8=22'
ct
L D
C= 16 1 , 0=19
E = 7.5, T66'
K (THREE GATES)=0.53
P,(GATES 1 &3)0.45
PROF! LE
8
qE9:= PLAN
ARKABUTLA*
(t 25 MODEL)
A= 10.5 , 8=16
C= 17 1 , 018.25
E= 8.5, T75'
K e (THREE GATES)=0.79
K 0 (GATES t &3)0.68
0
*CONDuIT NOT CIRCULAR
PROF LE
INTAKE HEAD LOSS I,_
v2
h -
e2 9
WHERE
LOSS COEFICIENT
V =GATE PASSAGE VELOCITY
g = ACCELERATON, GRAVITATIONAL
NOTE TO CONVERT LOSS COEFFICIENTS

INTO TERMS OF CONDUIT VELOCITY
HEAD MULTIPLY K. BV THE SQUARE OF
THE RATIO OF THE CONDUIT AREA TO
THE GATE PASSAGE F10W AREA t(CXE)
OF THE GATES CONCERNED.
Figure 65 Entrance losses of three gate conduits, from WES [43]

105
PLAN
EAST BRANCH
(1:25 MODEL)
A=23', 8=8'
T 23
:L FLO
H
PROF LE
FORT RANDALL
SEE CHART 221-1 FOR
Dl.4ENSIONS OF FORT (PROTOTYPE)
RANDALL INTAKE
K(TWO GATES)0 44
K0 (ONE GATE)=0.78
]jjM:::
EI
SARDIS
(1 25 MODEL)
PLAN
A= 6 8=235
- C= 12 1 ,D=17'
1/ H T
jR;
PROF LE
INTAPE HEAD LOSS h

v2
I = K-.-
WHERE
K.ZLOSS COEFFICIENT
V ZGATE PASSAGE VELOCITY
g = ACCELERATION, GRAVITATIONAL
NOTE TO CONVERT LOSS COEFFICIENTS

INTO TERMSOF CONOUIT VELOCIT'
HEAD, MULTIPLY KBY THE SQUARE OF
THE PATIO OF THE CONDUIT AREA TO
THE GATE PASSAGE FLOW AREA t(CE)
OF THE GATES CONCERNED
Figure 66 Entrance losses of two and four gate passages, from WES [43]
4.2 Friction losses
In Figure 67 the Moody diagram is presented. This diagram shows the wall
friction for pipe flow. The results can also be applied for free-surface
flow if the pipe diameter D is replaced by the hydraulic diameter Dh or the
hydraulic radius Rh.
D=Dh=4Rh (72)
Rh is the hydraulic radius which is defined as the ratio flow-section

wetted-perimeter.
The Moody diagram presents the friction coefficient related to the fully
developed boundary layer where the velocity profile harmonize with the pipe
diameter and the wall roughness.
The C f coefficient is related to the wall friction which equals the

boundary shear stress t in the flow:
t = C f (0.5 p 2 ) ( 73)
(2 represents the averaged velocity in the pipe).
The Reynoldsnumber TJ D/v is related to the average pipe velocity and the
pipe diameter (v represents the kinematic viscosity, with a magnitude
of about 10_6). For prototype structures the Reynolds numbers are so large
and the wails so rough that only the horizontal line parts of Figure 67 are
relevant. The friction formula for large Reynolds numbers is approximately:

Cf = 0.0606 [log (3 (74)
For instance, in Thijsse [36] the following k values are found:

old concrete k = 10-20 mm rusty steel k = 1-2 mm
plaster k = 0.5 mm riveted steel k = 0.5-2 mm
smoothed concrete k = 0.2 mm welded steel k = 0.1 min
rock excavation k = 0.2-2 m
sand in movement k = 10-100 mm stone revetment k = 5-20 mm

slopes (+ vegetation) k = 0.1-0.2 m riprap k = D90
107
Figure 68 illustrates how at the entrance of a rough pipe the local friction
developes with respect to the distance. Cfx is the friction coefficient
averaged from the mouth to the distance x. The higher C f at the mouth, is due
to the greater velocity gradients near the wall (at a contraction the velocity
profile tends to be block-shaped, without a boundary layer).
More general information about how friction developes with respect to the
distance can be found in Figure 69, from Harrison [15].
The friction coefficient has again been expressed in terms of a distance-

averaged value. Such friction tests are performed with a plate in an unbounded
flow. In practice, friction in a pipe or in a culvert will not fali below the
friction belonging to the fully developed flow in the culvert section over
long distances. The latter can be derived from the Moody diagram. Harrison
[15] shows that in the transition region between laminar and turbulent flow
there can exist lower C f values than presented in Figure 69. This is similar
to what is indicated in the Moody diagram of Figure 67 in the transition
region. Related velocity profiles are, among others, found in Kolkman [20].
Pugh [26] presents measurements in a rectangular prototype conduit. The

friction that was measured at the pipe inlet was greater than Tullis found
[ref. 371. Shear stress became constant at a distance x/D = 20 which is in the
range of what Tullis found. The other test results of Pugh confirmed what is
already known about velocity profiles and boundary layer development.
B. The developed flow
2
0
0.05 Ct)
1.5 z
0
0
0.02
0
z to 0.01
0
8
0 0.005
6 0.002
5 IkIi - 0.001
4 0.0005
0,0002
3 0.0001
2.10 e
10
B i6 10 7 io8
- ReD = Dh/V
Figure 67 The Moody diagram
108
Cfx FCfO
2.6 MEMO
BEEN
2.2 AIMMIN
31MIMEMMIMME
1.8 NMIMMMMIEMM
OMIMMMMIMMM
1.4 'a
MWNOMMMIMMM
1.0
IMMELEFFMOM
0 4 8 12 16 20
xlD
Fig. 68 Distance-averaged frjction factor at a pipe entrance from

Tullis/Wang [37]
0.016
xIk200
0.012
400
t-)
c
0 ._...S,. ROUGH/ TURBULE NT 103
0.008 04
L
c 10
- -
4
0.004
0-
10 10 6 10 7 10 8
Re: U 0 X/)
Figure 69 Distance-averaged friction factor at a plate (Harrison [151)

109
4.3 Slots
Loss coefficients of slots are presented as fractions of the velocity head

of the average flow velocity in the culvert section where the slot is
located.

LH = 1< 2/2g (75)
In general, the slot coefficients are small compared to other losses. A

possible additional effect of a slot is, when located just upstream of a
diverging section, that the flow spreads better over the divergence owing
to the induced turbulence. Then there is a greater recovery of potential
energy and reduction of energy loss.
Levin [22] tested a number of slot shapes located in a round tube (diameter
210 mm) along the whole circumference and in the walls of a squared
section, having a height and width of 150 mm. Figure 70 shows the
investigated slot sections.
The results of the test in the squared tunnel were not accurate because of
the deformation of the test section caused by hydrodynamic pressures.
Conclusions were:
Between Reynolds numbers of 3.5 10 and 8 10 (related to the tunnel
section) no effect of the Reynolds number was found
Low pressures locally appearing downstream of the slot faded out at a
distance of 12 diameters behind the slot
Rounded corners gave a decrease in loss in the order of 7-15%
The tests in a circular tube with the slot all around (only the slots with
sharp edges were tested again) show that the influence of the Reynolds
number was still observed till Re = 106 ; this is amazingly high. The losses
are shown in Figure 71.
110
v
0.050
0.045
31.5 form
J
[i1zes
mrn
0.040
JT1i JTL :.:::

0.035
{f9L j1 0.020
0.015
0.010
10 to Re

Figure 70 Slot sections of the Figure 71 Test resuits with the

tests by Levin [22] circular tube section
(D = 209.5 mm)
A theoretically plausible approximation is that losses are proportional to

the length (along the perimeter of the culvert) and inversed proportional
to the culvert section.
The further elaboration performed by Levin is not presented here; the
empirical formula are complex and slot losses as he found them are low
compared to other losses and their inaccuracies.
111
4.4 Exit losses and hydraulic jump relations
Exit losses are defined here as losses which occur at the transition of an
outlet into a downstream canal or a larger basin like a lake or the sea. In
certain cases the control section of the structure can alos be at the exit.
The exit can be an abrupt expansion or a gradual expansion with divergent
walis and/or a sloping bottom.
A reduction in exit losses gives an increase in discharge in the case of

submerged flow or intermediate flow (for definitions see Chapter 1).
However it has no effect on the discharge if the downstream water level is
so low that modular flow is involved. On the other hand, the modular
(maximum) flow will still occur at a higher downstream water-level. So it
is meaningful to calculate exit losses in case of a fully submerged flow
(at low velocities), and of intermediate flow till the limit of modular
flow.
The difference between submerged and intermediate flow causes that at
intermediate flow the exit losses with abrupt expansions should always be
calculated with the momentum equation whereas at fully submerged flow the
Carnot equation can be applied as well.
When no special measures are taken, the exit loss coefficient at an
(abrupt) large expansion is equal to 1, which means that all kinetic energy
is lost.
1f the outflow runs into a canal that having a limited section, it foliows
from the momentum equation that some of this energy is regained.
From the Carnot equation can be find for instance:
AH = (V1 - V 2 )2/2 g
This shows that a greater V 2 yields smaller losses.

Expansions that diverge gradually also regain potential energy, leading to
a further reduction of losses.
At modular flow conditions in the control section, super-critical flow will
occur downstream of this section, where energy loss takes place by bottom
friction and/or a hydraulic jump.
It has been decided to include in this section also the hydraulic jump
relations. They are important when studying the limit-conditions downstream,
where modular flow can still occur.
112
A. Abrupt expansions downstream of the control section
It has been discussed in Chapter 1, that at an abrupt (900) expansion the

momentum equation can be applied and losses can be calculated because then
the hydrodynamic forces, exerted on the back face of a weir, a gate or an
orifice are known to be hydrostatic. So the forces exerted on the down-
stream water by the rear side of the outlet are related to the local water
level. At the inflow of the downstream water, the pressures and the
advective momentum flux (fluid density x discharge x velocity pQV) are
often also known, being the pressure and momentum advection which belongs
to the modular flow condition of the control section.
1f there is a nappe over a shrap-crested weir, the highest downstream water
level at which modular flow in the control section can exist, is the water
level equal to the crest level. Together with the nappe momentum, the down-
water level h 3 at some distance further downstream can be determined. In
Figure 72 this is illustrated for a nappe with width w, flowing into a
downstream basin with width W. The coefficients m 1 and m 2 refer to
discharge coefficients. The aim of Figure 72 is to calculate h 3 when a
certain H is introduced. Solving these equation is laborious but the
procedure of Chapter 2 can be helpful.
A similar computation can be set up for a submerged hydraulic jump. The
discharge relation of flow under a gate can be treated sirnilarly to Q in
Figure 72.

113
Opening width w
Flume width W
H
07
h.2 C 2*~ d :
1 2
=IuI
1 '.-I I[-
Assumpt long:
1 atinoapheric pressure in upper part of II (side ventilation)
2 hydrostatic pressure in lower part of II
3 weli-distributed flow distribution at section 3
Q1 = m1w J 2 1fg(}1 1 - h 2 )) dz = m 1 w h 2 fg (H 1 - h2 )
(H -h )
2 =
= dz m 2 w (12g) (H 1 - h 2 )/ 2
m2w
0
Q = Q l (A)
11
2 =
Moentum-f1ux 1 = pm 1 w 2g (11 ) dz pm 1 w 2g (H 1 - h ) h2
1 h2 2
(11 h
1 )
-
Mornentum-flux 2 = pm2w 2 2g z dz = pm 2 w g 'I - h2 ) 2

0
Momentum-equation qpplied to the fluid between 1 and 3
t i+ Q 2 ) 2
pgW d + Momentum-fliix 1 Momentum-flux 2 = V?pgWd + p wd3 (B)

From I1 1 and Q (when m 1 and m 2 can be estimated) d 2 is found from A, and

from d 2 and Q, the value of d 3 js calculated with B (may be using the
program of Chapter 2).
Figure 72 Exit conditions with a submerged hydraulic jump caiculated with

the momenturn equation
114
B. Hydraulic jump relations (two-dimensional)
Only those hydraulic jump data which are related to the computation of
hydraulic losses are presented here. Not inciuded are considerations about
the jump stability, the design of stilling basins etc.
The application of the momentum equation for the hydraulic jump is well-
known, the results are presented in Figure 73.
Figure 73 is based on the momentum equation for a rectangular prismatic
canal:

'/ipgd 2 + pq(q/d 1 ) = pgd + pq (q/d 2 ) (76)
Division by pg and introducing the Froude number F for the upstream side
= gives
+ dF 2 = d 2 + ( d/d 2 ) F (77)
or:
(1 - d 2 /d 1 )(1 + d 2 /d 1 ) = 2 F(1 - d 2 /d 1 )/(d 2 /d 1 )
resulting in
d 2 /d 1 = (- 1 + 'v'l + 8Fi) (78)
This expression is found in many references. It is only valid at a well

distributed flow, because the momentum-flux does not contain a correctlon
factor.
In Figure 73 the following symbols are used. The index 1 refers to the
oncoming supercritical flow with depth d 1 , velocity V 1 , and index 2 to the
conditions behind the jump.
The term h refers to the difference in depth h. = (d 2 - d 1 ). The symbols
and H 2 refer to the energy levels.
Figure 73 can be well applied for the hydraulic jump at a horizontal bottom
when the starting conditions upstream of the jump are known. This is for
instance the case at free flow models of a gate.
115
For the design of stilling basins the author has designed a chart variant
from which directly the level of the bottom of the stilling basin can be
derived, see Figure 75. This is done by introducing the critical depth d
as the length by which all other length parameters are made dimensionless.
d 1 and d 2 again refer to upstream and related conjugate waterdepth. H 1 and
H 2 is respectively the upstream and downstream energy head.
The critical depth represents the discharge per unit of width:

d c = (q2/ g )1/3 (77)
The required energy loss AH of the hydraulic jump related to d determines

the Froude number. The bottom level of the stilling basin also follows
directly, because H 2 /d can be read in Figure 75. A second curve with 10%
extra reserve in conjugate depth (advised for the USBR stilling basins, see
Ven Te Chow [40]) is also included.
The use of the graphs, Figures 73 and 75, can be replaced by the
computation procedure of Chapter 2 (introducing directly the momentum
equation). In that case the momentum-flux can contain correction
coefficients for the uneven velocity distribution.
Particularly the momentum flux of the incoming flow will be sensitive to
values greater than 1. (see further Chapter 1, Figure 1.3).
116
1
1.0
CY Lr 0.9 0.1
0
0
0.8 0.2
1Hfl
0.7 0.3
mi
0.6 0.4
_ uSBR
IN
0.5e
IHIi
, 0.4 0.6
_
vII_ iIH!i
0.3 0.7
>
1110
0.2 0.8
0.1 VA 0.9
jump 1.0
n
2 3 4 5 6 7 8121620
Volues of F1
Figure 73 Relations for the hydraulic jump, from Ven Te Chow [40].
L
d2
II II
UndIa Iweok Oscillati fl9
Steody Jump Strong Jump
''N. jorop jump
s o rIot/ WOVY Best performance - Acceptob ExpensivesHlingbonond
orbojence ont performance rough sorfoce condilions
1 1 1 II .L
O 4 2 3 4 5 6 7 8 9 40 11 12 13 14 15 16 17 Id 1 O
Figure 74 The length of a hydraulic jump, from [40].

117
13 19
Hydraulic Jump Relations
12 18
r (dl/dcry2= /
4/dcr
- ---
11
- F 2" - F
4/3 17
Cjd cr =(d 1/d cr )(V 2 +1L1)

10 18
(10d 1/d cr)
-------/ ---/- - 15
14
7 13
/
Hi/d cr
(continued)
---
I 12
(continued)
5
4 Nlj 107 reser e in h2 H0/h er

V\
j
7-
el.
1
0
0 2 4 6 8 10 12 14
Index 1 refers to infiovring (supercritical) flow

Index 2 refers to outfiowing (subcritical) flow
Figure 75 Revised hydraulic jump relations for stilling basin design

118
C. Hydraulic jump on a slope
lof,
NEMN
OEM
M an
ORI
5uN1
Iui.iiiuiiiiiiiiiiiiiiiiii;i;;;;i;;J
v1//,
t
0
Figure 76 The hydraulic jump data sloping bottom (from Ven Te Chow [40])
For a sloping bottom the two-dimensional hydraulic jump relations are also
known, see Figure 76. It is important to know at which downstream water
level the modular flow gets (partially) submerged if the upstream jet is
the outfiow of a gate or a weir nappe . With Figure 76 this can only be
found with an iterative search. Therefore another graph has been prepared
in Figure 77, using the resuits of Figure 76. 1f tH, and d cr are known, the
dip in the water level can directly be read. AH the difference between the
energy level of the incoming flow and the downstream water level, and the
-
// ' - ---
119
dip is defined as h. 1f the dip is known the Froude number v 1 // d 1 can be

found, see the expression mentioned in Figure 77, and from that the length
of the jump can be read in Figure 76.
-
-
_.
h /d c
4
S:0
energy level
05
3 --
S 0 : slope
-- ..-
----
.10
S ..-;
1
2
5.3 :0
5:0.5
v 1i vga1
S c :O.2
0.25 F1 =[2AH+Ah
-
ir. 20 30 L dc
A H/dc
Figure 77 The dip in the water level at the hydraulic jump on a slope.
120
D. Gradual exoansions
For submerged flow at very low Froude number, losses in gradual expansions
are known from internal flow data.
The expression to define the losses is related to the Carnot equation at
abrupt expansions:

AH = F (V 1 - V 2 ) 2 /2g (78)
Free-surface flow tests were carried out by Formica [13]; Bos [7] compared
these resuits with Idelciks resuits for internal flow with various top
angles (see Figure 78). Other types of expansion investigated by Formica
are shown in Figure 79.
The scale of the Formica tests was quite small and presumably the Reynolds
number will stili have had influence. This, however, cannot be the explana-
tion of these losses which are twice as low as those of Idelcik. Idelcik
suggested that the expansion ratio is of little relevance. This is not a
priori evident. The curves of Idel'cik are for a = 90 a little higher than
the Carnot losses which were based on a hydrostatic pressure distribution
just after the expansion. This agrees with other references. So for this
moment it seems advisable to multiply all formica coefficients by 1.2.
In Figure 80 Idelcik data are presented for a gradual expansion with
outflow in a large canal, using an optimal diffusor angle (low losses).
VAL 1JES OF 4 itolO 1108 11o6 1105 004 1r3 002 1101

1.2
1.1
1.0
all vcllues
.9
.8
.7
.6
.5
4
-
two dImenJaoI
/
t
FORMICA
F1 0 79
W1/W2 :0.57
.3
.2
.1
0
2 3 4 5 9 7 8 9 10 20 30 40 50 90 70 80 90 100 180
). o'_
Figure 78 Expansion losses in closed and open conduit flow, from Bos [7]
and after Formica [13] and Idel'cik [17]
121
v
027
---
4
--- 07 0.29
06 0.45
- -
044
3 4
(v-v
VcIus of E in the qquation oHz
29
Figure 79 Losses at expansions of different shapes at free surface flow at

an expansion rate of 335 mm/205 mm, (taken from Bos [7] after
Formica [13])
122
2- dinierSiOfl0l ditlusor
d. opt
Loss of energy head: Constant depth:
H 011 = (V/2g) n a 2 /a 0 = A 2 /A0
n 1 a 1 /a 0 A 1 /A Q
Optimal diffusor angle a opt - when L/a 0 is given:
L /Oo
0.3 1,0 2.0 3.0 4,0 3,0 16,018.0 10 12 14
1.5 25 23
8 1 8,0 6.4 3.4 4.1 3,5 2.8 2.4 2.0 in degrees
2.0 33 IS 2 9,1 8,4 1,3 6.0 5.2 4,7 4.3
a
opt
2 S 37 26 8 14 12 10 9,4 8.0 1,0 6,3 5.6
3,0 39 27 20 16 13 12 II 9.1 8.0 1.2 6.4
4,0 42 30 21 II 5 3 12 10 9,0 8,2 7.4
6,0 45 II 23 8 16 14 II II 10 9,4 8,3 n 1 1 + (2 L/a)tg 0
opt
8.047 32 23 19 II 5 14 2 II 1091
10 48 33 24 20 II IS 14 12 II 10 95
14 49 34 25 20 II 16 14 13 12 II 9.9
20 50 35 25 21 18 16 IS 3 II II 10
L/Qo
0,5 1.0 2,0 3.0 4.0 59 6,0 8,0 10 12 14
1.5 004 004 0,04 0,04 0.05 0.06 0,06 0.08 0,10 0.11 0.13
2,0 012 0,09 0.01 0,01 0,06 0,01 0.07 0.02 008 0,10 0.12
2,3 0,18 0.14 0,11 0.10 0.09 0.09 0,09 0.09 0,09 0.10 0.10
30 0,23 0,18 0.14 0,12 0,11 0 II 0.10 0.10 0,10 0,I( 0.1 t Basic loss coefficient
4,0 0.30 0,24 0.19 0,16 0,13 0.14 0.13 0,12 0.12 0,12 0.12
60 0 1 38 031 0.25 0,21 0 19 0.18 0.17 0,16 0.13 0.14 0.14
80 043 036 0,28 0.25 0,22 020 0,19 0,11 0.16 0,16 013
00460380300260,240220,210.190,180.17016
4 050 041 033 0.29 016 014 0.22 020 019 0.18 0.14
20 0,33 0,44 0,35 0.31 0 28 0.25 0.24 0.22 0.20 0,19 0,19
0.4
+ r)
0.2
8 I
0 0.2 0.4 0.6 0.8 1.0
Figure 80 Determination of optimal diffusor outlet shape at constant depth
condition, from Idelcik [17]

123
E. Hydraulic jump in a gradual expansion
Downstream of a control section often divergent sections are applied. These

have two advantages:
- lower losses, particularly at the submerged flow conditions
- spreading of discharge over a larger width
Systematic test results are only available about conditions with a

horizontal bottom.
Two conditions can be considered:

- Free hydraulic jump, systematically investigated by Woodbury and
Padmanabhan, ref. [44]
- Submerged hydraulic jump (behind a gate with non-modular flow),
investigated by Abdel-Gawad and Carquodale, ref. [1].
The hydraulic jump in a divergent section is called "radial hydraulic

jump". The width at the throat w and the divergence angle 20 determine the
radius r 1 before the jump (being w/sin (e/2)) and the phenomena in the
divergent part are supposed to be a part of a circular hydraulc jump
condition.
For the free radial hydraulic jump the sequent depth can be derived from
Figure 81. There is still a slight dependence on the ratio of the radius r 2
and r 1 , so only a first estimate is obtained.
For the submerged hydraulic jump (definition scetch is Figure 84) the
relations found by Abdel-Gawed and Corquodal [1] are presented in Figure 85
and 86.
124
Nototion
r EGL-- jAE V/2g
V/2g 'y= d2 / d 1
FLOW
d = V,
Npid2
r0 = r,/r 1
r 2 -r %
1 2
P = pressure
F side pressure etc.
Figure 81 Definition scetch free radial jump

14
r0
1.20 - 1.40
12
o 1.41 - 1.60 / o
1.61 - 1.80
10 1.81 - 1.85
, 0
Eq. (75)
0-
east Sq.
8
.s
YO
0
1 2 3 4 5 6 7 S 9 10
Figure 82 Radial jump sequent depth ratio as function of approaching

Froude number and radius ratio, from Woodbury and
Padnahabnan [44]
125
6

S -
di 31
4 ;...
r2 - r 1 T
/ ..S
d2
3
ro
2 I2O_l4O
1 41 - 60
161180
181 - 185
0
1 2 3 4 5 6 7 8 9 10
F1
Figure 83 Free radial jump length [44]
H =energie heght
.7
.6
e
1. 0
- x
70
.5
AH
.3
.2
1
1 2 3 1 5 6 7 8 9 10
Figure 84 Free radial jump relative energy loss, also [44]

126
notat 10fl
4. :d/d
y0 = d2 / d1
P2 ro = R2/R0
1i
D0d1 F1 dl
Vg--
LR
r- DIr.ctlon
RO
Figure 85 Definition sketch submerged radial jump
iWel
model, S, radial
= free jump)
0.9 - (F
(S = submerged jump)
- r=1.8 7L:II_i.
0.8
0.7
0.6
T
-
0.5
--
0.4
/I'', /
0.3
/
0.2
/
____ H =energy head
0
1 2 3 4 5 6 7 8 9 10
--
Figure 86 Relative energy loss of the submerged radial jump, after

Abdel-Gawad and Corquodale [1]
127
0
XPRIMNTA1. DATA
2<W <4
70
4< y <6
6< Y <8
60 --F.RAOIAL --
- TH0RY
50 -
.8
40
30
10
0
2 3 4 5 6 7 e 9 10
Figure 87 Jump length of the submerged radial jump [1]

128
4.5 Trash racks
The losses at trash racks depend 0fl:
shape of the bar section,

blockage factor or solidity (S) which can also be written as Zw/W (w=
width of one single bar, W = width of the whole flow section if bars are
absent),
the angle of approach of the flow. 1f there is horizontal flow and the
bars stand inclined, the sin(e) is introduced, with e = 900 when the
bars are vertical.
when the bars stand vertical and the flow is horizontal but not
perpendicular to the plane of the rack, the skewness angle a is
introduced with 0 = 0 if the flow is perpendicular,
the other structure elements such as trusses, silis, slots which are
related to the rack.
Losses can be largely influenced by a skew approach angle, irregular flow

distribution, trash and ice.
1f Va is the approach velocity (or, when there are no bars, Va is the

average flow velocity in the trash rack section) the losses are expressed
in terms of:
AH = C (V 2 /2g) (79)
tr a
(this is method C of section 1.9 and Figure 8)

When the trash rack is inclined, the loss reduces to:
AH = C (V 2 /2g) sin
tr a
e (80)
Short bars (seen in flow direction) have higher losses than longer bars.
This is caused by the tendency of reattachment of the flow at the longer
bars.
1f the length is introduced as t, the increase in t is effective until
t/w > 5, see Figure 90. The effect of streamlining is important, see Figure
88, but the effect reduces when there is a skew approach flow. For further
data of flow resistance of different shape see also next section on bridge
piers.
129
Figure 89 shows C values for different sections of bars and various

tr
values of the solidity S.
Figure 90 contains information for greater blockage of flow. This
information is important because if the other trash-rack elements
(mentioned under point E) are taken into account the blockage of the flow
is much greater than from the bars only.
Trusses are sometimes badly shaped so that loss coefficients correspond to
bars with a small t value.
In section 1.9 the use of the Carnot formula for the computation of trash
rack losses has been discussed. Orsborn [26] also discusses this point,
referring to earlier work of Escande. For the contracted flow, a
contraction coefficient of .9 should be applied if S is in the order of 8%
and this coefficient is 0.65 if S is in the order of 50% (the contraction
coefficient is related to the net area between the bars). The contraction
coefficient is rather low and the important flow-area reduction can cause
that at an increased discharge (at free surface flow), modular flow
conditions occur causing a further decrease of discharge capacity.
In Figure 91 the losses are presented at skew approach velocities.

Allthough an empirical formula, constituted by multiplication of two
independent factors has its inherent inaccuracies, the losses correspond
within 20% with the other presented data. But in view of other inaccuracies
in the estimation of trashrack losses this is still acceptable.
INFLUENCE OF APPROACH-FLOW ANGLE ON BAR-LOSS COEFFIIENT

0
Conciitions: So1idiy = 0.37; Vertical angle e=90
Bar
Shape
Flow *
Bar Loss Coefficient, Ctr
Angles
a= 00 1.13 0.86 0.78 0.48 0.42 035
cx=30 1.46 0.76 0.71 0.43 0.68 0.22
cx=45 2.05 1.29 1.29 0.94 1.29 0.67
a=60 4.26 2.45 2.81 2.19 3.05 1.84
Ctr is loss coeficient A H=Ct r V 2/2g
Figure 88 Trash-rack losses at skew approach flow, from Orsborn [26].

130
w t'w=50
10
DEF1NIT1ON SKETCH -
LEGEND f
/
FELLENIUS TESTS (STOCHOLtl)
RC4 TESTS (MuNc 14)
L_
Li J /
- /
O _ -
1.0
C1
i70/J
THE NUM8ERS SHOWN
APEVALUESOFtw
/ (
>
1
0 4J SLIGHTLVROUNOEO
CORNERS
4
/
o
EEEEEE ?
It
10.9
t
I
50
O
00 01 03 04
Ar
NOTE HEAD LOSS THROUGH RACK IN FT

= VELOCITY AT SEC11ON
*1THOUT RACK N mis
Ct r HEAD LOSS COEFFICIENT
- AREA OF BARS
- AREA OF SECTION = SOLIDITY
OPEN CHANNEL FLOW
TRASH RACK LOSSES
Figure 89 Losses of different types of bars, after WES [43]

131
tr
measurement
point
20
10
measurement
points
5
Al1w
.,
r- -
2
meaurement '
points
1
.5 V J. 1
t/w
2
= 1
------ 3
- - 4 and greater
.2
0 .25 .5 .75 S = solidity
Figure 90 Loss coefficients for rectangular bars at greater solidity

after Orsborn [26]
132
rock
Lw]Iii:
,. V
Va /
Ctr; AH/( V 2/2g) 32
i 99JftI.LI
w w w w w w w w w w
llJ1Li
1 2 3 4 05w
5 6 7 8
10
rvai
1.4 3
8
1.2 3.
in]
2.
_v
1.0
0.8 2. d
M
5 F-P..PA
-
2
0.6
0.4
t' _ _
0.2 0
.-1
S = blockoge factor = w/W
Figure 91 Trashrack-bar losses at skew approach velocities, after

Idelcik {171
133
4.6 Resistance of bridge piers
In the following, only very schematised cases are considered. For example,
no deepening between the piers by the locally enlarged bedload capacity or
by scour and no rough bottom protection around the piers are taken into
account. The flow is assumed to be parallel to the canal or river axis.
Moreover, all experimental data were obtained by tests performed in
rectangular flumes and all data presented here relate to rectangular flow
sections. The ratio of water-depth over pier-width is not considered,
although there are indications that this parameter rnight have some
influence (Naudascher [251).
The above simplifications permit the application of one-dimensional
hydraulics, in the same way as discussed for other hydraulic structures in
section 1.4.
Three regimes of flow are distinguished (Figure 92); types A, B and C refer
to the WES classification [43].
A Subcritical flow condition (fully submerged and intermediate flow).
In this situation the water level upstream of the bridge depends on the
initial water depth (d.) and the discharge of the undisturbed river. The
initial water depth also corresponds to the depth downstream of the
bridge (d 2 in Figure 92) apart from a slight modification by friction
and head loss due to deceleration.
B Modular flow condition (critical flow).
A unique relation between the water level upstream of the bridge and the
discharge will occur here. It is usual to define the critical depth in
the control-section by using the full width of the flow section between
the bridge piers, neglecting the effect of side contraction of the flow.
The necessary correction coefficients are found from empirical data. The
energy level which belongs to the critical depth between the piers
determines the energy level upstream of the bridge and the corresponding
upstream water level can also be calculated. The water level downstream
of the bridge can be either supercritical or subcritical (see Figure
92). 1f the initial water level is in between these levels, the super-
critical water-level will occur in combination with a hydraulic jump
further downstream. 1f a downstream loss coefficient is estimated the
local downstream water-level can be derived, but this has no practical
importance.
134
Which condition - the supercritical or the subcritical one - will exist

is found from a further analysis of the flow downstream of the bridge.

w
C w
BRIDGE PEP
PLAN
CLASS A
d2
do> dc

d2 >d c
CLASS B
d1 > dcl
d2>dc
> dc d1 = dcl
d<d
d:dj T
d2<dcl i .
j
d,<dc d2<d
[LE VAT 10 N
Note: S = w/W = blockage ratio

w = total pier width
W = gross channel width
d = upstream depth
d 1 = depth within pier section
d 2 = downstream depth
d = critical depth within the unobstructed
channel section
d 1 = critical depth within the pier section
d. = initial depth
Figure 92 Water levels for the three flow regimes, A, B and C

135
C Super-critical flow condition.
At high river velocities the water level just upstream of the bridge as
it is found for type B flow is not sufficient to evoke a hydraulic jump
(upstream of the bridge) and the flow under the bridge remains thus also
supercritical. All water levels are now fully determined by the upstream
conditions and the upstream water level d 0 equals now the initial water
level d.. Regime C, which occurs at mountain rivers, is not of interest
for backwater effects; its interest is mainly in forces exerted by the
flow on the piers and the maximum water level under the bridge.
In table 4.6.1 the properties and needed calculations for each of the flow
regimes are summarized.
name per influence needed

definition initial depth calculations
A Subcritical d1 > d 1 d 2 = d. d = f (q, S, d 2 )
B Modular d 1 = d 1 no influence d = (q, S)

(critical)
C Supercritical d 1 < d 1 d 0 = d. d 2 = f (q, S, d 0 )
Table 4.6.1: Properties of flow regimes A, B and C
Figure 93, from WES [43], shows the limits of type A, B and C flow regimes
as function of the initial conditions and the blockage ratio.
Curve 1 is obtained by taking the critical depth in the section between the
piers and its corresponding energy head, from which the (equal) water
levels up- and downstream at the wider sectons are derived (subcriticai
condition without additional losses). The calculated water level applies
to the downstream limiting condition between modular flow and submerged
flow. Curve 2, now including hydraulic losses, follows from interpretation
of experiments; at the same critical flow condition between the piers the
downstream energy head and corresponding water level (which equals at the
limiting condition also the initial water level) is lower when losses are
involved. Curve 1 represents the subcritical depth upstream of the bridge
at modular flow. This depth is also the available sequent depth of a
hydraulic jump which can occur for supercritical flow in the upstream
region. From this depth the sequent supercritical flow conditions can be
derived, which results in curve 3.
136
In the following the three flow regimes will be considered in more

detail in the sequence C, B and A. For definitions see Figure 92. The
discharge per unit width, q, is defined in the initial situation without
bridge piers.
Flow type C (supercritical flow): No further experimental data are

available; WES advises to get insight in the flow behaviour by application
of the momentum equation. In Figure 94 a computation result for the water
level under the bridge is presented. The assumption was a constant energy
level without additional losses. However, this computation is on the unsafe
side so a factor a for flow contraction must be applied. It shows clearly
that an important set-up of the water level can be expected. The graph has
been obtained by application of the Bernoulli equation:
d 0 + q 2 /2g d 0 2 = d 1 + {q/(1-aS) 12/2g d 1 2 (81)
(d 0 and d 1 refer to the water depth upstream and between the piers)
In the figure d applies to critical depth at the section without the

bridge piers. Between the piers the discharge per unit width is higher than
upstream of the bridge

q 1 = q/(1-aS) (82)
Type B (modular flow): Figure 95 presents the upstream water level for two
nose shapes of a bridge pier. Contrary to subcritical flow where pier
losses are determined by the nose and the afterbody shape, only the con-
traction induced by the nose shape will determine the flow section at
modular flow.
Type A flow (fully-submerged and intermediate flow): the bridge pier losses
for low flow-velocities (due to the near-horizontal water level) are fully
comparable with the bar-losses at trashracks presented in the former
section. But at higher velocities a dip in the water level occurs, giving
an increase of losses. The higher velocity is expressed in terms of an
increased Froude number, a decreased relative water depth d/d or an
increased relative velocity head V 2 /2gd. (These numbers have in fact all
three the same significance, for instance F 2 = (did)3.)
For increasing velocity, the dip in the water level initially increases in
proportion to the velocity squared, but when the dip increases the velocity
137
in the contraction increases more rapidly due to the additional effect of a

decrease in water depth.
Two formula for bridge pier resistance are discussed below, the Rehbock
formula and the Yarnell formula (applied in the WES charts [43]). To have a
basis for comparison some uniform parameters are introduced in both
formula.
Rehbock (taken from Reh [29]) developed the following formula:

Cb V/2g (83)
C = [S - S(6 - 1fl (0.4 S + S 2 + 9 S 4 )(1 + F) (84)

br 0 0
= factor depending on the shape of the pier, S = blockage ratio = w/W,

F 2 = Froude number downstream = V 2 /I). The downstream Froude number is
also the Froude number for the initial condition. It is seen that Rehbock
took into account interaction between blockage and pier shape (the product
terms with 8 and S) but in his expression the amplification due to the
Froude numbers is the same for all bridge pier types. This indicates that
he tested the piers at relatively low velocities (near the modular flow
conditions it is mainly the flow contraction at the pier nose which is
important and not the total pier shape).
To enable comparison with other formulae for bridge pier losses Eqs. 83 and
84 are transformed.
Making use of the fact that the drag coefficient of one single bridge pier
equals (this will be discussed hereunder)
Cd = Cb/S (85)
the drag coefficient (derived from Equation 84) at low flow velocity
2 1 and zero blockage (S = 0) becomes:
F2
C do = 0.4 S 0 (86)

and Equation 84 transforms into
Cb= (Cd O _ CdS + 0.4 S)(S + 2.5 S 2 + 22.5 S)(1 + F 2 2 ) ( 87)

138
Remark:
Equation 85 is not formally derived here but when the drag of a pier is
defined as
D = Cd LL p V wd 0 (a)
and when the momentum equation is applied at a condition where there

are small velocities and a small head differences it can be found that
the pier drag is just compensated by the force induced by the
difference in water level
D = Ah (pg wd 0 ) (b)
And also can be found
Ah AH (c)
Now equation (b) in combination with (c) and (83) results, taking into
account that V 2 V o l in:
D = C (V 2 /2g) (pg wd0) = C b p V wd 0 (d)

br 0
Equation a and d combined result in eq. 85.
In Figure 96 the Cd of a number of piers as investigated by Rehbock is

presented. Although not all sizes are explicitly mentioned the drawing is
undistorted and all piers unless indicated otherwise, have a length of 6 to
6.7 times w (probably 6.66).
%.ES [43] proposes for Cb in eq. 83 the use of the Yarnell formula:
C br = (d 0 - d )/(V2 2 /2g) = 2K[K + 10(V 2 /2g d 2 ) - 0.6)(S + 15 S 4 ) ( 88)

2 2
The factor 1< can be transformed, using Eq. 85, into Cdo (V = 0 and S = 0)
C do = 2 1<2 - 1.2 1< (89)
from which:
1< = 0.3 + (0.09 + 0.5 C ) 0.5 (90)

do
139
The proposed Yarnell K factors are transformed in the next table into Cdo
and compared with the Rehbock values of corresponding shapes. Especially
for case a (the semicircular nose and tail), a great difference is seen.
Pier shape K derived Cd from Fig. 96
Semicircular nose and tail 0.90 0.54 0.84

Twin-cylinder piers with 0.95 0.67
connecting diaphragm
Twin-cylinder piers without 1.05 0.95
d iaphragm
90 deg. triangular nose and 1.05 0.95 1.15
tail
Square nose and tail 1.25 1.63 1.56
Table 4.6.2: The Yarnell K factors, the derived C factors therefrom and
do
corresponding values from Rehbock
The Yarnell formula is presented graphically in Figure 97. Figure 98 gives

a comparison with Rehbock results for shapes a, d and e (as indicated in
Table 4.6.1) for varying d 2 /d values.
The differences are probably due to the fact that the Rehbock formula is
here applied in the range beyond the experimental range from which It was
derived. Very small S values were probably not tested either, nor were the
very large d 2 /d values (to prevent too low values of the Reynolds number).
It is interesting that the best agreement is found for d 2 /d in the range
of 3 and for the square-nosed and -tailed piers.
In recent years an investigation of Al-Nassri and McBean [3] has been pu-
blished which covers a variety of shorter pier shapes. The range of F 2 was
between 0.25 and 0.6 the blockage varied between 12 and 35%. Figure 99
shows the test resuits. Instead of water depths the energy levels H with
respect to the bottom level were used.
A direct comparison with the Yarnell formula is difficult because the two
presentations are quite different; therefore only a comparison with the
Rehbock formula is presented. Shape 1 and 5 of Figure 96 was compared with
b and a of Figure 99; the comparison is presented in Figure 100. In general
the Rehbock formula gives higher water level differences (where the Al-
Nassri resistance is higher the diamond pier was also shorter and hence
less streamlined).
140
2.1
1.9
CLASS A
1.7 0:4
ASS A OR 8
1.5
t-)
1.3
t' 1
1.1
CLASS S
0.9
0.7
0.5
- --f::
Limitig
CLASS C
X by mom.rttum
CARSTANJEN)
metPtod dj =initial waterdepth
(without pier- s)
critical depth
without piers)
0.3
0.0 0.1 0.2 03 0.4
. blockage S:Zw/W
Figure 93 Flow regime as function of the initial flow conditions in a

rectangular canal, from WES [43]
141
HE
u
o
dcHticaI depth (unobstructed)
d 1 intIaI depth and d 0 d;
a > 1
pier
///// // ////////
0
0 05 10
d0 / d
Figure 94 Theoretical minimum change of the water level under the bridge
for flow type C (supercritical).
IIIIH
20
1.9
square nose pier
18
7 /
1.7
round riose piers

1.6 -
0
10
d: cr1 t. depth (unobstructed)
1.5
i '1IuiI d0LJE
1.4
t"
1.3
DEFINITION SKETCH
d 1 critdepth section 1
1.2
000 004 008 012 016 020 0.24 0.28 0.32
s= wIW
Figure 95 Upstream water depth for type B (modular) flow, from WES [43]

142
1 2 3 4 5 6
CL
C 156 1, 01 040 0,53 032

Do
7 8 9 10 11 12
CD O : 0,22 084 058 048 0,40 0,32
Figure 96 Drag coefficients of piers investigated by Rehbock after

Reh [29]
143
Ilulipil
'
0.12
D
oii IiIIw:IhIii 1IoiI IHl

-
cJ
L/ \
0
0 \
1
\
0.08 .I \
\
\

0.04
01)

1.5 2 2.5 3
- d/d
Round nose pters
At other shapes than round
nose the (d0 - d2 ) is a factor
A greciter
d0>d d2>d
DEFINITION SKETCN
Figure 97 Bridge pier losses for regime A (fully submerged and interme-
diate flow) based on the Yarnell equation, after WES [43]
144
\ \\\
VARNELL YARNELL
'S
S' \
_
u '-4
13 S0.3
z01 .oi
"4
005
0 0
\\
\\\\
s01 S0 1
\\\
.001 ______ _________ \\\ .001
REHCK REhbOCK
.0005
1 2 5 We 1 2 5 iL!J
d2 /dc d2 / dc
case a (tabie 4.62): case d:

s.mlrc
K nos. / toil 90 angle noe/taii
Rehbock Cr0.84 Rehbock C0 0 1.15
Yarneli K = (190 Yarriell K = 1,05
YARNELL
O.3
.o1
s'
"4 \
0
- a
.001
RENBOCK \\
.000 -
1 2 5
d 2 / dc
case e:
sqre nose / tOil
Retibock C 00 r1 56
YcirneU K = 1 25
Figure 98 Comparison between bridge pier losses computed from the Yarnell
and Rehbock equations

145
20,
LEGE1tD
:: &3 A B
S:0.333
S :0.278
o S :0.222
I
DIAMOMD - - 1
0 S)UIPE 5crrI 5c. F1ECTANGULAR
o S:0,166 z
9
118
8
7cm Hocto
---
--j
6
;3,:22;-PEAR SI)AP
7 cm
-_-___-( H-- 11cm
/
-
11=
c D
1 CIRCULAR
5 coij-._..
- .-
V2 /
1 L 1 t -- L. L I .t t J (equats F1 ( urtobstructed condltion))

01 2 0 4 .5 6.78910
r2 1 t S)
20 o --- --- --------- ----------
--
LEGENO fl79
rw -
S :0333
' 4 S :0.27e
:0,278 0 S :0.222
II
0 S :0,222 10
0
9 S :0,118
S :0,8
8 8
t / 1
6 /
105 (
6
..
0
.J0
-r S to
3 3
S
2 " 2 / ~.
1
01 2 3 4 5 6 7 8 9 10 01 2 3 4 .5 6 .1 6 9 10
F 2 0151 U2 F)I S)
20,
LE 0 EN t) - 1178
8 rIjK l
S :0,333 S :0.333
0 S :0.222
O S :0,278
o
o S :0.222
10 8 S :0,766
S:0.112
/
o S :0,166
..0 7
6 6
7: 5
:i '
0 0

/
o
2
II -L,-I-- L
1 j j I 01
F2 t(l-S) r2 (i . S)
Figure 99 Test resuits at regime A according to A1-Nassri and McBean [3]

146
20
20
III , . -
REHBOCP SHAPE 1 //j RHB0CK
-- rC00=053
SI-IAF'E S
1, /
t,
T..TJ /
10 * S S r0,3331
9
8 5=0,2781 II , 10
5=0333')
o S =0,222 > AL-NASSRI / 9
8 S =0 278
S AL-MASSPI
o S r0,166 rSHAPE B S o 5=0222 'SHAPE A
7 S =0,115) 7 o S= 0155
e S:011SJ
pier S
77//A per
'.5
1
- 4
1
1
III
0
13 ,o' 0
13
2
1/Je o
0,
1
//'.,I .. .
1/1
02 04 05 05 10 0.2 0.4 0.5 0 1.0
F21 (1-S) F2/ (1-S)
Figure 100 Comparison of pier resistance of the Rehbock formula and the
tests of A1-Nassri and McBean
147
REFERENCES
Abdel-Gawad, S.M. and Corquodale, J.A.Mc.,

"Analysis of the submerged radial hydraulic jump" in Canad. Journ. of
Civ. Eng. Vol 12, 1985 pp. 593-602
Abou-Seida, M.M. , Quraishi,A.A. ; "A flow equation for submerged

rectangular weirs", Proc.Inst.Civil Engineers, Vol.61, 1976.
A1-Nassri, S.A. and McBean, E.A. : "Bridge pier sharpe and contraction
ratios". Journal of Engineerign and Applied Sciences, Vol. 2, pp. 115-
129, 1983
Billor, J., Jaoui, A., Kolkman, P.A., Radu, M.T., de Vries, A.H.
"Hydraulic investigation for the diversion, the morning glory spillway
and the bottom outlet of the M'Dez dam in Marocco" (in French
language"). 13th Congress on large Dams, Q50, A62 ICOLD, 1979.
Boiten, W.; "Ice flap of the Borgharen barrage. Determination of

discharge coefficients." (original in Dutch), DELFT HYDRAULICS Rep.
M1318, 1980.
Boiten, W. ; "Vertical underflow gates discharge characteristics",

Research Report No 93. Agriculatural University of Wageningen, dept,
of Hydraulics and Catchment Hydrology, November 1989.
Bos, M.G. (ed); "Discharge measurement structures" pubi. DELFT

HYDRAULICS nr. 161 and Int. Inst. for Land Reclamation (ILRI),
Wageningen 1976.
Breadly, J.N.; "Rating curves for flow over Drum Gates", Transactions
ASCE Vol. 119, pp. 403-420, 1954.
Cozzo, G. ; "A formula to calculate the discharge coefficient of gates"

(original in Italian); lEnergia Elettrica nr 11-12, pp. 504-513,
1978.
DELFT HYDRAU1ICS " Discharge coefficients at overflow" (original in

Dutch), report M647, 1959.
148
REFERENCES (continued)
DELFT HYDRAULICS: "The rectangular broad-crested weir", report S170-

VIII, 1981.
DELFT HYDRAULICS: "The trapezoidal profile broad-crested weir",

report S170-XI, 1983.
Formica, G.,
"Preliminary experiments of head losses in canals with divergent or
convergent section" (in Italian).
Lenergia eletrica, Milano, Vol. 32 No 7, pp. 554-568, July 1955.
Gieseke, J. ;"Computation of discharge coefficient for Hower Bunger

valve" (original in German); "Die Wasserwirtschaft", 56, no. 10, pp.
315-323, 1966.
Harrison, A.J.M. , "Boundary layer displacement thickness on flat

plates", Proc. ASCE, Journal Hydraulic div. HY4, July 1967 paper 5339
pp. 79-91
Hay, N and Taylor, G.; "Performance and design of labyrinth weirs",

ASCE, Journ of Hydr. Div HYll, pp. 2337-2357, Nov. 1970.
Idel'cik, I.E.; "Momento des pertes de charge" Publ. Eyrolles, France

1969 (translated from Russian)
ISO, International Organization for Standardization: "Measurement of

liquid flow in open channels"; ISO Standards Handbook 16, Geneva 1983.
Kinsvater, C.E. and Carter, R.W.C.; "Discharge characteristics of

rectangular thin-plate weirs" Journ. Hydr.Div ASCE vol 83, HY6, paper
1453, Dec.1957.
Kolkxnan P.A. . ; "Considerations about the accuracy of discharge

relations of hydraulic structures and the use of scale models for
their calibration" IAHR symp. on "Modelling hydraulic structures"
Esslingen, Sept 1984.
REVERENCES (continued)
Kooman, D. and Korthof, R.M.; "Sili beam loads due to flow and waves";
In "Hydraulic Aspects of Coastal Structures; Developments in Hydraulic
Engineering related to the design of the Oosterschelde Storm Surge
Barrier in the Netherlands", part 2, pp 47-66, Delft University Press,
1980.
Levin, L.; "Theoretical and experimental study of fiow in slots"

(original in French); in "Le genie civil" No. 5, pp. 1-12, May 1971.
Nago, H. ; "Influence of Gate-Shapes on Discharge Coefficients";

Transactions Jap. Soc. of Civ. Engrs,Vol. 10, pp 116-119, Nov 1978.
Naudascher, E. ;"Hydraulics in open channels and hydraulic structures"

(original in German); Pubi. Springer, Vienna 1987.
Naudascher, E.; "Forces on hydraulic structures, Manual about

hydraulic structures design information"; Int. Ass. for Hydraulic
Research IAHR, publ. Balkema 1989 (in preparation).
Orsborn, J.F.,
"Rectangular-bar trashrack and baffie head losses", Proc. ASCE, journ.
Power Div., paper 6223, pp. 111-123, Nov. 1968.
Pickering, G.A.; "Lock culvert valve loss coefficients" Hydr.Lab IJS

Army Engin. Waterways Experiment Station WES. Sept 1981.
Pugh, C.A.,
"Turbulent flow in rectangular outlet conduit, Rend Lake dam, Big
Muddy river, Illinois"; Rep. H75-2, Hydr. Lab. of IJS Army Waterw. Exp.
Station WES, Jan. 1975.
Reh, U; "Simplified application of the Rehbock backwater formula for

bridges" (original in German). "Die Wasserwirtschaft", 1957 - 1958
pp. 240-242.
Rehbock, TH,; "The constantness of the discharge relation at sharp

edged weirs" (original in German), in "Der Bau-ingenieur" Vol.11,
No.48 pp. 821-826, 1930.
150
REFERENCES (cont inued)
Rouse, H.; "Engineering Hydraulics", publ. John Wiley 1950.
Rouv, G. and Inlekofer, H.; "Discharge over straight weirs with semi-
cylindrical crest" (original in German); "der Bauingenieur 49,
pp. 250-256, pubi. Springer, 1974.
Rouve, G. and Khader, M.H.A.; ,, Transition from a conduit to free

surface flow" IAHR Journal of Hydraulic Research, Vol. 7 No.3, 1969.
Sarginson, E.J.; "The influence of surface tension on weir flow"; IAHR

Journal of Hydraulic Research, Vol. 10, no.4, pp. 431-446, 1972.
Schnfeld, J.C. ; "Tidal computations by integration along

characteristics" (original in Dutch language); Note 55-14 of Central
Study Office, Rijkswaterstaat, 1935.
Thijsse, J.Th. ; "Hydraulica, in The Technical Vademecum" (originally

in Dutch the "technische vraagbaak") part W, pubi. Kluwer, the
Netherlands, 1951.
Tullis, J.P. and Wang, J.S.; "Turbulent flow in the entry region of a
rough pipe", Transactions ASME, pp. 62-68, March 1974.
USBR; "Design of small dams", US dept. of Interior, Bureau of

Reclamation USBR, revised edition 1977.
Varshney, R.S. and Mohanty, S.K.; "Discharge relations for submerged

weirs" Indian journ. of Power & River valley development, pp. 225-228,
July 1973.
Ven Te Chow; "Open Channel hydraulics", Publ. Mc Graw-Hill 1959.
Villemonte, J.R.; "Submerged weir discharge studies", Engin.News

Record, pp. 866-869, Dec. 25, 1947.
Wagner, W.E. ; "Determination of pressure-controled profiles of

Morning-Glory shaft spillways, "Symposium ASCE, transactions 1956.
151
REFEP.ENCES (continued)
WES; "Hydraulic Design Criteria", US Army Waterways Experiment

Station, Corps of Engineers, 17th issue 1977.
Woodbury, L.H. and Padmariabhan, G.; "Design implications of radial

hydraulic julnps" in "Civil Engineering for Practicing and Design
Engineers" Vol. 4, pp. 965-984, Pubi. Pergamon USA, 1983.
. location 'De Voorst'
S main office
main office
Rotterdamseweg 185
p.o. box 177
2600 MH Delft
The Netherlands
telephone (31) 15 - 56 93 53
telefax (31) 15 - 61 96 74
telex 38176 hydel-ni
location ' De Voorst'

Voorsterweg 28, Marknesse
p.o. box 152
8300 AD Emmeloord
The Netherlands
telephone (31) 5274 -29 22
telefax (31) 5274 - 35 73
telex 42290 hylvo-ni
North Sea
Amsterdam
Londen
Brussels .

Discharge Relations For Hydraulic Structures and Head Losses From Different Components

Uploaded by

Copyright:

Available Formats

Discharge Relations For Hydraulic Structures and Head Losses From Different Components

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Discharge Relations For Hydraulic Structures and Head Losses From Different Components

Uploaded by

Copyright:

Available Formats

QI

Discharge relations for hydraulic

Desk study for Rijkswaterstaat and

Q 953 delft hydraulics

GENERAL ASPECTS OF DISCHARGE RELATIONS AND LOSSES ..................1

1.1 Introduction ..................................................1

NUMERICAL PROCEDtffiE FOR COMPUTATION OF WATER LEVEL/DISCHARGE

2.1 Introduction ..................................................29

DISCHARGE RELATIONS FOR CONTROL SECTIONS ...........................36

3.1 The sharp-crested weir ........................................36

3.7 Labyrinth weir 69

4. LOSSES OF COMPONENTS UP- AND DOWNSTREAM OF THE CONTROL SECTION .....97

4.1 Orifice and intake losses .....................................97

a (gate) opening height (m)

A orrifice area (ni 2 )

Cb loss coefficient by bridge piers

Dh hydraulic diameter (4 x hydraulic radius) (in)

G0 height of gate opening (m)

K coefficient in the bridge-resistance formula of Yarnell

Lr length of radial JUMP (m)

Qg gate discharge (m 3 s' )

Q 0 discharge at modular flow (used as reference discharge) (m 3 s- 1 )

Rh hydraulic radius (m)

coordinate in vertical direction (m)

e divergence angle (= 2 a) (rad)

rate of submergence of submerged hydraulic jump (= d/d1)

1. GENERAL ASPECTS OF DISCHARGE RELATIONS AND LOSSES

Because internal flow is discussed in Part T of the IAHR manual, in the

The following type of equation is called a "discharge relation"

When in a structure one section controls either in total or mainly the

At a free-surface condition this equation becomes:

(w = culvert width, d = water depth in the culvert).

This relation differs markedly from Equation 1, and in equation 2b the

In this manual data are presented about discharge relations of control

Concerning the discharge characteristics of a control section with free

In Figure 1, a picture is shown of the generalized discharge/water-level

Figure 1 Generalized discharge/water-level relation for one geometry (and

In case the complete structure is considered, the characteristic property

Therefore, if the flow opening remains constant, the discharge increases

Figure 1 applies to one geometry only: at a structure with a gate such a

Resuits of computation of losses in free-surface flow at hydraulic

In Section 2 js introduced a computer algorithrn which can be used to handle

In one-dimensional far-field computational models (networks), and in two-

1f Ah < R, with R = f(h 0 ) then submerged flow: Q = f 2 (h 0 ) 1 2gAh

In dynamic networks or tidal computations one prefers to use discharge

The following type of expression might be suitable to represent the

is representative for the fully submerged flow condition (where Ah is

In this expression f 1 , f 2 and n can stili be functions of h 0 or Ah.

fully submerged regime, gets steeper in the intermediate flow regime. So a

Q = a + OAH + y AHb (with b > a > 0.5)

In two-dimensional far-field flow calculations the skew overfiow over dikes

Most computations of pressures and losses are based on one-dimensional

Figure 2 Correction coefficient a and 0 for a fully developed velocity

must be expressed correctly in the energy balance.

1.4 Some general remarks before starting the calculation of losses

In general an outlet work or a sluice must be approached as if they

There is a maximum downstream level (close to the control section)

Discharge-regulating structures are relatively short. The losses due to

Complicating factors in the analysis can be points F and G: