COURS DE PHYSIQUE

MÉCANIQUE DES FLUIDES

JIMMY ROUSSEL

2023

[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 © 2023 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 – 30 exercices et problèmes corrigés.
▶ Dynamique des fluides - Fluides parfaits et newtoniens – 40 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 CINÉMATIQUE 1
1.1 Le modèle continu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Description d’un fluide en écoulement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3 Conservation de la masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4 Caractéristiques d’un écoulement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

2 DYNAMIQUE DES FLUIDES PARFAITS 11

2.1 Bilan des forces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Relation fondamentale de la dynamique des fluides parfaits . . . . . . . . . . . . . . . . . 13
2.3 Statique des fluides . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.4 Écoulements permanents et incompressibles . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

3 FLUIDES VISQUEUX 23
3.1 Notion de viscosité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.2 Dynamique d’un écoulement visqueux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.3 Pertes de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.4 Traînée et portance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35

4 TENSION SUPERFICIELLE 39
4.1 Énergie de surface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.2 Forces capillaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.3 Mouillage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45

ANNEXES 49
A OPÉRATEURS DIFFÉRENTIELS 51
A.1 L’opérateur gradient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
A.2 L’opérateur divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
A.3 L’opérateur rotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
A.4 L’opérateur laplacien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
A.5 Accélération d’une particule de fluide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56

Pour en savoir plus 57

Grandeurs physiques et symboles mathématiques 59

Table des figures

1.1 Les trois états de la matière . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1

1.2 Les trois échelles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3 Ligne d’écoulement. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.4 Calcul du débit massique. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.5 Calcul du débit massique sortant d’une surface fermée. . . . . . . . . . . . . . . . . . . . . . . . 6
1.6 Particule de fluide dans un écoulement uniforme . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.7 Particule de fluide dans un écoulement radial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.8 Tube de courant. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.9 Particule de fluide dans un écoulement rotationnel . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.10 Quelques écoulement réels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

2.1 Notion de force pressante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

2.2 Mesures absolue et relative de pression. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.3 Nivellement barométrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.4 Cylindre en rotation. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.5 Calcul de la poussée d’Archimède sur un cube. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.6 Cas d’un volume quelconque. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.7 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. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.8 Vidange d’un récipient. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.9 Effet Venturi (©Wikimedia Commons) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.10 Tube de Pitot. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.1 Expérience de Couette. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

3.2 Contraintes internes dans le fluide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24

3.3 Mesures de la viscosité de polysaccharides . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.4 Bilan des forces de viscosité sur un élément de fluide. . . . . . . . . . . . . . . . . . . . . . . . . 26
3.5 Écoulement de Poiseuille. Position du problème. . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.6 Profil de vitesse. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.7 Diagramme de Moody. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
3.8 Ecoulement stationnaire dans un système de conduites. . . . . . . . . . . . . . . . . . . . . . . 34
3.9 Traînée sur un obstacle sphérique immobile. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.10 Portance et traînée. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.1 Interactions au sein d’un fluide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39

4.2 Coalescence de deux gouttes (©CNRS) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40

4.3 Production de particules dans la première chambre à bulles . . . . . . . . . . . . . . . . . . . . 41
4.4 Mise en évidence des forces capillaires. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.5 Densité linéique de forces superficielles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.6 Méthode de l’anneau. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
4.7 Existence d’une surpression dans une bulle. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
4.8 Lame de savon formant une caténoide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.9 Définition de l’angle de contact. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.10 Différents équilibres d’un liquide au contact d’un substrat solide. . . . . . . . . . . . . . . . . . 46
4.11 Expérience et schéma. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46

Liste des tableaux

1.1 Eau liquide à 25°C, 1 atm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1

1.2 Vapeur d’eau à 100°C, 1 atm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2

2.1 Autres unités de pression couramment utilisées. . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

3.1 Quelques valeurs de viscosité. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24

3.2 Quelques valeurs de temps de relaxation visco-élastique. . . . . . . . . . . . . . . . . . . . . . . 26
3.3 Quelques ordres de grandeur du nombre de Reynolds. . . . . . . . . . . . . . . . . . . . . . . . 28
3.4 Aérodynamique de quelques véhicules. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

4.1 Quelques valeurs de tension superficielle. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43

A.1 Expressions de l’opérateur gradient dans différents systèmes de coordonnées. . . . . . . . . 51

A.2 Expressions de la divergence dans différents systèmes de coordonnées. . . . . . . . . . . . . . 53
A.3 Expressions du rotationnel dans différents systèmes de coordonnées . . . . . . . . . . . . . . 54
A.4 Expressions du laplacien dans différents systèmes de coordonnées. . . . . . . . . . . . . . . . 55
 CINÉMATIQUE 1
Dans ce cours, nous étudions le fluide et son écoulement indépen- 1.1 Le modèle continu . . . . 1
damment des forces responsables de cet écoulement. L’état fluide . . . . . . . . 1
Le modèle continu . . . . 2
Version en ligne 1.2 Fluide en écoulement . . 3
Ligne d’écoulement . . . 3
[Link]
Notion de trajectoire . . . 4
[Link] Dérivée particulaire . . . 4
Accélération d’une particule 5
1.3 Conservation de la masse 6
1.1 Le modèle continu Densité de courant de
masse . . . . . . . . . . . 6
Équation de continuité . . 6
L’état fluide 1.4 Différents écoulements . 7
Interprétation de div 𝑣 ⃗ . . 7
Le terme fluide désigne un comportement qui s’oppose au comporte- Ecoulement incompressible 9
ment élastique ou plastique associé aux solides. Par définition, on dit Interprétation de rot 𝑣 ⃗ . . 9
que la matière est fluide lorsqu’elle se déforme aussi longtemps que
lui sont appliquées des contraintes tangentielles. En termes simples
on peut dire qu’un fluide coule quand un solide se déforme1 1 : Nous verrons plus tard qu’il existe
des milieux qui présentent un compor-
Fondamentalement, le comportement fluide est lié, au niveau mo- tement fluide lorsqu’on les observe
léculaire, à l’absence d’ordre à longue portée2 (contrairement aux sur de longues échelles de temps alors
cristaux) et à l’existence d’un chaos moléculaire (contrairement aux qu’ils ont un comportement élastique
voire plastique sur de petites échelles
solides). Ces propriétés se retrouvent notamment chez les gaz et les de temps.
liquides.
2 : Certains systèmes présentent un
ordre à longue portée suivant une
Liquide Gaz seule direction. Ils ont alors un carac-
Solide
tère cristallin selon cette direction et
fluide selon les autres. On les désigne
par le terme cristaux liquides.

Entropie FIG. 1.1 : Les trois états de la matière

Les liquides – Dans un liquide, les interactions (l’interaction de Van

der Waals, la liaison hydrogène, l’interaction électrostatique dans une
solution électrolytique etc ...) jouent un rôle clé. L’interaction est telle
que les molécules sont quasi en contact ce qui explique le carac-
tère quasi-incompressible des liquides : les liquides présentent un |Δ𝜌|
Δ𝑇 = 10 K = 0, 2%
𝜌
volume propre. Les variations du volume 𝑉 ou de la masse volumique
𝜌 avec la pression et la température se mesurent à l’aide du coeffi- |Δ𝜌|
Δ𝑝 = 1 bar = 0, 02%
𝜌
cient de dilatation 𝛼 et du coefficient de compressibilité 𝜒𝑇 :
TAB. 1.1 : Eau liquide à 25°C, 1 atm.
1 𝜕𝜌 1 𝜕𝜌
𝜒𝑇 = ∣ et 𝛼= − ∣
𝜌 𝜕𝑝 𝑇 𝜌 𝜕𝑇 𝑝

Pour l’eau, par exemple, la compressibilité vaut 𝜒𝑇 ≈ 4, 4.10−10 Pa−1

à 20 °C. Cela signifie qu’il faut augmenter la pression de 227 bars pour
 2 1 CINÉMATIQUE

voir la masse volumique augmenter de 1%. Les liquides ont égale-

ment un coefficient de dilatation très faible.

Le gaz – Dans un gaz, les particules interagissent peu, l’énergie est

avant tout cinétique. Les distances inter-atomiques sont grandes ce
|Δ𝜌|
Δ𝑇 = 10 K = 3% qui explique qu’à l’inverse des liquides, les gaz sont très compres-
𝜌
sibles. Pour un gaz, dans les conditions de pression et de température
|Δ𝜌|
Δ𝑝 = 1 bar = 100% raisonnables et loin de tout point critique, le modèle du gaz parfait
𝜌
est tout à fait suffisant.
TAB. 1.2 : Vapeur d’eau à 100°C, 1 atm.

Approximation courantes

▶ Cas des liquides : 𝜌 ≈ constante

▶ Cas des gaz : modèle du Gaz Parfait

𝑀𝑝
𝑝𝑉 = 𝑛 𝑅 𝑇 ⇒ 𝜌= avec 𝑅 = 8, 315 J.K−1 .mol−1
𝑅𝑇

Le modèle continu

Plusieurs approches sont possibles pour décrire un fluide :

▶ L’approche « Dynamique Moléculaire » : On peut chercher le

comportement de N molécules en résolvant de façon numé-
rique les équations de la mécanique du point appliquées aux
N corps. Bien entendu, la limitation des ordinateurs, impose 𝑁
petit. À l’heure actuelle le record est de l’ordre de 100 milliards
d’atomes.
▶ L’approche « milieu continu » : Lorsque le libre parcours moyen
ℓ est très petit devant l’échelle macroscopique, on choisit de
décrire le fluide à une échelle intermédiaire entre l’échelle ato-
mique et macroscopique : l’échelle mésoscopique.
▶ L’approche statistique : Lorsque le libre parcours moyen des
molécules est du même ordre de grandeur que l’échelle macro-
scopique on utilise les équations de physique statistique (équa-
tion de Boltzmann) pour décrire le fluide.

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 𝐿 est le diamètre du

tuyau quand on étudie l’écoulement dans un tuyau.
2. L’échelle des collisions ℓ ≪ 𝐿. ℓ est 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
 1.2 Description d’un fluide en écoulement 3

𝐿≫𝑎≫ℓ

• •
• • • •
•
•
• • • • • •
• •
• •• •
M •
•
• • •
• • • •
•
•
𝑣(M,
⃗⃗⃗⃗ 𝑡) • • •
•
•
• • • •
• •
• • •• •
• • •

FIG. 1.2 : Les trois échelles.

définir des grandeurs moyennes locales qui, elles, vont évoluer de

façon continue : la masse volumique locale 𝜌(M, 𝑡), la vitesse locale
𝑣(M,
⃗ 𝑡)... On donne à ce « sac de molécules » le nom de particule de Remarque : la vitesse 𝑣 en mécanique
fluide qu’il ne faut pas confondre avec la notion de molécule. des fluides désigne la norme du vec-
teur vitesse d’une particule de fluide.
En conséquence on peut avoir 𝑣(M) =
Conditions de validité du modèle 0 bien que la vitesse moyenne d’une
molécule soit non nulle :
Un milieu peut être considéré continu si le libre parcours moyen
𝑣𝑖 ∥ ≠ ‖⃗⃗⃗⃗⃗⃗⃗
𝑣 = ∥⃗⃗⃗⃗⃗⃗⃗ 𝑣𝑖 ‖
ℓ des molécules est petit devant la taille caractéristique 𝐿 du sys-
tème étudié. On définit le nombre de Knudsen
ℓ
𝐾𝑛 = ≪1
𝐿
Lorsque 𝐾𝑛 n’est pas petit devant 1, le modèle continu n’est plus
adapté.
Par exemple, dans la haute atmosphère, à l’altitude de 100 km, on a
ℓ ≈ 0, 3 m. Pour calculer l’écoulement autour d’un véhicule spatial
à cette altitude, le modèle continu ne conviendra pas.

1.2 Description d’un fluide en écoulement

Deux approches différentes existent. Le point de vue de Lagrange

consiste à s’intéresser à la trajectoire des particules de fluide. Celle
d’Euler se concentre sur l’évolution des propriétés du fluide en diffé-
rents points et au cours du temps.

Ligne d’écoulement

Adoptons l’approche d’Euler et supposons que l’on connaisse à chaque

instant 𝑡 le vecteur vitesse d’une particule de fluide située en M. Le ligne d’écoulement
vecteur vitesse 𝑣(M,
⃗ 𝑡) désigne alors un champ vectoriel.

Notion de ligne d’écoulement – Par définition, une ligne de courant

ou ligne d’écoulement, est une ligne de champ du vecteur vitesse, 𝑣(M,
⃗⃗⃗⃗ 𝑡)
c’est-à-dire une courbe 𝒞 telle qu’à un instant 𝑡 fixé, pour tout point
M ∈ 𝒞, 𝑣(M,
⃗ 𝑡) est tangente à 𝒞 en M. Lorsque le champ de vitesse M

ne dépend pas du temps, les lignes d’écoulement n’évoluent pas au

cours du temps : on dit que le régime d’écoulement est stationnaire
ou permanent.
FIG. 1.3 : Ligne d’écoulement.
4 1 CINÉMATIQUE

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
𝑢𝑥 + 𝑘𝑦 ⃗⃗⃗⃗⃗⃗⃗⃗

▶ Le régime est-il permanent ?

▶ Quelle est l’équation des lignes d’écoulement ?
▶ La vitesse est-elle constante le long d’une ligne de courant ?

Visualisation des 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 le sens de la vi-
tesse en différents points ce qui permet de reconstituer la carte du
champ de vitesse.

Notion de trajectoire

Dans la description de Lagrange, on s’intéresse à l’histoire de chaque

particule de fluide. Considérons une particule de fluide 𝒫 située en
(𝑥(𝑡), 𝑦(𝑡), 𝑧(𝑡)) à l’instant 𝑡. Par définition la trajectoire est la courbe
paramétrique 𝒞 d’équation

⎧𝑥(𝑡) ⎧𝑥̇ = 𝑣𝑥 (𝑥(𝑡), 𝑦(𝑡), 𝑧(𝑡), 𝑡)

{ {
𝒞 ⎨ 𝑦(𝑡) tel que ⎨
𝑦 ̇ = 𝑣𝑦 (𝑥(𝑡), 𝑦(𝑡), 𝑧(𝑡), 𝑡)
{ 𝑧(𝑡) { 𝑧 ̇ = 𝑣 (𝑥(𝑡), 𝑦(𝑡), 𝑧(𝑡), 𝑡)
⎩ ⎩ 𝑧

La trajectoire retrace l’histoire d’une particule alors que la ligne d’écou-

lement est un « instantanée » du champ de vitesse. De ce fait, ces
deux notions sont différentes. Par contre, lorsque que le régime d’écou-
lement est stationnaire, une particule suit nécessairement la ligne
d’écoulement sur laquelle elle se trouve puisque celle-ci est fixe.

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

Visualisation d’une trajectoire – On utilise des traceurs (colorants ou

fumées) et l’on prend une photo avec un long temps de pose.

Dérivée particulaire

Considérons une grandeur physique locale 𝐺(M, 𝑡) attachée à une

particule de fluide située en M à l’instant 𝑡. On peut penser à la tem-
pérature, la pression, la densité etc. Cherchons à calculer le taux de
 1.2 Description d’un fluide en écoulement 5

variation de cette grandeur lorsque l’on suit la particule. On appelle

cette grandeur la dérivée particulaire et on la note D𝐺
D𝑡 .

D𝐺 𝐺(𝑥 + 𝑣𝑥 𝛿𝑡, 𝑦 + 𝑣𝑦 𝛿𝑡, 𝑧 + 𝑣𝑧 𝛿𝑡, 𝑡 + 𝛿𝑡) − 𝐺(𝑥, 𝑦, 𝑧, 𝑡)

= lim
D𝑡 𝛿𝑡→0 𝛿𝑡
𝐺(𝑥, 𝑦, 𝑧, 𝑡) + 𝑣𝑥 𝜕𝐺 𝜕𝐺 𝜕𝐺
𝜕𝑥 𝛿𝑡 + 𝑣𝑦 𝜕𝑦 𝛿𝑡 + 𝑣𝑧 𝜕𝑧 𝛿𝑡 +
𝜕𝐺
𝜕𝑡 𝛿𝑡 − 𝐺(𝑥, 𝑦, 𝑧, 𝑡)
=
𝛿𝑡
ce qui donne la formule que l’on retiendra :

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

Accélération d’une particule

Calculons l’accélération d’une particule de fluide à partir du champ

de vitesse eulérien 𝑣(M,
⃗ 𝑡). L’accélération est le taux de variation du
champ de vitesse en suivant une particule de fluide. On a donc :

D𝑣𝑥 D𝑣𝑦 D𝑣
𝑎⃗ = 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗
𝑥+ ⃗⃗⃗⃗⃗⃗⃗⃗𝑦 + 𝑧 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢 𝑢
D𝑡 D𝑡 D𝑡 𝑧
ce qui donne
D𝑣𝑥 𝜕𝑣
𝑎𝑥 = = 𝑥 + (𝑣.⃗ ⃗⃗⃗⃗⃗⃗
∇)𝑣𝑥
D𝑡 𝜕𝑡
D𝑣𝑦 𝜕𝑣𝑦
𝑎𝑦 = = + (𝑣.⃗ ⃗⃗⃗⃗⃗⃗
∇)𝑣𝑦
D𝑡 𝜕𝑡
D𝑣𝑧 𝜕𝑣
𝑎𝑧 = = 𝑧 + (𝑣.⃗ ⃗⃗⃗⃗⃗⃗
∇)𝑣𝑧
D𝑡 𝜕𝑡
On pourra retenir le résultat sous forme compacte :

D𝑣 ⃗ 𝜕𝑣⃗
𝑎⃗ = = + (𝑣.⃗ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗ ♡ (1.2)
D𝑡 𝜕𝑡

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 bidimensionnel dont le champ de

vitesse est de la forme :

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

Soit P une particule de fluide située en 𝑥 = 1 et 𝑦 = 1 à l’instant 𝑡0 . Cal-

culer la vitesse scalaire, le vecteur accélération de P ainsi que le rayon de
courbure de sa trajectoire à l’instant 𝑡0 .
 6 1 CINÉMATIQUE

1.3 Conservation de la masse

Vecteur densité de courant de masse

Un écoulement est un phénomène de transport puisqu’il s’agit d’un

transfert de masse. C’est pourquoi il est naturel d’introduire la notion
⃗⃗⃗⃗ d𝑡
𝑣 de vecteur densité de courant de masse.

Pour cela, cherchons à exprimer la masse qui traverse une section

(S) lors d’un écoulement. Considérons une section infinitésimale d𝑆
2
autour d’un point M et calculons la masse d 𝑚 de fluide traversant
M
• 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 𝑚 = 𝜌(M, 𝑡)𝑣(M,
⃗ 𝑛 d𝑡 d𝑆
𝑡) ⋅ ⃗⃗⃗⃗
FIG. 1.4 : Calcul du débit massique.
où ⃗⃗⃗⃗
𝑛 est le vecteur normal à la section d𝑆. En sommant toutes les
contributions on obtient

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

On en déduit le flux de matière ou débit massique

def d𝑚
𝑄m = =∬ 𝜌(M, 𝑡)𝑣(M,
⃗ 𝑡) ⋅ ⃗⃗⃗⃗
𝑛d𝑆 en kg/s
d𝑡 M∈(𝑆)

Le débit massique est donc, au sens mathématique, le flux du vecteur

⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐽m = 𝜌𝑣 ⃗ :
Remarque : le débit volumique 𝑄V me-
sure le volume de fluide qui traverse la
surface (S) par unité de temps (unité : 𝑄m = ∬ ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐽m (M, 𝑡) ⋅ ⃗⃗⃗⃗
𝑛d𝑆 où ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐽m (M, 𝑡) = 𝜌(M, 𝑡) 𝑣(M,
⃗ 𝑡)
m3 .s−1 ) : M∈(𝑆)

2
1d 𝑚 Le vecteur ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐽m désigne le vecteur densité de courant de masse.
𝑄𝑉 = ∬ =∬ ⃗⃗⃗⃗ 𝑛 dS
𝑣.⃗⃗⃗⃗⃗
M∈(𝑆) 𝜌 d𝑡 M∈(𝑆)

É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 re-
𝑛⃗
pose sur une loi de conservation, celle de la masse.
𝑣⃗
d𝑆
Prenons un volume de contrôle fixe (𝑉 ) dans un fluide, délimité par
une surface fictive (𝑆). Soit 𝑚(𝑡) la masse contenue à l’intérieur de la
surface fermée à l’instant 𝑡. Par définition de la masse volumique,

𝑚(𝑡) = ∭ 𝜌(M, 𝑡) d𝜏
M∈(𝑉 )
(𝑆) surface fermée

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

FIG. 1.5 : Calcul du débit massique sor-
tant d’une surface fermée.
d𝑚(𝑡) ext
= −∬ 𝜌𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆
d𝑡 M∈(𝑆)
 1.4 Caractéristiques d’un écoulement 7

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

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

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.
ext
∬ ⃗⃗⃗ ⃗⃗
𝐴(M) ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 =∭ ⃗⃗⃗ ⃗⃗
div𝐴(M) d𝜏 ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗
avec div𝐴 ∇⋅𝐴 ⃗⃗⃗ ⃗⃗
M∈(𝑆) M∈𝑉

D’après le théorème de la divergence on obtient :

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

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

continuité

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

1.4 Caractéristiques d’un écoulement

Interprétation de la divergence de la vitesse

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
 8 1 CINÉMATIQUE

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

− =− = =
𝜌 D𝑡 𝜌 D𝑡 𝜌𝛿𝜏 2 D𝑡 𝛿𝜏 D𝑡

Finalement, on obtient

1 𝐷𝛿𝜏
div𝑣 ⃗ = ♡ (1.4)
𝛿𝜏 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. 1.6.

FIG. 1.6 : Particule de fluide dans un écoulement uniforme, capturé à différents instants.
simulation en ligne

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. 1.7, ce
qui prouve que l’écoulement est à divergence positive.

FIG. 1.7 : Particule de fluide dans un écoulement radial, capturé à différents instants.
simulation en ligne
 1.4 Caractéristiques d’un écoulement 9

Ecoulement incompressible
Définition
Un fluide est en écoulement incompressible quand les particules
de fluide ont un volume qui reste constant au cours de l’écou-
lement. Elles se déforment donc sans variation de masse volu-
mique :
D𝜌
=0
D𝑡

Par conséquent, d’après la relation (1.4), un fluide en écoulement in-

compressible vérifie la relation

div𝑣 ⃗ = 0 ⇔ ∬ 𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 ext = 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

∬ 𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 1 = ∬ 𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 2 soit 𝑄V1 = 𝑄V2 • ⃗⃗⃗⃗⃗⃗⃗⃗
𝑣2
(𝑆1 ) (𝑆2 )
(𝑆1 )

Le débit volumique se conserve le long d’un tube de courant. Si l’on (𝑆2 )

définit la vitesse moyenne dans le section S par :
FIG. 1.8 : Tube de courant.
𝑄def
𝑣 = V ♡ (1.5)
𝑆

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.

Interprétation du rotationnel de la vitesse

Exemple – Commençons par l’étude de l’exemple suivant. Soit un

écoulement bidimensionnel dont le champ de vitesse s’écrit

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

𝜕(−𝑘 𝑦)
1. Tout d’abord, on constate que div 𝑣 ⃗ = 𝜕𝑥 + 𝜕(𝑘 𝑥)
𝜕𝑦 = 0. L’écou-
lement est donc incompressible.
2. Le rotationnel de la vitesse vaut
𝜕(𝑘 𝑥) 𝜕(−𝑘 𝑦)
rot𝑣 ⃗ = ⃗⃗⃗⃗⃗⃗
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ ∇ ∧ 𝑣⃗ = ( − ) ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 = 2𝑘 ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧
𝜕𝑥 𝜕𝑦

3. Plaçons une particule carrée d’arête 2𝑎 en O. Quel est son mou-

vement ? La FIG. 1.9 montre que les sommets sont animés d’une
vitesse orthoradiale de sorte que la particule tourne autour de
 10 1 CINÉMATIQUE

√
O. La vitesse d’un sommet vaut 𝑣 = 2𝑘 𝑎 d’où une vitesse an-
gulaire de rotation 𝜔 = 𝑘.

FIG. 1.9 : Particule de fluide dans un écoulement rotationnel, capturé à différents instants (Simulation ©[Link])
simulation en ligne

Vecteur tourbillon – L’exemple précédent montre que pour un écou-

lement rotationnel, les particules de fluide tournent à une vitesse
angulaire égal à la moitié de la valeur du rotationnel. De manière
générale, on définit le vecteur tourbillon

⃗⃗⃗⃗⃗ 1 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
Ω = rot 𝑣⃗
2

Lorsque ⃗⃗⃗⃗⃗
Ω ≠ 0, l’écoulement est tourbillonnaire ce qui se traduit par
l’existence d’un mouvement de rotation des particules lors de l’écou-
lement.
Ecoulement potentiel – Lorsque ⃗⃗⃗⃗⃗
Ω = 0,⃗ le champ de vitesse est né-
cessairement un gradient.

rot𝑣 ⃗ = 0⃗
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ ⇒ 𝑣 ⃗ = ⃗⃗⃗⃗⃗⃗
∇𝜑

où 𝜑(𝑀 , 𝑡) désigne le potentiel des vitesses. Dans ce cas on parle

d’écoulement irrotationnel ou potentiel.

FIG. 1.10 : Quelques écoulement réels. À gauche, écoulement d’un fluide à la vitesse d’un millimètre par seconde autour d’un
obstacle cylindrique. L’écoulement est permanent, incompressible et irrotationnel caractérisé par l’équation △𝜑 = 0. À droite,
allée tourbillonnaire de Von Karman produit par un écoulement suffisamment rapide autour d’un obstacle cylindrique.
 DYNAMIQUE DES FLUIDES
PARFAITS 2
La dynamique des fluides relie l’écoulement d’un fluide aux actions 2.1 Bilan des forces . . . . . . 11
qui lui sont appliquées. Ce cours se limite aux écoulements pour les- Forces extérieures . . . . 11
quels les couches de fluide glissent les unes sur les autres sans dissi- Forces de pression . . . . 12
pation de chaleur. On parle alors de fluide parfait. On détaille particu- Force volumique pressante 12
lièrement le cas des fluides au repos ainsi que celui des écoulements 2.2 Relation fondamental . . 13
incompressibles et stationnaire. Équation d’Euler . . . . . 13
Résolution de l’eq. d’Euler 14
2.3 Statique des fluides . . . 15
Version en ligne
Cas des Liquide . . . . . . 15
[Link] Cas des gaz . . . . . . . . 16
[Link] Liquide en rotation . . . . 17
Poussée d’Archimède . . 17
2.4 Écoulements permanents
et incompressibles . . . . 18
2.1 Bilan des forces Théorème de Bernoulli . 18
Formule de Torricelli . . . 20
Effet Venturi . . . . . . . 20
Au sein d’un fluide parfait, on distingue deux types de forces : Tube de Pitot . . . . . . . 21
▶ Les forces dont l’origine est extérieure au fluide comme la pe-
santeur ;
▶ Les forces de contact entre particules de fluide que l’on appelle
forces internes.

Forces extérieures

Nous appelons forces extérieures les actions à distance qui agissent

sur toutes les molécules du fluide de sorte qu’elles sont proportion-
nelles au volume de fluide considéré. Pour une particule de fluide de
volume d𝜏 , la force extérieure se met sous la forme :

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 ext = ⃗⃗𝑓⃗ ⃗ ⃗ext d𝜏
v
[N] = [N.m−3 ] [m3 ]

La quantité 𝑓⃗⃗⃗ ⃗v⃗ext s’appelle force volumique extérieure. La résultante

des forces extérieures sur un volume (𝑉 ) se calcule ainsi :

⃗⃗⃗ ⃗⃗ ext = ∭
𝐹 ⃗⃗𝑓⃗ ⃗ ⃗ext d𝜏
v
M∈(𝑉 )

La pesanteur – Pour un fluide plongé dans un champ de pesanteur,

on a ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 ext = d𝑚⃗⃗𝑔⃗ ⃗ = 𝜌⃗⃗𝑔⃗ ⃗d𝜏 d’où

𝑓⃗⃗⃗ ⃗v⃗ext = 𝜌𝑔 ⃗ ♡ (2.1)

La force électromagnétique – Pour un plasma baignant dans un champ

électromagnétique on a

⃗⃗𝑓⃗ ⃗ ⃗ext = 𝜌 (𝐸 ⃗⃗⃗ ⃗⃗

⃗⃗⃗⃗⃗⃗ + 𝑣 ⃗ ∧ 𝐵) ♡ (2.2)
v e
 12 2 DYNAMIQUE DES FLUIDES PARFAITS

où 𝜌e désigne la densité volumique de charges.

Les forces d’inertie – Si le référentiel d’étude n’est pas galiléen, il faut

tenir compte des forces d’inertie. Une masse d𝑚 de fluide subit une
force d’inertie
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑎𝑒 − 2d𝑚⃗⃗⃗⃗⃗
d𝐹 ext = −d𝑚 ⃗⃗⃗⃗⃗⃗⃗ Ω ∧ 𝑣⃗
d’où
ext
𝑓⃗⃗⃗ ⃗v⃗ 𝑎𝑒 + 2⃗⃗⃗⃗⃗
= −𝜌(⃗⃗⃗⃗⃗⃗⃗ Ω ∧ 𝑣)⃗ ♡ (2.3)

Forces de pression

Dans un fluide au repos, une particule de fluide est soumise de la

part du fluide environnant à des actions dites internes au fluide. On
admettra que :
Fluide
▶ Ces forces sont proportionnelles à la surface sur laquelle elles
ext
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 s’exercent. En effet, microscopiquement, les actions inter-moléculaires
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 int sont de très courte portée (collision ou interaction de Van Der
Walls) de sorte qu’elles concernent uniquement les molécules
à la surface de la particule de fluide.
▶ Ces forces sont orientées perpendiculairement à la surface sur
Particule de fluide
laquelle elles s’exercent.
FIG. 2.1 : Notion de force pressante. ▶ Le coefficient de proportionnalité s’appelle la pression.

Mathématiquement, on écrira

int ext
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 = −𝑝(M) ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆
♡ (2.4)
−2 2
[N] = [N.m ][m ]

ext
où ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 représente le vecteur élément de surface dirigé vers l’exté-
rieur de la particule de fluide par convention.

1 : Blaise Pascal (1623–1662) : Mathé- Unités – Dans le Système International d’Unités, la pression s’exprime
maticien, physicien et philosophe fran- en pascal (symbole Pa), en hommage à Blaise Pascal1 .
çais, né à Clermont-Ferrand. Pascal
contribua au développement des ma-
thématiques (probabilités, invention 1 Pa = 1 N.m−2 ♡ (2.5)
du premier calculateur) de la philoso-
phie (« Les pensées ») et de l’hydrau- Exercice – À quelle masse répartie sur 1 cm2 correspond 1 bar si l’on prend
lique (principe de Pascal). Il met fin au 𝑔 = 10 m.s−2 ?
dogme aristotélicien selon lequel « la
nature a horreur du vide » . En effet, il Rép. 1 kg.
montre, avec la collaboration de son
beau frère Perrier, que l’ascension du
mercure dans l’expérience de Toricelli Il existe également d’autres unités encore très utilisées (TAB. 2.1).
est due à la pression atmosphérique.

Force volumique pressante

Dans un fluide parfait, les forces internes sont normales : elles se

résument aux forces de pression. Insistons sur le fait que les forces
de pression sont des actions de contact (interactions de courte por-
tée) et donc surfaciques. Ceci étant dit, montrons que la résultante
 2.2 Relation fondamentale de la dynamique des fluides parfaits 13

des forces de pression sur une particule de fluide peut s’exprimer à

l’aide d’une force volumique.

Dans un fluide, délimitons un volume (V) à l’aide d’une surface fer-

mée (𝑆) et calculons la résultante des forces de pression qu’exerce
le fluide extérieur sur ce volume. En un point M de la surface (𝑆), la Unités Équivalence
force de pression vaut en pascal
1 bar 105 Pa
1 atmosphère 1,013 ⋅ 105 Pa
⃗⃗⃗ ⃗⃗ = ∬
𝐹 −𝑝(M) ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 ext 1 torr = 1 mm Hg 133,3 Pa
M∈(𝑆)
Pressions rencontrées dans la Nature
Munissons nous d’un repère cartésien et exprimons la force suivant Centre du Soleil 1016 Pa
l’axe Ox : Centre de la Terre 1011 Pa
Record Haute pression en 1, 5 ⋅ 1010 Pa
𝐹𝑥 = ∬ −𝑝(M)⃗⃗𝑢⃗⃗𝑥 ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 ext
labo
M∈(𝑆)
Abysse des océans 108 Pa
Le théorème de Green-Ostrogradsky donne alors : Talons aiguilles 107 Pa
Atmosphère de la terre 105 Pa
𝜕𝑝 Trompe à eau 1000 Pa
𝐹𝑥 = ∭ −div (𝑝(M)⃗⃗𝑢⃗⃗𝑥 ) d𝜏 = ∭ − d𝜏 Ultra vide en labo 10−12 Pa
M∈(𝑉 ) M∈(𝑉 ) 𝜕𝑥
TAB. 2.1 : Autres unités de pression cou-
Ainsi la résultante des forces de pression s’écrit : ramment utilisées et quelques ordres
de grandeur.
𝜕𝑝 𝜕𝑝 𝜕𝑝
⃗⃗⃗ ⃗⃗ = ∭
𝐹 −( 𝑢
⃗⃗⃗⃗𝑥 + 𝑢
⃗⃗⃗⃗𝑦 + ⃗⃗⃗⃗ ) d𝜏 = ∭
𝑢 −⃗⃗⃗⃗⃗⃗
∇𝑝(M) d𝜏
M∈(𝑉 ) 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝑧 M∈(𝑉 )
(2.6)
On peut donc définir une force volumique ⃗⃗⃗⃗⃗⃗⃗ 𝑓𝑝 associée à la résultante
des forces de pression :

⃗⃗⃗⃗⃗⃗⃗
𝑓𝑝 = −⃗⃗⃗⃗⃗⃗
∇𝑝(M) ♡ (2.7)

Cette formule signifie que lorsque la pression est uniforme, la ré-

sultante des forces de pression sur un volume de fluide est nulle.
Lorsque la pression n’est pas uniforme, la résultante des forces de
pression sur une particule de fluide est opposée au gradient de pres-
sion c’est-à-dire orientée vers les valeurs décroissantes de pression :
La force de pression tend donc à déplacer les particules de fluide
vers les basses pressions.

L’équation (2.7) signifie aussi que la résultante des forces de pression

est conservative. On peut définir une énergie potentielle volumique de
pression 𝑒p = 𝑝.

2.2 Relation fondamentale de la dynamique

des fluides parfaits

Équation d’Euler

Supposons un fluide parfait en écoulement dans un référentiel gali- 2 : Si le référentiel n’est pas galiléen, il
léen2 . Appliquons le Principe Fondamentale de la Dynamique à chaque suffit d’ajouter dans le bilan des forces,
les forces volumiques d’inertie.
particule de fluide :
 14 2 DYNAMIQUE DES FLUIDES PARFAITS

▶ Système d’étude : une particule de fluide située en M à l’instant

𝑡 de masse d𝑚 = 𝜌(M, 𝑡) d𝜏
▶ Bilan des forces :

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹 = ⃗⃗𝑓⃗ ⃗v⃗ext d𝜏 + 𝑓⃗⃗⃗𝑝⃗ d𝜏 = (⃗⃗𝑓⃗ ⃗v⃗ext − ⃗⃗⃗⃗⃗⃗
∇𝑝(M, 𝑡)) d𝜏

2nde loi de Newton d𝑚 D⃗⃗

𝑣⃗⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
▶
D𝑡 = d𝐹 d’où l’équation d’Euler

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

Résolution de l’équation d’Euler

L’équation d’Euler est une équation aux dérivées partielles du pre-

mier 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
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 une constante
d’intégration par variable. On déterminera ces constantes d’intégra-
tion par les conditions aux interfaces (fluide1/fluide2 ou fluide/so-
lide).

Conditions aux limites d’un fluide parfait

▶ Condition sur 𝑣(M,

⃗ 𝑡) : à la traversée d’une interface, la com-
posante normale de la vitesse est continue.
▶ Condition sur 𝑝(M, 𝑡) : la pression est continue à la traversée
d’une interface fluide-fluide3
3 : condition valable quand on néglige
les effets de tension superficielle.
 2.3 Statique des fluides 15

2.3 Statique des fluides

La statique des fluides étudie les fluides au repos. Il existe alors un

référentiel dans lequel 𝑣(M,
⃗ 𝑡) = 0⃗ partout et à chaque instant. Dans
ce référentiel, l’équation d’Euler devient

⃗⃗𝑓⃗ ⃗ ⃗ext − ⃗⃗⃗⃗⃗⃗

∇𝑝 = 0⃗ (2.9)
v

Traitons quelques situations particulières.

Liquide dans un champ de pesanteur

Considérons un liquide au repos dans le champ de pesanteur.

▶ Supposons le liquide incompressible ⇒ 𝜌 = Cte

▶ Seule la pesanteur agit sur le fluide ⇒ 𝑓⃗⃗⃗ ⃗v⃗ext = 𝜌 𝑔 ⃗
▶ En projetant l’équation de la statique ( axe z descendant !) on
obtient
𝜕𝑝 𝜕𝑝 𝜕𝑝
=0 = 0 et = 𝜌𝑔 (2.10)
𝜕𝑥 𝜕𝑦 𝜕𝑧
La pression ne dépend donc que de 𝑧.
▶ Finalement, en intégrant la dernière relation

d𝑝
= 𝜌𝑔 ⇒ 𝑝(𝑧) = 𝑝0 + 𝜌 𝑔 𝑧 ♡ (2.11)
d𝑧

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

𝜌𝑔𝑧 représente la pression due au poids de la colonne de li-
quide 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 ?
105
Rép. 𝜌𝑔 ℎ = 105 ⇒ ℎ≃ 1000×9,8
= 10, 2 m

Conséquences –

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

(isobare 𝑝 = 𝑝ext ) aussi.
4 : Élève de Galilée, il découvre en 1644
▶ Mesure absolue de pression à l’aide du baromètre à mercure le principe du baromètre à mercure.
inventé par Evangelista Torricelli4 (FIG. 2.2). On a Le 19 septembre 1648, le beau-frère de
Blaise Pascal, Florin Périer, montre que
𝑝atm = 𝜌𝑔ℎ ⇒ 1 atm = 1, 01325.105 Pa → ℎ ≃ 760 mm de Hg la hauteur du mercure dans le baro-
mètre de Torricelli ne dépend que de
la pression atmosphérique. Il obtient à
Clermont Ferrand (460m) ℎ = 71, 2 cm
et sur le Puy de Dôme (1464m) ℎ =
Avant 1954, 1 atmosphère correspondait à une hauteur de mercure 62, 7 cm 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.

▶ Mesure relative de pression à l’aide d’un manomètre à liquide

en U (FIG. 2.2) : Δ𝑝 = 𝜌𝑔ℎ.
 16 2 DYNAMIQUE DES FLUIDES PARFAITS

Vide Atmosphère à la
pression 𝑝atm

Hgliq Gaz sous

ℎ pression
𝑝 > 𝑝atm
𝑔⃗ ℎ
Liquide
𝑝atm

𝑔⃗

FIG. 2.2 : Mesures absolue et relative de pression. À gauche, expérience de Torricelli (1644). À droite, Manomètre à liquide.

Gaz dans un champ de pesanteur

Les gaz étant compressibles, il faut utiliser l’équation d’état ainsi que
les principes de la thermodynamique pour résoudre le problème.

Modèle d’atmosphère isotherme – Considérons un gaz isotherme

(𝑇 = 𝑇0 ) en équilibre dans le champ de pesanteur terrestre. L’équa-
tion d’état du gaz parfait donne :

𝑀𝑝
𝜌=
𝑅𝑇0

En prenant l’axe 𝑧 ascendant, l’équation (2.10) devient

d𝑝 𝑀𝑝
=− 𝑔
d𝑧 𝑅𝑇0

équation différentielle linéaire d’ordre 1 dont la solution s’écrit

𝑅𝑇0
𝑝(𝑧) = 𝑝0 e−𝑧/𝐻 avec 𝐻 = (2.12)
𝑀𝑔

où 𝑝0 est la pression en 𝑧 = 0. La pression, comme la masse volu-

mique ou la densité de particules, décroît exponentiellement avec
l’altitude 𝑧 : c’est la loi de nivellement barométrique. La quantité 𝐻
𝑝 représente une distance caractéristique. Plus de 99% du gaz se re-
𝑝0 trouve sous l’altitude 𝑧 = 5 𝐻.

Si l’altitude est petite devant 𝐻, la loi (2.12) donne

37%𝑝0 𝛿𝑝 𝛿𝑧 𝛿𝑝
≃− ⇒ ∼ −10−4 pour 𝑧 = 1m
𝑝 𝐻 𝑝
5%𝑝0 𝑧
𝐻 3𝐻 En conclusion, on peut négliger la pesanteur dans un gaz sur une
FIG. 2.3 : Nivellement barométrique : échelle de quelques mètres.
décroissance exponentielle de la pres-
sion avec l’altitude.
Exercice – On considère que la troposphère et la stratosphère forment une
atmosphère de température moyenne 𝑇0 ≃ 250 K. Quelle est la hauteur
correspondant à 99% d’atmosphère en masse ?
Rép. 𝑧 = 4, 6 𝐻 ≃ 34 km.
 2.3 Statique des fluides 17

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.
Cependant, après un régime transitoire qui dépend de la viscosité du Liquide
𝑢𝑧
⃗⃗⃗⃗⃗⃗⃗⃗⃗
liquide, celui-ci tourne de façon solide à la même vitesse angulaire
𝑟 ⊗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝜃
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é ℛ′ ) M(𝑟, 𝜃, 𝑧)
𝑔⃗
munis d’un système de coordonnées cylindriques.

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

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

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

𝑓𝑖𝑐 = −2𝜌⃗⃗𝜔⃗ ⃗ ∧ 𝑣(M/ℛ
⃗ ′
) est nulle
car, 𝑣(M/ℛ
⃗ ′ ⃗
) = 0. Télescope à miroir liquide
▶ La force d’entraînement (ou force centrifuge) vaut ⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑢𝑟 .
𝑓𝑖𝑒 = 𝜌𝑟𝜔2 ⃗⃗⃗⃗⃗⃗⃗⃗
Des équipes de l’Université Laval,
▶ La force volumique de pesanteur vaut ⃗⃗⃗⃗⃗⃗⃗ 𝑓𝑔 = 𝜌 𝑔 ⃗ de l’University of British Columbia
et de l’Institut d’astrophysique de
L’équation de la statique des fluides (2.9) donne donc Paris, ont mis au point un téles-
cope Nommé Large Zénith Téles-
𝑢𝑟 + 𝜌𝑔 ⃗ − ⃗⃗⃗⃗⃗⃗
𝜌𝑟𝜔2 ⃗⃗⃗⃗⃗⃗⃗⃗ ∇𝑝 = 0⃗ cope (LZT) dont le miroir primaire
fait 6 mètres de diamètre. Contrai-
rement aux télescopes conven-
ce qui donne
tionnels dont le miroir est fait de
verre, le LZT a un miroir fait de
⎧ 𝜕𝑝 liquide réfléchissant, du mercure
{ = 𝜌𝑟𝜔2 → 𝑝 = 1/2𝜌𝜔2 𝑟2 + 𝑓(𝑧)
{ 𝜕𝑟 plus précisément qui adopte une
{ 𝜕𝑝 surface parabolique puisque mis
⎨ = 0 ↓ en rotation dans une cuve. Le mi-
{ 𝑟𝜕𝜑
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
𝜌𝑟2 𝜔2 La principale limitation du LZT,
𝑝(𝑟, 𝑧) = − 𝜌𝑔 𝑧 + Cte
2 et des autres miroirs liquides,
est qu’on ne peut le pointer
La surface libre étant une surface isobare, elle obéit à l’équation ailleurs qu’au zénith, sinon le li-
quide s’échappe de la cuvette.
2 Au-dessus du miroir, une étroite
𝜔2
𝑧 = 𝑧0 + bande de ciel défile à la vitesse de
2𝑔 la rotation terrestre. Un astre don-
né met environ une minute à tra-
La surface libre adopte une forme parabolique d’axe de révolution verser le champ de vision du téles-
(𝑂𝑧). cope.

Poussée d’Archimède

Imaginons un solide cubique d’arête 𝑎 immergé dans un liquide au

repos et calculons la résultante des forces de pression ⃗⃗Π.
⃗ ⃗⃗

▶ Tout d’abord, la pression ne dépend que de la profondeur 𝑝(𝑧) =

𝑝0 + 𝜌 𝑔 𝑧 ;
▶ Par symétrie, les forces de pression horizontales se compensent,
contrairement aux forces verticales du fait de la variation de
 18 2 DYNAMIQUE DES FLUIDES PARFAITS
Liquide

Position d’équilibre
⃗⃗⃗ ⃗⃗ est donc suivant O𝑧 :
pression avec la profondeur. Π
⃗⃗⃗⃗⃗⃗haut ⃗⃗⃗ ⃗⃗ = 𝑝(𝑧) 𝑎2 𝑢
𝐹 Π ⃗⃗⃗⃗𝑧 − 𝑝(𝑧 + 𝑎) 𝑎2 𝑢
⃗⃗⃗⃗𝑧
3
⃗⃗⃗⃗⃗⃗droite
𝐹 = −𝜌 𝑔 𝑎 𝑢 ⃗⃗⃗⃗𝑧
⃗⃗⃗⃗⃗⃗gauche
⃗⃗Π⃗ ⃗⃗ = −𝑚𝑑 𝑔 ⃗
𝐹
𝑔⃗
où 𝑚𝑑 désigne la masse de liquide déplacé.
⃗⃗⃗⃗⃗⃗bas
𝐹
On obtient donc une force ascendante égale au poids du volume de
FIG. 2.5 : Calcul de la poussée d’Archi- liquide déplacé.
mède sur un cube.
Ce résultat particulier se généralise sans difficulté. Si l’on considère
un solide de volume 𝑉 immergé dans un fluide quelconque au repos
Isobare
dans un champ de pesanteur, la résultante des forces de pression
s’écrit
⃗⃗⃗⃗⃗⃗
Π ⃗⃗Π⃗ ⃗⃗ = ∬ −𝑝(M) ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 ext = ∭ −⃗⃗⃗⃗⃗⃗
∇𝑝 d𝜏
M∈(𝑆) 𝑀∈(𝑉 )
d𝑆 𝑛⃗
Or, à l’équilibre 𝜌 𝑔 ⃗ = ⃗⃗⃗⃗⃗⃗
∇𝑝(𝑀 ) d’où
𝑔⃗

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

Fluide au repos M∈(𝑉 )

FIG. 2.6 : Cas d’un volume quelconque.

Cette force, opposée à la pesanteur, s’appelle la poussée d’Archi-
mède.

Théorème d’Archimède (250 av.J.C.)

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 ; ascension des ballons sondes ;

courants de convection etc.

2.4 Écoulements permanents et

incompressibles

Théorème de Bernoulli

Ce théorème énoncé en premier par Daniel Bernoulli, est une équa-

tion intégrale de l’équation d’Euler qui exprime la conservation de
l’énergie.
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é. On a donc
▶ Fluide parfait : les champs de vitesse et de pression vérifient
l’équation d’Euler.
▶ Écoulement permanent : 𝜕𝜕𝑡𝑣⃗⃗⃗⃗ = 0.⃗ Les trajectoires s’identifient
aux lignes de courant.
 2.4 Écoulements permanents et incompressibles 19

▶ Écoulement incompressible : le long d’une trajectoire 𝜌 reste

constant.
▶ De plus, nous supposerons que les forces volumiques extérieures
dérivent d’une énergie potentielle ⃗⃗𝑓⃗ ⃗v⃗ext = −⃗⃗⃗⃗⃗⃗
∇𝑒𝑝 . Par exemple,
dans le champ de pesanteur on a 𝑒𝑝 = 𝜌 𝑔 𝑧 (axe ascendant).

L’équation d’Euler devient donc

𝜌 (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗ ∇𝑝 + 𝑓⃗⃗⃗ ⃗v⃗ext

∇) 𝑣 ⃗ = −⃗⃗⃗⃗⃗⃗

or, le terme convectif peut aussi s’écrire (cf. Chapitre A)

2
𝑣
(𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗ = ⃗⃗⃗⃗⃗⃗
∇ + (⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣)⃗ ∧ 𝑣 ⃗
2
d’où l’équation :
2
𝑣
𝜌 (⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣)⃗ ∧ 𝑣 ⃗ = −⃗⃗⃗⃗⃗⃗
∇𝑝 − ⃗⃗⃗⃗⃗⃗
∇𝑒𝑝 − 𝜌⃗⃗⃗⃗⃗⃗
∇
2
Intégrons cette équation le long d’une ligne de courant entre deux
points A et B :
𝐵 𝐵 2 2 𝐵
𝑣 𝑣
∫ 𝜌 [(⃗⃗⃗⃗⃗⃗ ⃗ ⃗⃗⃗⃗⃗⃗⃗⃗
∇ ∧ 𝑣)⃗ ∧ 𝑣]⋅ ∇ (𝑝 + 𝑒𝑝 + 𝜌 )⋅⃗⃗⃗⃗⃗⃗⃗⃗
dℓ = − ∫ ⃗⃗⃗⃗⃗⃗ dℓ = − [𝑝 + 𝑒𝑝 + 𝜌 ]
𝐴 𝐴 2 2 𝐴

Le premier terme est nul puisque ⃗⃗⃗⃗⃗⃗⃗⃗

dℓ ∥ 𝑣.⃗ In fine,

𝑣(𝐴)2 𝑣(𝐵)2
⇒ 𝑝(𝐴) + 𝜌 + 𝑒𝑝 (𝐴) = 𝑝(𝐵) + 𝜌 + 𝑒𝑝 (𝐵)
2 2

Théorème de Bernoulli (1738)

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

fait, la quantité
𝑣2
𝑝 + 𝜌 + 𝑒𝑝 = Cte
2
le long d’une ligne de courant. Par exemple, dans le champ de
pesanteur on obtient

𝑣2
𝑝+𝜌 + 𝜌𝑔𝑧 = Cte
2
FIG. 2.7 : Daniel Bernoulli (1700 - 1782) :
médecin, physicien et mathématicien
Théorème de Bernoulli (1738) suisse. Il publia en 1738 son ouvrage
Hydrodynamica dans lequel il expose
D’une ligne de courant à l’autre, c’est la valeur de la constante qui le théorème fondamental de la méca-
change. nique des fluides qui porte son nom.

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.
 20 2 DYNAMIQUE DES FLUIDES PARFAITS

Formule de Torricelli
Liquide
section 𝑆
A• 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
avec la hauteur de liquide ℎ. On fera les hypothèses suivantes :

▶ 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 𝑠
•
B 𝑣(𝐵)
⃗
On cherche à calculer la vitesse d’écoulement 𝑣 à la sortie du trou.
L’application du théorème de Bernoulli sur une ligne de courant donne :
FIG. 2.8 : Vidange d’un récipient.
1 1
𝑝atm + 𝜌𝑔ℎ + 𝜌𝑣2 (A, 𝑡) = 𝑝atm + 𝜌𝑣2 (B, 𝑡)
2 2
Or, la conservation de la masse donne 𝑣(A, 𝑡) 𝑆 = 𝑣(B, 𝑡) 𝑠 d’où 𝑣(A, 𝑡) ≪
𝑣(B, 𝑡) car 𝑠 ≪ 𝑆. Finalement

𝑣(B, 𝑡) = √2𝑔ℎ(𝑡) [formule de Torricelli] (2.13)

On remarquera que la vitesse a la même expression que celle de la

chute libre d’un point matériel dans le champ de pesanteur. Le débit
volumique d’écoulement vaut donc :

𝑄𝑉 = 𝑠 𝑣 = 𝑠√2𝑔ℎ(𝑡)

Pour connaître l’évolution de la hauteur d’eau, il faut relier 𝑣 à ℎ(𝑡) :

dℎ 𝑄 𝑠
𝑣(A, 𝑡) = − = 𝑉 = √2𝑔ℎ
d𝑡 𝑆 𝑆
L’intégration de cette équation donne un temps de vidange

𝑆 2ℎ0
𝜏= √ avec ℎ0 = ℎ(𝑡 = 0)
𝑠 𝑔

En pratique, le jet de sortie est contractée. La section effective de sortie

est donc plus petite que la section de l’orifice. Si l’on veut tenir compte
de ce phénomène il faut remplacer 𝑠 par 𝛼𝑠 où 𝛼 est le coefficient de
contraction.

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.

Explications – Le long d’un tube horizontal, d’après le théorème de

Bernoulli on a
1 te
2𝜌 𝑣 + 𝑝 = C
2

𝑣 𝑆 = Cte
 2.4 Écoulements permanents et incompressibles 21

Ainsi, la conservation du débit impose une augmentation de vitesse

au niveau de l’étranglement et la relation de Bernoulli impose alors
une dépression au même niveau.

Cet effet peut être mis à profit pour les applications suivantes :

▶ Mesure de débit ou de vitesse ;

▶ Principe des trompes à eau montées sur les robinet des paillasses
de chimie ; FIG. 2.9 : Effet Venturi – By Gringotu-
▶ Douchette venturi produisant une économie d’eau ; madre (Own work) [Public domain], via
Wikimedia Commons
▶ Amélioration du tirage d’une cheminée, principe du carburateur,
vaporisateur, etc.

Exercice – Montrer que dans un tube de venturi, la dépression au niveau

du rétrécissement est lié au débit par la relation Δ𝑝 = 𝐾 𝑄2𝑉 où 𝐾 est un
facteur géométrique.
Rép. En appliquant la relation de Bernoulli entre un point situé au col du
tube et un point loin du tube on obtient Δ𝑝 = 1/2𝜌 (𝑣col 2
− 𝑣∞2
). Dans le
cas où l’écoulement est incompressible, le débit volumique se conserve :
𝑄𝑉 = 𝑣col 𝒜col = 𝑣∞ 𝒜∞ où 𝐴col désigne l’aire de la section au niveau du
col. Si l’on suppose un fort rétrécissement, on peut négliger 𝑣∞
2
devant 𝑣col
2

d’où
2 𝜌
Δ𝑝 = 𝐾 𝑄𝑉 avec 𝐾 ≃
2 𝐴2col

Tube de Pitot

Le tube de Pitot permet la mesure de la vitesse d’écoulement d’un gaz

subsonique (𝑣 ≪ 𝑐𝑠𝑜𝑛 ). On peut en effet le considérer incompressible
dans ce cas. On pratique dans un tube un orifice de prise de pression
en A et en B. Le point A est un point d’arrêt car la vitesse est nulle
(il n’ y a pas d’écoulement dans l’orifice , c’est juste une prise de
pression). Loin du tube de Pitot l’écoulement est supposé uniforme
de vitesse 𝑣∞ et de pression 𝑝0 . En B la pression vaut 𝑝0 car les lois
de l’hydrostatique s’appliquent dans une direction perpendiculaire à
un écoulement parallèle permanent incompressible.

B
•
A Tube de Pitot
• Manomètre

FIG. 2.10 : Tube de Pitot.

En A (point d’arrêt), en utilisant la relation de Bernoulli, la pression

vaut
1 2
𝑝𝐴 = 𝑝0 + 𝜌𝑣∞
2
En B, 𝑝𝐵 = 𝑝0 . d’où
2Δ𝑝
𝑣∞ = √
𝜌
avec Δ𝑝 = 𝑝𝐴 − 𝑝𝐵 différence de pression mesurable à l’aide d’un
manomètre.
22 2 DYNAMIQUE DES FLUIDES PARFAITS

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.
 FLUIDES VISQUEUX 3
Ce cours aborde l’étude des fluides réels, c’est-à-dire présentant de 3.1 Notion de viscosité . . . . 23
la viscosité. Nous voyons comment l’équation d’Euler et la relation Fluides newtoniens . . . . 23
de Bernoulli se transforment pour tenir compte des effets dissipatifs Mesure de viscosité . . . 24
dans le fluide. On introduit également les notions de perte de charge Fluides non newtoniens . 24
et de coefficients aérodynamiques. 3.2 Dynamique d’un écoule-
ment visqueux . . . . . . 26
Version en ligne Bilan des forces . . . . . . 26
Équation de Navier-Stokes 27
[Link] Le nombre de Reynolds . 27
[Link] 3.3 Pertes de charge . . . . . 29
Loi de Poiseuille . . . . . 29
Analogie électrique . . . 30
Notion de perte de charge 31
3.1 Notion de viscosité Théorème de Bernoulli
généralisé . . . . . . . . . 33
Nous avons vu au Chapitre 2 que dans un fluide parfait, la contrainte 3.4 Traînée et portance . . . 35
qui s’exerce sur une particule de fluide est toujours perpendiculaire Formule de Stokes . . . . 36
aux parois de celle-ci. Dans un fluide réel en écoulement, la contrainte Analyse dimensionnelle . 36
possède une composante tangentielle dite contrainte visqueuse. Coefficients aérodyna-
miques . . . . . . . . . . 37

Fluides newtoniens

Expérience de Couette – Considérons un fluide enfermé entre deux

cylindres, l’un mobile, l’autre fixé via un fil de torsion. On constate
que lorsque la cavité cylindrique extérieure est mise en rotation à
la vitesse angulaire 𝜔, le cylindre intérieur tourne d’un angle 𝛼 par
rapport à sa position d’équilibre.
Analysons en détail le phénomène. Fil de torsion

1. La torsion du fil conduit à l’existence d’un couple dont les forces α

de pression ne peuvent pas être responsables. On est donc obli-
gé d’admettre l’existence d’efforts tangentiels. Fluide visqueux
2. On observe que les particules de fluide adhèrent aux parois. Il
existe donc un gradient de vitesse au sein de l’écoulement.
3. Pour les fluides simples, l’angle 𝛼 augmente proportionnelle-
ment a 𝜔. Les efforts tangentiels augmentent donc proportion-
nellement au gradient de vitesse.
Interprétation – L’expérience montre que, lors de l’écoulement d’un
fluide, la pression ne suffit pas à expliquer les phénomènes et qu’il ω
convient d’introduire des forces tangentielles qui s’opposent au mou- FIG. 3.1 : Expérience de Couette.
vement du fluide. Ces forces, de type frottement, dues aux interac-
tions entre molécules du fluide, sont appelées forces de viscosité. La
contrainte (force par unité de surface) ⃗⃗⃗⃗
𝜎 qu’exerce une couche de
fluide supérieure sur un élément de surface d’une couche de fluide
inférieure, s’écrit :

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹
⃗⃗⃗⃗
𝜎1→2 = 𝑛 + 𝜎𝑡 𝑡 ⃗ avec 𝜎𝑛 = −𝑝
= 𝜎𝑛⃗⃗⃗⃗
d𝑆
 24 3 FLUIDES VISQUEUX

⃗⃗⃗⃗⃗
𝑛 𝑡⃗ Fluide newtonien – Entre deux couches successives de fluide en écou-
⃗⃗⃗⃗⃗⃗⃗⃗
𝜎𝑡
lement unidimensionnel à la vitesse 𝑣,⃗ il existe des contraintes tan-
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝜎𝑛
gentielles à l’écoulement qui accélèrent la couche la plus lente et
ralentissent la couche la plus rapide. Par définition d’un fluide new-
tonien, les forces visqueuses sont proportionnelles à la différence de
FIG. 3.2 : Contraintes internes dans le vitesse c’est-à-dire au gradient de vitesse.
fluide.

𝜕𝑣
𝜎𝑡 = 𝜂 = 𝜂 𝛾̇ ♡ (3.1)
𝜕𝑛

où 𝜕𝑣/𝜕𝑛 désigne le gradient de vitesse dans la direction normale à

la surface. De manière générale, la contrainte visqueuse varie comme
la vitesse de cisaillement 𝛾.̇ La constante de proportionnalité 𝜂 est ca-
ractéristique du fluide et désigne la viscosité dynamique du fluide.

Mesure de viscosité

L’analyse dimensionnelle de la relation (3.1) donne

[𝐹 ] L [𝐹 ]
[𝜂] = = 2 T
L2 [𝑣] L

Ainsi, la viscosité est homogène à une pression × temps. On l’ex-

1 : Jean-Louis Marie Poiseuille (1797– prime indifféremment en [Link] (Pa.s) ou en poiseuille (Pℓ)
1869) fut élève de l’École Polytech- en hommage à Jean-Louis Marie Poiseuille1 .
nique avant d’étudier la médecine.
Les recherches de Poiseuille ont porté Le viscosimètre est l’appareil de mesure de la viscosité. Différents
principalement sur l’hémodynamique,
types de viscosimètre existent suivant le type de fluide utilisé. Pour
c’est-à-dire l’étude de la circulation
sanguine et lui ont permis de dégager les liquides, on utilise essentiellement le viscosimètre de Couette ou
une loi sur l’écoulement des fluides le viscosimètre à tube capillaire.
visqueux dans des tubes capillaires.
Ordres de grandeur – Pour les liquides, la viscosité varie fortement
avec la température (elle diminue lorsque la température augmente).
Pour des liquides purs, elle suit une loi du type

𝜂 ∝ 𝑒𝑏/𝑇

Quant aux gaz, leur viscosité est plus difficile à mesurer car beaucoup
Liquide (20 °C) Viscosité (Pa.s)
plus faible que celle des liquides. Sa détermination peut se faire à
l’aide d’une :
Eau 1, 0.10−3
Huile d’olive 0,84 ▶ mesure de vitesse (viscosimètre à bille roulante, viscosimètre à
Glycérine pure 1,5
Mercure 1, 5.10−3
tube capillaire) ;
▶ mesure de fréquence de résonance d’une onde de cisaillement
Gaz Viscosité (Pa.s)
(viscosimètre à cristal piézo-électrique de torsion).
Vapeur d’eau (20 °C) 9, 7.10−6
Air sec (20 °C) 18, 2.10−6 Elle dépend peu de la pression
√ et augmente légèrement avec la tem-
He (25 °C) 19, 9.10−6 pérature (à peu près comme 𝑇 ).
N2 (25 °C) 17, 7.10−6

TAB. 3.1 : Quelques valeurs de viscosité.

Fluides non newtoniens

Le comportement newtonien (𝜎 = 𝜂𝛾)̇ s’observe :

▶ dans tous les gaz ;
 3.1 Notion de viscosité 25

▶ dans les liquides simples constitués de petites molécules (l’eau

par exemple) ;
▶ dans les solutions contenant des ions ou molécules à symétrie
sphérique.

Cependant la rhéologie2 montre qu’il existe des fluides pour lesquels 2 : Étude du comportement des
fluides en écoulement
la relation entre contrainte tangentielle et cisaillement est plus com-
plexe. Certains fluides vérifient la relation

𝜎𝑡 = 𝜂(𝛾)̇ 𝛾̇

où 𝜂(𝛾)̇ représente une viscosité apparente. Lorsque 𝜂(𝛾)̇ diminue

avec 𝛾,̇ le fluide coule d’autant plus facilement qu’il est cisaillé. On
parle alors de fluide rhéofluidifiant (sang, polymère fondu, etc.). Le
comportement inverse est désigné par le terme rhéoépaississant (ami-
don+eau). Il existe également des liquides, comme les peintures, qui
ne coulent que si la contrainte dépasse un valeur seuil.

FIG. 3.3 : Mesures de la viscosité de

polysaccharides de différentes masses
moléculaires en solution aqueuse à
3% en masse.

Comportement visco-élastique – Tout fluide se caractérise par un

temps de relaxation viscoélastique 𝜏ve . Lorsque un fluide est soumis à
une contrainte, on distingue trois types de comportement en fonction
du temps d’observation 𝑡.

▶ si 𝑡 ≪ 𝜏ve , le fluide adopte un comportement élastique (défor-

mation proportionnelle à la contrainte) ;
▶ si 𝑡 ≫ 𝜏ve le fluide adopte un comportement visqueux (vitesse
de cisaillement proportionnelle à la contrainte 𝜎 = 𝜂𝛾)̇ ;
▶ si 𝑡 ≃ 𝜏ve , le comportement est alors plus complexe ; il est dit
visco-élastique.

C’est pourquoi, du point de vue mécanique, la distinction entre un

solide et un liquide est artificielle. Ce que l’on appelle communé-
ment un liquide est un fluide de petit temps de relaxation (𝜏𝑣𝑒 = 1 ns
pour l’eau) et ce que l’on appelle un solide peut être vu comme un
fluide de grand temps de relaxation (𝜏𝑣𝑒 = 106 ans pour le manteau
 26 3 FLUIDES VISQUEUX

de la croûte terrestre). 𝜏𝑣𝑒 dépend fortement de la température ce

qui confère à certains systèmes un comportement fluide ou solide
suivant la température (bitume par exemple).
Matériaux 𝜏𝑣𝑒 (s)
Un exemple de fluide viscoélastique est la pâte de silicone connue
eau (20 °C) 1 ns
verre à vitre (400 °C) 32 ans sous le nom de “silly-putty”. Une boule de “silly-putty” rebondit sur
verre à vitre (20 °C) 1026 − 1030 𝑠 le sol comme une balle élastique (aux temps courts) mais s’étale
bitume (−5 °C) 10 s comme un fluide visqueux (aux temps longs) si on la pose sur une
bitume (40 °C) 1 ms
surface horizontale.
TAB. 3.2 : Quelques valeurs de temps
de relaxation visco-élastique.

3.2 Dynamique d’un écoulement visqueux

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.

Bilan des forces

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 :

⃗⃗⃗ ⃗⃗ = (⃗⃗𝑓⃗ ⃗v⃗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-

tante des forces visqueuses dans le cas particulier simple d’un écou-
lement suivant (𝑂𝑥) avec un gradient de vitesse suivant (𝑂𝑦) :
𝑦
⃗⃗⃗⃗⃗
𝑛 𝑣 ⃗ = 𝑣(𝑦)⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥
𝑡⃗
𝜎𝑡 (𝑦 + d𝑦)
⃗⃗⃗⃗⃗⃗⃗⃗
On remarque ici que div𝑣 ⃗ = 0. L’écoulement est donc bien incompres-
𝑡⃗
sible. Dans ce cas, la résultante des forces visqueuses s’exerçant sur
⃗⃗⃗⃗⃗⃗⃗⃗
𝜎𝑡 (𝑦)
une particule de fluide, s’écrit :
⃗⃗⃗⃗⃗
𝑛
2
FIG. 3.4 : Bilan des forces de viscosité d𝑣 d𝑣 d 𝑣⃗
sur un élément de fluide. 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

d⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝜂 = 𝜂Δ𝑣 ⃗ d𝜏 ♡ (3.2)

où Δ est l’opérateur laplacien (cf. Chapitre A).

 3.2 Dynamique d’un écoulement visqueux 27

Équation de Navier-Stokes

D’après la seconde loi de Newton appliquée à une particule de fluide,

on a :
D𝑣 ⃗ ⃗ d𝜏 + ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝜌d𝜏 ∇𝑝 d𝜏 + 𝑓⃗⃗⃗v,ext
= −⃗⃗⃗⃗⃗⃗ d𝐹𝜂 avec ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝐹𝜂 = 𝜂Δ𝑣 ⃗ d𝜏
D𝑡
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

𝜕𝑣⃗
𝜌[ + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗ ∇𝑝 + 𝑓⃗⃗⃗v,ext
∇)𝑣]⃗ = −⃗⃗⃗⃗⃗⃗ ⃗ + 𝜂Δ𝑣 ⃗ (3.3)
𝜕𝑡

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 3 : Voir [Link]
7 problèmes du millénaire mis à prix $ 1 000 000 par l’Institut Clay3 ! org/millennium-problems/
navier-stokes-equation
Conditions aux limites – L’équation de Navier-Stokes étant une équa-
tion du second ordre, sa résolution introduit deux constantes d’inté-
gration pour la pression 𝑝 et pour la vitesse 𝑣.⃗ On les détermine en
appliquant les conditions aux limites suivantes :

▶ continuité de la vitesse à la traversée d’une interface ;

▶ continuité de la contrainte normale et donc de la pression ;
▶ 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 convectif – et
d’un terme du second ordre – le terme de 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’importance du terme
convectif devant le terme de viscosité. On peut estimer l’ordre de
grandeur du terme convectif et du terme visqueux à partir de l’échelle
caractéristique 𝐷 du problème, de la vitesse moyenne d’écoulement
𝑣, de la masse volumique 𝜌 du fluide et de sa viscosité 𝜂.
2
𝑣 𝑣
∥𝜌 (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣∥⃗ ≈ 𝜌 et ‖𝜂Δ𝑣‖⃗ ≈ 𝜂
𝐷 𝐷2
D’où le nombre sans dimension appelé nombre de Reynolds

terme convectif 𝜌𝑣𝐷

𝑅e = = ♡ (3.4)
terme visqueux 𝜂
 28 3 FLUIDES VISQUEUX

TAB. 3.3 : Quelques ordres de grandeur

Type d’écoulement Nombre de Reynolds
du nombre de Reynolds.
Écoulement atmosphérique 𝑅e ≈ 1011
Écoulement sanguin dans l’aorte 𝑅e ≈ 104
Écoulement sanguin dans les capillaires 𝑅e ≈ 10−3
Domaine de l’aéronautique 𝑅e ≈ 107
Domaine de la microfluidique 𝑅e ≈ 10−3 − 1

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
𝜕𝑣⃗
𝜌 ∇𝑝 + 𝑓⃗⃗⃗ ⃗v⃗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
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

𝜕𝑣⃗
𝜌[ + (𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗ ∇𝑝 + ⃗⃗𝑓⃗ ⃗v⃗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 > 105 , l’écoulement devient turbulent, c’est-à-dire
que la vitesse en un point M varie dans le temps de façon erra-
tique. Dans ce cas, le problème étant analytiquement insoluble,
on utilise souvent des lois phénoménologiques associées à une
analyse dimensionnelle.

Exemple

Quel est l’ordre de grandeur du nombre de Reynolds associé à l’écoule-

ment autour d’une balle de tennis allant à la vitesse 𝑣 = 100 km.h−1 dans
l’air ?

Le diamètre d’une balle de tennis est de l’ordre de 7 cm, la masse volu-

mique de l’air de l’ordre de 1 kg.m−3 et sa viscosité de l’ordre de 2.10−5 de
sorte que
𝜌𝑣𝐷 1 × 100/3, 6 × 0, 07
𝑅e = ≃ ≃ 105
𝜂 2.10−5
L’écoulement est donc turbulent.
 3.3 Pertes de charge 29

3.3 Pertes de charge

Loi de Poiseuille

On s’intéresse à l’écoulement d’un fluide visqueux dans un long tube

cylindrique de rayon 𝑅 et de longueur 𝐿 ≫ 𝑅. Le tube est horizon-
tal (orienté suivant Oz) et l’écoulement est assuré grâce à l’existence
d’une différence de pression Δ𝑝 entre l’entrée du tube et la sortie du
tube.

M
• 𝑣(M)
⃗⃗⃗⃗
𝑧

FIG. 3.5 : Écoulement de Poiseuille. Po-

sition du problème.

Hypothèses de travail –

𝜕𝑣⃗
1. L’écoulement est permanent donc = 0⃗ ;
𝜕𝑡
2. L’écoulement est incompressible, par conséquent div𝑣 ⃗ = 0 ;
3. Le nombre de Reynolds est suffisamment petit pour supposer
un régime d’écoulement laminaire. En pratique, on considère
que c’est le cas, quand 𝑅e < 2000 ;
4. L’écoulement est parallèle à Oz et invariant par rotation autour
de l’axe Oz, d’où 𝑣 ⃗ = 𝑣(𝑟, 𝑧)⃗⃗𝑢⃗⃗𝑧 ;
5. Enfin, on néglige la pesanteur car 𝜌𝑔𝑅 ≪ Δ𝑝.

Calcul du champ de vitesse – Commençons par écrire l’équation de

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

𝜕𝑣⃗ 𝜕
𝑎(M,
⃗ 𝑡) = 𝑢𝑧 = 0⃗
+ 𝑣𝑧 𝑣𝑧 (𝑟)⃗⃗⃗⃗⃗⃗⃗⃗
𝜕𝑡 𝜕𝑧
L’accélération est nulle. En effet, les lignes de champ 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é.
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 donc à l’équation de Stokes :

 30 3 FLUIDES VISQUEUX

⃗⃗⃗⃗⃗⃗
∇𝑝 = 𝜂Δ𝑣.⃗ 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 équation apparemment paradoxale se résout
si les deux termes sont constants.
d𝑝 Δ𝑝 1 d d𝑣
=𝐾=− =𝜂 (𝑟 )
d𝑧 𝐿 𝑟 d𝑟 d𝑟

4 : Attention, ici Δ désigne une varia- où Δ𝑝 = 𝑝1 − 𝑝2 est la différence de pression4 entre l’entrée et la
tion et non un laplacien. sortie. En intégrant deux fois on obtient

Δ𝑝 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 ) (3.5)
4𝜂𝐿
𝑥
Le profil des vitesses est parabolique.

Calcul du débit volumique – Le débit volumique est le flux du vecteur

FIG. 3.6 : Profil de vitesse. vitesse à travers une section de la canalisation :
𝑅
𝜋𝑅4 Δ𝑝
𝑄𝑉 = ∬ 𝑣 ⃗ d𝑆⃗⃗⃗ ⃗ = ∫ 𝑣(𝑟) 2𝜋𝑟 d𝑟 =
0 8𝜂 𝐿

Ainsi, la différence de pression est directement reliée au débit volu-

mique par la formule

8𝜂𝐿
Δ𝑝 = 𝑄 [Formule de Poiseuille] ♡ (3.6)
𝜋𝑅4 𝑉

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 ?
Rép. 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 :
 3.3 Pertes de charge 31

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 Δ𝑝 = 𝑅𝐻 𝑄𝑉

Notion de perte de charge

Définition – La perte de charge est la pression supplémentaire qu’il

faut imposer entre les extrémités d’une canalisation pour assurer un
écoulement stationnaire et compenser le frottement visqueux. Deux
termes entrent dans le calcul des pertes de charge :

▶ La perte de charge en ligne dite perte de charge régulière due

aux frottements le long du trajet.
▶ La perte de charge singulière due à la présence d’obstacles lo-
calisés tels que les coudes, les robinets, les vannes, les modifi-
cations brutales de section etc.

La perte de charge sera notée Δ𝑝𝜂 et s’exprime en Pa.

Essayons de donner une forme générale à l’expression de la perte de

charge dans une canalisation à l’aide d’une analyse dimensionnelle.
Pour cela, utilisons le théorème Π.

Théorème Π
Le théorème Π ou théorème de Vashy-Buckingham est le théo-
rème fondamental de l’analyse dimensionnelle. Supposons que
nous cherchions une relation entre 𝑛 grandeurs physiques 𝑔𝑖=1..𝑛
que l’on considère pertinentes pour décrire un phénomène. No-
tons 𝑘 le nombre de dimensions fondamentales utilisées par ces
grandeurs (𝑘 ≤ 7).
Il existe alors (𝑛 − 𝑘) produits sans dimension notées 𝜋𝑖 tels que
𝑓(𝜋1 , ..., 𝜋𝑛−𝑘 ) = 0

Considérons une conduite de longueur 𝐿, de diamètre 𝐷 traversée

par un fluide de viscosité 𝜂 et de masse volumique 𝜌 circulant à la
vitesse moyenne 𝑣. Supposons qu’il existe une relation entre Δ𝑝𝜂 , 𝐿,
𝐷, 𝜌, 𝜂 et 𝑣.

Grandeurs 𝐿 𝐷 𝜌 𝜂 Δ𝑝𝜂 𝑣

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

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

théorème Π, il existe trois nombres sans dimensions 𝜋1 , 𝜋2 et 𝜋3 tels
que 𝑓(𝜋1 , 𝜋2 , 𝜋3 ) = 0. Fabriquons donc trois nombres indépendants
sans dimension :

▶ 𝐿 et 𝐷 étant de même dimension, leur rapport est adimension-

𝐿
né : 𝜋1 = ;
𝐷
 32 3 FLUIDES VISQUEUX

▶ Le théorème de Bernoulli nous enseigne que 12 𝜌𝑣 2 est homo-

Δ𝑝𝜂
gène à une pression. Ainsi 𝜋2 = ;
1/2𝜌𝑣 2
▶ Enfin, on sait que le nombre de Reynolds est sans dimension ;
𝜌𝑣𝐷
ce sera 𝜋3 = 𝑅e =
𝜂

Ces trois nombres sont liés par une loi.

Δ𝑝𝜂
𝜋2 = 1 2
= 𝑓(𝜋1 , 𝑅e )
2 𝜌𝑣

Par ailleurs, l’expérience montre que la perte de charge régulière Δ𝑝𝜂

est proportionnelle à 𝐿. Autrement dit, 𝑓(𝜋1 , 𝑅e ) = 𝜋1 𝜆(𝑅e ) d’où

1 𝐿
Δ𝑝𝜂 = 𝜆(𝑅e ) 𝜌𝑣 ̄ 2 ♡ (3.7)
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 ne
dépend que du nombre de Reynolds pour une canalisation lisse. Dans
le cas d’une canalisation rugueuse, un quatrième nombre sans di-
mension 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. La valeur de 𝜆 peut être obtenue à l’aide
d’abaque comme le diagramme de Moody(FIG. 3.7).

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

Que vaut 𝜆 en régime laminaire ? – On remarque sur le diagramme de

Moody que pour 𝑅e < 2000, le coefficient de perte de charge suit une
 3.3 Pertes de charge 33

loi de puissance (ce qui donne une droite avec une échelle logarith-
mique). En effet, en régime laminaire, la perte de charge est donnée
par la formule de Poiseuille

8𝜂𝐿 8𝜂𝐿 𝐿
Δ𝑝𝜂 = 𝑄𝑉 ̄ 2
= 𝑣𝜋𝑅 = 32𝜂𝑣 ̄ 2
𝜋𝑅4 𝜋𝑅4 𝐷
Or le nombre de Reynolds de cet écoulement laminaire s’écrit :

𝜌𝑣𝐷
̄ 𝜌𝑣𝐷
̄
𝑅e = ⇒ 𝜂=
𝜂 𝑅e

Finalement, en régime laminaire

1 𝐿 64
Δ𝑝𝜂 = 𝜆 𝜌𝑣 ̄ 2 avec 𝜆 = ( ) ♡ (3.8)
2 𝐷 𝑅e

Les pertes de charge en régime turbulent – 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. En conclusion, pour
diminuer l’ensemble des pertes de charge dans une canalisation, afin
de diminuer les coûts de fonctionnement dus aux pompes, il faut,
quand c’est possible :

▶ diminuer la longueur de canalisation ;

▶ diminuer le nombre d’accidents sur la canalisation ;
▶ diminuer le débit de circulation ;
▶ augmenter le diamètre des canalisations ;
▶ faire circuler des liquides le moins visqueux possible ;
▶ utiliser des matériaux de faible rugosité.

Dans l’expression du nombre de Reynolds pour une conduite non circu-

laire, il est d’usage d’utiliser le diamètre hydraulique 𝐷𝐻 = périmètre
4×aire
.

Pertes de charges singulières – De la même manière, on peut expri-

mer les pertes de charge singulières comme suit :

1 2
Δ𝑃𝑠 = 𝜉 𝜌𝑣inc
2
où 𝜉 est le coefficient de perte de charge singulière et 𝑣inc est la
vitesse moyenne incidente du fluide arrivant sur l’obstacle. Il existe
également des tables donnant 𝜉 en fonction du type d’obstacle.

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
 34 3 FLUIDES VISQUEUX

é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’ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑣𝐵

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

FIG. 3.8 : Ecoulement stationnaire dans

un système de conduites. d𝑚 = 𝜌𝑣𝐴 d𝑡 𝑆𝐴

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𝑡

Le débit massique est uniforme le long de la canalisation puisqu’en

régime permanent

div𝜌 𝑣 ⃗ = 0 ⇒ ∬ 𝜌𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 = ∬ 𝜌𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆
𝒮1 𝒮2

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

formes sur la section droite de la conduite (cf. remarque page 35). 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𝑡
 3.4 Traînée et portance 35

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 de la perte de charge,
on a
𝛿𝑊𝜂 = −𝑄𝑉 Δ𝑝𝜂 d𝑡 ♡ (3.9)

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
(3.10)
Notons au passage que l’on retrouve la relation de Bernoulli vu dans
le chapitre sur les fluides parfaits à condition de faire 𝒫 = 0 et Δ𝑝𝜂 =
0.

La relation (3.10) 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ù 𝑒𝑐1 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.

3.4 Traînée et portance

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
 36 3 FLUIDES VISQUEUX

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

ligne d’écoulement

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑣∞ ⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝑡

FIG. 3.9 : Traînée sur un obstacle sphé-

rique immobile.

nous contenter de la solution sans chercher à la justifier. Stokes

obtient qu’une sphère de rayon 𝑟, immobile, soumise à un écoule-
ment permanent incompressible et visqueux, ressent une force de
traînée ⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝑡 proportionnelle à la vitesse d’écoulement et à la taille de
la sphère.
⃗⃗⃗⃗⃗⃗⃗⃗ 𝑣∞ pour 𝑅e ≪ 1 ♡
𝐹𝑡 = 6𝜋𝜂𝑟 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ (3.11)

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 bille et d’autre part
Vitesse de sédimentation
aux forces visqueuses. Si l’on étudie la chute d’une bille sphérique
dans un fluide visqueux 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
autre liquide non miscible) utilise où 𝑣 ⃗ représente la vitesse de la bille dans le référentiel du labo-
la loi de Stokes car le nombre de ratoire. Cette loi est vérifiée avec une précision meilleure que 1%
Reynolds est assez petit. Lorsqu’on pour 𝑅e < 0, 3. Cette contrainte reste cependant assez forte. En ef-
laisse reposer un liquide conte-
nant de petites particules solides
fet, pour une bille de 1 cm de diamètre tombant dans l’air cela im-
(comme par exemple un mélange pose 𝑣 < 0, 5 mm.s−1 ce qui signifie que la loi du frottement linéaire
eau-argile), les particules vont dé- n’est pas valable (sauf au tout début) dans ce cas. Par contre si la
canter c’est-à-dire sédimenter au chute s’effectue dans un liquide visqueux tel le glycérol (grosso mo-
fond du récipient avec une vitesse
qui dépend de leur dimension ca-
do mille fois plus dense que l’air et un million de fois plus visqueux),
ractéristique. Le temps de décan- la contrainte devient 𝑣 < 0, 5 m.s−1 . Dans ce cas, la loi de Stokes peut
tation donne alors un renseigne- être utilisée si la bille n’est pas trop pesante.
ment sur la taille des grains. En
effet les grains tombent à une vi-
tesse constante pour laquelle le
poids apparent (poids moins la Analyse dimensionnelle
poussée d’Archimède) compense
la force de traînée :
L’analyse précédente n’est valable qu’à petit nombre de Reynolds et
4 3
6𝜋𝜂𝑟𝑣sed = 𝜋𝑟 (𝜌𝑠 − 𝜌𝑙 )𝑔 pour un obstacle sphérique. Plaçons un obstacle quelconque dans
3
2
un écoulement stationnaire de vitesse ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑣∞ loin de l’obstacle et cher-
d’où 𝑣sed = 9𝜂 (𝜌𝑠 − 𝜌𝑙 )𝑟2 .
chons la force 𝐹 due à l’écoulement. Le traitement analytique est pos-
sible pour des géométries simples et pour des valeurs faibles de 𝑅e
 3.4 Traînée et portance 37

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. Montrons par une analyse dimensionnelle
quelle forme doit prendre cette force 𝐹 .

Quelles sont les grandeurs physiques pertinentes du problème ? Il

faut tout d’abord préciser que la viscosité est essentielle pour justi-
fier l’existence d’une force de frottement. Si le fluide est parfait il n’y
a pas de traînée car les couches de fluide glissent sur l’obstacle. Les
paramètres pertinents sont donc 𝜂, 𝑣∞ la vitesse du fluide, 𝐷 une di-
mension caractéristique de l’obstacle, 𝜌 la masse volumique du fluide
et la force 𝐹 .

Grandeurs 𝑣∞ 𝐷 𝜌 𝜂 𝐹

Dimensions L.T−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 deux nombres sans dimensions 𝜋1 et 𝜋2 liés
entre eux. Choisissons le nombre de Reynolds comme premier facteur
adimensionné. Part ailleurs, la quantité 21 𝜌𝑣∞
2
est homogène à une
pression comme 𝐹 /𝐷 . Ainsi
2

𝜌𝑣∞ 𝐷 𝐹 /𝐷2
𝜋1 = 𝑅 e = et 𝜋2 = 2
𝜂 1/2𝜌𝑣∞

Finalement le théorème Π nous dit que :

1 2 2
𝜋2 = 𝑓(𝜋1 ) ⇒ 𝐹 = 𝜌𝑣 𝐷 𝑓(𝑅e ) (3.12)
2 ∞

Coefficients aérodynamiques

La force qu’exerce un fluide en écoulement autour d’un obstacle peut

se décomposer en deux composantes.

▶ 𝑣∞ : c’est la traînée ⃗⃗⃗⃗⃗⃗⃗⃗

Une composante parallèle à ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝐹𝑡 .
▶ Une composante perpendiculaire : c’est la portance ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝑝 .

Ces deux forces s’expriment comme le prévoit la formule (3.12). On

définit alors deux coefficients de frottement, le 𝐶𝑥 et le 𝐶𝑧 .

⃗⃗⃗ ⃗⃗
𝐹
⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝑝
ligne d’écoulement

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑣∞
⃗⃗⃗⃗⃗⃗⃗⃗
𝐹𝑡
FIG. 3.10 : Portance et traînée.

Coefficient de traînée – La formule (3.12) peut se réécrire :

1 2
𝐹𝑡 = 𝜌𝑣 𝑆 𝐶𝑥 (𝑅e )
2 ∞
 38 3 FLUIDES VISQUEUX

où S représente une surface caractéristique, en générale, la surface

frontale projetée. On constate expérimentalement que le 𝐶𝑥 est qua-
si constant en régime turbulent (𝑅e grand) ce qui correspond aux
situations courantes de l’aéronautisme, le nautisme, le cyclisme etc.
La traînée peut se décomposer en trois termes :
▶ la traînée visqueuse qui est liée aux frottements du fluide sur
l’obstacle ;
▶ la traînée de pression qui est liée à l’existence d’une dépression
dans le sillage de l’obstacle quand la couche limite se décolle ;
▶ la traînée induite par la portance.
La force de traînée que subit un véhicule (en l’air ou au sol) étant
opposée à sa vitesse, elle dissipe de l’énergie. La puissance dissipée
s’écrit :
1
𝒫 = − 𝜌𝑆𝐶𝑥 𝑣3
2
où 𝑣 est la vitesse du véhicule. La puissance dissipée est une fonction
cubique de la vitesse, elle ne devient donc importante qu’à haute vi-
5 : grosso modo, en dessous de tesse5 . Pour minimiser la consommation à grande vitesse, le concep-
60 km.h−1 pour une voiture, les frotte- teur aura intérêt à agir sur le produit 𝑆𝐶𝑥 . Le tableau ci-dessous
ments de roulement l’emportent sur le
frottement aérodynamique.
donne quelques exemples.

TAB. 3.4 : Aérodynamique de quelques véhicules.

Automobile 𝐶𝑥 𝑆𝐶𝑥 [m2 ] 𝐹𝑡 à 120 km/h [N] 𝒫 [kW] 𝒫 [Ch]

DAIHATSU UFE III (3 places) 0,168 0,235 160 5,33 7,2

CITROEN AX DIESEL (5 places) 0,31 0,570 389 13,0 17,6

« Automobile moyenne » 0,35 0,63 430 14,3 19,5

Coefficient de portance – Lorsque que l’obstacle solide présente trois

axes équivalents il ne peut pas exister de portance. C’est le cas de la
sphère et du cube. Dans le cas contraire, la portance fait intervenir
le coefficient de portance 𝐶𝑧 qui dépend de la forme du solide et de
l’écoulement :
1 2
𝐹𝑝 = 𝜌𝑣∞ 𝑆 𝐶𝑧 (𝑅e )
2
Par exemple, une aile d’avion présente un coefficient de portance qui
dépend :
▶ de l’angle d’attaque 𝛼. Lorsque cet angle augmente, la portance
augmente jusqu’à un angle 𝛼max pour lequel la portance est
maximum. Une fois cet angle dépassé, la portance s’effondre,
c’est le décrochage.
▶ du profil de l’aile, notamment de sa cambrure. Une aile symé-
trique d’angle d’attaque nul ( la corde est parallèle à ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑣∞ ) ne pré-
sente pas de portance. Notez que les hélices d’un hélicoptère
sont symétriques : leur portance est liée à leur inclinaison.
Pour un avion en vol, on cherche à avoir une faible traînée (pour
consommer moins de carburant) et un maximum de portance c’est-
𝐶𝑧
à-dire un rapport 𝐶 maximum. Ce rapport, appelé finesse de l’aile,
𝑥
est maximum pour un certain angle.
 TENSION SUPERFICIELLE 4
De nombreuses expériences de la vie quotidienne ne peuvent pas 4.1 Énergie de surface . . . . 39
trouver d’explication avec les lois vues jusqu’ici. C’est en tenant compte Tension de surface . . . . 39
des propriétés des interfaces que l’on peut les justifier notamment Conséquences . . . . . . 40
grâce au concept de tension superficielle. La capillarité est la science Rôle des tensio-actifs . . 41
qui s’intéresse à ces phénomènes et qui joue un rôle majeur dans de 4.2 Forces capillaires . . . . . 41
nombreux domaines scientifiques (climat, chimie de formulation, in- Mise en évidence . . . . . 41
Mesure de la tension super-
dustrie du verre etc.). On propose ici une présentation classique de la
ficielle . . . . . . . . . . . 42
capillarité ; pour les aspects dynamiques et une vision plus moderne
Théorème de Laplace . . 43
de cette science voir [1].
4.3 Mouillage . . . . . . . . . 45
Angle de contact . . . . . 45
Version en ligne Ascension capillaire . . . 46

[Link]
[Link]

4.1 Énergie de surface

Notion de tension de surface

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
simpliste (continuité de la vitesse et des contraintes pour un fluide
visqueux). En réalité, pour décrire correctement la physique aux in-
terfaces il faut tenir compte des interactions moléculaires à courte
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 • gaz ℱ2
superficielle.

Essayons de justifier cette notion par un modèle simple (voire sim-

pliste). Considérons un liquide ℱ1 en contact avec sa vapeur ℱ2 . Au
sein de chaque fluide les molécules subissent des interactions de • interface
très courte portée (interaction de Van Der Walls attractives). Ainsi,
dans ℱ1 , chaque molécule possède une énergie d’interaction 𝜖1 < 0
résultat de l’interaction attractive avec son voisinage immédiat. De
même on définira une énergie d’interaction 𝜖2 au sein du fluide ℱ2 .
• liquide ℱ1
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 mo-
léculaire est de l’ordre de la dimension 𝑎 d’une molécule. Appelons
alors 𝜖12 l’énergie d’interaction d’une ”molécule interfaciale”. Bien en- FIG. 4.1 : Interactions au sein d’un
fluide.
tendu, on a |𝜖12 | compris entre |𝜖1 | et |𝜖2 |. Si 𝑁 est le nombre de molé-
cules et 𝑁𝑠 le nombre de molécules à l’interface, l’énergie du liquide
ℱ1 vaut
𝐸1 = (𝑁 − 𝑁𝑠 )𝜖1 + 𝑁𝑠 𝜖12 = 𝑁 𝜖1 + 𝐸𝑠
 40 4 TENSION SUPERFICIELLE

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
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 𝛾 ≃ ♡ (4.1)
𝑎2

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

La tension superficielle est une grandeur positive qui caractérise une

interface. Ainsi la tension de surface de l’interface eau-air n’est pas
la même que celle de l’interface eau-huile. Elle s’exprime en J.m−2 .

Conséquences expérimentales

Forme des bulles et gouttes – Augmenter la surface d’un liquide

coûte donc de l’énergie. Ainsi un liquide adoptera une forme qui mi-
nimise 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 (on parle du phénomène de coalescence) et l’on
obtient après un certain temps de l’huile avec de l’eau au dessous.

FIG. 4.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.
Rép. Le volume des deux gouttes vaut 𝒱 = 8/3𝜋𝑟3 et l’interface eau-air
possède une aire totale 𝑆 = 8𝜋𝑟2 . Le volume √
se conservant, le rayon de la
goutte obtenu après coalescence vaut 𝑟′ = 3 2𝑟 et son interface mesure
√ 2
𝑆 ′ = 4 ( 3 2) 𝜋𝑟2 ≃ 6, 3𝜋𝑟2 < 𝑆. La transformation s’accompagne donc
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
une phase métastable en dessous de la pression de vapeur saturante.
On parle de retard à la vaporisation. Une simple perturbation locale
 4.2 Forces capillaires 41

peut suffire à déclencher la formation d’une première bulle de va-

peur. C’est ce phénomène qui fut employé dans les détecteurs de
particules du milieu du XXe siècle (chambres à bulles).
De manière analogue, il y a retard à la liquéfaction pour la vapeur sur-
saturante (phénomène utilisé dans un autre détecteur de particule :
la chambre à Wilson).

FIG. 4.3 : Production de particules dans

la première chambre à bulles à hydro-
Rôle des tensio-actifs gène liquide du CERN

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é-
diatement à la surface, avec la queue hydrophobe pointant à l’exté-
rieur de la surface. Ce processus s’accompagne d’une stabilisation
de la 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.

4.2 Forces capillaires

Mise en évidence

Expériences – Plongeons un contour métallique dans de l’eau savon-

neuse puis retirons 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 formant 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. 4.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 FIG. 4.4 : Mise en évidence des forces
capillaires.
tendent à minimiser l’aire du film d’eau savonneuse.
▶ Ces forces sont tangentes à l’interface. Membrane liquide
▶ 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 d𝑓⃗⃗⃗ ⃗
cadre rectangulaire dont un des coté est mobile, le liquide cherchant
à minimiser sa surface, il faut exercer une force sur la tige mobile
d𝑥
FIG. 4.5 : La tension superficielle peut
s’interpréter comme une densité li-
néique de forces.
 42 4 TENSION SUPERFICIELLE

pour maintenir la surface constante. On peut obtenir l’expression de

cette force à l’aide d’un raisonnement énergétique (méthode des tra-
vaux virtuels) : Supposons qu’un opérateur déplace de façon quasi-
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é à la
tige donne
d𝐸𝑐 = 0 = 𝛿𝑊 = 𝐹op d𝑥 + 𝛿𝑊𝑠
où 𝛿𝑊𝑠 est le travail résistant des forces capillaires. Or ces forces
dérivent d’une énergie potentielle :
Lorsque que l’on traite une membrane
d’eau savonneuse, il ne faut pas ou- 𝛿𝑊𝑠 = −d𝐸𝑆 = −2𝛾ℓ d𝑥 donc 𝐹op = 2𝛾ℓ
blier qu’il y a deux interfaces liquide-
gaz ce qui explique la présence du fac-
teur deux dans l’expression des forces D’après le principe des actions réciproques (troisième loi de Newton),
capillaires. cette force s’identifie à la résultante des forces capillaires s’exerçant
sur la tige.

À retenir
De façon général, on peut traiter une interface comme une mem-
brane tendue : chaque portion de surface est le siège de forces
capillaires réparties sur le contour 𝒞 délimitant la portion de sur-
face. Ces forces sont tangentes à l’interface, perpendiculaires en
tout point de 𝒞 et données par la relation

⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑓 = 𝛾dℓ ⃗⃗⃗⃗
𝑛

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

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

Mesure de la tension superficielle

dynamomètre

𝑔⃗⃗⃗⃗

anneau

boy

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

Méthode de l’anneau – 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 à étudier 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 force exercée sur l’an-

neau 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
 4.2 Forces capillaires 43

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𝜋𝑟𝛾 si 𝑟1 ≃ 𝑟2

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

liquide.

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

tension superficielle vaut quelques mN.m−1 . La tension superficielle
varie avec la température ; 𝛾 diminue quand la température augmente
jusqu’à s’annuler à la température critique.

TAB. 4.1 : Quelques valeurs de tension

Tension superficielle de quelques interfaces liquide-air
superficielle.
Liquides Température Tension superficielle [mN.m−1 ]
Mercure 18 °C 475
Eau 20 °C 73
Eau 80 °C 62
Huile d’olive 20 °C 32

Théorème de Laplace

Un petit contour circulaire pris dans la surface libre plane d’un li-
quide en équilibre 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 même contour ont
une résultante orientée vers l’intérieur de la sphère ; il faut donc une
surpression Δ𝑝 pour que l’équilibre existe. On voit immédiatement
que plus la courbure est importante et plus Δ𝑝 sera grand.

f2

R2 dx
2α
R1
2β
pint p ext
f1 dy

FIG. 4.7 : Existence d’une surpression

dans une bulle.

On peut facilement obtenir l’expression de la différence de pression à

l’aide d’un raisonnement énergétique. Considérons une bulle de gaz
contenue dans un liquide et appelons 𝛾 la tension superficielle de
l’interface gaz-liquide. Supposons que l’on fasse subir à la bulle une
transformation quasi statique en augmentant son rayon de d𝑅. L’aire
de l’interface augmente donc de d𝒮 = 8𝜋𝑅d𝑅 et le volume de la bulle
de d𝒱 = 4𝜋𝑅2 d𝑅. Appliquons le théorème de l’énergie cinétique en
choisissant comme système l’interface :

d𝐸c = 0 = 𝛿𝑊ext + 𝛿𝑊int

44 4 TENSION SUPERFICIELLE

Les forces extérieures sont les forces de pression 𝑝int et 𝑝ext de sorte
que
𝛿𝑊ext = −𝑝ext d𝒱 + 𝑝int d𝒱
Par ailleurs, les forces internes sont les forces de tension superfi-
cielles qui dérivent d’une énergie potentielle :

d𝑊int = −𝛾d𝑆 = −8𝛾𝜋𝑅d𝑅

On a donc
2𝛾
Δ𝑝4𝜋𝑅2 d𝑅 − 8𝛾𝜋𝑅 d𝑅 = 0 ⇒ Δ𝑝 = 𝑝int − 𝑝ext =
𝑅
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. Ainsi la surpression est d’autant plus impor-
tante que la courbure moyenne 𝐶 = 2/𝑅 est importante. La générali-
sation à une géométrie quelconque est donnée par la loi de Laplace-
Young :

Loi de Laplace-Young

La différence de pression entre deux milieux non miscibles sépa-

rés par une interface de tension superficielle 𝛾, est donnée par

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

où 𝑅1 et 𝑅2 sont les rayons de courbure principaux de la mem-

brane au point considéré. Ces rayons sont définis positifs quand
le rayon de courbure est du coté du milieu 1. Dans le cas d’une
interface sphérique, ces deux rayons de courbure s’identifient au
rayon de la sphère.

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 laisser
diffuser le gaz selon l’épaisseur de la membrane et la taille des mo-
lécules gazeuses. Le gaz emprisonné est en surpression par rapport
au liquide en vertu de la loi de Laplace-Young

2𝛾
Δ𝑝 =
𝑅
Ainsi la pression 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
 4.3 Mouillage 45

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.

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 → FIG. 4.8 : Lame de savon formant une
caténoide (©Berkeley Science Review).
∞). On peut aussi obtenir des lames avec deux rayons de courbures
opposées comme sur la FIG. 4.8 montrant une caténoïde.

4.3 Mouillage

Angle de contact
liq. ou gaz
Équilibre au contact de trois fluides – Déposons une petite quantité
de liquide 2 sur un autre liquide 1 plus dense. L’ensemble des points ⃗⃗𝛾⃗ ⃗23
liq. 2
en contact avec les deux liquides et l’air forme une ligne, dite ligne ⃗⃗𝛾⃗ ⃗13
𝜃
triple. Intéressons nous aux forces capillaires s’exerçant sur cette
ligne triple. L’équilibre n’est possible que si la résultante des forces ⃗⃗𝛾⃗ ⃗12
capillaires peut s’annuler ce qui définit l’angle de contact 𝜃 : liq. 1

⃗⃗𝛾⃗⃗12 + ⃗⃗𝛾⃗⃗23 + ⃗⃗𝛾⃗⃗31 = 0⃗

FIG. 4.9 : Définition de l’angle de
d𝑓⃗⃗⃗𝑖𝑗
⃗ contact.
où le vecteur ⃗⃗𝛾⃗⃗𝑖𝑗 =dℓ désigne la densité linéique de force capillaire.
Cet équilibre suppose que l’on puisse former un triangle (dit triangle
de Neumann) 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 contraire, le liquide 2 s’étale sur le li-
quide 1 : on dit qu’il y a mouillage 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.
Rép. Comme 73 > 18 + 32, la condition d’équilibre n’est pas respectée et
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 figure ci-contre, résultat d’un compromis entre le poids
qui tends à 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 𝜃 = 𝛾𝑆𝑉
 46 4 TENSION SUPERFICIELLE

Cette relation trouvée par Young en 1805 définit l’angle de contact

𝜃.
⃗⃗𝛾⃗ ⃗𝐿𝑉

On distingue trois cas de figure :

𝜃
⃗⃗𝛾⃗ ⃗𝑆𝑉 ⃗⃗𝛾⃗ ⃗𝑆𝐿
1. 𝜃 > 𝜋/2 : le liquide est non mouillant (exemple : verre-mercure-
absence de air)
mouillage 2. 𝜃 ∈ [0, 𝜋/2] : il y a mouillage partiel (exemple : verre-eau-air)
3. Lorsque l’angle de contact n’est pas défini, il y a mouillage com-
mouillage
partiel plet du liquide sur le substrat solide.

mouillage
complet Dans la direction perpendiculaire à la surface solide, les forces capillaires
ne se compensent pas ce qui signifie que le solide se déforme. On peut
FIG. 4.10 : Différents équilibres d’un li-
quide au contact d’un substrat solide. montrer que pour les solides courants, cette déformation est plus petite
que la taille d’un atome et donc négligeable.

Ascension capillaire

Quand on plonge un capillaire propre (tube étroit de rayon 𝑟) dans de

l’eau, on observe l’ascension d’une colonne d’eau dans le capillaire
malgré la pesanteur. 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. 4.11 : 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érique de rayon 𝑅.
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 −
𝑟
 4.3 Mouillage 47

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 𝜃 2𝛾 cos 𝜃 1
𝑝0 = 𝑝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).

Exemple

Pour l’eau dans un tube de verre propre l’angle de raccordement vaut 0

(l’eau mouille le tube). Si 𝑟 = 0, 01 mm on obtient
2𝛾 1
ℎ= = 1, 5 m
𝜌𝑔 𝑟

l’ascension peut donc être très importante. C’est ce qui explique par exemple
les remontées d’humidité par capillarité que l’on peut observer dans des
milieux poreux notamment dans certains bâtiments.

La loi de Jurin suppose que le ménisque est sphérique et donc que la

pression sous le ménisque est uniforme. Or rigoureusement, comme les
bords du ménisque sont plus haut que le centre, cette pression ne peut
pas être uniforme. Elle ne peut l’être qu’approximativement à condition
que l’élévation du ménisque soit négligeable devant la hauteur d’ascen-
sion ℎ. Il est facile de montrer que cette condition se traduit par

cos 𝜃 𝛾
𝑟≪ √ √
1 − sin 𝜃 𝜌𝑔

Pour les petits angles cela donne

𝛾
𝑟≪√ ≃ 2, 7 mm pour l’eau
𝜌𝑔
ANNEXES
 OPÉRATEURS DIFFÉRENTIELS A
A.1 L’opérateur gradient . . . 51
Complément de cours sur ce qu’il faut savoir à propos des opérateurs
Définition . . . . . . . . . 51
différentiels utilisés en physique. Propriétés . . . . . . . . . 52
A.2 L’opérateur divergence . 52
Définition . . . . . . . . . 52
Version en ligne Propriétés . . . . . . . . . 53
A.3 L’opérateur rotationnel . 54
[Link] Définition . . . . . . . . . 54
mecaflu_complement1.php Propriétés . . . . . . . . . 54
A.4 L’opérateur laplacien . . 55
Le laplacien scalaire . . . 55
Le laplacien vectoriel . . 56
A.5 Accélération d’une parti-
A.1 L’opérateur gradient cule de fluide . . . . . . . 56

Définition

L’opérateur gradient est un opérateur différentiel qui s’applique à un

champ scalaire (fonction scalaire dépendant de l’espace et du temps)
et le transforme en un champ vectoriel (vecteur dépendant de l’es-
pace et du temps). Il se lit « gradient » ou « nabla » et se note :

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑓(M, 𝑡) ou ⃗⃗⃗⃗⃗⃗
∇𝑓(M, 𝑡)

Dans le système de cordonnées cartésiennes le gradient s’exprime

ainsi :

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝜕𝑓(𝑥, 𝑦, 𝑧, 𝑡) 𝜕𝑓(𝑥, 𝑦, 𝑧, 𝑡) 𝜕𝑓(𝑥, 𝑦, 𝑧, 𝑡)

grad𝑓(𝑥, 𝑦, 𝑧, 𝑡) = ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥 + 𝑢𝑦 +
⃗⃗⃗⃗⃗⃗⃗⃗ ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧 ♡
𝜕𝑥 𝜕𝑦 𝜕𝑧
(A.1)
La TAB. A.1 donne les différentes expressions du gradient dans les
systèmes de coordonnées utilisés couramment en physique.

1 2
Exercice – Calculer le gradient des champs suivants : 𝑓(𝑥, 𝑦, 𝑧) = (𝑥 +
2
𝑦2 + 𝑧 2 ) et 𝑔(𝑟, 𝜃, 𝜑) = − 1𝑟 .

∇𝑓 = (𝑥, 𝑦, 𝑧) = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
Rép. ⃗⃗⃗⃗⃗⃗ OM et ⃗⃗⃗⃗⃗⃗ ∇𝑔 = 1
⃗⃗⃗ ⃗ ⃗⃗.
𝑢
𝑟2 𝑟

Système 𝑓(M, 𝑡) Expression de grad𝑓 TAB. A.1 : Expressions de l’opérateur

gradient dans différents systèmes de
coordonnées.
𝜕𝑓 𝜕𝑓 𝜕𝑓
Cartésien 𝑓(𝑥, 𝑦, 𝑧, 𝑡) 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ + 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗⃗ + ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢
𝜕𝑥 𝑥 𝜕𝑦 𝑦 𝜕𝑧 𝑧
𝜕𝑓 𝜕𝑓 𝜕𝑓
Cylindriques 𝑓(𝑟, 𝜃, 𝑧, 𝑡) 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗⃗ + ⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑢 + 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝜕𝑟 𝑟 𝑟𝜕𝜃 𝜃 𝜕𝑧 𝑧
𝜕𝑓 𝜕𝑓 𝜕𝑓
Sphériques 𝑓(𝑟, 𝜃, 𝜑, 𝑡) 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗⃗ + 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗⃗ + 𝑢
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝜕𝑟 𝑟 𝑟𝑑𝜃 𝜃 𝑟 sin 𝜃𝑑𝜑 𝜑
 52 A OPÉRATEURS DIFFÉRENTIELS

Propriétés

L’opérateur gradient est un opérateur linéaire et vérifie donc

⃗⃗⃗⃗⃗⃗
∇(𝛼𝑓 + 𝛽𝑔) = 𝛼⃗⃗⃗⃗⃗⃗
∇𝑓 + 𝛽 ⃗⃗⃗⃗⃗⃗
∇𝑔 avec (𝛼, 𝛽) ∈ ℝ2

Le gradient d’un produit de champs scalaires vaut

⃗⃗⃗⃗⃗⃗
∇(𝑓.𝑔) = 𝑓 ⃗⃗⃗⃗⃗⃗
∇𝑔 + 𝑔⃗⃗⃗⃗⃗⃗
∇𝑓

où 𝑓 et 𝑔 sont deux fonctions de l’espace et du temps.

Lien avec la différentielle – On peut définir le gradient à partir de

sa relation avec la différentielle. Soit M un point de l’espace et M’ un
point infiniment voisin, la différentielle d𝑓 représente la variation du
champ scalaire 𝑓 lorsque l’on se déplace de M à M’ à 𝑡 fixé :

def
∇𝑓(M, 𝑡) ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗
d𝑓 = 𝑓(M’, 𝑡) − 𝑓(M, 𝑡) = ⃗⃗⃗⃗⃗⃗ dℓ avec ⃗⃗⃗⃗⃗⃗⃗⃗
dℓ = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
MM’

En conséquence,

1 : La surface de niveau de 𝑓 est l’en-

semble des points M pour lesquels ▶ Le vecteur ⃗⃗⃗⃗⃗⃗
∇𝑓(M, 𝑡) est perpendiculaire à la surface de niveau1
𝑓(𝑀, 𝑡) conserve la même valeur à un de 𝑓 passant par M à l’instant 𝑡.
instant 𝑡 fixé. En dimension 𝑑 = 2, ▶ Le vecteur gradient est orienté vers les valeurs croissantes de 𝑓
cet ensemble donne une courbe de ni-
et sa norme mesure le taux de variation spatiale dans la direc-
veau.
tion de plus grande pente

d𝑓
∥⃗⃗⃗⃗⃗⃗
∇𝑓∥ =
dℓ

Exercice – Considérons le champ scalaire de l’espace bi-dimensionnel,

𝑓(𝑥, 𝑦) = 𝑥2 + 𝑦2 . Représenter les courbes de niveau puis calculer ⃗⃗⃗⃗⃗⃗
∇𝑓. Tra-
cer quelques vecteurs gradients.
Rép. Les courbes de niveau sont des cercles de centre O. On a ⃗⃗⃗⃗⃗⃗
∇𝑓 =
(2𝑥, 2𝑦) = 2⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
OM. Les vecteur gradients sont effectivement perpendiculaires
aux cercles.

A.2 L’opérateur divergence

Définition

L’opérateur divergence est un opérateur différentiel qui s’applique à

un champ vectoriel et qui renvoie un champ scalaire. Il se lit « diver-
gence » et se note :

⃗⃗⃗ ⃗⃗ 𝑡)
div𝐴(M, ou ⃗⃗⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗ 𝑡)
∇ ⋅ 𝐴(M,
 A.2 L’opérateur divergence 53

Cette notation permet de retenir l’expression de la divergence en co-

ordonnées cartésiennes :

𝜕/𝜕𝑥 𝐴𝑥
⃗⃗⃗ ⃗⃗ 𝑦, 𝑧, 𝑡) = ⎛
div𝐴(𝑥, ⎞⋅⎛ ⎞ = 𝜕𝐴𝑥 + 𝜕𝐴𝑦 + 𝜕𝐴𝑧
⎜ ⎟ ⎜
⎜ 𝜕/𝜕𝑦 ⎟ ⎜ 𝐴𝑦 ⎟
⎟ ♡
𝑑𝑥 𝑑𝑦 𝑑𝑧
⎝ 𝜕/𝜕𝑧 𝐴
⎠ ⎝ 𝑧 ⎠
(A.2)
La TAB. A.2 donne les différentes expressions de la divergence d’un
champ vectoriel exprimé dans différents systèmes de coordonnées.

⃗⃗⃗ ⃗ ⃗ = ∇ ⋅ 𝐴
⃗⃗⃗ ⃗ ⃗ TAB. A.2 : Expressions de la divergence
Système Expression de div𝐴
dans différents systèmes de coordon-
nées.
𝜕𝐴𝑥 𝜕𝐴𝑦 𝜕𝐴𝑧
cartésien + +
𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕(𝑟𝐴𝑟 ) 𝜕(𝐴𝜃 ) 𝜕𝐴𝑧
cylindriques + +
𝑟d𝑟 𝑟d𝜃 d𝑧
1 𝜕(𝑟2 𝐴𝑟 ) 1 𝜕(sin 𝜃 𝐴𝜃 ) 1 𝜕𝐴𝜑
sphériques + +
𝑟2 𝜕𝑟 𝑟 sin 𝜃 𝜕𝜃 𝑟 sin 𝜃 𝜕𝜑

⃗⃗⃗ ⃗⃗ 𝜃, 𝜑) = ⃗⃗𝑢⃗ ⃗𝑟⃗⃗ . Calculer la di-

Exercice – Considérons le champ vectoriel 𝐴(𝑟,
𝑟2
vergence de ce champ en tout point M autre que O.
⃗⃗⃗ ⃗⃗ = 0. On dit que 𝐴
Rép. On trouve div𝐴 ⃗⃗⃗ ⃗⃗ est un champ à flux conservatif (sauf
en O).

Propriétés

L’opérateur divergence est un opérateur linéaire et vérifie donc

⃗⃗⃗ ⃗⃗ + 𝛽 𝐵)
div(𝛼𝐴 ⃗⃗⃗ ⃗⃗ + 𝛽 div𝐵
⃗⃗⃗ ⃗⃗ = 𝛼 div𝐴 ⃗⃗⃗ ⃗⃗ avec (𝛼, 𝛽) ∈ ℝ2

La divergence d’un produit vaut

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

La divergence d’un champ est reliée au calcul du flux.

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.
ext
∬ ⃗⃗⃗ ⃗⃗
𝐴(M) ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 =∭ ⃗⃗⃗ ⃗⃗
div𝐴(M) d𝜏 ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗
avec div𝐴 ∇⋅𝐴 ⃗⃗⃗ ⃗⃗
M∈(𝑆) M∈𝑉

Sens physique – La divergence prend un sens bien précis en méca-

nique des fluides2 . Considérons une portion de fluide en mouvement 2 : Voir simulations à la page
[Link]/mecanique_des_-
dans un fluide décrit par le champ de vitesse 𝑣(M, ⃗ 𝑡). Au cours du fluides/[Link]
mouvement, le volume 𝒱 de cette portion varie suite aux déforma-
tions engendrées par l’écoulement. La divergence de la vitesse est
 54 A OPÉRATEURS DIFFÉRENTIELS

liée au taux de dilatation de la portion fluide par la relation

1 D𝒱
div𝑣 ⃗ =
𝒱 D𝑡

A.3 L’opérateur rotationnel

Définition

L’opérateur rotationnel est un opérateur différentiel qui transforme

un champ vectoriel en un autre champ vectoriel. Il se lit « rotation-
nel » et se note
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot 𝐴(M, ⃗⃗⃗ ⃗⃗ 𝑡) ou ⃗⃗⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗ 𝑡)
∇ ∧ 𝐴(M,
Cette notation permet de retenir l’expression du rotationnel en coor-
données cartésiennes :

𝜕 𝐴𝑥 𝜕𝐴𝑧 𝜕𝐴𝑦
⎛ −
⎛
⎜ 𝜕𝑥 ⎞
⎟
⎛
⎜ ⎞
⎟ ⎜
⎜ 𝜕𝑦 𝜕𝑧 ⎞
⎟
⎟
⎜ ⎟ ⎜
⎜ ⎟
⎟ ⎜ ⎟
⎜
⎜ ⎟
⎟ ⎜ 𝐴𝑦 ⎟ ⎜ 𝜕𝐴𝑥 𝜕𝐴𝑧 ⎟
𝜕
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot 𝐴 = ⎜ ⃗⃗⃗ ⃗⃗ ⎜ ⎟
⎟ ∧⎜
⎜ ⎟
⎟ =⎜
⎜
⎜ − ⎟
⎟ ♡ (A.3)
⎜
⎜ ⎟
𝜕𝑦 ⎟ ⎜⎜ ⎟
⎟ ⎜ 𝜕𝑧 𝑑𝑥 ⎟
⎟
⎜ ⎟ ⎜ ⎟
⎜ 𝜕 ⎟ ⎜ ⎜ 𝐴𝑧 ⎟
⎟ ⎜⎜ 𝜕𝐴𝑦 ⎟
𝜕𝐴𝑥 ⎟
⎝ 𝜕𝑧 ⎠ ⎝ −
⎠ ⎝ 𝑑𝑥 𝑑𝑦 ⎠

La TAB. A.3 donne les différentes expressions du rotationnel dans dif-

férents systèmes de coordonnées.

Propriétés

L’opérateur rotationnel étant linéaire, on a

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot (𝛼𝐴 ⃗⃗⃗ ⃗⃗ + 𝛽 𝐵)
⃗⃗⃗ ⃗⃗ = 𝛼 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot𝐴 ⃗⃗⃗ ⃗⃗ + 𝛽 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot𝐵 ⃗⃗⃗ ⃗⃗ avec (𝛼, 𝛽) ∈ ℝ2

Le rotationnel d’un gradient est nul.

rot ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ ∇ ∧ (⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑓 = ⃗⃗⃗⃗⃗⃗ ∇𝑓) = 0⃗

La divergence d’un rotationnel est nulle.

div (⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot𝐴) ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗
∇ ⋅ (⃗⃗⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗ = 0
∇ ∧ 𝐴)

TAB. A.3 : Expressions du rotationnel dans différents systèmes de coordonnées

Système ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot 𝐴 ⃗⃗⃗ ⃗ ⃗ = ⃗⃗⃗⃗⃗⃗⃗
∇∧𝐴 ⃗⃗⃗ ⃗ ⃗
𝜕𝐴𝑧 𝜕𝐴𝑦 𝜕𝐴𝑥 𝜕𝐴𝑧 𝜕𝐴𝑦 𝜕𝐴𝑥
cartésien ( − , − , − )
𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦
1 𝜕𝐴𝑧 𝜕𝐴𝜃 𝜕𝐴𝑟 𝜕𝐴𝑧 1 𝜕(𝑟𝐴𝜃 ) 1 𝜕𝐴𝑟
cylindrique ( − , − , − )
𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃
1 𝜕(sin 𝜃𝐴𝜑 ) 1 𝜕𝐴𝜃 1 𝜕𝐴𝑟 1 𝜕(𝑟𝐴𝜑 ) 1 𝜕(𝑟𝐴𝜃 ) 1 𝜕𝐴𝑟
sphérique ( − , − , − )
𝑟 sin 𝜃 𝜕𝜃 𝑟 sin 𝜃 𝜕𝜑 𝑟 sin 𝜃 𝜕𝜑 𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝑑𝜃
 A.4 L’opérateur laplacien 55

Le rotationnel d’un produit vaut

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot 𝑓 𝐴 ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗
∇ ∧ (𝑓 𝐴) ∇𝑓 ∧ 𝐴 ⃗⃗⃗ ⃗⃗ + 𝑓 ⃗⃗⃗⃗⃗⃗
∇∧𝐴 ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑓 ∧ 𝐴 ⃗⃗⃗ ⃗⃗ + 𝑓.⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot 𝐴 ⃗⃗⃗ ⃗⃗

Relation avec la circulation :

Théorème de Stokes
La circulation d’un champ vectoriel le long d’un contour 𝒞 fermé
et orienté est égal au flux du rotationnel de ce champ à travers
une surface 𝒮 délimité par 𝒞.

∮ ⃗⃗⃗ ⃗⃗
𝐴(M) ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗
dℓ = ∬ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot𝐴(M) ⃗⃗⃗ ⃗⃗ ⋅ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆
M∈𝒞 M∈𝒮

avec ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
d𝑆 orienté à partir du sens de parcours de 𝒞 et de la règle
du tire-bouchon.

Sens physique – En mécanique des fluides, le rotationnel du champ

de vitesse d’un fluide en écoulement est lié à la vitesse de rotation
Ω des particules de fluide au cours de leur mouvement.

⃗⃗⃗⃗⃗ 1
Ω = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot𝑣 ⃗
2

A.4 L’opérateur laplacien

Le laplacien scalaire

L’opérateur laplacien scalaire est un opérateur différentiel d’ordre

deux qui transforme un champ scalaire en un autre champ scalaire.
Le laplacien scalaire s’obtient en prenant la divergence du gradient
et se note △𝑓(M, 𝑡).

𝜕2𝑓 𝜕2𝑓 𝜕2𝑓

△𝑓(M, 𝑡) = div(⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑓) = ∇2 𝑓 = 2
+ 2 + 2 ♡ (A.4)
𝜕𝑥 𝜕𝑦 𝜕𝑧

La TAB. A.4 donne les expressions du laplacien scalaire dans différents

systèmes de coordonnées.

TAB. A.4 : Expressions du laplacien

Système Expression de △𝑓
dans différents systèmes de coordon-
nées.
𝜕2𝑓 𝜕2𝑓 𝜕2𝑓
cartésien 2
+ 2
+
𝜕𝑥 𝜕𝑦 𝜕𝑧2
1 𝜕 𝜕𝑓 1 𝜕2𝑓 𝜕2𝑓
cylindriques (𝑟 )+ 2 +
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃2 𝜕𝑧2
1 𝜕 𝜕𝑓 1 𝜕 𝜕𝑓 1 𝜕2𝑓
sphériques (𝑟2 )+ 2 (sin 𝜃 )+
2
𝑟 𝜕𝑟 𝜕𝑟 𝑟 sin 𝜃 𝜕𝜃 𝜕𝜃 2
𝑟2 sin 𝜃 𝜕𝜑
2
56 A OPÉRATEURS DIFFÉRENTIELS

Le laplacien vectoriel

Le laplacien s’applique également à un champ vectoriel. Dans ce cas

il renvoie un autre champ vectoriel et se note

⃗⃗⃗ ⃗⃗
△𝐴

Par définition, le laplacien vectoriel s’obtient à l’aide de l’identité

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
rot𝐴⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗
∇ ∧ (⃗⃗⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗
∇ ∧ 𝐴) ∇ (⃗⃗⃗⃗⃗⃗
∇ ⋅ 𝐴) ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
⃗⃗⃗ ⃗⃗ − ∇2 𝐴 grad(div𝐴) ⃗⃗⃗ ⃗⃗ − △𝐴
⃗⃗⃗ ⃗⃗

En coordonnées cartésiennes, les vecteur unitaires étant fixes, le la-

placien vectoriel d’un champ 𝐴 ⃗⃗⃗ ⃗⃗ est tout simplement, un vecteur dont
les composantes sont les laplaciens scalaires des composantes de
𝐴⃗⃗⃗ ⃗⃗ :
⃗⃗⃗ ⃗⃗ 𝑡) = (△𝐴 ) ⃗⃗⃗⃗⃗⃗⃗⃗
△𝐴(M, 𝑥 𝑢𝑥 + (△𝐴𝑦 ) ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑦 + (△𝐴𝑧 ) ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑧

A.5 Accélération d’une particule de fluide

On a vu au Chapitre 1 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 Ox de ce terme en utilisant l’égalité
⃗⃗⃗ ⃗⃗ ∧ (𝐵
𝐴 ⃗⃗⃗ ⃗⃗ ∧ 𝐶) ⃗⃗⃗ ⃗⃗ 𝐶)
⃗⃗⃗ ⃗⃗ = (𝐴. ⃗⃗⃗ ⃗⃗ 𝐵 ⃗⃗⃗ ⃗⃗ 𝐵)
⃗⃗⃗ ⃗⃗ − (𝐴. ⃗⃗⃗ ⃗⃗ 𝐶 ⃗⃗⃗ ⃗⃗ = 𝑣,⃗ 𝐵
⃗⃗⃗ ⃗⃗ avec 𝐴 ⃗⃗⃗ ⃗⃗ = ⃗⃗⃗⃗⃗⃗ ⃗⃗⃗ ⃗⃗ = 𝑢
∇𝑣𝑥 et 𝐶 ⃗⃗⃗⃗𝑥 :

(𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 ) 𝑢 ⃗⃗⃗⃗𝑥 ) ⃗⃗⃗⃗⃗⃗
⃗⃗⃗⃗𝑥 = (𝑣 ⃗ ⋅ 𝑢 ∇𝑣𝑥 − 𝑣 ⃗ ∧ (⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 ∧ 𝑢⃗⃗⃗⃗𝑥 )
1 ⃗⃗⃗⃗⃗⃗ 2
= 𝑣𝑥 ⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 − 𝑣 ⃗ ∧ (⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 ∧ 𝑢⃗⃗⃗⃗𝑥 ) = ∇𝑣 − 𝑣 ⃗ ∧ (⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 ∧ ⃗⃗⃗⃗⃗⃗⃗⃗
𝑢𝑥 )
2 𝑥
Ainsi en procédant de la même façon pour les deux autres compo-
santes, on obtient
1
(𝑣 ⃗ ⋅ ⃗⃗⃗⃗⃗⃗
∇) 𝑣 ⃗ = ⃗⃗⃗⃗⃗⃗ ⃗ ⃗⃗⃗⃗⃗⃗
∇ (𝑣𝑥2 + 𝑣𝑦2 + 𝑣𝑧2 )− 𝑣∧( ⃗⃗⃗⃗𝑥 + ⃗⃗⃗⃗⃗⃗
∇𝑣𝑥 ∧ 𝑢 ∇𝑣𝑦 ∧ 𝑢⃗⃗⃗⃗𝑦 + ⃗⃗⃗⃗⃗⃗
∇𝑣𝑧 ∧ 𝑢⃗⃗⃗⃗𝑧 )
2

On reconnait 𝑣2 dans le gradient et l’on voit apparaître ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗

rot𝑣 ⃗ dans le
dernier terme. On aboutit alors à une nouvelle expression de l’accé-
lération
𝜕 𝑣 ⃗ ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑣2
𝑎(M,
⃗ 𝑡) = + grad + (⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ rot𝑣)⃗ ∧ 𝑣 ⃗ ♡ (A.5)
𝜕𝑡 2
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. 39).
[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 −2
𝐺 Constante gravitationnelle 6,67430 × 10−11 m3 kg s

𝜖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 )
⃗⃗⃗ ⃗⃗
𝐵 Champ magnétique (T)

𝐶𝑥 , 𝐶𝑦 , 𝐶𝑧 Coefficients aérodynamiques (sans dimension)

ℰ Energie (J)
⃗⃗⃗⃗ ⃗
𝐸 Champ électrique (V.m−1 )

𝑒 Énergie volumique (J.m−3 )

⃗⃗⃗ ⃗⃗
𝐹 force (N)

𝑓⃗⃗⃗⃗ force volumique (N.m−3 )

𝑔⃗ Champ de pesanteur ([Link]−1 )

𝐽𝑚 Densité de courant massique (kg.m−2 )

𝐾𝑛 Nombre de Knudsen (sans dimension)

𝑀 Masse molaire ([Link]−1 )

𝑚 Masse (kg)

𝑛 Densité de particules (m−3 )

𝒫 Puissance (W)

𝑝⃗ Quantité de mouvement (kg.m.s−1 )

𝑝 Pression (Pa)

𝑄𝑚 Débit massique (kg.s−1 )

𝑄𝑉 Débit volumique (m3 .s−1 )

𝑅e Nombre de Reynolds (sans dimension)

𝑆 Surface (sa mesure en m2 )

𝑇 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, pulsation (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 temporelle de 𝑓(𝑡)

⟨𝑓⟩ Moyenne d’ensemble de 𝑓

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𝐴 ∇∧𝐴 ⃗⃗⃗ ⃗⃗ Rotationne d’un champ vectoriel

Δ𝑓 = ∇2 𝑓 Laplacien scalaire

∑ = ∑ ∑ Somme sur les couples (𝑖, 𝑗) avec 𝑖 ≠ 𝑗

couples (𝑖,𝑗) 𝑖 𝑗<𝑖
2023

[Link]/mecanique des fluides