ECOULEMENTS À SURFACE
LIBRE
OBJECTIVES
Savoir calculer les propriétés géométriques et hydrauliques d'un
écoulement à surface libre.
Définir les critères d'optimalité de la section et d'érosion pour le
dimensionnement des conduites et des canaux.
Être capable de classifier un écoulement selon son régime et la
variation des propriétés géométriques et hydrauliques.
Savoir calculer la courbe de remous pour un écoulement variant
graduellement.
Connaître les propriétés du ressaut hydraulique et les principes de
construction d'un bassin d'amortissement.
INTRODUCTION
Dans un écoulement à surface libre, le liquide s'écoule au contact de
l'atmosphère. Dans la nature, l'exemple typique d'un écoulement à surface
libre est celui des rivières. En milieu urbain, l'écoulement dans une conduite
d'égout pluvial et sanitaire doit normalement être à surface libre.
Une différence importante entre un écoulement en charge et un écoulement
à surface libre provient de la géométrie de l'écoulement.
Un écoulement en charge peut être considéré comme un état limite
d'écoulement à surface libre.
Les applications du présent chapitre sont nombreuses. Pour
évaluer une conduite d'égout existante, on doit vérifier si
l'écoulement se fait encore à surface. Pour une conduite
projetée, on doit trouver la dimension optimale qui permet la
circulation de l'eau à surface libre tout en gardant le projet
économiquement rentable. Dans le cas des canaux et des
rivières, on s’intéresse aussi bien à connaître la vitesse
d'écoulement pour protéger les berges contre l'érosion et
contre le transport des sédiments.
CLASSIFICATION DES ÉCOULEMENTS
Les écoulements à surface libre peuvent être permanents ou non permanents.
Dans un écoulement permanent, toutes les propriétés de l'écoulement sont
constantes dans le temps. En fait, un tel écoulement n'existe pas en pratique.
Cependant, l'hypothèse de la permanence de l'écoulement est souvent retenue.
Quand le débit d'écoulement varie significativement dans le temps, l'écoulement
est considéré non permanent.
L’écoulement à surface libre peut être uniforme, variant graduellement ou variant
brusquement.
Dans un écoulement uniforme, la profondeur et la vitesse
d'écoulement restent constantes sur une longue distance le
long de l'écoulement. Il faut alors que toutes les propriétés
géométriques et hydrauliques pertinentes restent constantes le
long de l'écoulement. Il s'agit, là encore, d'un écoulement idéal
qui, en toute rigueur, n'existe pas en pratique. L'hypothèse de
l'écoulement uniforme est cependant souvent retenue soit par
approximation soit dans un contexte de design.
Dans un écoulement variant graduellement, la profondeur
et la vitesse d'écoulement varient graduellement le long de
l'écoulement. Le cas typique est l'élévation du niveau d'eau
en amont d'une retenue. Il faut connaître l'effet d'un niveau
d'exploitation donné de la retenue sur l'inondation des terres
en amont.
Dans l'écoulement variant brusquement, la vitesse et le
niveau varient d'une manière sensible et brusque sur une
courte distance. L'exemple typique est celui du ressaut
hydraulique qui se produit au passage d'un écoulement
rapide a un écoulement lent dans le canal de fuite d'une
centrale hydroélectrique.
Finalement, un écoulement uniforme ou variant
graduellement peut encore être caractérisé selon son
régime : il peut être fluvial, critique ou torrentiel.
CARACTÉRISTIQUES GÉOMÉTRIQUES
DES ÉCOULEMENTS À SURFACE LIBRE
SECTION MOUILLÉE D'UN CANAL
La section d'un canal est la section perpendiculaire à son axe. Si le canal est de
forme cylindrique ou prismatique, comme tous les cas considérés ici, l'axe du
canal est parallèle à la génératrice du cylindre ou aux arêtes du prisme.
On appelle section mouillée A d'un canal la portion de la section du canal limitée
par les parois du canal et la surface libre .
PÉRIMÈTRE MOUILLÉ D'UN CANAL
On appelle périmètre mouillé P d'un canal, le périmètre de la section mouillée qui
inclut les parois solides mais ne comprend pas la surface libre
RAYON HYDRAULIQUE
On appelle rayon hydraulique RH le division de l'aire de la section mouillée A et du
périmètre mouillé P :
on peut interpréter le rayon hydraulique comme un indice de la performance de
la section d'écoulement.
PROFONDEUR OU TIRANT D'EAU
On appelle profondeur d'eau la hauteur d'eau au-dessus du point le plus bas de la
section perpendiculaire à l'écoulement. On la dénote souvent par h ou y
LARGEUR AU PLAN D'EAU, LARGEUR AU RADIER
La largeur de la section mouillée à la surface libre est appelée largeur au plan
d'eau. On la dénote généralement par B
La largeur de la section mouillée au niveau du fond du canal est appelée largeur
au radier, représentée par b
PENTE D'UN CANAL
On définit la pente d'un canal (S0) comme étant la
tangente de l'angle d'inclinaison du lit par rapport à
l'horizontale. Comme cet angle est généralement très
petit, la pente est, en pratique, numériquement égale
à l'angle donné en radians.
UN CANAL PRISMATIQUE
Un canal prismatique est un canal dont la pente et la
géométrie de la section restent constantes dans la
direction longitudinale du canal. Lorsqu'une de ces
conditions n'est pas respectée, le canal est dit non
prismatique.
ÉQUATION D' ÉNERGIE POUR LES
ÉCOULEMENTS À SURFACE LIBRE
FORMULATION GÉNÉRALE
l'équation de conservation d'énergie peut s'écrire pour les écoulements à
surface libre avec les notations habituelles :
α est un coefficient de correction qui
tient compte de la non uniformité de la
vitesse dans une section donnée
(comme généralement 1,04< α < 1,14,
on suppose souvent α = 1)
ÉCOULEMENT UNIFORME
Un écoulement est dit uniforme si les conditions suivantes sont satisfaites :
a) La section d'écoulement est prismatique et l'écoulement est permanent .
b) La profondeur de l'écoulement et la vitesse sont constantes le long du canal.
Dans un écoulement à surface libre en régime uniforme, la perte de charge par
unité de longueur ou pente hydraulique de l'écoulement Sf est égale à la pente
géométrique S0. Comme J est nécessairement positif (énergie absorbée par les
frottements), la relation Sf = S0 montre que le régime uniforme ne peut se produire
qu'en canal descendant.
En pratique, que la conduite soit enfouie ou en surface, sa pente sera a priori
égale à celle du terrain naturel, qui est une donnée du problème. Après avoir
choisi a priori le matériau de la conduite, On doit calculer le diamètre
économique pour véhiculer le débit de projet grâce à cette pente. Selon que
cette pente soit très abrupte ou très faible, pour produire des vitesses
respectivement très fortes ou très faibles d'un écoulement uniforme à surface libre,
des rajustements dans le design de la conduite seront nécessaires.
Il existe plusieurs formules pour relier la perte de charge unitaire aux propriétés de
l'écoulement en régime uniforme. La formule la plus utilisée est celle de Manning :
Le coefficient de Manning « n » ne dépend que de la nature des surfaces des parois
DIVERS TYPES DE PROBLÈMES POSÉS EN PRATIQUE
Calcul de la profondeur normale:
Lorsque l'écoulement est uniforme, la profondeur d'écoulement est appelée
profondeur normale.
Normalement, si on ne perturbe pas l'écoulement par l'introduction d'une singularité
quelconque, l'écoulement sera uniforme et la profondeur sera la profondeur
normale. L'équation de Manning peut être utilisée directement pour calculer la
profondeur normale. Cette profondeur sera comparée à la profondeur de
conception pour qualifier l'adéquation de l'ouvrage. La profondeur des écoulements
uniformes peut être calculée à l'aide de la figure :
Application 1
Un canal de section rectangulaire a une largeur B = 6 m et une pente de S0 = 0,0001. Son
coefficient de Manning est de n = 0,015.
De l'eau s'écoule dans ce canal à raison de Q = 6,0 mètres cubes par seconde.
Il faut trouver la profondeur normale y d'écoulement dans ce canal.
Calcul du débit et de la vitesse d’écoulement
Dans ce type de problèmes, on désire connaître la capacité hydraulique d'une
conduite ou d'un canal dans des conditions préétablies. Le débit calculé est comparé
au débit de conception pour vérifier l'adéquation de l'ouvrage face au refoulement
et aux besoins en eau. Le calcul de la vitesse est nécessaire pour vérifier si cette
dernière est assez élevée pour éviter la sédimentation de matières en suspension et
n'est pas trop forte pour causer l'érosion. Ces calculs sont directs et ne nécessitent ni
graphes ni itérations.
Application 2
Un canal ayant une section en forme de trapèze a une largeur au fond b = 5,0m et une pente
longitudinale S0 = 0,2%. La pente des berges est Z = 1,0. Le coefficient de Manning est n = 0,02.
Il faut calculer le débit Q d'écoulement si la profondeur d'eau ne doit pas dépasser y= 1,0m. Il
faut aussi déterminer la vitesse correspondante.
Application 3
Considérons une rivière à section composée, comme montrée sur la figure. La pente
longitudinale moyenne est S0 = 1,0 · 10-4 . Le coefficient de Manning est n = 0,02.
calculer le débit d'écoulement.
Calcul des dimensions d'un canal
Il s'agit d'un problème de conception proprement dit. En effet, le plus souvent le débit
est fixé par les besoins et la disponibilité de la ressource. La pente est imposée car elle
égale a priori la pente du terrain de pose. La rugosité de la conduite est dictée par le
choix du matériau relié au type de projet. En utilisant toutes ces données, On doit
trouver les dimensions économiques de l'ouvrage qui respectent certaines contraintes
de fonctionnement.
Application 4
Un canal à section rectangulaire doit transporter un débit Q = 0,40 m3/s sans dépasser une
profondeur y = 0,60m. Le canal, dont le coefficient de Manning est n = 0,016, est installé selon
une pente S0 = 0,05%. Il faut déterminer la largeur minimale B que l'on doit donner à ce canal.
Calcul du coefficient de Manning
Plusieurs paramètres ont été introduits pour caractériser le frottement de l'eau contre les parois et par
conséquent les pertes de charge.
Pour les écoulements à surface libre, l'équation de Manning qui permet de calculer la vitesse ou le
débit fait intervenir le coefficient de Manning n pour caractériser la rugosité des parois. Tous ces
coefficients ne sont pas directement accessibles à la mesure, même pour des conduites neuves.
Des expériences en laboratoire ou des mesures in situ permettent de calculer la valeur de ces
paramètres en mesurant simultanément toutes les propriétés hydrauliques qui apparaissent dans les
équations impliquées. Ceci est aussi le cas pour une conduite existante dont le vieillissement modifie
la rugosité. On sait que généralement la rugosité d'une conduite augmente avec l'âge mais il est
impossible de connaître, autrement que par la mesure, le coefficient de Manning actuel.
Application 5
Une conduite d'égout de longueur L = 100,0 m et de diamètre D =0,610 m a été installée selon
une pente S0 = 0,075%. Cette conduite véhicule un débit Q = 0,20 m3/s avec une profondeur
y = 0,61m, sans aucune mise en charge.
Trente années après son installation, la conduite n'est plus capable de véhiculer le débit initial de
conception sans provoquer une mise en charge H = 0,50 m à son entrée.
Il faut calculer les coefficients de Manning initial et final si l'on suppose que l'écoulement à la
sortie se fait toujours à surface libre et que les pertes de charge singulières sont négligeables.
Calcul de la pente ou de la perte de charge
Ce calcul consiste à déterminer la pente nécessaire pour assurer le fonctionnement
dans certaines conditions données d'écoulement à surface libre. Quand l'écoulement
est en charge, ce calcul consiste à déterminer la ligne piézométrique. la pente
s'exprime :
Pour une conduite pleine, on peut écrire une formule pour le calcul de la pente Sf de
la ligne d'énergie :
Application 6
Un entrepreneur promet de corriger le problème de la conduite considérée dans l'application 5,
en la réhabilitant par une méthode de tubage. Cette méthode consiste à introduire, à partir d'un
regard, une nouvelle conduite en PVC à l'intérieur de la conduite existante.
Malgré un diamètre plus petit (D = 0,59m), l'entrepreneur prétend qu'avec un coefficient de
Manning plus faible (n = 0,009 pour le PVC), la mise en charge sera sensiblement diminuée. A-t-il
raison?
CAS PARTICULIER DE LA FORME CIRCULAIRE
La section de forme circulaire a retenu plus particulièrement l'attention car elle est
pratiquement adoptée pour toutes les conduites de drainage fabriquées aujourd'hui.
L'usage des formules du tableau précèdent pour calculer les caractéristiques
géométriques d'une conduite de section circulaire pour différents degrés de remplissage
n'est cependant pas aisé. C'est pour cette raison que le tableau et la figure ci dessous
seront utilisés pour résoudre tous les types de problèmes des conduites partiellement
pleines. Ce tableau et cette figure expriment la relation qui existe entre la fraction de
remplissage de la conduite y/D et les autres
propriétés géométriques (A/AP et RH/ RHP) et hydrauliques V/VP et Q/ Qp)·
Toutes les propriétés sont normalisées par rapport aux caractéristiques de la conduite
pleine qui sont faciles à calculer. En effet on peut écrire :
Application 7
Une conduite d'égout de longueur L = 100m et de diamètre D = 0,59m a été installée selon une
pente S0 = 0,075%. Cette conduite, dont le coefficient de Manning est n = 0,009, doit drainer un
débit Q = 0,20m3s. Il faut calculer la profondeur d'eau et la vitesse d'écoulement .
Règles pratiques de dimensionnement des canaux
Pour choisir les dimensions et la pente d'un canal ou d'une galerie en
fonction de conditions d'écoulement spécifiées, on admet que
l'écoulement est uniforme.
Pour effectuer les calculs, on dispose de deux équations : l'équation de
calcul du débit Q = V.A, l'équation de Manning.
Avant d'accepter le choix final des dimensions d'un canal, il faut vérifier
si la vitesse normale d'écoulement n'est ni trop petite ni trop grande.
Les vitesses d'écoulement recommandées dépendent de la nature des
parois et du fond, ainsi que de la charge solide éventuelle que l'eau
doit transporter sans déposition. Le tableau suivant fournit l'ordre de
grandeur de la vitesse admissible en fonction de la nature des berges
et des caractéristiques de l'effluent.
Vitesses maximales acceptables pour canaux rectilignes, de pente modérée et stabilisée
Forme optimale de la section d'un canal
Afin de réduire les coûts d'excavation et de matériaux de construction des
canaux, on cherche à trouver la meilleure forme de la section pour véhiculer
un débit donné. L'équation de Manning montre que pour la même section
d'écoulement A, le débit d'écoulement est d'autant plus grand que le
périmètre mouillé P est faible. Pour les canaux ouverts, on peut montrer de ce
point de vue que la section la plus efficace est celle de forme semi-circulaire.
En prenant la forme semi-circulaire comme référence, les autres formes
usuelles se situent dans l'ordre suivant :
efficacité
section parabolique avec y/B = 0,5 0,986
section trapézoïdale avec z = 0,58; b =1,15y 0,967
section triangulaire avec z = 1 0,923
section rectangulaire avec b = 2y 0,923
NOTION D'ÉNERGIE SPÉCIFIQUE
Définition
L'énergie spécifique E est définie par :
Comme A est fonction de la profondeur, on peut dire que E = E(y) pour un
débit fixe donné.
Divers régimes d'écoulement
D'après la figure, on voit que pour une
énergie spécifique donnée E, il existe
deux profondeurs possibles, y1 et y2,
pour transporter un même débit Q.
L'écoulement se fait avec l'une ou
l'autre des profondeurs, y1 ou y2, selon
la rugosité et la pente du canal.
On voit aussi que pour assurer un débit
donné Q, l'écoulement requiert une
énergie minimale Emin à laquelle
correspond une profondeur critique
yc et une vitesse critique Vc.
L'énergie spécifique E est minimale lorsque
En remplaçant l'énergie par son expression on obtient :
Comme par définition dA/dy = B, on obtient la relation suivante :
Le nombre de Froude « Fr » est un nombre sans dimensions défini par :
Comme dE/ dy = 1 - Fr , trois cas peuvent se présenter :
l'écoulement est fluvial quand Fr < 1; dE/dy > 0 ;
ce cas se trouve dans la branche où V < Vc, y > yc ;
l'écoulement est critique quand Fr = 1, V = Vc, y = yc ;
ce cas se trouve au point le plus bas de la courbe où la pente est
nulle ;
l'écoulement est torrentiel quand Fr > 1 , dE/ dy < 0 ;
ce cas se trouve dans la branche où V > Vc, y < yc ;
Cas particulier de la section rectangulaire
Le calcul des propriétés hydrauliques critiques (yc, Vc, Ec, Sc) requiert
généralement des calculs itératifs ou l'utilisation d'abaques. La section
rectangulaire est une des rares exceptions où les calculs sont explicites.
Quand l'écoulement est critique, la section d'écoulement Ac se calcule
par la relation suivante :
On obtient :
Où q est le débit par unité de largeur du canal, appelé débit unitaire
On peut calculer l'énergie spécifique critique, la vitesse critique et la
hauteur critique :
Application 8
Dans un canal à section rectangulaire, l'eau s'écoule avec une
profondeur uniforme y = 0,9m. La largeur du canal est B = 2m. La vitesse
d'écoulement est V = 2 m/s.
a) calculer la profondeur critique de l'écoulement
b) calculer la vitesse critique de l'écoulement
c) qualifier l'écoulement
Écoulement critique : cas général
De façon générale, la profondeur critique se calcule par la solution de
l'équation :
Dans cette équation, A et B dépendent de y et la solution doit être
obtenue à l'aide d'une procédure itérative.
Il est impossible de trouver une relation explicite pour calculer la
profondeur critique en fonction des autres variables. On peut utiliser
aussi les courbes de la figure suivante pour calculer la profondeur
critique pour les sections géométriques courantes.
Application 9
Un canal ayant une section de forme trapézoïdale a un fond large de
b = 3,6m. Ses côtés ont une pente de 1. Le débit est Q = 10 m3/s.
Calculer la profondeur critique.
Notion de pente critique
Pour une certaine valeur de la pente, la profondeur normale yn est
égale à la profondeur critique. Cette valeur particulière de la pente du
fond du canal est appelée pente critique Sc.
Il est utile d'obtenir une relation explicite simple pour les canaux et
rivières très larges par rapport à leur profondeur. Dans ce cas, le
périmètre mouillé est :
Application 10
Soit un canal rectangulaire de la largeur B = 30,0m et de coefficient de
Manning n = 0,040. Il faut calculer la pente critique pour les débits Q
suivants : 30,0 m3/s, 60,0 m3/s et 90,0 m3/s.