COURS DE PHYSIQUE

MÉCANIQUE DES FLUIDES

JIMMY ROUSSEL

2025

[Link]/mecanique_des_fluides
Cours de mécanique des fluides – [Link]
JIMMY ROUSSEL, professeur agrégé à l’Ecole Nationale Supérieure de Chimie de Rennes

Copyright © 2025 Jimmy Roussel

cbn Ce document est sous licence Creative Commons « Attribution - Pas d’Utilisation Commerciale
4.0 International (CC BY-NC 4.0) ».
Pour plus d’informations : [Link]/licenses/by-nc/4.0/
Ce document est réalisé avec l’aide de KOMA - Script et LATEX en utilisant la classe kaobook.
1re édition – Janv. 2010
Version en ligne – [Link]/mecanique_des_fluides
Preface

Ce cours de mécanique des fluides est avant tout destiné à l’étudiant désirant maîtriser les bases de
la mécanique des fluides. Il est normalement accessible à un étudiant en fin de parcours de Licence
(L2-L3).
Ce cours couvre les aspects essentiels de la mécanique des fluides : notion de pression, tension super-
ficielle, écoulements parfaits, écoulements visqueux, notion de pertes de charges etc. En revanche la
notion de turbulence n’est pas abordée.
J’ai essayé le plus possible d’illustrer les différentes notions par des exemples ou de simples exercices.
Mais pour un entraînement plus poussé, j’invite le lecteur à se procurer les eBooks suivants :
• Statique des fluides - Pression et tension de surface – 45 exercices et problèmes corrigés.
• Dynamique des fluides - Fluides parfaits et newtoniens – 50 exercices et problèmes corrigés.

disponibles à l’adresse [Link]/femto

Jimmy Roussel
Table des matières
Preface iii

Table des matières v

1 PRESSION DANS UN FLUIDE 1

1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 La pression . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Fluide au repos dans le champ de pesanteur . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4 Généralisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

2 TENSION DE SURFACE 13
2.1 Énergie de surface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2 Forces capillaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.3 Mouillage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.4 Approche thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21

3 FLUIDE EN ÉCOULEMENT 25
3.1 Description d’un fluide en écoulement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.2 Conservation de la masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.3 Écoulement incompressible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31

4 DYNAMIQUE DES FLUIDES PARFAITS 35

4.1 Équations dynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2 Écoulements incompressibles stationnaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.3 Écoulements potentiels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41

5 FLUIDES VISQUEUX 45
5.1 Notion de viscosité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
5.2 Dynamique des écoulements visqueux incompressibles . . . . . . . . . . . . . . . . . . . . 47

6 PERTES DE CHARGE ET FRICTION 51

6.1 Pertes de charge dans une conduite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
6.2 Résistance au mouvement dans les fluides . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59

ANNEXES 63
A THÉORÈME DE LA DIVERGENCE 65

B ACCÉLÉRATION D’UNE PARTICULE DE FLUIDE 67

C THÉORÈME π 69
C.1 Énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
C.2 Applications en mécanique des fluides . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69

Pour en savoir plus 73

Grandeurs physiques et symboles mathématiques 75

Table des figures

1.1 Les trois échelles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1

1.2 Pression cinétique. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Le rapport d𝐹 /d𝑆 fluctue autour d’une valeur moyenne appelée pression. Les fluctuations
sont imperceptibles à l’échelle mésoscopique. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.4 Particule de fluide soumise aux forces pressantes du fluide l’entourant dans le cas où la
pression ne dépend que de 𝑥. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.5 Résultante des forces de pression sur une surface fermée. . . . . . . . . . . . . . . . . . . . . . 5
1.6 À gauche : manomètre à liquide. À droite : expérience de Torricelli (1644). . . . . . . . . . . . . 7
1.7 Nivellement barométrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.8 Évolution de la surface libre en cas d’accélération horizontale. . . . . . . . . . . . . . . . . . . . 9
1.9 Cylindre en rotation. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.10 Calcul de la poussée d’Archimède sur un cube. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.11 Cas d’un volume quelconque. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.1 Interactions au sein d’un fluide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

2.3 Production de particules dans la première chambre à bulles . . . . . . . . . . . . . . . . . . . . 14

2.2 Coalescence de deux gouttes (©CNRS) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.4 Mise en évidence des forces capillaires. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.5 Densité linéique de forces superficielles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.6 Stalagmométrie. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.7 Méthode de l’anneau. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.8 La résultante des forces capillaires qu’exerce la membrane sur un contour fermé est nulle si
la membrane est plane. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.9 À gauche : existence d’une surpression dans une goutte. La courbure entraîne l’existence
d’une résultante des forces capillaires dirigée vers le centre de courbure. À droite : forces
agissantes sur une portion hémisphérique d’une goutte liquide. Les flèches réparties sur le
contour circulaire représentent les forces capillaires. . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.10 Lame de savon formant une caténoïde . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.11 Définition de l’angle de contact. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.12 Différents équilibres d’un liquide au contact d’un substrat solide. . . . . . . . . . . . . . . . . . 20
2.13 Expérience et schéma. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.14 Bulle de savon enfermant de l’air à la pression 𝑝. . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

3.1 Notion de trajectoire. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

3.2 Notion de ligne d’écoulement. Notez que toutes les vitesses représentées concernent le même
instant. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.3 Visualisation des tourbillons derrière un obstacle plan. . . . . . . . . . . . . . . . . . . . . . . . 28
3.4 À gauche : lignes d’émission autour d’un corps fuselé à l’aide de la machine à fumées d’Etienne
Marey. À droite : ondes de choc lors d’un vol supersonique visualisées par strioscopie. . . . . 29
3.5 Calcul du débit massique. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.6 Calcul du débit massique sortant d’une surface fermée. . . . . . . . . . . . . . . . . . . . . . . . 30
3.7 Particule de fluide dans un écoulement uniforme . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
3.8 Particule de fluide dans un écoulement radial, capturé à différents instants. . . . . . . . . . . 32
3.9 Tube de courant. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4.1 Daniel Bernoulli (1700 - 1782) : médecin, physicien et mathématicien suisse. Il publia en 1738
son ouvrage Hydrodynamica dans lequel il expose le théorème fondamental de la mécanique
des fluides qui porte son nom. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.2 Effet venturi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.3 Vidange d’un récipient. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.4 Sonde de Pitot statique dite aussi sonde de Prandtl. . . . . . . . . . . . . . . . . . . . . . . . . . 40
4.5 Principe de la sonde de Pitot statique. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.6 Exemple d’écoulement potentiel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43

5.1 Rhéomètre de Couette. À gauche, vue de profil. À droite, vue de dessus. . . . . . . . . . . . . . 45

5.2 Contraintes internes dans le fluide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
5.3 Différents comportements rhéologiques. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
5.4 Bilan des forces de viscosité sur un élément de fluide. . . . . . . . . . . . . . . . . . . . . . . . . 48

6.1 Écoulement de Poiseuille. Position du problème. . . . . . . . . . . . . . . . . . . . . . . . . . . . 52

6.2 Profil de vitesse. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
6.3 Diagramme de Moody. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
6.4 Coefficients de perte de charge singulière sur quelques géométries. . . . . . . . . . . . . . . . 57
6.5 Ecoulement stationnaire dans un système de conduites. . . . . . . . . . . . . . . . . . . . . . . 57
6.6 Traînée sur un obstacle sphérique immobile. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
6.7 Portance et traînée. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
6.8 Quelques 𝐶𝑥 pour différents profils dans un écoulement turbulent avec 𝑅e ≃ 105 . . . . . . . 61
6.9 Évolution du 𝐶𝑥 d’une sphère en fonction du nombre de Reynolds. . . . . . . . . . . . . . . . . 61

A.1 Calcul du flux à travers un cube. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65

A.2 Flux à travers deux cubes adjacents . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66

C.1 Traînée sur un obstacle sphérique immobile. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69

C.2 Convection dans un échangeur cylindrique. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71

Liste des tableaux

1.1 Autres unités de pression couramment utilisées. . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

1.2 Eau liquide à 25°C, 1 atm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3 Vapeur d’eau à 100°C, 1 atm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

5.1 Quelques valeurs de viscosité. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46

5.2 Quelques ordres de grandeur du nombre de Reynolds. . . . . . . . . . . . . . . . . . . . . . . . 50
 PRESSION DANS UN FLUIDE 1
Dans ce cours, on s’intéresse au champ de pression au sein d’un 1.1 Introduction . . . . . . . . 1
fluide, notamment lorsqu’il est au repos dans un champ de force ex- 1.2 La pression . . . . . . . . 2
térieure. Origine de la pression . . 2
Unités . . . . . . . . . . . 3
Résultante des forces
Version en ligne
pressantes . . . . . . . . 3
[Link]/mecanique_des_fluides/[Link] Conséquences . . . . . . 5
1.3 Fluide au repos dans le
champ de pesanteur . . . 5
Cas des Liquides . . . . . 5
1.1 Introduction Cas des gaz . . . . . . . . 7
1.4 Généralisation . . . . . . 8
Équation fondamentale de
Le terme fluide désigne un comportement qui s’oppose au comporte-
la statique . . . . . . . . . 8
ment élastique ou plastique associé aux solides. Par définition, on dit
Liquide en translation
que la matière est fluide lorsqu’elle se déforme aussi longtemps que accélérée . . . . . . . . . 9
lui sont appliquées des contraintes tangentielles. En termes simples Liquide en rotation . . . . 10
on peut dire qu’un fluide coule quand un solide se déforme. Poussée d’Archimède . . . 10

Fondamentalement, le comportement fluide est lié, au niveau molé-

culaire, à l’absence d’ordre à longue portée (contrairement aux cris- 1 : Certains systèmes présentent un
taux1 ) et à l’existence d’un chaos moléculaire (contrairement aux so- ordre à longue portée suivant une
lides). Ces propriétés se retrouvent notamment chez les gaz et les seule direction. Ils ont alors un carac-
tère cristallin selon cette direction et
liquides. fluide selon les autres. On les désigne
par le terme cristaux liquides.
Plusieurs approches sont possibles pour décrire un fluide :

• L’approche « Dynamique Moléculaire » consiste à simuler le mou-

vement moléculaire en résolvant de façon numérique les équa-
tions mécaniques appliquées à un grand nombre de molécules.
• L’approche « milieu continu » est une description à une échelle
où les grandeurs thermodynamiques ont un sens.
• L’approche statistique traite la dynamique moléculaire en adop-
tant un point de vue statistique.

𝐿≫𝑎≫ℓ

• • • • •

• • • • • •
• •
• •
• • ••
M • ••
• • • •
•
•
• •
𝑣(M,
⃗⃗⃗⃗ 𝑡) •
• • • •
• • • •
• •
• • • ••
• •
•
FIG. 1.1 : Les trois échelles.

La mécanique des fluides repose sur la deuxième approche. En effet,

dans les situations courantes on peut, en général, distinguer trois
échelles :

1. L’échelle macroscopique 𝐿. Par exemple, le diamètre d’un tuyau

d’arrosage caractérise l’échelle macroscopique d’un tel système.
 2 1 PRESSION DANS UN FLUIDE

2. L’échelle des collisions ℓ ≪ 𝐿. Généralement on choisit le libre

parcours moyen, c’est-à-dire la distance moyenne parcourue
par une molécule entre deux collisions successives. À cette échelle,
les grandeurs varient de façon discontinue et imprévisible.
3. L’échelle mésoscopique 𝑎 telle que ℓ ≪ 𝑎 ≪ 𝐿. À cette échelle,
les fluctuations sont lissées de sorte que l’on peut définir des
grandeurs locales continues.
Particule de fluide – On choisit alors comme échelle d’observation,
l’échelle mésoscopique. On considère, autour d’un point M, un vo-
lume mésoscopique 𝛿𝜏 . Typiquement un volume de 1 µm3 convient.
Ce volume contient un grand nombre de particules ce qui permet de
définir des grandeurs moyennes locales qui, elles, vont évoluer de
façon continue : la masse volumique locale 𝜌(M, 𝑡), la vitesse locale
𝑣(M,
⃗ 𝑡)...

1.2 La pression

Origine de la pression

Le chaos moléculaire qui règne au sein d’un fluide est à l’origine de la

pression cinétique. Les collisions sur un élément de surface mésosco-
fluide paroi
pique ou macroscopique produisent une force moyenne qui présente
les propriétés suivantes :
d𝑆
⃗⃗⃗⃗⃗⃗
d𝐹 • cette force est orthogonale à l’élément de surface ;
⃗⃗⃗⃗⃗⃗⃗⃗
𝑣1 • son intensité est proportionnelle à l’aire de l’élément de sur-
⃗⃗⃗⃗⃗⃗⃗⃗
𝑣2 face.
Si l’on note d𝐹 l’intensité de la force et d𝑆 l’aire de l’élément de
FIG. 1.2 : Pression cinétique.
surface, on a
d𝐹 = 𝑝c d𝑆
où 𝑝c désigne la pression cinétique2 .
L’interaction intermoléculaire est aussi responsable d’une pression
2 : On montre en théorie cinétique des
gaz que moléculaire 𝑝m . Les interactions entre molécules peuvent amplifier
1 ⋆
(𝑝m > 0) ou diminuer (𝑝m < 0) cette pression cinétique. Les collisions
𝑝c = 𝑛 𝑚⟨𝑣2 ⟩ inter-moléculaires l’augmentent alors que les interactions attractives
3
où 𝑛⋆ est la densité moléculaire (mo-
de van der Waals la diminuent.
lécules par m3 ), 𝑚 la masse d’une mo-
En résumé, dans un fluide tout élément de surface subit des forces
lécule et ⟨𝑣2 ⟩ le carré moyen de la vi-
tesse des particules. pressantes qui s’écrivent

⃗⃗⃗ ⃗⃗ = 𝑝 d𝑆 ⃗⃗⃗⃗
d𝐹 𝑛 avec 𝑝 = 𝑝c + 𝑝m ♡ (1.1)

d𝐹
d𝑆
microscopique mésoscopique macroscopique
FIG. 1.3 : Le rapport d𝐹 /d𝑆 fluctue au-
tour d’une valeur moyenne appelée 1 nm 1 µm 1 mm
pression. Les fluctuations sont imper-
ceptibles à l’échelle mésoscopique. Échelle
 1.2 La pression 3

𝑝 désigne la pression qui règne au point où se trouve l’élément de

surface. La pression est donc un champ scalaire que nous noterons
𝑝(M) avec M un « point mésoscopique » du fluide.
Dans un gaz, la température est di-
• Dans les gaz, la pression d’origine cinétique est prépondérante rectement reliée à l’énergie cinétique
(𝑝m ≪ 𝑝c ). moyenne des molécules, via la rela-
• Dans le cas limite où l’on néglige l’interaction moléculaire, on tion
1 3
obtient un gaz parfait pour lequel 𝑚⟨𝑣2 ⟩ = 𝑘B 𝑇
2 2
où 𝑘B désigne la constante de Boltz-
𝑛𝑅𝑇 mann. Il en découle
𝑝GP = 𝑝c = ♡ (1.2)
𝑉 1 𝑛𝑁A 𝑛𝑅𝑇
𝑝c = 𝑚⟨𝑣2 ⟩ =
3 𝑉 𝑉
• Dans les liquides, les pressions moléculaire et cinétique sont à condition de poser 𝑅 = 𝑁A 𝑘B .
du même ordre de grandeur (en valeur absolue) de sorte qu’en
général 𝑝 = 𝑝c + 𝑝m ≪ 𝑝c avec 𝑝m < 0. On peut même rencon-
trer des états liquides à pression négative, mais ces états sont
métastables.

Unités
3 : Blaise Pascal (1623–1662) : Mathé-
Dans le Système international d’unités, la pression s’exprime en pas- maticien, physicien et philosophe fran-
çais. Il contribua au développement
cal (symbole Pa), en hommage à Blaise Pascal3 . des probabilités, de la philosophie
(« Les pensées ») et de l’hydraulique. Il
1 Pa = 1 N.m−2 ♡ (1.3) montra, avec la collaboration de son
beau frère Perrier, que l’ascension du
mercure dans l’expérience de Torricelli
Il existe également d’autres unités encore très utilisées (TAB. 1.1). est due à la pression atmosphérique.
TAB. 1.1 : Autres unités de pression cou-
Unités Équivalence en pascal
ramment utilisées et quelques ordres
1 bar 105 Pa de grandeur.
1 atmosphère 1,01325 ⋅ 105 Pa
1 torr = 1 mm Hg 133,3 Pa
Pressions rencontrées dans la Nature
Centre du Soleil 1016 Pa
Centre de la Terre 1011 Pa
Record Haute pression en labo 1, 5 ⋅ 1010 Pa
Abysse des océans 108 Pa
Talons aiguilles 107 Pa
Atmosphère de la terre 105 Pa
Trompe à eau 103 Pa
Ultra vide en labo 10−12 Pa

Exercice – À quelle masse répartie sur 1 cm2 correspond 1 bar si l’on prend
𝑔 = 10 m.s−2 ?
Solution : 1 kg.

Résultante des forces pressantes

Comme on vient de le voir, les forces pressantes sont des actions de

contact (interactions de courte portée) qui agissent sur des surfaces.
Si l’on considère une particule de fluide, son environnement produit
à sa surface des forces pressantes dont on cherche à exprimer la
résultante.
 4 1 PRESSION DANS UN FLUIDE

particule de fluide
d𝜏 = d𝑥 × d𝑦 × d𝑧
𝑧

𝑦
𝑝(𝑥) d𝑦d𝑧 𝑝(𝑥 + d𝑥) d𝑦d𝑧 𝑥
FIG. 1.4 : Particule de fluide soumise
aux forces pressantes du fluide l’en-
tourant dans le cas où la pression ne
dépend que de 𝑥.

Considérons pour cela une particule de fluide en forme de parallélé-

pipède située en (𝑥, 𝑦, 𝑧) de côtés d𝑥, d𝑦 et d𝑧. Supposons dans un
premier temps que la pression varie seulement avec 𝑥. On voit alors
par symétrie, que la résultante des forces est nécessairement selon
⃗⃗⃗ ⃗⃗ cette résultante. On a
𝑢𝑥 . Appelons d𝐹
⃗⃗⃗⃗⃗⃗⃗⃗

⃗⃗⃗ ⃗⃗ = 𝑝(𝑥) d𝑦d𝑧 ⃗⃗⃗⃗⃗⃗⃗⃗

d𝐹 𝑢𝑥 − 𝑝(𝑥 + d𝑥) d𝑦d𝑧 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥

d𝑥 caractérisant l’échelle mésoscopique, on peut utiliser le dévelop-

pement de Taylor d’ordre un :

𝜕𝑝
𝑝(𝑥 + d𝑥) = 𝑝(𝑥) + d𝑥
𝜕𝑥
On aboutit alors au résultat

⃗⃗⃗ ⃗⃗ = − 𝜕𝑝
d𝐹 d𝑥d𝑦d𝑧 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥
𝜕𝑥
Si l’on avait considéré une pression dépendant de la variable 𝑦 nous
aurions obtenu une force selon ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑦 faisant intervenir la dérivée par-
tielle 𝜕𝑝/𝜕𝑦, etc. Ainsi, dans le cas où la pression dépend des trois
variables d’espace 𝑥, 𝑦, 𝑧, on obtient

⃗⃗⃗ ⃗⃗ = − ( 𝜕𝑝 ⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 𝑢𝑥 +
𝜕𝑝
⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑦 +
𝜕𝑝
𝑢 ) d𝑥d𝑦d𝑧
⃗⃗⃗⃗⃗⃗⃗⃗
𝜕𝑥
⏟⏟⏟⏟⏟⏟⏟⏟⏟⏟⏟ 𝜕𝑦 𝜕𝑧 𝑧 ⏟
d𝜏
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑝

Résultante des forces pressantes

Une particule de fluide de volume d𝜏 subit de la part du fluide qui

l’entoure, une force de pression totale

⃗⃗⃗ ⃗⃗ = −⃗⃗⃗⃗⃗⃗
d𝐹 ∇𝑝(M) d𝜏 (1.4)

4 : Toutes les forces volumiques se- On peut donc définir une force volumique4 de pression (en N ⋅
ront notées avec un 𝑓 minuscule. m−3 ) :
⃗⃗⃗ ⃗⃗
d𝐹
⃗⃗⃗⃗⃗⃗⃗
𝑓𝑝 = = −⃗⃗⃗⃗⃗⃗
∇𝑝(M) (1.5)
d𝜏

En sommant toutes ces forces sur un volume macroscopique V déli-

mité par une surface S, on obtient la résultante des forces pressantes
qui agit sur sa surface. On a donc deux moyens d’exprimer cette ré-
 1.3 Fluide au repos dans le champ de pesanteur 5

sultante :

𝑛ext
⃗⃗⃗⃗⃗
⃗⃗⃗ ⃗⃗ = ∯
𝐹 𝑛ext = ∭
−𝑝(M) d𝑆 ⃗⃗⃗⃗ −⃗⃗⃗⃗⃗⃗
∇𝑝(M) d𝜏 ♡ (1.6)
M∈𝑆 M∈𝑉

Volume V

Conséquences
Fluide extérieur

• Comme on l’a vu, on peut définir une force volumique associée FIG. 1.5 : Résultante des forces de pres-
sion sur une surface fermée.
à la résultante des forces pressantes :

⃗⃗⃗ ⃗⃗
d𝐹
⃗⃗⃗⃗⃗⃗⃗
𝑓𝑝 = = −⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad(𝑝)
d𝜏
On constate donc que cette force volumique est une force conser-
vative qui dérive d’une énergie potentielle volumique

𝑒𝑝 (J ⋅ m−3 ) = 𝑝 (Pa) ♡ (1.7)

La mécanique nous enseigne alors que dans un fluide les forces

de pression tendent à déplacer le fluide vers les zones de basse
pression pour minimiser l’énergie potentielle de pression.
• Par ailleurs, l’une des propriétés du gradient d’une fonction
𝑓(𝑥, 𝑦, 𝑧) est de donner un vecteur perpendiculaire à la sur-
face de niveau définie par 𝑓(𝑥, 𝑦, 𝑧) = Cte . Il en découle que les
lignes de force de pression sont orthogonales aux isobares.
• Enfin, si la pression est uniforme on a ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad(𝑝) = 0⃗ et donc

𝑛ext = 0⃗
∯ 𝑝d𝑆 ⃗⃗⃗⃗ (pression uniforme) ♡ (1.8)

Exercice – Soit un solide hémisphérique de rayon 𝑅 plongé dans l’air de

⃗⃗⃗ ⃗⃗ qui s’exerce sur la sur-
pression uniforme 𝑝0 . Calculer la force pressante 𝐹
face sphérique du solide en utilisant la propriété (1.8).

1.3 Fluide au repos dans le champ de

pesanteur

Dans un fluide au repos par rapport à un référentiel, les forces in-

ternes se résument aux forces de pression. Voyons dans un premier
temps comment le champ de pesanteur conditionne l’évolution de la
pression au sein du fluide.

Liquide dans un champ de pesanteur

Les liquides – Dans un liquide, les interactions intermoléculaires5 5 : L’interaction de van der Waals, la
jouent un rôle important. Elles sont telles que les molécules sont liaison hydrogène, l’interaction élec-
trostatique dans une solution électro-
quasi en contact ce qui explique le caractère quasi-incompressible lytique etc.
 6 1 PRESSION DANS UN FLUIDE

de ces fluides : les liquides présentent un volume propre. En effet,

les variations de densité avec la pression et la température restent
faibles comme le montre la TAB. 1.2. On adopte alors le modèle sui-
vant :

Modèle de l’état liquide

|Δ𝜌|
Δ𝑇 = 10 K → = 0, 2% Sauf indication contraire, on considère les liquides incompressibles
𝜌
et indilatables. Autrement dit
|Δ𝜌|
Δ𝑝 = 1 bar → = 0, 02%
𝜌
𝜌 = Cte
TAB. 1.2 : Eau liquide à 25°C, 1 atm.

Pression dans un liquide incompressible – Considérons un liquide au

repos dans le champ de pesanteur 𝑔.⃗ Isolons une particule de fluide
de volume d𝜏 et effectuons un bilan des forces :

⃗⃗⃗ ⃗⃗ = 𝜌d𝜏 𝑔 ⃗ ;
• Pesanteur : d𝑃
⃗⃗⃗ ⃗⃗ = −⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
• Pression du fluide environnant : d𝐹 grad𝑝 d𝜏

Appliquons la seconde loi de Newton, sachant que le liquide est sup-

posé au repos dans le référentiel terrestre considéré galiléen.

𝜌d𝜏 𝑔 ⃗ − ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑝 d𝜏 = 0⃗

ce qui donne en projection sur les axes cartésiens :

Attention, l’axe (𝑂𝑧) est descendant ! 𝜕𝑝 𝜕𝑝 𝜕𝑝

=0 =0 et = 𝜌𝑔
𝜕𝑥 𝜕𝑦 𝜕𝑧

La pression ne dépend que de 𝑧 d’après les deux premières relations.

En intégrant la dernière relation on aboutit à

𝑝(𝑧) = 𝑝0 + 𝜌 𝑔 𝑧 ♡ (1.9)

La pression augmente linéairement avec la profondeur. Le terme 𝜌𝑔𝑧

représente la pression due au poids de la colonne de liquide qui
surplombe un point à la profondeur z.

Exercice – De quelle profondeur faut-il s’enfoncer dans l’océan pour voir la

pression augmenter de 1 bar si l’on prend 𝑔 = 9, 8 m.s−2 ?
Solution : 10,2 m.

Conséquences –

• Les isobares sont horizontales, par conséquent la surface libre

(isobare 𝑝 = 𝑝atm ) aussi.
6 : Élève de Galilée, il découvre en 1644 • Mesure relative de pression à l’aide d’un manomètre à liquide
le principe du baromètre à mercure. en U (cf. FIG. 1.6) :
Le 19 septembre 1648, le beau-frère de 𝑝int − 𝑝atm = 𝜌𝑔ℎ
Blaise Pascal, Florin Périer, montre que
la hauteur du mercure dans le baro- • Mesure absolue de pression à l’aide du baromètre à mercure
mètre de Torricelli ne dépend que de inventé par Evangelista Torricelli6 (cf. FIG. 1.6). On a
la pression atmosphérique. Il obtient à
Clermont Ferrand (460m) ℎ = 71,2 cm
et sur le Puy de Dôme (1464m) ℎ = 𝑝atm = 𝜌𝑔ℎ ⇒ 1 atm = 1, 01325.105 Pa → ℎ ≃ 760 mm de Hg
62,7 cm de mercure.
 1.3 Fluide au repos dans le champ de pesanteur 7

Vide
Atmosphère à la
pression 𝑝atm 𝑔⃗⃗⃗⃗
Hgliq
Gaz sous ℎ
pression Liquide
𝑝int > 𝑝atm ℎ 𝑝atm
FIG. 1.6 : À gauche : manomètre à li-
quide. À droite : expérience de Torricel-
li (1644).

Avant 1954, 1 atm correspondait à une hauteur de mercure de 76 cm pris

à 0 °C en un lieu de champ de pesanteur standard 𝑔 = 9, 806650 m.s−2
ce qui donnait 𝑝atm = 1, 013249.105 Pa. À l’heure actuelle on a 1 atm =
1,013 25 × 105 Pa par définition.

Gaz dans un champ de pesanteur

Les gaz – Dans un gaz, les molécules interagissent peu, et leur énergie |Δ𝜌|
Δ𝑇 = 10 K → = 3%
est avant tout cinétique. Les distances inter-atomiques sont grandes 𝜌
ce qui explique que, contrairement aux liquides, les gaz sont très com- |Δ𝜌|
Δ𝑝 = 1 bar → = 100%
pressibles. Par ailleurs, dans les conditions de pression et de tempé- 𝜌
rature raisonnables et loin de tout point critique, le modèle du gaz TAB. 1.3 : Vapeur d’eau à 100°C, 1 atm.
parfait est tout à fait suffisant.

Modèle de l’état gazeux

Sauf indication contraire, on considère les gaz parfaits. Par consé-

quent, si l’on suppose l’équilibre thermodynamique à l’échelle lo-
cale, on aura en chaque point du fluide

𝑚 𝑛𝑀gaz 𝑀gaz 𝑝(M)

𝜌(M) = = =
𝑉 𝑉 𝑅𝑇 (M)
−1
où 𝑀gaz est la masse molaire du gaz et 𝑅 = 8,314 J ⋅ K−1 ⋅ mol la
constante des gaz parfaits.

Pression dans une atmosphère isotherme – Considérons un gaz iso-

therme (𝑇 = 𝑇0 ) en équilibre dans le champ de pesanteur terrestre
𝑢𝑧 et procédons comme dans la section précédente. L’équi-
𝑔 ⃗ = −𝑔 ⃗⃗⃗⃗⃗⃗⃗⃗ Attention ici on prend un axe (O𝑧) as-
libre mécanique se traduit par la même équation : cendant !

𝜌d𝜏 𝑔 ⃗ − ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑝 d𝜏 = 0⃗

Là encore, la pression ne dépend que de 𝑧, et la projection selon l’axe

vertical donne l’équation

d𝑝
= −𝜌𝑔
d𝑧
Ici 𝜌 dépend de 𝑧 car la pression en dépend. En utilisant l’équation
d’état du gaz parfait, on obtient l’équation différentielle

d𝑝 𝑀gaz 𝑔
+ 𝑝=0 (1.10)
d𝑧 𝑅𝑇0
 8 1 PRESSION DANS UN FLUIDE

𝑝
Cette équation différentielle d’ordre un sans second membre a pour
𝑝0 solution :
𝑅𝑇0
𝑝(𝑧) = 𝑝0 e−𝑧/𝐻 avec 𝐻 = (1.11)
𝑀gaz 𝑔
où 𝑝0 est la pression en 𝑧 = 0. La pression, comme la masse volu-
37%𝑝0
mique décroît exponentiellement avec l’altitude 𝑧.

5%𝑝0 𝑧 • Appliquée à l’atmosphère isotherme elle donne la loi de nivelle-

𝐻 3𝐻 ment barométrique. Si l’on prend une température moyenne de
FIG. 1.7 : Nivellement barométrique : 15°C, la distance caractéristique 𝐻 vaut alors (pour l’air 𝑀air =
décroissance exponentielle de la pres- −1
sion avec l’altitude.
29 g ⋅ mol )

8,314 × 288
𝐻= = 8,4 km
29 × 10−3 × 9,81

• Pour une petite variation d’altitude, la loi (1.10) donne

𝛿𝑝 𝛿𝑧 𝛿𝑝
≃− ⇒ ∼ −10−4 pour 𝛿𝑧 = 1 m
𝑝 𝐻 𝑝

En conclusion, on peut négliger la pesanteur dans un gaz sur

une échelle de quelques mètres.

Exercice – Considérons une colonne cylindrique verticale d’atmosphère de

section droite 𝑆 = 1 m2 . Calculer sa masse à partir de la loi (1.11) sachant
que 𝑝0 = 1013,25 hPa et 𝑔 = 9,81 m ⋅ s−2 .
Solution : 10,3 × 103 kg.

1.4 Généralisation

Équation fondamentale de la statique

La démarche que nous avons suivi pour étudier l’équilibre d’un fluide
dans un champ de pesanteur, peut se généraliser à un fluide dans un
champ de forces extérieures quelconque. L’équation que l’on obtient
est l’équation fondamentale de la statique des fluides.

Supposons un fluide au repos dans un référentiel galiléen et soumis à

un champ de forces extérieures. Chaque parcelle de fluide de volume
d𝜏 subit, en plus des forces de pressions, une force
ext
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 = 𝑓⃗⃗⃗⃗ext d𝜏

de sorte que l’équilibre se traduit par l’équation

𝑓⃗⃗⃗⃗ext d𝜏 − ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑝 d𝜏 = 0⃗

d’où l’on tire

𝑓⃗⃗⃗⃗ext − ⃗⃗⃗⃗⃗⃗
∇𝑝 = 0⃗ ♡ (1.12)

Cette équation associée à l’équation d’état du fluide permet de dé-

terminer l’évolution de la pression. La résolution de l’équation diffé-
 1.4 Généralisation 9

rentielle fait apparaître une constante d’intégration que l’on obtient

à l’aide des conditions aux limites.

Liquide en translation accélérée

𝑧
Imaginons un verre d’eau initialement au repos dans un train immo-
𝑥
bile. La surface libre de l’eau est alors horizontale. À partir de 𝑡 = 0,
le train adopte un mouvement uniformément accéléré, d’accélération 𝑔⃗⃗⃗⃗
𝑢𝑥 où ⃗⃗⃗⃗⃗⃗⃗⃗
𝑎⃗ = 𝑎 ⃗⃗⃗⃗⃗⃗⃗⃗ 𝑢𝑥 est horizontal. Après un court instant pendant lequel
le liquide oscille, un état de repos est atteint. Cherchons comment la
pression varie dans le verre d’eau.
𝑎 ⃗⃗⃗⃗
⃗⃗⃗⃗ = 0
Tout d’abord, raisonnons dans le référentiel lié au verre : ce dernier
n’est pas galiléen puisqu’il est accéléré par rapport au référentiel ter- 𝑧
restre. Dans ce référentiel, il faut alors tenir compte des forces d’iner- 𝑥
tie. Faisons le bilan des forces extérieures que subit un élément de
𝑔⃗⃗⃗⃗
fluide de volume d𝜏 et de masse volumique 𝜌.

⃗⃗⃗ ⃗⃗
d𝑃
⃗⃗⃗ ⃗⃗ = 𝜌d𝜏 𝑔 ⃗
• Pesanteur : d𝑃 d’où ⃗⃗⃗⃗⃗⃗⃗
𝑓𝑔 = = 𝜌 𝑔.⃗ 𝑎
⃗⃗⃗⃗
d𝜏
d⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹ie FIG. 1.8 : Évolution de la surface libre
𝐹ie = −𝜌d𝜏 𝑎⃗ d’où ⃗⃗⃗⃗⃗⃗⃗⃗⃗
• Force d’inertie : d⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑓ie = = −𝜌𝑎.⃗ en cas d’accélération horizontale.
d𝜏

Appliquons l’équation fondamentale de la statique (1.12) :

𝜌 𝑔 ⃗ − 𝜌𝑎⃗ − ⃗⃗⃗⃗⃗⃗
∇𝑝 = 0⃗

Puisque le gradient de pression a deux composantes, l’une selon la

verticale (orientée selon ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 ) et l’autre selon l’horizontale (selon ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥 ),
on en déduit que la pression dépend des coordonnées 𝑥 et 𝑧. Proje-
tons tout d’abord l’équation de la statique selon ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 :

𝜕𝑝
= −𝜌𝑔 puis intégrons 𝑝(𝑥, 𝑧) = −𝜌𝑔𝑧 + 𝑓(𝑥)
𝜕𝑧
Dérivons la relation obtenue par rapport à 𝑥 et comparons à ce que
donne la projection de l’équation de la statique selon ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥 :

𝜕𝑝 𝜕𝑝
= 𝑓 ′ (𝑥) et = −𝜌 𝑎
𝜕𝑥 𝜕𝑥
ce qui mène à 𝑓 ′ (𝑥) = −𝜌 𝑎, soit 𝑓(𝑥) = −𝜌𝑎𝑥 + 𝑓(0). Finalement le
champ de pression s’écrit

𝑝(𝑥, 𝑧) = −𝜌𝑔𝑧 − 𝜌𝑎𝑥 + 𝑝(0, 0)

Les surfaces isobares sont donc des plans d’équation −𝜌𝑔𝑧 − 𝜌𝑎𝑥 =
Cte . La surface libre étant une surface isobare particulière elle vérifie
la même équation. On constate alors que la surface libre est inclinée
par rapport à l’horizontal d’un angle 𝛼 donné par
𝑎
tan 𝛼 = −
𝑔
 10 1 PRESSION DANS UN FLUIDE

Liquide en rotation

𝜔 Un flacon cylindrique ouvert, contient un liquide de masse volumique

𝜌. On fait tourner le flacon autour de son axe à la vitesse angulaire 𝜔.
Ce liquide n’est donc pas au repos dans le référentiel du laboratoire.
Liquide
⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 Cependant, après un régime transitoire qui dépend de la viscosité du
𝑟 ⊗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝜃 liquide, celui-ci tourne de façon solide à la même vitesse angulaire
•
M ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑟 que le cylindre. Ainsi, dans le référentiel lié au cylindre, le liquide est
𝑔⃗⃗⃗⃗ au repos. On raisonnera donc dans ce référentiel tournant (noté ℛ′ )
munis d’un système de coordonnées cylindriques.

Bilan des forces – Le référentiel n’étant pas galiléen, il faut tenir

FIG. 1.9 : Cylindre en rotation.
compte des forces d’inertie :

• Tout d’abord, la force de Coriolis ⃗⃗⃗⃗⃗⃗⃗⃗⃗

𝑓𝑖𝑐 = −2𝜌⃗⃗𝜔⃗ ⃗ ∧ 𝑣(M/ℛ
⃗ ′
) est nulle
Télescope à miroir liquide car, 𝑣(M/ℛ
⃗ ′
) = 0.⃗
• La force d’entraînement (ou force centrifuge) vaut ⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑢𝑟 .
𝑓𝑖𝑒 = 𝜌𝑟𝜔2 ⃗⃗⃗⃗⃗⃗⃗⃗
Des équipes de l’Université Laval,
de l’University of British Columbia
⃗⃗⃗⃗⃗⃗⃗
• La force volumique de pesanteur vaut 𝑓𝑔 = 𝜌𝑔.⃗
et de l’Institut d’astrophysique de
Paris, ont mis au point un téles- L’équation de la statique des fluides (1.12) donne donc
cope Nommé Large Zénith Téles-
cope (LZT) dont le miroir primaire
𝑢𝑟 + 𝜌 𝑔 ⃗ − ⃗⃗⃗⃗⃗⃗
𝜌𝑟𝜔2 ⃗⃗⃗⃗⃗⃗⃗⃗ ∇𝑝 = 0⃗
fait 6 mètres de diamètre. Contrai-
rement aux télescopes conven-
tionnels dont le miroir est fait de ce qui donne
verre, le LZT a un miroir fait de
liquide réfléchissant, du mercure ⎧ 𝜕𝑝
plus précisément qui adopte une { = 𝜌𝑟𝜔2 → 𝑝 = 1/2𝜌𝜔2 𝑟2 + 𝑓(𝑧)
{ 𝜕𝑟
surface parabolique puisque mis { 𝜕𝑝
en rotation dans une cuve. Le mi- = 0 ↓
⎨ 𝑟𝜕𝜑
roir parabolique obtenu permet de {
faire l’image d’une étoile au foyer { 𝜕𝑝 𝜕𝑝
{ = −𝜌𝑔 et = 𝑓 ′ (𝑧) ⇒ 𝑓(𝑧) = −𝜌𝑔 𝑧 + Cte
de la parabole avec précision. Sa ⎩ 𝜕𝑧 𝜕𝑧
focale vaut

𝑓=
𝑔 Finalement la pression dépend de la distance à l’axe 𝑟 et de 𝑧 :
2𝜔2
La principale limitation du LZT, 𝜌𝑟2 𝜔2
et des autres miroirs liquides, 𝑝(𝑟, 𝑧) = − 𝜌𝑔 𝑧 + Cte
2
est qu’on ne peut le pointer
ailleurs qu’au zénith, sinon le li- La surface libre étant une surface isobare, elle obéit à l’équation
quide s’échappe de la cuvette.
Au-dessus du miroir, une étroite
𝜔2 2
bande de ciel défile à la vitesse de 𝑧 = 𝑧0 + 𝑟
la rotation terrestre. Un astre don- 2𝑔
né met environ une minute à tra-
verser le champ de vision du téles- La surface libre adopte une forme parabolique d’axe de révolution
cope. (O𝑧).

Poussée d’Archimède

Imaginons un solide cubique d’arête 𝑎 complètement immergé dans

un liquide au repos dans le champ de pesanteur. Calculons la résul-
tante des forces de pression ⃗⃗Π⃗ ⃗⃗ qui s’exercent sur le cube.

• Tout d’abord, la pression ne dépend que de la profondeur :

𝑝(𝑧) = 𝑝0 + 𝜌 𝑔 𝑧
 1.4 Généralisation 11

• Par symétrie, les forces de pression horizontales se compensent, Position d’équilibre

contrairement aux forces verticales du fait de la variation de
pression avec la profondeur. ⃗⃗Π⃗ ⃗⃗ est donc suivant O𝑧 : ⃗⃗⃗⃗⃗⃗haut
𝐹

⃗⃗Π⃗ ⃗⃗ = 𝑝(𝑧) 𝑎2 𝑢
⃗⃗⃗⃗𝑧 − 𝑝(𝑧 + 𝑎) 𝑎2 𝑢
⃗⃗⃗⃗𝑧 ⃗⃗⃗⃗⃗⃗gauche
𝐹 ⃗⃗⃗⃗⃗⃗droite
𝐹
3
= −𝜌 𝑔 𝑎 𝑢 ⃗⃗⃗⃗𝑧
𝑔⃗⃗⃗⃗
⃗⃗⃗ ⃗⃗ = −𝑚d 𝑔 ⃗
Π ⃗⃗⃗⃗⃗⃗bas
𝐹

où 𝑚d désigne la masse de liquide déplacé. FIG. 1.10 : Calcul de la poussée d’Archi-

mède sur un cube.
On obtient donc une force ascendante opposée au poids du volume
de liquide déplacé.
Ce résultat particulier se généralise sans difficulté quelle que soit la
forme du corps immergé et quel que soit le fluide. Appelons V0 le
volume d’un solide et S0 l’aire de sa surface. Immergeons ce corps
dans un fluide quelconque au repos dans un champ de pesanteur.
La résultante des forces de pression qu’il subit s’écrit

⃗⃗⃗ ⃗⃗ = ∯
Π 𝑛ext d𝑆 = ∭
−𝑝(M) ⃗⃗⃗⃗ −⃗⃗⃗⃗⃗⃗
∇𝑝 d𝜏
M∈(S0 ) M∈(V0 )

La deuxième intégrale étant issue du résultat mathématique (1.6). Il

faut bien avoir conscience qu’il n y a pas, ici, de gradient de pression
dans le volume 𝑉0 puisqu’il n y pas pas de fluide. ⃗⃗⃗⃗⃗⃗
∇𝑝 représente le
gradient de pression qui règnerait si le volume 𝑉0 était occupé par
isobare
le fluide. Dans ce cas, on aurait 𝜌 𝑔 ⃗ = ⃗⃗⃗⃗⃗⃗
∇𝑝(M) puisque l’on suppose ⃗⃗⃗⃗⃗⃗
Π
l’équilibre. On en déduit 𝑛ext
⃗⃗⃗⃗⃗

⃗⃗Π⃗ ⃗⃗ = ∭ −𝜌(M)𝑔 ⃗ d𝜏 = −𝑚d 𝑔 ⃗

M∈(V0 )
𝑔⃗⃗⃗⃗
Volume 𝑉0
Cette force, opposée à la pesanteur, s’appelle la poussée d’Archi-
mède.
Fluide au repos
Théorème d’Archimède (250 av.J.C.) FIG. 1.11 : Cas d’un volume quelconque.

Tout corps immergé partiellement ou totalement dans un fluide

subit de la part de celui-ci une poussée verticale, dirigée vers le
haut, appelée poussée d’Archimède, dont l’intensité est égale au
poids du fluide déplacé. Le point d’application de cette force est le
centre de poussée ; il est différent, en général, du centre de gravité.

Applications –
• Flottaison des bateaux et autres corps flottants ;
• Ascension des ballons sondes ;
• Courants de convection etc.

Exercice – Je monte sur une balance qui m’indique que je « pèse 67,2 kg ».
Imaginons que la pression atmosphérique passe de 1 bar à 2 bar. Comment
évolue l’indication de la balance ? Pouvez-vous donner un ordre de gran-
deur ?
 TENSION DE SURFACE 2
De nombreuses expériences de la vie quotidienne ne peuvent pas 2.1 Énergie de surface . . . . 13
trouver d’explication avec les lois vues jusqu’ici. C’est en tenant compte Origine microscopique . . 13
des propriétés des interfaces qu’on peut les justifier notamment grâce Conséquences expérimen-
au concept de tension superficielle. La capillarité est la science qui tales . . . . . . . . . . . . 14
s’intéresse à ces phénomènes et qui joue un rôle majeur dans de Rôle des tensio-actifs . . 15
nombreux domaines scientifiques (climat, chimie de formulation, in- 2.2 Forces capillaires . . . . . 15
dustrie du verre etc.). Mise en évidence . . . . . 15
Application à la mesure de
On propose ici une présentation classique de la capillarité ; pour les tension de surface . . . . 16
aspects dynamiques et une vision plus moderne de cette science Théorème de Laplace . . 18
voir [1]. 2.3 Mouillage . . . . . . . . . 20
Angles de contact . . . . 20
Version en ligne Ascension capillaire . . . 21
[Link]/mecanique_des_fluides/[Link] 2.4 Approche thermodyna-
mique . . . . . . . . . . . 21
Identité thermodynamique 22
Exemple de calcul . . . . 22
2.1 Énergie de surface

Origine microscopique

Jusqu’ici nous avons considéré le fluide comme un milieu continu

contraint par des conditions aux limites que l’on a traité de façon
1 : Jusqu’ici on a admis la continui-
simpliste1 . En réalité, pour décrire correctement la physique aux in- té de la pression aux interfaces par
terfaces il faut tenir compte des interactions moléculaires à courte exemple.
portée qui s’exercent de part et d’autre de l’interface. La théorie clas-
sique de la capillarité consiste à modéliser une interface comme une
surface mathématique auquel on associe une certaine élasticité re-
présentée par la propriété physique que nous appellerons tension
superficielle.
• gaz ℱ2
Essayons de justifier cette notion par un modèle simple. Considérons
un liquide ℱ1 en contact avec sa vapeur ℱ2 . Au sein de chaque fluide
les molécules subissent des interactions de très courte portée (inter-
actions de van der Waals attractives). Ainsi, dans ℱ1 , chaque molé-
• interface
cule possède une énergie d’interaction 𝜖1 < 0 résultat de l’interac-
tion attractive avec son voisinage immédiat. De même on définira
une énergie d’interaction 𝜖2 au sein du fluide ℱ2 . Par contre, il existe
une couche de fluide dans laquelle les molécules sont soumises à
l’action des deux fluides. L’épaisseur de ce film moléculaire est de • liquide ℱ1
l’ordre de la dimension 𝑎 d’une molécule. Appelons alors 𝜖12 l’éner-
gie d’interaction d’une « molécule interfaciale ». Bien entendu, on a
|𝜖12 | compris entre |𝜖1 | et |𝜖2 |. Si 𝑁 est le nombre de molécules et 𝑁𝑠
FIG. 2.1 : Interactions au sein d’un
le nombre de molécules à l’interface, l’énergie du liquide ℱ1 vaut fluide.

ℰ1 = (𝑁 − 𝑁𝑠 )𝜖1 + 𝑁𝑠 𝜖12 = 𝑁 𝜖1 + ℰ𝑠

où 𝐸𝑠 = 𝑁𝑠 (𝜖12 − 𝜖1 ) représente l’énergie de l’interface. On voit donc

que l’on peut associer à l’interface une énergie liée à l’anisotropie
 14 2 TENSION DE SURFACE

des forces d’interaction moléculaire. Le terme d’anisotropie 𝜖12 − 𝜖1

est positif. Enfin, le nombre de molécules à l’interface est proportion-
nel à l’aire 𝑆 de l’interface. On a approximativement 𝑁𝑠 ≃ 𝑆/𝑎2 de
sorte que l’on peut écrire

𝜖12 − 𝜖1
ℰ𝑠 = 𝛾𝑆 avec 𝛾 ≃ >0 ♡ (2.1)
𝑎2

où 𝛾 est par définition la tension superficielle.

Ainsi, la tension superficielle est une grandeur positive qui caracté-

rise une interface et représente l’énergie interfaciale par unité de sur-
face. Elle s’exprime en J.m−2 et dépend notamment des interactions
moléculaires.

Conséquences expérimentales

Augmenter l’aire d’une interface de d𝑆, présente donc un coût éner-

gétique qui s’élève à
dℰ = 𝛾d𝑆 ♡ (2.2)

Forme des bulles et gouttes – Ainsi un liquide adoptera une forme

qui minimise la surface compte tenu des contraintes. On montre que
pour un volume donné la surface qui minimise l’énergie est une sphère.
Par exemple une goutte d’huile dans un mélange eau-alcool de den-
sité identique sera sphérique. De la même manière, les bulles de gaz
carbonique dans le champagne sont sphériques.

Coalescence – On montre aussi que deux gouttes sphériques auront

intérêt à former une goutte plus grosse. Ainsi quand on agite énergi-
quement un mélange eau-huile on obtient une émulsion de petites
gouttes d’huile dans l’eau. Cette émulsion est instable : les petites
gouttes coalescent et l’on obtient après un certain temps de l’huile
avec de l’eau au dessous.

FIG. 2.2 : Coalescence de deux gouttes

(©CNRS)

Exercice – Considérons deux gouttes d’eau sphériques de rayon 𝑟 qui coa-

lescent pour ne former qu’une seule goutte de rayon 𝑟′ . Montrer que cette
transformation s’accompagne d’un gain d’énergie.

Retard des changements d’état – Lorsque l’on détend un liquide de

façon isotherme, la thermodynamique prévoit qu’en dessous d’une
pression dite pression de vapeur saturante, le liquide change de phase
pour se vaporiser. Cependant la formation de la première bulle de
vapeur coute de l’énergie de sorte que le liquide peut exister dans
FIG. 2.3 : Production de particules dans une phase métastable en dessous de la pression de vapeur saturante.
la première chambre à bulles à hydro-
gène liquide du CERN On parle de retard à la vaporisation. Une simple perturbation locale
peut suffire à déclencher la formation d’une première bulle de vapeur.
C’est ce phénomène qui fut employé dans les détecteurs de parti-
2 : Ce phénomène fut utilisé dans cules du milieu du XXe siècle (chambres à bulles). De manière ana-
un autre détecteur de particule : la
logue, il y a retard à la liquéfaction pour la vapeur sursaturante2 .
chambre à Wilson.
 2.2 Forces capillaires 15

Rôle des tensio-actifs

Les tensioactifs sont constituées de molécules amphiphiles c’est-à-

dire munies d’un pôle hydrophile et d’une longue chaîne hydrophobe.
Lorsqu’un tensioactif est ajouté à de l’eau il vient se placer immédia-
tement à la surface, avec la queue hydrophobe pointant à l’extérieur
de la surface. Ce processus s’accompagne d’une diminution de l’éner-
gie de surface et donc d’une chute de la tension superficielle. Ce n’est
qu’une fois la surface saturée, et n’offrant plus d’espace disponible
à de nouvelles molécules amphiphiles, que les tensioactifs vont for-
mer des structures organisées au sein du liquide : ce sont les micelles.
Les micelles sont des structures sphériques ou ellipsoïdales dont la
surface est constituée des têtes hydrophiles des tensioactifs, alors
que les queues hydrophobes de ces derniers sont regroupées à l’in-
térieur. L’effet nettoyant des tensioactifs découle du fait que les sub-
stances hydrophobes, telles que les matières grasses, peuvent être
contenues à l’intérieur des micelles. La concentration de tensioac-
tifs au-dessus de laquelle les micelles commencent à se former est
connue comme la concentration micellaire critique (CMC). Elle s’ob-
tient en déterminant la concentration à partir de laquelle la tension
superficielle cesse de chuter.

2.2 Forces capillaires

Mise en évidence
Expérience

Plongeons un contour métallique dans de l’eau savonneuse puis reti-

rons le. Il se forme alors une membrane liquide plane qui s’appuie sur
le contour. Auparavant, nous avons pris soin de fixer un fil souple for-
mant une boucle qui adopte une forme quelconque. Lorsque l’on crève au
moyen d’une aiguille la membrane située dans la boucle, celle-ci adopte
aussitôt une géométrie circulaire (FIG. 2.4). Cette expérience montre que
la boucle est soumise à des forces dites forces capillaires aux propriétés
suivantes :

• Ces forces sont perpendiculaires en chaque point du contour et

tendent à minimiser l’aire du film d’eau savonneuse.
FIG. 2.4 : Mise en évidence des forces
• Ces forces sont tangentes à l’interface.
capillaires.
• Elles sont réparties de façon uniforme. On peut donc définir une
densité linéique de force d𝐹 /dℓ.

De la même façon, si l’on forme une lame d’eau savonneuse sur un Membrane liquide
cadre rectangulaire dont un des côtés est mobile, le liquide cherchant
à minimiser sa surface, il faut exercer une force sur la tige mobile pour ⃗⃗⃗⃗⃗⃗
𝐹 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹op
maintenir la surface constante.
d𝑥
On peut obtenir l’expression de cette force à l’aide d’un raisonnement FIG. 2.5 : La tension superficielle peut
énergétique : supposons qu’un opérateur déplace de façon quasi- s’interpréter comme une densité li-
néique de forces.
statique la tige de longueur ℓ en produisant une force ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹op . Si l’on
note d𝑥 le déplacement, le théorème de l’énergie cinétique appliqué
 16 2 TENSION DE SURFACE

à la tige T donne

dℰ𝑐 = 0 = ∑ travaux = 𝐹op d𝑥 + 𝛿𝑊liq→T

où 𝛿𝑊liq→T est le travail des forces capillaires sur la tige. Or, lors du
Lorsque que l’on traite une membrane déplacement de la tige la membrane reçoit une énergie
d’eau savonneuse, il ne faut pas ou-
blier qu’il y a deux interfaces liquide- dℰ = 2 𝛾 × (ℓd𝑥) = 𝑊T→liq = −𝑊liq→T
gaz ce qui explique la présence du fac-
teur deux.
En réinjectant 𝑊liq→T dans le théorème de l’énergie, on trouve 𝑓op =
2𝛾ℓ. Et comme 𝐹op = 𝐹 (évolution quasi-statique), on obtient

𝐹 = 2𝛾ℓ (2 interfaces) ♡ (2.3)

De façon général, on peut traiter une interface comme une membrane

tendue : chaque portion de surface est le siège de forces capillaires
réparties sur le contour 𝒞 délimitant la portion de surface. Ces forces
sont tangentes à l’interface, perpendiculaires en tout point de 𝒞 et
données par la relation

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 = 𝛾dℓ ⃗⃗⃗⃗
𝑛 avec ⃗⃗⃗⃗
𝑛⟂𝒞 ♡ (2.4)

Finalement, la tension superficielle est une force par unité de lon-

gueur, c’est pourquoi on l’exprime couramment en N.m−1 .

Application à la mesure de tension de surface

Mesure relative par stalagmométrie – La mesure consiste à compter

les gouttes sortant d’un compte-gouttes. Chaque goutte se détache
au niveau du col de rayon 𝑟 lorsque le poids du volume de la goutte
dépasse les forces capillaires. Si l’on note 𝑚 la masse de la goutte au
moment où elle se détache, on a

⃗⃗⃗⃗⃗⃗
𝐹 𝑚𝑔 = 2𝜋𝑟𝛾
2𝑅
2𝑟 En général on connaît le volume 𝑉 du compte-goutte et l’on compte
le nombre 𝑛 de gouttes. La relation précédente devient

𝜌𝑉
𝑚⃗⃗𝑔⃗⃗ 𝑔 = 2𝜋𝑟𝛾
𝑛
FIG. 2.6 : Stalagmométrie.
où 𝜌 est la masse volumique du liquide. Si l’on suppose que 𝑟 ne dé-
pend pas du liquide utilisé, mais seulement du diamètre du compte-
goute, alors on peut appliquer le même raisonnement pour un liquide
de référence (eau à 20°C en général) :

𝜌ref 𝑉ref
𝑔 = 2𝜋𝑟𝛾 (étalonnage)
𝑛ref

Le rapport des deux relations précédentes mène à la relation

𝜌/𝑛
𝛾 = 𝛾ref
𝜌ref /𝑛ref
 2.2 Forces capillaires 17

La connaissance de la masse volumique et du nombre de gouttes per-

met donc de comparer la tension de surface avec celle d’un liquide
de référence.

Mesure absolue par la méthode d’arrachement – Historiquement la

méthode de l’anneau a été la première à être développée. Il s’agit
de plonger un anneau (en platine en général) dans le liquide à étu-
dier puis de le remonter délicatement de façon à étirer un film au-
dessous de l’anneau. Au cours de l’étirement du film de liquide, la

dynamomètre

𝑔⃗⃗⃗⃗

anneau

boy

FIG. 2.7 : Méthode de l’anneau.

force exercée sur l’anneau est mesurée à l’aide d’un dynamomètre et

le système passe par un seuil où la force est maximale : dans ce cas
les forces capillaires sont verticales. Si l’on note 𝑟1 le rayon intérieur
de l’anneau et 𝑟2 son rayon extérieur, on a la relation.

𝐹max = 2𝜋(𝑟1 + 𝑟2 )𝛾 ≃ 4𝜋𝑟𝛾 avec 𝑟1 ≃ 𝑟2 = 𝑟

Sa mesure permet donc de déterminer la tension superficielle du

liquide.

Ordre de grandeur – À température ordinaire, pour les liquides mo-

léculaire, la tension superficielle vaut quelques dizaines de mN.m−1 .
Les liquides métalliques ont des tensions superficielles un ordre de
grandeur au dessus.

Tension superficielle de quelques interfaces liquide-air

Liquides Température Tension superficielle [mN.m−1 ]
Benzène 20 °C 29
Eau 20 °C 73
Eau 80 °C 62
Glycérine 20 °C 63
Hexadécane 20 °C 27
Huile d’olive 20 °C 32
Mercure 18 °C 475
Plomb 327 °C 462

La tension superficielle varie avec la température. Pour les corps purs,

𝛾 diminue linéairement avec la température jusqu’à s’annuler à la
 18 2 TENSION DE SURFACE

température critique

𝑇
𝛾 = 𝛾0 (1 − ) (formule d’Eötvös)
𝑇c

Théorème de Laplace

Membrane liquide Un petit contour pris dans la surface libre plane d’un liquide en équi-
libre est soumis à des forces de tension superficielle situées dans
son plan et dont la résultante est nulle.

Par contre, pour une surface sphérique, les forces exercées sur ce
FIG. 2.8 : La résultante des forces ca-
pillaires qu’exerce la membrane sur même contour ont une résultante orientée vers l’intérieur de la sphère ;
un contour fermé est nulle si la mem- il faut donc une surpression Δ𝑝 pour que l’équilibre existe. On voit
brane est plane. immédiatement que plus la courbure est importante et plus Δ𝑝 sera
grand.

Système
𝑝ext

𝐹𝑝 ext
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
FIG. 2.9 : À gauche : existence d’une
surpression dans une goutte. La cour- ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 𝐹𝑝 int
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
bure entraîne l’existence d’une résul-
tante des forces capillaires dirigée vers ⃗⃗⃗⃗⃗⃗cap
𝐹
le centre de courbure. À droite : forces
agissantes sur une portion hémisphé-
rique d’une goutte liquide. Les flèches 𝑝int 𝑝ext 𝑝int
réparties sur le contour circulaire re-
présentent les forces capillaires.

Calculons la surpression Δ𝑝 qui règne à l’intérieur d’une goutte sphé-

rique de rayon 𝑅 en faisant un bilan des forces sur une demi-goutte :

• La pression extérieure produit une force pressante

⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝑝 ext = − ∬ 𝑝ext⃗⃗⃗⃗
𝑛 d𝑆 = −𝑝ext 𝜋𝑅2 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧

• La pression intérieure produit une force pressante

⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝑝 int = 𝑝int 𝜋𝑅2 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧

• Enfin, l’hémisphère inférieure exerce sur le contour du système

des forces capillaires orientées vers le bas dont la résultante
s’écrit
⃗⃗⃗ ⃗⃗cap = ∮ 𝛾dℓ(−⃗⃗⃗⃗⃗⃗⃗⃗
𝐹 𝑢𝑧 ) = −2𝜋𝑅𝛾 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧

L’équilibre du système se traduit par

(𝑝int − 𝑝ext ) 𝜋𝑅2 − 2𝜋𝑅𝛾 = 0

ce qui donne une surpression

2𝛾
Δ𝑝 = 𝑝int − 𝑝ext = ♡ (2.5)
𝑅
 2.2 Forces capillaires 19

De la même manière, à l’intérieur d’une bulle de savon, il règne une

surpression
4𝛾
Δ𝑝 = 𝑝int − 𝑝ext =
𝑅
où le facteur 4 est dû au fait que la bulle de savon présente deux
interfaces liquide-gaz.

Ordre de grandeur – Pour une bulle de savon de rayon 𝑅 ≃ 1 cm

et 𝛾 ≃ 25.10−3 N.m−1 on obtient Δ𝑝 ≃ 10 Pa. Pour faire des grosses
bulles il faut fournir beaucoup d’énergie (surface importante) et gé-
nérer une faible surpression ; il faut donc souffler tout doucement et
longtemps.

Mûrissement d’une mousse – Dans une mousse humide (mousse à

raser, mousse de bière, etc.), du gaz est enfermé dans des bulles
sphériques séparées par un film liquide pouvant plus ou moins lais-
ser diffuser le gaz selon l’épaisseur de la membrane et la taille des
molécules gazeuses. Le gaz emprisonné est en surpression par rap-
port au liquide en vertu de la loi (2.5). Par ailleurs la surpression est
plus importante dans les petites bulles. C’est ce qui explique le phé-
nomène de mûrissement d’une mousse : Le gaz contenu dans les
petites bulles traversent la membrane liquide par diffusion pour se
diriger dans les zones de moins grande pression, c’est-à-dire, dans
les grosses bulles. Les petites bulles se vident donc dans les grosses
et la mousse s’enrichit en grosses bulles.

La généralisation à une géométrie quelconque est donnée par la loi

de Laplace-Young :

Loi de Laplace-Young

La différence de pression à la traversée d’une interface séparant

un fluide intérieur et un fluide extérieur est donnée par

1 1
𝑝int − 𝑝ext = 𝛾 ( + )
𝑅1 𝑅2

où 𝑅1 et 𝑅2 sont les rayons de courbure de l’interface selon deux

directions orthogonales. Par convention, ces rayons sont positifs
quand le centre de courbure est dans le fluide intérieur. Dans le
cas d’une interface sphérique, ces deux rayons de courbure s’iden-
tifient au rayon de la sphère.

Surface minimale – Lorsque l’on trempe une structure métallique

dans une eau de savon, on obtient une surface minimale (l’interface
va chercher a minimiser l’énergie superficielle) qui a la propriété sui-
vante : si la surface est ouverte, Δ𝑝 = 0 et donc

1 1
+ =0
𝑅1 𝑅2

On dit que la courbure moyenne est nulle. Dans la plupart des cas on
obtient des lames planes qui forment une surface minimale (𝑅1 , 𝑅2 →
∞). On peut aussi obtenir des lames avec deux rayons de courbures FIG. 2.10 : Lame de savon formant une
caténoïde (©Berkeley Science Review).
opposées comme sur la FIG. 2.10 montrant une caténoïde.
 20 2 TENSION DE SURFACE

2.3 Mouillage

Angles de contact

Équilibre au contact de trois fluides – Déposons une petite quantité

de liquide 2 sur un autre liquide 1 plus dense et non miscible. L’en-
semble des points en contact avec les deux liquides et l’air forme un
contour appelé ligne triple. Intéressons nous aux forces capillaires
s’exerçant sur cette ligne triple. Notons ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝛾𝑖𝑗 la force capillaire par uni-
té de longueur due à l’interface entre les fluides 𝑖 et 𝑗. L’équilibre n’est
gaz ou liq.3
possible que si la résultante des forces capillaires peut s’annuler ce
qui définit l’angle de contact 𝜃 :
⃗⃗𝛾⃗ ⃗23
liq.2
⃗⃗𝛾⃗ ⃗13
𝜃 ⃗⃗𝛾⃗⃗12 + ⃗⃗𝛾⃗⃗23 + ⃗⃗𝛾⃗⃗31 = 0⃗
⃗⃗𝛾⃗ ⃗12
liq.1 Cet équilibre suppose que l’on puisse former un triangle avec les trois
vecteurs ⃗⃗𝛾⃗⃗𝑖𝑗 ce qui n’est possible qu’à condition que chaque tension
de surface soit inférieure à la somme des deux autres. Dans le cas
FIG. 2.11 : Définition de l’angle de contraire, le liquide 2 s’étale sur le liquide 1 : on dit qu’il y a mouillage
contact.
total.

Exercice – Déposons une goutte d’huile d’olive sur de l’eau. Sachant que
𝛾eau-air = 73 mN.m−1 , 𝛾huile-air = 32 mN.m−1 et 𝛾huile-eau = 18 mN.m−1 , dire s’il
y a étalement ou non.
Solution : l’huile d’olive s’étale.

Équilibre d’un liquide au contact d’un solide – Déposons une goutte

de liquide sur un support plan. En général, le liquide adopte la forme
décrite sur la FIG. 2.12, résultat d’un compromis entre le poids qui
tend à diminuer la position du centre de gravité de la goutte et des
forces capillaires qui tendent à minimiser l’aire de la surface libre. À
l’équilibre, la résultante des forces capillaires en un point de la ligne
triple s’annule. En projection sur le plan de la surface solide, on a
donc
𝛾𝑆𝐿 + 𝛾𝐿𝐺 cos 𝜃 = 𝛾𝑆𝐺
Cette relation trouvée par Young en 1805 définit l’angle de contact 𝜃.
On distingue trois cas de figure :

⃗⃗𝛾⃗ ⃗𝐿𝐺

𝜃
⃗⃗𝛾⃗ ⃗𝑆𝐺 ⃗⃗𝛾⃗ ⃗𝑆𝐿 𝜃 > 𝜋/2 𝜃 < 𝜋/2 𝜃→0
Support absence de mouillage mouillage partiel mouillage complet

FIG. 2.12 : Différents équilibres d’un liquide au contact d’un substrat solide.

• 𝜃 > 𝜋/2 : le liquide est non mouillant (exemple : verre-mercure-

air)
• 𝜃 ∈ [0, 𝜋/2] : il y a mouillage partiel (exemple : verre-eau-air)
• Lorsque l’angle de contact n’est pas défini, il y a mouillage com-
plet du liquide sur le substrat solide.
 2.4 Approche thermodynamique 21

Ascension capillaire

Quand on plonge un tube fin3 de rayon 𝑟 dans de l’eau, on observe 3 : Ces tubes ont la finesse d’un
l’ascension d’une colonne d’eau dans le capillaire malgré la pesan- cheveu d’où l’appellation tubes capil-
laires.
teur. Cette ascension est d’autant plus importante que le rayon est
petit.
Cte
ℎ= (loi de Jurin)
𝑟
où la constante dépend du liquide et de l’angle de contact. On peut

𝜃
2𝑟

FIG. 2.13 : Expérience et schéma.

utiliser la loi de Laplace-Young pour démontrer la loi de Jurin. En effet,

considérons l’interface liquide raccordée à la paroi du tube avec un
angle 𝜃 et supposons que le ménisque est sphérique4 de rayon 𝑅. 𝛾
4 : Ce qui suppose 𝑟 ≪ ℓc = √ 𝜌𝑔
L’air étant à la pression 𝑝0 , la pression qui règne dans le liquide au
voisinage du ménisque vaut 𝑝0 − Δ𝑝 avec Δ𝑝 = 2𝛾/𝑅 en vertu de la
loi de Laplace-Young. Le rayon de courbure vaut 𝑅 = 𝑟/ cos 𝜃 d’où

2𝛾 cos 𝜃
𝑝ℎ = 𝑝0 −
𝑟
Or, si l’on applique les lois de l’hydrostatique au niveau de la surface
libre du récipient dans lequel plonge le tube, on trouve

2𝛾 cos 𝜃 1
𝑝0 = 𝑝ℎ + 𝜌𝑔ℎ ⇒ ℎ=
𝜌𝑔 𝑟

Il y a donc ascension capillaire si le liquide mouille la paroi (𝜃 < 𝜋/2).

Par contre pour un liquide non mouillant il y a descente capillaire
(cas du mercure dans un capillaire en verre).

Exercice – On plonge côte à côte deux tubes capillaires de diamètres inté-

rieurs 𝑑1 = 0,20 mm et 𝑑2 = 0,40 mm dans de l’huile d’olive. On mesure une
différence de hauteur d’ascension de 31 mm. En déduire la tension super-
ficielle de l’huile d’olive en supposant que l’huile mouille parfaitement les
parois du tube capillaire. On donne 𝜌huile = 918 g ⋅ L−1
Solution : 𝛾huile = 28 mN ⋅ m−1 .

2.4 Approche thermodynamique

Dans toute la suite on note 𝐴 l’aire de l’interface au lieu de 𝑆, notation
que l’on réserve à l’entropie.
 22 2 TENSION DE SURFACE

Identité thermodynamique

Soit un système liquide que l’on décrit par son volume 𝑉 , son aire 𝐴
sa température 𝑇 et sa pression 𝑝. Appelons également 𝑈 son éner-
gie interne, et 𝑆 son entropie. Imaginons une transformation infini-
tésimale réversible pendant laquelle le système reçoit un transfert
thermique 𝛿𝑄rev et un travail 𝛿𝑊 rev = −𝑝 d𝑉 + 𝛾d𝐴 où 𝛾 est la ten-
sion superficielle.

En vertu du second principe, on a 𝛿𝑄rev = 𝑇 d𝑆, et le premier principe

donne l’identité thermodynamique

d𝑈 = 𝛿𝑄rev + 𝛿𝑊 rev soit d𝑈 = 𝑇 d𝑆 − 𝑝d𝑉 + 𝛾d𝐴

Introduisons maintenant l’énergie libre, une grandeur notée 𝐹 (comme

free energy) et définie par

def
𝐹 = 𝑈 − 𝑇𝑆

Sa différentielle s’écrit

d𝐹 = d𝑈 − 𝑇 d𝑆 − 𝑆d𝑇 = −𝑆d𝑇 − 𝑝d𝑉 + 𝛾d𝐴

identité qui montre que l’énergie libre est une fonction de 𝑇 , 𝑉 et

𝐴.

Si maintenant on fixe la température et l’on suppose le liquide in-

compressible (𝑉 est donc constant), on obtient d𝐹 = 𝛾d𝐴.

À retenir
Le terme 𝛾d𝐴 représente la variation d’énergie libre du liquide
quand on augmente son aire à température constante.

Exemple de calcul thermodynamique : la loi de Laplace

Essayons de retrouver la formule de Laplace relative à une bulle de

𝑝0 , 𝑇0 savon à l’aide d’une analyse thermodynamique.
membrane 𝑅
Considérons une bulle de savon de rayon 𝑅, enfermant un volume 𝑉
air 𝑝, 𝑇0 d’air et présentant une interface d’aire 𝐴 = 2×4𝜋𝑅2 . Prenons comme
système l’ensemble {volume d’air + membrane de savon}.
FIG. 2.14 : Bulle de savon enfermant de
l’air à la pression 𝑝. On impose une température extérieure 𝑇0 et une pression 𝑝0 . Dans
ce cas, on sait qu’à l’équilibre l’enthalpie libre 𝐺 = 𝑈 + 𝑝0 𝑉 − 𝑇0 𝑆
sera minimale.
Rappel On détermine donc l’équilibre en résolvant l’équation d𝐺 = 0. On a
Pour une évolution monotherme ici
et monobare, le second principe d𝐺 = d𝑈 + 𝑝0 d𝑉 − 𝑇0 d𝑆
stipule que Δ𝑆 ≥ 𝑇𝑄 et le premier
0
principe Δ𝑈 = 𝑄 − 𝑝0 Δ𝑉 si les avec, en vertu de l’additivité de l’entropie et de l’énergie
seules forces extérieures sont les
forces de pression. Il est assez fa- d𝑈 = d𝑈air + d𝑈membrane et d𝑆 = d𝑆air + d𝑆membrane
cile de montrer alors que Δ𝐺 ≤ 0 :
l’enthalpie libre décroît jusqu’à at-
teindre un minimum à l’équilibre.
 2.4 Approche thermodynamique 23

Le système est à la température 𝑇0 , et l’on note 𝑝 la pression de l’air

dans la bulle que l’on cherche à déterminer. On a donc

d𝑈air = 𝑇0 d𝑆air − 𝑝d𝑉 et d𝑈membrane = 𝑇0 d𝑆membrane + 𝛾d𝐴

En rassemblant tous les termes et en remplaçant 𝑉 et 𝐴 par leur

expression, on trouve

d𝐺 = [(𝑝0 − 𝑝)4𝜋𝑅2 + 16𝜋𝑅𝛾] d𝑅

Ainsi, l’équilibre est obtenu quand d𝐺 = 0 soit lorsque

4𝛾
𝑝 − 𝑝0 = (CQFD)
𝑅
 FLUIDE EN ÉCOULEMENT 3
Dans ce cours, nous étudions le fluide et son écoulement indépen- 3.1 Description d’un fluide en
damment des forces qui en sont responsables. écoulement . . . . . . . . 25
Deux points de vue . . . . 25
Écoulement stationnaire . 26
Version en ligne
Accélération . . . . . . . . 27
[Link]/mecanique_des_fluides/[Link]
Dérivée particulaire . . . 28
Techniques de visualisation28
3.2 Conservation de la masse 29
3.1 Description d’un fluide en écoulement Débit massique . . . . . . 29
Débit volumique . . . . . 30
Équation de continuité . . 30
Deux points de vue 3.3 Écoulement incompressible 31
Dilatation . . . . . . . . . 31
Tout d’abord, comme on l’a déjà dit, on décrit le fluide comme un en- Conservation du débit
semble d’une infinité de particules de fluide. On considère, autour volumique . . . . . . . . . 32
Sa dimension doit être grande devant
d’un point P, un volume mésoscopique d𝜏 . Ce volume contient un le libre parcours moyen des molécules.
grand nombre 𝑁 de molécules ce qui permet de définir des gran- Dans les conditions standards, un vo-
deurs moyennes locales qui, elles, vont évoluer de façon continue. lume de 1 µm3 convient.
On définit alors la vitesse locale 𝑣 ⃗ d’une particule de fluide P par

1 𝑁
𝑣P =
⃗⃗⃗⃗⃗⃗⃗ ∑ ⃗⃗⃗⃗⃗⃗
𝑣 = ⟨⃗⃗⃗⃗⃗⃗
𝑣𝑖 ⟩
𝑁 𝑖=1 𝑖

où ⃗⃗⃗⃗⃗⃗
𝑣𝑖 est le vecteur vitesse d’une molécule située dans le volume d𝜏 .
𝑣P = 0⃗ bien que ⃗⃗⃗⃗⃗⃗
Ainsi, lorsque le fluide est au repos, ⃗⃗⃗⃗⃗⃗⃗ 𝑣𝑖 ≠ 0.⃗

Deux approches différentes existent pour décrire le fluide en mouve-

ment.
Trajectoire de P
Point de vue de Lagrange – On s’intéresse à l’histoire des particules
P(𝑡0 )
de fluide en suivant leur trajectoire. Chaque trajectoire est donc as-
sociée à une particule donnée (que l’on notera P). Au cours du trajet,
le vecteur vitesse est tangent à la trajectoire, et l’accélération de la P(𝑡1 )
𝑣
⃗⃗⃗⃗⃗⃗⃗⃗P (𝑡0 )
particule s’écrit P(𝑡2 )
𝑣
⃗⃗⃗⃗⃗⃗⃗⃗P (𝑡1 )
d⃗⃗⃗⃗⃗⃗⃗
𝑣
⃗⃗⃗ ⃗P⃗⃗ = P
𝑎 𝑣
⃗⃗⃗⃗⃗⃗⃗⃗P (𝑡2 )
d𝑡
conformément aux lois vues en mécanique du point.
FIG. 3.1 : Notion de trajectoire.
Point de vue d’Euler – Dans cette approche, on s’intéresse au fluide
dans son ensemble à un instant 𝑡. Pour chaque point du fluide (no-
tons le M) on observe le vecteur vitesse à l’instant 𝑡. Comme ce vecteur
dépend du temps 𝑡 et de la position du point considéré, on note cette M1

grandeur 𝑣(M
⃗⃗⃗⃗ 1 , 𝑡)
𝑣(M,
⃗ 𝑡) ou 𝑣(𝑥, ⃗ 𝑦, 𝑧, 𝑡)

M2
Du point de vue mathématique, il s’agit d’un champ vectoriel spatio-
temporel. Pour représenter le champ de vitesse, on utilise la notion 𝑣(M
⃗⃗⃗⃗ 2 , 𝑡)
de ligne de champ : à un instant 𝑡 fixé, on trace différentes courbes qui
ont la particularité d’avoir leur tangente en un point colinéaire à la
vitesse en ce même point, ceci pour tous les points de la courbe. Ces Ligne d’écoulement

FIG. 3.2 : Notion de ligne d’écoulement.

Notez que toutes les vitesses représen-
tées concernent le même instant.
26 3 FLUIDE EN ÉCOULEMENT

courbes sont appelées lignes de courant ou lignes d’écoulement. No-

tez que la ligne d’écoulement est fondamentalement différente d’une
trajectoire, car tous les points de la ligne sont observés au même ins-
tant, contrairement à la trajectoire.
Pour un problème à deux dimensions, l’équation 𝑓(𝑥, 𝑦) = 0 d’une
ligne d’écoulement s’obtient en résolvant l’équation différentielle

d𝑦 𝑣𝑦 (𝑥, 𝑦, 𝑡)
= avec 𝑡 fixe
d𝑥 𝑣𝑥 (𝑥, 𝑦, 𝑡)

Exercice – On considère un écoulement bidimensionnel dont le champ de

vitesse s’écrit :

𝑣(M,
⃗ 𝑡) = −𝑘𝑥 ⃗⃗⃗⃗⃗⃗⃗⃗ 𝑢𝑦 avec 𝑘 = Cte
𝑢𝑥 + 𝑘𝑦 ⃗⃗⃗⃗⃗⃗⃗⃗

• Quelle est l’équation des lignes d’écoulement ?

• La vitesse est-elle constante le long d’une ligne de courant ?

Écoulement stationnaire
Définition
1 : On parle aussi d’écoulement per- Un écoulement est stationnaire1 si le champ de vitesse ne dépend
manent pas explicitement du temps 𝑡. Mathématiquement cela se traduit
pas
𝜕 𝑣(𝑥,
⃗ 𝑦, 𝑧, 𝑡)
= 0⃗
𝜕𝑡

Cela signifie que :

• l’expression du champ de vitesse ne dépend que de 𝑥, 𝑦 et 𝑧
(en coordonnées cartésiennes) ;
• la carte du champ de vitesse est invariable au cours du temps ;
• les lignes de courant sont fixes dans le temps.

Attention, cela ne signifie pas que les particules se déplacent à vitesse

constante.

Exemple

Supposons un fluide en écoulement présentant un champ de vitesse

𝑥
𝑣(M,
⃗ 𝑡) = 𝑢 avec 𝜏 > 0
⃗⃗⃗⃗⃗⃗⃗⃗
𝜏 𝑥

L’écoulement est stationnaire puisque 𝜕 𝑣(M,

⃗ 𝑡)/𝜕𝑡 = 0.⃗ De plus, en tout
point, le vecteur vitesse à la même direction. Par conséquent, les lignes
d’écoulement sont des droites parallèles de vecteur directeur ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥 .

Intéressons-nous à la vitesse d’une particule de fluide en particulier. Ap-

pelons P la particule située en 𝑥0 à l’instant 𝑡 = 0 et notons 𝑥P (𝑡) et
𝑣P (𝑡) respectivement sa position et sa vitesse au cours du temps. Le lien
⃗⃗⃗⃗⃗⃗⃗
entre la « vitese eulérienne » 𝑣(M,
⃗ 𝑡) et la « vitesse lagrangienne » ⃗⃗⃗⃗⃗⃗⃗
𝑣P (𝑡) se
 3.1 Description d’un fluide en écoulement 27

traduit par
𝑥P d𝑥P
𝑣(𝑃 𝑣P (𝑡) soit
⃗ (𝑡), 𝑡) = ⃗⃗⃗⃗⃗⃗⃗ =
𝜏 d𝑡
équation différentielle d’ordre un dont la solution s’écrit

d𝑥P 𝑥
𝑥P = 𝑥0 e𝑡/𝜏 d’où 𝑣P = = 0 e𝑡/𝜏
d𝑡 𝜏
Comme on le voit, la vitesse de P varie au cours du temps bien que l’écou-
lement soit stationnaire.

Sur l’exemple précédent, on constate également que la trajectoire

d’une particule se confond avec la ligne d’écoulement. Cette propriété
se généralise pour tous les écoulements stationnaires.

À retenir
En régime stationnaire, les trajectoires tracent les lignes d’écoule-
ment.

Accélération d’une particule

Supposons que l’on connaisse le champ de vitesse eulérien 𝑣(M, ⃗ 𝑡)

et que l’on veuille déterminer l’accélération 𝑎(M,⃗ 𝑡), c’est-à-dire l’ac-
célération d’une particule de fluide située à l’instant 𝑡 au même point
M.
La difficulté du calcul réside dans le fait que si une particule se situe
en M à l’instant 𝑡, elle ne situe plus au même point à l’instant 𝑡 + 𝛿𝑡.
C’est pourquoi, on ne peut pas écrire 𝑎(M,⃗ 𝑡) = 𝜕 𝑣(M,
⃗ 𝑡)/𝜕𝑡.
Considérons une particule P située en M(𝑥, 𝑦, 𝑧) à l’instant 𝑡. Elle se
trouve alors en M’(𝑥+𝛿𝑥, 𝑦 +𝛿𝑦, 𝑧 +𝛿𝑧) à l’instant 𝑡+𝛿𝑡. Par définition,
la particule présente une accélération

𝑣(M’,
⃗ 𝑡 + 𝛿𝑡) − 𝑣(M,
⃗ 𝑡)
𝑎⃗ = lim
𝛿𝑡→0 𝛿𝑡
𝑣(𝑥
⃗ + 𝛿𝑥, 𝑦 + 𝛿𝑦, 𝑧 + 𝛿𝑧, 𝑡 + 𝛿𝑡) − 𝑣(𝑥,
⃗ 𝑦, 𝑧, 𝑡)
= lim
𝛿𝑡→0 𝛿𝑡

En utilisant le développement de Taylor à l’ordre un, on obtient

𝜕𝑣⃗ 𝜕𝑣⃗ 𝜕𝑣⃗ 𝜕𝑣⃗

𝑣(𝑥+𝛿𝑥,
⃗ 𝑦+𝛿𝑦, 𝑧+𝛿𝑧, 𝑡+𝛿𝑡) = 𝑣(𝑥,
⃗ 𝑦, 𝑧, 𝑡)+𝛿𝑥 +𝛿𝑦 +𝛿𝑧 +𝛿𝑡
𝛿𝑡→0 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡

De sorte que

𝛿𝑥 𝜕 𝑣 ⃗ 𝛿𝑦 𝜕 𝑣 ⃗ 𝛿𝑧 𝜕 𝑣 ⃗ 𝜕 𝑣 ⃗
𝑎⃗ = lim + lim + lim +
⏟
𝛿𝑡→0 𝛿𝑡 𝜕𝑥 ⏟
𝛿𝑡→0 𝛿𝑡 𝜕𝑦 𝛿𝑡→0 𝛿𝑡 𝜕𝑧
⏟ 𝜕𝑡
𝑣𝑥 𝑣𝑦 𝑣𝑧

Expression que l’on peut mettre sous forme compacte

𝜕𝑣⃗ 𝜕 𝜕 𝜕
𝑎⃗ = + (𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧 ) 𝑣 ⃗
𝜕𝑡 ⏟⏟⏟⏟⏟⏟⏟⏟⏟⏟⏟
𝜕𝑥 𝜕𝑦 𝜕𝑧
⃗⃗⃗⃗ ⃗⃗⃗⃗⃗⃗⃗
𝑣⋅ ∇
 28 3 FLUIDE EN ÉCOULEMENT

Champ d’accélération

Si une particule de fuide présente un vecteur vitesse 𝑣(M,⃗ 𝑡) lors-

qu’il est en M à l’instant 𝑡, alors son accélération vaut

𝜕𝑣⃗
𝑎(M,
⃗ 𝑡) = + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗ (3.1)
𝜕𝑡

Le premier terme est lié au caractère non permanent de l’écoulement

alors que le second au fait que la particule, en se déplaçant, visite des
endroits où la vitesse change. On l’appelle le terme convectif.

Exercice – On considère un écoulement dont le champ de vitesse s’écrit

𝑥
𝑣(M,
⃗ 𝑡) = 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗
𝑡+1 𝑥
• L’écoulement est-il stationnaire ?

Dérivée particulaire

Nous venons de calculer le taux de variation du vecteur vitesse lorsque

l’on suit une particule. Une telle grandeur est appelée dérivée parti-
culaire de la vitesse. Nous la notons
D𝑣 ⃗
D𝑡
Le raisonnement effectuée pour la vitesse peut s’appliquer à toute
grandeur physique 𝐺(M, 𝑡) attachée à une particule de fluide située
en M à l’instant 𝑡. Lorsque l’on suit la particule, le taux de variation
de la grandeur 𝐺 vaut

D𝐺 𝜕𝐺
= + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇)𝐺 ♡ (3.2)
D𝑡 𝜕𝑡

Techniques de visualisation

Lignes d’écoulement – On utilise des particules réfléchissantes que

l’on photographie avec un court temps de pose. On a accès ainsi à
des segments brillants qui donnent la direction du vecteur vitesse
en différents points ce qui permet de reconstituer la carte du champ
de vitesse.
Lignes d’émission – Courbe décrivant la position à un instant donné
de l’ensemble des particules d’un fluide qui sont passées antérieure-
ment en un même point donné. Expérimentalement on les visualise
FIG. 3.3 : Visualisation des tourbillons
derrière un obstacle plan.
en émettant en continue un traceur en différentes points et en pre-
nant une photographie instantanée.
Carte de densité – La strioscopie est une technique interférométrique
permettant d’accéder aux variations d’indice de réfraction et donc
aux variations de densité.
 3.2 Conservation de la masse 29

FIG. 3.4 : À gauche : lignes d’émission autour d’un corps fuselé à l’aide de la machine à fumées d’Etienne Marey. À droite : ondes
de choc lors d’un vol supersonique visualisées par strioscopie.

3.2 Conservation de la masse

Débit massique

Un écoulement est un phénomène de transport puisqu’il se traduit

par un transfert de masse. Et comme pour tout phénomène de trans-
port, on est amené à définir la notion de flux et de courant.

Commençons par le flux de masse, que l’on appelle plus communé-

ment débit massique. Par définition, le débit massique, que l’on note
𝐷𝑚 , est la masse qui traverse une surface par unité de temps. Plus
précisément, si d𝑚 est la masse qui traverse la surface (S) entre 𝑡 et
𝑡 + d𝑡, on a
d𝑚 = 𝐷𝑚 d𝑡
♡ (3.3)
(kg) (kg ⋅ s−1 ) (s)

Le débit massique peut se calculer à l’aide d’une intégrale de surface2 2 : appelée aussi flux mathématique
sous la forme
𝐷𝑚 = ∬ ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑗𝑚 ⋅ ⃗⃗⃗⃗ 𝑛 d𝑆 ♡ (3.4)
(𝑆)

où ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑗𝑚 désigne le vecteur courant de masse. Cherchons à exprimer le
vecteur courant en fonction des caractéristiques du fluide en écoule-
ment.

Pour cela, considérons une section infinitésimale d𝑆 autour d’un point

2
M et calculons la masse d 𝑚 de fluide traversant d𝑆 pendant d𝑡. Cette
masse se trouve dans le prisme de base d𝑆 et de génératrice 𝑣 ⃗ d𝑡. On
a donc 𝑡
2 𝑣⃗⃗ ⃗d
d 𝑚 = 𝜌 × (d𝑆 × 𝑣 d𝑡 cos 𝛼) , 𝑡) 𝛼
𝑣⃗⃗ ⃗(M
𝑛 d𝑡 d𝑆
= 𝜌𝑣 ⃗ ⋅ ⃗⃗⃗⃗ • ⃗⃗⃗⃗⃗
𝑛
M
où ⃗⃗⃗⃗
𝑛 est le vecteur normal à la section d𝑆. Pour une surface macro-
scopique (S), il suffit de sommer toutes les contributions : dS

d𝑚 = (∬ 𝜌(M, 𝑡)𝑣(M,
⃗ 𝑛d𝑆) d𝑡
𝑡) ⋅ ⃗⃗⃗⃗
M∈(𝑆)

FIG. 3.5 : Calcul du débit massique.

 30 3 FLUIDE EN ÉCOULEMENT

On en déduit

def
𝐷𝑚 = ∬ ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑗𝑚 (M, 𝑡) ⋅ ⃗⃗⃗⃗
𝑛d𝑆 avec ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑗𝑚 (M, 𝑡) = 𝜌(M, 𝑡)𝑣(M,
⃗ 𝑡) ♡
(𝑆)

(3.5)

Débit volumique

Divisons l’équation (3.1) par la masse volumique :

d𝑚 𝐷
= d𝑉 = 𝑚 d𝑡
𝜌 𝜌

où d𝑉 représente le volume transportée pendant d𝑡. On en déduit le

débit volumique

def d𝑉
𝐷𝑉 = d𝑡 = ∬M∈(𝑆) 𝑣(M,
⃗ 𝑡) ⋅ ⃗⃗⃗⃗
𝑛 d𝑆
♡ (3.6)
m3 ⋅ s−1 m ⋅ s−1 m2

Ainsi le flux du vecteur vitesse à travers une surface, donne le débit

volumique qui le traverse.

Équation de continuité

Établissons la première équation fondamentale de la mécanique des

fluides. Il s’agit d’une contrainte imposée à 𝑣(M,
⃗ 𝑡) et 𝜌(M, 𝑡) qui tra-
duit la conservation de la masse.

Prenons un volume de contrôle fixe (𝑉 ) dans un fluide, délimité par

une surface (𝑆). Soit 𝑚(𝑡) la masse contenue à l’intérieur de la surface
𝑛ext
⃗⃗⃗⃗⃗ fermée à l’instant 𝑡. Par définition de la masse volumique,
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑗𝑚 d𝑆

𝑚(𝑡) = ∭ 𝜌(M, 𝑡) d𝜏
M∈(𝑉 )

(S)
Cette masse varie à cause du flux de matière à travers (𝑆) :

FIG. 3.6 : Calcul du débit massique sor- d𝑚(𝑡) sortant

tant d’une surface fermée. = −𝐷𝑚 = −∯ 𝑛ext d𝑆
𝜌𝑣 ⃗ ⋅ ⃗⃗⃗⃗
d𝑡 M∈(𝑆)

𝑛ext est dirigé vers l’extérieur de la surface fermée (𝑆) ce qui ex-
où ⃗⃗⃗⃗
plique l’origine du signe – devant l’intégrale. Or, on a également

d𝑚(𝑡) 𝜕𝜌(M, 𝑡)
=∭ d𝜏
d𝑡 M∈𝑉 𝜕𝑡
D’après le théorème de la divergence on obtient :

𝜕𝜌
∭ [div(𝜌𝑣)⃗ + ] d𝜏 = 0 ∀𝑉
M∈(𝑉 ) 𝜕𝑡
 3.3 Écoulement incompressible 31

d’où l’équation de conservation de la masse, dite aussi équation de

continuité

𝜕𝜌
div(𝜌𝑣)⃗ + =0 partout et à chaque instant ♡ (3.7)
𝜕𝑡

Théorème de Green-Ostrogradsky ou théorème de la divergence

Le flux d’un champ vectoriel 𝐴(M) ⃗⃗⃗ ⃗⃗ à travers une surface fermée
(𝑆) est égal à l’intégrale sur le volume 𝑉 limité par (𝑆) de la diver-
gence du champ vectoriel.

∯ ⃗⃗⃗ ⃗⃗
𝐴(M) 𝑛ext d𝑆 = ∭
⋅ ⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗
div𝐴(M) d𝜏
M∈(𝑆) M∈𝑉

avec
𝜕𝐴𝑦
⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗
div𝐴 ∇⋅𝐴 ⃗⃗⃗ ⃗⃗ = 𝜕𝐴𝑥 + +
𝜕𝐴𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

Exercice – Que devient l’équation de continuité pour un liquide incompres-

sible ?

3.3 Écoulement incompressible

Dilatation d’une particule de fluide

La quantité div(𝑣)⃗ prend une signification bien précise en mécanique

des fluides. Partons de la relation

⃗⃗⃗ ⃗⃗ = 𝑓div𝐴
div(𝑓.𝐴) ⃗⃗⃗ ⃗⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
⃗⃗⃗ ⃗⃗ + 𝐴. grad𝑓

Appliqué au vecteur densité de courant de matière ⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝐽m = 𝜌𝑣 ⃗ cela

donne
div(𝜌𝑣)⃗ = 𝜌div𝑣 ⃗ + 𝑣.⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝜌
En utilisant l’équation de continuité on obtient

𝜕𝜌
+ 𝑣.⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝜌 = −𝜌div𝑣 ⃗
𝜕𝑡
On reconnaît dans le terme de gauche, la dérivée particulaire de 𝜌.

1 D𝜌
div𝑣 ⃗ = −
𝜌 D𝑡

Si l’on note 𝛿𝑚 et 𝛿𝜏 la masse et le volume d’une particule de fluide

en mouvement on peut écrire

1 D𝜌 1 D (𝛿𝑚/𝛿𝜏 ) 𝛿𝑚 D𝛿𝜏 1 D𝛿𝜏

− =− = =
𝜌 D𝑡 𝜌 D𝑡 𝜌𝛿𝜏 2 D𝑡 𝛿𝜏 D𝑡
 32 3 FLUIDE EN ÉCOULEMENT

Finalement, on obtient

1 𝐷𝛿𝜏
div𝑣 ⃗ = ♡ (3.8)
𝛿𝜏 D𝑡

La divergence de la vitesse d’écoulement représente ainsi la vitesse

de dilatation de la particule de fluide.

Exemple 1 : écoulement unidimensionnel uniforme – Considérons

l’écoulement décrit par le champ de vitesse

𝑣(M,
⃗ 𝑢𝑥 avec 𝑣 = Cte
𝑡) = 𝑣 ⃗⃗⃗⃗⃗⃗⃗⃗

Les lignes de courant sont des droites parallèles et l’écoulement est à

divergence nulle. Les particules de fluide se déplacent sans se dilater
comme le montre la FIG. 3.7.

FIG. 3.7 : Particule de fluide dans un écoulement uniforme, capturé à différents instants.

Exemple 2 : écoulement radial – Considérons l’écoulement décrit en

coordonnées polaires par le champ de vitesse

𝑣(M,
⃗ 𝑢𝑟 avec 𝑣 = Cte
𝑡) = 𝑣 ⃗⃗⃗⃗⃗⃗⃗⃗

Les lignes de courant sont des droites issues de O et les particules

de fluide se déplacent en se dilatant comme le montre la FIG. 3.8, ce
qui prouve que l’écoulement est à divergence positive.

FIG. 3.8 : Particule de fluide dans un écoulement radial, capturé à différents instants.

Conservation du débit volumique

Définition
Un fluide est en écoulement incompressible quand les particules
de fluide se déplacent sans subir, ni dilatation, ni contraction.
Donc,
Écoulement incompressible ⇒ div𝑣 ⃗ = 0

Par conséquent, d’après le théorème de la divergence, un fluide en

 3.3 Écoulement incompressible 33

écoulement incompressible vérifie la relation

𝑛ext d𝑆 = ∭ div𝑣 ⃗ d𝜏 = 0
∯ 𝑣 ⃗ ⋅ ⃗⃗⃗⃗
(𝑆) (𝑉 )

la vitesse est à flux conservatif.

Tube de courant – Toutes les lignes de courant qui s’appuient sur
une courbe 𝒞 fermée constituent un tube de courant. Dans ce cas, la
conservation du flux de vitesse s’exprime par Tube de
couran
• 𝑣
⃗⃗⃗⃗⃗⃗⃗⃗1 t
−∬ 𝑛ext d𝑆 = ∬
𝑣 ⃗ ⋅ ⃗⃗⃗⃗ 𝑛ext d𝑆
𝑣 ⃗ ⋅ ⃗⃗⃗⃗ soit 𝐷𝑉 1 = 𝐷𝑉 2
(𝑆1 ) (𝑆2 ) (𝑆1 ) • 𝑣
⃗⃗⃗⃗⃗⃗⃗⃗2

Le débit volumique se conserve le long d’un tube de courant.

(𝑆2 )
Si l’on définit la vitesse moyenne dans la section 𝑆 par
FIG. 3.9 : Tube de courant.
def 𝐷𝑉
⟨𝑣⟩ = ♡ (3.9)
𝑆

On obtient
⟨𝑣1 ⟩ 𝑆1 = ⟨𝑣2 ⟩ 𝑆2
Autrement dit, dans un tube de courant, le resserrement des lignes
de courant provoque une augmentation de la vitesse moyenne.
 DYNAMIQUE DES FLUIDES
PARFAITS 4
La dynamique des fluides relie l’écoulement d’un fluide aux actions 4.1 Équations dynamiques . 35
qui lui sont appliquées. Ce chapitre se limite aux écoulements pour Équation d’Euler . . . . . 35
lesquels les couches de fluide glissent les unes sur les autres sans Résolution de l’équation
dissipation de chaleur. On parle alors de fluide parfait. On détaille d’Euler . . . . . . . . . . . 35
particulièrement le cas des écoulements incompressibles et station- 4.2 Écoulements incompres-
naires. sibles stationnaires . . . 36
Théorème de Bernoulli . 36
Interprétation énergétique 37
Version en ligne Effet Venturi . . . . . . . 38
[Link]/mecanique_des_fluides/[Link] Formule de Torricelli . . . 39
Sonde de Pitot . . . . . . 40
4.3 Écoulements potentiels . 41
4.1 Équations dynamiques Vorticité . . . . . . . . . . 41
Écoulement potentiel . . 43

Équation d’Euler

Supposons un fluide en écoulement dans un champ de forces exté-

rieures que l’on sait exprimer. On admettra que :
1 : Si le référentiel n’est pas galiléen, il
• le fluide est en écoulement dans un référentiel galiléen1 ; suffit d’ajouter dans le bilan des forces,
• le fluide est parfait : les forces internes se résument aux forces les forces d’inertie.
de pression2
Appliquons le Principe fondamentale de la dynamique à une parti- 2 : Nous verrons dans le chapitre sui-
cule de fluide située en M à l’instant 𝑡, et de masse d𝑚 = 𝜌(M, 𝑡) d𝜏 . vant qu’il y a en fait des forces de fric-
tion, dites forces de viscosité, que l’on
Effectuons un bilan des forces. néglige ici.
ext
• Forces extérieures : ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 = 𝑓⃗⃗⃗⃗ext d𝜏 avec 𝑓 ext en N ⋅ m−3 .
int
• Forces pressantes internes : ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 = −⃗⃗⃗⃗⃗⃗ ∇𝑝(M, 𝑡) d𝜏
ext int
La 2nde loi de Newton d𝑚 D⃗⃗
𝑣⃗⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ + ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
D𝑡 = d𝐹 d𝐹 donne

𝜕𝑣⃗
𝜌(M, 𝑡) [ + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗ ∇𝑝(M, 𝑡) + 𝑓⃗⃗⃗⃗ext
∇) 𝑣]⃗ = −⃗⃗⃗⃗⃗⃗ ♡ (4.1)
𝜕𝑡

Cette relation est appelée équation d’Euler.

Résolution de l’équation d’Euler

La relation (4.1) est une équation aux dérivées partielles du premier

ordre. On remarque qu’elle est non linéaire à cause de la présence du
terme convectif (𝑣.⃗ ⃗⃗⃗⃗⃗⃗
∇)𝑣 ⃗ ; c’est ce qui rend les problèmes de mécanique
des fluides mathématiquement redoutables...
Regardons si nous disposons d’assez d’équations pour traiter un pro-
blème de mécanique des fluides parfaits.
Le fluide est incompressible – Dans ce cas la masse volumique est
fixée. Le problème présente donc 4 inconnues scalaires : le champ de
 36 4 DYNAMIQUE DES FLUIDES PARFAITS

pression 𝑝(M, 𝑡) et le champ de vitesse 𝑣(M,

⃗ 𝑡) (3 composantes). Il faut
donc 4 équations scalaires ! L’équation d’Euler n’en donne que 3. La
quatrième est donnée par l’équation de continuité div𝑣 ⃗ = 0

Le fluide est compressible – La masse volumique peut varier sous

l’effet de la pression mais aussi sous l’effet de la chaleur. En géné-
ral le fluide possède une équation d’état locale 𝜌(𝑝, 𝑇 ). Le problème
présente donc 6 inconnues scalaires : le champ de pression 𝑝(M, 𝑡),
les trois composantes du champ de vitesse 𝑣(M, ⃗ 𝑡), la masse volu-
mique 𝜌(M, 𝑡) et la température 𝑇 (M, 𝑡). Il faut donc 6 équations sca-
laires. L’équation d’Euler en donne 3, la quatrième est donnée par
l’équation de continuité div(𝜌𝑣)⃗ + 𝜕𝜌 𝜕𝑡 = 0, la cinquième par l’équa-
tion d’état du fluide 𝜌(𝑝, 𝑇 ) et la dernière par le premier principe de
la thermodynamique. Dans ce cas une bonne modélisation des trans-
ferts thermiques est nécessaire ce qui rend le problème très ardu. Par
exemple, étudier une étoile ou la combustion d’une flamme nécessite
ces 6 équations et surtout de gros ordinateurs...

L’équation d’Euler et de continuité sont des équations différentielles

du premier ordre ; leur intégration va donc produire des constantes
d’intégration. On déterminera ces constantes d’intégration par les
conditions aux interfaces (fluide1/fluide2 ou fluide/solide).

Conditions aux limites d’un fluide parfait

• Condition sur 𝑣(M,

⃗ 𝑡) : à la traversée d’une interface, la com-
3 : Condition valable si les effets
de tension superficielle sont négli- posante normale de la vitesse est continue.
geables. Si tel n’est pas le cas, il faut • Condition sur 𝑝(M, 𝑡) : la pression est continue à la traversée
utiliser la relation de Laplace qui relie d’une interface fluide-fluide3
la variation de pression entre deux in-
terfaces en fonction de la courbure de
l’interface et de la tension interfaciale.

4.2 Écoulements incompressibles

stationnaires

Théorème de Bernoulli

Ce théorème fut énoncé en premier par Daniel Bernoulli, en 1738.

Hypothèses – Le théorème de Bernoulli dans sa formulation clas-

sique ne s’applique qu’aux écoulements stationnaires incompres-
sibles et sans viscosité. Détaillons :

• Le fluide est parfait ; il obéit donc à l’équation d’Euler.

𝜕𝑣⃗
• L’écoulement est stationnaire. = 0⃗ et les trajectoires s’iden-
𝜕𝑡
tifient aux lignes de courant.
• L’écoulement est également incompressible. Le long d’une tra-
FIG. 4.1 : Daniel Bernoulli (1700 - 1782) : jectoire 𝜌 reste constant.
médecin, physicien et mathématicien • De plus, nous supposerons que les forces volumiques extérieures
suisse. Il publia en 1738 son ouvrage se résument à la pesanteur. Si l’on munit l’espace d’un axe O𝑧
Hydrodynamica dans lequel il expose
le théorème fondamental de la méca- vertical ascendant, on a
nique des fluides qui porte son nom.
𝑢𝑧 = −𝜌⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑓⃗⃗⃗⃗ext = 𝜌 𝑔 ⃗ = −𝜌𝑔 ⃗⃗⃗⃗⃗⃗⃗⃗ grad(𝑔 𝑧)
 4.2 Écoulements incompressibles stationnaires 37

L’équation d’Euler devient donc

𝜌 (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗ = −⃗⃗⃗⃗⃗⃗
∇𝑝 − 𝜌⃗⃗⃗⃗⃗⃗
∇(𝑔 𝑧)

2
𝑣
L’identité4 (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗ = ⃗⃗⃗⃗⃗⃗
∇ + (⃗⃗⃗⃗⃗⃗∇ ∧ 𝑣)⃗ ∧ 𝑣 ⃗ permet de réécrire l’équa- 4 : cf. Annexe B.
2
tion :
𝑣2
𝜌 (⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣)⃗ ∧ 𝑣 ⃗ = −⃗⃗⃗⃗⃗⃗
∇𝑝 − 𝜌⃗⃗⃗⃗⃗⃗
∇(𝑔𝑧) − 𝜌⃗⃗⃗⃗⃗⃗
∇
2
Intégrons cette équation le long d’une ligne de courant entre deux
points A et B : Rappelons que le long d’une trajec-
toire (qui se confond ici avec une
B 𝐵 2 ligne de courant), la masse volu-
𝑣
∫ 𝜌⏟ ⃗⃗⃗⃗⃗⃗
[(⏟
∇⏟⏟ ∧ 𝑣)⃗ ⏟
∧ 𝑣]⃗ ⋅⃗⃗⃗⃗⃗⃗⃗⃗
⏟⏟ ∇ (𝑝 + 𝜌𝑔𝑧 + 𝜌 ) ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗
dℓ = − ∫ ⃗⃗⃗⃗⃗⃗ dℓ mique conserve une valeur constante
A 𝐴 2 puisque l’écoulement est incompres-
⟂⃗⃗𝑣⃗⃗
sible.

Le premier membre est nul puisque ⃗⃗⃗⃗⃗⃗⃗⃗

dℓ ∥ 𝑣.⃗ Le second membre s’écrit
𝐵
𝑣2 𝑣2 𝑣2
− ∫ d (𝑝 + 𝜌𝑔𝑧 + 𝜌 ) = (𝑝A + 𝜌𝑔𝑧A + 𝜌 A )−(𝑝B + 𝜌𝑔𝑧B + 𝜌 B )
𝐴 2 2 2

In fine, on aboutit à la relation

𝑣A2 𝑣2
𝑝A + 𝜌𝑔𝑧A + 𝜌 = 𝑝B + 𝜌𝑔𝑧B + 𝜌 B
2 2

Théorème de Bernoulli (1738)

Pour un écoulement incompressible et permanent d’un fluide par-

fait dans le champ de pesanteur, la quantité

𝑣2
𝑝+𝜌 + 𝜌𝑔𝑧 = Cte le long d’une ligne de courant (4.2) Attention, l’axe O𝑧 doit être ascen-
2
dant !
NB : d’une ligne de courant à l’autre, c’est la valeur de la constante
qui change.

Interprétation énergétique
2
La conservation de la quantité 𝑝 + 𝜌 𝑣2 + 𝜌𝑔𝑧 exprime la conservation
de l’énergie le long d’une ligne de courant.
2
• 𝜌 𝑣2 représente l’énergie cinétique volumique ;
• 𝜌𝑔𝑧 l’ énergie potentielle volumique de pesanteur ;
• et la pression 𝑝 représente l’énergie potentielle volumique as-
sociée aux forces de pression.

Démonstration
Considérons un fluide incompressible en écoulement permanent plongé
dans un champ de pesanteur. Considérons, un filet de fluide de section in-
finitésimale, et appliquons le théorème de l’énergie cinétique entre deux
38 4 DYNAMIQUE DES FLUIDES PARFAITS

instants voisins 𝑡 et 𝑡 + d𝑡 :

dℰc = 𝛿𝑊

d𝑆2
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹2 𝑣2
⃗⃗⃗⃗⃗⃗⃗⃗
𝑝2
𝑣2 d𝑡
𝑣1 d𝑡

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑣1
⃗⃗⃗⃗⃗⃗⃗⃗ 𝑧2
𝐹1

𝑝1
d𝑆1 𝑧1

À un instant 𝑡 le système étudié est une portion du tube de courant com-

pris entre les sections d𝑆1 (𝑡) et d𝑆2 (𝑡). À l’instant 𝑡 + d𝑡 la section d𝑆1 (𝑡)
s’est déplacée de 𝑣1 d𝑡 et la section d𝑆2 (𝑡) de 𝑣2 d𝑡. Chaque surface balaye
le même volume d𝜏 car le fluide est incompressible :

d𝜏 = 𝑣1 d𝑆1 d𝑡 = 𝑣2 d𝑆2 d𝑡

Commençons par exprimer l’énergie cinétique du système :

1
ℰc (𝑡) = 𝜌d𝜏 𝑣1 2 + ℰ0c (𝑡)
2
où ℰ0c (𝑡) désigne l’énergie cinétique de la portion de fluide autre que le
volume d𝜏 . De même, à l’instant 𝑡 + d𝑡, on a
1
ℰc (𝑡 + d𝑡) = ℰ0c (𝑡 + d𝑡) + 𝜌d𝜏 𝑣2 2
2
Supposons l’écoulement stationnaire. Dans ce cas ℰ0c (𝑡) = ℰ0c (𝑡 + d𝑡) et
l’on aboutit à
1
dℰc = 𝜌d𝜏 (𝑣2 2 − 𝑣1 2 )
2
Quant au travail mécanique, il se résume à deux termes :

1. le travail de pesanteur 𝛿𝑊𝑔 = −dℰp = −𝜌d𝜏 𝑔(𝑧2 − 𝑧1 ).

2. le travail des forces de pression 𝛿𝑊𝑝 = −𝑝ext d𝑉 = (𝑝1 − 𝑝2 )d𝜏

Notez qu’il n’y a pas de travail des forces internes pour deux raisons :

1. au sein du système, les forces de pression ne peuvent ni dilater ni

comprimer les particules de fluide car l’écoulement est incompres-
sible ;
2. les forces de friction dues à la viscosité sont négligés ici puisque
l’on traite le cas d’un fluide parfait.

Le théorème de l’énergie cinétique donne finalement

1 1
𝜌𝑣 2 + 𝜌𝑔𝑧1 + 𝑝1 = 𝜌𝑣2 2 + 𝜌𝑔𝑧2 + 𝑝2
2 1 2
On obtient la relation de Bernoulli, montrant ainsi qu’elle exprime tout
simplement la conservation de l’énergie mécanique.

Effet Venturi

Dans un tube horizontal de section 𝑆 variable, l’écoulement d’un

fluide en écoulement incompressible et permanent s’accompagne
d’une dépression là où il y a rétrécissement : c’est l’effet Venturi.
 4.2 Écoulements incompressibles stationnaires 39

• 𝑣
⃗⃗⃗⃗⃗⃗⃗⃗A B • 𝑣B
⃗⃗⃗⃗⃗⃗⃗⃗
A 𝑆B

𝑆A
FIG. 4.2 : Effet venturi

Explications – Le long d’un tube horizontal, d’après la conservation

du débit on a
𝐷𝑉 = 𝑣A 𝑆A = 𝑣B 𝑆B
Ainsi, la conservation du débit impose une augmentation de vitesse
là ou il y a rétrécissement. Le théorème de Bernoulli impose

1 1
𝜌 𝑣 2 + 𝑝A +  

𝜌𝑔𝑧 2
A = 𝜌 𝑣 B + 𝑝B + 


𝜌𝑔𝑧
2 A 2 B

Il existe alors une dépression au niveau du rétrécissement donnée

par
1
𝑝B − 𝑝A = 𝜌 (𝑣A 2 − 𝑣B 2 ) < 0
2
Cet effet peut être mis à profit pour les applications suivantes :

• Mesure de débit (débitmètre à effet venturi) ;

• Principe des trompes à eau montées sur les robinet des paillasses
de chimie ;
• Douchette venturi produisant une économie d’eau ;
• Amélioration du tirage d’une cheminée, principe du carburateur,
vaporisateur, etc.

Exercice – Montrer que la hauteur manométrique ℎ de la FIG. 4.2 s’écrit

2
ℎ = 𝐾𝐷𝑉 où 𝐾 est un facteur géométrique que l’on exprimera.

Formule de Torricelli

Considérons un réservoir cylindrique rempli d’un liquide dans lequel

on perce un orifice. La formule de Torricelli relie le débit d’écoulement
section 𝑆
avec la hauteur de liquide ℎ. On fera les hypothèses suivantes : A•
𝑔⃗⃗⃗⃗

• La section 𝑆 du cylindre est très grande devant la section de

l’orifice : 𝑠 ≪ 𝑆 ;
• On considère le liquide incompressible et parfait ;
ℎ(𝑡)
• Enfin, on considère que l’écoulement est en régime stationnaire.

section 𝑠
On cherche à calculer la vitesse d’écoulement 𝑣(B, 𝑡) à la sortie du
trou. L’application du théorème de Bernoulli sur une ligne de courant
donne : •
B 𝑣(B,
⃗⃗⃗⃗ 𝑡)
1 1
𝑝atm + 𝜌𝑔ℎ + 𝜌𝑣2 (A, 𝑡) = 𝑝atm + 𝜌𝑣2 (B, 𝑡) FIG. 4.3 : Vidange d’un récipient.
2 2
 40 4 DYNAMIQUE DES FLUIDES PARFAITS

Or, la conservation du débit volumique donne

𝑣(A, 𝑡) 𝑆 = 𝑣(B, 𝑡) 𝑠

d’où 𝑣(A, 𝑡) ≪ 𝑣(B, 𝑡) car 𝑠 ≪ 𝑆. Finalement

On peut noter que la vitesse du jet sor- 𝑣(B, 𝑡) = √2𝑔ℎ(𝑡) [formule de Torricelli] ♡ (4.3)
tant dépend du temps ce qui montre
que l’écoulement n’est pas rigoureuse-
ment stationnaire. On peut cependant
La vitesse a la même expression que celle de la chute libre d’un point
montrer que 𝜕 𝑣/𝜕𝑡
⃗⃗⃗⃗ est négligeable si matériel dans le champ de pesanteur.
𝑠 ≪ 𝑆. On parle d’écoulement quasi-
stationnaire. Cherchons à déterminer le temps de vidange 𝑇 . Pour cela, reprenons
la conservation du débit volumique :
𝑠
𝐷𝑉 = 𝑠√2𝑔ℎ(𝑡) = 𝑆𝑣(A, 𝑡) d’où 𝑣(A, 𝑡) = √2𝑔ℎ(𝑡)
𝑆
Par ailleurs, le point A situé sur la surface libre a pour côte 𝑧 = ℎ(𝑡)
𝑣𝐴 = dℎ
d’où une vitesse ⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑢𝑧 . On a donc
d𝑡 ⃗⃗⃗⃗⃗⃗⃗⃗
Attention aux signes. La vitesse en A a
pour norme 𝑣𝐴 = |ℎ|̇ = −ℎ,̇ car ℎ(𝑡) dℎ 𝑠
décroît.
𝑣(A, 𝑡) = − = √2𝑔ℎ
d𝑡 𝑆
Il s’agit d’une équation différentielle que l’on peut résoudre en sépa-
rant les variables :
0 √ 𝑇
dℎ 𝑠 2𝑔
∫ √ =− ∫ d𝑡
ℎ0 ℎ 𝑆 0

où ℎ0 désigne la hauteur initiale du liquide. On obtient le résultat

suivant :
𝑆 2ℎ
𝑇 = √ 0
𝑠 𝑔

En pratique, le jet de sortie est contracté de sorte que la section effective

de sortie 𝑠 est légèrement plus petite que la section de l’orifice.

Sonde de Pitot

La sonde de Pitot fut inventé en 1732 par l’ingénieur français Henri

5 : On parle souvent de tube de Pitot.
Pitot. Il s’agit d’un tube5 destiné à mesurer la vitesse d’écoulement
d’un fluide ; il est largement employé dans le génie chimique et en
aéronautique sous une version améliorée que l’on appelle sonde de
Prandtl.

La sonde est constituée d’un tube de diamètre 𝑑 au bout arrondi et

percée d’un trou cylindrique parallèle au tube. Sur la face latérale du
FIG. 4.4 : Sonde de Pitot statique dite tube à une distance de l’ordre de 3𝑑 se situe une ou plusieurs prises
aussi sonde de Prandtl. de pression. On place cette sonde parallèlement à un écoulement
de fluide et l’on mesure la différence de pression entre les pressions
axiale et latérale.

6 : Pour les gaz, l’approximation est Admettons que le fluide soit en écoulement stationnaire et incom-
justifiée si l’écoulement est subso- pressible6 . Le point A est un point d’arrêt car la vitesse est nulle (il n’
nique (𝑣 ≪ 𝑐son )
 4.3 Écoulements potentiels 41

B’
•
B∞ •
• B manomètre
• •
A∞ 𝑣
⃗⃗⃗⃗∞ A différentiel

Δ𝑝

FIG. 4.5 : Principe de la sonde de Pitot statique.

y a pas d’écoulement dans l’orifice , c’est juste une prise de pression).

Le théorème de Bernoulli appliqué entre A∞ et A donne

1 2
𝑝A = 𝑝∞ + 𝜌𝑣∞
2
Le même théorème appliqué entre B∞ et B’ donne

1 2 2
𝑝B’ = 𝑝∞ + 𝜌 (𝑣∞ − 𝑣B’ )
2
Par ailleurs, si l’on admet que l’écoulement est parallèle autour de la
À retenir
sonde, alors il en découle
𝑝B = 𝑝B’ Quand un fluide parfait ou vis-
queux présente un écoulement pa-
ce qui conduit à rallèle, la pression obéit aux lois
1 de la statique dans une direction
Δ𝑝 = 𝑝A − 𝑝B = 𝜌𝑣 2 perpendiculaire à l’écoulement.
2 B’
La différence de pression indiquée par le manomètre permet de re-
monte à la vitesse en B’. De surcroît si l’on suppose que la section
du tube de Pitot est suffisamment faible devant la section du tube
de courant pour ne pas influencer la vitesse d’écoulement alors on a
𝑣B’ = 𝑣∞ .

Exercice – Un tube de Pitot dans un écoulement d’air mesure une différence

de pression Δ𝑃 = 0,6 mbar. Quelle est la vitesse d’écoulement ?
Rép. 10 m/s.

4.3 Écoulements potentiels

Vorticité

On appelle vorticité le rotationnel de la vitesse. Il s’agit d’un vecteur

qui se note 𝜔
⃗⃗⃗ ⃗ :

𝜔(M,
⃗⃗⃗ ⃗ 𝑡) = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot 𝑣(M, ⃗ 𝑡) = ⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣(M,
⃗ 𝑡) ♡ (4.4)

À l’instar des lignes de courant, on définit les lignes de vorticité comme

les lignes du champ vectoriel 𝜔(M,
⃗⃗⃗ ⃗ 𝑡).
42 4 DYNAMIQUE DES FLUIDES PARFAITS

Exemple

Considérons l’écoulement bidimensionnel dont le champ de vitesse s’écrit

𝑣(M,
⃗ 𝑡) = −𝑘 𝑦 ⃗⃗⃗⃗⃗⃗⃗⃗ 𝑢𝑦 avec 𝑘 = Cte
𝑢𝑥 + 𝑘 𝑥 ⃗⃗⃗⃗⃗⃗⃗⃗

Tout d’abord, on constate que

𝜕(−𝑘 𝑦) 𝜕(𝑘 𝑥)
div 𝑣 ⃗ = + =0
𝜕𝑥 𝜕𝑦

Il s’agit donc d’un écoulement incompressible. La vorticité de cet écoule-

ment est donnée par

𝜕/𝜕𝑥 −𝑘𝑦
⎛ ⎞ ⎛ ⎞
⃗⃗⃗ ⃗ = ⃗⃗⃗⃗⃗⃗
𝜔 ∇ ∧ 𝑣⃗ = ⎜
⎜ 𝜕/𝜕𝑦 ⎟ ⎜ 𝑘𝑥 ⎟
⎟ ∧ ⎜ ⎟ = 2𝑘 ⃗⃗𝑢⃗ ⃗𝑧⃗ ⃗
⎝ 𝜕/𝜕𝑧 ⎠ ⎝ 0 ⎠

Cet écoulement présente donc une vorticité uniforme.

La vorticité a une signification cinématique simple : elle donne le

double du vecteur rotation local des particules de fluide. Ce vecteur,
noté ⃗⃗⃗⃗⃗
Ω, est parfois appelé vecteur tourbillon.

⃗⃗⃗⃗⃗ 1
Ω = ⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣⃗
2

Exemple

Reprenons l’exemple précédent et plaçons une particule carrée d’arête

2𝑎 en O. Quel est son mouvement ?

On a vu qu’au cours de l’écoulement, cette particule conserve son aire

puisque div𝑣 ⃗ = 0.⃗ Sur la simulation ci-dessus on constate que les som-
mets sont animés d’une vitesse orthoradiale de sorte que √ la particule
tourne autour de O. La vitesse d’un sommet vaut 𝑣 = 2𝑘, 𝑎 d’où une
vitesse angulaire de rotation Ω = 𝑘 = 21 ‖⃗⃗𝜔‖.
⃗⃗

Si l’on prend le rotationnel de l’équation d’Euler, on obtient l’équa-

tion d’évolution de la vorticité. Dans le cas d’un fluide incompressible
siège d’un champ de force extérieur conservatif, on aboutit à l’équa-
tion
𝜕⃗⃗𝜔⃗ ⃗ ⃗⃗⃗⃗⃗⃗
= ∇ ∧ (𝑣 ⃗ ∧ 𝜔)
⃗⃗⃗ ⃗
𝜕𝑡
En conséquence, si un fluide parfait est irrotationnel à un instant
donné (⃗⃗𝜔(M,
⃗⃗ 𝑡0 , ) = 0⃗ partout), alors il le reste aux instants ultérieurs
(𝜕⃗⃗𝜔/𝜕𝑡
⃗⃗ ⃗
= 0 pour 𝑡 > 𝑡0 ).
 4.3 Écoulements potentiels 43

Écoulement potentiel

Le phénomène de viscosité – absent ici – est source de vorticité no-

tamment dans le sillage des obstacles et dans une couche, dite couche
limite. Mais en général, en dehors ce ces zones le fluide peut être
traité comme un fluide parfait. Et si en plus l’écoulement est initiale-
ment irrotationnel, alors il le restera ultérieurement. C’est pourquoi
ce type d’écoulement revêt une certaine importance en dynamique
des fluides.
Lorsque 𝜔⃗⃗⃗ ⃗ = 0,⃗ le champ de vitesse s’écrit nécessairement comme un
gradient. On dit qu’il dérive d’un potentiel des vitesses 𝜑(M, 𝑡) :

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot 𝑣(M, ⃗ 𝑡) = 0⃗ d’où 𝑣(M,
⃗ 𝑡) = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝜑(M, 𝑡)

C’est pourquoi, on parle dans ce cas d’écoulement irrotationnel ou

potentiel. Dans ces conditions, l’équation d’Euler devient

𝜕 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝜑 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑣2
𝜌(M, 𝑡) [ ∇𝑝(M, 𝑡) + 𝑓⃗⃗⃗⃗ext
+ grad ( )] = −⃗⃗⃗⃗⃗⃗
𝜕𝑡 2

Si l’on suppose en plus le fluide incompressible dans un champ de

force conservatif d’énergie potentielle volumique 𝑒ext
p , on aboutit à

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝜕𝜑 1 2 ⃗
grad [𝜌 + 𝜌𝑣 + 𝑒ext
p + 𝑝] = 0
𝜕𝑡 2

soit, après intégration :

𝜕𝜑 1 2
𝜌 + 𝜌𝑣 + 𝑒ext
p + 𝑝 = 𝑓(𝑡) ♡ (4.5)
𝜕𝑡 2

où 𝑓(𝑡) est une fonction arbitraire du temps. Notons que le potentiel

des vitesses est indéterminée : si 𝜑(M, 𝑡) convient, alors 𝜑(M, 𝑡) + 𝑔(𝑡)
convient aussi, de sorte que l’on peut toujours choisir le potentiel
des vitesses de sorte que 𝑓(𝑡) = 0.

FIG. 4.6 : Écoulement potentiel d’un

fluide à la vitesse d’un millimètre par
seconde autour d’un obstacle cylin-
drique. L’écoulement est permanent,
incompressible et irrotationnel.

Bien entendu, dans le cas d’un écoulement stationnaire, 𝜑 ne dépend

plus explicitement du temps, et l’on retrouve la relation de Bernoul-
li :
1 2
𝜌𝑣 + 𝑒ext
p +𝑝 =C
te
2
à une nuance près : la constante est la même pour toutes les lignes
d’écoulement.
 FLUIDES VISQUEUX 5
Nous avons vu au chapitre précédent que dans un fluide parfait, la 5.1 Notion de viscosité . . . . 45
contrainte qui s’exerce sur une particule de fluide est toujours per- Expérience de Couette . . 45
pendiculaire aux parois de celle-ci. Dans un fluide réel en écoule- Fluide newtonien . . . . . 46
ment, la contrainte possède aussi une composante tangentielle dite Rhéologie . . . . . . . . . 47
contrainte visqueuse. Sa prise en compte permet d’établir l’équation 5.2 Dynamique des écoule-
de Navier-Stokes qui représente une version améliorée de l’équation ments visqueux incompres-
d’Euler pour un fluide réel. sibles . . . . . . . . . . . 47
Force volumique visqueuse 48
Équation de Navier-Stokes 48
Version en ligne Le nombre de Reynolds . 49
[Link]/mecanique_des_fluides/[Link]

5.1 Notion de viscosité

Expérience de Couette

Considérons un fluide enfermé entre deux cylindres, l’un mobile, l’autre

−1
fixé via un fil de torsion de constante de torsion 𝐶 (en N ⋅ m ⋅ rad ).
L’ensemble forme ce qu’on appelle un rhéomètre. On constate que
lorsque la cavité cylindrique extérieure est mise en rotation à la vi-
tesse angulaire 𝜔, le cylindre intérieur tourne d’un angle 𝜃 par rapport
à sa position d’équilibre.

Fil de torsion

Liquide

𝐿 équilibre pour 𝜔 = 0
•
𝑒 0
𝜔 >
𝜃
o ur
𝑅1 p
re
𝑅2 u ilib
éq
𝜔 FIG. 5.1 : Rhéomètre de Couette. À
𝜔
gauche, vue de profil. À droite, vue de
dessus.

Ce que l’on observe :

• Les particules de fluide adhèrent aux parois. Il existe donc un
gradient de vitesse au sein de l’écoulement.
• Plus la vitesse angulaire 𝜔 est importante, plus l’angle de tor-
sion 𝜃 l’est aussi.
• Pour les fluides simples, cette relation est linéaire.

Conséquences :
• La torsion du fil conduit à l’existence d’un couple produit par le
fluide. Ce couple, noté ℳΔ équilibre le couple de torsion −𝐶𝜃
de sorte que ℳΔ = 𝐶𝜃.
 46 5 FLUIDES VISQUEUX

• Les forces de pression ne peuvent pas expliquer ces actions.

On est donc obligé d’admettre l’existence d’efforts tangentiels.
On définit la contrainte tangentielle ⃗⃗𝜎⃗ ⃗t⃗ qui s’applique sur un
élément de surface d𝑆 par

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹t = ⃗⃗𝜎⃗ ⃗t⃗ × d𝑆
♡ (5.1)
N = Pa × m2

• Ces forces dépendent du gradient de vitesse qui règne au sein

du rhéomètre. Ce gradient de vitesse, aussi appelé vitesse de
cisaillement, est noté 𝛾.̇
Le rhéomètre permet de mesurer la contrainte 𝜎t et la vitesse de ci-
saillement 𝛾̇ grâce à la mesure de l’angle de torsion 𝜃 et la vitesse
angulaire 𝜔. En effet, on a

ℳΔ /𝑅1 𝐶𝜃 Δ𝑣 𝑅2 𝜔
𝜎t = = 2
et 𝛾̇ = =
𝑆lat 2𝜋𝐿𝑅1 Δ𝑅 𝑅2 − 𝑅1

La courbe 𝜎t = 𝑓(𝛾)̇ est appelée rhéogramme et décrit le comporte-

ment dynamique du liquide en écoulement.

Fluide newtonien
Définition
Un fluide est newtonien quand la contrainte visqueuse est propor-
tionnelle à la vitesse de cisaillement :

𝜎t = 𝜂 𝛾 ̇ (5.2)

avec 𝜂 une constante appelée viscosité dynamique.

L’analyse dimensionnelle de la relation (5.2) donne

1 : Jean-Louis Marie Poiseuille (1797–
1869) fut élève de l’École Polytech-
nique avant d’étudier la médecine.
[𝐹 ] L [𝐹 ]
[𝜂] = = 2 T
Les recherches de Poiseuille ont porté L2 [𝑣] L
principalement sur l’hémodynamique,
c’est-à-dire l’étude de la circulation Ainsi, la viscosité est homogène à une pression × temps. On l’ex-
sanguine et lui ont permis de dégager prime indifféremment en [Link] (Pa.s) ou en poiseuille (Pℓ)
une loi sur l’écoulement des fluides
visqueux dans des tubes capillaires.
en hommage à Jean-Louis Marie Poiseuille1 .
Ordre de grandeur – Tous les gaz sont des fluides newtoniens avec
une viscosité de l’ordre de 1×10−5 Pa⋅s dans les conditions normales.
Liquide (20 °C) Viscosité (Pa.s) Elle varie peu avec la température. Quant aux liquides, la viscosité
varie sur plusieurs ordres de grandeurs (TAB. 5.1) et reste très sensible
Acétone 3,2 × 10−4
Eau 1,0 × 10−3 à la température (elle diminue lorsque la température augmente).
Mercure 1,5 × 10−3
Huile d’olive 0,84
Pour résumer, au sein d’un fluide newtonien, il existe deux types d’ac-
Glycérine pure 1,5 tions de contact :
Gaz Viscosité (Pa.s) • les actions pressantes dont les contraintes sont normales :
Vapeur d’eau (20 °C) 9,7 × 10−6
Air sec (20 °C) 18,2 × 10−6 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 n
He (25 °C) 19,9 × 10−6 ⃗⃗⃗⃗⃗⃗⃗⃗
𝜎n = = −𝑝 ⃗⃗⃗⃗
𝑛
d𝑆
N2 (25 °C) 17,7 × 10−6

TAB. 5.1 : Quelques valeurs de viscosité. Rappelons que ses forces sont conservatives.
 5.2 Dynamique des écoulements visqueux incompressibles 47

𝑣
⃗⃗⃗⃗ ⃗⃗⃗⃗⃗
𝑛
• les actions visqueuses dont les contraintes sont tangentielles 𝑡⃗
⃗⃗𝜎⃗ ⃗ t⃗⃗
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 t 𝜕𝑣
⃗⃗𝜎⃗ ⃗t⃗ = =𝜂 𝑡⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝜎n
d𝑆 𝜕𝑛
où 𝜕𝑣/𝜕𝑛 désigne le gradient de vitesse dans la direction nor-
FIG. 5.2 : Contraintes internes dans le
male à la surface. Ces forces sont dissipatives. fluide.

Rhéologie

La rhéologie est la science qui étudie le comportement des fluides en

écoulement. Cette science repose sur l’emploi de rhéomètres2 , appa- 2 : On dit aussi viscosimètres.
reils destinés à mesurer la relation entre la contrainte tangentielle 𝜎t
et la vitesse de cisaillement 𝛾.̇ Le comportement newtonien s’observe
dans tous les gaz et dans les liquides simples constitués de petites
molécules (l’eau par exemple).

Toutefois, certains liquides ont des comportements rhéologiques qui

s’écartent du comportement newtonien. La FIG. 5.3 en donnent quelques
exemples.
contrainte

viscosité

Fluide Newtonien 𝜂 = Cte

Exemples : Eau, huiles, Gaz…

vitesse de cisaillement vitesse de cisaillement

Fluide Pseudoplastique
contrainte

viscosité

Exemples : Peintures, vernis,

produits cosmétiques, sang,
suspensions minérales…

vitesse de cisaillement vitesse de cisaillement

Fluide plastique (à seuil)

contrainte

viscosité

Exemples : pate dentifrice, cho-

colat, moutarde, graisse, sus- 𝜎0
pensions d’argile…

vitesse de cisaillement vitesse de cisaillement

FIG. 5.3 : Différents comportements rhéologiques.

5.2 Dynamique des écoulements visqueux

incompressibles

Lorsque le fluide est newtonien et incompressible, les équations de

Newton appliquées à chaque particule de fluide prennent la forme
des équations de Navier-Stokes.
 48 5 FLUIDES VISQUEUX

Force volumique visqueuse

Plaçons nous dans un référentiel galiléen et effectuons un bilan des

forces sur une particule de fluide située en M à l’instant 𝑡, de masse
d𝑚 = 𝜌(M, 𝑡) d𝜏 . En plus des forces de pression et des forces exté-
rieures volumiques, il faut ajouter la résultante des forces visqueuses :

⃗⃗⃗ ⃗⃗ = (𝑓⃗⃗⃗⃗ext − ⃗⃗⃗⃗⃗⃗

d𝐹 ∇𝑝) d𝜏 + d⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝜂

L’expression de d𝐹 ⃗⃗⃗ ⃗𝜂⃗ est en général assez compliquée mais elle se

simplifie dans le cas des fluides newtoniens et incompressibles.

𝑦 𝑣
⃗⃗⃗⃗ ⃗⃗⃗⃗⃗
𝑛
Cas d’un écoulement parallèle unidimensionnel – Calculons la résul-
𝑡⃗
⃗⃗⃗ ⃗ t⃗⃗(𝑦 + d𝑦) tante des forces visqueuses dans le cas particulier simple d’un écou-
𝜎
lement suivant (𝑂𝑥) avec un gradient de vitesse suivant (𝑂𝑦) :
𝑡⃗
⃗⃗𝜎⃗ ⃗ t⃗⃗(𝑦)
⃗⃗⃗⃗⃗
𝑛 𝑣 ⃗ = 𝑣(𝑦) ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥

On remarque ici que div𝑣 ⃗ = 0. L’écoulement est donc bien incompres-

FIG. 5.4 : Bilan des forces de viscosité
sur un élément de fluide. sible. Dans ce cas, la résultante des forces visqueuses s’exerçant sur
une particule de fluide, s’écrit :
2
d𝑣 d𝑣 d 𝑣⃗
d⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝜂 = 𝜂 [ (𝑦 + d𝑦) − (𝑦)] d𝑥d𝑧 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥 = 𝜂 2 d𝜏
dy dy dy

On voit apparaître une force volumique qui s’exprime comme le la-

placien de la vitesse. Cette formule obtenue dans un cas particulier
se généralise aux écoulements incompressibles des fluides newto-
niens. On admettra que pour un fluide newtonien incompressible, la
résultante des forces visqueuses s’écrit

𝜕2𝑣 ⃗ 𝜕2𝑣 ⃗ 𝜕2𝑣 ⃗

d⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝜂 = 𝜂Δ𝑣 ⃗ d𝜏 = 𝜂 ( 2 + 2 + 2 ) d𝜏 ♡ (5.3)
𝜕𝑥 𝜕𝑦 𝜕𝑧

où Δ est l’opérateur laplacien.

Équation de Navier-Stokes

Appliquons la seconde loi de Newton à une particule de fluide :

D𝑣 ⃗ ⃗⃗⃗⃗ext d𝜏
𝜌d𝜏 − ⃗⃗⃗⃗⃗⃗
= ⏟ ∇𝑝 d𝜏 + 𝜂Δ
⏟ 𝑣 ⃗ d𝜏 + 𝑓⏟
D𝑡
pression viscosité autres

En divisant par d𝜏 , on obtient l’équation de Navier-Stokes.

Équation de Navier-Stokes

Pour un fluide incompressible newtonien, la dynamique de l’écou-

lement vérifie l’équation

𝜕𝑣⃗
𝜌[ + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗ ∇𝑝 + 𝑓⃗⃗⃗⃗ext + 𝜂Δ𝑣 ⃗
∇)𝑣]⃗ = −⃗⃗⃗⃗⃗⃗ (5.4)
𝜕𝑡
 5.2 Dynamique des écoulements visqueux incompressibles 49

Il s’agit donc d’une équation aux dérivées partielles du second ordre

et non linéaire. Cette équation recèle encore quelques mystères qui
résistent à la sagacité de nos meilleurs mathématiciens puisque l’existence
et l’unicité d’une solution de l’équation de Navier-Stokes est l’un des
7 problèmes du millénaire mis à prix $ 1 000 000 par l’Institut Clay3 ! 3 : Voir [Link]
org/millennium-problems/
Conditions aux limites – L’équation de Navier-Stokes étant une équa- navier-stokes-equation
tion aux dérivées partielles, sa résolution introduit des constantes
d’intégration. On les détermine en appliquant les conditions aux li-
mites suivantes :

• continuité de la vitesse à la traversée d’une interface (adhé-

rence) ;
• continuité de la pression si la surface est plane ;
• continuité de la contrainte tangentielle.

Le nombre de Reynolds

La complexité provient essentiellement de la présence, dans l’équa-

tion de Navier-Stokes, d’un terme non linéaire – le terme d’inertie –
et d’un terme du second ordre – le terme de viscosité.
𝜕𝑣⃗
+ 𝜌 (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗ = −⃗⃗⃗⃗⃗⃗ ⏟𝑣 ⃗ +𝑓⃗⃗⃗⃗
ext
𝜌 ∇𝑝 + 𝜂Δ
𝜕𝑡 ⏟⏟⏟⏟⏟
inertie viscosité

Dans de nombreux cas, on peut négliger l’un des deux termes devant
l’autre. On définit alors un facteur sans dimension, qui estime l’impor-
tance du terme d’inertie devant le terme de viscosité. On peut estimer
l’ordre de grandeur de ces termes à partir de la masse volumique 𝜌
du fluide, de sa viscosité 𝜂, d’une vitesse caractéristique de l’écoule-
ment 𝑣 et d’une distance 𝑑 qui caractérise l’échelle de variation de la
vitesse. En effet,

𝑣2 𝑣
∥𝜌 (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣∥⃗ ∼ 𝜌 et ‖𝜂Δ𝑣‖⃗ ∼ 𝜂
𝑑 𝑑2
D’où le nombre sans dimension appelé nombre de Reynolds

terme d’inertie 𝜌𝑣𝑑

𝑅e = = ♡ (5.5)
terme visqueux 𝜂

Ce nombre joue un rôle très important en mécanique des fluides car

il permet de distinguer trois types d’écoulement.

• L’écoulement à petit nombre de Reynolds 𝑅e ≪ 1 – L’écoule-

ment est laminaire et essentiellement gouverné par la viscosité.
Le terme d’inertie est négligeable et l’équation de Navier-Stokes
devient
𝜕𝑣⃗
𝜌 ∇𝑝 + 𝑓⃗⃗⃗⃗ext + 𝜂Δ𝑣 ⃗
= −⃗⃗⃗⃗⃗⃗
𝜕𝑡
équation qui a le bon goût d’être linéaire. Si l’écoulement est
permanent, on obtient le régime de Stokes.
• L’écoulement à grand nombre de Reynolds 𝑅e ≫ 1 – On montre
dans ce cas que les effets visqueux sont concentrés sur les
 50 5 FLUIDES VISQUEUX

bords, dans une fine couche appelée couche limite, et dans le

sillage des obstacles. Hors de ces zones, le terme visqueux est
négligeable et l’on retrouve l’équation d’Euler

𝜕𝑣⃗
𝜌[ + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗ ∇𝑝 + 𝑓⃗⃗⃗⃗ext
∇) 𝑣]⃗ = −⃗⃗⃗⃗⃗⃗
𝜕𝑡

• Écoulement turbulent – La viscosité stabilise et régularise les

écoulements de façon générale. Ainsi, quand le nombre de Rey-
nolds augmente le régime laminaire devient instable voire tur-
bulent. La transition entre le régime laminaire et turbulent se
produit dans une certaine gamme de valeur du nombre de Rey-
nolds qui dépend du problème. On peut retenir qu’en général,
lorsque 𝑅e ≳ 103 − 104 , l’écoulement devient turbulent : la vi-
tesse en un point M varie dans le temps de façon erratique. Dans
ce cas, le problème étant analytiquement insoluble, on utilise
souvent des lois phénoménologiques associées à une analyse
dimensionnelle.
TAB. 5.2 : Quelques ordres de grandeur
Type d’écoulement Nombre de Reynolds Régime
du nombre de Reynolds.
Écoulement atmosphérique 𝑅e ≈ 1011 turbulent
Écoulement sanguin dans l’aorte 𝑅e ≈ 104 turbulent
Écoulement sanguin dans les capillaires 𝑅e ≈ 10−3 laminaire visqueux
Domaine de l’aéronautique 𝑅e ≈ 107 turbulent
Domaine de la microfluidique 𝑅e ≈ 10−3 − 1 laminaire visqueux

Exercice – De l’eau (𝜂 = 1 × 10−3 Pa ⋅ s) à 20 °C et 1 × 105 Pa coule dans une

−1
canalisation de diamètre 𝑑 = 1 cm avec un débit 𝐷𝑉 = 10 L ⋅ min . Calculer
le nombre de Reynolds.
 PERTES DE CHARGE ET FRICTION 6
Nous abordons dans ce chapitre deux conséquences aux caractère 6.1 Pertes de charge . . . . . 51
visqueux de la plupart des fluides : l’existence de pertes de charge Observations expérimen-
dans les canalisations et, la résistance au mouvement d’un corps tales . . . . . . . . . . . . 51
dans un fluide. Écoulement de Poiseuille 52
Loi de Poiseuille . . . . . 53
Analogie électrique . . . 54
Version en ligne
Formule générale . . . . . 54
[Link]/mecanique_des_fluides/[Link] Pertes de charge singulière 56
Th. de Bernoulli généralisé 57
6.2 Résistance au mouvement 59
6.1 Pertes de charge dans une conduite Formule de Stokes . . . . 59
Coefficient de traînée Cx . 60

Observations expérimentales
Expérience

Considérons un récipient de diamètre 𝐷 dans lequel on branche un tube

horizontal de diamètre 𝑑 ≪ 𝐷. Ce tube dispose de 3 tubes verticaux en
contact avec l’air.

Dans un premier temps, on bouche la sortie du tube horizontal et on

remplit le récipient d’un fluide. On observe ceci :

Diamètre 𝐷

Liquide
visqueux
ℎ0
𝑔⃗⃗⃗⃗
𝐿

Dans un deuxième temps, on ouvre le tube horizontal : le niveau d’eau

descend lentement, mais la hauteur des colonnes d’eau diffère d’un tube
à l’autre.

A
•

ℎ(𝑡)

B
•

Analyse – Lorsque la sortie du tuyau est obturée, le liquide est en

équilibre hydrostatique. La pression dans le tuyau est uniforme et
vaut 𝑝atm + 𝜌𝑔ℎ0 . Lorsque le liquide coule, les colonnes d’eau verti-
cales sont en quasi-équilibre (le niveau descend très lentement) de
sorte que la hauteur d’eau nous renseigne sur la pression qui règne
 52 6 PERTES DE CHARGE ET FRICTION

dans le tuyau. On en déduit que la pression diminue le long de la

canalisation. On dit qu’il y a une perte de charge en ligne.

Ce résultat est en contradiction avec le théorème de Bernoulli. En

effet, si on l’applique entre le point A situé à la surface libre du large
récipient, et un point B situé dans la canalisation, on obtient

1
𝑝B = 𝑝atm + 𝜌 (𝑣A 2 − 𝑣B 2 ) + 𝜌𝑔ℎ(𝑡)
⏟⏟
2 ⏟⏟⏟⏟⏟⏟⏟⏟⏟
𝑓(𝑡)

On prévoit donc que la pression en B dépend du temps mais pas de

sa position, puisque ℎ, 𝑣A = −ℎ̇ et 𝑣B = 𝑣A (𝑑2 /𝐷2 ) ne dépendent
que du temps. Ce qui est en contradiction avec nos observations. On
en déduit que l’on ne peut pas appliquer le théorème de Bernoulli,
car l’on ne peut pas considérer le fluide parfait : c’est la viscosité du
fluide qui explique le phénomène de perte de charge.

Écoulement de Poiseuille

Cherchons tout d’abord à caractériser l’écoulement qui s’établit le

long d’une conduite cylindrique horizontale dans laquelle s’écoule
un fluide newtonien de viscosité 𝜂. On note 𝑑 le diamètre du tuyau,
𝑅 son rayon et 𝐿 sa longueur.

Hypothèses de travail –

1. L’écoulement est supposé stationnaire, d’où 𝜕 𝑣/𝜕𝑡 ⃗ = 0.⃗

2. L’écoulement est laminaire, parallèle à O𝑧 et invariant par rota-
tion autour de l’axe Oz, d’où 𝑣 ⃗ = 𝑣(𝑟, 𝑧) 𝑢
⃗⃗⃗⃗𝑧 .
3. Le fluide est incompressible.
4. Enfin, la pesanteur est négligeable sur l’épaisseur de la conduite.

M• 𝑣(M)
⃗⃗⃗⃗ = 𝑣(𝑟, 𝑧) ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧
𝑧
FIG. 6.1 : Écoulement de Poiseuille. Po-
sition du problème. fluide visqueux (𝜂, 𝜌)

Commençons par écrire l’équation de continuité :

𝜕(𝑟𝑣𝑟 ) 𝜕(𝑣𝜃 ) 𝜕𝑣𝑧 𝜕𝑣

div𝑣 ⃗ = + + =0= 𝑧 ⇒ 𝑣 ⃗ = 𝑣(𝑟) ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧
𝑟𝜕𝑟 𝑟𝜕𝜃 𝜕𝑧 𝜕𝑧
La vitesse ne dépend que de 𝑟. Calculons l’accélération :

𝜕𝑣⃗ 𝜕
𝑎(M,
⃗ 𝑡) = +𝑣 𝑢𝑧 ) = 0⃗
(𝑣(𝑟) ⃗⃗⃗⃗⃗⃗⃗⃗
𝜕𝑡 𝜕𝑧
L’accélération est nulle. En effet, les lignes de courant sont des droites
horizontales et se confondent avec la trajectoire des particules (ré-
gime stationnaire). Or, si la vitesse ne dépend pas de z cela signifie
que les particules de fluide se déplacent avec une vitesse constante
en direction et en intensité. L’accélération est donc nulle. On peut aus-
si ajouter que chaque particule de fluide est soumise à deux forces
qui se compensent : les forces de pression et les forces de viscosité.
 6.1 Pertes de charge dans une conduite 53

Sans force de pression, c’est-à-dire sans différence de pression, il ne

peut pas avoir d’écoulement stationnaire.

L’équation de Navier-Stokes se réduit à ⃗⃗⃗⃗⃗⃗

∇𝑝 = 𝜂Δ𝑣.⃗ Projetons cette
relation dans la base cylindrique :

𝜕𝑝
∣ = 0
∣ 𝜕𝑟
∣ 𝜕𝑝 d𝑝 1 d d𝑣
∣ 𝑟𝜕𝜃 = 0 ⇒ =𝜂 (𝑟 )
d𝑧 𝑟 d𝑟 d𝑟
∣
∣ 𝜕𝑝 1 d d𝑣
= 𝜂 (𝑟 )
∣ 𝜕𝑧 𝑟 d𝑟 d𝑟

Ainsi, la pression ne dépend que de 𝑧. Le terme de gauche de la der-

nière équation ne dépend donc que de 𝑧 alors que celui de droite ne
dépend que de 𝑟. Cette relation apparemment paradoxale se résout
si les deux termes sont constants.
d𝑝 1 d d𝑣
=𝐾 et 𝜂 (𝑟 ) = 𝐾
d𝑧 𝑟 d𝑟 d𝑟

Intéressons-nous au champ de vitesse. En intégrant deux fois on ob-

tient
𝐾 2
𝑣(𝑟) = 𝑟 + 𝐶1 ln 𝑟 + 𝐶2
4𝜂
où 𝐶1 et 𝐶2 sont deux constantes d’intégration. La vitesse doit être
définie en 𝑟 = 0 ce qui implique 𝐶1 = 0. Enfin, les conditions aux
limites imposent 𝑣(𝑅) = 0 d’où

−𝐾 2 𝑟 2
𝑣(𝑟) = (𝑅 − 𝑟2 ) = 𝑣max [1 − ( ) ]
4𝜂 𝑅

Le profil des vitesses est parabolique et l’écoulement est le plus

FIG. 6.2 : Profil de vitesse.

rapide le long de l’axe O𝑧.

Loi de Poiseuille

Revenons sur la relation relative à la pression

d𝑝
=𝐾 avec 𝐾 < 0
d𝑧

où le signe de 𝐾 est dicté par le sens de l’écoulement (ici selon O𝑧).

On en déduit
𝑝(𝑧) = 𝑝(0) + 𝐾𝑧
La pression diminue linéairement le long de la conduite comme on
l’a observé expérimentalement.
54 6 PERTES DE CHARGE ET FRICTION

Définition
On appelle perte de charge, la chute de pression entre l’aval et
l’amont d’une canalisation. Elle s’exprime en Pa et sera notée Δ𝑝
dans ce cours.

Ici la perte de charge s’écrit

Δ𝑝 = 𝑝(0) − 𝑝(𝐿) = −𝐾𝐿 > 0

Déterminons maintenant la relation entre la perte de charge et le

débit volumique. Le débit volumique est le flux du vecteur vitesse à
travers une section de la canalisation :
𝑅 4
𝐾𝜋𝑅
𝐷𝑉 = ∬ 𝑣 ⃗ ⋅ d𝑆⃗⃗⃗ ⃗ = ∫ 𝑣(𝑟) 2𝜋𝑟 d𝑟 = −
0 8𝜂

En remplaçant 𝐾 par −Δ𝑝/𝐿 on aboutit à la loi de Poiseuille :

8𝜂𝐿
Δ𝑝 = 𝐷 [Loi de Poiseuille] ♡ (6.1)
𝜋𝑅4 𝑉

On note que la perte de charge double quand la longueur de tuyau

double. En revanche, elle est très sensible au diamètre puisqu’elle
est divisée par 16 dans le diamètre double.

Exercice – On veut perfuser un patient en 1 H avec un flacon de 0,5 L de

plasma de viscosité 𝜂 = 1,4 × 10−3 Pa ⋅ s et de densité proche de l’eau. Si
l’aiguille utilisée a une longueur de 3 cm et un diamètre de 0,4 mm, à quelle
hauteur minimale faut-il installer le flacon ?
Solution : 95 cm.

Analogie électrique

On peut faire une analogie avec la conduction électrique et définir

une résistance hydraulique 𝑅𝐻 analogue de la résistance électrique :

Concepts électriques Concepts hydrauliques

Potentiel 𝑉 Pression 𝑝
ddp 𝑈 = 𝑉1 − 𝑉2 Différence de pression Δ𝑝
Intensité du courant électrique 𝐼 Débit volumique 𝐷𝑉
Loi d’Ohm 𝑈 = 𝑅𝐼 Loi de Poiseuille Δ𝑝 = 𝑅𝐻 𝐷𝑉

Formule générale

La formule de Poiseuille est valable tant que l’écoulement est lami-

naire. L’expérience montre qu’il faut pour cela que le nombre de Rey-
nolds ne dépasse pas 2 × 103 :

8𝜂𝐿 𝜌⟨𝑣⟩𝑑
Δ𝑝 = 𝐷 si 𝑅e = < 2 × 103 ♡ (6.2)
𝜋𝑅4 𝑉 𝜂
 6.1 Pertes de charge dans une conduite 55

Dans le cas contraire, l’écoulement est instable et la turbulence se

développe soit de façon intermittente soit de façon totale. Il faut alors
utiliser des lois empiriques. Ces lois reposent dans un premier temps
sur une analyse dimensionnelle.

Admettons que la perte de charge dépende seulement de la longueur

𝐿 de la conduite, de son diamètre 𝑑, de la viscosité 𝜂 du fluide, de
sa masse volumique 𝜌 ainsi que de la vitesse moyenne ⟨𝑣⟩ d’écou-
lement. Avec ces 5 grandeurs, on peut former deux nombres sans
dimensions :

• Le nombre de Reynolds 𝑅e = 𝜌⟨𝑣⟩𝑑/𝜂 ;

• le rapport 𝑥 = 𝐿/𝑑.

Ainsi, la perte de charge Δ𝑝 dépend de ces deux grandeurs adimen-

sionnées et de 3 autres grandeurs indépendantes. Choisissons 𝜌, ⟨𝑣⟩
et 𝑑 de sorte que Δ𝑝 doit pouvoir se mettre sous la forme
𝛽
Δ𝑝 = 𝑓(𝑅e,
⏟ 𝑥) 𝜌𝛼 ⟨𝑣⟩ 𝑑𝛾
sans dimension

Il nous reste à déterminer les exposants 𝛼, 𝛽 et 𝛾 pour que la relation

1 : Il suffit d’écrire l’équation aux di-
soit dimensionnellement correcte. On trouve1 𝛼 = 1, 𝛽 = 2 et 𝛾 = 0, mensions :
d’où
2 [𝐹 ] MLT−2
Δ𝑝 = 𝑓(𝑅e, 𝑥) × 𝜌⟨𝑣⟩ [Δ𝑝] = = = ML−1 T−2
[𝑆] L2
Par ailleurs, l’expérience montre que la perte de charge en ligne est et
proportionnelle à 𝐿, donc proportionnelle à 𝑥. C’est pourquoi l’ex- 𝛽 𝛼 𝛽
[𝜌𝛼 ⟨𝑣⟩ 𝑑𝛾 ] = (ML−3 ) (LT−1 ) L𝛾
pression de la perte de charge peut se mettre sous la forme
= M𝛼 L−3𝛼+𝛽+𝛾 T−𝛽

1 2𝐿 Par identification, on obtient

Δ𝑝 = 𝜆(𝑅e ) 𝜌⟨𝑣⟩ ♡ (6.3)
2 𝑑 ⎧ 𝛼 = 1
{
⎨ −3𝛼 + 𝛽 + 𝛾 = −1
Cette relation est appelée équation de Darcy-Weisbach. Le facteur adi- {
⎩ −𝛽 = −2
mensionné 𝜆 désigne le coefficient de perte de charge régulière et ce qui mène au résultat.
ne dépend que du nombre de Reynolds pour une canalisation lisse.

Dans le cas d’une canalisation rugueuse, un troisième nombre sans

dimension intervient : 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. Il suffit donc de réaliser des expériences
avec des écoulement de différents nombre de Reynolds et des tuyaux
de différentes rugosité pour déterminer de façon empirique la rela-
tion 𝜆(𝑅e , 𝜖/𝑑). Les résultats empiriques se lisent sur le diagramme
de Moody (FIG. 6.3).

Quelques remarques sur ce diagramme :

• Il s’agit d’un diagramme logarithmique. Pour le lire, on déter-

mine la rugosité relative pour savoir quelle courbe il faut sé-
lectionner (axe de droite). Ensuite, il suffit de suivre la courbe
sélectionnée jusqu’à la valeur du nombre de Reynolds désirée,
puis de lire à gauche la valeur du coefficient de perte de charge.
• En régime laminaire les pertes de charge sont données par la
loi de Poiseuille qui correspond à un coefficient de pertes de
charge
64
𝜆= pour 𝑅e < 2 × 103
𝑅e
 56 6 PERTES DE CHARGE ET FRICTION

FIG. 6.3 : ©Donebythesecondlaw at the English language Wikipedia [GFDL]

Cette loi se matérialise par une droite décroissante dans le dia-

gramme de Moody (« laminar flow »)

Exercice – Un tube en acier de 25 mètres de long, de diamètre intérieur

39,5 mm, de rugosité 𝜖 = 0,1 mm achemine de l’eau froide à 15°C (𝜌 =
998 g/L et 𝜂 = 1,14 mPa ⋅ s) avec un débit de 2,00 L/s. Calculer la perte de
charge en Pa.
Solution : 74 kPa

En régime turbulent, 𝜆 augmente brutalement et pour les grands nombres

de Reynolds, le coefficient de perte de charge conserve une valeur constante
qui ne dépend que de la rugosité relative de la conduite. Dans ce cas, les
pertes de charge varient comme le carré du débit. C’est pourquoi, l’ana-
logie avec la loi d’ohm comme la notion de résistance hydraulique n’ont
plus de sens en régime turbulent.

Pertes de charge singulière

Les pertes de charges en ligne étudiées jusqu’ici sont appelées pertes

de charge régulières.

Lorsque dans une conduite, l’écoulement subit de brusques varia-

tions de section ou de direction, il se produit des pertes de charges
dites singulières. Elles sont données par la formule

1 2
Δ𝑝s (en Pa) = 𝐾 𝜌⟨𝑣avant ⟩
2
 6.1 Pertes de charge dans une conduite 57

où ⟨𝑣avant ⟩ est la vitesse d’écoulement avant le passage de l’obstacle,

et 𝐾 le coefficient de perte de charge singulière. Ce dernier dépend
de la nature de l’obstacle (FIG. 6.4).

𝑆2
𝑆1
Divergent
2
𝑆1
𝐾 = (1 − ) si 𝛼 > 22,5°
𝑆2 𝛼

𝑆1 𝑆2
Convergent
2√
1 𝑆
𝐾= (1 − 2 ) sin 𝛼 si 𝛼 > 22,5°
2 𝑆1 𝛼

Coude brusque 𝛼
2 𝛼 4
𝐾 = sin 𝛼 + 2 sin ( )
2

Coude arrondi
𝑑 𝛼
7/2
2𝛼 𝑑 FIG. 6.4 : Coefficients de perte de
𝐾= [0,131 + 1,847 ( ) ] 𝑅
𝜋 2𝑅 • charge singulière sur quelques géomé-
tries.

Théorème de Bernoulli généralisé

Nous avons vu dans le chapitre consacré aux fluides parfaits, que

dans le champ de pesanteur, un fluide incompressible en écoulement
stationnaire voit la quantité 21 𝜌𝑣2 + 𝑝 + 𝜌𝑔𝑧 se conserver le long d’une
ligne de courant ce qui traduit la conservation de l’énergie. Voyons
comment s’exprime le bilan d’énergie dans le cas d’un fluide réel en
écoulement stationnaire dans une conduite en tenant compte des
échanges mécaniques avec des machines hydrauliques.

Considérons un fluide en écoulement stationnaire et incompressible

dans un système de conduites où il traverse des machines hydrau-
liques avec lesquelles il peut échanger de l’énergie :

• des pompes donneront de l’énergie mécanique au fluide ;

• des turbines recevront de la part du fluide de l’énergie méca-
nique.

Si l’on note 𝒫 la puissance échangée avec le fluide, on a 𝒫 > 0 pour

les pompes et 𝒫 < 0 pour les turbines.

Machine hydraulique

Puissance fournie : 𝒫
𝑝𝐵
B B’ ⃗⃗⃗⃗⃗⃗⃗⃗
𝑣 B

A A’
𝑝𝐴
𝑣A Pertes de charge : Δ𝑝
⃗⃗⃗⃗⃗⃗⃗⃗

FIG. 6.5 : Ecoulement stationnaire dans

d𝑚 = 𝜌𝑣A d𝑡 𝑆A un système de conduites.
58 6 PERTES DE CHARGE ET FRICTION

Bilan d’énergie cinétique – Considérons comme système, le fluide

situé dans la conduite entre A et B à l’instant 𝑡 et entre A’ et B’ à
l’instant 𝑡 + d𝑡. Pendant d𝑡, la masse transférée est

d𝑚 = 𝐷𝑚 d𝑡 = 𝜌𝐷𝑉 d𝑡

Rappelons que le débit volumique est uniforme le long de la canali-

sation puisque le fluide est incompressible.

Pour simplifier on considère que les grandeurs physiques sont uni-

formes sur la section droite de la conduite. Le théorème de l’énergie
cinétique donne

dℰc = ℰc (𝑡 + d𝑡) − ℰc (𝑡) = 𝛿𝑊

avec 𝛿𝑊 le travail de toutes les forces (extérieures et intérieures)

s’exerçant sur le système.

Tout d’abord, le régime étant permanent, la portion de fluide située

entre A’ et B conserve son énergie cinétique de sorte que

1
ℰc (𝑡 + d𝑡) − ℰc (𝑡) = ℰBB’ AA’
c − ℰc = 𝐷 d𝑡 (𝑣𝐵
2 2
− 𝑣𝐴 )
2 𝑚
Par ailleurs, les forces de pesanteur travaillent d’où le transfert mé-
canique

d𝑊𝑔 = −𝐷𝑚 d𝑡 𝑔 (𝑧𝐵 − 𝑧𝐴 ) [axe ascendant]

Les machines hydrauliques fournissent une puissance mécanique 𝒫

de sorte que le transfert mécanique effectué pendant la durée d𝑡
s’écrit
𝛿𝑊méca = 𝒫 d𝑡
Quant aux forces de pression, leur travail s’exprime par

𝛿𝑊𝑝 = −𝑝ext d𝑉 = 𝑝𝐴 𝐷𝑉 d𝑡 − 𝑝𝐵 𝐷𝑉 d𝑡

Le fluide est également le siège de forces intérieures. Or, l’écoule-

ment étant incompressible, les particules de fluide se déforment sans
changer de volume ce qui explique que les forces de pression interne
ne travaillent pas. Par contre, le fluide est le siège d’un travail résis-
tant 𝛿𝑊𝜂 des forces visqueuses. Par définition, la perte de charge ex-
primée en Pa correspond au travail des forces visqueuses par unité
de volume :
𝛿𝑊𝜂 = −𝐷𝑉 d𝑡 × Δ𝑝 ♡ (6.4)

Finalement le théorème de l’énergie cinétique donne la relation de

Bernoulli généralisée

1 2 1 2
𝜌 𝑣 + 𝜌 𝑔 𝑧 𝐵 + 𝑝𝐵 = 𝜌 𝑣 𝐴 + 𝜌 𝑔 𝑧𝐴 + 𝑝𝐴 + 𝒫/𝑄𝑉 − Δ 𝑝 ♡
2 𝐵 2
(6.5)
Notons au passage que l’on retrouve la relation de Bernoulli vu au
Chapitre 4 à condition de faire 𝒫 = 0 et Δ𝑝 = 0.
 6.2 Résistance au mouvement dans les fluides 59

La relation (6.5) utilise l’approximation des écoulements unidimension-

nels, ce qui revient à confondre la vitesse avec la vitesse moyenne sur
une section droite de la conduite. Cette approximation produit des er-
reurs sur l’expression de l’énergie cinétique. En effet, le bilan rigoureux
de l’énergie donne

𝒫
𝑝1 + ⟨𝑒𝑐1 ⟩ + 𝜌𝑔𝑧1 + = 𝑝2 + ⟨𝑒𝑐2 ⟩ + 𝜌𝑔𝑧2 + Δ𝑝𝜂
𝑄𝑉

où ⟨𝑒𝑐 ⟩ est la moyenne de l’énergie cinétique volumique défini par

𝑣2 1 1 1
⟨𝑒c ⟩ = ⟨𝜌 ⟩= ∬ ( 𝜌𝑣2 ) 𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 = 𝛼 𝜌⟨𝑣⟩2
2 𝑄𝑉 2 2

où 𝛼 est un coefficient correctif qui dépend du profil de vitesse dans la

canalisation. En pratique les valeurs de 𝛼 sont les suivantes :

• 𝛼 = 1 écoulement piston (vitesse uniforme dans toute la section)

• 𝛼 = 2 écoulement laminaire visqueux (newtonien)
• 𝛼 = 1, 01 − 1, 1 écoulement turbulent.

Comme le terme d’énergie cinétique est souvent faible devant les autres
termes, une erreur sur celle ci a peu d’effet sur le résultat. C’est pourquoi
nous prendrons systématiquement 𝛼 = 1.

6.2 Résistance au mouvement dans les fluides

Si l’on met de coté la poussée d’Archimède, la force que ressent un

solide plongé dans un écoulement stationnaire tridimensionnelle est
nulle si le fluide n’est pas visqueux. En revanche, l’écoulement vis-
queux autour d’un obstacle solide produit une force qui présente
deux composantes : la composante dans le sens de l’écoulement est
appelée force de traînée, la composante perpendiculaire est la force
de portance.

Formule de Stokes

Stokes s’est intéressé à la force de traînée qu’un écoulement visqueux

produit autour d’une sphère. Il s’est placé dans le cas où l’écoulement
𝑣
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
est gouverné par la viscosité c’est-à-dire pour les petits nombres de ∞
⃗⃗⃗⃗⃗⃗
𝐹
Reynolds.

La résolution complète est assez longue et nous allons nous conten-

FIG. 6.6 : Traînée sur un obstacle sphé-
ter de la solution sans chercher à la justifier. Stokes obtient qu’une rique immobile.
sphère de rayon 𝑟, immobile, soumise à un écoulement permanent
incompressible et visqueux, ressent une force de traînée 𝐹 ⃗⃗⃗ ⃗⃗ propor-
tionnelle à la vitesse d’écoulement et à la taille de la sphère.

⃗⃗⃗ ⃗⃗ = 6𝜋𝜂𝑟 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗

𝐹 𝑣∞ pour 𝑅e ≪ 1 ♡ (6.6)

où ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑣∞ représente la vitesse de l’écoulement par rapport à la sphère
et loin de la sphère. Cette force de trainée est liée d’une part à un
champ de pression plus important en avant de la sphère et d’autre
part aux forces visqueuses.
 60 6 PERTES DE CHARGE ET FRICTION

Si l’on étudie la chute d’une bille sphérique dans un fluide visqueux

Vitesse de sédimentation au repos (loin de la bille), il faut écrire
La physique des suspensions (par-
ticules solides mélangées à un li- ⃗⃗⃗ ⃗⃗ = −6𝜋𝜂𝑟𝑣 ⃗
𝐹
quide) et des émulsions (goutte-
lettes liquides dispersées dans un où 𝑣 ⃗ représente la vitesse de la bille dans le référentiel du laboratoire.
autre liquide non miscible) utilise
Cette loi est vérifiée avec une précision meilleure que 1% pour 𝑅e <
la loi de Stokes car le nombre de
Reynolds est assez petit. Lorsqu’on 0, 3.
laisse reposer un liquide conte-
nant de petites particules solides
(comme par exemple un mélange
eau-argile), les particules vont dé-
Coefficient de traînée Cx
canter c’est-à-dire sédimenter au
fond du récipient avec une vitesse L’analyse précédente n’est valable qu’à petit nombre de Reynolds et
qui dépend de leur dimension ca-
pour un corps sphérique.
ractéristique. Le temps de décan-
tation donne alors un renseigne-
Pour un corps de géométrie différente, la force qu’exerce un fluide
ment sur la taille des grains. En
effet les grains tombent à une vi- en écoulement stationnaire unidirectionnel présente deux compo-
tesse constante pour laquelle le santes.
poids apparent (poids moins la
poussée d’Archimède) compense • Une composante parallèle à ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑣∞ : c’est la traînée (on dit aussi
la force de traînée : résistance) ⃗⃗⃗⃗⃗⃗⃗⃗
𝐹∥ .
6𝜋𝜂𝑟𝑣sed
4
= 𝜋𝑟3 (𝜌𝑠 − 𝜌𝑙 )𝑔 • Une composante perpendiculaire : c’est la portance ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹⟂ .
3
2
d’où 𝑣sed = 9𝜂 (𝜌𝑠 − 𝜌𝑙 )𝑟2 . ligne d’écoulement ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹⟂ ⃗⃗⃗⃗⃗⃗
𝐹

𝑣
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑧
∞
⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹∥ ⊗ 𝑥
𝑦
FIG. 6.7 : Portance et traînée.

Le traitement analytique est possible pour des géométries simples et

pour des valeurs faibles de 𝑅e comme nous venons de le voir. Pour
des grands nombres de Reynolds, on procède en général à des expé-
riences sur maquette en soufflerie pour avoir accès à la force. Mon-
trons par une analyse dimensionnelle quelle forme doit prendre cette
force 𝐹 .
Admettons que la force dépende de la vitesse 𝑣∞ , de la densité 𝜌 du
fluide, de sa viscosité 𝜂 et d’une longueur caractéristique ℓ associé
2 : Pour une sphère on choisit géné- à la taille du corps2 . Avec ces 4 grandeurs indépendantes, on peut
ralement son diamètre, pour une aile
construire un nombre sans dimension, à savoir le nombre de Rey-
d’avion sa corde.
nolds
𝜌𝑣 ℓ
𝑅e = ∞
𝜂
On choisit alors d’exprimer la force en fonction de 𝑅e , 𝜌, 𝑣∞ et ℓ :

𝐹 = 𝑓(𝑅e )𝜌𝛼 𝑣∞ 𝛽 ℓ𝛾

L’analyse dimensionnelle aboutit à

𝐹 = 𝑓(𝑅e )𝜌 𝑣∞ 2 ℓ2

On définit alors deux coefficients de frottement – le 𝐶𝑥 et le 𝐶𝑧 –

associés respectivement à la résistance et à la portance :

1 2 1 2
𝐹∥ = 𝜌𝑣 𝑆 𝐶𝑥 (𝑅e ) et 𝐹⟂ = 𝜌𝑣 𝑆 𝐶𝑧 (𝑅e )
2 ∞ 2 ∞
 6.2 Résistance au mouvement dans les fluides 61

où 𝑆 désigne une surface caractéristique de référence. On choisit gé-

néralement la surface projetée dans le plan (𝑦𝑧) pour la force de traî-
née, et dans le plan (𝑥𝑦) pour la force de portance. Ainsi, il suffit de
mesurer les forces aérodynamiques en soufflerie sur une maquette
en échelle réduite pour différentes valeurs de nombre de Reynolds
pour prévoir les forces en échelle réelle.

Sphère 1/2 Sphère Cône 60° Cube Cube à 45° Cylindre2x1 Corps profilé
sens de l’

écoulement
𝐶𝑥 = 0,47 0,42 0,50 1,05 0,80 0,85 0,04

FIG. 6.8 : Quelques 𝐶𝑥 pour différents profils dans un écoulement turbulent avec 𝑅e ≃ 105 .

Dans le cas d’une sphère (sans rotation propre) de rayon 𝑟, on ob-

tient

𝜋𝑟2 2 𝜌𝑣 (2𝑟)
𝐹 = 𝐹∥ = 𝜌𝑣∞ 𝐶𝑥 (𝑅e ) avec 𝑅e = ∞ ♡ (6.7)
2 𝜂

FIG. 6.9 : Évolution du 𝐶𝑥 d’une sphère

en fonction du nombre de Reynolds.

En général, on constate que le 𝐶𝑥 est quasi constant en régime tur-

bulent ce qui correspond aux situations courantes de l’aéronautisme,
le nautisme, le cyclisme etc. Pour 𝑅e ≪ 1 on a 𝐶𝑥 = 24/𝑅e et l’on
retrouve la loi de Stokes (6.6).

Exercice – Déterminer la vitesse limite de chute d’une goutte d’eau sphé-

rique de rayon 𝑟 = 1,0 mm dans l’air (20 °C et 1,0 bar) sachant que 𝐶𝑥 ≃ 0,5
pour 500 < 𝑅e < 105 .
Solution : 6,6 m ⋅ s−1 .
ANNEXES
 THÉORÈME DE LA DIVERGENCE A
Version en ligne
[Link]/mecanique_des_fluides/[Link]

Le théorème de la divergence relie le flux d’un champ vectoriel à tra-

vers une surface fermée S, à la somme d’un scalaire en tout point
du volume enfermé par S. Ce théorème fait appel à l’opérateur diver-
gence, d’où son nom. La divergence est un opérateur qui s’applique
à un champ vectoriel et retourne un champ scalaire. Il se note

⃗⃗⃗ ⃗⃗
div𝐴 ou ⃗⃗⃗⃗⃗⃗
∇⋅𝐴 ⃗⃗⃗ ⃗⃗

La dernière notation permet de retrouver son expression en coordon-

nées cartésiennes :

𝜕/𝜕𝑥 𝐴𝑥 (𝑥, 𝑦, 𝑧)
⃗⃗⃗⃗⃗⃗ ⎛𝜕/𝜕𝑦 ⎟
⃗⃗⃗ ⃗⃗ 𝑦, 𝑧) = ⎜ ⎛
⎞⋅⎜ ⎞ = 𝜕𝐴𝑥 + 𝜕𝐴𝑦 + 𝜕𝐴𝑧
∇ ⋅ 𝐴(𝑥, ⎜ ⎟ ⎜𝐴𝑦 (𝑥, 𝑦, 𝑧)⎟ ⎟
𝜕𝑥 𝜕𝑦 𝜕𝑧
⎝ 𝜕/𝜕𝑧 ⎠ ⎝ 𝐴𝑧 (𝑥, 𝑦, 𝑧) ⎠
⃗⃗⃗ ⃗ ⃗
𝐴(P)
𝑎
Imaginons un petit cube d’arête 𝑎 centré en M(𝑥, 𝑦, 𝑧), et dont les
faces sont des plans cartésiens. Voyons comment s’écrit le flux d’un ⃗⃗⃗⃗⃗
𝑛
•
champ vectoriel 𝐴⃗⃗⃗ ⃗⃗ à travers ce cube. P(𝑥′ , 𝑦 + 𝑎2 , 𝑧′ )
(2)
•
Commençons par exprimer le flux 𝜙1 à travers la face (1) perpendicu- M(𝑥, 𝑦, 𝑧) (1)
𝑧
laire à l’axe (O𝑦) et située en 𝑦 + 𝑎/2 :

⃗⃗⃗ ⃗⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗ 𝑦
𝜙1 = ∬ 𝐴 𝑢𝑦 d𝑆 = ∬ 𝐴𝑦 (𝑥′ , 𝑦 + 𝑎/2, 𝑧 ′ ) d𝑥′ d𝑧′ 𝑥
(1) (1)
FIG. A.1 : Calcul du flux à travers un
cube.
De même, le flux à travers la face située en 𝑦 − 𝑎/2 vaut

⃗⃗⃗ ⃗⃗ ⋅ (−⃗⃗⃗⃗⃗⃗⃗⃗ Par convention, la normale ⃗⃗⃗⃗⃗𝑛 est per-

𝜙2 = ∬ 𝐴 𝑢𝑦 ) d𝑆 = ∬ −𝐴𝑦 (𝑥′ , 𝑦 − 𝑎/2, 𝑧 ′ ) d𝑥′ d𝑧 ′ pendiculaire à la surface et orientée
(2) (2)
vers l’extérieur ; ce qui explique la pré-
sence du signe −.
Appelons 𝜙𝑦 le flux à travers ces deux faces, et faisons tendre 𝑎 → 0.
On peut alors considérer l’intégrand constant et égale à sa valeur au
centre de la face :

𝜙𝑦 = 𝜙1 + 𝜙2 = [𝐴𝑦 (𝑥, 𝑦 + 𝑎/2, 𝑧) − 𝐴𝑦 (𝑥, 𝑦 − 𝑎/2, 𝑧)] 𝑎2

𝑎 étant un infiniment petit, on peut légitimement remplacer

𝐴𝑦 (𝑥, 𝑦 + 𝑎/2, 𝑧) − 𝐴𝑦 (𝑥, 𝑦 − 𝑎/2, 𝑧) 𝜕𝐴𝑦 (𝑥, 𝑦, 𝑧)

par
𝑎 𝜕𝑦

𝜕𝐴𝑦 (𝑥, 𝑦, 𝑧) 3
ce qui donne 𝜙𝑦 = 𝑎 . Ce même raisonnement réitéré sur
𝜕𝑦
les faces perpendiculaires aux axes (O𝑥) et (O𝑧) aboutit à

𝜕𝐴𝑥 (𝑥, 𝑦, 𝑧) 3 𝜕𝐴𝑧 (𝑥, 𝑦, 𝑧) 3

𝜙𝑥 = 𝑎 et 𝜙𝑧 = 𝑎
𝜕𝑥 𝜕𝑧
 66 A THÉORÈME DE LA DIVERGENCE

Finalement, le flux 𝜙 du champ vectoriel 𝐴 ⃗⃗⃗ ⃗⃗ à travers un cube infinité-

simal centré en (𝑥, 𝑦, 𝑧) de volume infinitésimal d𝜏 = 𝑎3 vaut

𝜕𝐴𝑥 (𝑥, 𝑦, 𝑧) 𝜕𝐴𝑦 (𝑥, 𝑦, 𝑧) 𝜕𝐴𝑧 (𝑥, 𝑦, 𝑧) ⃗⃗⃗ ⃗⃗ d𝜏

𝜙=( + + ) d𝜏 = div𝐴
𝜕𝑥 𝜕𝑦 𝜕𝑧

Mettons dorénavant côte à côte deux cubes infinitésimaux. Lorsque

𝜙1 = −𝜙2

(2) ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑛1 • 𝑛2 (1)
⃗⃗⃗⃗⃗⃗⃗⃗⃗
FIG. A.2 : Le flux à travers deux cubes
adjacents se réduit au flux à travers la
surface qui délimite le volume consti-
tué par la réunion des deux cubes.

l’on calcule le flux à travers ces deux cubes réunis, on s’aperçoit que
la contribution due aux surfaces adjacentes se compensent, car les
normales à ces faces sont opposées. Aussi, le flux total se réduit-il au
flux à travers la surface frontière. Dès lors, on conçoit qu’en empilant
de tels cubes en nombre infini, on puisse reconstituer un volume
fini, de sorte que le flux à travers la surface frontière soit égal à la
somme des flux élémentaires produits à travers chaque petit cube
constituant le volume. C’est le sens du théorème de la divergence.

Théorème de la divergence

L’intégrale de la divergence d’un champ vectoriel sur un volume

V est égal au flux de ce champ à travers la surface fermée qui
délimite le volume.

∯ ⃗⃗⃗ ⃗⃗
𝐴(M) 𝑛 d𝑆 = ∭
⋅ ⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗
div𝐴(M) d𝜏
M∈S M∈V

Ce théorème est aussi appelé théorème de Green-Ostrogradsky

ACCÉLÉRATION D’UNE PARTICULE
DE FLUIDE B
Ce complément donne une expression alternative de l’accélération
faisant intervenir les opérateurs gradient et rotationnel.

Version en ligne
[Link]/mecanique_des_fluides/acceleration-
[Link]

On a vu dans le Chapitre 3 consacré à la cinématique des fluides,

que l’accélération d’une particule de fluide située en M à l’instant 𝑡
pouvait s’obtenir à l’aide du champ de vitesse 𝑣(M,
⃗ 𝑡) :

𝜕𝑣⃗
𝑎(M,
⃗ 𝑡) = + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗
𝜕𝑡
où le dernier terme désigne la partie convective de l’accélération. Ex-
plicitons la composante suivant O𝑥 de ce terme en utilisant l’identité
du double produit vectoriel :

⃗⃗⃗ ⃗⃗ ∧ (𝐵
𝐴 ⃗⃗⃗ ⃗⃗ ∧ 𝐶) ⃗⃗⃗ ⃗⃗ ⋅ 𝐶)
⃗⃗⃗ ⃗⃗ = (𝐴 ⃗⃗⃗ ⃗⃗ 𝐵 ⃗⃗⃗ ⃗⃗ ⋅ 𝐵)
⃗⃗⃗ ⃗⃗ − (𝐴 ⃗⃗⃗ ⃗⃗
⃗⃗⃗ ⃗⃗ 𝐶

⃗⃗⃗ ⃗⃗ = 𝑣,⃗ 𝐵
Prenons 𝐴 ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗ = 𝑢
∇𝑣𝑥 et 𝐶 ⃗⃗⃗⃗𝑥 :

𝑣 ⃗ ∧ (⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 ∧ ⃗⃗⃗⃗⃗⃗⃗⃗ 𝑢𝑥 ) ⃗⃗⃗⃗⃗⃗
𝑢𝑥 ) = (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗ ∇𝑣𝑥 − (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 ) ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥

Le premier membre s’écrit aussi 𝑣 ⃗ ∧ (⃗⃗⃗⃗⃗⃗

∇ ∧ 𝑣𝑥 𝑢 𝑥 ). Quant au premier
⃗⃗⃗⃗⃗⃗⃗⃗
terme du second membre, on a
1
𝑢𝑥 ) ⃗⃗⃗⃗⃗⃗
(𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗ ∇𝑣𝑥 = 𝑣𝑥 ⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 = ⃗⃗⃗⃗⃗⃗
∇𝑣 2
2 𝑥

Enfin le dernier terme est équivalent à (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗ 𝑢𝑥 . Ainsi on obtient

∇) 𝑣𝑥 ⃗⃗⃗⃗⃗⃗⃗⃗

1
𝑣 ⃗ ∧ (⃗⃗⃗⃗⃗⃗ 𝑢𝑥 ) = ⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣𝑥 ⃗⃗⃗⃗⃗⃗⃗⃗ ∇𝑣 2 − (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣𝑥 𝑢⃗⃗⃗⃗⃗⃗⃗⃗
2 𝑥 𝑥

En procédant de la même façon pour les deux autres composantes,

on obtient, après avoir sommé les trois relations,

1 ⃗⃗⃗⃗⃗⃗
𝑣 ⃗ ∧ (⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣)⃗ = ∇ (𝑣𝑥 2 + 𝑣𝑦 2 + 𝑣𝑧 2 ) − (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗
2
où l’on reconnait 𝑣2 dans le gradient. En substituant dans l’expression
de l’accélération, on aboutit à une nouvelle expression de l’accéléra-
tion :
𝜕 𝑣 ⃗ ⃗⃗⃗⃗⃗⃗ 𝑣2
𝑎(M,
⃗ 𝑡) = + ∇ ( ) + (⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣)⃗ ∧ 𝑣 ⃗
𝜕𝑡 2
L’opérateur ⃗⃗⃗⃗⃗⃗
∇∧ s’appelle le rotationnel ; il s’applique à un champ vec-
toriel et retourne également un champ vectoriel. Il se note indifférem-
ment
rot ou ⃗⃗⃗⃗⃗⃗
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ ∇∧
68 B ACCÉLÉRATION D’UNE PARTICULE DE FLUIDE

On voit ici que l’accélération fait intervenir le rotationnel de la vitesse,

que l’on appelle vorticité en mécanique des fluides. Pour résumer,

𝜕 𝑣 ⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑣2
𝑎(M,
⃗ 𝑡) = + grad ( ) + ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot (𝑣)⃗ ∧ 𝑣 ⃗ (B.1)
𝜕𝑡 2
 THÉORÈME π C
Version en ligne C.1 Énoncé . . . . . . . . . . 69
[Link]/mecanique_des_fluides/[Link] C.2 Applications en mécanique
des fluides . . . . . . . . 69
Force de traînée . . . . . 69
Pertes de charge . . . . . 70
C.1 Énoncé Convection forcée . . . . 71

Le théorème π ou théorème de Vaschy-Buckingham est le théorème

fondamental de l’analyse dimensionnelle.

Supposons que nous cherchions une relation entre 𝑛 grandeurs phy-

siques 𝑔𝑖=1..𝑛 que l’on considère pertinentes pour décrire un phéno-
mène. Notons 𝑘 le nombre de dimensions fondamentales utilisées
par ces grandeurs (𝑘 ≤ 7).
1 : Dans l’article original de Buckin-
Le théorème stipule qu’il existe alors (𝑛−𝑘) produits sans dimension gham, ces nombres sans dimension
notées1 𝜋𝑖 tels que 𝑓(𝜋1 , ..., 𝜋𝑛−𝑘 ) = 0. sont notés Π, d’où le nom du théo-
rème
Ainsi, procéder à une analyse dimensionnelle consiste à déterminer
ces nombres sans dimensions, ce qui ne pose pas de problème en gé-
néral. Une fois ces nombres définies, le problème se ramène alors à
une relation entre grandeurs adimensionnées. Cette dernière fait ap-
paraître une constante ou une fonction à déterminer empiriquement.
Ainsi, lorsque le problème est trop complexe pour une modélisation,
l’analyse dimensionnelle et un guide précieux avant toute analyse
empirique.

C.2 Applications en mécanique des fluides

Force de traînée sur une sphère

Plaçons un obstacle sphérique dans un écoulement stationnaire de

vitesse 𝑣∞ ⃗⃗⃗ ⃗⃗ que
⃗ loin de l’obstacle et cherchons la force de trainée 𝐹
𝑣
⃗⃗⃗⃗∞
produit l’écoulement sur l’obstacle. Que donne l’analyse dimension- ⃗⃗⃗⃗⃗⃗
𝐹
nelle ?

L’expérience montre que cette force est influencée par les facteurs FIG. C.1 : Traînée sur un obstacle sphé-
suivants : rique immobile.

• le diamètre 𝑑 de la sphère ;
• la masse volumique 𝜌 du fluide ;
• sa viscosité 𝜂 ;
• ainsi que la vitesse d’écoulement 𝑣∞ .

Rassemblons dans un tableau la dimension des cinq grandeurs que

l’on cherche à relier :
70 C THÉORÈME π

Grandeurs 𝑣∞ 𝑑 𝜌 𝜂 𝐹

Dimensions LT−1 L ML−3 ML−1 T−1 MLT−2

On a 𝑛 = 5 grandeurs et 𝑘 = 3 dimensions différentes. D’après le

théorème 𝜋, il existe 𝑛 − 𝑘 = 2 nombres sans dimensions 𝜋1 et 𝜋2 liés
entre eux. Il est facile de construire de tels nombres :
• On sait que le nombre de Reynolds est un facteur adimension-
né. On choisit donc
𝜌𝑣∞ 𝑑
𝜋1 = 𝑅e =
𝜂
1
• Par ailleurs, la quantité est homogène à une pression,
2
2 𝜌𝑣∞
comme le rapport 𝐹 /𝑑 . On pose alors
2

𝐹 /𝑑2
𝜋2 = 2
1/2𝜌𝑣∞

En vertu du théorème 𝜋 nous savons que

1 2 2
𝜋2 = 𝑓(𝜋1 ) soit 𝐹 = 𝜌𝑣 𝑑 𝑓(𝑅e ) (C.1)
2 ∞
Des études en soufflerie permettent d’accéder à la fonction inconnue
𝑓(𝑅e ).

Pertes de charges régulières

Considérons une conduite cylindrique horizontale traversée par un

fluide newtonien circulant avec un certain débit. On observe une dif-
férence de pression –dite perte de charge– entre l’amont et l’aval du
tuyau. Notons-là Δ𝑝𝜂 . L’expérience montre que les pertes de charge
sont influencées par les paramètres suivants :
• la longueur 𝐿 de la conduite ;
• son diamètre 𝑑 ;
• la masse volumique du fluide ;
• sa viscosité 𝜂 ;
• la vitesse moyenne ⟨𝑣⟩ de l’écoulement ;
• ainsi que l’état de surface de la canalisation que l’on peut ca-
ractériser par la rugosité 𝜖.
Déterminons les dimensions de ces septs grandeurs que l’on cherche
à relier.

Grandeurs 𝐿 𝑑 𝜌 𝜂 Δ𝑝𝜂 ⟨𝑣⟩ 𝜖

Dimensions L L ML−3 ML−1 T−1 ML−1 T−2 L.T−1 L

On a 𝑛 = 7 grandeurs et 𝑘 = 3 dimensions différentes impliquées

dans ces grandeurs. D’après le théorème 𝜋, il existe 𝑛−𝑘 = 4 nombres
sans dimension –notons-les 𝜋1 , 𝜋2 , 𝜋3 et 𝜋4 – qui doivent être reliées
par une loi. Cherchons donc à construire quatre nombres indépen-
dants sans dimension.
 C.2 Applications en mécanique des fluides 71

• 𝐿 et 𝑑 étant de même dimension, leur rapport est adimension-

né : 𝜋1 = 𝐿/𝑑.
• 𝜖 et 𝑑 sont aussi de même dimension : 𝜋2 = 𝜖/𝑑 (rugosité rela-
tive).
• On sait que le nombre de Reynolds est sans dimension ; ce sera
𝜌⟨𝑣⟩𝑑
𝜋3 = 𝑅 e = .
𝜂
• Enfin, le théorème de Bernoulli nous enseigne que 21 𝜌⟨𝑣⟩ 2 est
Δ𝑝𝜂
homogène à une pression. Ainsi 𝜋4 = .
1/2𝜌⟨𝑣⟩ 2
D’après le théorème 𝜋, il existe une relation entre ces quatre nombre
sans dimension. Par exemple :

𝜋4 = 𝑓(𝜋1 , 𝜋2 , 𝜋3 )

Autrement-dit,
Δ𝑝𝜂
1 2
= 𝑓(𝐿/𝑑, 𝜖/𝑑, 𝑅e )
2 𝜌⟨𝑣⟩

L’analyse dimensionnelle est ici terminée. La fonction 𝑓(𝐿/𝑑, 𝜖/𝑑, 𝑅e )

est à déterminer, soit de façon empirique, soit à l’aide d’une modéli-
sation2 2 : cf. la loi de Poiseuille par exemple
au Chapitre 6.
Empiriquement, on constate que la perte de charge Δ𝑝𝜂 est propor-
tionnelle à 𝐿. Par conséquent, on peut écrire

𝐿
𝑓(𝐿/𝑑, 𝜖/𝑑, 𝑅e ) = × 𝜆(𝜖/𝑑, 𝑅e )
𝑑
où 𝜆(𝜖/𝑑, 𝑅e ) est une fonction inconnue à déterminer. Pour résumer,
on a
1 𝐿
Δ𝑝𝜂 = 𝜆(𝜖/𝑑, 𝑅e ) 𝜌⟨𝑣⟩ 2
2 𝑑
Cette relation est appelée équation de Darcy-Weisbach. Le facteur adi-
mensionné 𝜆 désigne le coefficient de perte de charge régulière et
cette analyse prévoit qu’il ne dépend que du nombre de Reynolds et
de l’état de surface de la canalisation via la rugosité relative.

Convection forcée

La convection est un mode de transfert thermique qui implique un

déplacement collectif de fluide. La matière fluide chaude, en se dé-
plaçant, cède de l’énergie aux parties plus froides. Lorsque le déplace-
ment du fluide est forcé par des pompes ou des turbines par exemple,
on parle de convection forcée.

Prenons le cas d’un simple échangeur cylindrique de diamètre inté- 𝑑 ⃗⃗ fluide(𝜆, 𝑐𝑝 , 𝜌, 𝜂)

⟨⃗⃗𝑣⟩
rieur 𝑑 dans lequel circule un fluide de masse volumique 𝜌, de viscosi-
té 𝜂, de conductivité thermique 𝜆 et de capacité thermique massique
𝑐𝑝 . Lorsque le fluide s’écoule à une vitesse moyenne ⟨𝑣⟩, un échange FIG. C.2 : Convection dans un échan-
convectif a lieu avec l’intérieur du tuyau. En vertu de la loi de New- geur cylindrique.

ton, le courant thermique (en W ⋅ m−2 ) est proportionnel à l’écart de

température entre le fluide et le solide. On a

𝑗th = ℎΔ𝑇
72 C THÉORÈME π

où ℎ est le coefficient de transfert. Ce dernier dépend essentielle-

ment de la nature du fluide et de son écoulement. Montrons sur la
géométrie de la FIG. C.2 comment l’analyse dimensionnelle permet
de prévoir des lois d’échelle concernant le coefficient de transfert ℎ.
Comme d’habitude, dressons une table avec toutes les grandeurs a
priori impliquées.

Grandeur 𝜌 ⟨𝑣⟩ 𝑑 𝜂 𝜆 𝑐𝑝 ℎ

Dimension ML−3 LT−1 L ML−1 T−1 MLT−3 Θ−1 L2 T−2 Θ−1 MT−3 Θ−1

On cherche une relation entre 𝑛 = 7 grandeurs mettant en jeu 𝑘 = 4

dimensions indépendantes. En vertu du théorème 𝜋, il existe 𝑛−𝑘 = 3
grandeurs adimensionnées liées entre elles. On peut effectivement
construire trois nombres sans dimension :
𝜌⟨𝑣⟩𝑑
• Le nombre de Reynolds 𝑅e = qui caractérise l’écoule-
𝜂
ment.
𝜂 𝑐𝑝
• Le nombre de Prandtl 𝑃r = . Plus le nombre de Prandtl est
𝜆
élevé, plus la diffusion de masse domine devant la diffusion de
chaleur.
ℎ𝑑
• Le nombre de Nusselt 𝑁u = . Plus le nombre de Nusselt
𝜆
est important plus le transfert par convection domine devant
le transfert par conduction.
D’après le théorème 𝜋, on a

𝜆
𝑁u = 𝑓(𝑅e , 𝑃r ) ⟹ ℎ= 𝑓(𝑅e , 𝑃r )
𝑑

La fonction 𝑓 peut être déterminée de façon empirique ou numérique.

Par exemple, pour les fluides usuels et lorsque le régime turbulent
est établi, on obtient
1/3
𝑁u = 0, 023 𝑅e0,8 𝑃r si 𝑃r ≥ 0, 5
Pour en savoir plus
[1] Pierre-Gilles de GENNES, Françoise Brochard WYART et David QUÉRÉ. Gouttes, bulles, perles et ondes.
Belin, 2002 (cf. p. 13).
[2] Jean-Pierre HULIN, Etienne GUYON et Emmanuel GUYON. « Pilotes et tourbillons ». In : bulletin de la
société française de physique (mai 2002), p. 4.
[3] S. CANDEL. Mécanique des fluides 2ème cycle. Paris : Dunod, 2001.
[4] P. FLEURY et J.P. MATHIEU. Mécanique Physique. Sous la dir. d’EYROLLES. 1961.
[5] Ludovic KEISER et al. « Fragmentation Marangoni ». In : Reflets de la Physique 59 (sept. 2018), p. 33-
35.
[6] D. LYNCH. « Les mascarets ». In : Pour la Science (déc. 1982).
[7] C. CLANET et L. BOCQUET. « Les ricochets ». In : Bulletin de la SFP 152 (déc. 2005).
[8] David QUÉRÉ, José BICO et Denis RICHARD. « Le mouillage nul (ou presque) ». In : Bulletin de la SFP
125 (juill. 2000), p. 8.
[9] Dominique LANGEVIN et al. « Les mousses ». In : Bulletin de la Société française de physique 134
(mai 2002), p. 11-16.
[10] G. DUPEU et al. « Le football et ses trajectoires ». In : Reflets de la physique 28 (mars 2012), p. 10-14.
[11] Marc RABAUD et Frédéric MOISY. « Du neuf dans les sillages ». In : Reflets de la physique 39 (2014),
p. 10-13.
[12] Tom SCHELLER. « Micro-écoulements dans des films liquides : application aux mousses aqueuses ».
In : Mémoire,(Université de Liège, Année académique 2004-2005) ().
[13] Lydéric BOCQUET. « The physics of stone skipping ». In : American Journal of Physics 71.2 (2003),
p. 150-155.
[14] Ian BRUCE. « Concerning drops ». In : American Journal of Physics 52.12 (1984), p. 1102-1105.
[15] Feredoon BEHROOZI et al. « The profile of a dew drop ». In : American Journal of Physics 64.9 (1996),
p. 1120-1125.
[16] Isabelle GALLAGHER. Autour des équations de Navier-Stokes — Images des Mathématiques, CNRS.
2010. URL : [Link]
[17] G. LIGER-BELAIR. « Histoire illustrée d’une bulle de champagne ». In : Bulletin de la SFP 127 (déc.
2000).
[18] D. QUÉRÉ et J-M Di MEGLIO. « Sur un thème de Rayleigh ». In : Pour la Science 152 (juin 1990).
[19] B. DARBOIS-TEXIER et al. « La caléfaction ». In : Reflets de la physique 37 (déc. 2013).
[20] A. Amon P. PANIZZA et L. COURBIN. « Comment circulent des gouttes dans un laboratoire sur puce ? »
In : Reflets de la physique 36 (oct. 2013).
[21] D. LANGEVIN. « La coalescence ». In : Bulletin de la SFP 115 (juill. 1998).
[22] D. LANGEVIN. « Les tensioactifs ». In : BUP, Hors série 875 (2005).
[23] D. QUÉRÉ et al. « Les surfaces superhydrophobes ». In : Images de la physique (2005), p. 239-244.
Grandeurs physiques et symboles
mathématiques

Constantes physiques définies par le SI (valeurs exactes)

−1
ℎ Constante de Planck 6,626 070 15 × 10−34 J ⋅ Hz

𝑐 Vitesse de la lumière dans le vide 299 792 458 m ⋅ s−1

Δ𝜈Cs Fréquence hyperfine du 133 Cs 9 192 631 770 Hz

𝑒 Charge élémentaire 1,602 176 634 × 10−9 C

𝑘B Constante de Boltzmann 1,380 649 × 10−23 J ⋅ K−1

−1
𝑁A Nombre d’Avogadro 6,022 140 76 × 1023 mol
−1
𝑅 = 𝑘 B 𝑁A Constante des gaz parfaits 8,314 462 618 J ⋅ K−1 ⋅ mol

𝐾cd Efficacité lumineuse 683 lm ⋅ W−1

Autres constantes physiques

−1
𝐺 Constante gravitationnelle 6,67430 × 10−11 m3 ⋅ kg ⋅ s−2

𝜖0 Permittivité diélectrique du vide 8,85418781 × 10−12 F ⋅ m−1

𝜇0 Perméabilité magnétique du vide 1,256637062 × 10−6 H ⋅ m−1

𝑚e Masse de l’électron au repos 9,10938370 × 10−31 kg

𝑚p Masse du proton au repos 1,672621923 × 10−27 kg

𝑚n Masse du neutron au repos 1,674927498 × 10−27 kg

Grandeurs physiques

𝑎⃗ Accélération (m.s−2 )

𝐶𝑥 , 𝐶𝑦 , 𝐶𝑧 Coefficients aérodynamiques (sans dimension)

ℰc Energie cinétique (J)

ℰp Energie potentielle (J)

𝑒c Énergie cinétique volumique (J.m−3 )

𝑒p Énergie potentielle volumique (J.m−3 )

⃗⃗⃗ ⃗⃗
𝐹 force (N)

𝑓⃗⃗⃗⃗ force volumique (N.m−3 )

𝑔⃗ Champ de pesanteur ([Link]−1 )

𝐽𝑚 Densité de courant massique (kg.s−1 .m−2 )

𝑀 Masse molaire ([Link]−1 )

𝑚 Masse (kg)

𝑛⋆ Densité de particules (m−3 )

𝒫 Puissance (W)

𝑝 Pression (Pa)

𝐷𝑚 Débit massique (kg.s−1 )

𝐷𝑉 Débit volumique (m3 .s−1 )

𝑅e Nombre de Reynolds (sans dimension)

𝑇 Température (K)

𝑡 Temps (s)

𝑉 Volume (sa mesure en m3 )

𝑣⃗ Vitesse (m.s−1 )

𝑊 Travail (J)

𝛾 Tension superficielle (J.m−2 )

𝜆 Coefficient de perte de charge (sans dimension)

𝜂 Viscosité (Pa.s)

𝜌 Masse volumique (kg.m−3 )

𝜎 Contrainte (N.m−2 )

𝜔, Ω Vitesse angulaire (rad.s−1 )

Symboles mathématiques
def
= Relation de définition

∼ Égal en ordre de grandeur

𝐴≫𝐵 𝐴 très grand devant 𝐵

𝐴≪𝐵 𝐴 très petit devant 𝐵

⟨𝑔⟩ Moyenne d’une grandeur 𝑔

d𝑓
d𝑡 Dérivée première par rapport au temps
d𝑛 𝑓
d𝑡𝑛 Dérivée n-ième par rapport au temps

(⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥 , ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑦 , ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 ) Base cartésienne

(𝑟, 𝜃, 𝑧) Coordonnées cylindriques

(⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑟 , ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝜃 , ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 ) Base cylindrique

(𝑟, 𝜃, 𝜑) Coordonnées sphériques

(⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑟 , ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝜃 , ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝜑 ) Base sphérique
𝐴𝑧 ⃗⃗⃗ ⃗⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗
Composante suivant l’axe (O𝑧) : 𝐴𝑧 = 𝐴 𝑢𝑧

∫𝒟 Intégration sur un domaine 𝒟

⃗⃗⃗ ⃗⃗
∫C 𝐴(M) ⋅ dℓ ⃗ ⃗⃗⃗ ⃗⃗ le long du circuit C
Circulation de 𝐴

⃗⃗⃗ ⃗⃗
∬S 𝐴(M) 𝑛 d𝑆
⋅ ⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗
Flux d’un champ vectoriel 𝐴

∭V 𝑓(M) d𝜏 Intégrale de volume

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑓 ou ⃗⃗⃗⃗⃗⃗ ∇𝑓 Gradient d’un champ scalaire
⃗⃗⃗ ⃗⃗ ou ⃗⃗⃗⃗⃗⃗
div𝐴 ∇⋅𝐴 ⃗⃗⃗ ⃗⃗ Divergence d’un champ vectoriel
⃗⃗⃗ ⃗⃗ ou ⃗⃗⃗⃗⃗⃗
rot𝐴 ∇∧𝐴 ⃗⃗⃗ ⃗⃗ Rotationnel d’un champ vectoriel

Δ𝑓 = ∇2 𝑓 Laplacien scalaire
2025

[Link]/mecanique_des_fluides