Bul81 - Spillway Shockwaves and Air Entrainement

This Bulletin has been prepared by Subcommittee No.
4
of the Committee on Hydraulics for Dams,
with the assistance from W. H. Hager (Switzerland)
Original text in English-French translation par Y. Le May
Ce Bulletin a ete prepare par le Sous-Comite n° 4

du Comite de l'Hydraulique des Barrages,
avec la collaboration de W H. Hager (Suisse)
Texte original en anglais - Traduction en franc;ais par Y. Le May
SPILLWAYS. SHOCKWAVES
AND AIR ENTRAINMENT
Review and recommendations
,
EVACUATEURS. ONDES DE CHOC
ET ENTRAINEMENT D'AIR
Synthese et recommandations
Commission Internationale des Grands Barrages - 151, bd Haussmann, 75008 Paris

Tel. : (33-1) 40 42 67 33 - Telex : 641320 ICOLD F - Fax : (33-1) 40 42 60 71
AVERTISSEMENT – EXONERATION DE RESPONSABILITE :
Les informations, analyses et conclusions auxquelles cet ouvrage renvoie

sont sous la seule responsabilité de leur(s) auteur(s) respectif(s) cité(s).
Les informations, analyses et conclusions contenues dans cet ouvrage

n'ont pas force de Loi et ne doivent pas être considérées comme un
substitut aux réglementations officielles imposées par la Loi. Elles sont
uniquement destinées à un public de Professionnels Avertis, seuls aptes à
en apprécier et à en déterminer la valeur et la portée et à en appliquer avec
précision les recommandations à chaque cas particulier.
Malgré tout le soin apporté à la rédaction de cet ouvrage, compte tenu de

l'évolution des techniques et de la science, nous ne pouvons en garantir
l'exhaustivité.
Nous déclinons expressément toute responsabilité quant à l'interprétation et

l'application éventuelles (y compris les dommages éventuels en résultant ou
liés) du contenu de cet ouvrage.
En poursuivant la lecture de cet ouvrage, vous acceptez de façon expresse

cette condition.
NOTICE – DISCLAIMER :
The information, analyses and conclusions referred to herein are the sole
responsibility of the author(s) thereof.
The information, analyses and conclusions in this document have no legal

force and must not be considered as substituting for legally-enforceable
official regulations. They are intended for the use of experienced
professionals who are alone equipped to judge their pertinence and
applicability and to apply accurately the recommendations to any particular
case.
This document has been drafted with the greatest care but, in view of the
pace of change in science and technology, we cannot guarantee that it
covers all aspects of the topics discussed.
We decline all responsibility whatsoever for how the information herein is

interpreted and used and will accept no liability for any loss or damage
arising therefrom.
Do not read on unless you accept this disclaimer without reservation.

COMMITIEE ON HYDRAULICS FOR DAMS
COMITE DE L'HYDRAULIQUE DES BARRAGES
( 1987 -1993)
Chairman/President
USA/Etats-Unis J. J. CASSIDY
Sub-Committee 1/Sous-Comite 1
Spain/Espagne J. A. HERRERAS
(coordinatorI coordinateur)
Brazil/Bresil A. RICOY DE OLIVEIRA
Chine/China SHEN CHONGGANG
France/France S. ALAM
Norway/Norvege D.K.LYSNE
Portugal/Portugal F. OLIVEIRA LEMOS
USSR/URSS V. M. SEMENKOV
USA/Etats-Unis J. J. CASSIDY
Netherlands/Pays-Bas P. KOLKMAN
Japan/ Japon M. HINO
USA/Etats-Unis A. H. MAKARECHIAN
Finland/Finlande M. AUTII (coordinator/coordinateur)
Canada/Canada M. DROUIN
Czechoslovakia/Tchecoslovaquie V. MATOUSEK
Denmark/Danemark H. MAI
Sweden/Suede L. BILLFALK
Germany/ Allemagne J. GIESECKE
Belgium/Belgique F. ZWAENEPOEL
Czechoslovakia/Tchecoslovaquie K. HAINDL
Switzerland/Suisse D. VISCHER
Finland/Finlande P. REITER (coordinator/coordinateur)
Australia/ Australie J. WINTON
Great Britain/Grande-Bretagne P.A. BACK
Spain/Espagne L. BERGA
South Africa/ Afrique du Sud A. ROOSEBOOM
2
SOMMAIRE CONTENTS
AVANT-PROPOS FOREWORD
l. INTRODUCTION l. INTRODUCTION
2. ONDES DE CHOC 2. SHOCK WA YES IN CHUTES

SUR LES COURSIERS
3. ENTRAINEMENT D'AIR 3. NATURAL AIR ENTRAINMENT

NATUREL
4. AERATEURS D'EVACUATEURS 4. SPILLWAY AERATORS

DE CRUE
3
TABLE DES MATIERES
AVANT-PROPOS .................................................................................................................... 11
I. INTRODUCTION ............................................................................................................ 12
References ...... .. ........ .. .. ...... ...... .. .. ........ .. .. ...... ...... .. .... .. ...... .. ...... .. .. ...... .. ...... .. ...... .. .. ...... .. .. .. 22
2. ONDES DE CHOC SUR LES COURSIERS .............................................................. 24

2.1. Introduction .. .. .. ...... ...... .. .. ........ .. ........ .. .... ...... .. ........ .. .. .... .. .. ...... .. ........ .. ...... .. .. ...... .. 24
2.2. Deviation brusque d'un mur .................................................................................. 26
2.3. Reduction des ondes stationnaires ........................................................................ 32
2.4. Deviation infiniment petite d'un mur .................................................................. 34
2.5. Retrecissement d'un canal ...................................................................................... 36
2.6. Elargissement d'un canal ........................................................................................ 48
2. 7. Courbes .. .. .. ...... .. .. .... .. .... .. .. .. .. ........ .. .... .. ...... ...... .. .. ...... .. .. ...... .. ........ .. ...... .. ...... .. .. .. .. .. 52
2.8. Reduction des on des a val ........ .. ...... ...... ...... .. .......... .. ...... .. ........ .. .. ...... .. .. .... .. .. .. .... 56
2.9. Notations .................................................................................................................. 64
2.10. References .. .. .. ........ .... .. ...... .. .. .. .. ...... .. .... .. .... .. ...... .. .. ...... .. .. ...... .. ...... .. .. .... .. .. .. .... .. .. .. .. 66
3. ENTRAINEMENT D'AIR NATUREL ........................................................................ 70

3.1. Introduction ................................................................................................................ 70
3.2. Mecanisme de l'entrainement d'air .......................................................................... 72
3.3. Formation de l'entrainement d'air .......................................................................... 74
3.4. :Ecoulement aere uniforme ........................................................................................ 76
3.5. :Ecoulement aere graduellement varie .................................................................... 82
3.6. Aper9u historique ...................................................................................................... 86
3.7. Notations .................................................................................................................... 88
3.8. References .................................................................................................................... 90
4. AERATEURS D'EVACUATEURS DE CRUE ............................................................ 94

4.1. Introduction .. .. ...... .... .. ...... .. .. ........ .. ...... ...... .. .... .. .. ...... .. ........ .. ...... .. ........ .. ........ .. .. ...... 94
4.2. Aeration en radier .... .. .... .. .. .......... .. .. .... ............ ........ .. ........ .. .. ...... .. .. ...... .. ...... .. .. ........ 96
4.3. Echappement de I' air ................................................................................................ 106
4.4. Conception des aerateurs ........ .. .. ...... ...... ...... ...... .. ........ .. ...... .. ...... .. .. ...... .. .. ...... .. .. .. .. I 08
4.5. Notations .................................................................................................................... 112
4.6. References .................................................................................................................... 114
4
TABLE OF CONTENTS
FOREWORD .......................................................................................................................... 11
I. INTRODUCTION ............................................................................................................ 13
References .... .. .. .. ............ ...... ...... .................... ...... .................................................... ........... 22
2. SHOCK WA YES IN CHUTES ...................................................................................... 25

2.1. Introduction .... ...... .................. ........ ...... ...... ...... ........ .. .. ...... .. ...... .. ........ ........ ........ .. .. 25
2.2. Abrupt Wall Deflection .......................................................................................... 27
2.3. Reduction of Standing Waves ................................................................................ 33
2.4. Infinitesimal Wall Deflection ................................................................................ 35
2.5. Channel Contraction ................................................................................................ 37
2.6. Channel Expansion .. ........ .. .. .......... .... .. .... .. .... .. ........ .. ........ .. ...... .. .. ...... .. ........ .. ........ 49
2.7. Bends .......................................................................................................................... 53
2.8. Reduction of Tailwater Waves .............................................................................. 57
2.9. Notation .................................................................................................................... 65
2.10. References .................................................................................................................. 66
3. NATURAL AIR ENTRAINMENT................................................................................ 71

3.1. Introduction ................................................................................................................ 71
3.2. Mechanism of Air Entrainment .............................................................................. 73
3.3. Incipient Air Entrainment ........................................................................................ 75
3.4. Uniform Aerated Flow .............................................................................................. 77
3.5. Gradually Varied, Aerated Flow .............................................................................. 83
3.6. Historical Review ...................................................................................................... 87
3.7. Notation ...................................................................................................................... 89
3.8. References .................................................................................................................... 90
4. SPILLWAY AERATORS.................................................................................................. 95
4.1. Introduction ................................................................................................................ 95
4.2. Bottom Aeration .. ...... ...... .. .. ........ .. ...... .. .... ...... .. .......... ........ ...... .. ...... .. ...... .. ........ .. .. .. 97
4.3. Air Detrainment .......................................................................................................... 107
4.4. Design of Aerators .... ........ .. .. ........ .. ...... ...... ........ .. ...... .. ...... .. ...... .. ...... .. ...... .. .. .... .... .. I 09
4.5. Notation ...................................................................................................................... 113
4.6. References .................................................................................................................... 114
5
LISTE DES FIGURES
Fig. I. - Conception des coursiers d'evacuateurs de crue.

Fig. 2. - Ondes transversales sur les coursiers du barrage d'Itaipu (Bresil/Paraguay).
Fig. 3. - Ondes transversales dues aux piles et aux chambres des vannes.
Fig. 4. - Surface d'un ecoulement aere sur un coursier.
Fig. 5. - Ecoulement aere sur un coursier a forte pente (barrage Norris).
Fig. 6. - Aerateur d'un coursier en canal vitre.
Fig. 7. - Vue de l'evacuateur de crue du barrage Grand Coulee.
Fig. 8. - a) Deviation d'un mur. b) Onde transversale dans un canal en laboratoire.
Fig. 9. - Deviation d'un mur et lignes de courant dans un canal.
Fig. 10. - Principe de l'interference d'ondes applique a un canal.
Fig. 11. - Angle de deviation d'un mur en fonction du nombre de Froude local F.
Fig. 12. - Retrecissement de canal.
Fig. 13. - Ecoulement dans un retrecissement de canal.
Fig. 14. - Conception d'un retrecissement en forme d'entonnoir.
Fig. 15. - Retrecissement en forme d'eventail a l'aval d'un evacuateur de crue.
Fig. 16. - Ecoulement dans un retrecissement plan en forme d'eventail.
Fig. 17. - Retrecissement de coursier en forme d'eventail. Notation.
Fig. 18. - Vue en perspective de la geometrie courbe du radier d'un retrecissement en forme
d'eventail.
Fig. 19. - Differences de surface d'ecoulement. a) Radier plan. b) Radier courbe.
Fig. 20. - Ondes transversales dans un elargissement brusque de canal.
Fig. 21. - Canal avec elargissement infini.
Fig. 22. - Graphique general de la surface h/h,.
Fig. 23. - Representation geometrique generate d'un elargissement.
Fig. 24. - Ecoulement torrentiel dans une courbe.
Fig. 25. - Hauteurs des ondes extremes en fonction du nombre de Froude F et de la courbure
relative b/R de la courbe.
Fig. 26. - Relevement d'une courbe au dechargeur de Saint-Chamas.
Fig. 27. - Association d'un retrecissement et d'une courbe dans un canal a ecoulement
torrentiel.
Fig. 28. - Ecoulement dans une courbe d'un canal trapezoidal.
Fig. 29. - Ecoulement torrentiel a surface libre.
Fig. 30. - Evacuateur de crue avec retenue et saut de ski. Differents ecoulements.
Fig. 31. - Projection de gouttes d'eau au-dessus de la surface d'ecoulement.
Fig. 32. - Formation d'eau blanche sur un coursier.
Fig. 33. - Profits de concentration d'air dans un ecoulement aere uniforme pour differentes
pentes de radier.
Fig. 34. - Ecoulement aere uniforme. a) Effets de la concentration d'air sur le coefficient de
frottement. b) Repartition de la concentration d'air.
6
LIST OF FIGURES
Fig. I. - Concepts of spillway chutes.

Fig. 2. - Cross-waves in chutes of Itaipu dam (Brazil/Paraguay).
Fig. 3. - Cross-waves due to piers and gate chambers.
Fig. 4. - Surface pattern of aerated chute flow.
Fig. 5. - Aerated flow on steep chute (Norris dam).
Fig. 6. - Chute aerator in a laboratory channel.
Fig. 7. - View at Grand Coulee dam spillway.
Fig. 8. - a) Wall deflection. b) Cross-wave in laboratory channel.
Fig. 9. - Wall deflection and streamline pattern in a channel.
Fig. 10. - Principle of wave interference applied to a channel.
Fig. 11. - Angle of wall deflection as a function of local Froude number F.
Fig. 12. - Channel contraction.
Fig. 13. - Flow pattern in channel contraction.
Fig. 14. - Design of funnel-shaped contraction.
Fig. 15. - Fan-shaped contraction downstream from a spillway.
Fig. 16. - Flow features in planar fan-type contraction.
Fig. 17. - Fan-type chute contraction. Notation.
Fig. 18. - View of curved bottom geometry of fan-type contraction.
Fig. 19. - Differences in flow surfaces. a) Plane bottom. b) Curved bottom.

Fig. 20. - Cross-wave formation at abrupt channel expansion.
Fig. 21. - Infinitely expanding channel.
Fig. 22. - Generalized surface plot h/h 1•
Fig. 23. - Generalized geometry of expansion.
Fig. 24. - Supercritical flow in bend.
Fig. 25. - Extreme wave heights as a function of Froude number F and relative bend
curvature b/R.
Fig. 26. - Curve banking at the " dechargeur de Saint-Chamas ".
Fig. 27. - Combination of contraction and bend in a channel with supercritical flow.
Fig. 28. - Bend flow in trapezoidal channel.

Fig. 29. - Supercritical free-surface flow.
Fig. 30. - Spillway with reservoir and ski jump. Different flows.
Fig. 31. - Droplets above the water surface.
Fig. 32. - Occurrence of white water in chute.
Fig. 33. - Air concentration profiles for uniform aerated flow and various bottom slopes.
Fig. 34. - Uniform aerated flow. a) Effect of air concentration on friction factor. b) Air
concentration distribution.
7
Fig. 35. --- Profil de la surface libre pour un ecoulement aere graduellement varie.
Fig. 36. - Zones de basse pression dues a des irregularites de surface et degats de surface.
Fig. 37. - Ecoulement a proximite d'un aerateur.
Fig. 38. - Dispositifs d'alimentation en air pour des deflecteurs.
Fig. 39. - Vue laterale d'un aerateur de coursier.
Fig. 40. - Aerateur de fond avec interface air-eau.
Fig. 41. - Representation graphique des equations (4.2) et (4.7).
Fig. 42. - Profil en long d'un coursier s'adaptant a la topographie.
Fig. 43. - Aerateur pour une forme concave de radier.
8
Fig. 35. - Free surface profile for gradually varied aerated flow.
Fig. 36. - Local low pressure zones due to surface irregularities and surface damage.
Fig. 37. - Flow in .the vicinity of aerator.
Fig. 38. - Air supply systems for chute flow deflectors.
Fig. 39. - Side view at chute aerator.
Fig. 40. - Bottom aerator with air-water interface.
Fig. 41. - Graphical display of eq. (4.2) and (4.7).
Fig. 42. - Longitudinal chute profile with adaptation to topography.
Fig. 43. - Aerator for concave bottom geometry.
9
A VANT-PROPOS
Les evacuateurs de crue et les vidanges de fond representent un pourcentage

eleve du cout total des grands barrages. Leur conception revet ainsi une grande
importance du point de vue de l'econornie et de Ia securite. Ce Bulletin qui traite de
l'etude hydraulique des ecoulernents a grande vitesse, des ecoulernents aeres, de la
cavitation et des ondes de choc, donne des recornrnandations aux ingenieurs charges
du projet des vidanges de fond et des coursiers a forte pente.
La synthese complete presentee dans ce Bulletin a ete etablie par le Sous-Cornite
no 4 du Cornite de l'Hydraulique des Barrages, cree en 1987, et plus specialernent par
un groupe de travail constitue de representants des Cornites Nationaux allernand et
suisse.
Le texte de base des quatre chapitres a ete prepare soigneusernent par
W. H. Hager (Suisse) et revu par !es rnernbres du Sous-Cornite. Les cornrnentaires
presentes par huit Cornites Nationaux ont permis d'arneliorer la redaction initiale.
Qu'ils soient tous vivement remercies.
J. J. Cassidy J. Giesecke
President du Cornite Coordinateur
de l'Hydraulique des Barrages du Sous-Comite n° 4
FOREWORD
Spillways and bottom outlets constitute a significant percent of the total cost of
major dams. Their proper design is, thus, very important from the standpoint of
economy as well as safety. This Bulletin, which considers the hydraulic concepts of
design for high-velocity flow, aeration, cavitation and shock waves, has been
prepared as a general design guideline for use of engineers engaged in the design of
bottom outlets and steep spillway chutes.
The comprehensive review, presented in this Bulletin, was performed by
Subcommittee No. 4 of the Committee on Hydraulics which was established in 1987.
The working group of the Subcommittee, which developed the successive drafts for
this Bulletin, was formed by the representatives of the German and Swiss National
Committees.
We wish to thank W. H. Hager, Switzerland, who carefully wrote the basis of
the four chapters; thanks are also due to the members of the Subcommittee who
contributed strongly to the development of the Bulletin. In addition we wisk to thank
the eight National Committees who provided comments on the initial draft, the
incorporation of which provided definite improvements.
J. J. Cassidy J. Giesecke
Chairman, Committee Coordinator,
on Hydraulics for Darns Subcommittee No. 4
11
1. INTRODUCTION
Un barrage comprend deux elements importants : l'evacuateur de crue et le

dissipateur d'energie. L'evacuateur controle le deversement, en particulier !ors des
crues, et achemine l'eau a l'aval du barrage. Le dissipateur d'energie re9oit le debit
evacue au pied du barrage et Je restitue sans danger a la riviere OU dans Un lac.
Le raccordement entre le seuil de l'evacuateur et le dissipateur est assure au
moyen de dispositifs appropries. Deux points importants sont examines dans le
Bulletin: l) en raison des grandes differences de niveau entre le seuil de l'evacuateur
et le lit de la riviere, l'energie cinetique au pied du barrage peut etre enorme et les
vitesses peuvent etre de l'ordre de 30 a 50 ms~l; 2) la distance horizontale entre les
deux elements precites est faible lorsque le deversement sur la crete tombe en chute
libre dans un bassin de reception; par contre, cette distance peut atteindre plusieurs
centaines de metres lorsque l'ecoulement se fait dans un canal d'evacuation; ce
dernier ouvrage est aussi designe sous le nom d'evacuateur avec coursier, en
soulignant le mot coursier ou a lieu l'ecoulement. On traitera ici de cet ouvrage de
raccordement.
La conception, les types, les ouvrages, le materiel d'equipement, l'exploitation,
l'entretien et Jes reparations sont examines de fa9on detaillee dans le Bulletin CIGB
no 58 ( 1987).
II y a deux types de conception possibles des evacuateurs de crue (Fig. I) :
- Le plus souvent, la rugosite des parois du coursier est faible de sorte que
l'ecoulement sera turbulent dans le regime hydraulique transitoire. II en resulte que
Jes pertes de charge par frottement sont petites et les vitesses grandes.
- Une autre conception a ete recemment recommandee par Sorensen
(1985) [8] (*), Vittal et Porey (1987) [13], et Stephenson (1988) [9] : il s'agit d'un
coursier en forme de gradins. Une partie importante de l'energie hydraulique peut
etre dissipee le long du coursier et le dissipateur d'energie peut, de ce fait, etre reduit.
La presente etude ne porte que sur les coursiers avec ecoulement a grande
vitesse, representes schematiquement sur la Fig. I a). Une grande attention a ete
portee aux coursiers de ce type au cours de ces vingt dernieres annees qui ont vu les
digues en terre de grande hauteur prendre une place importante dans la technique
des barrages.
Bien qu'il existe un grand nombre d'evacuateurs en galerie (par exemple, gaJerie
raccordee a un evacuateur de surface en tulipe), le present rapport ne traitera pas
specialement des ecouJements en conduit ferme, et n'abordera done pas Jes proble-
mes specifiques de vibration, d'hydraulique a basse pression et d'ondes de pres-
sion [2]. En outre, on n'examinera ni Jes organes de l'evacuateur, tels que vannes ou
pertuis [l], ni les dissipateurs d'energie [6] ou les questions particulieres relatives aux
modeles [7]. Sauf indication contraire, la discussion sera limitee aux evacuateurs
avec coursiers a ciel ouvert, d'axe rectiligne.
(*) Les chiffres entre crochets renvoient a la liste des references de chaque chapitre.
12
1. INTRODUCTION
Two major parts of a dam are the spillway and the energy dissipator. The
spillway controls the overflow, particularly during periods of floods and conveys the
water away from the dam. The energy dissipator, on the other hand, receives the flow
at the foot of the dam and returns it safely back to the adjacent river or lake.
Between the spillway crest and the dissipator lies a certain region to be
connected by appropriate means. Two major points are of concern here : 1) due to
large differences of elevation between spillway crest and original river the kinetic
energy at the foot of the dam may be large and velocities can be of the order of 30
to 50 ms~ 1 ; 2) the horizontal reach between the two elements is small when the flow
is guided over the dam crest and falls through the air on an impact basin. The
distance can amount to hundreds of meters when the flow is in a spilling channel.
The latter is also referred to as spillway with stress on the word way where the
spilling occurs, or better to a chute. This connecting element will be considered
herein.
The principles for spillway design, including types, structures, equipment,

operation, maintenance and repairs are comprehensively described in the
ICOLD Bulletin 58 (1987).
Spillways may be designed according to two concepts (Fig. 1) :
- Usually, the boundary roughness of the chute is kept small such that the flow
will be turbulent in the hydraulic transitional regime. The resulting friction losses are
thus small and the velocities large.
- A recent design recommended by Sorensen (1985) [8] (*), Vittal and Porey
(1987) [13) as well as Stephenson (1988) [9] is based on cascades of drop structures.
As a result, a significant portion of hydraulic energy may be dissipated along the
chute, and the dissipator size may be considerably reduced.
The present study accounts only for high velocity chutes, as shown in Fig. 1a).
Such chutes have received considerable attention in the last twenty years when the
high earth dam became an efficient type of structure in dam engineering.
Although there exist a number of tunnel spillways such as those discharging

from morning glory spillways, this report does not refer to closed conduit flow. Nor
does it deal with specific questions relating to vibration, low pressure hydraulics, and
pressure surges [2]. Moreover, neither the spillway elements, such as gates or
openings [l], nor the energy dissipator [6] or particular questions relating to
model-prototype [7] will be considered. The discussion will be restricted to the open
chute spillway of straight axis unless otherwise stated.
(*) Numbers between brackets refer to references of each chapter.
13
a)
- ~-J~ -
~-
Le chapitre 2 est consacre aux ondes de choc qui constituent des perturbations
stationnaires a la surface libre d'un ecoulement torrentiel. Le terme « onde de choc »
tire son origine de la dynamique des gaz et, en hydraulique, le terme « ondes
transversales » est plus approprie. La Fig. 2 donne des exemples caracteristiques de
formation d'ondes transversales. Dans Jes deux cas, ces ondes prennent naissance a
l'extremite aval des piles de l'evacuateur de crue et se manifestent sur une grande
longueur du coursier. Comme on le verra plus loin, ces ondes stationnaires sont
provoquees par toutes perturbations situees sur le coursier, telles que changement de
direction, retrecissements et elargissements, OU meme piles (Fig. 3). Le chapitre 2
aborde brievement la question des ondes transversales et indique Jes possibilites de
reduire Jes ondes de hauteur excessive.
14
Fig. 1
Concepts of spillway chutes.
Conception des coursiers d'evacuateurs de crue.
a) Smooth boundary roughness, combined with a} Faible rugosite des parois : ecoulement a
high-velocity chute flow and terminal energy grande vitesse sur le coursier et dissipateur
dissipator. d'energie ii /'extremite du coursier.
b) Rough chute bottom, combined with highly b) Coursier en forme de gradins (forte rugosite) ·
turbulent cascade flow and partial terminal ecoulement tres turbulent et dissipation
energy dissipation. d'energie reduite au pied.
Chapter 2 focuses on shock waves, which are standing free surface

perturbations of supercritical flow. The term shock wave originates from gas
dynamics, and one may speak more appropriately of cross-waves in hydraulics. Fig. 2
shows typical examples of cross-wave generation. These are generated in both cases
at the rear end of spillway piers and extend far down the chute. As will be seen, these
standing waves are due to any perturbations located in the chute, such as change of
direction, contractions and expansions, or even piers (Fig. 3). In chapter 2, a brief
outline on cross-waves is provided and reference is made to possibilities for reducing
waves of excessive height.
15
a)
Fig. 2
Development of cross-waves in chutes of Itaipu Dam (Brazil/Paraguay).
Formation d'ondes transversales sur /es coursiers de /'evacua/eur de crue du harrage d'/taipu
( Bresi l/ Paraguayj.
a) Downstream view (VA W D 6200). a) Vue d'aval (VAW D 6200).
b) Upstream view (YAW D 4249). h) Vue d'amon/ (VAW D4249}.
16
'
Fig. 3
Development of cross-waves due to piers and gate chambers
(ASCE Task Committee, 1961) [II].
Formation d'ondes transversales dues aux piles et aux chamhres des vannes
(ASCE, 1961) [llj
17
Le chapitre 3 traite de l'entrainement d'air nature! dans Jes ecoulements sur Jes
coursiers. Etant donne que Jes publications sur ce sujet sont tres abondantes et
qu'aucune synthese n'existe jusqu'a maintenant, un bref aperyu general historique
sera presente sur cette question. Une distinction est faite entre la zone evolutive,
)'aeration naissante, l'ecoulement aere graduellement varie et I'ecoulement uni-
forme, par analogie avec Jes courbes de remous d'un ecoulement non aere dans un
canal. Des termes en seront tires pour Jes caracteristiques d'un ecoulement varie
intervenant dans I'etude d'un coursier : profondeur du melange air-eau, concentra-
tion moyenne de l'ecoulement air-eau, courbes de remous modifiees, ... Les meca-
nismes de I'entrainement d'air sont egalement examines. La Fig. 4 montre un
ecoulement entierement aere sur un coursier et indique la configuration de surface
agitee et tres aeree. Les photographies de la Fig. 5, publiees par Hickox en 1945 [4],
donnent un exemple d'ecoulement aere sur un coursier.
Fig. 4
Surface pattern of aerated chute flow, exposure time I 0 '
(according to Straub and Lamb, 1953) [IO].
Surface d'un ecoulement aere sur w1 coursier, temps d'exposition JO ' s
(d'apres Straub et Lamb, 1953) [JO}.
Enfin, le chapitre 4 est consacre aux aerateurs d'evacuateurs de crue. De tels

organes luttent efficacement contre Jes degats par cavitation. L'aerateur fournit de
I'air supplementaire a Ia couche inferieure de I'ecoulement sur le coursier et
fonctionne ainsi comme un amortisseur de choc [12]. La surface du beton est, de ce
fait, protegee. On examine egalement le probleme de I'augmentation de hauteur des
18
In chapter 3, the natural air entrainment in chute flow is discussed. Since the
literature on this topic is quite extensive, and no systematic review exists, a historical
summary is presented. Distinction between the developing region, incipient aeration,
gradually varied aerated flow and uniform flow is made, much in analogy to
backwater curves in non-aerated channel flow. Expressions are derived for various
flow characteristics involved in chute design, such as the mixture flow depth, the
average concentration of air-water flow and the modified backwater curves. The
mechanisms of air entrainment are also discussed. Fig. 4 shows a fully aerated chute
flow and reveals the rough and highly porous surface pattern. An example of aerated
chute flow is presented with reference to photographs originally published by
Hickox in 1945 [4] (Fig. 5).
Fig. 5
Aerated flow on steep chute. Norris Dam Spillway
where the effect of incipient aeration is well seen both for :
Ecoulement aere sur un coursier a forte pente. Evacuateur de crue du barrage Norris,
ou /'on voit bien la naissance de /'aeration dans /es deux cas :
a) Small discharges. a) Debit faible.

b) Large discharges b) Debit e/eve
(according to Hickox, 1945) [41. (d'apres Hickox, 1945) [4].
Chapter 4 concludes the present study with spillway aerators. Such elements
were found adequate to counter cavitation damage. The aerator supplies the bottom
layer of chute flow with additional air, resulting in a sort of shock-absorber [12]. The
concrete surface is thus protected from damage. The problem of increasing the
channel side walls to account for the bulkage of the air-water mixture is also
19
murs lateraux du canal pour tenir compte du gonflement du melange air-eau. Le
chapitre presente des methodes pour l'etude des conduits d'alimentation en air et
pour !'estimation de la quantite d'air entrainee suivant la geometrie de l'aerateur. On
donne egalement des informations sur la perte d'air le long du couriser, ainsi que des
criteres pour determiner !'emplacement des aerateurs. La Fig. 6 represente un
aerateur dans un canal vitre et montre que l'entrainement d'air a partir d'aerateurs
de fond est un phenomene de « cisaillement ».
La vue d'un ecoulement sur un coursier est tres impressionnante (Fig. 7).
L'ecoulement est contr6le sur le seuil de l'evacuateur puis evacue par le coursier. A
l'extremite des piles des vannes, l'ecoulement commence souvent a entrainer de I' air,
l'eau devenant « blanche ». Suivant la pente du coursier et la valeur du debit evacue,
ce phenomene peut se produire plus a l'aval. Une fois l'ecoulement aere, un bruit
effrayant, des embruns et des rafales de vent se manifestent. L'energie de cette masse
d'air et d'eau, tres turbulente, doit etre efficacement dissipee dans un bassin
d'amortissement situe au pied du barrage.
Ow= 107.21/s
Fr= 7.36
,:\p=0
13 =fimax
Fig. 6
Chute aerator in a laboratory channel (VAW, 33173/12).
Aerateur d'un coursier en canal vitre (VAW 33173/12).
20
discussed. The chapter introduces procedures for designing both the air supply ducts
and determining the quantity of air entrained by a certain aerator geometry.
Guidance is given with regard to the loss of air along the chute, and criteria for
aerator positioning. Fig. 6 shows a lab aerator and stresses the fact that air
entrainment by bottom aerators is a shearing phenomenon.
A tail water view of a chute is very impressive (Fig. 7). The flow is controlled by
a spillway crest and then directed along the chute. At the end of gate walls, the flow
often starts to entrain air, that is the water becomes " white ". Depending on the
slope of chute and the discharge, this process may be shifted more downstream.
Once the flow is aerated, it is accompanied by a tremendous noise, spray
development and gusty winds. The energy of this highly turbulent mass of air-water
mixture must then be calmed down and be effectively dissipated by an appropriate
stilling basin located at the toe of dam.
........
Fig. 7
View at Grand Coulee Dam spillway (Gumensky, 1949) [3].

Vue de f'evacuateur de crue du barrage Grand Coulee (d'apres Gumensky. 1949) (3/.
21
REFERENCES
l. CASSIDY J. J. and ELDER R. A. (1984): "Spillways of high dams". Developments
in Hydraulic Engineering 2, 159-182. P. Novak ed. Elsevier Applied Science
Publishers, London and New York.
2. FALVEY H. T. (1990) : "Cavitation in chutes and spillways". Water Resources
Technical Publication Engineering Monograph 42. US. Dept. Interior, Bureau of
Reclamation, Denver.
3. GUMENSKY D. B. (l 949) : " Air entrained in fast water affects design of training
walls and stilling basins ". Civil Engineering 19 ( 12), 35-37; 93.
4. HICKOX G. H. ( 1945) : " Air entrainment on spillway faces ". Civil Engineering
15 ( 12), 562-563.
5. ICOLD (1987) : "Spillways for dams". !COLD Bulletin 58, Paris.
6. LOCHER F. A. and Hsu S. T. (1984) : " Energy dissipation at high dams ".
Developments in Hydraulic Engineering, Vol. 2, 183-238, P. Novak ed. Elsevier
Applied Science Publishers, London and New York.
7. NOVAK P. and CABELKA J. (1981): "Models in Hydraulic Engineering". Pitman,
Bosten.
8. SORENSEN R. M. (l 985) : " Stepped spillway hydraulic model investigation ".
Journal of Hydraulic Engineering ASCE l l l (12), 1461-1472; 113 (8), 1095-1097.
9. STEPHENSON D. (1988): "Stepped energy dissipators ".The Internal. Symposium
on Hydraulics for High Dams, Beijing, 1228-1235.
JO. STRAUB L. G. and LAMB 0. P. (1953). " Experimental studies of the air
entrainment in open channel flow ". V. JAHR Congress Minneapolis, 425-437.
11. Task Committee on Air Entrainment in Open Channels (1961)." Aerated flows
in open channels ". Proc. ASCE, J. Hydraulics Division 87 (3), 73-86; 87 (6),
221-231.
12. VISCHER D. ( 1985). " Die Schu13rinnenbeliiftung als Ma13nahme gegen
Kavitationserosion ". Osterreichische Wasserwirtschaft 37 (3/4), 117-123.
13. VITTAL N. and POREY P. D. (1987)." Design of cascade stilling basins for high
dam spillways ". Journal of Hydraulic Engineering ASCE, 113 (2), 225-237.
22
2. ONDES DE CHOC SUR LES COURSIERS
2.1. INTRODUCTION
Le nombre de Froude local
IF = V/c (2.1)
ou : V est la vitesse locale,
C, Ja vitesse de propagation OU ceJerite des ondes de surface, joue Un role
important dans Jes ecoulements a surface libre en canal. Pour un ecoulement
unidimensionnel, la celerite depend uniquement de la geometrie de la section
transversale et peut etre exprimee comme suit :
C = Vg A/Bs (2.2)
A, surface de la section transversale,
B" largeur en surface,
g, acceleration de la pesanteur.
Dans le cas d'un canal rectangulaire de largeur b et de surface transversale
A = bh, ou h est la profondeur de l'ecoulement, !'equation (2.2) donne :
C=Vg·h
Dans le cas d'un ecoulement peu profond bidimensionnel, la celerite peut
changer dans Jes directions transversale et longitudinale suivant la surface locale de
l'ecoulement.
Lorsque la vitesse locale est egale a la celerite, l'ecoulement est qualifie de
critique et IF= I. Lorsque IF< I, l'ecoulement est fluvial; lorsque IF> I, l'ecou-
lement est torrentiel. Le nombre de Froude dans la dynamique des fluides est
analogue au nombre de Mach en aerodynamique.
Les ecoulements fluvial et torrentiel ont un comportement tout a fait different,
en ce qui concerne l'etude mathematique du probleme. L'ecoulement avec IF < I a
un comportement elliptique, c'est-a-dire que des problemes aux limites doivent etre
resolus. Au contraire, l'ecoulement avec IF > I a un comportement hyperbolique
entrainant des particularites du type ondes. Ces dernieres constituent des problemes
caracteristiques le long du canal dans Jes deux directions, amont et aval. Chaque
obstacle dans le canal - par exemple, pile, deviation ou jonction de mur - cree une
perturbation distincte le long de courbes caracteristiques bien definies. Les perturba-
tions provoquent des ondes stationnaires a la surface de l'eau : creux ou - ce qui
est plus dangereux - surelevations. De telles ondes sont appelees ondes transversa-
les ou - par analogie avec la dynamique des gaz - ondes de choc. Dans Jes
ecoulements permanents, ces ondes ne se deplacent pas sauf pour Jes variations
turbulentes, et necessitent une revanche supplementaire. II y a done lieu de recher-
cher des moyens permettant de reduire de telles ondes.
Dans ce qui suit, la formation d'ondes transversales sera examinee, tout
d'abord, pour Jes deviations de mur de faible convergence dans un canal infiniment
large. La notion d'onde transversale sera introduite et Jes caracteristiques a l'amont
et a l'aval de l'onde seront determinees. Le cas du canal de largeur finie sera ensuite
examine et la forme des ondes de surface sera etudiee.
24
2. SHOCK WAVES IN CHUTES
2.1. INTRODUCTION
Open channel flows are significantly influenced by the local Froude number :
IF = V/c (2.1)
where : V is the local velocity and
c the propagation velocity or celerity of shallow water waves. The celerity
of a one-dimensional flow depends solely on the cross-sectional geometry
and may be expressed as
c = Vg AIB, (2.2)
A is the cross-sectional area, B, the surface width and g the gravitational
acceleration.
For a rectangular channel of width b and cross-sectional area A bh, where

h is the flow depth, Eq. (2.2) obtains :
c=Vg·h
For a two-dimensional shallow flow, the celerity may change both m the
transverse and longitudinal directions according to the local flow surface.
Flows where the local velocity is equal to the celerity are referred to as critical,
and IF = I. When IF < I, a flow is subcritical, and supercritical flows occur if
IF > I. The Froude number in fluid dynamics is analogous to the Mach number in
aerodynamics.
Subcritical and supercritical flows behave quite differently as regards the
mathematical solution process. Flows with IF < I behave elliptically, that is
boundary value problems must be solved. In contrast, flows for which IF > I are
hyperbolic problems that involve wave-type features. The latter are typical marching
problems both in the forward and backward directions. Each obstacle in the channel
- such as a pier, a wall deflection or a junction - creates a distinct perturbation
along well defined curves which are referred to as characteristics. The perturbations
result in standing water surface waves that are either troughs or - more
dramatically - wave peaks. Such waves are also termed cross-waves or - according
to gas dynamics - shock waves. For steady flow, these waves do not move except
for turbulent fluctuations, and demand additional freeboard. Means to reduce
cross-waves are thus of engineering concern.
In the following, the development of cross-waves will be first outlined for small
converging wall deflections in an infinitely wide channel. The notion of cross-wave
is introduced and quantities up- and downstream of a cross-wave shall be related.
Then, a channel of finite width is considered, and the pattern of surface waves is
studied.
25
On poursuivra par l'etude des retrecissements des canaux, en forme d'enton-
noir, d'eventail ou de tuyere. Cette derniere forme de retrecissement ne sera
d'ailleurs pas examinee. Les retrecissements en forme d'entonnoir sont efficaces
pour des rapports de contraction relativement eleves, c'est-a-dire lorsque la largeur
aval n'est pas inferieure a la moitie de la largeur amont (a l'entonnement). Pour des
rapports de contraction plus petits, comme dans le cas des entrees d'evacuateurs de
crue, le retrecissement en forme d'eventail est recommande.
Une autre disposition a examiner est l'elargissement du canal, bien qu'il attire
beaucoup mains !'attention que le retrecissement. La methode d'etude classique sera
presentee et !es resultats d'essais sur modele seront examines. Des ecoulements dans
des elargissements peuvent se produire sur des coursiers !ors de manreuvre de vannes
ma! equilibree ou dans des zones de transition entre pente raide et pente mains raide
d'un coursier.
2.2. DEVIATION BRUSQUE D'UN MUR
Soit un ecoulement torrentiel parallele dans un large canal rectangulaire dont

la pente du radier est telle que le frottement sur le mur est compense (ecoulement
uniforme) [6, 28]. En uncertain point, lemur est devie vers l'ecoulement, d'un angle
+ 8 (Fig. 8 a). La deviation du mur cause une perturbation de l'ecoulement le long
d'une ligne droite faisant un angle ~ par rapport a la direction initiale de l'ecoule-
ment. L'angle ~ est appele « angle de choc » et la ligne correspondante « onde
transversale ». L'angle de choc ~est toujours plus grand que !'angle de deviation du
mur e.
...
V1
-7
"" ."-/ /v1n
~t
I< --
-,
~
--......
....-...
I
a)
th1
-V1
i - --
V2
--.-
-l\1
26
Next, channel contractions of funnel-shaped, fan-shaped, and nozzle-shaped
geometries are discussed. The latter will not be considered further. Funnel-shaped
contractions are efficient for relatively large contraction ratios, that is when the
downstream width is not less than half of the approaching width. For smaller
contraction ratios, as may typically occur at spillway entrances, the fan-shaped
contraction is recommended.
The next structure to be considered is the channel expansion, although it

received much less attention than the channel contraction. The classical approach is
presented, and the results of model tests will be discussed. Flows in channel
expansions may occur in chutes with unbalanced gate operation or at transitions
from steep to less steep chute portions.
2.2. ABRUPT WALL DEFLECTION
Rouse and Chow [6, 28] considered parallel supercritical flow in a wide
rectangular channel whose bottom slope is such that wall friction is compensated
(uniform flow). At a certain point, the wall is deflected towards the flow by the angle
+ 0 (Fig. 8 a). The wall deflection causes a perturbation of flow along a straight
curve whose deviation from the original direction of flow is p. P is referred to as
shock angle, and the corresponding curve as cross-wave. The shock angle Pis always
larger than the wall deflection 0.
Fig. 8
a) Wall deflection and notation for cross-wave. a) Deviation d'un mur et notation relative a
l'onde transversale.
b) Cross-wave in laboratory channel. b) Onde transversale dans un canal en
laboratoire.
27
En designant respectivement par Jes indices (I) et (2) Jes positions a l'amont et
a l'aval de l'onde stationnaire, !'application du principe des quantites de mouvement
et de !'equation de continuite donne [6J :
(2.4)
h2 = tg ~ (2.5)
h1 tg (~ - 8)
(2.6)
OU:
IF1 = V1/(gh1)112 (2.7)
est le nombre de Froude de l'ecoulement entrant et
y = h2lh1 (2.8)
est le rapport des profondeurs de l'ecoulement.
Le systeme d'equations a ete etabli en supposant que l'eau etait peu profonde,
c'est-a-dire que la pression etait hydrostatique et que Jes distributions des vitesses
etaient uniformes.
Pour des parametres donnes h 1, V1et 8, on peut calculer h 2, ~et IF2. Le systeme
d'equations ne conduit pas a des solutions precises, mais ii existe un graphique [15)
permettant d'obtenir les valeurs recherchees.
Dans le cas OU IF1 sin~ > I, I' equation (2.4) peut etre remplacee par !'equation
approchee suivante [10] :
(2.9)
et Hager [9) a montre que :
~ = e+ 3_ IF11 (2.10)
2i/2
Pour~ < n/4 (45°) et IF 1 > 2, Jes valeurs caracteristiques tirees de !'equation
(2.10) presentent un ecart inferieur a 2° par rapport a la valeur exacte. De plus, on
peut representer le nombre de Froude aval IF2 par !'equation approchee :
(2.11)
28
Let indices (I) and (2) denote locations upstream and downstream of the
standing wave, respectively. The application of the momentum principle and the
continuity equation yields [6] :
(2.4)
hz = tg ~ (2.5)
h1 tg (~ - 0)
IF}= y- 1
[ IF 12 - /y (Y - l)(Y + 1) 2] (2.6)
in which
IF 1 = V/(gh 1) 112 (2.7)
is the approaching Froude number and
(2.8)
is the ratio of flow depths.
The set of equations was derived assuming shallow water, that is hydrostatic
pressure and uniform velocity distributions.
For given parameters h 1, V1 and 0, one may compute h 2, ~and IF 2 • The system
of equations cannot be solved explicitly but Ippen [15] prepared a diagram from
which the unknown quantities may be determined graphically.
For the case IF 1 sin ~ > l, Eq. (2.4) may be approximated by [IO] :
,;- l
Y = v 2 IF1 ·sin ~ - - (2.9)
2
and Hager [9] showed that :
~ = e + _l_ IF1 1 (2. l 0)

2{2
For ~ < n/ 4 (45°) and IF 1 > 2, typical values of Eq. (2. l 0) deviate less than 2°
from the " exact " value. Further, an approximation for the downstream Froude
number IF2 is :
(2.1 l)
29
Les equations (2.4) a (2.6) ont ete verifiees experimentalement par Ippen et
Harleman [17]. Pour I < Y < 2, on a observe des ressauts ondules tandis que Jes
fronts des ondes transversales etaient raides pour y > 2 et le passage de h1 a h2
presque brutal. Pour Y > 2, le present systeme d'equations concordait bien avec Jes
observations.
Exemple 2.1
Soit un coursier dont la pente du radier est S = 7 %, le debit unitaire (par metre
0
de largeur) q = 15 m 2 s- 1, le coefficient de rugosite n = 0,011 sm- 113 • Quelle est la

profondeur de l'ecoulement a l'aval d'un mur de deviation d'angle e = 4°, si
l'ecoulement entrant est uniforme?
1) La profondeur de l'ecoulement uniforme d'apres !'equation de Manning-
Strickler est :
et !'on a :
VI = q/hl = 20 ms-I
IF 1 = Vil(gh 1) 112 = 7,37
2) L'angle de l'onde transversale p = 9,5°, d'apres !'equation (2.10).
3) La profondeur de l'ecoulement a l'aval de l'onde transversale est determinee
comme suit :
Y = 1,22 d'apres !'equation (2.9)
d'ou h 2 = 1,22 x 0,75 m = 0,915 m
4) Le nombre de Froude IF 2 = 6,6, d'apres !'equation (2.11); l'ecoulement aval
reste done torrentiel.
Une application du mur de deviation est le canal presentant un changement de
direction brusque mais faible. La Fig. 9 represente un canal dont le nombre de
Froude de l'ecoulement a l'entree IF 1 > 2. Aux points A et B, le canal change de
direction d'un angle 0. Au point A, le changement de direction rentre dans
l'ecoulement : la deviation du mur est qualifiee de « positive » et la profondeur de
l'ecoulement augmente. L'ecoulement est identique a celui precedemment examine
(Fig. 8). Au contraire, la deviation du mur au point B est « negative » et une onde
negative se forme avec diminution de la profondeur de l'ecoulement. Cette onde est
differente de l'onde positive du fait que le changement de profondeur est progressif
et non brutal.
30
Eqs. (2.4) to (2.6) were experimentally verified by lppen and Harleman [17]. For
< Y < 2, undular jumps were observed whereas the fronts of cross-waves for
Y > 2 were sharp, and the transition from h 1 to h 2 almost abrupt. For Y > 2, the
present system of equations compared well with observations.
Example 2.1
Given a chute of bottom slope S0 = 7 %, discharge per unit width q = 15 m2 s- 1,
roughness coefficient n = 0.011 sm- 113 • What is the flow depth downstream of a wall
deflection of e = 4°, if the approaching flow is uniform?
I) Uniform flow depth according to Manning-Strickler equation
and thus
VI = q/hl = 20 ms-I
and IF 1 = V/(gh 1) 112 = 7.37
2) Angle of cross-wave p = 9.5° according to Eq. (2.10).
3) Flow depth downstream from cross-wave :
Y = 1.22 according to Eq. (2.9), that is

h2 = 1.22 x 0.75 m = 0.915 m
4) Froude number IF2 = 6.6 according to Eq. (2.11 ), thus the tail water flow
remains supercritical.
An application of the wall deflection as previously considered is the channel
with an abrupt but small change of direction. Fig. 9 shows a channel with an
approaching Froude number IF 1 > 2. At points A and B, the channel changes
direction by 8. At point A, the change is into the flow and reference is made to a
positive wall deflection associated with an increase of flow depth. The flow is similar
to that previously discussed (Fig. 8). In contrast, the wall deflection at point B is
negative, and a negative wave develops associated with a decrease in flow depth.
This wave is different from the positive wave as the depth change is gradual and not
abrupt.
Fig. 9
Wall deflection and streamline pattern in a channel of finite width.
Deviation d'un mur et representation des /ignes de courant dans un canal de largeur finie.
( - ) Positive wave. ( - ) Onde positive.
(- - -) Negative wave. (- - -) Onde negative.
31
La zone a l'amont des ondes transversales, c'est-a-dire a gauche de ACB dans
la Fig. 9, n'est pas influencee par le changement de direction et Jes lignes de courant
restent paralleles au canal d'entree. Dans Jes zones ACE et BCD, Jes lignes de
courant sont paralleles aux murs du canal aval et Jes profondeurs de l'ecoulement
augmentent et diminuent respectivement.
Le point C est !'intersection des deux ondes transversales. Toutes deux pene-
trent dans des zones perturbees, mais aboutissent - en raison de !'inversion - a la
zone d'etat initial. II en resulte que Jes caracteristiques dans la zone CDEF sont
egales a h 1 et IF 1• Cependant, Jes lignes de courant ont devie de 2 8 par rapport a
la direction de l'ecoulement initial.
A l'aval de la ligne DFE, une inversion se produit de nouveau : Jes lignes de
courant deviennent paralleles aux murs. La profondeur de l'ecoulement augmente
dans la zone DFH et diminue dans la zone EFG. Le long du mur AEG, la surface
s'eleve tout d'abord brusquement puis tombe au-dessous de la profondeur initiale h 1•
L'ecoulement devient ainsi agite dans le canal, ce qui necessite des murs plus hauts
que ceux de la zone d'entree. Une partie importante de l'etude des coursiers porte
sur la reduction des sommets des ondes en vue d'obtenir un ecoulement presque
uniforme.
2.3. REDUCTION DES ONDES STATIONNAIRES

Au cours des cinquante dernieres annees, un certain nombre de dispositifs
relativement simples ont ete mis au point pour reduire Jes sommets des ondes locales.
II est evident que ces dispositifs ne peuvent supprimer totalement Jes ondes dans
n'importe quelles conditions d'ecoulement, pour Jes raisons suivantes :
- Jes couches limites et Jes effets de courbure peuvent localement modifier le
modele mathematique;
- Jes effets de la pente et le frottement sur lemur ne sont pas pris en compte;
- Jes debits differents du debit de projet ne satisfont pas aux conditions pour
!'elimination des ondes.
Ces facteurs ne peuvent etre pris en compte que dans des grands modeles
physiques OU dans une etude sur modeJe numerique de l'ecouJement, tenant COmpte
des effets de courbure des lignes de courant. Les methodes fondamentales de
reduction des ondes guideront le choix des techniques de modelisation.
Le moyen le plus efficace pour reduire Jes ondes transversales est base sur le
principe de !'interference d'ondes. Selon ce principe, une onde negative est produite
pour remplacer une onde positive identique. Par exemple, une deviation de mur
positive dans un canal de largeur finie, telle que representee sur la Fig. I 0 a), conduit
a une augmentation de la profondeur d'eau h 2 > h 1• Dans la Fig. 10 b), !'inverse se
produit conduisant a une diminution de la profondeur d'eau a l'aval.
a)µv
32
The domain upstream of the cross-waves, that is on the left of ACB in Fig. 9,
is not influenced by the change of direction, and the streamlines remain parallel to
the approaching channel. In the two zones ACE and BCD, the streamlines are
parallel to the walls of the downstream channel, and the flow depths have increased
and decreased, respectively.
Point C is the intersection of two cross-waves. Both waves continue into
perturbed domains, but lead - due to inversion - to the domain of original state.
As a result, the characteristics in CDEF are equal to h 1 and IF 1• However, the
streamlines are deflected by 2 8 relative to the original flow direction.
Downstream from the curve DFE, inversion occurs again, and the streamlines
become parallel to walls. The flow depth increases in domain DFH, and decreases
in domain EFG. Along the wall AEG, the surface first rises suddenly, and then falls
below the initial depth h 1• The flow in the channel thus becomes tumbling and walls
higher than in the approaching zone are needed. An important consideration in the
design of chutes is the reduction of wave peaks to obtain almost uniform flow.
2.3. REDUCTION OF STANDING WAVES

During the last fifty years, a number of relatively simple means to reduce local
wave peaks were developed. Clearly, these may not be considered as a complete
removal of waves for any flow condition, since :
- boundary layers and curvature effects may locally modify the computational
model;
- slope effects and wall friction are not considered, and
- discharges other than the design discharge do not satisfy the conditions for
wave elimination.
These implications can only be accounted for when large physical models are
used or the flow is modelled numerically by an approach accounting for streamline
curvature effects. The basic methods of wave reduction should guide the more
involved modelling techniques.
The most effective means in reducing cross-waves is based on the principle of
wave interference. According to this principle, a negative wave is produced to
compensate for an identical positive wave. For example, a positive wall deflection
in a channel of finite width according to Fig. 10 a) is associated with an increase of
flow depth h 2 > h 1• In Fig. 10 b) the inverse is shown which leads to a decrease of
downstream flow depth.
Fig. 10
Principle of wave interference applied to channel of finite width.
Principe de /'interference d'ondes applique a un canal de largeur finie.
a) Increasing flow depth and IF 1 > IF,. a) Augmentation de la profondeur de /'ecoule-
ment et IF, > IF,.
b) Decreasing flow depth and IF 1 < IF,. b) Diminution de la profondeur de /'ecoulement et
IF, < IF,.
33
2.4. DEVIATION INFINIMEMENT PETITE D'UN MUR
Precedemment, la deviation du mur 8 avait une valeur finie. Lorsque 8 - 0,

Jes fonctions h et V peuvent etre considerees comme continues. Selon von Karman
[37], on peut ecrire la relation suivante entre la profondeur de l'ecoulement et la
vitesse locale :
tg~ d8
2
dh = (: ) (2.12)
Si, de plus, la charge de l'ecoulement H = h + V2/(2 g) est consideree com me

constante, rapportee au fond rectiligne du canal, la relation :
1 (2.13)
8 - 8 1 = v'3·arctg [ v'3 !)112] -
(IF2 _
arctg [
(IF2 _
]
!)1/2
est obtenue pour la fonction 8 (IF). lei, 8 1 = 8 (IF 1) en [rad] correspond a la

constante d'integration dependant du nombre de Froude IF 1 de l'ecoulement entrant.
Pour IF > 2, !'approximation
8 - 81 = 2 (IF2 - 1,8) (2.14)

(IF2 - 1)312
est meilleure que 5 % (Fig. 11 ). Pour IF > 4, !'equation (2.14) peut etre simplifiee sous
la forme suivante :
2
(2.15)
IF
60°
40°
e-a1
20°
IF
0
1 3 5 7 9 11
34
2.4. INFINITESIMAL WALL DEFLECTION
The wall deflection 0 was previously assumed finite. For 0 -. O one may
consider the functions h and V as continuous. According to von Karman (37] the
following relation holds between flow depth and local velocity :
(~ ) tg~ d0
2
dh = (2.12)
Further, if the energy head H = h + V2I (2 g) is assumed constant relative to the

straight channel bottom, the relation :
(2.13)
e- e I = i/3. arctg [
(IF2 -
i/3 1)112] - arctg [
(IF2 -
I
1)112
]
is obtained for the function 0 (IF), where 0 1 = 0 (IF 1) in [rad] corresponds to the
integration constant depending on the approaching Froude number IF 1. For IF > 2
the approximation
e- 01 = 2 (IF2 - 1.8) (2.14)

(IF2 _ l )312
holds better than 5 % (Fig. 11 ). For IF > 4, Eq. (2.14) can even be simplified to :
2
= - (2.15)
IF
Fig. 11
Angle of wall deflection (8 - 8 1) as a function
of local Froude number IF according to Eq. (2.13).
Angle de deviation d'un mur (e - e 1) en fonction
du nombre de Froude local IF, d'apres /'equation (2.13).
Approximations : (- - -) Eq. (2.14), Approximations : (- - -) equation (2.14),
( ... ) Eq. (2.15). (. ..) equation (2.15).
35
Pour de grandes valeurs de IF, l'effet de la variation de la profondeur de
l'ecoulement est faible et on peut admettre la vitesse constante V = V 1 au lieu de
H = H 1• Dans ce cas, la profondeur locale d'ecoulement h varie comme suit (22]
(2.16)
, IF 1 = -V1-
OU
Vgh1
Une deviation de mur positive + 8 conduit a une augmentation de profondeur
d'ecoulement h/h 1 > I. Des applications des methodes de von Karman et de Knapp
sont donnees dans le paragraphe 2.7.
2.5. RETRECISSEMENT D'UN CANAL
2.5.1. Types de retrecissement

Un retrecissement de canal dans un ecoulement torrentiel se presente aux endroits
OU
la pente du radier est augmentee;
la rugosite est reduite, ou
le debit augmente.
Le premier cas est un probleme technique. Un retrecissement bien corn;:u se
caracterise par un ecoulement presque uniforme dans la partie de canal retrecie.
Sinon, de fortes ondes transversales se formeront et se manifesteront loin a l'aval;
des murs de canal de grande hauteur seront necessaires. En outre, un tel ecoulement
peut entrainer beaucoup d'air et un bassin d'amortissement situe a l'extremite du
coursier peut souffrir d'un ecoulement entrant dissymetrique.
Un probleme important est l'etude du retrecissement pour passer d'un canal
d'entree de largeur b 1 a un canal de largeur b 3 < b 1• Taubert [33, 34] a distingue Jes
retrecissements en forme d'entonnoir, en forme d'eventail et en forme de tuyere
(Fig. 12). Le premier type a une forme polygonale entre la partie large et la partie
etroite du canal. Dans la forme en eventail, on utilise des arcs circulaires dont la
courbure est orientee vers l'ecoulement, alors que dans la forme en tuyere deux arcs
circulaires sont disposes en S.
l
I
I
I I
a) b) c)
+ +
36
For large IF the effect of variation in flow depth is small and one may assume
constant velocity V = V1 instead of H = H 1• Then, the local flow depth h varies as
(22] :
h
- =
IF2 · 2
I ·Sill (2.16)
h1
·in which IF 1 = V/(gh 1) 112 •

For positive wall deflection + e the relative flow depth is h/h1 > I.
Applications of the approaches by von Karman, and Knapp are given in section 2.7.
2.5. CHANNEL CONTRACTION
2.5.1. Types of contraction

A channel contraction in supercritical flow occurs at locations where
the bottom slope is increased;

the roughness is reduced, or
the discharge decreases.
Usually, the first case is of engineering concern. A well-designed contraction is
characterised by an almost uniform flow in the contracted channel reach. Otherwise,
heavy cross-waves may develop which continue far downstream, and channel walls
of considerable height are needed. Moreover, such flow may entrain much air, and
a stilling basin located at the chute extremity might suffer of asymmetric approach
flow.
The design of a contraction from an approaching channel of width b 1 to a
channel of width b 3 < b 1 is important. Taubert (33, 34] distinguished between the
funnel-shaped, the fan-shaped and the nozzle-shaped contractions (Fig. 12). The first
type is of polygonal shape from the wide to the narrow channel portions. For
fan-shaped contractions, circular arcs curved against the flow are used, while two
circular arcs arranged in an S-shape yield a nozzle-shaped contraction.
Fig. 12
Channel contraction.
Rerrecissement de canal.
a) Funnel. a) Entonnoir.
b) Fan. b) Eventail.
c) Nozzle. c) Tuyere.
37
Le retrecissement en forme d'eventail presente le mm1mum de deviations
elementaires et correspond ainsi a la forme optimale du point de vue geometrique.
En ce qui concerne le retrecissement en forme de tuyere, !'arc amont aura un rayon
plus petit que l'arc aval. Les retrecissements en forme d'entonnoir et en forme
d'eventail seront examines ci-apres.
2.5.2. Retrecissement en forme d'entonnoir

La Fig. 13 represente un retrecissement, de trace rectiligne, avec lignes de
courant et ondes transversales. L'indice « 1 » correspond aux conditions d'ecoule-
ment a l'entree, h 1, Y1 et IF 1 = Y1/(gh 1) 112 etant respectivement la profondeur, la
vitesse et le nombre de Froude de l'ecoulement a l'entree. Les largeurs du canal
d'entree et du canal de sortie sont respectivement b 1 et b 3 • Pour un angle de
retrecissement arbitraire 8, le schema de l'ecoulement est represente sur la Fig. 13.
L'onde transversale positive, ayant comme origine le point A, est reflechie sur lemur
oppose au point C et penetre dans le canal aval. Simultanement, des ondes negatives
se forment au point D. La superposition des ondes positives et negatives conduit a
un ecoulement agite (voir Fig. 10). II en resulte que l'ecoulement dans le canal aval
connait des perturbations.
Fig. 13
Schematic flow pattern in channel contraction according
to lppen and Dawson (1951) [ 16].
Representation schematique de /'ecoulement dans un retrecissement de canal,
d'apres lppen et Dawson (1951) [16}
38
The fan-shaped contraction has the mm1mum sum of elementary wall
deflections, and thus corresponds to the optimum shape from the geometric point
of view. As regards the nozzle-shaped contraction, the radius of the upstream arc
should be chosen much smaller than the second radius. In what follows, the
funnel- and the fan-shaped contractions will be considered.
2.5.2. Funnel-shaped contraction

Fig. 13 shows the straight-lined contraction including streamlines and
cross-waves. Index " I " denotes approaching flow conditions, such that h 1, V 1 and
IF 1 = V/(gh 1) 112 are the approaching flow depth, velocity and Froude number.
Further b 1 and b 3 are the widths in the approaching and tailwater channels. For any
arbitrary contraction angle 8 the flow pattern is as shown in Fig. 13. The positive
cross-wave with origin at point A is reflected at the opposite wall at point C, and
continues into the tailwater. Simultaneously, negative waves are generated at point
D. The superposition of both positive and negative waves leads to tumbling flow as
already shown by Fig. I 0. As a result, the flow in the tail water channel is far from
being without perturbations.
A C,D
G)
®
a) + L1 L2--
L
•h2 ,----
b)
+h1
~-
____ J
Fig. 14
b
Design of funnel-shaped contraction according
to lppen and Dawson (1951) [16].
Conception d'un retrecissement de canal en forme d'entonnoir,
d'apres lppen et Dawson (1951) [16]
a) Plan view. a) Vue en plan.
b) Longitudinal section with surface profile b) Coupe longitudinale avec pro/ii de surface le
along(-) walls and (- - -) channel axis. long des murs ( - ) et dans /'axe (- - -) du
canal.
39
La Fig. 14 represente une nette amelioration de l'ecoulement, resultant d'un
choix correct de !'angle de retrecissement. Si l'onde transversale positive est dirigee
exactement vers le point D, !'interference des ondes se produit comme decrit
precedemment. L'ecoulement qui en resulte dans le canal aval n'a pratiquement pas
de perturbations de surface. II s'ensuit que seul l'ecoulement dans la zone (2) n'est
pas axial; la Fig. 14 b) represente des profils caracteristiques de surface, le long des
murs lateraux et le long de l'axe du canal.
Le cas represente par la Fig. 14 a fait l'objet d'etudes complementaires par
Harrison [l l] et Sturm [32]. D'apres Hager et Bretz [1 O], !'angle de retrecissement 8
peut etre determine par la formule approchee suivante
arctg 8 = -1 -(~ - 1) (2.17)

2 IF1 b3
si 8 < lOo.
Pour de plus grandes valeurs de 8, Jes resultats devront etre verifies vis-a-vis de
!'engorgement, etant donne que IF 3 peut etre proche de l'ecoulement critique. Un
ressaut hydraulique se produit alors dans le retrecissement. De tels phenomenes sont
a ecarter dans tous les cas et on doit assurer une transition progressive de IF 1 a IF3 •
La conception sera basee sur une valeur minimale du nombre de Froude ava1 IF 3 >
2, de fa<;:on que 2 < IF 3 < IF 1.
L'etude se base sur le debit pour lequel la hauteur maximale de l'onde se
produit. Ce debit n'est pas connu a priori. On considere, en general, le debit maximal
comme debit de projet. Pour tous Jes autres debits, la forme du retrecissement ne sera
pas Conforme a )'equation (2.17). Le trace des ondes doit alors etre calcule suivant
Jes methodes presentees par Bagge et Herbich [4] OU dans le chapitre 2.8, en incluant
le frottement sur le mur. Des resultats preliminaires tires d'une etude sur modele
furent presentes par Cotta et Albina [7] pour Jes debits s'ecartant du debit de projet.
Exemple 2.2
Soit Jes caracteristiques de l'ecoulement entrant de l'exemple 2.1 : h 1 = 0,75 m
et IF 1 = 7,37. Quelle est !'angle optimal de convergence 8 si la largeur du canal
d'entree est b 1 = 20 met celle du canal aval b 3 = 15 m?
Avec r- 1 = b 1/b 3 = 1,333 et IF 1 = 7,37, )'equation (2.17) donne arctg 8 = 1,3°.
D'apres Hager et Bretz [10], le rapport des nombres de Froude est IF3 /IF 1 =
(l/2) (l + r) = 0,88, d'ou IF3 = 6,45 > 2. L'engorgement de l'ecoulement est ainsi
empeche.
2.5.3. Retrecissement en forme d'eventail

Le retrecissement examine ci-dessus correspond a des convergences relative-
ment faibles r = b 3 /b 1 > 113 a 1/2. Si le rapport rest plus petit, par exemple de
l'ordre de 1/5, IF 1 doit etre superieur a 10 afin que IF3 > 2. Dans ces cas, la surface
en plan devient excessive. II en resulte que Jes canaux fortement retrecis doivent etre
traites differemment.
Hartung et Knauss [ 12] ont attribue une surface minimale aux retrecissements
de canaux, en forme d'eventail. Ce type d'ouvrage sera, en particulier, con9u pour
Jes evacuateurs de crue dont le seuil est argue. La section d'entree a l'aval du seuil
de l'evacuateur a un rayon R 1 (Fig. 15). Les caracteristiques de ce retrecissement
40
Fig. 14 shows a significant improvement of flow pattern based on a correct
choice of the angle of contraction 8. If the positive cross-wave is exactly directed
to point D, wave interference occurs as previously described. The resulting tailwater
flow is then practically free from surface perturbations. As a consequence, only the
flow in domain (2) is non-axial, and Fig. 14 b) shows typical surface profiles both
along the side walls and the channel axis.
The design case as shown in Fig. 14 was further analysed by Harrison [I I] and
Sturm (32]. According to Hager and Bretz [10) an approximate expression for
contraction angles 8 < I 0° is
arctg 8 (2.17)
For larger 8, the result should be checked against choking as IF3 may be close
to critical flow. Then, a hydraulic jump occurs in the contraction. Such phenomena
must be avoided to provide a continuous transition from IF 1 to IF3 • The design should
be based on a minimum downstream Froude number IF 3 > 2, such that 2 < IF 3 <
IF1.
The design is based on the discharge with the maximum wave height. This
discharge is not known a priori. Usually, the maximum discharge is assumed as the
design discharge. For all other discharges, the contraction shape will not conform
to Eq. (2.17). The wave pattern must then be computed according to approaches
presented by Bagge and Herbich [4) or in section 2.8 by including wall friction.
Preliminary results from a model study for the off-design case were presented by
Cotta and Albina [7].
Example 2.2
Consider the approaching characteristics of example 2.1, that is h 1 = 0.75 m and
IF 1 = 7.37. What is the optimum angle of convergence 0 if the approaching width
b 1 = 20 m and the downstream width b 3 = 15 m are considered?
With r- 1 = b/b3 = l.333 and IF 1 = 7 .37, one may obtain from Eq. (2.17)
arctg 0 = l.3°. According to Hager and Bretz [I OJ, the ratio of Froude number is
IF/IF 1 = (1/2) (1 + r) = 0.88, thus IF 3 = 6.45 > 2. Choking flow is thus inhibited.
2.5.3. Fan-shaped contraction

The preceeding referred to relatively modest contraction ratios r = b 3/b 1 > 1/3
to 112. However, if the ratio r is in the range 115, then IF 1 must be larger than JO in
order that IF3 > 2. In such cases, the area covered in plan view is excessive and the
highly contracting channel must be treated differently.
Hartung and Knauss (I 2] attributed a minimum area to fan-shaped channel

contractions. This type of structure is typically designed for arc-crested spillways.
The entrance section downstream from the spillway crest has a radius R 1 (Fig. 15).
The features of this contraction are a fan-type intake, a contracting channel and a
41
---• .~-rB-
~- ®
tbd
b) 1
Fig. 15
Fan-shaped contraction downstream from a spillway.
Retrecissement en forme d'eventail a /'aval d'un evacuateur de crue.
a) Sectional view. a) Coupe lonfriludinale.

b) Plan view. b) Vue en plan.
(I) Entrance section. (/) Entree.
(2) Tailwater channel. (2) Sortie.
sont : une prise d'eau en forme d'eventail, un canal convergent et un coursier aval
prismatique. Le rapport de convergence est m = bd/bu et !'angle au centre y. Le
radier a l'aval de la zone d'entree est une surface plane et a generalement une pente
longitudinale Constante.
Les conditions d'ecoulement a l'entree (!) sont : la profondeur de l'ecoulement
h 1, la vitesse moyenne V1 et la charge H 1• En supposant un retrecissement entiere-
ment symetrique sans obstacles, tels que des murs guideaux, l'ecoulement tend a se
concentrer au centre de courbure M, la zone DE le long du mur lateral etant privee
d'ecoulement (Fig. 16). Cet effet est du a la composante de la force centrifuge
agissant sur Jes particules se depla9ant le long des lignes de courant courbes. La
profondeur maximale de l'ecoulement hM est situee au point C. Plus a l'aval, le profil
de surface dans !'axe diminue jusqu'au point F, puis atteint un deuxieme maximum
au point H.
L'onde transversale venant du point M intersecte lemur au point E. Sa hauteur
maximale hm se situe au point I, legerement a l'aval de E, et determine la hauteur
du mur lateral. Plus a l'aval, un ecoulement torrentiel avec surface agitee s'etablit.
42
A,
Fig. 16
Flow features in planar fan-type contraction.
Caracteristiques de /'ecoulement dans un ri!trecissement en forme d'eventai/.
b) Axial cross-section. b) Coupe /ongitudinale dans /'axe.
c) Surface profile along the wall. c) Pro/ii de la surface le long du mur.
prismatic tailwater chute. The contraction ratio is m = bd/bu and the sector angle
y. The invert downstream from the entrance section is planar and normally of
constant longitudinal slope.
The flow conditions at the entrance section (1) are flow depth h 1, average
velocity V, and energy head H,. Assuming a fully symmetrical contraction without
any obstacles such as guiding walls, the flow tends to concentrate to the center of
curvature M, associated with a domain of zero flow along the side wall DE (Fig. 16).
This effect must be attributed to the centrifugal force component acting on particles
moving along curved streamlines. The maximum flow depth hM occurs at point C.
Further downstream, the axial surface profile decreases to a minimum at point F and
reaches a second maximum at point H.
The cross-wave from point M intersects the wall at point E. Its maximum hm is
slightly downstream at point I and determines the height of the side wall. Further
downstream the well-known and previously described tumbling surface pattern of
supercritical flow is established.
43
Les deux valeurs relatives YM = hM/H 1 et Ym = hm/H 1 dependent de IF 1 et de
la pente du radier S0 • Pour S0 "" IO% et y "" 40° ( = 0,70 rad), on obtient d'apres
Hartung et Knauss [12] :
YM = 3,66 y (2 + IF1 2)- 112 , l < IF, < 3,5 (2.18)
YM = 0,915 [I - (1,278 - y 2 (2 + IF/)- 112) 112
] IF 1 ~ 3,5 (2.19)
y etant exprimee en radian.
Pour d'autres pentes de radier, on peut ecrire Ym = a YM• OU a < l depend du
rapport de convergence m = bd/bu et du debit a l'entree. Des valeurs de m entre l/3
et 1/2 seront adoptees pour empecher la formation d'ondes transversales dans un
retrecissement en forme d'eventail. Les equations (2.18) et (2.19) peuvent etre
utilisees dans une etude preliminaire. Une etude complete devra, cependant, com-
prendre une· optimisation de la pente du coursier S0 , de !'angle au centre y et du
rapport de retrecissement m. Les resultats de Hartung et Knauss sont preliminaires
et s'appliquent exclusivement a une surface plane de radier. Le coefficient a depend
de plusieurs parametres et une etude sur modele est necessaire pour le projet final
des grands ouvrages.
Le Corps des Ingenieurs de I' Armee des Etats-Unis [35] a mis au point un
evacuateur de crue en forme d'eventail pour le barrage Genegantslet, New York. Le
projet definitif conduisait a une distribution reguliere de l'ecoulement au moyen de
!'augmentation de la longueur des murs d'entree. Les conditions d'ecoulement sur
le coursier etaient egalement ameliorees par un seuil triangulaire avec face amont
verticale, a l'entree du coursier, et par le realignement des murs lateraux. Michels [26]
a presente une etude dans laquelle un seuil avec avant busque est associe a un
coursier qui converge dans une auge.
Une surface de radier plane est une geometrie tres simple du point de vue de
la construction, mais est sujette a des restrictions en ce qui concerne l'hydraulique
de l'ecoulement. Anastasi [3] a etudie un retrecissement de coursier permettant de
supprimer, en grande partie, !es ondes transversales au moyen d'un radier ayant une
courbure transversale. La Fig. 17 represente un retrecissement de coursier, en forme
d'eventail, ayant une largeur b 0 a l'amont et une largeur bct a l'aval. La courbure du
mur a l'amont est I/Ru et la courbure varie le long du retrecissement suivant une
clotho.ide, c'est-a-dire que R- 1 = 0 a l'aval. La Fig. 17 indique egalement les lignes
de courant caracteristiques et !es lignes equipotentielles en forme d'arc circulaire de
longueur b;. ' est !'angle de convergence local et •u !'angle a l'entree.
Anastasi [2, 3] a presente une methode permettant de definir la geometrie de la
surface du radier du retrecissement; cette methode suppose une variation lineaire du
rayon de courbure transversal; au point P :
R/R; = b/n; (2.20)

OU R; R; (x) represente la geometrie de la ligne limite du retrecissement, n; est la
coordonnee transversale et x la coordonnee longitudinale mesuree a partir de la
section aval. En supposant une charge transversale horizontale et une profondeur
d'ecouJement transversale Constante h, on obtient le resultat suivant
y2- S ·sm
-dz = - . (Rn;;)
0 (2.21)
dn g R
44
Both relative values YM = hM/H 1 and Ym = hm/H 1 depend on IF 1, and the
bottom slope S0 • For S0 " " 10 % and y "" 40° ( = 0.70 rad) Hartung and Knauss [ 12]
obtained :
YM = 3.66 y (2 + IF/)- 112 , 1 < IF 1 < 3.5 (2.18)
YM = 0.915 [l - (i.278 - y - 2 (2 + IF/)- 112) 112

] IF 1 ~ 3.5 (2.19)
in which y [rad].
For other bottom slopes Ym may be expressed as Ym = a YM wherein a < I
depends on the contraction ratio m = bd/bu and the approaching discharge.
Typically, m values between 1/3 and I /2 should be considered to prevent the
formation of cross-waves in the fan. Eqs. (2.18) and (2.19) may be suitable for
preliminary design. A definite design should however include an optimization of
chute slope S0 , sector angle y and contraction ratio m. The results of Hartung and
Knauss are preliminary and apply exclusively to a plane bottom geometry. The
coefficient a depends on many parameters and a model study is necessary for the
final design of larger structures.
The US Army, Corps of Engineers [35] developed a fan-shaped spillway for

Genegantslet Dam, New York. The final design involved an even distribution of flow
by increasing the lengths of the approaching walls. The flow conditions in the chute
were also improved by a triangular sill with vertical front face at the chute entrance
and by realigned side walls. A unique design was presented by Michels [26] which
involved a combination of mitre-prowed weir with a convergent chute into a trough.
A plane bottom geometry is quite simple in terms of construction but restrictive

in terms of flow features. Anastasi [3] studied a chute contraction in which
cross-waves were largely suppressed by a transversally curved bottom geometry. In
Fig. 17 a fan-type chute contraction is shown with upstream width bu and
downstream width bct. The wall curvature of the upstream section is 1/Ru, and the
wall curvature along the contraction varies according to a clothoid, that is R 1 = 0
at the downstream section. Fig. 17 shows also typical streamlines and circular
arc-shaped equipotential lines of length b,. The local angle of convergence is T and
Tu is the angle at the entrance section.
Anastasi's [2, 3] approach for the bottom surface of contraction geometry

assumes a linear variation of transverse radius of curvature, that is at point P :
R/R; = b/n; (2.20)

Herein, R; = R; (x) is the geometry of contraction boundary, n, the transverse
coordinate and x the longitudinal coordinate measured from the downstream
section. Assuming a horizontal transverse energy head, and a constant transverse
flow depth h, the transverse bottom slope obtains :
y2 .
-dz = --So·sm (Rn;;) (2.21)
dn g R
45
Fig. 17
Fan-type chute contraction. Notation.
Retrecissement de coursier en forme d'eventail. Notation.
/
I
Fig. 18
View of curved bottom geometry
of fan-type contraction (Anastasi, l 982) [3].
Vue en perspective de la geometrie courbe du radier
d'un retrecissement en forme d'eventail (d'apres Anastasi, 1982) [3]
46
a) b)
Pig. 19
Difrerences in flow surfaces - side view (lop) and plane (bouom).
Di!fere11ces de surface d'eco11/eme111 - en liaut : "" e de cote: cu bas : 1•11e e11 pla11.
a) Plane bouom. 01 Radier plan.
b) Curved bouom. h) Radier courbe.
47
Suivant !'equation (2.21), la pente transversale du radier est egale a la somme
du rapport entre l'acceleration centrifuge et !'acceleration de la pesanteur, et de la
composante de pente transversale due a la pente longitudinale S0 • La Fig. 18 donoe
une configuration caracteristique de radier.
Les vues ·de la Fig. 19 comparent des ouvrages de prise d'eau, en forme
d'eventail, pour deux geometries : a) radier plan, b) radier convexe. Les differences
de surface d'ecoulement sont significatives. Des essais ont demontre que les radiers
convexes amelioraient beaucoup l'ecoulement meme pour des debits differents du
debit de projet. Par rapport a la s uppression des ondes basee sur le principe de
!'interference (voir ci-dessus), le radier convexe peut ainsi etre considere com me une
adaptation souple a un ecoulement meme torrentiel.
Le principe de la compensation de I'acceleration centrifuge par I' acceleration de
la pesanteur fut applique par Vischer [36) aux elargissements, aux retrecissements en
forme d'eventail et aux raccordements. Actuellement, ce moyen de reduire les ondes
transversales est Jimite aux courbes de canal.
2.6. ELARGISSEMENT D'UN CANAL
Les ecoulcments sortant d ' une conduite, d'un evacuateur de crue ou d'un
coursier a pente raide peuvent deboucher dans un canal aval de largeur plus grande.
Fig. 20
Cross-wave formation at abrupt channel expansion.
Fom101io11 d'ondes transversales dans rm elargisseme111 hn1sq11e d'un canal.
48
According to Eq. (2.21) the transverse slope of the bottom is equal to the sum
of the ratio between centrifugal to gravitational accelerations, and the transverse
slope component due to longitudinal slope S Fig. 18 shows a typical bottom
0 •
configuration.
Fig. 19 compares fan-type intake structures for the plane and the convex bottom
geometries. The differences in flow surfaces are significant. Experiments verified
that convex bottoms greatly improve flows even at non-design discharge. Compared
to the interference-type wave suppression as discussed previously in this section, the
convex bottom geometry may thus be regarded as a flexible adaptation to even
supercritical flow.
The principle of compensating the centrifugal acceleration with gravitation was

further applied by Vischer [36] to expansions, funnel-type contractions and
junctions. This method of reducing cross-waves is actually restricted to channel
curves.
2.6. CHANNEL EXPANSION
Flow from a closed conduit, spillway or steep chute may emerge to a

downstream channel with a larger width. If the transition from the narrow to the
a}
-,- - -- ~~-- -- - - - - . --- - - · - - --· ---
b)
Fig. 21
Infinitely expanding channel.
Canal avec elargissement infini.
b) Side view along channel axis. b) Profil le long de /"axe du canal.
49
Si la transition entre le canal etroit et le canal plus large est trop brusque, des ondes
transversales importantes peuvent se former comme le montre la Fig. 20. Par contre,
une transition trop progressive conduit a un ouvrage de grande longueur et couteux.
Rouse et al. [28] ont etudie des elargissements de canal dans lesquels la pente
de frottement est egale a la pente du radier. Pour un canal avec elargissement infini
(Fig. 21 ), on a demontre que la surface libre h (x, y) dependait seulement de la
profondeur h 1, de la largeur relative b 1/h 1et du nombre de Froude IF 1 = Q/(gb~hi) 112
de l'ecoulement entrant. Pour des canaux d'entree relativement larges (h/b 1 ~ I),
la distribution des pressions est presque hydrostatique sauf dans la zone laterale de
l'elargissement. La profondeur relative est donnee par
E_ = f (~. L, IF 1) (2.22)
h1 b1 b1
En outre, les deux parametres independants x/b 1 et IF 1 peuvent etre associes
lineairement sous la forme [(x/b 1)1F 1- 1]. La Fig. 22 donne une representation
graphique generale pour 2 < IF 1 < 8. Rouse et al. [28] ont determine une courbe
limite YL(x) contenant 90 % du debit, pour fe = I, et pouvant etre definie par
('equation suivante :
YL
b1
= !
2
[1 + f (-x-)112]
b1 IF1
e
(2.23)
1.5 -------~------,..-.-----------"T--------,
1 /
0.3
x
0.9 0.8 0.6 0.5 0.4 b11F1
0
0 0.5 1 1.5 2
Fig. 22
Generalized surface plot h/h 1 for 2 < IF, < 8
(Rouse et al., 1951) [28]. (- - -) Eq. (2.23).
Graphique general de la surface h/h 1 pour 2 < IF, < 8
(d'apres Rouse et al., 1951) [28]. (- - -) equation (2.23).
50
wider channel is too abrupt, large cross-waves occur as shown in Fig. 20. On the
other hand, if the transition is too gradual, the design is long and uneconomical.
Rouse et al. [28] considered the design of channel expansions in which the
friction slope equals the bottom slope. For an infinitely expanding channel (Fig. 21 ),
the free surface h (x, y) was found to be a function of the inflow depth h 1, the relative
inflow width b/h 1 and the inflow Froude number IF 1 = Q/(gbrhi) 112 • For relatively
wide approaching channels (h 1/b 1 < 1), the pressure distribution is almost
hydrostatic except for the lateral expansion area. Thus, the relative depth is given
by:
(2.22)
Further the two independent parameters x/b 1 and IF 1 may be combined linearly
into [(x/b 1)1Fj 1]. Fig. 22 shows a generalized contour plot for 2 < IF 1 < 8. Rouse,
et al. [28] determined a limit curve yJx) which enclosed 90 % of discharge,
for fe = I, and may be expressed as :
~~ ~ [I (b ~J ]
2
= + fe
1
(2.23)
y
2 b1
.
... ~ -~-·~-3 //"
-
.. ;,~/o~/ -~ bib1
1 -
-
~--/2
-
~-/
.... ··~--------
~
Fig. 23
Generalized geometry of expansion
with(. .. ) Eq. (2.23) with f, = 114 and (- - -) Eq. (2.24).
Representation geometrique generale d'un elargissement,
suivant /'equation (2.23) avec f.. = //4 (. .. )et /'e4uation (2.24) (- - -).
51
II ya lieu de noter que la Fig. 22 est seulement applicable aux faibles pentes de
radier. Pour Jes fortes pentes, l'ecoulement a tendance a moins s'etendre laterale-
ment.
Une zone de transition depuis le canal d'entree de largeur b 1 jusqu'au canal aval
de largeur b 2 est souvent necessaire. Pour des valeurs finies de b 2/b 1, on peut utiliser
!'equation (2.23) avec fe = 114 pour representer la courbe d'extension limite a
laquelle on doit ajouter une courbe convergente. La Fig. 23 presente une etude
generale. La longueur de la zone de transition L, est determinee comme suit :
(2.24)
Par rapport aux retrecissements de canal, un elargissement de canal conduit a

une geometrie a lignes courbes. D'apres Rouse et al. [28], la conception representee
sur la Fig. 23 donne un ecoulement sortant essentiellement uniforme.
Des etudes sur l'ecoulement torrentiel dans un elargissement ont ete egalement
effectuees par Koch [23] !ors de !'analyse de deux jets se heurtant contre un canal
horizontal, et par Herbich et Walsh [14] avec utilisation d'une methode numerique.
Cependant, leurs observations s'ecartaient beaucoup des resultats des calculs, ce qui
fut attribue au fait d'avoir neglige Jes effets du frottement sur le mur et de la
courbure. L'etude de Rouse peut servir d'approche grossiere dans un projet, bien que
la longueur de la zone de transition L, soit tres grande.
2.7. COURBES
Les courbes d'un canal a ecoulement libre perturbent un ecoulement torrentiel

de la meme fa9on que Jes retrecissements. Le mur courbe exterieur correspond a une
deviation positive et le mur courbe interieur a une deviation negative. Ces deux
deviations conduisent a une surface d'ecoulement complexe, comme l'ont deja
montre Ippen et Knapp [18], et Knapp [22] proposa une methode d'etude. Comme
represente sur la Fig. 24 a), l'ecoulement torrentiel courbe dans le canal est caracte-
rise par un type d'ondes stationnaires. Un couple de niveaux d'ondes, l'un maximal
et l'autre minimal, se manifeste a 0, 2 0, ... a partir de la section d'entree. Chaque
maximum sur le mur courbe exterieur est suivi d'un minimum.
Avec ~o = arcsin CIFo-l) comme angle d'onde elementaire de l'ecoulement a
l'entree, 0 est relie a la courbure relative du canal par !'equation suivante
0 = arctg [ b/R ] (2.25)

(I + b/2R) tg ~
La Fig. 24 b) donne une representation graphique ou !'on constate que 0

augmente quand le nombre de Froude IF0 de l'ecoulement entrant et la courbure
52
Note that Fig. 22 is applicable only to small bottom slopes. For large slopes the
flow has a tendency to expand less laterally.
A transition from the approaching channel of width b 1 to the downstream

channel of width b 2 is often required. For finite values of b 2 /b 1 Eq. (2.23) with fe =
114 can be used as the limit expanding curve, to which a converging curve must be
added. Fig. 23 presents a generalized design. The length of transition L, is :
(2.24)
Compared to channel contractions, a channel expansion involves a curved

expansion geometry. According to Rouse, et al. [28], the design according to Fig. 23
gives an essentially uniform outflow. The approach was not checked for discharges
other than the design discharge.
Supercritical expanding flow were also studied by Koch [23] when analyzing
two impinging jets on a horizontal channel, and by Herbich and Walsh [14] using
a numerical approach. Their observations were in poor agreement with the
computations, however, and the deviation was attributed to the neglect of wall
friction and curvature effects. Presently, the design of Rouse can be taken as a rough
design guide although the length of transition L, is extremely long.
2.7. BENDS
Open channel bends perturb a supercritical flow similarly to contractions. The

outer bend wall corresponds to a positive deflection, and the inner bend wall to a
negative deflection. Both yield a complex surface pattern, as described by lppen and
Knapp [18] and Knapp [22] proposed a design procedure. Curved supercritical
channel flow is characterized by a standing wave pattern (Fig. 24 a)). A pair of
maximum and minimum wave elevations occurs at locations 8, 2 8, ... from the
inflow section. Each maximum at the outer bend wall is followed by a minimum.
The elementary wave angle of the approaching flow P0

= arcsin (IF0 - 1
) is
related to the angle 8 and the relative channel curvature as
e = arctg [ b/R ] (2.25)

(1 + b/2R) tg p
According to Fig. 24 b) the angle 8 is seen to increase with both increasing

inflow Froude number IF0 and relative curvature b/R. The domain of small IF0
53
IF- 1
o --
\
\
\
b/R "" " \
\
0.2
"2", CO\
\
\ \
\ \
R e \ \
\ I
t -----
IF-1
0
0
a) b) 0 0.2 0.4 0.6 0.8 1
Fig. 24
Supercritical flow in bend.
Ecoulement torrentiel dans une courbe.
a) Schematic flow pattern with • wave a) Representation schematique de /'ecoulement :
maximum and o wave minimum. • onde maximale, o onde minimale.
b) Angle 0 as a function of IF" 1
and b/R (Sinni- b) Angle e enfonction de IF;, I et de b/R (d'apres
ger and Hager, 1989) [31]. Sinniger et Hager, 1989) [31]
relative b/R augmentent. La zone des faibles valeurs de IF0 (lignes-tirets de la

Fig. 24 b)) est une zone ou se forment des ressauts ondules. Pour les grandes valeurs
de IF0 et b/R < I /2, I' equation (2.25) peut s'ecrire sous une forme approchee
(2.26)
Les hauteurs minimale et maximale (indice « e ») des on des peuvent etre

calculees a partir des equations (2.16) et (2.25). Dans la Fig. 25, Ye= ± he/h est 0
represente graphiquement en fonction du nombre de Froude IF0 de l'ecoulement a

l'entree et de la courbure relative b/R de la courbe. Si (b/2R)IF 2 < I, Jes hauteurs 0
des ondes extremes peuvent s'exprimer sous la forme suivante tiree de l'eq. (2.26) :
(2.27)
Suivant I' equation (2.27), Ye est toujours du cote prudent et ne sera applique que
pour Jes faibles valeurs de b/R < 0,5.
54
3
Ye
·b/R o.os 0.2
2
0
0 0.2 0.4 0.6 0.8 1
Fig. 25
Extreme wave heights y, as a function of IF,, and b/R
(Sinniger and Hager, 1989) [31].
Hauteurs des ondes extremes y,. en fonction de IF,, et de b/R
(d'apres Sinniger et Hager, 1989) {Jlj
shown dashed in Fig. 24 b)) is a region of undular jumps. For large IF and 0
b/R < I 12, Eq. (2.25) may be approximated as :
(2.26)
The minimum and maximum (index « e »)wave heights can be computed with
Eqs. (2.16) and (2.25). In Fig. 25, Ye = ± h/h 0 are plotted as a function of
approaching Froude number IF and relative bend curvature b/R. Provided
0
(b/2R)IF 2 < I, the extreme wave heights Ye may be approximated based on

0
Eq. (2.26) as :
(2.27)
Ye according to Eq. (2.27) is always on the safe side and should be applied only
for small values of b/R < 0.5.
55
Exemple 2.3
Soit un canal de largeur b = 15 m, de profondeur d'ecoulement h = 0,75 met
de nombre de Froude IF = 7,0. Quelles sont les caracteristiques de J'ecoulement
0
dans une courbe de rayon de courbure R = 75 m et d'angle de deviation i; = 70°?

Avec ~ 0 = arcsin (IF 1) = 8,2° et b/R = 0,2, ['equation (2.25) donne
0
-
8 = 51,5°; d'ou ~ ± 112 8 = 8,2° ± 25,8°. Pour l'onde maximale, on a Ymax = 9,28
d'apres )'equation (2.16); la profondeur maximale de l'ecoulement est done dix fois
environ superieure a la profondeur de I'ecoulement a J'entree. A l'evidence, ce projet
est mauvais et on doit rechercher des moyens pour reduire J'onde maximale.
2.8. REDUCTION DES ONDES AV AL
D'apres I' equation (2.27), les extremes des ondes augmentent de fa<;:on significa-
tive avec le nombre de Froude IF de I'ecoulement entrant; Jes ondes peuvent aussi
0
continuer a se manifester sur un long parcours dans le canal aval. Cela n'est pas
souhaitable des points de vue hydraulique et economique. Knapp [22] a propose un
certain nombre de moyens permettant de reduire les ondes dans des canaux courbes.
Ces moyens comprennent :
- Relevement de la courbe, en donnant au radier une pente transversale S,
pour s'opposer a la pente centrifuge. D'apres Jes principes de la mecanique, on a
St -- v~ (2.28)
gR
Le relevement doit etre progressif, en partant d'une pente transversale nulle au
debut de la courbe pour atteindre la valeur definitive a )'emplacement de la premiere
onde maximale. Le relevement presente un serieux inconvenient, a savoir que S,
depend de la vitesse V0 a I'entree et done du debit Q. Cependant, par rapport a
S, = 0, son effet pour un debit s'ecartant du debit de projet est une amelioration de
J'ecoulement.
- Ailettes courbes multiples, en inserant des elements concentriques dans Ia
courbe, ce qui conduit a une courbe plus etroite.
- Interference, comme cela est applique dans le cas des retrecissements.
Parmi ces trois methodes, le relevement s'adapte le mieux aux canaux ou la
vitesse d'entree ne varie pas beaucoup. Afin que le niveau maximal de la surface
reste bas, le mur interieur est abaisse (Fig. 26), ce qui peut conduire a un cout
d'excavation supplementaire. L'utilisation d'ailettes est restreinte a cause des trans-
ports solides. Egalement, Jes ailettes perturberont I'ecoulement et creeront une
nouvelle configuration d'ondes transversales. Un grave inconvenient des ailettes et
des seuils est le risque de cavitation pour des vitesses relativement grandes. Des
ailettes peuvent etre incorporees comme des piles de pont. Les courbes en spirale ne
sont pas recommandees du fait qu'elles n'ameliorent pas l'ecoulement. L'utilisation
de canaux non rectangulaires sera egalement evitee si possible. Les resultats
presentes s'appliquent a des nombres de Froude superieurs a 1,5 afin d'exclure Jes
ressorts ondules.
56
Example 2.3
Given a channel of width b = 15 m, flow depth h = 0.75 m and Froude number
IF = 7.0. What are the flow features in a curve of which the radius is R = 75 m and
0
the angle of deviation <; = 70°?

With ~ 0 = arcsin (IF0 - 1 ) = 8.2° and b/R = 0.2, Eq. (2.25) yields 8 = 51.50,
and thus ~ ± 112 8 = 8.2° ± 25.8°. The wave maximum thus obtains Ymax = 9.28
according to Eq. (2.16). The maximum flow depth is thus almost ten times as large
as the approaching depth. Clearly, this design would be poor, and means to reduce
the wave maximum must be sought.
2.8. REDUCTION OF TAILWATER WAVES
According to Eq. (2.27) the wave extrema increase significantly with the
approaching Froude number IF Also, the waves may persist over a long reach in the
0
•
the tail water channel. This feature is undesirable both from hydraulic and economic
points of view. Knapp [22] suggested a number of means by which the wave pattern
in curved channels may be reduced. These include :
- Curve banking by providing a transverse bottom slope S, to counter the

centrifugal slope. According to basic mechanics, the slope is given by
S t -- v; (2.28)
gR
The banking must be gradual, starting from zero transverse slope at the
beginning of curve to the final value at about the location of the first wave maximum.
A serious disadvantage of banking is that the slope S, depends on the approach
velocity V0 , and therefore on the discharge Q. Yet, its effect for off-design discharge
is improving the flow as compared to S, = 0.
- Multiple curved vanes, by inserting concentric elements in the curve and

thus narrowing a vaned curves.
- Interference treatment as was applied to contractions.
Of these three methods, banking is most adaptable to channels where the
approaching velocity does not change significantly. In order that the maximum
surface elevation remains low, the inner wall is depressed (Fig. 26). This may involve
extra excavation cost. The use of vanes is restricted because of debris. Also, vanes
will disturb flow and create a new pattern of cross-waves. A serious disadvantage of
vanes and sills is the danger of cavitation for relatively high velocities. Vanes may
be incorporated as bridge piers. Spiral curves are not recommended since they do
not improve the flow pattern. Also, the use of nonrectangular channels should be
avoided whenever possible. The results presented apply to Froude numbers larger
than 1.5 to exclude undular jumps.
57
Fig. 26
Curve banking at the " Dechargeur de Saint-Chamas ", Q = 60 m 3/s
(La Houille Blanche 3, 1966).
Relevement d'une courbe au dechargeur de Saint-Chamas, Q = 60 m 3/s
(La Houille Blanche 3, 1966).
D'autres etudes sur Jes ecoulements torrentiels dans des canaux courbes ont ete
realisees par Poggi [27J qui a compare Jes methodes de von Karman et de Ippen-
Knapp a ses experiences. Les conclusions furent les suivantes :
- Pour des courbes de faible courbure relative (b/R < 0,04) et de faible pente
de radier (S 0 < 5 %), les deux methodes indiquent de fa~on acceptable Jes ondes
extremes.
- Pour des courbures plus grandes et des pentes de radier plus fortes, la
methode de Ippen-Knapp peut etre recommandee. Ainsi, I' equation (2.16) au lieu de
!'equation (2.13) sera utilisee.
L'influence de la pente S du radier a ete analysee par Sananes et Acatey [29].
0
En se basant sur les experiences de Poggi [27], ils trouverent que la methode
d'Ippen-Knapp n'etait en accord avec les observations que si I'ecoulement le long
de la courbe etait accelere, c'est-a-dire lorsque la pente le long de la courbe etait au
moins deux fois plus grande que la pente S a I'entree. Le premier maximum d'onde
0
depasse egalement tous Jes autres maxima. Ces resultats obtenus egalement par
Krause [24] soutiennent grandement le concept de reduction des ondes par releve-
ment de la courbe.
58
rh
I I f
~ __111 _,;:;,;:;/ 3.0
=-- I 1---
l,I k
---7.5 I 0
t' -~-+~
1- 4.2-J 1@
2
Fig. 27
Combination of contraction and bend in a channel
with supercritical flow (Hartung and Knauss, 1967) [12].
Association d'un retrecissement et d'une courbe dans un canal
a ecoulement 10rrentiel (d'apres Hartung et Knauss, 1967) [12].
a) With dividing wall. a) Avec mur separateur.

b) Without dividing wall. b) Sans mur separateur.
W, ® Transverse surface profiles at sections 5 W, ® Coupes transversa/es de /'ecou/ement
and 9. aux sections 5 et 9.
Further studies on supercritical flows in curved channels were due to Poggi [27]
who compared the approaches of both von Karman and Ippen-Knapp with his
experiments. He concluded that :
- For curves with small relative curvature (b/R < 0.04) and small bottom
slope (S 0 < 5 %) both approaches reasonably predict the wave extrema.
- For larger curvature and steeper bottom slope the approach of lppen-Knapp
as presented herein may be recommended. Thus, Eq. (2.16) instead of Eq. (2. 13)
should be used.
The effect of bottom slope S0 was analysed by Sananes and Acatey [29]. Based
on the experiments of Poggi [27], they found that the approach by lppen-Knapp is
in agreement with observations only if the flow along the curve is accelerated, that
is when the slope along the curve is at least two times as large as the approaching
slope S0 • Also, the first wave maximum is higher than all other maxima. These
findings which were also stated by Krause [24] greatly favour the concept of wave
reduction by curve banking.
59
Harrisson [11] a conduit des experiences sur Jes ecoulements a grande vitesse
dans des canaux de diverses formes. II constata des effets importants de la vitesse
et du rayon de la courbe sur la formation des ondes. C'est dans le canal rectangulaire
que Jes ondes sont le moins actives.
Hartung [13] a examine Jes aspects pratiques des ecoulements dans des coursiers
et a recommande des courbes equipees d'ailettes comme dispositif pour reduire Jes
concentrations d'ecoulement. La Fig. 27 donne un exemple d'association d'un
retrecissement avec une courbe dans un canal, tel que l'ont propose Hartung et
Knauss [12]. Le canal est divise en deux parties par un mur separateur, ce qui
permet de reduire de moitie l'effet de la force centrifuge. L'extremite aval du mur
separateur peut etre determinee suivant le principe des interferences d'ondes. Un tel
dispositif se limite a un certain domaine de debits au-dessous du debit de projet. Le
mur separateur peut aussi etre utilise comme support pour la couverture du canal.
La conception optimale doit etre etudiee sur modele hydraulique.
Un calcul de la surface d'un ecoulement torrentiel dans une courbe d'un canal
trapezoidal, Jes cotes lateraux ayant une pente m, a ete effectue par Lenau [25]. Le
cas h/R ~ I et IF0 > 2 fut considere, pour lequel Jes cretes des ondes deviennent
presque periodiques, comme cela est suppose dans la methode classique. La hauteur
maximale d'onde fut determinee par la formule suivante :
Ym = I + IF} (~ + 2~ho) (2.29)
qui est analogue a l'equation (2.27) pour de faibles valeurs de b/R. La Fig. 28
presente des photographies d'ecoulement dans une courbe d'un canal trapezoidal,
d'apres Kleinschroth et Abelein [20]. Ils ont recommande :
- des courbes de transition, a double rayon de courbure, a l'amont et a l'aval
du coude,
un mur lateral exterieur raide (pente horiz./vert. = 0,5),
- une surface de radier rugueuse,
afin d'ameliorer l'ecoulement. L'effet de la pente transversale du radier sur la
reduction des ondes stationnaires etait faible et un canal rectangulaire devrait etre
normalement adopte. Une etude sur modele des courbes dans un canal trapezoidal
fut conduite par Kleinschroth et Siegerstetter [21]. L'equation (2.29) fut verifiee et
une revanche supplementaire se justifia pour l'entrainement d'air et Jes projections
d'eau.
Ellis et Pender [8] ont propose une solution au probleme des ecoulements
torrentiels, en se basant sur Jes observations de modeles. Ils ont resolu Jes equations
bidimensionnelles de l'ecoulement torrentiel permanent. Le programme de calcul du
systeme d'equations caracteristiques a donne des resultats satisfaisants en tenant
compte des effets de la pente, du frottement et d'autres elements relatifs a la
dissipation de l'energie. Les calculs furent appliques a des courbes simples, des
courbes en forme de S et des elargissements. II est reconnu que de tels calculs ne
peuvent remplacer les etudes sur modeles, mais on peut reduire la duree et le cout
des etudes sur modeles physiques en effectuant des etudes numeriques.
Un important developpement du calcul bidimensionnel des ecoulements a
surface libre est du a Jimenez et a Chaudhry [19], qui ont utilise des programmes de
calculs aux differences finies, explicites, concernant les ondes de choc. Si l'on tient
60
Harrisson [11] made experiments on high velocity flows in channels of various
shapes. He found appreciable effects of both velocity and curve radius on wave
generation. The lowest wave action occurred in the rectangular channel.
Hartung [13] discussed practical aspects of chute flows and recommended

vaned curves as a means to reduce flow concentrations. Fig. 27 shows an example
of combining a channel contraction with a channel bend due to Hartung and
Knauss [12]. The vane divides the channel into two parts and reduces the effect of
centrifugal force in half. The downstream end of vane may be determined according
to the principle of wave interference. The use of vanes is then limited to a certain
percentage of discharge below the design discharge. The dividing wall may also be
the foundation for a channel cover. The optimum design must be sought with a
hydraulic model.
Lenau [25] computed the surface pattern of supercritical flow in a trapezoidal

channel bend of side slope m. The case h/R ~ 1 and IF > 2 was considered for
0
which the pattern of crests becomes almost periodic, as is assumed in the

conventional approach. The maximum wave height was predicted as :
Ym = J + 1Fo2 (* + 2';ho) (2.29)
analogous to Eq. (2.27) for small b/R. Fig. 28 shows photographs of bend flow in
a trapezoidal channel according to Kleinschroth and Abelein [20]. They
recommended :
- transition curves of double radius of curvature both up- and downstream
from the bend,
- steep outer side wall of slope 0.5 (horizontal) to I (vertical), and
- rough outer side wall,
to improve the flow. The effect of transverse bottom slope on the reduction of the
standing wave pattern was small and a rectangular channel should be normally used.
A model study on bends in a trapezoidal channel was conducted by Kleinschroth
and Siegerstetter [21]. Eq. (2.29) was verified and an additional freeboard should
account for air entrainment and splash water.
Ellis and Pender [8] introduced an alternative to the design of supercritical flows
based on model observations. They solved the two-dimensional equations of steady
supercritical flow. The computational scheme for the system of characteristic
equations was effective in including effects of slope, friction and other energy
dissipating elements. The computations were applied to simple bends, S-shaped
bends, and expansions. Although these computations cannot replace model studies,
the time and cost of physical modelling could be reduced by performing numerical
studies.
A significant development of the computation of 20 free surface flows was

achieved by Jimenez and Chaudhry [19] using explicit, shock-capturing
finite-difference schemes. Provided the boundary conditions are properly accounted
61
a)
b)
Fig. 28
Bend flow in trapezoidal channel according
to Kleinschroth and Abelein ( 1985) [20].
Ecoulement dans une courbe d'un canal trapezofdal
(d'apres Kleinschroth et Abelein, 1985) [20}.
a) Smooth outer side wall. a) Mur lateral exterieur /isse
b) Rough outer side wall. b) Mur lateral exterieur rugueux.
62
2.5
2.0
1.5
1.0
0.5
9. 87
Y/Ho
a) 2. 5 0. 00
NOTE CHANNEL WALL::i

H/H() NOT SHOWN
2 10
b)
X/HO -------- 59.61
-10
Fig. 29
Supercritical free-surface flow, flow direction left to right
(Jimenez and Chaudhry, 1988) (19].
Ecoulement torrentiel a surface fibre, direction de /'ecoulement
de la gauche vers la droite (d'apres Jimenez et Chaudhry, 1988) {19]
a) Computation of shock-reflection. a) Ca/cul de la reflexion de /'onde de choc.
b) 3 D-plot of water surface in circular-arc b) Representation tridimensionnelle de la surface
contraction - both for IF 4 and de /'ecoulement dans un retrecissement
Maccormack number v = 0.98. presentant une courbe circulaire - pour IF = 4
et une valeur v (nombre de MacCormack)
0,98.
63
compte correctement des conditions aux limites, on peut obtenir des resultats
satisfaisants. On a constate que l'hypothese de repartition hydrostatique des pres-
sions imposait des restrictions pour Jes ecoulements presque critiques.
La Fig. 29 donne deux solutions caracteristiques : calcul de la reflexion de
l'onde de choc et surface de l'ecoulement dans un retrecissement ayant une courbe
circulaire. Le rapport presente egalement en detail Jes techniques de calculs des
ecoulements a surface libre, dans le domaine bidimensionnel.
Actuellement, Jes modeles mathematiques relatifs aux ecoulements torrentiels
dans Jes canaux evoluent rapidement. Ence qui concerne l'ecriture des programmes,
ii semble important que des systemes explicites, avec une operation courante de
prevision-correction, fonctionnent avec succes. La definition des conditions aux
limites, la modelisation des sommets des ondes relativement raides et des effets du
frottement, posent encore quelques problemes.
Un excellent compte rendu des calculs numeriques des ecoulements a surface
libre est donne par Chaudhry [5].
2.9. NOTATIONS ADOPTEES DANS LE CHAPITRE 2
A [m2] surface de la section transversale.

a [-) rapport des profondeurs maximales d'ecoulement.
b [m] largeur du canal rectangulaire.
Bs [m] largeur en surface.
~ [-] angle de choc
~o [-] angle de choc elementaire.
c [ms- 1] celerite de l'onde.
0 [-] angle au centre du secteur.
IF [-] nombre de Froude.
8 [-) angle de deviation du mur.
g [ms- 2] acceleration de la pesanteur.
h [m] profondeur de l'ecoulement.
he [m] hauteur de l'onde extreme.
hm [m] profondeur maximale de l'ecoulement au mur.
hM [m] profondeur maximale de l'ecoulement au centre.
L [m] longueur.
L, [m] longueur de la zone de transition.
m [-) rapport de retrecissement (de convergence).
mh/R [-] caracteristiques du coursier trapezoidal.
n [sm-113] coefficient de rugosite de Manning.
ni [m] coordonnee curviligne transversale.
q [m2s-1] debit par unite de largeur.
R [m] rayon de courbure.
r [-] rapport de largeurs.
so [-] pente du radier (du fond).
64
for, favorable results may be obtained. The assumption of hydrostatic pressure
distribution was found to impose restrictions for nearly critical flows.
Fig. 29 shows two typical solutions which include the computation of

shock-reflexion and the water surface in a circular-arc contraction. The paper also
contains a thorough state of the art on computational methods for 20 free surface
flows.
Actually, the numerical modelling of supercritical channel flow is rapidly
evolving. It is important from the point of view of program writing that explicit
schemes with a predictor-corrector routine seem to work successfully. There are still
some problems with the definition of the boundary conditions, the modelling of
relatively sharp wave crests and effects of friction.
An excellent account on the numerical computation of free surface flows is

given by Chaudhry [5].
2.9 NOTATION CHAPTER 2
A [m2] cross-sectional area.

a [-] ratio of maximum flow depths.
b [m] width of rectangular channel.
Bs [m] surface width.
~ [-] shock angle.
~o [-] elementary shock angle.
c [ms- 1] celerity.
0 [-] sector angle.
IF [-1 Froude number.
e [-] wall deflection angle.
g [ms- 2] gravitational acceleration.
h [m] flow depth.
he [m] extreme wave height.
hm [m] maximum wall flow depth.
hM [m] maximum center flow depth.
L [m] length.
L, [m] transition length.
m [-] contraction ratio.
mh/R [-] characteristics of trapezoidal chute
n [sm-113] Manning roughness coefficient.
ni [m] transverse curved coordinate.
q [m2s-1] discharge per unit width.
R [m] radius of curvature.
r [-] width ratio.
so [-] bottom slope.
65
st (-] pente transversale du radier.
't (-] angle de convergence.
v [ms- 1] vitesse dans une section transversale.
x [m] coordonnee longitudinale.
c; [-] angle de la courbe.
y (-] rapport de profondeurs d'ecoulement.
y [m] coordonnee cartesienne transversale.
Ye [-] valeur relative de he.
z [m] cote du radier (du fond).
Indices
ecoulement a l'entree (amont).
2 ecoulement a la sortie (aval).
d a val.
u amont.
2.10. REFERENCES
I. ABBOTI M. B. ( 1966) : "Introduction to the Methods of Characteristics ".American

Elsevier, New York.
2. ANASTASI G. (1980) : " L'Eliminazione delle onde d'urto negli scivoli a
contrazione ".XVII Convegno di Idrau/ica e Costruzioni Idrauliche(Palermo) A 2,
1-12.
3. ANASTASI G. (1982) Die Aufhebung von StoBwellen in
SchuBrinnenverengungen ". Mitteilung No. 59, Versuchsanstalt filr Wasserbau,
Hydrologie und Glaziologie, ETH Ztirich, ed. D. Vischer, Zurich.
4. BAGGE G. and HERBICH J.B. (1967): "Transitions in supercritical open-channel
flow". Journal of Hydraulics Division ASCE, 93 (5), 23-41; 1968, 94 (3), 803-804;
1969, 95 (I), 453-454.
5. CHAUDHRY M. H. (1992): "Computational Open Channel Flow". Van Nostrand
Reinhold, New York.
6. CHOW V. T. (1959). "Open channel hydraulics ".McGraw-Hill, New York.
7. COTIA R. D. and ALBINA H. C. (1972): "Contribucion al diseno de transiciones
para canal es con regimen supercritico ". V. Congreso Latinamericano de
Hidraulica (Lima) A 12, 1-8.
8. ELLIS J. and PENDER G. (1982) : " Chute spillway design calculations ". Proc.
Institution Civil Engineers (London) 73, 299-312.
9. HAGER W. H. (1989) : " Supercritical flow in channel junctions ". Journal of
Hydraulic Engineering ASCE, 115 (5), 595-616.
IO. HAGER W. H. and BRETZ N. V. (1987) : Discussion to " Simplified design of
contractions in supercritical flow". Journal of Hydraulic Engineering ASCE, 113
(3), 422-425.
66
st [-] transverse bottom slope.
't [-] convergence angle.
v [ms-'] cross-sectional velocity.
x [m] longitudinal coordinate.
~ [-] bend angle.
y [-] ratio of flow depths.
y [m] transverse cartesian coordinate.
Ye [-] relative value of he.
z [m] bottom elevation.
Suscripts
approaching flow.
2 tailwater flow.
d downstream.
u upstream.
11. HARRISON A. J.M. (1966) : " Design of channels for supercritical flow ". Proc.
Institution Civil Engineers (London) 35, 475-490; 36, 557-565.
12. HARTUNG F. and KNAUSS J. (1967) : " Developments to improve economy,
capacity and efficiency of structures controlling the passage of flood water
through reservoirs ". International Commission on Large Dams ICOLD, IX
Congress (Istanbul); Q. 33, R. 14, 227-249.
13. HARTUNG F. (1972) : " Gestaltung von Hochwasserentlastungsanlagen bei
Talsperrendammen ". Wasserwirtschaft 62 (1/2), 39-51.
14. HERBICH J.B. and WALSH P. (1972) : " Supercritical flow in rectangular
expansions ". Journal of Hydraulics Division ASCE, 98 (9), 1691-1700; 99 (8),
1276-1277.
15. IPPEN A. T. (1951) : "Mechanics of supercritical flow - a symposium". Part 1,
Trans. ASCE 116, 268-295.
16. IPPEN A. T. and DAWSON J. H. (1951):" Design of channel contractions". Trans.
ASCE 116, 326-346.
17. IPPEN A. T. and HARLEMAN D.R. F. (1956): "Verification of theory for oblique
standing waves ". Trans. ASCE 121, 678-694.
18. IPPEN A. T. and KNAPP R. T. (1936) : "A study of high-velocity flow in curved
channels of rectangular cross-section ". Trans. American Geophysical Union VII
Annual Meeting, 516-521.
19. JIMENEZ 0. F. and CHAUDHRY M. H. (1988) : " Computation of supercritical
free-surface flows ". Journal of Hydraulic Engineering ASCE, 114 (4), 377-395.
20. KLEINSCHROTH A. and ABELEIN R. (1985) : " Modelluntersuchungen iiber die
Wasserspiegellage in der Kriimmung einer Schu13rinne ". Mitteilung No. 44, Inst.
Hydraulik und Gewasserkunde, TU Munchen, P.-G. Franke ed., Munchen.
67
21. KLEINSCHROTH A. and SIEGERSTETTER L.A. (1986) : " Investigations of
supercritical flow in a curved section of an extremely rough channel ". XX
Convegno di Idraulica e Costruzioni Idrauliche (Padova) 873-884.
22. KNAPP R. T. (1951) : "Design of channel curves for supercritical flow". Trans.
ASCE, 116, 296-325.
23. KocH K. (1968) : " Die gegenseitige Strahlablenkung auf horizontaler Soble ".
Bericht Nr. 15, Versuchsanstalt fiir Wasserbau, Oscar v. Miller lnstitut, ed. F.
Hartung, Munich/Obernach.
24. KRAUSE D. ( 1970) : " EinfluB der Trassierungselemente auf den Spiegelverlauf
in gekriimmten SchuBrinnen ". Technischer Bericht No. 6, Institut Hydraulik und
Hydrologie, TH Darmstadt, R. Schroder ed., Darmstadt.
25. LENAU C. W. (1979) : " Supercritical flow in bends of trapezoidal section ".
Journal of Engineering Mechanics Division ASCE, 105 (EM I), 43-54.
26. MICHELS V. (1985) : " Rosslynne spillway : mitre-prowed overfall weir with
sharply convergent walls and skew-impact trough apron". XX/ /AHR Congress
(Melbourne) 3, 628-632.
27. POGGI B. (1956): "Correnti veloce nei canali in curva ". L'Energia Elettrica 34,
465-480.
28. ROUSE H., BHOOTHA, B. V. and Hsu E. Y. (1951) : " Design of channel
expansions ". Trans. ASCE 116, 1369-1385.
29. SANANES F. and ACATAY T. (1962) : " Contribution a l'etude des ecoulements
68
supercritiques dans des canaux coudes ". Comptes Rendus de l'Academie des
Sciences (Paris), 255, 1-3.
30. SANANES F. and XANTHOPOULOS T. (1965) : " Amelioration des conditions
d'ecoulement dans un canal courbe de section rectangulaire a forte pente ''.
Comptes Rendus de l'Academie des Sciences (Paris) 260, 801-803.
31. SINNIGER R. 0. and HAGER W. H. (1989) : " Constructions Hydrauliques -
Ecoulements Stationnaires ". Presses Polytechniques Romandes, Lausanne.
32. STURM T. W. (1984) : "Simplified design of contractions in supercritical flow".
Journal of Hydraulic Engineering ASCE, 111 (5), 871-875.
33. TAUBERT U. (1971) : " Der Abflul3 in Schu13rinnen-Verengungen ''. Der
Bauingenieur 46 (11), 385-392.
34. TAUBERT U. (1974) : " The design of spillway contractions using computer
simulation ". Water Power 26 (8), 282-287.
35. United States Army, Corps of Engineers ( 1952) : " Spillway for Genegantslet
dam Genegantslet creek, New York". Technical Memorandum 2-351. Waterways
Experiment Station, Vicksburg.
36. VISCHER D. L. (1988) : "A design principle to avoid shockwaves in chutes". The
International Symposium on Hydraulics for High Dams (Beijing), 391-396.
37. VON KARMAN T. (1938) : "Eine praktische Anwendung der Analogie zwischen
Oberschallstromung in Gasen und iiberkritischer Stromung in offenen
Gerinnen ''. Zeitschrift far Angewandte Mathematik und Mechanik 18, 49-56.
69
3. ENTRAINEMENT D'AIR NATUREL
3.1. INTRODUCTION
Le phenomene des ecoulements aeres revet un grand inten~t dans le projet des
canaux a forte pente, appeles aussi coursiers. L'entrainement d'air dans J'ecoulement
conduit a un melange air-eau et le gonflement de J'ecoulement aere necessite des
murs lateraux plus hauts que dans le cas d'un ecoulement non aere. L'inconvenient
du gonflement est compense par une suppression des degats par cavitation.
L'entrainement d'air dans un ecoulement sur un coursier concerne toujours un
ecoulement a grande vitesse, sous reserve que la rugosite aux parois presente les
caracteristiques habituelles. L'entrainement d'air dans un ecoulement sur coursier
est qualifie de nature), etant donne qu'aucun appareil externe, tel qu'un ventilateur,
n'est necessaire pour l'alimentation en air. Le processus d'entrainement d'air est lie
aux caracteristiques de turbulence. La presence d'air sous forme de bulles dans le
melange air-eau influe sur la turbulence. Comme Jes bulles montent vers la surface
libre, le melange est influence par Jes composantes longitudinale et transversale.
L'entrainement d'air a la surface sur Jes coursiers peut s'expliquer en se
reportant a la Fig. 30. A la sortie de la retenue, sur la crete de l'evacuateur de crue,
J'ecoulement s'accelere. Lorsque le niveau de turbulence a la surface libre est
suffisamment eleve, de l'air commence a etre entraine. Dans la zone d'acceleration,
l'ecoulement est designe sous le nom d'ecoulement non uniforme. Si la geometrie et
la rugosite du canal ne varient pas a I'aval, un ecoulement aere uniforme peut
s'etablir, ressemblant beaucoup a un ecoulement uniforme classique. Dans cette
zone, Jes distributions de la vitesse et de la teneur en air ne varient pas. Cet etat
d'equiJibre peut etre perturbe par des Variations de Ja pente du radier OU de )a Jargeur
du canal. II en resulte qu'un ecoulement aere graduellement varie se forme jusqu'a
ce qu'un etat d'equilibre soit retabli. L'ouvrage de sortie, tel que le saut de ski sur
la Fig. 30, ne sera pas decrit ici. Sauf indication contraire, tous Jes resultats presentes
concernent un coursier rectiligne.
---@---
70
3. NATURAL AIR ENTRAINMENT
3.1. INTRODUCTION
Aerated flow is a significant phenomenon in the design of steep channels, also

referred to as chutes. Air entrainment in the water flow leads to an air-water mixture
and the bulking of aerated flow requires higher side walls than for non-aerated flow.
The disadvantage of bulking is balanced by a suppression of cavitation damage.
Air entrainment in chute flow is always related to high velocity flow, provided
the usual boundary roughness is considered. The air entrainment in chute flow is
said to be natural, as no external mechanism such as a fan is needed for the air
supply. The process of air entrainment is associated with the turbulence
characteristics. The presence of air in the form of bubbles in the air-water mixture
affects the turbulent structure. As the bubbles rise to the free surface, the mixture is
influenced by a combination of longitudinal and transverse components.
The surface air entrainment in chutes can be understood by referring to Fig. 30.
As water leaves the reservoir over the spillway crest, the flow accelerates. When the
turbulence level at the free surface is sufficiently high, air begins to be entrained. In
the accelerating region, the flow is referred to as non-uniform. If the channel
geometry and roughness do not vary further downstream, a uniform aerated flow
may be established, much in the manner of conventional uniform flow. In this
domain, both the velocity and air concentration distributions remain invariant. This
state of equilibrium may be perturbed by changes of bottom slope, or channel width.
As a result, gradually varied aerated flow occurs until a state of equilibrium is
re-established. The terminal structure such as the ski jump in Fig. 30 will not be
described herein. All results to be developed refer to a straight chute unless otherwise
stated.
Fig. 30
Schematic spillway with reservoir (left)
and ski jump (right).
Representation schematique d'un evacuateur de crue avec retenue (a gauche)
et saut de ski (a droite).
@ Regions of uniform flow. @ Zones d'ecoulement uniforme.
@ Regions of non-uniform flow. @ Zones d'ecoulement non uniforme.
@ Aerated flow (right). 0 Ecoulement aere (a droite).
(vj) Non-aerated flow (downstream from reser- (vj) Ecoulement non aere (a l'aval de la retenue).
voir).
71
Les paragraphes 3.2 et 3.3 traitent respectivement du mecanisme de l'entraine-
ment d'air et de !'emplacement de la formation de l'entrainement d'air. Dans le
paragraphe 3.4, l'ecoulement aere uniforme sera examine, tandis que le paragraphe
3.5 concerne l'ecoulement aere non uniforme. Un apen;:u historique sur Jes dudes
d'ecoulement accelere dans Jes coursiers termine le chapitre.
3.2. MECANISME DE L'ENTRAINEMENT D' AIR
Une premiere description du mecanisme d'entrainement d'air fut presentee par

Straub et Anderson [49]. On demontra que !'aeration de l'ecoulement commern;:ait
dans une zone ou la couche limite atteignait la surface de l'eau. En consequence,
!'aeration depend beaucoup de l'intensite de la turbulence pres de la surface libre.
Ainsi que l'ont decrit Straub et Anderson, des« paquets d'eau » jaillissent au-dessus
de la surface de l'ecoulement puis retombent par gravite, de !'air etant ainsi
« insuffle » dans l'ecoulement. L'air est reparti dans celui-ci par turbulence. Deux
zones furent distinguees : I) la zone superieure pres de la surface libre, 2) la zone
inferieure pres du fond du canal.
Volkart [57] considere que !'air est entraine par Jes gouttes d'eau ejectees de
l'ecoulement, qui, en retombant, entrainent des bulles d'air (Fig. 31). Falvey et
Ervine [13] attribuent l'entrainement d'air sur Jes coursiers a la turbulence entretenue
dans l'ecoulement le Jong du coursier.
Le transport de bulles a l'interieur de l'ecoulement necessite une composante de
la vitesse, vers le bas, dans l'ecoulement turbulent, superieure a la vitesse d'ascension
des bulles. Pour des bulles d'air de I a JO mm de diametre, la vitesse d'ascension est
d'environ 0,30 ms- 1• Falvey et Ervine [13] ont montre que des bulles de diametre
superieur a JO mm avaient une forte probabilite d'etre entrainees, mais une faible
probabilite d'etre transportees vers le bas dans l'ecoulement. C'est le contraire pour
72
Sections 3.2 and 3.3 refer to the mechanism of air entrainment and the location
of incipient air entrainment. In 3.4, the uniform aerated flow will be discussed,
whereas section ;LS refers to non-uniform aerated flow. A historical account on the
development of accelerated chute flow concludes the chapter.
3.2. MECHANISM OF AIR ENTRAINMENT
A first description of the air entrainment mechanism was provided by Straub

and Anderson [49]. The aeration of flow was shown to begin in a region where the
boundary layer reached the water surface. Accordingly, the aeration depends
significantly on the turbulence intensity near the free surface. Then, as described by
Straub and Anderson," clumps of water" break through it and fall back by gravity,
whereby the air is " insufflated " into the stream. The air is distributed in the flow
by turbulent transfer. A distinction was made between I) an upper region near the
free surface and 2) the lower region near the channel bottom.
Volkart [57] postulated that air is entrained by droplets of water ejected from the
flow which, when falling back, entrain air bubbles (Fig. 31 ). Falvey and Ervine [ 13]
attributed the air entrainment on chutes to turbulence being continuously generated
within the flow along a chute.
The transport of bubbles into the flow requires a turbulent downward velocity
component larger than the bubble rise velocity. For typical air bubbles of I to I 0 mm
diameter, the rise velocity is about 0.30 ms- 1• Falvey and Ervine [13] were also able
to show that bubbles larger than 10 mm have a high probability of being entrained
but a small probability of being transported down into the flow. The opposite was
found for bubbles smaller than I mm in diameter. Bubbles of diameter between I
Fig. 31
Droplets above the water surface
(Volkart, 1980) [57].
Surface velocity 4.67 ms ',height of water column 40 mm.
Projection de goulles d'eau au-dessus de la sw:face de l'i!coulement
(d'apres Vo/kart, 1980) [57].
Vitesse a la surface 4,67 m/s, hauteur de colonne d'eau 40 mm.
73
!es bull es de diametre inferieur a 1 mm. Ce sont done !es bull es de 1 a I 0 mm que
!'on trouve le plus souvent, !es bulles de petit diametre etant reparties au-dessous de
celles de grand diametre.
3.3. FORMATION DE L'ENTRAINEMENT D' AIR
Le depart de l'ecoulement aere peut facilement etre mis en evidence par I' « eau
blanche », tel que le montre Ia Fig. 5. Pour que l'entrainement d'air se produise,
Gangadharaiah et al. [ 15] donnent deux conditions : 1) I' energie cinetique des
remous en surface doit etre plus grande que l'energie de tension en surface, 2)
l'ecoulement doit avoir atteint un etat entierement turbulent.
Le niveau de turbulence minimal pour que l'entrainement d'air prenne nais-
sance peut s'exprimer sous la forme d'une vitesse de courant V [ms- 1) > 0,24/lf
pour des bulles de 5 mm de diametre, lf etant le coefficient de turbulence [13]. Le
coefficient de turbulence lr est defini comme le rapport de la racine carree
moyenne des variations des vitesses transversales a la vitesse moyenne de
l'ecoulement. La condition 1) pour qu'il y ait entrainement d'air est normalement
satisfaite sur Jes ouvrages reels.
Wood et al. [60) ont etabli la formule suivante pour l'epaisseur o de Ia couche
limite :
0
- = 00212 ( -x.)0.11 (x.)-0.10
- (3.1)
x ' H, k,
cette epaisseur etant independante de la forme de l'evacuateur pour des coursiers a
peu pres rectangulaires. Dans cette formule, H, est la difference de niveau d'eau
entre le niveau du reservoir et le point de depart de l'entrainement d'air, x la distance
par rapport a la crete de l'evacuateur, k, la hauteur de rugosite d'apres Nikuradse
(Fig. 32). L'equation (3.1.) indique que o augmente presque lineairement avec x, et
tres peu avec k,, mais cette acceleration (exprimee par HJ a un effet retardateur.
Le point de depart x = xi survient Iorsque la profondeur de l'ecoulement h est
egale a l'epaisseur de la Couche limite 0 (h = hi) = oi. II en resulte que le profil de
la surface libre doit etre defini par un calcul de ligne de remous et I' equation (3.1)
donnera o; (x;). Une autre methode par approximations successives fut egalement
presentee par Falvey [12] pour la localisation du point de depart de !'aeration.
74
to 10 mm are thus most commonly found, with small diameter bubbles distributed
below large diameter bubbles.
3.3. INCIPIENT AIR ENTRAINMENT
The onset of aerated flow may easily be visualised by the so-called white water
as shown in Fig. 5. Gangadharaiah et al. [15] quoted two conditions for inception :
I) the kinetic energy of the surface eddies must be greater than the surface tension
energy; and 2) the flow must have become fully turbulent.
1]
The minimum turbulence level for incipient air entrainment may be expressed
as streamwise velocity V [ms- > 0.24/ lf for bubbles of 5 mm diameter, with
lf as turbulence number [13]. The turbulence number lf is defined as the ratio of
root-mean square (rms) value of the transverse velocity fluctuations to the mean
flow velocity. Condition l) for incipient air entrainment is normally satisfied
for prototype structures.
Wood et al. [60] found for the thickness o of boundary layer :
Q= 0.0212 (~)0.11 (~)-0.10 (3. l)

x H, k,
independent of the spillway shape for approximately rectangular chutes. Here, H, is
the difference of water surface elevations between inception point and reservoir
level, x the distance from the spillway crest and k, the roughness height according
to Nikuradse (Fig. 32). Eq. (3. l.) shows that o increases almost linearly with x, and
slightly with k,, but that acceleration (expressed by H,) has a retarding effect.
The point of inception x = X; occurs if the flow depth h is equal to the boundary
layer thickness o(h = h;) = o;. As a result, the free surface profile must be computed
by a backwater procedure and Eq. (3. l) must be applied to determine o; (x;). An
alternative trial and error method for locating the point of incipient aeration was also
presented by Falvey [12].
Fig. 32
Occurrence of white water in chute.
Formation d"eau blanche sur un coursier.
75
Exemple 3.1
Soit un evacuateur de forme prismatique de pente 8 = 15° (S = 0,259) et de 0
rugosite equivalente ks = 0,8 mm. Trouver le point de depart de l'entrainement d'air

pour un debit de 25 m 3 s-- 1 par unite de largeur (q = 25 m 2 s- 1).
La charge sur la crete est H = 312 (qlg 112 ) 213 = 6 m suivant la formule de
l'ecoulement critique. A partir d'un calcul de ligne de remous, le profil de la surface
libre peut etre determine comme indique dans le Tableau ci-dessous; x est l'abscisse
longitudinale ayant comme origine la crete de l'evacuateur. Figurent egalement dans
ce Tableau : z, cote du radier du coursier au-dessous de Ia crete de I'evacuateur;
Hs = H-z-h; 8, epaisseur de la couche Ii mite (Fig. 32), donnee par !'equation (3.1 ).
On constate que le point de depart de l'entrainement d'air est situe a xi = 120 m ou
hi = 1,05 m.
x [m] 0 55 110 165 225 280 335
h [m] 4,0 1,30 1,08 0,96 0,86 0,78 0,73
z[m] 0 - 14,2 - 28,5 - 42,7 - 58,2 - 72,5 - 86,7
H,[m] 2,0 18,9 33,4 47,7 63,3 77,8 92,0
8 [m] 0 0,49 0,98 1,48 2,02 2,53 3,03
3.4. ECOULEMENT AERE UNIFORME
Dans un ecoulement aere uniforme un equilibre entre la force de masse et la

force retardatrice s'est etabli. Un tel ecoulement peut se produire asymptotiquement
a l'extremite des longs coursiers prismatiques, de pente de radier constante, mais
physiquement s'assimile mieux a une quantite de regularisation.
Soit C la concentration moyenne d'air dans le temps, en un point, definie
comme le volume d'air par volume de melange. La profondeur equivalente d'eau
pure est :
hw = Joo (1 -
0
C) dy (3.2)
avec y, ordonnee perpendiculaire a !'axe du canal. La vitesse moyenne de l'eau

(indice « w ») est al ors :
(3.3)
ou q est le debit d'eau par unite de largeur. Une profondeur caracteristique du

melange air-eau est la profondeur h 90 ou la concentration moyenne d'air est de 90 %.
La concentration moyenne d'air, prise sur la profondeur, et designee par C, peut etre
definie par !'equation suivante :
(3.4)
76
Example 3.1
Given a wide prismatic spillway of bottom slope 8 = 15° (S 0 = 0.259) and
equivalent roughness k, = 0.8 mm. Find the point of inception for a unit discharge
q = 25m 2 s- 1•
The energy head on the crest is H = 312 (q/g 112) 213 = 6 m according to the
critical flow formula. Based on a backwater computation, the free surface profile
may be given as indicated in the Table below. The longitudinal coordinate x has its
origin at the spillway crest. Also included are the location z of chute bottom below
the crest, the difference head H, = H-z-h and the thickness of boundary layer o
(Fig. 32) according to Eq. (3.1). The point of inception is located at xi = 120 m where
hi = 1.05 m.
x [m] 0 55 I IO 165 225 280 335
h[m] 4.0 1.30 1.08 0.96 0.86 0.78 0.73
z [m] 0 - 14.2 - 28.5 - 42.7 - 58.2 - 72.5 - 86.7
H,[m] 2.0 18.9 33.4 47.7 63.3 77.8 92
5 [m] 0 0.49 0.98 1.48 2.02 2.53 3.03
3.4. UNIFORM AERATED FLOW
In uniform aerated flow an equilibrium between body and retarding forces has
become established. Equilibrium flow may occur asymptotically at the end of long
prismatic chutes of constant bottom slope, but is physically more relevant as a
normalising quantity.
Let C be the local time-averaged air concentration defined as volume of air per
volume mixture. The equivalent clear water depth is then :
hw = f~0
(J - C) dy (3.2)
with y as coordinate perpendicular to the channel axis and with q as water discharge
per unit width. The average water (index " w ") velocity is :
Yw = _9_ (3.3)
hw
A characteristic depth of the air-water mixture is the depth h 90 where the average
air concentration is 90 %. The depth-averaged mean air concentration C may then
be defined as :
(3.4)
77
Rao et al. [35] ont etabli la relation suivante entre le rapport de la densite du
melange a la densite de l'eau, et la concentration moyenne d'air c :
_e_ = I - 1,1 C (3.5)

Pw
Wood (61) a demontre que le profil de la concentration d'air et la repartition des

vitesses dependaient exclusivement de la coordonnee Y90 = y/h 90 et de la pente du
coursier 8. La Fig. 33 indique Jes profils de concentration Ce (Y90 ) pour differentes
valeurs de 8. L'indice « e » correspond a l'ecoulement en equilibre.
Wood [61) a egalement demontre que le coefficient de frottement fe de l'ecoule-
ment en equilibre dependait de la rugosite relative k/h, du rayon hydraulique Rh et
de la concentration moyenne d'air Ce. D'apres Jes resultats obtenus par Straub et
Anderson, le rapport fe/fw depend seulement de C., fw etant le coefficient de
frottement pour un ecoulement d'eau pure. Pour Ce < 20 %, l'effet du melange
air-eau est negligeable. Pour des concentrations moyennes d'air Ce plus grandes, le
rapport f/fw_ diminue lorsque Ce augmente. Wood a etabli une courbe f/fw en
fonction de Ce, pour une rugosite de coursier type (Fig. 34 a). La relation peut etre
exprimee approximativement comme suit :
(3.6)
Cain et Wood (5, 6] ont demontre que !'air entraine n'influern;:ait pas la
repartition des vitesses de l'ecoulement aere uniforme. Avec µ = v/V90, OU V90
q/h 90 , le profil des vitesses suivait la relation [7] :
µ = Y90116
(3.7)
independante de Ce. La profondeur relative d'ecoulement est Y90 = y/h 90 , comme
defini precedemment.
0.2
0.2 0.4 0.6 0.8 1
78
Rao et al. [35] correlated the density ratio of mixture to water flow and the mean
air concentration C as follows :
2-=l-l.lC (3.5)
Pw
Wood [61] was able to show that the air concentration profile and the velocity
distribution depend exclusively on the coordinate Y90 = y/h 90 and on the chute
slope 8. Fig. 33 shows the concentration profiles Ce (Y90 ) for various values 8. Index
" e " accounts for equilibrium flow.
Wood [61] also showed that the friction factor fe of equilibrium flow depends
on the relative roughness k/h, the hydraulic radius Rh and the mean air
concentration Ce. Using Straub and Anderson's data, the ratio f/fw depends only on
Ce. Herein, fw is the friction factor for pure water flow. For Ce < 20 %, the effect
of air-water mixture is negligible. For larger mean air concentrations Ce, the ratio
fe/fw decreases with increasing Ce. Wood presented a curve fe/fw as a function of Ce
for typical chute roughness (Fig. 34 a). It may be approximated as
(3.6)
Cain and Wood [5, 6] demonstrated that entrained air does not affect the
velocity distribution of uniform aerated flow. With µ = V/V90 where V90 = q/h 90 ,
the velocity profile followed the relation [7] :
(3.7)
independent of Ce. The relative flow depth is Y90 = y/h 90 , as previously defined.
Fig. 33
Concentration profiles C, (Y,,,) for uniform aerated flow
and various bottom slopes 8 according to Wood (1985) [61]
based on the data of Straub and Anderson ( 1960) [49].
Profils de concentration d"air C,. (Y",) dans un ecoulement aere uniforme
pour differentes pentes de radier e
(d'apres Wood, 1985 {61} et /es resu/tats obtenus par Straub et Anderson, 1960 [49}).
79
Le profil de concentration d'air d'un ecoulement aere uniforme peut etre
analyse en se referant aux etudes de Wood [59] et de Hager [17]. La Fig. 33 montre
que Ce (Y90) depend seulement de la pente du fond 8, en considerant la rugosite aux
limites du canal d'essai de Straub et d'Anderson. En designant par C 0 Ia concentra-
tion asymptotique au fond avec :
C 0 = 1,25 8 3 [rad], 0° < 8 <40°
C 0 = 0,65 sin 8 [rad], 40° < 8 <80° (3.8)

on peut demontrer que la concentration d'air normalisee c = (Ce - C 0)/(0,90 - C 0)
varie comme suit [17] :
c = 1,05 exp [- 0,54 (y* + 0,30) 2] (3.9)
ou y* = 2 (l - Y90 )/(sin 8) 112 correspond a une position mesuree au-dessus du fond
du canal. Ainsi, toutes Jes courbes de concentration Ce (Y90 ) pour differents 8
peuvent etre transformees en un graphique unique c (y*), comme represente sur la
Fig. 34 b).
La concentration moyenne d'air Ce dans une section transversale peut etre
calculee par integration sur Ia section. On obtient :
Ce = 0,75 (sin 8) 0 •75 (3.10)

indiquant une augmentation presque lineaire de Ce avec la pente 8 du coursier.
En s'appuyant sur un grand nombre d'experiences, Gangadharaiah et al. [15]
ont propose pour la concentration moyenne d'air Ce l'equation:
(3. 11)
Le parametre de coursier Q tient compte de la forme du canal et de la hauteur
de rugosite, et IFw = Vw/(ghw 3 ) 112 est le nombre de Froude de l'ecoulement d'eau
pure. Pour des canaux rectangulaires, on a obtenu le resultat - qui est inexact du
point de vue dimensionnel - n = 1,35 n, n etant le coefficient de rugosite de
Manning.
Avec hw = (l - Ce) h 99 , le rapport Y99 = h 99 /hw peut s'ecrire [44] :
y 99 = I + __l_,3_5_TJ~ (3. 12)
I + 2 (hwlb)
1 I 1
c
0.5 - 0.5
y*
I I
0
0.2 0.4 0.6 0.5 1 1.5
b) 0
80
The concentration profile of uniform aerated flow may further be analyzed
when following Wood [591 and Hager [17]. According to Fig. 33, Ce (Y90 ) depends
only on the bottom slope 0 for the boundary roughness existing in Straub and
Anderson's test channel. If C 0 is the asymptotic bottom concentration with
3
C0 = 1.25 0 [rad], 0° < 0 <40°
C0 = 0.65 sin 0 [rad], 40° < 0 <80° (3.8)

the normalized air concentration c = (Ce C 0)/(0.90 - C 0 ) varies as [17] :
c = 1.05 exp [- 0.54 (y* + 0.30) 2] (3.9)

112
The parameter y* = 2 (I - Y90)/(sin 0) corresponds to a scaled position over
the channel bottom. All concentration curves Ce (Y90 ) for various 0 may thus be
transformed in a single plot c (y*) as shown in Fig. 34 b).
The average cross-sectional air concentration Ce may be computed by

integration over the cross-section and obtains :
Ce = 0.75 (sin 0)075 (3.10)

This indicates an almost linear increase of Ce with the chute slope 0.
Based on a large number of experiments, Gangadharaiah et al. [151 proposed for
the mean air concentration Ce :
(3.11)
The chute parameter Q accounts for the channel shape and roughness height,
and IFw = Vw/(gh/) 112 is the Froude number of pure water flow. For rectangular
channels, the dimensionally incorrect result Q = 1.35 n with n as Manning's
roughness coefficient was found.
With hw = (l - C.) h 99 , the ratio Y 99 = h 99 /hw obtains [441 :
Y99 = I + t. 35 TJ (3.12)
I + 2 (hwlb)
Fig. 34
Uniform aerated flow.
Ecoulement aere uniforme.
a) Effect of mean air concentration on friction a) Effet de la concentration moyenne d'air sur le
factor. coefficient de frottement.
b) Air concentration distribution c (y*). h) Repartition de la concentration d"air c (r*).
81
ou T] = [hw sin 3 8/(g3 n 2)] 114 represente la caracteristique de rugosite dimensionnelle
du coursier. Pour les coursiers larges, la surface moyenne du melange eau-air Y 99
augmente lineairement avec T], c'est-a-dire essentiellement avec la pente e du
coursier, et tres peu avec le coefficient de rugosite. L'effet de la profondeur
d'ecoulement est tres faible. L'equation (3.12) est une base simple pour le calcul de
la profondeur h 99 de l'ecoulement aere uniforme sur le coursier, en considerant la
profondeur h" de l'ecoulement d'eau pure.
Les equations (3.11) et (3.12) peuvent etre transformees en rempla~ant le debit
par la formule de Manning. En se limitant aux coursiers larges (hw/b ~ I), la
concentration moyenne d'air Ce peut s'exprimer comme suit [17]
ce = 1,35 ri (3.13)
La concentration moyenne d'air a l'equilibre Ce depend principalement de la
pente e du coursier, comme l'indiquait aussi !'equation (3.10), et egalement du
coefficient de rugosite de Manning.
Exemple 3.2
Quelles sont Jes caracteristiques de l'ecoulement aere uniforme de l'exemple
precedent 3.1? Avec ks 116 /(ng 112 ) = 8,2 selon Strickler pour un ecoulement rugueux
turbulent, on obtient pour la rugosite de Manning la valeur n = 0,0119. En
outre, avec Q = 1,35 n = 0,016 et le nombre de Froude IFw = q/(ghw 3 ) 112 =
25/(9,81 x 0,72 3 ) 112 = 13,06, I' equation (3.11) donne pour Ce, concentration
moyenne d'air, la valeur de 43 % . L'equation (3.8) donne comme concentration au
fond la valeur approximative de 2,2 %. De plus, avec T] = 0,55 et hw/b ~ I, on obtient
pour la profondeur de l'ecoulement du melange air-eau Y99 = 1,75, et done h 99 =
1,28 m.
3.5. ECOULEMENT AERE, GRADUELLEMENT VARIE
La zone situee entre le point de naissance de l'entrainement d'air et l'ecoule-

ment aere uniforme est graduellement variee et a ete designee sous le nom de zone
d'ecoulement aere, graduellement varie (Fig. 30). Cette zone est caracterisee par une
variation locale de concentration d'air, de vitesse et de profondeur d'ecoulement;
cette zone peut etre assimilee a un ecoulement graduellement varie dans un canal,
excepte que la concentration d'air C et la densite p sont des inconnues supplemen-
taires. En se basant sur la similitude de l'ecoulement graduellement varie et de
l'ecoulement uniforme pour ce qui concerne la pente de frottement [8], on a stipule
le concept de similitude pour les profils de concentration d'air dans Jes ecoulements
aeres [61]. On obtient ainsi, pour le profil de surface, !'equation suivante
dh
sin e- Sr - h sin e (~)
(3.14)
2
dx cos 8 - IF
ou d8/dx correspond a la variation de la pente du fond et IF q/(gh 3) 112 est le
2
nombre de Froude. La pente de frottement Sr = (fa/ 4 Rh) (V /2 g) est calculee a
partir de !'equation (3.6) en tenant compte de la concentration d'air comme indique
precedemment.
82
in which TJ = [hw sin 3 E>/(g 3 n 2)] 114 is the dimensional chute roughness characteristics.
For wide chutes the mean surface of the air-water mixture Y99 increases linearly with
TJ, that is mainly with the chute slope E> and moderately with the roughness
coefficient. The effect of flow depth is very small. Eq. (3.12) is a simple basis for
computing the uniform aerated flow depth h 99 in chute flow when accounting for the
pure water flow depth hw.
Eqs. (3.11) and (3.12) may be developed further when substituting Manning's
formula for discharge. Restricting considerations to wide chutes (hJb < 1), the
mean air concentration Ce may be expressed as [17]
(3.13)
The average equilibrium air concentration Ce depends mainly on the chute slope
E> as was also predicted with Eq. (3 .10) and also on the Manning roughness
coefficient.
Example 3.2
What are the uniform aerated flow characteristics of example 3.1? With
k, 116 /(ng 112 ) = 8.2 according to Strickler for turbulent rough flow, Manning's
roughness value obtains n = 0.0119. Further, with Q = 1.35 n = 0.016 and the
Froude number IFw = q/(ghw 3) 112 = 25/(9.81 x 0.72 3) 112 = 13.06, the average air
concentration is Ce = 43 % according to Eq. (3.11 ). The bottom concentration C 0 is
2.2 %, roughly, from Eq. (3.8). Further, with T] = 0.55 and hw/b < 1 the mixture flow
depth obtains Y 99 = 1.75, thus h 99 = 1.28 m.
3.5. GRADUALLY VARIED, AERATED FLOW
The region between the inception point and the uniform aerated flow is
gradually varied, and is referred to as gradually varied, aerated flow zone (Fig. 30).
This domain is characterized by a local variation of air concentration, velocity and
flow depth. The region is similar to gradually varied channel flow, except that both
air concentration C and density p are additional unknowns. Based on the similarity
of gradually varied and uniform flow as regards the friction slope [8] the concept of
similarity for air concentration profiles in aerated flow was stipulated [61]. The
equation for the surface profile then obtains :
dh
sin E> - Sr - h sin E> (~)
dx cos E> - IF2 (3.14)
where dE>/dx corresponds to the change of bottom slope and IF = q/(gh 3) 112 is the
Froude number. The friction slope Sr = (f/ 4 Rh) (V 2/2 g) is computed from Eq. (3.6)
by accounting for the effect of air concentration as previously discussed.
83
Une methode approchee considere la profondeur hum de l'ecoulement aere
uniforme comme la profondeur d'ecoulement asymptotique, mais utilise !'equation
relative a l'ecoulement non aere graduellement varie. lei, l'indice (( u » correspond
a « uniforme » et « m » au « melange air-eau ». Avec Jes parametres non
dimensionnels :
x = sin e. x/hum•
on obtient [44] :
(3.15)
La Fig. 35 represente la solution Ym(X) pour diverses valeurs de Ye- En partant

du point de naissance de l'entrainement d'air (xi, hJ, et connaissant la profondeur
h de l'ecoulement aere uniforme et la profondeur hem de l'ecoulement critique
a6;e, on peut localiser sur le graphique le point de naissance non dimensionnel
(Xi = sin 8 · x/hum; Y mi = h/hum).
La courbe avec la valeur Ye souhaitee est egalement choisie. On peut alors
proceder, a l'aval le long du coursier, c'est-a-dire en allant de la gauche vers la droite
a)-20 -15 -10 -5

Fig. 35 a)
Free surface profile Ym(X) for gradually varied aerated flow.
Profil de la surface fibre Y"' (X) pour un ecoulement aere graduellement varie.
( ... ) Critical flow. (. .. ) Ecoulement critique.
a) - 20 ,,;;; X ,,;;; 0. a) - 20,,;;; X,,;;; 0.
84
An approximate procedure accounts for the uniform aerated flow depth hum as
the asymptotic flow depth but the equation for non-aerated gradually varied flow.
Index " u " refers to uniform and " m " to the air-water mixture. With the
nondimensional parameters :
x I
= sin e. x/hum,
the result obtains [44) :
dhm sin 0 - Srm
(3.15)
dx 1 - IF,;
Fig. 35 shows the solution Ym (X) for various values of Ye. Starting at the
inception point (xi; h) and with known uniform aerated flow depth hum and aerated
critical flow depth hem, the nondimensional inception point (Xi = sin 0 · x/hum;
y mi = h/hum) may be located on the plot.
Also, the curve with the desired Ye value is selected. Then one may proceed
downstream along the chute that is from left to right in the diagram to obtain the
... . . . 5
x..__ ym
4
4.5 ....__ 4 3.5

Ye 3
\
2
b)-100 -80 -40
Fig. 35 b)
Free surface profile Y m(X) for gradually varied aerated flow.
Profil de la surface fibre Y"' (X) pour un ecoulement aere graduellement varie.
( ... ) Critical flow. (. ..) Ecoulement critique.
b) - 100 ,,;;; x , ; ; - 20. b) - JOO ,,;;; X,,;;; - 20.
85
du graphique, a l'obtention du profil de la surface de transition aux points x choisis,
OU Jes profondeurs du melange sont hm = y mhum·
D'autres etudes relatives a la zone OU se developpe l'entrainement d'air ont ete
presentees par Keller et al. [28], et par Thandaveswara et Rao [51]. Leurs resultats
experimentaux peuvent etre utilises pour l'etalonnage des modeles numeriques.
3.6. APER(:U HISTORIQUE
Parmi Jes premieres etudes sur le phenomene d'entrainement d'air nature!, on

peut citer celles d'Ehrenberger (1926, 1930) [9, 10] et de Ludin (1928) [41]. Ehrenber-
ger put enregistrer la distribution verticale de la teneur en air et etablit que la vitesse
dans une section transversale variait avec (sin a) 0·4, tg a etant la pente du canal. La
Vitesse minimaJe pour qu'un eCOUJement aere existe fut trOUVee egaJe a 4 ms- 1•
Aux Eta ts- u nis, Jes premieres etudes furent effectuees par Lane ( 1936, 1939) [38,
39] et Stevens (1938) [45], bien que Houk (1927) [27] ait fait mention de problemes
d'entrainement d'air au barrage de derivation de Willwood. Une etude systematique
de l'entrainement d'air dans Jes canaux a ecoulement libre fut presentee par Hall
(1943) [22]. II fut etabli que l'entrainement d'air dependait principalement de la
rugosite aux parois, du nombre de Froude local et de la geometrie locale du canal.
Hickox ( 1945) [24] a decrit la naissance de !'aeration sur Jes coursiers a forte pente.
Gumensky (1949) [16] trouva que l'eau aeree avait un effet negligeable sur le
comportement des bassins d'amortissement. D'autres resultats experimentaux furent
obtenus par Viparelli (1951, 1953, 1954, 1958) [52, 53, 54, 56], Halbronn (1951, 1952,
1955) [18, 19, 21] et Halbronn et al. (1953) [20], Straub et Lamb (1953, 1956) [47, 48],
Straub et al. ( 1954) [46], Yevdjevich et Levin ( 1953) [62], et Levin ( 1955) [40]. Cette
serie de donnees se termina en partie par la contribution de Straub et Anderson
( 1960) [49], analysant la distribution de la teneur en air dans des conditions
d'equilibre, c'est-a-dire d'ecoulement aere uniforme.
En dehors des etudes en laboratoire effectuees par Annemiiller (1961) [3],
Rajaratnam (1962) [43], Anderson (1965) [2], Killen et Anderson (1969) [31], et Lai
( 1968) [33], et des observations in situ faites par Okada et al. ( 1959) [42], Harrison et
Owen (1967) [23], et Aki (1969) [!], une seconde periode d'intense experimentation
debuta dans Jes annees 70. Lakshmana Rao et al. (1968, 1970 et 1971) [34, 35, 36]
mirent au point une methode theorique pour determiner Jes distributions de la teneur
en air et de la vitesse dans Jes zones turbulentes libres et Jes zones turbulentes le long
des murs. Le debut d'entrainement d'air et des formules pour le coefficient moyen
d'entrainement furent etudies par Gangadharaiah et al. (1968, 1970) [14, 15]. La zone
ou se developpe l'entrainement d'air fut analysee par Keller et al. (1974) [28]; Keller
et Rastogi (1975, 1977) [29, 30] examinerent Jes problemes concernant Jes evacuateurs
de crue. Parallelement, Thandaveswara et Lakshmana Rao ( 1978) [51] menerent une
etude experimentale sur la zone ou I' entrainement se developpe. Volkart (1980, 1982)
[57, 58] expliqua le mecanisme de l'entrainement d'air et etendit Jes resultats aux
conduites partiellement remplies, tandis que Falvey (1979) [l l] proposa une relation
pour la teneur moyenne en air dans une section transversale d'un ecoulement aere
uniforme. Cain et Wood (1981 a, b) [5, 6] presenterent un nouvel appareillage de
mesure et continuerent Jes experiences au barrage Aviemore. Wood et al. (1983) [60]
etablirent une methode modifiee pour localiser le debut de l'entrainement d'air;
Wood (1983) [59] reexamina l'ecoulement aere uniforme et presenta un rapport sur
86
transition surface profile at selected locations x, where the mixture depths are hm =
Ymhum·
Further contributions to the developing region were offered by Keller et al. [28],
as well as Thandaveswara and Rao [51 ]. Their experimental data may be used for
calibration of numerical models.
3.6. HISTORICAL REVIEW
Early contributions to the natural air entrainment phenomenon were provided

by Ehrenberger (1926, 1930) [9, 10] and Ludin (1928) [41]. Ehrenberger was able to
record the vertical distribution of air concentration and found that the
cross-sectional velocity varied with (sin a)04 . Herein tg a is the slope of the channel.
The minimum velocity for the presence of aerated flow was set at 4 ms- 1•
First studies in the United States were conducted by Lane (1936, 1939) [38, 39]
and Stevens 1938 [45], although Houk (1927) [27] referred to problems of air
entrainment for the Willwood diversion dam. A systematic analysis of air
entrainment in open channels was presented by Hall ( 1943) [22]. It was found that
air entrainment depends mainly on the boundary roughness, the local Froude
number and the local channel geometry. Hickox (1945) [24] described incipient
aeration on steep chutes. Gumensky (1949) [16] found the effect of aerated water
negligible on the performance of stilling basins. Further experimental results were
obtained by Viparelli (1951, 1953, 1954, 1958) [52, 53, 54, 56], Halbronn (1951, 1952,
1955) [18, 19, 21] and Halbronn et al. (1953) [20], Straub and Lamb (1953, 1956) [47,
48], Straub et al. ( 1954) [46], Yevdjevich and Levin (I 953) [62], and Levin ( 1955) [40].
These series of data was partially concluded by the contribution of Straub and
Anderson ( 1960) [49], who analysed the distribution of air concentration in
equilibrium conditions, that is uniform aerated flow.
Apart from lab studies by Annemiiller ( 1961) [3], Rajaratnam ( 1962) [43],
Anderson (1965) [2], Killen and Anderson ( 1969) [31 ], and Lai (1968) [33], and field
observations by Okada et al. (1959) [42], Harrison and Owen (1967) [23] and Aki
(1969) [l], a second intense period of experimentation started in the Seventies.
Lakshmana Rao et al. (1968, 1970 and 1971) [34, 35, 36] developed a theoretical
approach for the concentration and velocity distributions in the free turbulent and
the wall turbulent zones. Incipient air entrainment and expressions for the average
entrainment coefficient were analysed by Gangadharaiah et al. ( 1968, 1970) [14, 15].
The developing region was investigated by Keller et al. (1974) [28] and Keller and
Rastogi (1975, 1977) [29, 30] applied further results on spillways. In parallel,
Thandaveswara and Lakshmana Rao (1978) [51] conducted an experimental study
on the developing region. Vol kart ( 1980, 1982) [57, 58] explained the mechanism of
air entrainment, and extended the results on partially filled pipes, whereas Falvey
( 1979) [ 11] proposed a relation for the average cross-sectional air concentration in
uniform aerated flow. Cain and Wood (1981 a, b) [5, 6] presented a novel
instrumentation and continued the experiments on the Aviemore dam. Wood et al.
(1983) [60] introduced a modified procedure for locating the incipient air
entrainment, and Wood (1983) [59] re-examined the uniform aerated flow and
presented a state-of-the-art report in I 985 [61], as did Rao and Kobus (1975) [37] and
Kobus (1985) [32], as well as Falvey (1980) [12] and Falvey and Ervine (1988) [13].
87
I'etat de la question en 1985 [61], ainsi que le firent egalement Rao et Kobus (1975)
[37], Kobus (1985) [32], Falvey (1980) [12], Falvey et Ervine (1988) [13]. Vers la fin des
annees 80, l'etude, sur modele mathematique, de l'ecoulement aere dans un canal a
surface libre fut entreprise. Des rapports furent etablis par Hopping et Hoopes
( 1988) [26) et par Chanson ( 1989) [7].
b [m] Iargeur du coursier.

c [-] concentration d'air normalisee.
c [-] concentration.
c [-] concentration moyenne dans une section transversale.
Co [-) concentration d'air au fond (en radier).
0 [m] epaisseur de la couche limite.
f [-] coefficient de frottement.
0 [-) pente du radier (du fond).
H [m] charge energetique.
H, [m] difference de niveaux d'eau en un point.
TJ [sm-113] parametre de rugosite du coursier.
k, [m] hauteur de rugosite equivalente.
µ [-) vitesse normalisee.
n [m-113 s+1] coefficient de rugosite de Manning.
n [sm- 113] parametre du coursier.
q [m2 s- 11 debit par unite de largeur.
Rh [m] rayon hydraulique.
p [kgm- 3) densite.
Sr [-) pente de frottement.
so [-] pente du radier (du fond).
11 [-) coefficient de turbulence.
v [ms- 1) vitesse dans une section transversale.
x [-] coordonnee longitudinale normalisee.
x [m] coordonnee longitudinale.
y [m] cote au-dessus du lit.
y [-) profondeur relative de l'ecoulement.
y* [-] parametre de position.
Ye [-) rapport de Ia profondeur de l'ecoulement critique a la profondeur
de l'ecoulement uniforme.
Ym [-] profondeur d'ecoulement normalisee.
z [m] coordonnee verticale.
88
In the late 1980s, the numerical modelling of aerated open channel flow was
initiated. Contributions were made by Hopping and Hoopes (1988) [26] and
Chanson ( 1989) [7].
3.7. NOTATION CHAPTER 3
b [m] chute width.

c [-] normalised air concentration.
c [-] concentration.
t [-] cross-sectional average concentration.
Co [-] bottom air concentration.
0 [m] thickness of boundary layer.
IF [-] Froude number.
f [-] friction factor.
e [-] bottom slope.
H [m] energy head.
h [m] flow depth.
H, [m) difference of local water elevations.
Tl [sm-113] chute roughness parameter.
k, [m] equivalent sandroughness height.
µ [-] normalised velocity.
n [m- 113 s+ 1] Manning roughness coefficient.
Q [sm-113] chute parameter.
q [m 2s- 1] discharge per unit width.
Rh [m] hydraulic radius.
p [kgm- 3] density.
Sr [-] friction slope.
so [-] bottom slope.
ll [-] turbulence number.
v [ms- 1] cross-sectional velocity.
x [-] normalised length coordinate.
x [m] longitudinal coordinate.
y [m] elevation above bed.
y [-] relative flow depth.
y* [-] position parameter.
Ye [-] ratio of critical to uniform flow depths.
Ym [-] normalised flow depth.

z [m] vertical coordinate.
89
Indices
a air.
c ecoulement critique.
e equilibre.
naissance (debut).
m melange.
u ecoulement uniforme.
w eau.
3.8 REFERENCES
1. AKI S. (1969) : " Field measurements of velocity and pressure on spillway
chute ". XI II /AHR Congress Kyoto, 5 ( 1), 167-172.
2. ANDERSON A. G. ( 1965) : " Influence of channel roughness on the aeration of
high-velocity open-channel flow". XI JAHR Congress Leningrad, I (37), 1-13.
3. ANNEMULLER H. (1961) : "Die Berechnung der Abfluptiefen in Schuprinnen ".
Der Bauingenieur 36, 222-226.
4. BAUER W. J. (1954) : "Turbulent boundary layer on steep slopes". Trans. ASCE
119, 1212-1242.
5. CAIN P. and WOOD I. R. (l 981 a) : " Instrumentation for aerated flow on
spillways ". Proc. ASCE, J. Hydraulics Division 107 (l 1), 1407-1424.
6. CAIN P. and WooD I. R. ( 1981 b) : " Measurement of self-aerated flow on a
spillway ". Proc. ASCE, J. Hydraulics Division 107 (11), 1425-1444; 109 (I),
145-146.
7. CHANSON H. (1989) : " Flow downstream of an aerator - aerator spacing ".
J. Hydraulic Research 27 (4), 519-536.
8. CHOW V. T. (1959) : "Open channel hydraulics ".McGraw-Hill Book Company,
New York.
9. EHRENBERGER R. (1926) : "Wasserbewegung in steilen Rinnen mit besonderer
Beriicksichtigung der Selbstbeli.iftung ". Oesterreichischer Ingenieur- und Archi-
tektenverein 78 (15/16), 155-160; 78 (17/18), 175-179.
10. EHRENBERG ER R. ( 1930) : " Eine neue Geschwindigkeitsformel fiir kilnstliche
Gerinne mit starken Neigungen (Schusstennen) und Berechnung der
Selbstbelilftung des Wassers ". Die Wasserwirtschaft 23 (28), 573-575; 23 (29),
595-598.
11. FALVEY H. T (1979) : " Mean air concentration of self-aerated flows ". Proc.
ASCE, J. Hydraulics Division 105 (I), 91-96; 105 (11), 1469; 106 (8), 1407-1408.
12. FALVEY H. T (1980) : " Air-water flows in hydraulic structures ". US Dept.
Interior, Water Resources Technical Publication, Engineering Monograph 41. US
Printing Office, Denver Col.
13. FALVEY H. T. and ERVINE D. A. (1988) : " Aeration in jets and high velocity
flows ". Model - prototype correlation ofhydraulic structures Colorado Springs, ed.
P. H. Burgi, 25-55.
90
Subcripts
a air.
c critical flow.
e equilibrium.
inception.
m mixture.
u uniform flow.
w water.
14. GANGADHARAIAH T., SEETHARAMIAH K. and LAKSHMANA RAo N. S. (1968): .. On

the mean velocity distribution in turbulent flow in a chute ". Conf Hydraulics
and Fluid Mechanics, The Institution of Engineers, Australia, 124-128; 248-249.
15. GANGADHARAIAH T., LAKSHMANA RAo N. s. and SEETHARAMIAH K. (1970) :
" Inception and entrainment in self-aerated flows". Proc. ASCE, J. Hydraulics
Division 96 (7), 1549-1565.
16. GUMENSKY D. B. (1949 : "Air entrained in fast water affects design of training
walls and stilling basins". Civil Engineering 19 (12), 35-37; 93.
17. HAGER W. H. (1991) : "Uniform aerated chute flow". Proc. ASCE, J. Hydraulic
Engineering, 117 (4), 528-533.
18. HALBRONN G. (1951) : " Mesure des concentrations et des vitesses dans un
courant mixte d'air et d'eau ". La Houille Blanche 6 (3), 394-405.
19. HALBRONN G. (1952) : " Etude de la mise en regime des ecoulements sur !es
ouvrages a forte pente ". La Houille Blanche 7 (l), 21-40; 7 (3), 347-371; 7 (5)
702-722.
20. HALBRONN G., DURAND R. and COHEN DE LARA G. (1953) : "Air entrainment
in steeply sloping flumes ". V JAHR Congress Minneapolis, 455-466.
21. HALBRONN G. (1955) : Discussion to " The development of the turbulent
boundary layer on steep slopes". Proc. ASCE, J. Engineering Mechanics Division
81 (606), l-9.
22. HALL L. S. (1943) : "Open channel flow at high velocities". Trans. ASCE 108,
1394-1434.
23. HARRISON A. J.M. and OWEN M. W. (1967)." A new type of structure for flow
measurement in steep streams ". Proc. Institution Civil Engineers (London), 36,
273-296; 37, 811-818.
24. HICKOX G. H. (1945) : "Air entrainment on spillway faces". Civil Engineering,
15 (12), 562-563.
25. HINO M. (1961) : " On the mechanism of self-aerated flow on steep slope
channels". IX JAHR-Congress Dubrovnik, 123-132.
26. HOPPING P. N. and HOOPES, J. A. (1988) : "Development of a numerical model
to predict the behaviour of air-water mixtures in open channels ". Model -
prototype correlation of hydraulic structures Colorado Springs, ed. P. H. Burgi,
419-428.
91
27. HOUK I.E. (1927) : " New design features in Willwood diversion dam ".
Engineering News-Record 99 (17), 660-664.
28. KELLER R. J., LAI K. K. and WOOD I. R. (1974) : .. Developing region in
self-aerated flow ". Proc. ASCE, J. Hydraulics Division l 00 (4), 553-568; 100 (I I),
1729-I730; 101 (2), 319-320; 101 (9), 1278-I280; 102 (4), 529-530.
29. KELLER R. J. and RASTOGI A K. (1975) : " Prediction of flow development on
spillways ". Proc. ASCE, J. Hydraulics Division I 0 I (9), 1171-1184; I 02 (9),
I401-1404; 103 (6), 664.
30. KELLER R. J. and RASTOGI AK. (I977) : " Design chart of predicting critical
point on spillways". Proc. ASCE, J. Hydraulics Division 103 (12), 1417-1429;
104 (12), 1678.
31. KILLEN J.M. and ANDERSON A.G. (1969) : "A study of the air-entrained flow
in open channels ". XIII JAHR-Congress (Kyoto) 2, 339-347.
32. KoBus H. E. (1985) : " An introduction to air-water flows in hydraulics ".
Mittei/ung Nr. 61, Institut fiir Wasserbau, Universitat Stuttgart, Eigenverlag des
Instituts.
33. LAI K. K. (1968) : " A study of air entrainment in steep open channels ".
Conference on Hydraulics and Fluid Mechanics, The Institution of Engineers,
Australia, 41-44; 241-242.
34. LAKSHMANA RAO N. s., SEETHARAMIAH K. and GANGADHARIAH T. (1968) : "A
new electrical probe for measurement of air concentration in air-water flows ".
J. Institution of Engineers (India) 48 (11), 1677-1693.
35. LAKSHMANA RAO N. s., SEETHARAMIAH K. and GANGADHARAIAH T. (1970) :
" Characteristics of self-aerated flow". Proc. ASCE, J. Hydraulics Division 96 (2),
331-355; 96 (12), 2622-2627; 97 (12), 2077.
36. LAKSHMANA RAO N. s. and GANGADHARAIAH T. (1971) : " Self-aerated flow
characteristics in wall region ". Proc. ASCE, J. Hydraulics Division 97 (9),
1285-1303.
37. LAKSHMANA RAo N. S. and KOBUS H. E. (1975): "Characteristics of self-aerated
free-surface flows ". Water and wastewater current research and practise Nr. I 0,
Erich Schmidt-Verlag, Berlin-Bielefeld-Miinchen.
38. LANE E. W. (1936) : " Recent studies on flow conditions in steep chutes ".
Engineering News-Record 116 (I), 5-7.
39. LANE E.W. (1939) : " Entrainment of air in swiftly flowing water ". Civil
Engineering 9 (2), 89-91; 9 (6), 371.
40. LEVIN L. (1955) : " Quelques reflexions sur la mecanique de l'ecoulement des
melanges d'eau et d'air ". La Houi/le Blanche l l (8), 555-558.
41. LUDIN A (1928) Modellversuche fiir die Schussrinne des
Shannon-Kraftwerkes ". Der Bauingenieur 9 (11 ), J 79-184.
42. OKADA A, KUDO M. and FUKUHARA K. (1959) : " Field tests of high velocity
flow and its air context ». VIII JAHR-Congress Montreal, 1D, l-20.
43. RAJARATNAM N. (1962):" Energy loss in a short-smooth-steep chute". J. Central
Board of Irrigation and Power, 19 (9) 725-730.
92
44. SINNIGER R. 0. and HAGER W. H. ( 1989) : " Constructions Hydrauliques -
Ecoulements stationnaires ". Presses Polytechniques Romandes, Lausanne.
45. STEVENS J.C. (1938) : "Hydraulic knowledge". Engineering News-Record, 120,
236-238.
46. STRAUB L. G., KILLEN J. M. and LAMB 0. P. (1954) : " Velocity measurement of
air-water mixtures". Trans. ASCE 119, 207-220.
47. STRAUB L. G. and LAMB 0. P. (1953) : " Experimental studies of the air
entrainment in open channel flow ". V /AHR-Congress Minneapolis, 425-437.
48. STRAUB L. G. and LAMB 0. P. (1956) : " Studies of air entrainment in
open-channel flow ". Trans. ASCE 121, 30-44.
49. STRAUB L. G. and ANDERSON A.G. (1960) : "Experiments on self-aerated flow
in open channels". Trans. ASCE 125, 456-486.
50. Task Committee on Air Entrainment in Open Channels (1961) : " Aerated flows
in open channels ". Proc. ASCE, J. Hydraulics Division 87 (3), 73-86; 87 (6),
221-231.
51. THANDAVESWARA B. S. and LAKSHMANA RAo N. S. (1978) : " Developing zone
characteristics in aerated flows ". Proc. ASCE, J. Hydraulics Division 104 (3),
385-396; 105 (3), 279-281; 105 (11), 1451-1452.
52. VIPARELLI M. (1951) : " Premessa ad una ricerca sulle correnti rapide ".
L'Energia Elettrica 28 (1), 13-22.
53. VIPARELLI M. (1953): "The flow in a flume with 1 : l slope". V /AHR-Congress
Minneapolis, 415-423.
54. VIPARELLI M. (1954) : "Correnti rapide - risultati in canaletta a 450 ". L 'Energia
Elettrica 31(6),393-405.
55. VIPARELLI M. (1957):" Fast water flow in steep channels". VI /AHR Congress
Lisbon, 039, 1-12.
56. VIPARELLI M. (1958) : " Correnti rapide ". L'Energia Elettrica 35 (7), 633-649.
57. VoLKART P. ( 1980) : " The mechanism of air bubble entrainment in self-aerated
flow ". International J. Multiphase Flow 6, 411-423.
58. VoLKART P. U. (1982):" Self-aerated flow in steep, partially filled pipes". Proc.
ASCE, J. Hydraulics Division l 08 (9), l 029-1046.
59. WooD I. R. (1983) : " Uniform region of self-aerated flow ". Proc. ASCE,
J. Hydraulic Engineering I 09 (3), 447-461.
60. WooD I. R., ACKERS, P. and LOVELESS J. (1983) : " General method of critical
point on spillways ". Proc. ASCE, J. Hydraulic Engineering 109 (2), 308-312.
61. WOOD I. R. (1985) : "Air water flows". Keynote Address XX/ /AHR Congress
Melbourne, 18-29.
62. YEVDJEVICH V. and LEVIN L. (1953) : " Entrainment of air in flowing water and
technical problems connected with it". V JAHR Congress Minneapolis, 439-454.
93
4. AERATEURS D'EVACUATEURS DE CRUE
4.1. INTRODUCTION
Tous les evacuateurs de crue presentent des defauts de construction, le plus
souvent de faible importance, tels que irregularites aux joints de beton, ou decro-
chements dus a des tassements locaux (Fig. 36). De telles irregularites peuvent
provoquer une separation de l'ecoulement, associee a une reduction de pression
locale derriere l'irregularite [2, 38]. Les ecoulements sur Jes evacuateurs de crue ont
des vitesses de 20 a 40 ms- 1• Comme la reduction depression peut s'exprimer sous
la forme d'un multiple de la charge cinetique a l'approche, la pression dans la zone
de separation peut tomber au-dessous de la pression locale de vapeur Pv·
Le phenomene de vaporisation resultant d'une baisse de pression - au lieu
d'une augmentation de temperature jusqu'a la temperature d'ebullition - est appele
« cavitation ». Lorsque la pression augmente de nouveau au-dessus de la pression
de vapeur, Jes bulles de vapeur deviennent instables et eclatent, produisant ainsi des
micro-jets de vitesse elevee avec des pointes de pression tres fortes et de tres courte
duree. La cavitation en elle-meme n'est pas souvent dangereuse, a moins qu'elle ne
se produise pres des parois de I'ecoulement. Ces pointes depression, atteignant des
charges de IO 000 m de hauteur d'eau, peuvent alors endommager la surface par Jes
phenomenes de fatigue et d'impact [12]. Le beton ou I'acier qui a subi une cavitation
presente I'aspect d'un materiau soumis a une erosion. Les degiits par cavitation
resultent essentiellement d'un processus de fracturation dans des ecoulements d'eau
claire a partir de vitesses comprises entre 12 et 15 ms- 1, en fonction de la geometrie
de l'ecoulement.
On doit empecher par tous Jes moyens Jes degiits par cavitation, car ceux-ci
deteriorent progressivement l'ouvrage. On peut le faire en eloignant Jes zones de
cavitation des parois, en profilant Jes surfaces pour empecher Jes separations dans
J'ecouJement, OU en utiJisant des materiaux resistant a )a cavitation [39, 40).
Cependant, de tels procedes peuvent ne pas etre applicables aux evacuateurs de crue
ayant de grandes surfaces de contact.
L'adjonction d'air dans l'ecoulement a grande vitesse peut constituer un moyen
efficace et economique d'empecher Jes degiits par cavitation. La presence d'air pres
des parois de l'ecoulement agit comme un coussin permettant de reduire l'effet de
l'eclatement des bulles d'air. Cet effet resulte de la diminution de la vitesse sonique
de l'eau due a la presence d'air. Pour des vitesses d'ecoulement atteignant 35 ms- 1,
une teneur en air entre 6 et 8 % peut empecher l'erosion par cavitation [21].
a)
94
4. SPILLWAY AERATORS
4.1. INTRODUCTION
All spillways contain minor construction defects, such as irregularities at
concrete joints or changes due to local settlements (Fig. 36). Such irregularities may
cause separation of flow, associated with a local pressure reduction behind the
irregularity [2, 38]. For spillway flow, the velocities are typically in the range of 20
ms-I to 40 ms- 1• As the pressure reduction may be expressed as a multiple of the
approaching velocity head, the pressure in the separation zone may fall below the
local vapor pressure Pv·
The phenomenon of vaporization due to a pressure decrease - instead of a

temperature increase until boiling temperature - is known as cavitation. As the
pressure increases again above the vapor pressure, vapor bubbles become unstable
and collapse, thereby generating high speed micro-jets associated with intense
pressure peaks of extremely short duration. Cavitation in itself is often not
dangerous, unless it occurs near the boundaries of flow. Then, pressure peaks up to
heads of 10 000 m water column may damage the surface by fatigue and impact
phenomena (12]. Concrete or steel exposed to cavitation looks as if it had
experienced erosion. Cavitation damage is essentially a fracturing process and may
start in clear water flows at velocities between 12 to 15 ms- I, depending on the flow
geometry.
Compared to cavitation, cavitation damage must be inhibited by all means as

it progressively damages a structure. Cavitation damage may be countered by
guiding all cavitational zones away from the boundaries, streamlining the surfaces
to inhibit separation or using cavitation resistant materials [39, 40]. Yet, these
methods may be unsuitable on spillways with large contact surfaces.
An effective and economic means to prevent cavitation damage is the addition

of air to the high velocity flow. The presence of air near the boundaries of flow acts
as a cushion by which the effect of bubble collapse may be rerluced. This effect is
due to the reduction of the sonic velocity of the water by the presence of air. For
approaching velocities up to 35 ms- 1, an air concentration between 6 and 8 % may
effectively inhibit cavitation erosion [21]. When the natural air entrainment is
Fig. 36
Local low pressure zones due to surface irregularities
and typical surface damage.
Zones de basse pression dues a des irregularites de surface
et degats de surface caracteristiques.
(· • ·) Separating streamline. (· • ·) Ligne de courant de separation.
95
Lorsque l'entrainement d'air nature! ne suffit pas a assurer la concentration
d'air necessaire au fond, on peut installer le long du coursier des dispositifs speciaux
d'entrainement d'air, appeles aerateurs de coursier.
Parmi Jes premiers rapports consacres aux aerateurs de coursier, on peut citer
Jes Russes (Galperin, Semenkov et autres) et Jes Bresiliens (Pinto). Eccher et
Siegenthaler [9] presenterent un compte rendu des rapports publies au cours des
annees 70 sur cette question. Des details sur la cavitation dans Jes ouvrages
hydrauliques furent fournis par Hamilton [14]. Les effets d'echelle furent examines
au Symposium d'Esslingen en 1984 (Marcano et Castillejo [17]; Pan et Shao [20];
Wood [43)). Michel [l 8] ainsi que Laali et Michel [16] ont analyse le mecanisme
d'entrainement d'air. L'attention fut attiree sur la pulsation dans la cavite et ses
caracteristiques physiques. En ce qui concerne l'entrainement d'air, on indiqua le
role important du nombre de Froude et du coefficient de cavitation. Des etudes sur
modeles furent effectuees par de Fazio et Wei [8], Wood [43] et Pinto et al. [22, 23,
24, 25, 26, 27], avec comparaison des resultats avec Jes observations sur les ouvrages
reels. Une attention particuliere fut portee aux dispositifs d'aeration du barrage de
Guri (Marcano et al. [17]; Castillejo et al. [3]; Zagustin et al. [46, 47)), du barrage
d'Alicura (Vernet et al. [35)) et du barrage d'Uribante (Coleman et al. [6]). Aksoy et
al. [I], Prusza et al. [28] et Vischer [37) examinerent les aspects pratiques de projet,
tandis que Wood [44) resuma les connaissances actuelles en matiere d'aerateurs de
coursier.
4.2. AERATION EN RADIER

Les degats par cavitation dans Jes ecoulements a grande vitesse sur Jes coursiers
peuvent etre empeches en ajoutant de l'air a l'ecoulement, si l'entrainement d'air
nature! ne suffit pas. Un dispositif relativement simple, permettant d'augmenter la
teneur en air sur le radier d'un coursier, est l'aerateur de fond, ou simplement
l'aerateur. La Fig. 37 represente l'ecoulement a proximite d'un aerateur, en distin-
guant quatre zones [40) :
l) une zone d'approche, caracterisee generalement par un ecoulement aere de
surface;
2) une zone de transition, ou l'ecoulement est devie par une rampe;
3) une zone comportant une aeration en radier, ou de l'air entre dans la cavite
au moyen de tuyaux ou de rampes sur la paroi;
4) une zone d'ecoulement aere en radier resultant de l'aerateur; le long de cette
zone l'ecoulement est graduellement varie.
Pour empecher Jes degats par cavitation, un aerateur supplementaire est
necessaire lorsque la concentration C 0 de l'air au fond descend au-dessous de 6 a 8 %.
L'equation (3.8) donne C 0 ;;.. 6 % si 8 > 21°, pour un ecoulement uniforme. II en
resulte que Jes coursiers de pente relativement faible 8 < 30° en general, et Jes zones
d'acceleration de l'ecoulement sur les coursiers de pente raide sont particulierement
exposes au risque de degats par cavitation si la vitesse est suffisamment elevee pour
causer de la cavitation.
Le long de la cavite a l'aval de l'aerateur, un melange d'eau et d'air pulverise
a forte turbulence existe. Comme la pression dans la cavite est subatmospherique,
!'air est fourni par le dispositif d'aeration et entraine dans l'ecoulement du fait d'une
difference de pression Lip entre !'atmosphere et la cavite. La teneur en air dans la
96
insufficient for the bottom air concentration required, special a1r entraining
structures called chute aerators must be provided along the chute.
Early contributions to chute aerators were given by the Russians Galperin,

Semenkov and others, and the Brazilians around Pinto. A review of papers in the
Seventies was due to Eccher and Siegenthaler (9). Details on cavitation in hydraulic
structures were supplied by Hamilton [14). Effects of scale were discussed in
Esslingen (1984) (Marcano and Castillejo [17]; Pan and Shao [20]; Wood [43]).
Michel [18] as well as Laali and Michel [16] analysed the mechanism of air
entrainment. Emphasis was placed on cavity pulsation and its physical
characteristics. Regarding the entrainment of air, a dominant role of Froude number
and cavitation number was found. Model studies were provided by de Fazio and
Wei [8], Wood [43] and Pinto et al. [22, 23, 24, 25, 26, 27], by comparing the results
with the prototype. The aeration devices of the Guri dam (Marcano, et al. [ 17];
Castillejo et al. [3]; Zagustin et al. [46, 471), the Alicura dam (Vernet et al. [35]) and
the Uribante dam (Coleman et al. [6]) received particular attention. Aksoy et al. [l],
Prusza et al. [28] and Vischer [37] referred to practical design aspects, whereas
Wood [44] summarised the modern knowledge on chute aerators.
4.2. BOITOM AERATION

Cavitation damage in high velocity chute flow may be prevented when adding
air to the flow provided the natural air entrainment does not suffice. A relatively
simple structure to increase the bottom air concentration is the bottom aerator, or
simply the aerator. As shown in Fig. 37 the flow in the vicinity of an aerator consists
of four zones which include [40] :
l) approaching zone, usually a surface aerated flow;
2) transition zone where the approaching flow is deflected by a ramp;

3) bottom aeration zone where air enters the cavity through air ducts or wall
ramps; and
4) bottom aerated flow due to aerator, along which gradually varied flow
occurs.
To prevent cavitation damage, a further aerator is required when the bottom air
concentration C 0 decreases below 6 to 8 %. From Eq. (3.8), C 0 ;;;., 6 % if 0 > 21° for
uniform flow. Therefore, chutes of relatively small slope 0 < 30° in general, and
accelerating flow portions of steep chutes are particularly in danger of cavitation
damage provided the velocity is sufficiently high to cause cavitation.
Along the cavity downstream from the aerator, a highly turbulent spray-air
mixture exists. As the cavity pressure is subatmospheric, air is supplied by the
aeration system and entrained in the flow due to a pressure difference ~p between
the atmosphere and the cavity. Thus, the cavity air concentration at zone is ®
97
cavite de la zone Q) est done constante pour une valeur ~p au-dessous de la pression
atmospherique (Fig. 37). Dans la zone d'impact du jet (zone de transition entre la
zone Q) et la zone @), la pression au fond augmente nettement jusqu'a une valeur
maximale. Plus a l'aval, la pression en radier tend asymptotiquement vers la
profondeur de l'ecoulement, et la concentration d'air au fond diminue progressive-
ment par suite de depart d'air.
A partir d'etudes sur modeles et sur ouvrages reels en Russie et au Bresil,
Vischer et al. [36) et Volkart [39) ont propose trois types de dispositif d'entrainement
d'air : le deflecteur, le decrochement et la rainure. Le dernier type n'a pas ete retenu,
etant donne que les rainures sont sujettes a submersion, ce qui entraine une rupture
de l'entrainement d'air par suite de penetration d'eau a l'interieur des tuyaux
d'alimentation en air.
Un aerateur doit conduire a un entrainement d'air tres uniforme sur la largeur
du coursier et avec le moins de perturbation possible dans l'ecoulement sur le
coursier. Comme dispositif optimal d'entrainement d'air, on peut recommander une
combinaison du deflecteur et du decrochement. Des dispositifs d'alimentation en
air, tels que ceux representes sur la Fig. 38, ont ete proposes par Volkart et
Rutschmann [40). Ces dispositifs ne provoquent que de petites ondes transversales
et permettent une distribution d'air sur toute la section transversale du coursier.
Rutschmann [27), Rutschmann et Volkart [32] ont propose une methode pour
determiner l'ecoulement d'air dans un conduit de distribution (Fig. 38 b)). L'ecoule-
ment sur le coursier, a l'aval de l'aerateur, reste a peu pres plan.
Le fonctionnement hydraulique d'un aerateur de coursier peut etre decrit
comme suit. Une difference de pression ~p entre !'atmosphere et la cavite entraine
un ecoulement d'air Qa dans le dispositif d'alimentation en air. Le debit d'air depend
principalement de ~p et de la geometrie du dispositif d'alimentation en air. Le debit
Qa peut etre calcule par extension de !'equation de Bernoulli, tenant compte de la
somme des pertes de charge depuis !'atmosphere a travers le conduit d'alimentation
et le dispositif eventuel de distribution jusqu'a la cavite sous le jet.
L'entrainement d'air dans l'ecoulement sur le coursier depend principalement
de la difference de pression ~p dans la cavite, des caracteristiques du jet et de la
geometrie de l'aerateur. Dans ce qui suit, seul le dernier phenomene sera examine,
etant donne que le dispositif d'alimentation en air depend beaucoup des caracteris-
tiques geometriques du canal d'approche [29, 31).
98
constant at a value ~p below atmospheric pressure (Fig. 37). At the jet impact zone
(transition from zone ® to @)) the bottom pressure significantly increases to a
maximum value. Further downstream, the bottom pressure head tends
asymptotically to the flow depth, and the bottom air concentration decreases
gradually due to air detrainment.
Based on both model and prototype analyses in Russia and Brazil, Vischer et
al. [36] and Volkart [39] proposed three types of air entrainment devices, namely the
deflector, the offset and the groove. Since the latter are susceptible to submergence,
that is the breakdown of the air entrainment as water is flowing back into the air
supply conduit, grooves are excluded.
An aerator should entrain the air most uniformly over the width of chute and
with the least disturbance in the chute flow. The combination of deflector and offset
may be recommended as the optimum air entrainment device. Air supply systems,
such as shown in Fig. 38, were proposed by Volkart and Rutschmann [40]. These
involve only a small development of cross-waves and are particularly able to
distribute the air across the entire chute. Rutschmann [27] and Rutschmann and
Volkart [32] presented an approach for the determination of air flow in a distribution
duct as shown in Fig. 38 b). The chute flow downstream from this aerator remains
nearly plane.
The hydraulics of a chute aerator may be described as follows. A driving

pressure ~p between the atmosphere and the cavity provides an air flow Qa across
the air supply system. The air discharge is mainly dependent on ~p and the geometry
of the air supply system. Q. may be computed by an extended Bernoulli equation
which accounts for the sum of head losses from the atmosphere across the duct and
eventual distribution system to the cavity below the jet.
The air entrainment in the chute flow, on the other hand, depends mainly on
cavity pressure difference ~p, the characteristics of the mixture jet and the aerator
geometry. In the following, only the latter phenomenon will be studied, as the air
supply supply system depends strongly on the entrance geometry [29, 31].
Fig. 39 shows schematically the flow at a chute aerator. The chute axis is
Fig. 37
Schematic flow in the vicinity of aerator.
Representation schematique d'un ecoulement a proximite d'un aerateur.
(... ) Bottom pressure head curve. (. ..) Courbe de la hauteur piezometrique en radier.
(For details see text.) (Se reporter au texte pour /es details.)
99
•·>·,:::: . ·•·•••:ta·.·.·.·.Eit·u ,,,,.
.••'•· •••.•:. ··:.·;···:•·····}·'•.+•····1 ;:;:;:;:;•C!:•·•'····:·.········

a)
Fig. 38
Air supply systems for chute flow deflectors.
Dispositifs d'alimentation en air pour des deflecteurs d'ecoulement sur un coursier.
a) Without distribution duct. a) Sans conduit de distribution.
b) With distribution duct. b) Avec conduit de distribution.
La Fig. 39 represente schematiquement l'ecoulement sur un aerateur de cour-

sier. L'axe du coursier est suppose rectiligne dans la section longitudinale et dans la
vue en plan. h et V sont la profondeur et la vitesse moyenne de l'ecoulement
d'approche, 8 la pente du coursier, a !'angle du deflecteur, t, et t, les hauteurs du
deflecteur et du decrochement ( « marche »), Li la longueur du jet depuis l'aerateur
jusqu'au point de pression maximale au fond (Fig. 37).
Les caracteristiques de l'entrainement d'air d'un aerateur de fond peuvent etre
exprimees par le taux d'entrainement d'air P = Q.!Qw ou les indices « a » et « w »
correspondent respectivement a !'air et a l'eau. Si l'on peut negliger !es effets de la
viscosite, de la tension de surface et de la turbulence [I I], p peut s'exprimer comme
suit :
p= f (IF, IE, 8, a, ~· ~) (4.1)
Dans cette equation, IF = V (gh)- 112 est le nombre de Froude a l'approche et

IE = V2 (t.p/p)- 1 le nombre d'Euler de la cavite. Rutschmann et Hager [33, 44] ont
demontre que l'effet de la sous-pression de la cavite pouvait etre represente comme
suit :
-p = [2- arctg (3. I 0- 3 t.IE)1~ (4.2)

Pm ax 1t
ou Pmax = P (t.p = 0) est le coefficient maximal d'entrainement d'air pour une

pression de nappe nulle. t.IE = IE - IEmin est la difference du nombre d'Euler, avec
100
e
a)
Fig. 39
Schematic side view at chute aerator.
Vue laterale schematique d'un aerateur de coursier.
a) Detail. a) Detail.
b) Overall view. b) Vue generale.
assumed straight both in the longitudinal section and in plan view. h and V are the
approaching flow depth and average velocity, 8 the slope of chute, a the deflector
angle, t, and ts the height of deflector and offset (« step ») and Lj the length of jet
from aerator to the point of maximum bottom pressure (Fig. 37).
The air entrainment characteristics of a bottom aerator may be expressed by the

air entrainment rate p = Q/Qw in which indices « a » and « w » refer to air and
water, respectively. If the viscosity, surface tension and turbulence effects can be
overlooked [I I], p may be expressed as :
p= f (IF, IE, e, a, ~· ~) (4.1)
Herein, IF
= V (gh)- 112 is the approaching Froude number and
IE = V (~p/p)- the cavity Euler number. Rutschmann and Hager [33, 44] were able
2 1
to show that the effect of cavity subpressure is :
- p = [2- arctg (3 . I 0- 3 ~IE)]~ (4.2)

~max 1t
Pmax = P(~p= 0) is the maximum air entrainment coefficient for zero

nappe pressure, further ~IE = IE - IEmin is a difference Euler number with
101
IEmin = IE (~ = 0). La derniere expression correspond a la pression relative de la
cavite pour laquelle l'entrainement d'air disparait. Emin peut etre ajuste par un
ensemble limite de donnees :
1
10- 3 IE min= - -(tg a) 1•15 exp [l,15 (tg 8)2] + _!_T 2 (4.3)
2,3 3 s
ou T, = t/h. Pour une geometrie donnee d'aerateur, on peut calculer IEmin et I'effet
de la sous-pression lip de la cavite depend exclusivement de ~max-
Suivant que le coursier a l'amont de l'aerateur est court (comme c'est le cas dans
Jes laboratoires d'hydraulique) ou long (comme sur Jes ouvrages reels), on peut faire
une distinction entre l'ecoulement d'approche en cours d'evolution et l'ecoulement
d'approche etabli. Dans le dernier cas, l'ecoulement graduellement varie (aere)
domine et on peut appliquer ce qui suit.
A partir des etudes faites par Pinto et al. [22, 23], Minor [ 19], Pinto [27], et des
resultats experimentaux obtenus par Rutschmann [30, 31], le coefficient maximal
d'entrainement d'air ~max peut etre relie a la longueur relative du jet \ = L/h sous
la forme suivante [33, 34] :
~max = 0,030 (Ai - 5) pour IF > 6 (4.4)
Ce resultat est interessant dans la mesure OU le taux maximal d'entrainement
d'air depend lineairement de la longueur de la surface de separation entre l'eau et
!'air. L'entrainement d'air depuis un aerateur est done un phenomene avec domi-
nance de cisaillement, ainsi que cela peut etre demontre par d'autres analyses
(Fig. 40).
Dans le cas d'un ecoulement d'approche etabli, la longueur relative du
jet A,i = L/h pour un aerateur en forme de rampe peut etre calculee comme suit [33] :
:~s": [1 (1 + 2~~;;2 8)
112
A,i = + ] (4.5)
et pour un aerateur en forme de decrochement (a = 0)
"-i = w(~)112 (4.6)

cos 8
Dans ces equations, T = (t, + tJ/h est egal a la hauteur relative de l'aerateur et
a est l'angle de tir effectif du jet (Fig. 39 b). On a :
a
(4.7)
a
102
IEmin = IE (~ = 0) as the relative cavity pressure at which the air entrainment
vanishes. Emin was adjusted by a limited set of data as :
1
10- 3 IE min = - -(tg a) 115 exp [l.15 (tg 8)2) + ! T s2 (4.3)
2.3 3
where Ts= t/h. For any given aerator geometry, IEmin may be computed, and the
effect of cavity subpressure ~p depends exclusively on the maximum air entrainment
rate ~max-
Depending on whether the aerator approaching chute is short (as typically in
hydraulic labs) or long (as in prototypes), a distinction between so-called developing
and developed approaching flow may be made. In the latter case, gradually varied
(aerated) flow prevails and the following may be applied.
Based on considerations by Pinto et al. [22, 23), Minor [19), Pinto [27) and
experimental data by Rutschmann [30, 31), the maximum air entrainment coefficient
~max may be related to the relative length of jet, Ai = L/h as [33, 34] :
~max = 0.030 (Ai - 5) for IF> 6 (4.4)

This result is interesting in so far as the maximum air entrainment rate depends
on the length of interface between water and air. Aerator air entrainment therefore
is a shear-dominated phenomenon, as may also be corroborated by other evidence
(Fig. 40).
For developed approaching flow, the relative length of jet Ai = L/h for ramp
aerators may be computed as [33) :
A_ <i·IF
J - cos e
2
[i + (i + 2T·cos8)
(aIF)2
112
] (4.5)
and for offset aerators (a = 0)

11
Ai= IF ( -2T-) 2
(4.6)
cos e
Herein, T = (t, + tJ/h is equal to the relative height of aerator and a is the
effective take-off angle (Fig. 39 b).
a (4.7)
a
Fig. 40
Bottom aerator with shear-dominated air-water interface.
Aerateur de fond avec inteiface de cisaillement air-eau.
103
ou Tr= t/h est la hauteur de la rampe et Tgh est la fonction tangente hyperbolique.
Le calcul du taux d'entrainement d'air devient ainsi direct si Jes conditions suivantes
sont satisfaites : T < I, tga < 0, 15 et 6 < IF < 12. Le coursier ne devra pas
egalement avoir une pente transversale importante, son axe sera rectiligne, et un
ecoulement a peu pres bi-dimensionnel devra s'etre etabli. A ce jour, aucune etude
experimentale systematique n'a ete effectuee pour verifier le present modele
mathematique. Rutschmann et Hager [34] ont pu verifier la presente methode a partir
de resultats obtenus sur des ouvrages reels.
Une autre methode d'etude, faisant intervenir le niveau de turbulence sur la
rampe, !'angle de sortie du jet et Jes pertes a travers Jes orifices d'aeration, a ete
proposee par Glazov [13].
Exemple 4.1.
Soit un coursier de pente de radier 8 = 28°, avec un ecoulement de profondeur
h = 1,2 m et de vitesse V = 27 ms- 1 a l'approche. Que! est le taux maximal
d'entrainement d'air Pmax pour un aerateur de pente a= 4°, de hauteur tr= 0,2 m et
de hauteur de decrochement ts= 0,2 m?
Le nombre de Froude de l'ecoulement d'eau pure a l'approche est
IF= 27/(9,81 x 1,2) 112 = 7,9 et l'aerateur peut etre caracterise geometriquement
par a = 0,070, T = 0,333 et Tr = 0, 167. On a done a = 0,983 x 0,070 = 0,069
d'apres !'equation (4.7). De plus, \ = 13,3 d'apres !'equation (4.5) et, done,
Pmax = 0,25. On a done pour l'entrainement maximal d'air par unite de Jargeur
qa = 0,25 x 27 x 1,2 = 8,1 m3 s- 1•
Exemple 4.2.
Que! est l'effet d'une difference depression ~p/(pg) = 0,3 m de hauteur d'eau
sur l'entrainement d'air?
Le nombre d'Euler est egal a IE = (V2 /g)/(~p/pg) = 247,7. En outre, avec
a = 4°, 8 = 28° et T, = 0,167, on obtient : IEmin = 65,2 (eq. 4.3). II en resulte :
~IE = 247,7 - 65,2 = 182,5 et PIPmax = 0,449. Le coefficient d'entrainement d'air
est done (eq. 4.2) : p = 0,449 x 0,25 = 0,112.
D'autres recommandations pour le projet d'aerateur ont ete donnees par
Rutschmann et Hager [34]. La pente de l'aerateur devrait etre comprise entre 5° et
10°. De plus, on a demontre que l'efficacite optimale de l'aerateur etait atteinte pour
1 1
W~max a/a
0.5 0.5
10- 3·.11E Trf a

0
0.5 1 1.5 b)o 0.5 1 1.5
104
where T, = t,/h is the height of ramp and « Tgh » the hyberbolic tangent function.
Computation of the air entrainment rate thus becomes straightforward if the
following conditions are satisfied : T < 1, tga < 0.15 and 6 < IF < 12. Also, the
chute should have no significant transverse slope, its axis is straight and nearly
two-dimensional flow should have established. Actually, no systematic experimental
studies were conducted to verify the present computational model. Rutschmann and
Hager [34] were able to verify the present approach by prototype data.
An alternative design procedure involving the turbulence level at the ramp, the
exit angle of jet and the losses through the air vent was proposed by Glazov [13].
Example 4.1.
Given a chute with approach flow depth h = 1.2 m, velocity V = 27 ms - 1 and
bottom slope 8 = 28°. What is the maximum air entrainment rate Pmax for an aerator
of slope a= 4°, height t,= 0.2 m and height of step t, = 0.2 m?
The approach Froude number is IF = 27 /(9.81 x l.2) 112 = 7.9 and the aerator
geometry may be described by a = 0.070, T = 0.333 and T, = 0.167. Thus
a = 0.983 x 0.070 = 0.069 according to Eq. (4.7). Further, \ = 13.3 according to
Eq. (4.5) and thus Pmax = 0.25. The maximum air entrainment per unit width is thus
qa = 0.25 X 27 x 1.2 = 8.1 m 3 s- 1•
Example 4.2.
What is the effect of a pressure difference ~p/(pg) = 0.3 m water column on the
air entrainment?
The Euler number obtains : IE = (V2 /g)/(~p/pg) = 247.7. Further, with
a = 4°, 8 = 28° and T, = 0.167, one obtains : IEmin = 65.2 from Eq. (4.3). As a
result, ~IE = 247.7 - 65.2 = 182.5 and PIPmax = 0.449. The air entrainment
coefficient thus is from Eq. (4.2) : p = 0.449 x 0.25 = 0.112.
Additional guidelines for aerator design were given by Rutschmann and
Hager [34]. Regarding the slope of aerator, values between 5° and I 0° should be
adopted. It was further shown that the optimum aerator efficiency occurs for
Fig. 41
Graphical Display.
Representation graphique.
a) Of Eq. (4.2). a) De /'equation (4.2).
b) Of Eq. (4.7). b) De /'equation (4. 7).
105
~opt = (1/3) ~max· Cette relation permet d'estimer la surface de la section des conduits
d'alimentation en air. Rutschmann et Hager [34] ont egalement etabli une relation
entre la hauteur optimale de l'aerateur et le nombre de Froude a l'entree, et Jes pentes
du coursier et de l'aerateur. Avec ces indications et une verification de projets
d'aerateurs existants [45], le projet d'aerateur standard ne doit pas poser de proble-
mes.
4.3. ECHAPPEMENT DE L' AIR
A partir de la methode de Wood [43], Chanson [5] a admis qu'au-dela du point

d'impact du jet, c'est-a-dire dans la zone d'ecoulement aere graduellement varie
(Fig. 37), le taux d'entrainement d'air etait constant et egal au taux a l'equilibre Ce.
En outre, le taux d'echappement d'air a ete suppose proportionnel a la concentration
locale d'air C = C (x). La variation locale du debit d'air qa est done :
dq" = (Cc - C) u, cos e (4.8)

dx
ou u, est la vitesse d'ascension des bulles.
L'equation de continuite (I - C) qa C qw conduit a
dC = U,·COS e (Cc - C) (I - C)2 (4.9)

dx qw
On peut ecrire la condition aux Ii mites comme suit : c (0) = C; + c*. Pour une
concentration d'air nulle a l'entree C; = 0, on a C (0) = C*, tandis que pour un
ecoulement aere uniforme a l'entree la condition aux limites s'ecrit : c (0) =
Ce + C*. Considerant la condition aux limites C (0) = C*, avec :
c* = -~- ( 4.10)
1+ ~
!'equation (4.2) conduit avec Jes parametres :
J_ J- Ce
I -- C (4.11)
c= - '
1 _1--Ce
I - C*
a la relation suivante :
x* = - [I - c+ Inc] (4.12)
ou X (u, cos8/qw) x. On a trouve comme vitesses d'ascension des bulles d'air :

u, = 0,16 ms- 1 sur des modeles reduits et 0,40 ms- 1 sur des ouvrages reels. La
106
Popt = (1/3) Pmax· This relation provides an estimate for the area of air supply
conduits. Rutschmann and Hager [34) also related the optimum height of aerator to
the approach Froude number, and the slopes of chute and aerator. With these
indications, and a check on existing aerator designs [45], the design of the standard
aerator becomes straightforward.
4.3. AIR DETRAINMENT
Based on Wood's [43) approach, Chanson [5] assumed that beyond the jet
impact point, that is in the gradually varied aerated flow zone (Fig. 37), the rate of
air entrainment is constant and equal to the equilibrium rate Ce. Further, the rate of
air escape was assumed proportional to the local air concentration C = C (x).
Therefore, the local change of air flow qa is with u, as the bubble rise velocity
dq.
_a= - -C)
(Ce - Ur COS
e
H (4.8)
dx
Using the continuity equation (I - C) qa C qw yields
dC = U,·COS e (Ce -C) (I - C)2 (4.9)

dx qw
The boundary condition may be expressed as c (0) Ci + c*. For
approaching zero air concentration Ci = 0 this simplifies to C (0) = C*, whereas
for uniform aerated approaching flow the boundary condition reads C (0) =
cc + C*. Accounting for the boundary condition c (0) = c* with :
(4.10)
according to Eq. (4.2) yields with the parameters
I_ I - Ce
1-C (4.11)
the two parameter set :
x* = - [I - c + Inc] (4.12)
where X = (u, cos8/qw) x. The bubble rise velocity was found as u, = 0.16 ms- 1
for scale models and 0.40 ms- 1 for prototypes. The difference in u, may be explained
107
difference de u, peut s'expliquer par la difference de diametre db des bulles d'air [7].
En general, on peut admettre que la vitesse d'ascension des bulles d'air est constante.
A I'aval du point d'impact du jet, Jes bulles d'air soot redistribuees et tendent
asymptotiquement vers la concentration d'air a l'equilibre, telle qu'elle est definie
par )'equation (3.9). Une fois u, defini, le profil de la concentration d'air C (x) dans
la zone de transition peut etre calcule suivant )'equation (4.12). De plus, en suppo-
sant un profil de teneur en air C (y) defini suivant )'equation (3.9), une teneur
moyenne minimale cmin = 30 % est necessaire pour assurer une teneur en air de 5 %
a 10 % pres du fond. II en resulte que Jes aerateurs de fond soot absolument
necessaires pour Jes coursiers de pente 8 < 20°, comme deja vu au chapitre 3, a
condition que l'ecoulement soit presque uniformement aere. Pour Jes evacuateurs OU
8 est superieur a 30°, avec un ecoulement aere presque uniforme, aucun aerateur ne
sera necessaire puisque l'entrainement d'air nature! est suffisant pour le proteger le
fond contre l'erosion par cavitation.
Exemple 4.3
Considerons I'exemple 4.2. Quel est le profil longitudinal de la concentration
moyenne d'air C (x) si l'ecoulement a l'amont de l'aerateur ne contient pas d'air et
si la vitesse d'ascension des bulles d'air est u, = 0,40 ms- 1 ?
Avec ~ = 0,112, la concentration d'air initiale correspondante est C* = 0,101
d'apres )'equation (4.10) si l'ecoulement d'approche n'est pas aere. L'equation (3.11)
donne pour la concentration d'air a l'equilibre Ce = 0,259. Le Tableau ci-dessous
indique !'evolution de la concentration moyenne d'air C dans la section transversale.
On constate que la concentration d'air Ce est atteinte asymptotiquement au bout de
I 00 m environ.
x [m] 0 10 20 40 60 100
x* 0 0,34 0,68 1,36 2,04 3,40
c 1 0,385 0,236 0,105 0,050 0,012
c 0,101 0,205 0,227 0,245 0,252 0,257
4.4. CONCEPTION DES AERATEURS

La conception et !'emplacement des aerateurs de coursier dependent beaucoup
des caracteristiques particulieres du site : topographie, geologie et methode de
construction. Selon Rutschmann [30, 31 ], on doit tenir compte des points suivants :
- Le profil en long du coursier sera constitue de parties de pente constante,
s'adaptant a la topographie. Les points de changement de pente peuvent etre
consideres comme des emplacements possibles pour Jes aerateurs de fond. II en
resulte que Jes ouvrages de transition peuvent etre evites. On peut aussi placer Jes
aerateurs aux courbes dans le plan vertical, le centre de courbure etant eloigne de
l'eau.
108
by the various sizes of bubble diameter d 0 [7]. Normally, the bubble rise velocity is
assumed constant.
Downstream from the jet impact point, the air bubbles are redistributed and
tend asymptotically to the equilibrium air concentration as given in Eq. (3.9). The
transition concentration profile C (x) may be computed with Eq. (4.12) once u, is
determined. Further, upon assuming an air concentration profile C (y) according to
Eq. (3.9), a minimum average concentration Cmin = 30 % is needed to secure a near
bottom air concentration of 5 % to I 0 %. This is to say that bottom aerators are
definitely needed for chutes with slopes 8 < 20° as already discussed in chapter 3
provided the flow is almost uniformly aerated. For spillways where 8 > 30° with
almost uniform aerated flow, no aerator will be required since the natural air
entrainment is sufficient to protect the bottom from cavitational erosion.
Example 4.3
Consider example 4.2. What is the longitudinal profile of mean air
concentration C (x) if no air is contained upstream from the aerator, and the air
bubble rise velocity is u, = 0.40 ms - i?
With r) = 0.112, the corresponding initial air concentration is C* = 0.1 OI
according to Eq. (4.10) if the approaching flow is not aerated. From Eq. (3.11) the
equilibrium air concentration obtains Ce = 0.259. The Table below gives the
development of cross-sectional mean air concentration C. The asymptotic air
concentration Ce is reached within 100 m, roughly.
x [m] 0 JO 20 40 60 JOO
x* 0 0.34 0.68 1.36 2.04 3.40
c I 0.385 0.236 0.105 0.050 0.012
c 0.101 0.205 0.227 0.245 0.252 0.257
4.4. DESIGN OF AERATORS

The design and location of chute aerators greatly depends on site-specific
characteristics, such as topography, geology and the method of construction.
According to Rutschmann [30, 31] the following items need consideration :
- The longitudinal chute profile should be composed of constant sloping
portions in agreement with the topography. Changes of chute slope may then be
considered as possible locations for bottom aerators. As a result, transitional
structures can be avoided. Aerators may also be located at curves in the vertical
section with the center of curvature away from the water.
109
,. 152
(mJ
Fig. 42
Longitudinal chute profile with adaptation
to topography (Rutschmann, 1988) [30, 31].
Pro/ii en long d'un coursier s'adaptant a la topographie
(d'apres Rutschmann, 1988) [30, 31]
- A partir d'une courbe de remous de l'ecoulement d'eau pure le long de tout

le coursier, le coefficient de cavitation :
(J = p/(pg) - pJ(pg) (4.13)

yz/(2 g)
peut etre calcule. lei, p/(pg) est la charge piezometrique sur le radier du coursier, et
p/(pg) ~ - 10 m la pression de vapeur. La pression sur le radier se compose d'une
partie statique et d'une partie dynamique resultant de la courbure des lignes de
courant. On peut ainsi tracer une courbe cr (x) du coefficient de cavitation et le
comparer au coefficient de debut de cavitation cri. Ce dernier depend de la qualite
du beton et de la teneur en air de l'eau. Pour cr < cri, la cavitation se produit et une
aeration au fond est necessaire pour empecher Jes degats par cavitation.
- Selon Falvey [10], on peut adopter cri = 0,20 pour Jes ouvrages en beton
classique, ce qui correspond a une vitesse Iimite de 30 a 40 ms- 1 environ. II ya lieu
de noter que le debit maximal n'est pas necessairement la condition determinante.
Une application de cette methode a ete presentee par Volkart et Rutschmann [41].
Exemple 4.4
La Fig. 42 donne un exemple OU )'emplacement des aerateurs a ete determine
uniquement par rapport a la topographie du coursier. Dans cet exemple, un plus
petit nombre d'aerateurs aurait suffi pour la protection du coursier, mais cela aurait
110
26.3°
~
oy
-+----------------~~----~- ---
Fig. 43
Aerator for concave bottom geometry.
Aerateur pour une forme concave de radier.
- Based on a backwater curve for pure water flow along the entire chute, the
cavitation number :
(J = p/(pg) - pj(pg) (4.13)

y2/(2 g)
may be computed. Herein, p/(pg) is the pressure head at the chute bottom, and
p/(pg) """ - 10 m the vapor pressure head. The bottom pressure is composed of the
static and dynamic portions due to streamline curvature. One may thus plot a
cavitation parameter curve cr (x) and compare it with the incipient cavitation number
cr;. The latter depends on the quality of concrete and air content in the water. For
a < cr;, cavitation occurs and bottom aeration is needed to prevent cavitation
damage.
- According to Falvey [IO] one may set cri = 0.20 for usual concrete structures,
corresponding to a limit velocity of roughly 30 to 40 ms- 1• Note that the maximum
discharge is not necessarily the determining condition. An application of this
method was presented by Volkart and Rutschmann [41].
Example 4.4
Fig. 42 shows an example in which the location of aerators was chosen only
relative to the topography of chute. In this example, less aerators would have been
sufficient for chute protection, but this would have required relatively costly
111
necessite des ouyrages de transition relativement couteux. L'aerateur@a une forme
tres particuliere, etant situe dans une auge, et une submersion pouvant se produire
pour un deflecteur trop petit. Dans ce cas, ii y aurait penetration d'eau dans le
dispositif d'alimentation en air et coupure du circuit d'air. D'autre part, le detlecteur
ne doit pas avoir une pente trop raide, ni etre trop long, car ii pourrait caviter a son
extremite aval; l'ecoulement pourrait egalement etre rejete trop loin du coursier et
des murs de grande hauteur seraient necessaires. La Fig. 43 donne une solution pour
l'exemple considere. Elle consiste en un deflecteur courbe ressemblant a une cuillere
de dissipation.
L'aeration en radier est tres sensible a un exces d'aeration conduisant a des murs
lateraux du coursier de hauteur excessive et a une augmentation de la vitesse due a
la diminution du frottement sur Jes parois (Fig. 34 a). Le bassin de dissipation
d'energie peut se comporter differemment de ce qui a ete prevu, par suite d'une
insuffisance d'eau a l'aval, ou meme de la chasse du ressaut. En conclusion, on peut
indiquer que le projet d'aerateurs doit comprendre des etudes sur modele a une
echelle appropriee.
Ce qui precede concerne exclusivement Jes canaux rectilignes dont la pente
transversale du radier est nulle. Pour Jes canaux presentant une courbe, le projet des
aerateurs est beaucoup plus complique et des effets nuisibles causes par des ondes
transversales peuvent se produire. Actuellement, la plupart de tels organes soot
etudies sur modeles reduits. Une importante contribution a ete presentee par Wei et
de Fazio [42] qui ont resolu Jes equations regissant l'ecoulement potentiel par la
methode des elements finis. Les projets futurs pourraient s'appuyer sur cette
approche lorsque Jes conditions aux limites soot correctement prises en compte. Le
modele reduit restera le moyen le plus commode pour etudier Jes ecoulements sur
Jes coursiers, a condition que l'echelle du modele soit bien choisie pour permettre
une transposition acceptable des resultats a l'ouvrage reel.
a [-] angle du detlecteur.

a [-] angle de tir effectif.
~ [-] taux d'entrainement d'air.
c [-] concentration moyenne d'air.
c* [-] concentration d'air pour un ecoulement aere uniforme.
ci [-] concentration moyenne initiale d'air.
Co [-] concentration d'air au fond (en radier).
c [-] concentration relative d'air.
db [m] diametre d'une bulle d'air.
~IE [-] difference du nombre d'Euler.
~p [Nm- 2] difference de pression.
IE [-] nombre d'Euler.
e [-] pente du coursier.
112
trans1t10n structures. Aerator @ has a quite specific shape as it is located in a
trough and submergence might occur for a too small deflector. Then, water would
run back into the air supply system and cut the air stream. On the other hand. the
deflector must not be too steep and long as it could cavitate at its downstream
extremity. Also, the flow could be deflected too far from the chute and impose high
chute walls. Fig. 43 shows a solution as proposed for the example considered. It
consists of a curved deflector and resembles a flip bucket.
Bottom aeration is quite susceptive to overaeration of which effects are

excessive height of chute side walls and particularly increase of velocity due to
decreasing wall friction (Fig. 34 a). The stilling basin could then behave differently
as designed due to insufficient tailwater or even sweepout of the jump. As a general
conclusion, the design of aerators should thus include model studies at an
appropriate scale.
The previous discussion was limited to straight channels with zero transverse
bottom slope. For curved channels, the aerator design is much more involved and
adverse effects from cross-waves may occur. Actually, most of such elements are
analysed by scale models. An outstanding contribution was forwarded by Wei and
de Fazio [42], who solved the governing equations for potential flow numerically by
the finite element method. Future designs could be based on this approach when
correctly incorporating the boundary conditions. The scale model will remain the
most convenient means to design chute flows provided its scale permits a reasonable
up-scaling to the prototype.
4.5. NOTATION CHAPTER 4
a [-J deflector angle.

i'i [-] effective take-off angle.
~ [-] air entrainment rate.
c [-] average air concentration.
c* [-] air concentration for uniform aerated flow.
ci [-] initial average air concentration.
Co [-] bottom air concentration.
l: [-] relative air concentration.
db [m] bubble diameter.
liIE [-] difference Euler number.
,lip [Nm- 2] difference pressure.
IE [-] Euler number.
IF [-] Froude number.
e [-] chute slope.
l 13
Li [m] longueur du jet.
\ [-] longueur relative du jet.
p/(pg) [m] charge piezometrique.
Q [m3s-1] debit.
q [m2s-1] debit unitaire.
p [kgm- 3] densite.
(J [-] coefficient de cavitation.
(Ji [-] coefficient de debut de cavitation.
T [-] hauteur relative de l'aerateur.
T, [-] hauteur relative du deflecteur.
T, [-] hauteur relative du decrochement.
t, [m] hauteur du deflecteur.
ts [m] hauteur du decrochement.
u, [ms- 1] vitesse d'ascension des bulles d'air.
v [ms- 1] vitesse.
x [-] coordonnee normalisee.
x* [-] coordonnee longitudinale normalisee.
Indices
a air.
e equilibre.
max maximum.
min minimum.
v vapeur.
w eau.
4.6. REFERENCES
I. AKSOY S. and ETHEMBABAOGLU S. (1979) : "Cavitation damage at the discharge

channels of Keban dam ". XIII /COLD-Congress (New Delhi) Q. 50, R. 21,
369-379.
2. BALL J. W. (1976) : "Cavitation from surface irregularities in high velocity".
J. Hydraulics Division ASCE, I 02 (9) 1283-1297.
3. CASTILLEJO N. and MARCANO A. (1988) : " Aeration at Guri final stage
spillways ". Model-Prototype Correlation of Hydraulic Structures Colorado
Springs (Col.) 102-109.
4. CHANSON H. (1989a) : " Study of air entrainment and aeration devices ".
5. CHANSON H. ( 1989b) : " Flow downstream of an aerator - aerator spacing ".
114
h [m] flow depth.
Lj [m] jet length.
A,j [-] normalised jet length.
p/(pg) [m] pressure head.
Q [m3s .. 1] discharge.
q [m2s-1] unit discharge.
p [kgm- 3] density.
(J [-] cavitation number.
CT; [-] incipient cavitation number.
T [-] relative aerator height.
T, [-] relative deflector height.
T, [-] relative offset height.
t, [m] deflector height.
ts [m] offset height.
u, [ms- 1] bubble rise velocity.
v [ms- 1] velocity.
x [-] normalised coordinate.
x* [-] normalised length coordinate.
Subscripts
a air.
e equilibrium.
max maximum.
mm mm1mum.
v vapor.
w water.
6. COLEMAN H. W., SIMPSON A. R. and DE GARCIA L. M. (1983) : " Aeration for

cavitation protection of Uribante spillway". ASCE sponsored Conf. on Frontiers
in Hydraulic Engineering (Cambridge, Mass.) 438-443.
7. COMO LET R. (1979) : " Gas bubble motion in a liquid medium ". (in French). La
Houille Blanche 34 ( 1), 31-42.
8. DE FAZIO F. G. and WEI C. Y. (1983):" Design of aeration devices on hydraulic
structures ". ASCE sponsored Conf. on Frontiers in Hydraulic Engineering
(Cambridge, Mass.) 426-43 l.
9. EccHER L. and SIEGENTHALER A. (1982) : " Spillway aeration of the San Roque
project ". Water Power and Dam Construction 33 (9), 37-40.
10. FALVEY H. T. (1983):" Preventing cavitation in tunnel spillways". Water Power
& Dam Construction 34 (8), 13-15.
11. FALVEY H. T. and ERVINE D.A. (1988) : " Aeration in jets and high velocity
flows ''. Model-Prototype Correlation of Hydraulic Structures, Colorado Springs,
ed. P. H. Burgi, 25-55.
115
12. FALVEY H. T (1990) : "Cavitation in chutes and spillways". Water Resources
Technical Publication. US Printing Office, Engineering Monograph 42, US Dept.
of Interior, Bureau of Reclamation, Denver, Col.
13. GLAZOV A. I. ( 1985) : " Calculation of the air-capturing ability of a flow behind
an aerator ledge ". Hydrotechnical Construction 18 (11), 554-558.
14. HAMILTON W. S. (1983) : " Preventing cavitation damage to hydraulic
structures". Water Power 35 (11), 40-43; 35 (12), 48-53; 36 (1), 42-45.
15. KOSCHITZKY H.P. and KOBUS H. (1988) : " Hydraulics and design of spillway
aerators for cavitation prevention in high speed flows ". International
Symposium on Hydraulics for High Dams Beijing, 724-733.
16. LAALI A. R. and MICHEL J.M. (1984) : "Air entrainment in ventilated cavities :
case of a fully developed half-cavity. ". J. Fluids Engineering 106 (9), 327-335.
17. MARCANO A. and CASTILLEJO N. (1984) : " Model-prototype comparison of
aeration devices of Guri dam spillway ". Symp. Scale Effects in Modelling
Hydraulic Structures (Esslingen) 4.6, 1-5.
18. MICHEL J.M. (1984) : "Some features of water flows with ventilated cavities".
J. Fluids Engineering I 06 (9), 319-326.
19. MINOR H.-E. (1984) Erfahrungen mit Schussrinnenbeliiftung ".
Wasserwirtschaft 77 (6), 292-295.
20. PANS. and SHAO Y. (1984) : "Scale effects in modelling air demand by a ramp
slot". Symp. on Scale Effects in Modelling Hydraulic Structures (Esslingen), 4.7,
1-5.
21. PETERKA A. J. ( 1955) : " The effect of entrained air on cavitation pitting ".
V. !AHR Congress Minnesota, 507-518.
22. PINTO N. L. de S. and NEIDERT S. H. (1982):" Model - prototype conformity in
aerated spillway flow ". Hydraulic Modelling of Civil Engineering Structures
Coventry (England), 273-284.
23. PINTO N. L. de S., NEIDERT S. H. and OTA J. J. (1982): "Aeration at high velocity
flows ". Water Power and Dam Construction 33 (2) : 34-44.
24. PINTO N. L. de S. and NEIDERT S. H. (1983) : " Evaluation entrained air flow
through aerators ". Water Power and Dam Construction 34 (8), 40-42.
25. PINTO N. L. de S. (1984) : " Model evaluation of aerators in shooting flow".
Symp. on Scale Effects in Modelling Hydraulic Structures 4.2, 1-6.
26. PINTO N. L. de S. (1986):" Basic hydraulics of shooting flows over aerators".
ASCE Specialty Conf. Advancements in Aerodynamics, Fluid Mechanics and
Hydraulics I 007-1017.
27. PINTO N. L. de S. (1989) : " Designing aerators for high velocity flow". Water
Power & Dam Construction 40 (7), 44-48.
28. PRUSZA Z. V., MANTELLINI T. P. and SEMENKOV v. ( 1983) : " Remedial measures
against spillway cavitation ". XX !AHR Congress (Moscow) 3, 468-476.
29. RUTSCHMANN P. (1987) Die transversale Druckverteilung unter
Sprungstrahlen von Schussrinnenbeliiftern - Folgen filr den Lufteintrag ".
Wasserwirtschaft 77 (5), 219-223.
30. RUTSCHMANN P. (1988a) : "Calculation and optimum shape of spillway chute
116
aerators ". Model - Prototype Correlation of Hydraulic Structures Colorado
Springs (Col.), 118-127.
31. RuTSCHMANN P. ( l 988b) : " Beliiftung von Schu~rinnen - Ausfilhrung und
Bemessung ". Internat. Symposium Erosion, Abrasion und Kavitation.
Mitteilung No. 99, 149-169, Versuchsanstalt filr Wasserbau, Hydrologie und
Glaziologie, ETH-Zurich, ed. D. Vischer, Zurich.
32. RUTSCHMANN P. and VOLKART P. (1988) : " Spillway chute aeration ". Water
Power and Dam Construction 40 ( 1), l 0-15.
33. RUTSCHMANN P. and HAGER w. H. (1990a) : " Air entrainment by spillway
aerators". J. Hydraulic Engineering ASCE, 116 (6), 765-782; 117 (4), 545.
34. RUTSCHMANN P. and HAGER w. H. (1990b) : " Design and performance of
spillway chute aerators ". Water Power and Dam Construction, 42 (1), 36-42.
35. VERNET G. F., ANGELACCIO C. M. and CHIVIDINI M. F. (1988)
" Model-prototype comparison in the Alicura chute spillway air system "
Model-Prototype Correlation of Hydraulic Structures Colorado Springs (Col.)
110-117.
36. V1sCHER D., VoLKART P. and SJEGENTHALER A. (1982) : " Hydraulic modelling
of air slots in open chute spillways ". Hydraulic Modelling of Civil Engineering
Structures Coventry (England), 239-252.
37. VISCHER D. ( 1985) : " Die Schussrinnenbeliiftung als Ma~nahme gegen
Kavitationserosion ". Oesterr. Wasserwirtschaft 37 (3/4), 117-123.
38. VISCHER D. (1987) : " Kavitation an Schussrinnen ". Wasserwirtschaft 77 (6),
288-291.
39. VOLKART P. (1983) : " Air slots for flow aeration ". Mitteilung No. 66,
Versuchsanstalt filr Wasserbau, Hydrologie und Glaziologie, Swiss Federal
Institute of Technology, Zurich, ed. D. Vischer, Zurich.
40. VoLKART P. and RuTSCHMANN P. (1984):" Air entrainment devices". Mitteilung
No. 72, Versuchsanstalt filr Wasserbau, Hydrologie und Glaziologie, Swiss
Federal Institute of Technology, Zurich, ed. D. Vischer, Zurich.
41. VOLKART P. and RUTSCHMANN P. (1986) : " Aerators on spillway chutes :
fundamentals and applications ". Advancements in Aerodynamics, Fluid
Mechanics and Hydraulics Minneapolis (Minnesota), 15 pages.
42. WEI C. Y. and DE FAZIO F. G. ( 1982) : " Simulation of free jet trajectories for the
design of aeration devices on hydraulic structures ". Finite Elements in Water
Resources Proc. IV Int. Conference, Hannover. Springer, Wien-Berlin.
43. WooD I. R. (1984) : " Air entrainment in high speed flows ". Symp. on Scale
Effects on Modelling Hydraulic Structures (Esslingen), 4.1, 1-7.
44. WOOD I. R. (1988) : " Aerators - the interaction of nappe and duct air
entrainment ". Int. Symp. on Hydraulics for High Dams (Beijing) 611-618.
45. WooD I. R. (1991) : "Air entrainment in free-surface flows". JAHR Hydraulic
Structures and Design Manual, 4. A. A. Balkema : Rotterdam/Brookfield.
46. ZAGUSTIN K. and CASTILLEJO N. (1983) : " Model-prototype correlation for flow
aeration in the Guri dam spillway ". XX JAHR-Congress (Moscow) 3, 443-452.
47. ZAGUSTIN K., MANTELLINI T. and CASTILLEJO N. (1982) : "Some experience on
the relationship between a model and prototype for flow aeration in spillways".
Hydraulic Modelling of Civil Engineering Structures Coventry (England), 285-295.
117
Imprimerie de Montligeon
61400 La Chapelle Montligeon
Depot legal : janvier 1992
No 15833
ISSN 0534-8293
Couverture : Olivier Magna
Copyright © ICOLD - CIGB
Archives informatisées en ligne Computerized Archives on line
The General Secretary / Le Secrétaire Général :

André Bergeret - 2004
_______________________________________________________________
International Commission on Large Dams
Commission Internationale des Grands Barrages
151 Bd Haussmann -PARIS –75008
http://www. icold-cigb.net ; http://www. icold-cigb.org

Bul81 - Spillway Shockwaves and Air Entrainement

Transféré par

Informations du documentcliquez pour développer les informations du document

Informations du documentcliquez pour développer les informations du document

Droits d'auteur :

Formats disponibles

Bul81 - Spillway Shockwaves and Air Entrainement

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Bul81 - Spillway Shockwaves and Air Entrainement

Transféré par

Droits d'auteur :

Formats disponibles

This Bulletin has been prepared by Subcommittee No.

Ce Bulletin a ete prepare par le Sous-Comite n° 4

Commission Internationale des Grands Barrages - 151, bd Haussmann, 75008 Paris

Les informations, analyses et conclusions auxquelles cet ouvrage renvoie

Les informations, analyses et conclusions contenues dans cet ouvrage

Malgré tout le soin apporté à la rédaction de cet ouvrage, compte tenu de

Nous déclinons expressément toute responsabilité quant à l'interprétation et

En poursuivant la lecture de cet ouvrage, vous acceptez de façon expresse

The information, analyses and conclusions in this document have no legal

We decline all responsibility whatsoever for how the information herein is

Do not read on unless you accept this disclaimer without reservation.

2. ONDES DE CHOC 2. SHOCK WA YES IN CHUTES

3. ENTRAINEMENT D'AIR 3. NATURAL AIR ENTRAINMENT

4. AERATEURS D'EVACUATEURS 4. SPILLWAY AERATORS

2. ONDES DE CHOC SUR LES COURSIERS .............................................................. 24

3. ENTRAINEMENT D'AIR NATUREL ........................................................................ 70

4. AERATEURS D'EVACUATEURS DE CRUE ............................................................ 94

2. SHOCK WA YES IN CHUTES ...................................................................................... 25

3. NATURAL AIR ENTRAINMENT................................................................................ 71

Fig. I. - Conception des coursiers d'evacuateurs de crue.

Fig. I. - Concepts of spillway chutes.

Fig. 19. - Differences in flow surfaces. a) Plane bottom. b) Curved bottom.

Fig. 28. - Bend flow in trapezoidal channel.

Les evacuateurs de crue et les vidanges de fond representent un pourcentage

Un barrage comprend deux elements importants : l'evacuateur de crue et le

The principles for spillway design, including types, structures, equipment,

Although there exist a number of tunnel spillways such as those discharging

(*) Numbers between brackets refer to references of each chapter.

Chapter 2 focuses on shock waves, which are standing free surface

Enfin, le chapitre 4 est consacre aux aerateurs d'evacuateurs de crue. De tels

a) Small discharges. a) Debit faible.

View at Grand Coulee Dam spillway (Gumensky, 1949) [3].

For a rectangular channel of width b and cross-sectional area A bh, where

2.2. DEVIATION BRUSQUE D'UN MUR

Soit un ecoulement torrentiel parallele dans un large canal rectangulaire dont

The next structure to be considered is the channel expansion, although it

2.2. ABRUPT WALL DEFLECTION

et Hager [9) a montre que :

and Hager [9] showed that :

~ = e + _l_ IF1 1 (2. l 0)

de largeur) q = 15 m 2 s- 1, le coefficient de rugosite n = 0,011 sm- 113 • Quelle est la

I) Uniform flow depth according to Manning-Strickler equation

Y = 1.22 according to Eq. (2.9), that is

2.3. REDUCTION DES ONDES STATIONNAIRES

2.3. REDUCTION OF STANDING WAVES

Precedemment, la deviation du mur 8 avait une valeur finie. Lorsque 8 - 0,

Si, de plus, la charge de l'ecoulement H = h + V2/(2 g) est consideree com me

est obtenue pour la fonction 8 (IF). lei, 8 1 = 8 (IF 1) en [rad] correspond a la

8 - 81 = 2 (IF2 - 1,8) (2.14)

Further, if the energy head H = h + V2I (2 g) is assumed constant relative to the

e- 01 = 2 (IF2 - 1.8) (2.14)

2.5. RETRECISSEMENT D'UN CANAL

2.5.1. Types de retrecissement

·in which IF 1 = V/(gh 1) 112 •

2.5. CHANNEL CONTRACTION

2.5.1. Types of contraction

the bottom slope is increased;

2.5.2. Retrecissement en forme d'entonnoir