MECANIQUE DES FLUIDES

APPLIQUEE

Dr YONLI H. Fabien
1
Sommaire

AVANT-PROPOS .................................................................................................................... 5
CHAPITRE I : INTRODUCTION A LA MECANIQUE DES FLUIDES ......................... 7
I. Définition de la mécanique des fluides ............................................................................... 7
II. Définition d’un fluide ........................................................................................................ 7
III. Propriétés des fluides ........................................................................................................ 8
1. Propriétés communes des liquides et des gaz ..................................................................... 8
2. Différences entre les propriétés des fluides ........................................................................ 8
3. Densité d’un fluide ............................................................................................................. 8
4. Viscosité d’un fluide ......................................................................................................... 10
5. Pression de vapeur ............................................................................................................ 11
CHAPITRE II : STATIQUE DES FLUIDES...................................................................... 13
I. Grandeur de base : La pression ......................................................................................... 13
1. Définition de la pression ................................................................................................... 13
2. Isotropie de la pression ..................................................................................................... 13
3. Unités de pression ............................................................................................................. 14
4. Modèles de comportement................................................................................................ 15
II. Equation générale de la statique des fluides .................................................................... 15
1. Cas des fluides incompressibles ....................................................................................... 17
2. Cas où les forces de volume se réduisent à la force de pesanteur .................................... 18
3. Statique des fluides compressibles ................................................................................... 19
III. Théorème de Pascal : Transmission des pressions ......................................................... 19
IV. Théorème d’Archimède .................................................................................................. 20
V. Action de l’eau sur une surface ....................................................................................... 21
1. Action de l’eau sur une surface droite (en position inclinée ou verticale) ....................... 21
2. Action de l’eau sur une surface plane en position horizontale ......................................... 22
3. Action de l’eau sur une surface gauche ............................................................................ 23
4. Action de l’eau sur une surface fermée ............................................................................ 24
5. Détermination du centre de poussée ................................................................................. 25
CHAPITRE III : Cinématique des fluides ........................................................................... 26
I. Vitesse et accélération d’une particule fluide ................................................................... 26
II. Ecoulements permanents et non permanents ................................................................... 26
III. Description lagrangienne et eulérienne du mouvement ................................................. 26

2
 1. Description de Lagrange ................................................................................................... 27
2. Description d’Euler........................................................................................................... 27
IV. Débit volumique – Débit massique ................................................................................ 29
V. Principe de conservation de la masse : équation de continuité ........................................ 30
1. Cas d’un écoulement permanent ....................................................................................... 31
2. Cas d’un fluide homogène incompressible en écoulement permanent ............................. 31
VI. Ecoulement potentiel et écoulement plan ....................................................................... 32
1. Ecoulement potentiel ........................................................................................................ 32
2. Ecoulement plan ............................................................................................................... 32
3. Ecoulement potentiel plan ................................................................................................ 33
4. Lignes de courant et équipotentielles ............................................................................... 33
5. Ecoulements dans les milieux poreux .............................................................................. 34
CHAPITRE IV : Dynamique des fluides incompressibles ................................................. 37
I. Equations générales du mouvement .................................................................................. 37
1. Equations de Navier Stokes .............................................................................................. 37
2. Equations d’Euler : Cas des fluides non visqueux ........................................................... 39
3. Théorème de Bernoulli ..................................................................................................... 40
II. Interprétations physiques du théorème de Bernoulli........................................................ 41
1. Interprétation en énergie ................................................................................................... 41
2. Interprétation en graphique ............................................................................................... 42
III. Régimes d’écoulement – Nombre de Reynolds ............................................................. 42
CHAPITRE V : Notions d’hydraulique à surface libre ........................................................ 47
I. Définitions ......................................................................................................................... 47
I. 2. Paramètres géométriques............................................................................................... 47
I.3. Paramètres hydrauliques ................................................................................................ 50
I.4. Notion de charge hydraulique, théorème de Bernoulli .................................................. 50
II. Ecoulement uniforme ...................................................................................................... 52
II.1. Définition ...................................................................................................................... 52
II.2. Etude des écoulements uniformes ................................................................................ 52
II.2.1. Principe fondamental ................................................................................................. 52
II.2.2. Calcul de la pente d’énergie ....................................................................................... 52
II.2.3. Equations de l’écoulement uniforme ......................................................................... 53
II.2.4. Estimation des valeurs de 𝐊𝐒 ..................................................................................... 54

3
II.3. Problèmes concernant les écoulements uniformes ....................................................... 55
II.3.1. Problème de type 1..................................................................................................... 55
II.3.2. Problème de type 2..................................................................................................... 55
II.3.2.1. Résolution par les abaques ...................................................................................... 56
III.3.2.2. Résolution par la méthode de la débitance ............................................................ 58
II.3.3. Problème de type 3..................................................................................................... 59
Section hydrauliquement favorable de forme trapézoïdale .................................................. 59
CHAPITRE VI : Notions d’hydraulique en charge .............................................................. 61
I. Etude des pertes de charge ................................................................................................ 61
II. Problèmes courants en HEC ............................................................................................ 62

4
AVANT-PROPOS

Les fluides sont présents partout. De l’air dans l’atmosphère, à l’eau dans les lacs, fleuves et
rivières en passant par les huiles et le sang dans les veines. Les fluides ont été et demeurent
étudiés par une vaste communauté de chercheurs. Que ce soit au niveau des fondements ou des
applications, beaucoup de problèmes y relatifs nécessite une meilleure compréhension de ces
milieux.

Ce programme de mécanique des fluides appliquée est destiné aux élèves ingénieurs de niveau
« L » du cycle LMD des écoles d’ingénieurs. Il a pour objectif d’établir et de résoudre les
équations de l’hydrodynamisme pour des écoulements courants rencontrés dans le métier
d’ingénieur (génie civil, hydraulique, génie minier). Les applications sont nombreuses, à titre
d’exemple on pourrait citer la détermination des forces de poussée hydrostatique sur un
ouvrage tel qu’un barrage, un parc à résidus, le débit de fuite à travers les digues et les
fondations de ces barrages, l’écoulement dans les conduites (exhaure d’eau souterraine par
exemple), etc.

Le programme comporte quatre chapitres, assortis de séries de problèmes, permettant d'aborder

les bases de la discipline.

Le premier chapitre concerne la définition de la mécanique des fluides ainsi que les propriétés
des fluides.

Le chapitre 2 est consacré à la statique des fluides. Après avoir introduit la notion de pression,
la loi fondamentale de la statique est formulée. Sa définition conduit à celle de l’hydrostatique
ainsi qu’à celle de la poussée d’Archimède.

Le chapitre 3 porte sur la cinématique des fluides qui vise la description du mouvement.
Successivement sont évoquées les deux types de représentation utilisées en mécanique des
fluide, le principe de conservation de la masse et les quelques types d’écoulements (potentiel et
plan).

Le quatrième chapitre aborde les lois de conservations (quantité de mouvement et énergie)

pouvant être très utiles dans certaines situations que rencontre l’ingénieur. Un accent est porté
sur le théorème de Bernoulli appliqués aux fluides parfaits puis aux fluides visqueux.

5
Les chapitres cinq et six introduisent des notions d’hydraulique générale (hydraulique à surface
libre et hydraulique en charge). Concernant l’hydraulique à surface libre, l’accent est porté sur
les écoulements uniformes.

6
CHAPITRE I : INTRODUCTION A LA MECANIQUE DES FLUIDES

I. Définition de la mécanique des fluides

La mécanique des fluides est l’étude du comportement des fluides en repos ou en mouvement.
Elle étudie l’action des forces sur les fluides. Lorsqu’elle s’intéresse à l’étude des fluides au
repos, on parle de statique des fluides. Si elle s’intéresse à l’étude des fluides en mouvement, il
s’agit de la dynamique des fluides. L’objet de la cinématique des fluides est de décrire ce
mouvement tandis que celui de la dynamique des fluides est d’établir la relation entre les forces
agissantes et ce mouvement.

II. Définition d’un fluide

La matière peut se présenter sous trois (3) états : l’état solide, l’état liquide et l’état gazeux.

A l’état solide, la matière présente un aspect cohésif. Les espacements entre les molécules d’un
solide sont bien définis et les liaisons intermoléculaires y sont très élevées : les molécules sont
arrangées sous forme de lattices spécifiques qui limitent leur mouvement.

Les fluides (liquides ou gaz) quant à eux sont des substances susceptibles de s’écouler plus ou
moins facilement sous l’action de forces extérieures. La structure moléculaire d’un fluide est
telle que les espaces qui séparent les différentes molécules sont essentiellement constants
(quelques changements mineurs dus à la pression et à la température) si bien qu’une masse
donnée de fluide occupe un volume fini dans l’espace.

Un fluide est défini comme un milieu continu ne pouvant être maintenu au repos lorsque des
efforts de cisaillement lui sont appliqués. Il est constitué de « particules fluides ». Ces particules
sont des domaines que l’on considère petits à l’échelle macroscopique mais qui comportent un
grand nombre de molécules, ce qui permet de les estimer grands à l’échelle microscopique.

7
III. Propriétés des fluides

Tous les fluides possèdent des caractéristiques permettant de décrire leurs conditions physiques
dans un état donné. On essaie d’exprimer ces caractéristiques, qu’on appelle propriétés du
fluide, au moyen d’un nombre limité d’unités de base.

1. Propriétés communes des liquides et des gaz

Les fluides présentent tous un certain degré de compressibilité.

Ils n’offrent aucune résistance au changement de forme : ils s’écoulent sous l’effet de forces de
cisaillement si petites soient-elles.
La structure moléculaire d’un fluide est telle que les espaces qui séparent les différentes
molécules sont essentiellement constants si bien qu’une masse donnée de fluide occupe un
volume fini dans l’espace. Il en résulte que les fluides n’ont pas de forme propre et prennent
surtout celle du récipient les contenant.

2. Différences entre les propriétés des fluides

Les liquides sont incompressibles contrairement aux gaz. Leur masse volumique est plus
importante que celle des gaz. Les forces intermoléculaires sont extrêmement faibles chez les
gaz, ce qui engendre que les gaz n’ont pas de volume propre. Leurs molécules se déplacent
jusqu’à ce qu’elles heurtent une autre molécule de gaz ou les parois du vase. Par contre les
forces intermoléculaires sont faibles chez les liquides mais suffisamment fortes pour leur
conférer un volume propre.

3. Densité d’un fluide

1° La densité d’un fluide est sa caractéristique principale qui forme un champ de densité scalaire
dans l’espace occupé par ce fluide.

8
2° La masse volumique (𝑘𝑔/𝑚3 ), ou brève densité, est définie ainsi :
𝑑𝑚
𝜌=
𝑑𝑉
où 𝑑𝑚 est la masse totale de toutes les molécules d’une particules fluide de volume 𝑑𝑉.
C’est une mesure de la concentration de la matière, exprimée comme une masse 𝑑𝑚, par
volume unitaire 𝑑𝑉.

3° Le poids volumique (N/m3 ) est défini de la façon suivante :

𝛾 = 𝜌𝑔
Il représente la force de gravité agissant sur la masse par volume unitaire.

4° L’inverse de la masse volumique (m3 /kg) est appelée volume spécifique

𝑣𝑠 = 1/𝜌
Il représente le volume par masse unitaire.

5° La densité spécifique d’un fluide est définie comme étant le rapport entre la masse volumique
d’un fluide et la masse volumique de l’eau dans des conditions standards (𝑝𝑎 = 1 𝑎𝑡𝑚, 𝑇 =
39,8 °𝐶) :
𝑠𝑠 = 𝜌/𝜌𝑒𝑎𝑢

6° Pour un liquide, la masse volumique 𝜌 dépend de la température et de la pression. La

variation de masse volumique 𝑑𝜌 en fonction de la variation de pression dp est :
𝑑𝜌
= 𝜒𝑑𝑝
𝜌
où 𝜒 est le coefficient de compressibilité isothermique. Cette relation montre qu’une
augmentation de pression dp correspond à une augmentation de la masse volumique 𝑑𝜌.
Dans les conditions normales, on a pour l’eau 𝜒 = 5. 10−10 𝑃𝑎−1 .

7° Pour un gaz, la masse volumique 𝜌 dépend de la température et de la pression, ce qui est

exprimé par la relation suivante, qu’on appelle loi du gaz :
𝑝′
= 𝑅𝑇
𝜌𝑛
où 𝑝′ est la pression absolue, 𝑇 la température absolue et 𝑅 la constante du gaz. Pour un gaz
parfait, l’exposant vaut 𝑛 = 1. Pour l’air à la pression atmosphérique, 𝑅 = 287 𝑚2 /(𝑠 2 𝐾°).

9
4. Viscosité d’un fluide

1° La viscosité est l’aptitude d’un fluide à s’écouler plus ou moins facilement. Comme les
frottements dans un solide, la viscosité s’oppose à l’écoulement du fluide. La viscosité est une
propriété intensive qui permet de distinguer les fluides parfaits (dépourvus de viscosité) des
fluides réels/visqueux.

répartition
de vitesse
dz

z
z dV
x

2° Dans un fluide, la force de frottement (force tangentielle) par unité de surface ou tension de
frottement 𝜏𝑧𝑥 est proportionnelle au gradient de vitesse de l’écoulement 𝑑𝑉/𝑑𝑧. La loi de la
viscosité de Newton ou de Stokes s’exprime donc ainsi :
𝑑𝑉
𝜏𝑧𝑥 ≈
𝑑𝑧
Lorsque le fluide se déplace en couches parallèles (écoulement dit laminaire), le facteur de
proportionnalité est le coefficient de viscosité  ; on écrit alors :
𝑑𝑉
𝜏𝑧𝑥 = 
𝑑𝑧
Telle est la définition de la viscosité dynamique ou absolue . La tension de frottement 𝜏𝑧𝑥 agit
dans le sens opposé à l’écoulement.
La viscosité cinématique 𝜈 est définie comme étant le rapport entre la viscosité dynamique  et
la masse volumique 𝜌 :

𝜈=
𝜌

10
4° Un fluide visqueux présente des forces intermoléculaires plus élevées que celle d’un fluide
moins visqueux.
Ex : L’eau coule plus vite que l’huile (l’huile est plus visqueuse).

5° Dans le système SI, la viscosité dynamique  a pour unité 𝑘𝑔/𝑚. 𝑠 soit encore l’unité 𝑃𝑎. 𝑠.
Il existe également la poise et le Poiseuille.
1 𝑃𝑜𝑖𝑠𝑒 = 1 𝑔/𝑐𝑚. 𝑠.
1 𝑃𝑜𝑖𝑠𝑒𝑢𝑖𝑙𝑙𝑒 (𝑃𝑖) = 1 𝑃𝑎. 𝑠 = 10 𝑃𝑜𝑖𝑠𝑒(𝑃𝑜).
Ces anciennes unités sont données en hommage à Jean-Louis Marie Poiseuille (1797-1869),
physicien et médecin ayant travaillé sur les fluides visqueux dans les tubes capillaires.
La viscosité cinématique a pour unité dans le système SI 𝑚2 /𝑠 ; dans le système CGS, l’unité
est le stokes.
1 𝑠𝑡𝑜𝑘𝑒𝑠 = 1 𝑐𝑚2 /𝑠.

6° La viscosité est fonction de la température. L’eau à 20° C à la pression atmosphérique a une

viscosité cinématique 𝜈 = 10−6 𝑚2 /𝑠.

7° Les fluides qui obéissent à la loi de viscosité de Newton

𝑑𝑉
𝜏𝑧𝑥 =  𝑑𝑧 sont appelés fluides newtoniens. L’eau et d’autres liquides ainsi que l’air

appartiennent à cette catégorie. Certains fluides tels que les plastiques liquides ou les
suspensions concentrées ne se comportent pas selon cette loi et sont appelés fluides non
newtoniens. L’étude des propriétés de ces fluides constituent le domaine de la rhéologie.

5. Pression de vapeur

1° Tous les liquides ont tendance à s’évaporer, la phase liquide se transforme en phase gazeuse.
Les molécules de vapeur (gaz) exercent une pression, appelée pression de vapeur pv . Les unités
dans le système SI sont le Pa. Cette pression varie pour chaque fluide et croît en fonction de la
température.

2° L’eau bout lorsque la pression absolue ou imposée devient inférieure ou égale à la pression
de vapeur. Le point d’ébullition dépend donc de la pression imposée ainsi que de la température.

11
3° Si dans un liquide en mouvement la pression absolue devient inférieure à la pression de
vapeur pv , il se produit une évaporation. Le liquide en question contient des bulles de vapeur
en suspension ; c’est le phénomène de cavitation.

Propriétés physiques de quelques fluides

Fluide 𝝆 (𝒌𝒈/𝒎𝟑 ) 𝒔𝒔 (−) 𝝂 (𝒎𝟐 /𝒔) 𝒑𝒗 (𝑵/𝒎𝟐 )

Eau (1 atm, 20° C) 998,2 0,9982 1,004.10−6 2339,0
Eau de mer (1 atm, 1025,0 1,0250 1,2.10−6 -
15 °C)
Mercure (1 atm, 13555,0 13,5550 1,1.10−7 0,17
15,6° C)
Alcool éthylique 788,6 0,7910 1,5.10−6 5860,0
(96% dans l’eau)
Glycérine (1atm, 1257,6 1,2600 1,2.10−3 0,014
20 °C)
Air (20 °C) 1,21 0,0012 1,49.10−5 -
Hélium (1atm, 0,178 0,0002 1,04.10−4 -
0°C)

12
CHAPITRE II : STATIQUE DES FLUIDES
La statique des fluides est la science qui étudie les fluides au repos.

I. Grandeur de base : La pression

1. Définition de la pression

La pression est le rapport de l’intensité de la force F exercée perpendiculairement à une

surface sur l’aire S de cette surface. Ainsi, la pression P est donnée par : P = F / S.
Autour d’un point M de la surface S, isolons un élément de surface ΔS sur lequel s’exerce
perpendiculairement la force ΔF. La pression au point M est la limite de ΔF/ΔS quand ΔS
tend vers 0.
∆𝐹 𝑑𝐹
𝑃 = lim ;𝑃=
𝛥𝑆→0 ∆𝑆 𝑑𝑆
Considérant la normale extérieure unitaire 𝑛⃗ d’une surface, la force pressante s’écrit
𝐹 = −𝑃𝑆𝑛⃗.

2. Isotropie de la pression

Un liquide au repos est isotrope pour ce qui concerne la pression. Cela signifie qu'en un point
quelconque d'un fluide, la pression en ce point est identique dans toutes les directions.
En effet, soit un élément de surface dS’ de centre M’ et un cylindre infiniment petit de section
droite dS’. L’autre base dS, de centre M, est d’orientation quelconque définie par l’angle 𝛼. Les
dimensions linéaires de ce cylindre, donc MM’, sont infiniment petites du premier ordre.

Nous avons par définition :

⃗⃗⃗ = −𝑝′𝑑𝑆′𝑛′
𝑑𝐹 = −𝑝𝑑𝑆𝑛⃗ et 𝑑𝐹′ ⃗⃗⃗

13
En écrivant l’équilibre des forces agissant sur ce cylindre élémentaire, on peut négliger les
forces de volume (infiniment petits d’ordre 3) devant les forces de surface (infiniment petits
d’ordre 2) et on aura :
⃗⃗⃗ + ∑ 𝑑𝐹𝑙𝑎𝑡 = ⃗0
𝑑𝐹 + 𝑑𝐹′

où 𝑑𝐹𝑙𝑎𝑡 sont les forces de surface latérale du cylindre.

La projection de cette équation dans la direction ⃗⃗⃗

𝑛′ donne :
𝑝𝑑𝑆𝑐𝑜𝑠𝛼 = 𝑝′𝑑𝑆′
⃗⃗⃗ , elles s’éliminent.
Puisque les forces de surface agissant sur les faces latérales sont normales à 𝑛′
Comme 𝑑𝑆𝑐𝑜𝑠𝛼 = 𝑑𝑆’, on aura 𝑝 = 𝑝′.
Cette égalité est vraie quelque soit 𝛼 ; la pression est donc indépendante de l’orientation 𝑑𝑆.

3. Unités de pression

Dans le système international d’unités, la pression s’exprime en pascal (Pa): 𝟏𝑷𝒂 = 𝟏 𝑵/𝒎𝟐 .
Toutefois, il s’agit d’une unité très petite qui s'utilise surtout dans les calculs en raison de son
appartenance au système SI.
Il existe d’autres unités de pression, les plus courantes étant l’atmosphère (atm), le bar, le
« millimètre de mercure ».
Multiples du Pascal :
1𝑘𝑃𝑎 = 103 𝑃𝑎
{1𝑀𝑃𝑎 = 106 𝑃𝑎 .
1𝐺𝑃𝑎 = 109 𝑃𝑎
Autres unités
Le bar
1𝑏𝑎𝑟 = 105 𝑃𝑎.
L’atmosphère
1𝑎𝑡𝑚 = 1,01325. 105 𝑃𝑎.
Le kilogramme force /cm2
1𝑘𝑔. 𝑓/𝑐𝑚2 = 9,81 × 104 𝑃𝑎.
Le millimètre de mercure
105 𝑃𝑎 = 760 𝑚𝑚 𝑑𝑒 𝑚𝑒𝑟𝑐𝑢𝑟𝑒.

14
4. Modèles de comportement

Lorsqu'on modifie la pression d'un fluide, sa masse volumique en est affectée, dans une mesure
plus ou moins grande suivant la nature du fluide. Cette variation est traduite par la loi de
comportement, c'est-à-dire en définitive par l'équation d'état.
Les liquides très peu sensibles aux variations de pression ont leur comportement en général
bien décrit par le modèle incompressible 𝜌 = 𝑐𝑡𝑒.
Au contraire, les gaz très sensibles aux variations de pression ont leur comportement dépendant
de l'environnement thermique. Deux modèles extrêmes sont envisagés :
1° Le modèle isotherme, dans lequel le contact du fluide avec son environnement est
supposé parfait ; alors le gaz reste à une température constante.
2° Le modèle adiabatique, qui suppose au contraire que le gaz n'échange aucune chaleur avec
l'extérieur.

II. Equation générale de la statique des fluides

Dans cette partie de l’étude nous appliquons les principes fondamentaux de la mécanique à un
fluide considéré comme un milieu continu que nous découpons en éléments de volume
infiniment petits.

Considérons donc un élément de fluide de forme parallélépipédique dont les arêtes dx, dy et
dz sont parallèles aux axes. Cet élément est en équilibre sous l’action des forces suivantes :

- Les forces de volume ou forces massiques ρdxdydzΓ où Γ est la force extérieure agissant sur
la masse fluide et rapportée à l’unité de masse (densité massique des forces de volume). Le
fluide étant un milieu matériel d’élément de volume 𝑑𝑣, il possède une masse volumique et
donc une masse 𝑑𝑚 = 𝜌𝑑𝑣.

Les forces de volume sont dues au fait que la masse fluide est située dans un champ caractérisé
en chaque point par un vecteur généralement variable d’un point à l’autre. Il peut s’agir des
champs électriques, magnétiques, d’inertie. Mais le cas le plus fréquent est sans doute celui du
champ de pesanteur caractérisé par le vecteur 𝑔 que l’on peut considérer uniforme pour les
applications courantes.

15
- Les forces de contact ou forces de surface qui se réduisent aux forces de pression en
hydrostatique (cas des fluides parfaits ; pour les fluides réels, les forces dues au frottement
encore appelés tensions tangentielles doivent être prises en compte). Ces forces de contact sont
exercées par les particules fluides extérieures à l’élément sur la surface ds qu’il délimite. Elles
s’exercent donc aux frontières du domaine. Ces forces sont normales aux 6 faces du
parallélépipède donc parallèles aux axes.

Soit P=P(x,y,z) la pression au centre de l’élément de volume de coordonnées M(x,y,z).

Considérons les forces de pression perpendiculaires aux faces dydz : ces grandeurs
𝑃(𝑥 − 𝑑𝑥/2, 𝑦, 𝑧) et 𝑃(𝑥 + 𝑑𝑥/2, 𝑦, 𝑧) sont les pressions en deux points infiniment voisins de
𝑀.

𝜕𝑃 𝑑𝑥
A gauche, elles valent 𝑃(𝑥 − 𝑑𝑥/2, 𝑦, 𝑧)𝑑𝑦𝑑𝑧𝑖 = (𝑃 − 𝜕𝑥 ) 𝑑𝑦𝑑𝑧𝑖.
2
𝜕𝑃 𝑑𝑥
A droite, elles valent −𝑃(𝑥 + 𝑑𝑥/2, 𝑦, 𝑧)𝑑𝑦𝑑𝑧𝑖 = − (𝑃 + 𝜕𝑥 ) 𝑑𝑦𝑑𝑧𝑖.
2

En additionnant ces deux quantités, on obtient la composante des forces de pressions

𝜕𝑃
Suivant x : − 𝜕𝑥 𝑑𝑥𝑑𝑦𝑑𝑧

16
Par un raisonnement semblable suivant les autres axes, on obtient :
𝜕𝑃
Suivant y : − 𝜕𝑦 𝑑𝑥𝑑𝑦𝑑𝑧
𝜕𝑃
Suivant z : − 𝜕𝑧 𝑑𝑥𝑑𝑦𝑑𝑧

La résultante des forces de contact s’écrit donc :

𝜕𝑃 𝜕𝑃 𝜕𝑃
(− 𝜕𝑥 𝑑𝑥𝑑𝑦𝑑𝑧; − 𝜕𝑦 𝑑𝑥𝑑𝑦𝑑𝑧; − 𝜕𝑧 𝑑𝑥𝑑𝑦𝑑𝑧).

L’équilibre statique ∑ 𝐹 = ⃗0 donne

𝜕𝑃 𝜕𝑃 𝜕𝑃
𝜌𝑑𝑥𝑑𝑦𝑑𝑧𝛤 − ( 𝑖+ 𝑗+ ⃗ ) 𝑑𝑥𝑑𝑦𝑑𝑧 = ⃗0
𝑘
𝜕𝑥 𝜕𝑦 𝜕𝑧
Cette équation peut s’écrire en reconnaissant l’expression du gradient de P.

⃗ = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝝆𝜞 𝒈𝒓𝒂𝒅𝑷 (1)

Il s’agit là de l’équation fondamentale de la statique des fluides.

Il est également nécessaire de tenir compte des conditions aux limites. Pour un fluide présentant
une surface libre sur laquelle la pression donnée est égale à la pression atmosphérique, alors il
faut que la pression à l’intérieur du fluide P(x, y, z) se raccorde tout au long de cette surface
avec la pression atmosphérique donnée. Dans le cas contraire, il y aurait discontinuité.
Remarque : Aucune hypothèse sur la nature du fluide, ni sur les forces de volume n’a encore
été émise.

1. Cas des fluides incompressibles

Pour les fluides incompressibles tels que l’eau, la masse volumique ne varie pas.
L’équation (1) devient alors :
𝑃
𝛤 = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑔𝑟𝑎𝑑 (2)
𝜌
Dans le cas fréquent où la densité massique des forces de volume dérive d’un potentiel U,
on a :
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑈 (3)
𝛤 = −𝑔𝑟𝑎𝑑
𝑃
Une identification avec (2) donne – 𝑈 = 𝜌 + 𝑐𝑡𝑒 (4).

17
2. Cas où les forces de volume se réduisent à la force de pesanteur

Dans ce cas 𝛤 = 𝑔

0
𝛤=( 0 )
−𝑔
L’équation (3) devient
𝜕𝑈
−
𝜕𝑥
0 𝜕𝑈
( 0 )= − ⇒ 𝑈 = 𝑔𝑧
𝜕𝑦
−𝑔
𝜕𝑈
−
( 𝜕𝑧 )
𝑃
L’équation de la statique nous donne maintenant −𝑔𝑧 = 𝜌 + 𝑐𝑡𝑒

Soit 𝑃 + 𝜌𝑔𝑧 = 𝑐𝑡𝑒

En hydrostatique donc, la quantité 𝑃 + 𝜌𝑔𝑧 est une constante : théorème fondamental de
l’hydrostatique.
La quantité 𝑃 + 𝜌𝑔𝑧 est souvent appelée la pression étoilée 𝑃 ∗ .
𝑃∗ = 𝐶𝑡𝑒 .

Les hypothèses de base de ce théorème sont les suivantes :

Le fluide considéré :
Est isotherme ;
Est incompressible (sa masse volumique est indépendante de la pression mais surtout de la
température) ;
Occupe un domaine d’étendue réduite dans lequel il est en équilibre et soumis uniquement à
l’action de la pesanteur.

Conséquences
Les surfaces équipotentielles sont des plans horizontaux (parce que 𝑈 = 𝑔𝑧 = 𝑐𝑡𝑒 ⇔ 𝑧 =
𝑐𝑡𝑒).
Les équipotentielles et les isobares sont confondues (U = cte ⇔ 𝑃 = 𝑐𝑡𝑒). La pression ne
varie pas dans un plan horizontal mais elle varie sur le plan vertical.

18
La surface libre d’un liquide au repos est donc toujours horizontale (parce que la surface libre
est une isobare établie à la pression atmosphérique).

Application : Etablir que la différence de pression entre 2 points ne dépend que de la différence
de leurs côtes et de la masse volumique ou du poids volumique.

Déterminer la pression en tout point d’altitude z du réservoir suivant. Il contient trois

fluides non miscibles de masse volumiques ρ1, ρ2, ρ3.

3. Statique des fluides compressibles

Dans le cas des fluides compressibles, la masse volumique n’est pas constante. L’équation de
la statique à elle seule est insuffisante pour déterminer les inconnues 𝑃 et 𝜌. Si la
thermodynamique fournit une relation directe et biunivoque entre P et 𝜌 (hypothèse de fluide
barotrope), alors 𝑑𝑃 = −𝜌𝑔𝑑𝑧 et 𝑃 = 𝑓(𝜌) sont les équations permettant de résoudre le
problème.

III. Théorème de Pascal : Transmission des pressions

Toute variation de pression en un point d’un liquide en équilibre est transmise intégralement à
tous les autres points du même liquide. Cette proposition est connue sous le nom du principe
de Pascal. En effet, soient deux points A et B d’un liquide de masse volumique ρ en équilibre.
Si on modifie la pression au point A de ∆𝑃, celle-ci devient 𝑃𝐴 ′ = 𝑃𝐴 + ∆𝑃.

19
Ayant 𝑃𝐵 ′ − 𝑃𝐴 ′ = 𝜌𝑔ℎ = 𝑃𝐵 − 𝑃𝐴 ,
On obtient 𝑃𝐵 ′ = 𝑃𝐵 + (𝑃𝐴 ′ − 𝑃𝐴 ) = 𝑃𝐵 + ∆𝑃.

Blaise Pascal (1623 –1662) est un mathématicien, physicien,

inventeur, philosophe, moraliste et théologien français.
Le théorème de Pascal ou principe de Pascal avait déjà été
connu de Simon Stevin. On lui donne aussi le nom de
paradoxe hydrostatique. Pascal donna cependant aux travaux
de Stevin une forme belle et ordonnée à travers ses exposés.

Blaise Pascal (1623-1662)

IV. Théorème d’Archimède

Tout corps plongé dans l’eau subit de la part de ce fluide une poussée d’Archimède verticale
dirigée de bas en haut et égale au poids du fluide déplacé. Cette force est appliquée au centre
de gravité du fluide déplacé.

Mathématicien, physicien, ingénieur, Archimède (287 av. J.-

C. – 212 av. J.C.) est un scientifique grec considéré comme
l’un des plus grands savants de l’antiquité.
C’est la poussée d’Archimède qui « permet » à une
montgolfière de flotter dans l’air, de « flotter » sur la mer
morte, qui explique qu’un plongeur coule à une certaine
profondeur.

Archimède (-287 à -212) Le traité des corps flottants d’Archimède jeta les bases de ce
qui sera plus tard la science nommée hydrostatique.

20
V. Action de l’eau sur une surface

Soit une surface de géométrie quelconque immergée dans un liquide. En général on demande
de répondre aux trois questions suivantes :
1° Quelle est l’intensité de la force qui s’applique sur cette surface ?
2° Où est situé le point d’application de cette force ?
3° Quelle est sa direction ?
La réponse à ces questions diffère selon la géométrie de ces surfaces.

1. Action de l’eau sur une surface droite (en position inclinée ou verticale)

Considérons une surface plane quelconque (vanne par exemple) dont la trace sur une coupe
verticale est le segment AB. Cette surface est inclinée d’un angle α.

α
Cy
y

F1 F2

F
u

La résultante des forces élémentaires exercées par le fluide est appelée poussée.

Soit 𝐹 cette force. Elle est la somme des forces agissant sur la face 1 et sur la face 2.

Sur la face 1, il s’exerce la force ⃗⃗⃗⃗ ⃗ 𝒅𝑺 = ∫𝑺 𝒑𝒅𝑺 𝒏

𝑭𝟏 = ∫𝑺 𝒑𝒏 ⃗
Si la surface libre est l’origine des z, la loi de l’hydrostatique donne :

21
𝑷 + 𝝆𝒈𝒛 = 𝒄𝒕𝒆 = 𝑷𝒂𝒕𝒎

Il en résulte que ⃗⃗⃗⃗ ⃗

𝑭𝟏 = ∫𝑺 (𝑷𝒂𝒕𝒎 − 𝝆𝒈𝒛)𝒅𝑺 𝒏

Sur la face 2, il s’exerce la force ⃗⃗⃗⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝒅𝑺 = − ∫ 𝒑𝒂𝒕𝒎 𝒅𝑺 𝒏

𝑭𝟐 = ∫𝑺 𝒑𝒂𝒕𝒎 (−𝒏) ⃗
𝑺

Il en résulte que la force nette ou poussée hydrostatique vaut :

⃗𝑭 = ⃗⃗⃗⃗
𝑭𝟏 + ⃗⃗⃗⃗ ⃗ = −𝝆𝒈 ∫ 𝒛𝒅𝑺 𝒏
𝑭𝟐 = − ∫ 𝝆𝒈𝒛𝒅𝑺𝒏 ⃗⃗
𝑺 𝑺

On a ∫𝑺 𝒛𝒅𝑺 = 𝒛𝑮 𝑺 où 𝑧𝐺 désigne la profondeur du centre d’inertie de la surface plane.

Connaissant 𝒖𝑮 on peut retrouver 𝒛𝑮 par la formule 𝒛𝑮 = −𝒖𝑮 𝒔𝒊𝒏𝜶.
En définitive la poussée s’écrit :
⃗ = −𝝆𝒈𝒛𝑮 𝑺𝒏
𝑭 ⃗

Cette relation a été établie dans le cas où la surface plane est à l’interface entre l’eau et
l’atmosphère.
Pour une surface plane immergée, de façon générale, il faudra tenir compte de la pression
atmosphérique et l’action de l’eau vaut :
⃗𝑭 = (𝑷𝒂𝒕𝒎 − 𝝆𝒈𝒛𝑮 )𝑺𝒏
⃗⃗

2. Action de l’eau sur une surface plane en position horizontale

Considérons une paroi de surface S immergée horizontalement à une profondeur h. La pression

est la même sur la surface car on se situe à z = cste.
Dans ce cas zG = −h en se référant à la démonstration faite en 1 et donc :
⃗ = ρghSn
F ⃗.
La force s’applique normalement à la surface.

22
 h

S
⃗
n

Quelle que soit la forme des réservoirs, s’ils sont remplis du même liquide, jusqu’à la même
hauteur ℎ, les fonds de même surface 𝑆 sont soumis à la même force de pression.

S S S

3. Action de l’eau sur une surface gauche

Lorsqu’on a une telle surface, les principes énoncés précédemment sont valables pour les
projections de cette surface suivant les plans verticaux et suivant z la poussée vaut :

𝑭𝒛 = 𝝆𝒈𝑽∗

23
 z

V*

4. Action de l’eau sur une surface fermée

𝑭𝒛 = 𝝆𝒈𝑽

Où V est le volume d’eau déplacé : Principe d’Archimède.

24
5. Détermination du centre de poussée

Pour certains calculs en structure (barrages, vannes, etc.), on a besoin du point d’application de
la poussée appelé centre de poussée. Sa position s’obtient en écrivant que le moment des forces
élémentaires par rapport à un point ou une droite est égal au moment de la poussée.

Dans le cas d’une surface plane à l’interface entre l’eau et l’atmosphère, en considérant le
moment par rapport à la droite 𝑂𝑦.
Soit P la position du centre de poussée :
Le moment résultant par rapport à cet axe

⃗⃗⃗⃗⃗ ˄𝐹 = ∫ 𝑂𝑀
𝑂𝑃 ⃗⃗⃗⃗⃗⃗ ˄𝑑𝐹
𝑆

⃗⃗⃗⃗⃗⃗ ˄𝑑𝐹 = ∫ 𝑢(−𝜌𝑔𝑧𝑑𝑆) = −𝜌𝑔 ∫ 𝑢𝑧𝑑𝑆

∫ 𝑂𝑀
𝑆 𝑆 𝑆

Il doit être équilibré par le moment de la poussée qui vaut

−𝜌𝑔𝑧𝐺 𝑆𝑢𝑃 ;

Après simplification, on obtient donc 𝑧𝐺 𝑆𝑢𝑃 = ∫𝑆 𝑢𝑧𝑑𝑆

En remarquant que 𝑧 = −𝑢𝑠𝑖𝑛𝛼 𝑒𝑡 𝑧𝐺 = −𝑢𝐺 𝑠𝑖𝑛𝛼, on peut encore écrire :

𝑢𝐺 𝑆𝑢𝑃 = ∫ 𝑢2 𝑑𝑆
𝑆

où le terme de droite est le moment d’inertie 𝐼𝑜𝑦 de la surface par rapport à l’axe 𝑂𝑦.

Le théorème de Huygens permet d’obtenir le moment d’inertie par rapport à 𝐺𝑦.

𝐼 = 𝐼𝐺𝑦 + 𝑆𝑢𝐺 2 .

Et par suite 𝑢𝐺 𝑆𝑢𝑃 = 𝐼𝐺𝑦 + 𝑆𝑢𝐺 2.

𝐼𝐺𝑦
𝑢𝑃 − 𝑢𝐺 =
𝑆𝑢𝐺

Cette dernière relation montre que le centre de poussée n’est pas confondu avec le centroïde G
il est toujours situé en dessous. On peut aussi écrire en utilisant z,

𝐼𝐺𝑦
𝑧𝑃 − 𝑧𝐺 = 𝑠𝑖𝑛2 𝛼
𝑆𝑧𝐺

25
CHAPITRE III : Cinématique des fluides

La cinématique est l’étude du mouvement des particules fluides sans tenir compte des forces
qui causent le mouvement.

I. Vitesse et accélération d’une particule fluide

Les lois de la physique newtonienne reposant sur les notions de vitesse et d’accélération de
particules, on introduit la notion de particule fluide qui est un petit élément de masse (petit à
l’échelle de l’ingénieur mais très grand à l’échelle moléculaire).
La vitesse d’une particule fluide est le taux de variation temporelle du vecteur position de cette
particule.
⃗⃗⃗⃗⃗⃗
𝑑𝑂𝑀
⃗ =
𝑉
𝑑𝑡
⃗ : vecteur vitesse
V
⃗⃗⃗⃗⃗⃗
𝑂𝑀: vecteur position de la particule fluide situé au point M.
L’accélération d’une particule fluide est le taux de variation temporelle du vecteur vitesse de
cette particule.
𝑑𝑉⃗
𝛾=
𝑑𝑡
𝛾 : vecteur accélération.

II. Ecoulements permanents et non permanents

Un mouvement est dit permanent lorsque la vitesse ne dépend du temps de façon explicite.
C'est-à-dire que la vitesse de la particule qui se trouve en un point donné reste la même à des
instants différents. Il s’en suit en général que les autres variables de l’écoulement sont
indépendantes du temps.

Dans le cas contraire, le mouvement est dit non permanent.

III. Description lagrangienne et eulérienne du mouvement

Le mouvement d’un fluide est repérable au moyen de variables adéquates ; deux méthodes sont
considérées : les variables de Lagrange et les variables d’Euler.

26
1. Description de Lagrange

La description lagrangienne du mouvement consiste à suivre individuellement une particule

fluide dans son mouvement.
Soient (𝑥0 ,𝑦0 ,z0 ) les coordonnées de cette particule à un instant initial t0. A l’instant t, sa
position est définie par ses coordonnées x,y,z en fonction de 𝑥0 ,𝑦0 ,z0 et t.
𝑥 = 𝑓(𝑥0 , 𝑦0 , z0 , t) 𝑦 = 𝑔(𝑥0 , 𝑦0 , z0 , t) 𝑒𝑡 𝑧 = ℎ(𝑥0 , 𝑦0 , z0 , t) (∗)

Le lieu géométrique des positions successives d’une particule au cours du temps s’appelle sa
trajectoire.
Les équations paramétriques (en t) de la trajectoire sont donc données par les relations (*).
On peut en déduire la vitesse de cette particule à tout instant t :
𝑑𝑥 𝜕𝑓
𝑉𝑥 (𝑥0 , 𝑦0 , z0 , 𝑡) = =
𝑑𝑡 𝜕𝑡
𝑑𝑦 𝜕𝑔
𝑉𝑦 (𝑥0 , 𝑦0 , z0 , 𝑡) = =
𝑑𝑡 𝜕𝑡
𝑑𝑧 𝜕ℎ
𝑉𝑧 (𝑥0 , 𝑦0 , z0 , 𝑡) = =
𝑑𝑡 𝜕𝑡

Pt (x, y, z)
trajectoire

Pt0 (x0 , y0 , z0 )

2. Description d’Euler

La description eulérienne du mouvement consiste à considérer un point fixe M de l’espace et à

étudier, en fonction du temps, ce qui se passe en ce point. Nous déterminerons ainsi, en fonction
⃗ de la particule qui se trouve en ce point M de coordonnées x, y, z. Nous
du temps, la vitesse 𝑉
aurons alors les composantes de la vitesse en chaque instant t :
𝑉𝑥 = 𝐹(𝑥, 𝑦, 𝑧, 𝑡)

27
𝑉𝑦 = 𝐺(𝑥, 𝑦, 𝑧, 𝑡)
𝑉𝑧 = 𝐻(𝑥, 𝑦, 𝑧, 𝑡)
Les x,y,z,t sont appelées variables d’Euler.
On appelle ligne de courant une courbe tangente en chacun de ses points au vecteur vitesse en
ce point. Mathématiquement, cela se traduit par :
⃗⃗⃗⃗⃗⃗ = 𝜆𝑉
𝑑𝑂𝑀 ⃗ ou 𝑑𝑂𝑀
⃗⃗⃗⃗⃗⃗ ⌃𝑉
⃗ = ⃗0
𝑑𝑥 𝑑𝑦 𝑑𝑧
= =
𝑉𝑥 𝑉𝑦 𝑉𝑧
Les équations des lignes de courant s’obtiennent donc par intégration des équations ci-dessus.

⃗ (Vx , Vy , Vz )
V
ligne de
courant
M(x,y,z)
)

D’une façon générale, la représentation d’Euler, plus simple que celle de Lagrange, est utilisée
plus fréquemment.

28
 Jean Louis Lagrange (1736 - 1813) est un mathématicien,
mécanicien et astronome italien naturalisé français. Il a établi
de nombreuses correspondances avec Leonhard Euler (1707
–1783) certainement le plus grand mathématicien de son
époque et physicien.
Jean Louis Lagrange (1736 Bien en avance sur son temps, en mécanique des fluides,
- 1813) Lagrange introduit le concept de potentiel de vitesses pour
tout écoulement fluide réel, le potentiel de courant, pour un
fluide incompressible.
Euler fut le premier à poser les équations désormais connues
sous le nom d'équations d'Euler des fluides parfaits, dans
« Mémoires de l'Académie royale des sciences et des belles
lettres de Berlin ».

Leonhard Euler (1707 –

1783)

IV. Débit volumique – Débit massique

Le débit volumique est la quantité de volume de fluide qui traverse une section donnée par
unité de temps.
Considérons un élément de surface dS traversé par le fluide et normal à la direction de la vitesse.
Dans un intervalle de temps dt, la quantité de fluide passée à travers dS est contenue dans un
cylindre de longueur Vdt et de base dS.

dS

Vdt

29
Son volume vaut alors :
𝑑𝑣 = 𝑉𝑑𝑡𝑑𝑆
Le débit entrant à travers la surface élémentaire dS est donc égal à :
𝑑𝑣
𝑞𝑣 = 𝑑𝑄 = = 𝑉𝑑𝑆
𝑑𝑡
Ceci est valable dans le cas où 𝑉 ⃗ et 𝑛⃗ sont colinéaires.
Dans le cas contraire, la longueur du cylindre évoqué vaut VdtcosΘ (Ѳ étant l’angle entre la
vitesse et le vecteur unitaire 𝑛⃗ normal à la surface 𝑑𝑆) et
𝑑𝑣
𝑞𝑣 = 𝑑𝑄 = = 𝑉𝑐𝑜𝑠𝛩𝑑𝑆 = 𝑉 ⃗ . 𝑛⃗𝑑𝑆
𝑑𝑡
Le débit massique ou débit en masse 𝑑𝑚̇ est la quantité de masse qui traverse dS par unité de
temps :
𝑞𝑚 = 𝑑𝑚̇ = 𝜌𝑞𝑣 = 𝜌𝑉 ⃗ . 𝑛⃗𝑑𝑆
A travers une section de grandeur finie S, le débit s’obtient par intégration de vitesse.
⃗ . 𝑛⃗𝑑𝑆
𝑄 = ∫ 𝑑𝑄 = ∫ 𝑉
𝑆 𝑆

⃗ . 𝑛⃗𝑑𝑆
𝑚̇ = ∫ 𝑑𝑚̇ = ∫ 𝜌𝑉
𝑆 𝑆

V. Principe de conservation de la masse : équation de continuité

L’équation de continuité, qui est une équation fondamentale de la mécanique des fluides,
exprime le principe de conservation de la masse.

Etablissement de l’équation
Soit un milieu fluide continu sans source ni puits. On considère l’élément de volume de
dimension dx, dy et dz. La projection sur le plan est donnée par la figure ci-dessous.

Le volume de fluide entrant par la face AB dans le temps dt est Vxdydzdt correspondant à la
masse ρVxdydzdt.

30
 𝜕(𝜌𝑉𝑥 )
De même, la masse de fluide qui sort par BC est [𝜌𝑉𝑥 + 𝑑𝑥] 𝑑𝑦𝑑𝑧𝑑𝑡.
𝜕𝑥

La perte de masse par ces deux faces est la différence entre la masse sortante et la masse
entrante.
𝜕𝜌𝑉𝑥
Elle est égale à 𝑑𝑥𝑑𝑦𝑑𝑧𝑑𝑡.
𝜕𝑥

En faisant le calcul pour les autres directions et pour l’ensemble des six (6) faces, la masse
𝜕(𝜌𝑉𝑥 ) 𝜕(𝜌𝑉𝑦 ) 𝜕(𝜌𝑉𝑧 )
perdue est 𝑑𝑚 = [ + + ] 𝑑𝑥𝑑𝑦𝑑𝑧𝑑𝑡.
𝜕𝑥 𝜕𝑦 𝜕𝑧

Or la masse que le parallélépipède contient à l’instant t est bien 𝜌𝑑𝑥𝑑𝑦𝑑𝑧 et sa diminution est
𝜕𝜌
𝑑𝑚 = − 𝜕𝑡 𝑑𝑡𝑑𝑥𝑑𝑦𝑑𝑧.

En égalant ces deux expressions de dm, on obtient

𝜕𝜌 𝜕(𝜌𝑉𝑥 ) 𝜕(𝜌𝑉𝑦 ) 𝜕(𝜌𝑉𝑧 )
+ + + = 0 ∶ Equation de continuité
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
Cette équation peut s’écrire encore :
𝜕𝜌
⃗)=0
+ 𝑑𝑖𝑣(𝜌𝑉
𝜕𝑡

1. Cas d’un écoulement permanent

𝜕𝜌
La masse volumique est indépendante du temps donc = 0 et l’équation devient
𝜕𝑡
⃗)=0
𝑑𝑖𝑣(𝜌𝑉

2. Cas d’un fluide homogène incompressible en écoulement permanent

La masse volumique est constante en tout point de l’espace.

⃗ ) = 0 devient 𝑑𝑖𝑣(𝑉
𝑑𝑖𝑣(𝜌𝑉 ⃗)=0

Si l’on considère un volume élémentaire 𝑑𝑉 = 𝑑𝑥𝑑𝑦𝑑𝑧, qu’on le multiplie par l’équation

précédente et qu’on intègre cette quantité par rapport au volume, on obtient :

⃗ 𝑑𝑉 = 0
∫ 𝑑𝑖𝑣 𝑉
𝑉

31
Selon le théorème de Gauss, on peut transformer une intégrale de volume en une intégrale de
surface (fermée) :

⃗ 𝑑𝑉 = ∫ 𝑉
∫ 𝑑𝑖𝑣 𝑉 ⃗ 𝑛⃗𝑑𝑆 = 0
𝑉 𝑆

Donc pour un fluide incompressible, l’interprétation physique de cette équation est la suivante :
les débits entrant et sortant à travers une surface quelconque fermée doivent être égaux.

Cette équation est d’une grande utilité dans toute l’hydrodynamique.

VI. Ecoulement potentiel et écoulement plan

1. Ecoulement potentiel

⃗ dérive d’un potentiel Ф tel que 𝑉

Un écoulement est potentiel lorsque le champ de vitesse 𝑉 ⃗ =
⃗ = ⃗0.
−𝑔𝑟𝑎𝑑Ф ⇔ 𝑟𝑜𝑡𝑉
Il existe une grande classe d’écoulements où la vitesse dérive d’un potentiel Ф. C’est le cas des
écoulements en milieux poreux où le potentiel hydraulique est la charge hydraulique (au
𝑉2
coefficient K près) qui se résume à deux termes seulement car le terme d’énergie cinétique 2𝑔

est négligeable (voir chapitre suivant).

𝑃
Ф = 𝐾𝐻 = 𝐾(𝑧 + )
𝜌𝑔
La présence du signe moins dans l’expression du vecteur vitesse, oriente cette dernière dans le
sens des potentiels décroissants.
⃗ = ⃗0 est dit irrotationnel.
Tout écoulement pour lequel 𝑟𝑜𝑡𝑉
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ Ф) = 0 soit ∆Ф = 0.
On en déduit que 𝑑𝑖𝑣(𝑔𝑟𝑎𝑑
Ф est une fonction harmonique car elle vérifie l’équation de Laplace.

2. Ecoulement plan

⃗ = (𝑉𝑥 , 𝑉𝑦 , 0).
Un écoulement est plan si le champ de vecteur vitesse s’écrit 𝑉

32
 ⃗ = 𝑟𝑜𝑡𝜎 et
Pour un écoulement plan d’un fluide incompressible, on peut trouver 𝜎 tel que 𝑉
0
𝜎 = (0)
𝜓

La fonction 𝜓 ainsi définie est appelée fonction de courant.

3. Ecoulement potentiel plan

Pour un écoulement potentiel et plan, ψ est aussi une fonction harmonique ∆ψ = 0.

Ceci parce qu’on aura les équations de Cauchy-Riemann:
𝜕ψ 𝜕Ф
𝑉𝑥 = =−
𝜕𝑦 𝜕𝑥
𝜕ψ 𝜕Ф
𝑉𝑦 = − =−
𝜕𝑥 𝜕𝑦
Il s’ensuit que le Laplacien de la fonction courant est nul :
𝜕 2Ф 𝜕 2Ф
∆𝜓 = − =0
𝜕𝑥𝜕𝑦 𝜕𝑦𝜕𝑥

4. Lignes de courant et équipotentielles

Les lignes de courant sont les lieux des points où 𝜓 = 𝑐𝑡𝑒.

Les équipotentielles sont les lieux des points où Ф = 𝑐𝑡𝑒.
Lorsqu’une ligne équipotentielle Ф et une ligne de courant 𝜓 se coupent en un point donné 𝑃,
elles sont orthogonales. Ces lignes forment donc un réseau orthogonal.

33
5. Ecoulements dans les milieux poreux

Equation de Darcy
Un milieu poreux est un milieu composé d’une phase solide et d’une (ou plusieurs) phase
liquide.
L’écoulement à travers un milieu poreux isotrope est donné par l’équation de Darcy :
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝐻
⃗ = −𝐾𝑔𝑟𝑎𝑑
𝑉
⃗ est le vecteur vitesse moyenne sur la section, appelée vitesse superficielle ou de filtration.
Où 𝑉
𝐾 est la conductivité hydraulique ou perméabilité.
Dans la direction selon 𝑥 on écrit :
𝜕𝐻
𝑉𝑥 = −𝐾
𝜕𝑥
𝑃
Ici 𝐻 se réduit à la charge piézométrique 𝑧 + 𝜌𝑔.

Ingénieur général des ponts et chaussées,

Henri Darcy (1803 - 1858) est un
hydraulicien français. Ses travaux en
adduction d’eau et sur le chemin de fer de la
ville de Dijon ont contribué fortement au
développement de sa ville.
Dans son célèbre ouvrage Les fontaines
publiques de la ville de Dijon, il établit la loi
empirique de Darcy permettant d’exprimer le
débit d’un fluide incompressible en régime
stationnaire en écoulement au travers d’un
milieu poreux.
La loi de Darcy et ses dérivées sont

Henry Darcy (1803 – 1858) aujourd’hui constamment utilisées dans des

domaines à enjeux forts pour la sécurité des
travaux publics (transport, construction...),
ainsi que pour l'alimentation en eau, en gaz et
en pétrole et de nombreux domaines de la
géotechnique et de l'industrie.

34
On rencontre des écoulements poreux dans différentes circonstances i) les écoulements
souterrains – mouvements des eaux, du pétrole ou du gaz dans les sols ou les roches -, ii) les
écoulements dans les fondations et les digues en terre, iii) le flux d’humidité dans les
constructions, iv) la filtration, v) l’irrigation et le drainage.

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝐻), on obtient l’équation de continuité :

⃗ = 𝑑𝑖𝑣(-𝐾𝑔𝑟𝑎𝑑
En prenant 𝑑𝑖𝑣𝑉
𝜕𝑉𝑥 𝜕𝑉𝑦 𝜕𝑉𝑧
⃗ =
𝑑𝑖𝑣𝑉 + + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧
qui a pour conséquence
0 = −𝐾∆𝐻
Donc
𝜕²𝐻 𝜕²𝐻 𝜕²𝐻
∆𝐻 = + + =0
𝜕𝑥² 𝜕𝑦² 𝜕𝑧²
Les problèmes des écoulements dans les milieux poreux sont donc bien ramenés à ceux des
écoulements potentiels.

La cinématique dans les milieux poreux peut permettre de déterminer par exemple, le débit de
fuite sous la digue de ce barrage :

De déterminer de façon locale (autour d’un puits de pompage), la conductivité hydraulique de

la formation aquifère ci-dessous et par conséquent la productivité de l’ouvrage ou bien encore
le débit que l’on peut évacuer par pompage.

35
36
CHAPITRE IV : Dynamique des fluides incompressibles

La dynamique des fluides est l’étude de leur mouvement en tenant compte des forces qui
engendrent le mouvement.

I. Equations générales du mouvement

En statique, on a démontré que
1
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃 + 𝛤 = ⃗0
− 𝑔𝑟𝑎𝑑
𝜌
En dynamique on aura
1
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃 + 𝛤 = 𝛶
− 𝑔𝑟𝑎𝑑 ⃗
𝜌
⃗.
Il s’agit de la traduction de loi de la mécanique newtonienne selon laquelle ∑ 𝐹𝑒𝑥𝑡 = 𝑚𝛶

1. Equations de Navier Stokes

Lorsque le fluide est newtonien et incompressible, les équations de Newton appliquées à

chaque particule fluide prennent la forme des équations de Navier-Stokes.

Bilan des forces pour un fluide visqueux

Si l’on se place dans un référentiel galiléen et qu’on effectue un bilan de forces sur une particule
fluide située en M à l’instant 𝑡 de masse 𝑑𝑚 = 𝜌(𝑀, 𝑡)𝑑𝑣, il faut tenir compte :
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃𝑑𝑣
- des forces de pression −𝑔𝑟𝑎𝑑
- des forces de volume 𝜌𝛤 𝑑𝑣
- et des forces visqueuses 𝑑𝐹𝜇
La résultante des forces vaut donc :
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃)𝑑𝑣 + 𝑑𝐹𝜇
𝑑𝐹 = (𝜌𝛤 − 𝑔𝑟𝑎𝑑

L’expression 𝑑𝐹𝜇 est assez compliquée mais se simplifie dans le cas particulier des fluides
newtoniens et incompressibles.

37
Expression des forces visqueuses

Cas d’un écoulement parallèle unidimensionnel

Calculons la résultante des forces visqueuses dans le cas d’un écoulement suivant 𝑜𝑥 avec un
gradient de vitesse suivant 𝑜𝑦 :
𝑣 = 𝑣(𝑦)𝑖
On remarque qu’il s’agit bien d’un écoulement incompressible :
𝑑𝑣 𝑑𝑣 𝑑²𝑣
𝑑𝐹𝜇 = 𝜇 [ (𝑦 + 𝑑𝑦) − (𝑦)] 𝑑𝑥𝑑𝑧𝑖 = 𝜇 𝑑𝑣
𝑑𝑦 𝑑𝑦 𝑑𝑦²
Cette formule obtenue dans le cas particulier unidimensionnel, se généralise de la façon
suivante :
𝑑𝐹𝜇 = 𝜇∆𝑣𝑑𝑣

Equation de Navier-Stokes
D’après la seconde loi de Newton appliquée à une particule fluide :
𝑑𝑣
𝑑𝐹 = 𝜌𝑑𝑣 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃𝑑𝑣 + 𝜌𝛤 𝑑𝑣 + 𝜇∆𝑣𝑑𝑣
= −𝑔𝑟𝑎𝑑
𝑑𝑡
Pour un fluide newtonien incompressible, la dynamique de l’écoulement vérifie l’équation :
⃗
𝑑𝑉
𝜌 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃 + 𝜌𝛤 + 𝜇∆𝑣
= −𝑔𝑟𝑎𝑑
𝑑𝑡

38
 La « mystérieuse » équation de Navier Stokes a été trouvée quasi
simultanément par l’ingénieur, mathématicien et économiste français
Henri Navier (1785 –1836) et le mathématicien britannique Georges
Stokes (1819 –1903).
Puisqu’elle décrit le mouvement des fluides (gaz ou liquides), son
domaine d’application est vaste. Elle peut être utilisée pour étudier la
circulation du sang dans nos artères, la trajectoire de l’air autour
Henri Navier (1785
d’une aile d’avion, le mouvement des grandes masses d’air dans
–1836)
l’atmosphère, etc. Elle est un must pour les ingénieurs de toutes les
disciplines.
Cette équation est jugée compliquée car c’est une équation
différentielle non-linéaire à cause du terme 𝑣. ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑔𝑟𝑎𝑑 𝑣. Cette non-
linéarité explique certains comportements chaotique des fluides
appelée turbulence qu’il est difficile d’appréhender.
Aucune solution explicite de cette équation n’a été trouvée ; pour des

Georges Stokes conditions initiales définies, l’unicité d’une solution pour tout temps t

(1819 –1903) n’est pas démontrée.

L’équation de Navier Stokes fait partie des 7 problèmes du millénaire
mis à prix 1 000 000 de dollars par la fondation Clay.

2. Equations d’Euler : Cas des fluides non visqueux

L'équation d'Euler découle de l'équation de Navier-Stokes. Si l’on ne considère que l’action des
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑔𝑧 et si on néglige les forces de
forces de gravité comme forces de volume : 𝛤 = −𝑔𝑟𝑎𝑑
frottement visqueux :
𝑑𝑣
𝜌 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃 − 𝜌𝑔𝑟𝑎𝑑
= −𝑔𝑟𝑎𝑑 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑔𝑧
𝑑𝑡

Et pour un fluide incompressible

⃗
𝑑𝑉
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ (𝑃 + 𝜌𝑔𝑧) = 𝜌
−𝑔𝑟𝑎𝑑
𝑑𝑡

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑃∗ = 𝜌 𝑑𝑉⃗
Soit −𝑔𝑟𝑎𝑑 𝑑𝑡

39
3. Théorème de Bernoulli

Cas d’un fluide parfait

𝝏
En considérant que lors de l’écoulement, le mouvement est permanent (𝝏𝒕 ≡ 𝟎), il n’y a pas de

force de frottement (fluide parfait), la masse volumique est constante (fluide incompressible),
l’action de la gravité est la seule force extérieure;
Bernouilli a établi que le long d’une ligne de courant, la quantité définie par
𝑃 𝑉2
𝐻= +𝑧+
𝜌𝑔 2𝑔
est constante.
𝐻 appelée charge hydraulique et est définie à une constante additive près car l’origine de l’axe
vertical z est arbitraire. La dimension de la charge hydraulique est une longueur.

On peut démontrer le théorème de Bernoulli à partir de celui d’Euler et en considérant les

hypothèses simplificatrices citées ci-dessus.

La charge hydraulique représente une énergie par unité de poids qui est la somme de trois
termes :
𝑉2
: énergie cinétique par unité de poids
2𝑔

z: énergie potentielle de position par unité de poids

𝑃
: énergie des forces de pression par unité de poids
𝜌𝑔

Il s’agit donc d’un théorème de conservation de l’énergie.

Cas des fluides réels (visqueux)

Pour de tels fluides, il y a des forces de frottement entre les particules et les parois et entre
particules elles-mêmes, dues à la viscosité du fluide, qui vont dissiper l’énergie. La charge
hydraulique diminue donc le long de la ligne de courant dans le sens de l’écoulement.
Ainsi, sur une ligne de courant entre 2 points A et B, on a cette fois-ci
𝑃𝐴 𝑉𝐴 2 𝑃𝐵 𝑉𝐵 2
+ 𝑧𝐴 + = + 𝑧𝐵 + + ∆𝐻
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔
∆𝐻 est appelé perte de charge.

40
Pour résoudre un problème d’écoulement d’un fluide réel, on fait appel à des résultats
expérimentaux, en particulier ceux de l’ingénieur et physicien britannique Osborne Reynolds.

Daniel Bernoulli (1700 – 1782) est un médecin, physicien et

mathématicien suisse, ami d’Euler.
L’ouvrage Hydrodynamica où il publie le principe de Bernoulli est
considéré comme l’œuvre fondatrice de l’hydrodynamique.
Initialement utilisé pour des fluides en circulation dans une
conduite, ce principe a trouvé un important champ d'application en

Bernoulli (1700 – aérodynamique.

1782)

II. Interprétations physiques du théorème de Bernoulli

1. Interprétation en énergie

Par unité de masse Par unité de volume Par unité de poids

(𝐽/𝑘𝑔) ou 𝑚2 /𝑠 2 (𝐽/𝑚2 ) ou 𝑃𝑎 (𝐽/𝑁 𝑜𝑢 𝑚)
Energie potentielle 𝑔𝑧 𝜌𝑔𝑧 𝑧
de situation
Energie cinétique 𝑉2 𝑉2 𝑉2
𝜌
2 2 2𝑔
Energie potentielle 𝑃/𝜌 𝑃 𝑃
de pression 𝜌𝑔
Domaine Dynamique des Dynamique des Dynamique des
d’utilisation fluides fluides fluides
compressibles incompressibles incompressibles

𝑉2
L’expression relative à l’unité de volume est 𝑃 + 𝜌𝑔𝑧 + 𝜌 = 𝐶𝑡𝑒 est souvent utilisée en
2
aérodynamique ou en ventilation.
𝑃 𝑉2
Celle relative à l’unité de poids 𝜌𝑔 + 𝑧 + 2𝑔 = 𝐶𝑡𝑒 est souvent utilisée en hydraulique.

41
2. Interprétation en graphique

Dans l’expression relative à l’unité de poids, on a :

𝑧 = côte ou hauteur de position du fluide considéré.
𝑉2
ℎ𝑣 = hauteur cinétique ou hauteur d′ où doit tomber dans le vide une particule fluide
2𝑔
pour acquérir la vitesse V = √2gh.
𝑃
ℎ𝑃 = hauteur d′ une colonne du fluide considéré qui mesure la pression P.
𝜌𝑔

Plan de charge
V²/2g

P/ρg
Ligne piézométrique

Tube de courant
Plan de référence

𝑃
On appelle hauteur piézométrique la somme 𝜌𝑔 + 𝑧.
𝑃 𝑉2
Sa courbe représentative est appelée ligne piézométrique. La somme 𝑧 + 𝜌𝑔 + 2𝑔 est la

hauteur correspondant à la charge totale de l’écoulement.

III. Régimes d’écoulement – Nombre de Reynolds

Reynolds par ses expériences menées en 1883 lors de l’écoulement d’un liquide dans une
conduite cylindrique rectiligne dans laquelle arrive un filet de liquide coloré, a montré
l’existence de deux régimes d’écoulement :
- le régime d’écoulement laminaire
- le régime d’écoulement turbulent

42
Régime laminaire :
Les filets fluides sont des lignes régulières, sensiblement parallèles entre elles.

Filet coloré

Régime turbulent
Les filets s’enchevêtrent, s’enroulent sur eux-mêmes.

Filet coloré

En utilisant des fluides de diverses viscosités, et en faisant varier les conditions expérimentales
(débit et diamètre de la canalisation d’écoulement), Reynolds a montré qu’il existait un
paramètre 𝑅𝑒 , sans dimension, qui permet de classer ces différents comportements.
𝑉. 𝑑 𝜌. 𝑉. 𝑑
𝑅𝑒 = =
𝜈 𝜇
- 𝜌: masse volumique du fluide en 𝑘𝑔/𝑚3
- 𝑉 : Vitesse moyenne d’écoulement à travers la section considérée en 𝑚/𝑠
- 𝑑 : diamètre de la conduite ou largeur de la veine fluide en 𝑚
- 𝜈 : viscosité cinématique du fluide en 𝑚2 /𝑠
- 𝜇 : viscosité dynamique du fluide en 𝑃𝑎. 𝑠

Si 𝑅𝑒 < 2000, l’écoulement est laminaire

Si 𝑅𝑒 > 2000, l’écoulement est turbulent :
- lisse si 2000 < 𝑅𝑒 < 10000
- rugueux si 𝑅𝑒 > 10000

43
 Osborne Reynolds est un ingénieur et physicien irlandais qui fit
d’importantes contributions à l’hydrodynamique et à la dynamique
des fluides.
Léonard de Vinci fut le premier à étudier la formation de tourbillons
dans l’air ou dans l’eau qui s’écoule autour d’un obstacle :
phénomène auquel il donna le nom actuel de turbulence. Il n’en
chercha pas l’explication théorique mais en fit de très belles
descriptions.

Osborne Reynolds Il faut attendre 1883 pour que Reynolds établisse une règle pour

(1842 – 1912) définir si un écoulement sera laminaire ou turbulent. Le nombre

adimensionnel de Reynolds traduit le rapport du transfert par
convection au transfert par diffusion de la quantité de matière et
permet de distinguer les régimes d’écoulement.

Notion de perte de charge

Lors du déplacement d’un fluide réel, les frottements visqueux dissipent son énergie entraînant
des pertes de charges. Celles-ci correspondent à la pression supplémentaire qu’il faut exercer
aux extrémités d’une canalisation pour assurer un écoulement stationnaire et compenser le
frottement visqueux.
Deux termes entrent en jeu dans le calcul des pertes de charge :
- la perte de charge en ligne ou perte de charge régulière due aux frottements le long du trajet
- la perte de charge singulière due à la présence d’obstacles localisés tels que les coudes, les
robinets, les vannes, les modifications brutales de section etc.

L’expérience montre que si l’on considère un écoulement fluide dans une canalisation de
longueur 𝐿 et de diamètre 𝐷, la perte de charge ∆𝐻 vaut :
1 𝑉² 𝐿
∆𝐻 = 𝜆(𝑅𝑒)
2 𝑔𝐷
Cette relation s’appelle équation de Darcy-Weisbach.
Le facteur sans dimension 𝜆 désigne le coefficient de perte de charge régulière et ne dépend
que de Reynolds pour une canalisation lisse.

44
Pour une conduite rugueuse, 𝜆 dépend aussi de la rugosité relative є/𝐷 qui mesure le rapport
de la hauteur moyenne des aspérités de la paroi interne de la conduite sur son diamètre interne.
1 є 2,51
= −2𝑙𝑜𝑔 [ + ]
√𝜆 3,7𝐷 𝑅𝑒√𝜆
𝜆 peut s’obtenir grâce au diagramme de Moody page suivante.

En régime laminaire, la formule de Poiseuille donne également la perte de charge

8𝜇𝐿 8𝜇𝐿 𝐿
∆𝐻 = 𝑄𝑣 4
= 𝑉𝜋𝑅² 4
= 32𝜇𝑉
𝜌𝑔𝜋𝑅 𝜌𝑔𝜋𝑅 𝜌𝑔𝐷²

𝜌. 𝑉. 𝐷 𝜌. 𝑉. 𝐷 32 𝑉² 𝐿
𝑃𝑢𝑖𝑠𝑞𝑢𝑒 𝑅𝑒 = , 𝑎𝑙𝑜𝑟𝑠 𝜇 = 𝑒𝑡 ∆𝐻 =
𝜇 𝑅𝑒 𝑅𝑒 𝑔 𝐷
Finalement en régime laminaire :
1 𝑉² 𝐿 64
∆𝐻 = 𝜆 𝑎𝑣𝑒𝑐 𝜆 =
2 2𝑔 𝐷 𝑅𝑒

En régime turbulent, 𝑅𝑒 étant grand,

1 є
= −2𝑙𝑜𝑔 [ ]
√𝜆 3,7𝐷

la perte de charge ne dépend que de la rugosité relative.

45
46
CHAPITRE V : Notions d’hydraulique à surface libre

Les écoulements à surface libre sont ceux pour lesquels la surface libre est en contact avec
l’atmosphère. Pour ces écoulements, la surface libre est soumise à la pression atmosphérique.
On rencontre ce type d’écoulement dans les rivières naturelles ou aménagés, dans les canaux
d’irrigation, de drainage, dans les conduites dont le liquide n’occupe pas toute la section, dans
les réseaux d’assainissement urbain (eaux pluviales, eaux usées).
L’hydraulique à surface libre est la partie de l’hydraulique générale qui étudie les canaux à ciel
ouvert, les cours d’eau naturels aménagés, les canalisations ne fonctionnant pas à pleine section.
L’hydraulique générale utilise les principes fondamentaux de la mécanique comme la relation
fondamentale de la dynamique, le théorème de l’énergie cinétique, le théorème des quantités
de mouvement, le théorème de la conservation de l’énergie totale…
Les équations générales de l’hydraulique sont établies par la mécanique des fluides. Les
simplifications de ces équations obéissent à des hypothèses qui dépendent du type
d’écoulement.

I. Définitions

I. 2. Paramètres géométriques
Les paramètres géométriques essentiels sont :
- le tirant d’eau 𝑦 qui est la hauteur d’eau entre la surface libre du liquide et le fond du chenal ;
- la largeur en miroir 𝐿 ou largeur en gueule ou encore largeur du plan d’eau qui est largeur de
la surface mouillée au contact de la surface libre ;

47
 - la section mouillée 𝑆 est la surface de la section du chenal depuis la surface libre jusqu’au
fond du chenal ;
- le périmètre mouillé 𝑃 est la longueur de la ligne de contact, pour la section donnée, de l’eau
avec les parois et le fond du chenal ;
- le rayon 𝑅𝐻 hydraulique est le rapport entre section mouillée et périmètre mouillé.(𝑅𝐻 = 𝑆/𝑃)
- le diamètre hydraulique 𝐷𝐻 = 4𝑅𝐻

Pour un canal artificiel de forme trapézoïdale, on peut compléter ces définitions de paramètres
géométriques par les suivantes :

L
Rev
Rev Ligne d’eau
dz y
y 1
m Z1 Z2 fond
Θ
dx
b

- la largeur en radier ou base 𝑏

- les fruits des berges 𝑚 permettent de définir l’inclinaison du canal par rapport à la verticale
(ou l’horizontale)
𝑚 = 𝑐𝑜𝑡𝑔𝛩
- un des paramètres géométriques indispensable à l’étude des chenaux, artificiels comme
naturels, siège d’un écoulement à surface libre est la pente du fond du canal, généralement notée
𝑑𝑧
i. Elle est définie par 𝑖 = 𝑑𝑥, z étant la côte du fond du chenal à la section d’abscisse x.

Dans le cas d’une section trapézoïdale, on a les relations suivantes entre les paramètres
géométriques :
𝐿 = 𝑏 + 2𝑚𝑦
𝑆 = 𝑦(𝑏 + 𝑚𝑦)

48
𝑃 = 𝑏 + 2𝑦√1 + 𝑚²

Pour un canal rectangulaire 𝛼 = 90° ⇒ 𝑚 = 0

𝐿=𝑏 𝑆 = 𝑦𝑏 𝑃 = 𝑏 + 2𝑦

Pour un canal triangulaire 𝑏 = 0

𝐿 = 2𝑚𝑦 𝑆 = 𝑚𝑦² 𝑃 = 2𝑦√1 + 𝑚²

Dans le cas d’une section circulaire, on a les relations suivantes :

1
𝑆1 = 𝑅²𝛩
2
2𝜋 − 𝜃 𝜃
𝛼= =𝜋−
2 2
ℎ
𝑐𝑜𝑠𝛼 = ⇒ ℎ = 𝑅𝑐𝑜𝑠𝛼
𝑅
𝛩 𝛩
ℎ = 𝑅𝑐𝑜𝑠 (𝜋 − ) = −𝑅𝑐𝑜𝑠
2 2
𝐿 𝛩 𝜃
= 𝑅𝑠𝑖𝑛𝛼 = 𝑅𝑠𝑖𝑛 (𝜋 − ) ⇒ 𝐿 = 2𝑅𝑠𝑖𝑛
2 2 2
𝛩 𝛩
𝐿ℎ −2𝑅𝑠𝑖𝑛 2 𝑅𝑐𝑜𝑠 2 −𝑅²
𝑆2 = = = 𝑠𝑖𝑛𝛩
2 2 2
𝑅2
𝑆 = 𝑆1 + 𝑆2 = (𝛩 − 𝑠𝑖𝑛𝛩)
2
𝑃 = 𝑅𝛩

49
 𝛩
𝑦 = 𝑅 + ℎ = 𝑅 − 𝑅𝑐𝑜𝑠
2
𝛩
𝑦 = 𝑅(1 − 𝑐𝑜𝑠 )
2

I.3. Paramètres hydrauliques

Profil de vitesse

Berge à berge Surface au fond

Vitesses moyenne
1 1
𝑈 = 𝑆 ∑𝑖 𝑣𝑖 𝑑𝑆𝑖 𝑈 = 𝑆 ∬ 𝑣 𝑑𝑠

Débit
∆𝑣
𝑄= = 𝑆. 𝑈
∆𝑡

𝑄 = ∬ 𝑣 𝑑𝑠

I.4. Notion de charge hydraulique, théorème de Bernoulli

La charge hydraulique représentant l’énergie totale par unité de poids s’exprime en longueur,
l’unité est donc le mètre (𝑚). Géométriquement, la charge hydraulique est comptée
verticalement et a de ce fait pour expression :
𝑣2
𝐻 = 𝑧 + 𝑦. 𝑐𝑜𝑠(𝜂) +
2𝑔

50
 𝜂 = 𝐼 est la pente de fond du canal (en réalité 𝐼 = tan(𝜂). Dans le cadre de ce cours, cette pente
est supposée faible, ce qui permet d’écrire l’expression de la charge moyenne en un point
donné :
𝑣2
𝐻 =𝑧+𝑦+
2𝑔

𝑣2
2𝑔 𝐿𝑖𝑔𝑛𝑒 𝑑𝑒 𝑐ℎ𝑎𝑟𝑔𝑒 𝑜𝑢 𝑙𝑖𝑔𝑛𝑒 𝑑′é𝑛𝑒𝑟𝑔𝑖𝑒

𝑦. 𝑐𝑜𝑠(𝜂) 𝑦 𝑆𝑢𝑟𝑓𝑎𝑐𝑒 𝑙𝑖𝑏𝑟𝑒 𝑜𝑢 𝑙𝑖𝑔𝑛𝑒 𝑑′𝑒𝑎𝑢

𝑧 𝐹𝑜𝑛𝑑 𝑑𝑢 𝑐𝑎𝑛𝑎𝑙

𝜂
𝑅é𝑓é𝑟𝑒𝑛𝑐𝑒
Charge hydraulique en un point d'une section donnée

On peut définir une charge moyenne dans une section donnée. Dans ce cas, on introduit un
coefficient α, appelé coefficient de Coriolis ou encore coefficient d’énergie cinétique, qui
permet de tenir compte de la répartition des vitesses dans la section considérée. L’expression
de la charge moyenne est
𝑈2
𝐻 = 𝑧 + 𝑦 + 𝛼 2𝑔

Le coefficient de Coriolis est défini comme suit :

1
𝛼= ∫ 𝑣 3 𝑑𝑠
𝑈3𝑆 𝑆
En écoulement laminaire 𝛼 vaut 2.
Il vaut 1,02 en écoulement turbulent.
En HSL, le régime d’écoulement est presque toujours turbulent, ainsi on devrait prendre 𝛼 =
1,02 dans l’expression de la charge moyenne. Mais vu que des erreurs de lecture et de mesure
peuvent être commises facilement, on prendra une fois pour toute 𝛼 = 1.

D’une manière générale, la charge entre totale varie entre deux sections distinctes à cause des
frottements sur les parois (perte d’énergie). La différence entre les deux charges est appelée
51
perte de charge. La pente de la ligne d’énergie n’est rien d’autre que la perte de charge par unité
de longueur ou encore perte de charge unitaire. Elle est notée 𝐽 et a pour expression
𝑑𝐻
𝐽=−
𝑑𝑥
Cette définition traduit en même temps le théorème de Bernoulli, à savoir que la variation de
charge entre deux sections distinctes est égale à la perte de charge entre ces deux sections :
𝐻1 − 𝐻2 = 𝐽∆𝑥

II. Ecoulement uniforme

II.1. Définition
Un écoulement permanent est dit uniforme lorsque les valeurs de la hauteur d’eau, de la vitesse
moyenne et du débit ne changent pas d’une section à l’autre. Pour un tel écoulement, les filets
liquides sont rectilignes, parallèles entre eux et parallèles aux parois, les trajectoires ne se
distinguent pas des lignes de courant.

II.2. Etude des écoulements uniformes

II.2.1. Principe fondamental

De la définition de l’écoulement uniforme, on arrive à une définition synthétique qui résume
les définitions et hypothèse d’un écoulement uniforme : Dans un écoulement uniforme à surface
libre, la perte de charge unitaire est égale à la pente du fond du canal elle-même parallèle à la
surface libre : 𝑰 = 𝑱 = 𝑰𝟎. Tous les calculs, dans le cas d’un tel type d’écoulement, se feront à
partir de ce principe fondamental.
𝐽 représente la pente de la ligne de charge, 𝐼0 celle de la surface libre et 𝐼 celle du canal.

II.2.2. Calcul de la pente d’énergie

Dans le principe fondamental énoncé ci-dessus, on voit que la connaissance de la perte de
charge suffit à connaitre la pente du fond ou celle de la surface libre. La relation entre le
frottement contre la paroi s’exprime en fonction de la pente de la ligne d’énergie par 𝜏0 =
𝜌𝑔𝑅𝐻 𝐽, ou 𝑅𝐻 est le rayon hydraulique dans la section considérée. En outre 𝑅𝐻 𝐽 est fonction
de la vitesse moyenne U de l’écoulement, mesurée dans cette section. 𝑅𝐻 𝐽 = 𝜑(𝑈).

52
Cette fonction peut être expliquée par un coefficient de rugosité, appelé souvent par abus de
langage « rugosité ».

II.2.3. Equations de l’écoulement uniforme

La relation 𝑅𝐻 𝐽 = 𝜑(𝑈) devient 𝑅𝐻 𝐼 = 𝜑(𝑈) dans le cas des écoulements uniformes en vertu
du principe fondamental énoncé ci-dessus. La connaissance 𝜑(𝑈) permet d’obtenir une relation
qui permet de résoudre les problèmes qui peuvent se poser dans un écoulement uniforme (voir
plus loin).
𝑈2
Selon Chézy 𝜑(𝑈) peut s’écrire sous la forme 𝜑(𝑈) = ou C est appelé coefficient de Chézy
𝐶2

qui tient compte de la rugosité du canal, de sa forme, mais aussi des conditions d’écoulement.
On a ainsi, ce qui s’appelle la formule de Chézy qui s’écrit :
𝑈 = 𝐶√𝑅𝐻 𝐼

Manning et Strickler ont reformulé la relation de Chézy sous la forme :

2/3
𝑈 = 𝐾𝑆 𝑅𝐻 √𝐼, pour Strickler, 𝐾𝑆 étant le coefficient de rugosité de Strickler. Il est fonction de
la nature de la paroi ; il est d’autant plus grand que celle-ci est lisse.
1 2/3
𝑈 = 𝑛 𝑅𝐻 √𝐼, pour Manning n étant le coefficient de rugosité de Manning
1
On voit que ces deux formules sont identiques, avec la relation 𝐾𝑆 = 𝑛. On peut également

déduire une relation entre l’un de ces coefficients et le coefficient de Chézy :

1/6 1 1/6
𝐶 = 𝐾𝑆 𝑅𝐻 ou encore 𝐶 = 𝑛 𝑅𝐻 .

Ces différentes formules permettent aussi l’estimation du débit. L’estimation du débit selon la
formule de Strickler pourra s’écrire.

2/3 𝑆 5/3
𝑄= 𝐾𝑆 𝑆𝑅𝐻 √𝐼 = 𝐾𝑆 2/3 √𝐼
𝑃

Selon Darcy-Weisback, la perte de charge unitaire s’écrit :

𝜆 𝑈² 𝜆 𝑄²
𝐽= =
𝐷 2𝑔 𝐷 2𝑔𝑆²

Selon Calmon-Lechapt

53
 𝑄𝑛
𝐽=𝑎
𝐷𝑚

II.2.4. Estimation des valeurs de 𝐊 𝐒

Le coefficient de Strickler est le paramètre le plus difficile à connaitre et la plus grande source
d’erreur. Les méthodes d’estimation du coefficient de rugosité de Strickler 𝐾𝑆 , et donc de
Manning 𝑛, ainsi que des tables permettant de choisir des valeurs de ces coefficients en
fonctions des types de parois, pour des matériaux cohérents, sont présentées dans Mar 2004.
Dans les cas des matériaux non cohérents, on peut estimer ce coefficient à l’aide de la formule
𝐾𝑆 = 26(𝑅𝐻 /𝑑35 )1/6
Où d35 est le diamètre des grains tel que 35% en poids des particules ont un diamètre inférieur
à d35.
Pour les sols cohérents, le coefficient de rugosité de Strickler est fourni par des tables dédiées.

Estimation des valeurs de 𝐊 𝐒 dans le cas des sections composées

Considérons une section mouillée composée de n autres sections ayant des rugosités propres.

A’ H’
A H
1
3
2

B C F
G

D E

On considère que la vitesse est identique dans chaque section :

𝑈𝑖 = 𝑈𝑗 = 𝑈
2/3
2/3 𝑆1
𝑈𝑖 = 𝐾𝑆1 𝑅𝐻 √𝐼 = 𝐾𝑆1 2/3
√𝐼
𝑃1
On déduit que
3/2
𝑈
𝑆𝑖 = 𝑃𝑖 ( )
𝐾𝑆1 √𝐼

54
On écrit enfin que la section totale est la somme des sections élémentaires :

𝑆 = ∑ 𝑆𝑖
𝑖

Soit :
3/2 3/2
𝑈 𝑈
𝑃( ) = ∑ 𝑃𝑖 ( )
𝐾é𝑞 √𝐼 𝑖
𝐾𝑆1
√𝐼
Donc :
𝑃 𝑃𝑖
=∑
𝐾𝑒𝑞 3/2 𝑖
𝐾𝑆𝑖 3/2

II.3. Problèmes concernant les écoulements uniformes

II.3.1. Problème de type 1

C’est le cas le plus facile à étudier. Lorsqu’on connait la géométrie de notre canal et qu’on
désire calculer un des paramètres, connaissant les deux autres (comme présenté dans le tableau
ci-dessous), il suffit d’appliquer directement la formule de Manning-Strickler.

On connait On cherche
I 𝑲𝑺 Q
Q 𝑲𝑺 I
Q I 𝑲𝑺

II.3.2. Problème de type 2

Ce type de problème se pose lorsque pour un débit connu, et pour une hauteur normale d’eau
donnée, on cherche les dimensions du canal en mesure de contenir ce débit OU lorsque la forme
du canal est connue et qu’on cherche la position de la ligne d’eau (on dit qu’on cherche la
hauteur d’eau normale ou profondeur normale).

Notion de profondeur normale

La profondeur normale 𝑦𝑛 pour un débit Q donné et une pente I donnée est la profondeur de
l’écoulement en régime uniforme. Dans ce cas, la surface libre représente le niveau normal. Le

55
rayon hydraulique, la section mouillée et le périmètre mouillé, correspondant à la profondeur
normal, sont appelés les éléments normaux.

Il existe 3 méthodes de résolution de ce problème itératif :

- Résolution par méthode numérique
- Résolution par abaques
- Résolution par la méthode de la débitance

II.3.2.1. Résolution par les abaques

Les canaux les plus courants ont souvent des formes trapézoïdales (de préférence isocèle s’il
n’y a pas de contrainte particulière) ou circulaires, à la limite parabolique.
Pour des canaux ayant des sections trapézoïdales ou circulaires, des abaques permettent de lire
directement la profondeur normale à partir des caractéristiques de la section.
- Pour les sections circulaires, il suffit de calculer la quantité :
𝑄
8
𝐾𝑆 √𝐼𝐷 3

et de la projeter sur la courbe (graduation logarithmique) correspondant à la section circulaire

𝑦𝑛
pour obtenir la valeur de . Connaissant le diamètre du canal, on en déduit la profondeur
𝐷

normale
- Pour les sections trapézoïdales, on calcule la quantité
𝑄
8
𝐾𝑆 √𝐼𝑏 3

On projette le résultat sur la courbe correspondant au fruit de berge de la section trapézoïdale

𝑦𝑛
considérée pour obtenir la valeur de on en déduit la profondeur normale.
𝑏

56
 𝑄
8
𝐾𝑆 √𝐼𝐷 3

𝑦𝑛 D
𝐷 2 𝑦𝑛 𝑦𝑛
 y 𝑏

1 𝑦𝑛
m
b

𝑄 57
8
𝐾𝑆 √𝐼𝑏 3
III.3.2.2. Résolution par la méthode de la débitance
La méthode de la débitance, part de l’égalité suivante, issue de la formule de Manning-
Strickler :
𝑄 2/3
= 𝐾𝑆 𝑆𝑅𝐻
√𝐼
𝑄
On peut remarquer que la partie gauche de l’égalité est fixe, puisque le débit et la pente sont
√𝐼
2/3
entièrement connus. La deuxième partie, 𝐷 = 𝐾𝑆 𝑆𝑅𝐻 , appelée débitance, en revanche est
variable.
Le but de la méthode est de construire la courbe de la débitance en fonction du tirant d’eau 𝑦.
𝑄
La profondeur normale est obtenue en projetant la débitance particulière sur la courbe
√𝐼

obtenue.
𝑦

𝑦𝑛

2/3
𝑄 𝐷 = 𝐾𝑆 𝑆𝑅𝐻
√𝐼

On peut aussi procéder par incréments de 𝑦 de telle sorte à aller dans le sens de la convergence
𝑄
de la débitance vers le rapport . Cette méthode peut être illustrée par le tableau suivant :
√𝐼
𝑄
y S P R D √𝐼
1 21 22.8284271 0.91990569 1390.42064
2 44.00 25.66 1.71 4412.83
3 69.00 28.49 2.42 8711.63
2.5 56.25 27.07 2.08 6411.58
2.25 50.06 26.36 1.90 5373.76
2.251 50.09 26.37 1.90 5377.76
2.2505 50.07 26.37 1.90 5375.76 5375.76

58
On peut écrire à l’ordinateur ou sur une calculatrice programmable de petits programmes
capables de résoudre des équations implicites telles 𝑦 = 𝑓(𝑦).

II.3.3. Problème de type 3

Ce problème se pose lorsque l’on dispose d’un débit fixé et que l’on désire dimensionner le
canal (largeur et hauteur d’eau correspondant à ce débit) de manière à évacuer ce débit, pour
une pente 𝐼 et une rugosité données.

Section hydrauliquement favorable de forme circulaire

De tous les canaux, les sections circulaire et demi-circulaire, de par leur géométrie, sont celles
qui offrent pour une section mouillé donnée, le plus petit périmètre. Elles sont donc
naturellement hydrauliquement favorables. Mais les difficultés de mise en œuvre ne militent
pas en leur faveur.

Section hydrauliquement favorable de forme trapézoïdale

Pour déterminer les relations existant entre les dimensions (largeur et tirant d’eau) permettant
de réaliser un tel type de problème, on écrira que le périmètre mouillée est minimum, sachant
que la section mouillée est fixée. Cela revient à écrire que les dérivées partielles de ces deux
grandeurs sont nulles.
Pour une section trapézoïdale, on a 𝑆 = 𝑦(𝑏 + 𝑚𝑦) et 𝑃 = 𝑏 + 2𝑦√1 + 𝑚2 , on a ainsi

𝑑𝑆 = 0 ⇒ 𝑦𝑑𝑏 + (𝑏 + 2𝑚𝑦)𝑑𝑦 = 0
{
𝑑𝑃 = 0 ⇒ 𝑑𝑏 + 2√1 + 𝑚2 𝑑𝑦 = 0

Ce système n’a de solution différente de la solution triviale (0,0) que si le discriminant du

système est nul. Ainsi, on a

2𝑦√1 + 𝑚2 − (𝑏 + 2𝑚𝑦) = 0
Soit

𝑏 = 2𝑦 (√1 + 𝑚2 − 𝑚)

A partir de cette relation entre la largeur au radier et la hauteur, on peut déduire celles de la
surface mouillée, du périmètre mouillé et du rayon hydraulique :

59
𝑆 = 𝑦(𝑏 + 𝑚𝑦) = 𝑦 (2𝑦 (√1 + 𝑚2 − 𝑚) + 𝑚𝑦) = 𝑦 2 (2√1 + 𝑚2 − 𝑚)

Posons 𝜆 = 2√1 + 𝑚2 − 𝑚
𝑆 s’écrit tout simplement
𝑆 = 𝜆𝑦 2
Par ailleurs,

𝑃 = 𝑏 + 2𝑦√1 + 𝑚2 = 2𝑦 (√1 + 𝑚2 − 𝑚) + 2𝑦√1 + 𝑚2 = 2𝑦 (2√1 + 𝑚2 − 𝑚)

Soit finalement :
𝑃 = 2𝜆𝑦
Le rayon hydraulique s’en déduit immédiatement :
𝑦
𝑅𝐻 =
2
Remarque : En adoptant la même démarche pour une section rectangulaire et pour une section
𝑦
demi-circulaire, on arrive au même résultat, à savoir que 𝑅𝐻 = 2.

60
CHAPITRE VI : Notions d’hydraulique en charge

I. Etude des pertes de charge

Il existe deux types de pertes de charge : les pertes de charge linéaires et les pertes de charge
singulières.

Les premières sont celles produites par la conduite sur le profil en long, les secondes sont
localisées. On appelle les premières pertes de charges régulières ou encore linéaires vu que la
ligne de charge est approximativement une droite.

Pertes de charges linéaires

L étant une longueur imposée, l’estimation de la perte de charge revient à estimer la pente de
la ligne de charge.

Plusieurs formules permettent d’estimer cette pente :

- Maning-Strickler
2/3
𝑄 = 𝐾𝑆 𝑆𝑅𝐻 √𝐽

- Chézy

𝑄 = 𝐶𝑆√𝑅𝐻 𝐽

- Darcy-Weisback

𝜆𝐿 𝑈 2 𝜆𝐿 𝑄 2
𝐽𝐿 = =
𝐷 2𝑔 𝐷 2𝑔𝑆 2
8𝜆𝐿 𝑄²
𝐽=
𝜋²𝑔 𝐷5

- Calmon-Lechapt

𝑄𝑛
𝐽=𝑎
𝐷𝑚

- William-Hazen

61
 10,675 𝑄1,852
𝐽= 1,852 𝐷 4,87
𝐶𝐻𝑊

On remarquera que la formule de Calmon-Lechapt est une généralisation de toutes les autres
formules de pertes de charge. Il suffit de déterminer les coefficients 𝑎, 𝑛 et 𝑚.

Pour Darcy-Weisback

8𝜆
𝑎= , 𝑛 = 2, 𝑚 = 5
𝜋²𝑔

Pour William-Hazen

10,675
𝑎= 1,852 , 𝑛 = 1,852, 𝑚 = 4,87
𝐶𝐻𝑊

II. Problèmes courants en HEC

Dans la résolution des problèmes courants susceptibles d’être rencontrés en HEC on fera
constamment recours au théorème de Bernoulli.

Le principe fondamental de la résolution de ces problèmes est simple : autant d’équations

que d’inconnues.

Les problèmes que nous serons amenés à résoudre sont : calcul du débit 𝑄, calcul du diamètre
𝐷 (dimensionnement de la conduite), calcul de la côte d’eau dans un réservoir, une bâche, un
lac, un barrage ou tout autre plan d’eau, calcul de la puissance d’une pompe ou d’une turbine,
association série/parallèle des conduites, calcul de butée.

62
Application

Une conduite de diamètre 𝐷, de longueur 𝐿, de rugosité 𝑘 sert à effectuer la vidange d’un

réservoir vers un autre. Sur la conduite sont localisés un débitmètre et ainsi qu’une vanne. Les
coefficients de pertes de charge de toutes les singularités (vanne, débitmètre, rétrécissement
brusque à la sortie du réservoir amont, élargissement brusque à l’entrée du réservoir aval) sont
répertoriés sur la figure. Les niveaux d’eau des réservoirs étant maintenus constants quel est
le débit véhiculé par la conduite ?

L’estimation de pertes de charge linéaires s’effectuera par la formule de Darcy-Weisback.

𝑧1 = 1260 𝑚

𝑘𝐷 = 1,2
𝑘𝑅𝐵 = 0,8
D
𝑧2 = 1245 𝑚
V

𝑘𝑉 = 0,52
𝑘𝐸𝐵 = 1
𝐷 = 600 𝑚𝑚
Conduite 𝐿 = 6350 𝑚
𝜀 = 0,1 𝑚𝑚

63