La mécanique des uides

Hassan ADOCH
Professeur agrégé au Lycée Ibn Timiya
Classes préparatoires aux grandes écoles, Marrakech
 TABLE DES MATIÈRES

1
1.1

O C H
1 Étude phénoménologique des uides
Outils mathématiques . . . . . . . . . . .
Surface et volume . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
1
1
1

D
1.2 Opérateurs vectoriels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
2 Modélisation d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.1 Dénition d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

A
2.2 Modèle de milieu continu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3 Caractéristiques d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3.1 Pression d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3.2 Notion de viscosité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3.3 Condition aux limites cinématiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2 Cinématique des uides 7
1 Description du mouvement d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1 Description Lagrangienne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2 Description Eulérienne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3 Dérivée particulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 Conservation de la masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1 Débit massique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2 Équation de conservation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3 Écoulements particuliers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1 Écoulement incompressible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.2 Écoulement irrotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.3 Écoulement stationnaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4 Exemple d'application : Tourbillon de Rankin . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3 Dynamique des uides 15
1 Généralités . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.1 Rappels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2 Les forces usuels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2 Équation d'Euler - Théorème de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1 Équation d'Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.2 Théorème de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

i
 TABLE DES MATIÈRES

2.3 Exemple d'application : eet Venturi . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

3 Équation de Navier-Stockes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.1 Démonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.2 Nombre de Reynolds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.3 Exemple d'application : écoulement de Poiseuille . . . . . . . . . . . . . . . . . . . . . 20
4 Pertes de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4.1 Présentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4.2 Pertes de charge singulières . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4.3 Pertes de charge régulières . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
5 Force de traînée et de portance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

H
5.1 Présentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
5.2 Force de traînée et de portance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

C
4 Bilans dynamiques et thermodynamiques 30
1 Présentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1.1 Cadre d'étude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

O
1.2 Bilan de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2 Bilans dynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.1 Bilan de quantité de mouvement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31

D
2.2 Bilan de moment cinétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.3 Bilan d'énergie cinétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3 Bilans thermodynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34

A
3.1 Bilan énergétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.2 Bilan entropique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35

H.ADOCH ii / 35 adoch.hasaan@gmail.com
CHAPTER 1

ÉTUDE PHÉNOMÉNOLOGIQUE DES FLUIDES

Introduction

O C H
L'étude des uides est une partie très intéressante en physique. elle permet d'expliquer plusieurs

D
phénomènes à partir de comment vole un avion jusqu'à prédire la météo. Le but de ce chapitre est de
préciser le cadre d'étude et de dénir quelques vocabulaires.

A
1 Outils mathématiques
1.1 Surface et volume
En mécanique de point, on a vue trois systèmes des coordonnées à savoir cartésien, cylindrique et sphérique.
On rappelle que le vecteur déplacement est donné par :
En sphérique
En cartésien En cylindrique
−−−→
−−−→ −−−→ dOM = dr→−e r + rdθ→
−
eθ
dOM = dx→
−
e x + dy →
−
e y + dz →
−
ez dOM = dr→
−
e r + rdϕ→
−
e ϕ + dz →
−
ez
+r sin(θ)dϕ→−e ϕ
On peut dénir trois surfaces On peut dénir trois surfaces
élémentaires : élémentaires : On peut dénir trois surfaces
élémentaires :
−
→ −
→
dS x = dydz →
−
ex dS r = rdϕdz →
−
er −
→
dS r = r2 sin(θ)dϕdθ→
−
er
−
→ −→
dS y = dxdz →
−
ey dS ϕ = dr dz →
−
eϕ −
→
−
→ −
→ dS θ = r sin(θ)drdϕ→
−eθ
dS z = dxdy →
−
ez dS z = dr rdϕ→
−ez −
→
dS ϕ = dr r dθ→
−
eϕ
Le volume élémentaire est : Le volume élémentaire est :
Le volume élémentaire est :
dτ = dx dy dz dτ = r dr dϕ dz
dτ = r2 sin(θ) drdϕdθ
−
→
Pour déterminer la direction d'orientation de vecteur dS d'une surface, on a deux cas :
−
→
• Cas d'une surface fermée : dans ce cas le vecteur dS est toujours perpendiculaire sur la surface et orienté
vers l'extérieur.
• Cas d'une surface ouverte : Dans ce cas on remarque que cette surface appuie sur un contour fermé.

1
 CHAPTER 1. ÉTUDE PHÉNOMÉNOLOGIQUE DES FLUIDES

−
→
L'orientation de ce contour (généralement choisie arbitrairement) donne l'orientation de vecteur dS par la
règle de tire bouchon.

1.2 Opérateurs vectoriels

Lorsqu'on manipule les grandeurs vectorielles ou scalaires, on remarque l'apparition de quelques opérateurs
mathématiques qu'on propose de dénir dans cette partie.
• Le gradient : est un opérateur linéaire qui s'applique à un champ de scalaire et décrit un champ de
vecteurs qui représente la variation de la valeur du champ scalaire dans l'espace. Son expression dépend de
la base choisie. En eet on a :

H
Base cartésienne :
−−→ ∂f →
− ∂f →
− ∂f →
−
grad(f ) = ex+ ey+ ez

C
∂x ∂y ∂z
Base cylindrique :
−−→ ∂f →
− 1 ∂P →
− ∂f →
−
grad(f ) = er+ eϕ+ ez

O
∂r r ∂ϕ ∂z
Base sphérique :
−−→ ∂f →
− 1 ∂f →
− 1 ∂f →
−
grad(f ) = er+ eθ+ eϕ

D
∂r r ∂θ rsin(θ) ∂ϕ

• La divergence : est un opérateur linéaire qui s'applique à un champ vectoriel. Elle exprime la tendance

A
de ce champ à uer localement hors d'un petit volume entourant le point M où est calculée la divergence. En
eet on a :
Base cartésienne :
→
− ∂Ax ∂Ay ∂Az
div( A ) = + +
∂x ∂y ∂z
Base cylindrique :
→
− 1 ∂rAr 1 ∂Aϕ ∂Az
div( A ) = + +
r ∂r r ∂ϕ ∂z
Base sphérique :
→
− 1 ∂r2 Ar 1 ∂sin(θ)Aθ 1 ∂Aϕ
div( A ) = 2 + +
r ∂r r sin(θ) ∂θ r sin(θ) ∂ϕ

• Le rotationnel : est un opérateur mathématique qui transforme un champ de vecteurs en un autre champ
de vecteurs. Il exprime la tendance de ce champ à tourner autour du point étudié. Ses expressions sont :
Base cartésienne :
→
−
ex →
−
ey →
−
ez
     
−→ →− ∂Az ∂Ay →
− ∂Ax ∂Az →
− ∂Ay ∂Ax →
− ∂ ∂ ∂
rot( A ) = − ex+ − ey+ − ez=
∂y ∂z ∂z ∂x ∂x ∂y ∂x ∂y ∂z
Ax Ay Az

Base cylindrique :
     
−→ →− 1 ∂Az ∂Aϕ →
− ∂Ar ∂Az →
− 1 ∂(rAϕ ) ∂rAr →
−
rot( A ) = − er+ − eϕ+ − ez
r ∂ϕ ∂z ∂z ∂r r ∂r ∂ϕ

H.ADOCH 2 / 35 adoch.hasaan@gmail.com
 CHAPTER 1. ÉTUDE PHÉNOMÉNOLOGIQUE DES FLUIDES

Base sphérique :
     
−→ →− 1 ∂ (sin(θ)Aϕ ) ∂Aθ →
− 1 ∂Ar 1 ∂(rAϕ ) →
− 1 ∂(rAθ ) ∂Ar →
−
rot( A ) = − er+ − eθ+ − eϕ
r sin(θ) ∂θ ∂ϕ r sin(θ) ∂ϕ r ∂r r ∂r ∂θ

• Le laplacien : c'est un opérateur mathématique de deuxième ordre qui exprime la variation d'un champ
scalaire ou vectoriel dans l'espace. Son expression est donnée par les relations suivante
Base cartésienne :
∂ 2f ∂ 2f ∂ 2f
∆(f ) = + +
∂x2 ∂y 2 ∂z 2

H
Base cylindrique :
1 ∂ 2f ∂ 2f
 
1 ∂ ∂f
∆(f ) = r + +

C
r ∂r ∂r r2 ∂ϕ2 ∂z 2
Base sphérique :

O
∂ 2f
   
1 ∂ 2 ∂f 1 ∂ ∂f 1
∆(f ) = 2 r + 2 sin(θ) + 2 2
r ∂r ∂r r sin(θ) ∂θ ∂θ r sin (θ) ∂ϕ2

D
Remarques :
→
−
∗ En coordonnées cartésiennes on peut utiliser le vecteur nable 5 pour écrire les autres opérateurs.

A
En eet :
−−→ →
− →
− →
− →
− −→ →− →
− → − →
−
grad(f ) = 5(f ) ; div( A ) = 5. A ; rot( A ) = 5 ∧ A et ∆f = 5 2 f

∗ Le laplacien peut être appliquer à un scalaire ou à un vecteur généralement on utilise la même

→
−
notation dans les deux cas mais on peut trouver aussi la notation ∆ pour le laplacien vectoriel.
∗ Les opérateurs sont reliés entre eux par plusieurs relations. On mentionne ici :
−→ −−→  → − −→ → − 
rot grad(f ) = 0 div rot( A ) = 0
−−→  −→ −→ → −  −−→  →
−  →
−
div grad(f ) = ∆f rot rot( A ) = grad div( A ) − ∆ A

−−→ →
− → − −−→ →
−  → − → − −−→ →− → − −→ → − →
− −→ → −
grad A . B = A .grad B + B .grad A + A ∧ rot( B ) + B ∧ rot( A )
−→ →
− → − →− −−→ →
− → − −−→ → − → − →
− →
− →
−
rot A ∧ B = − A .grad B + B .grad A + A div( B ) − B div( A )
→
− →− → − −→ →− →
− −→ → −
div A ∧ B = B .rot( A ) − A .rot( B )
∗ Ces opérateurs permet de convertir une intégrale simple en une autre double et une double en une
triple. On mentionne ici deux théorèmes à savoir :
‹ ˚
− −
→ → →
−
• Théorème de Green-Ostrogradsky : A .dS = div( A ).dτ
Σ V
˛ ¨
− →
→ − −→ →−
• Théorème de Stockes : A . dl = rot( A ).dS
L S

H.ADOCH 3 / 35 adoch.hasaan@gmail.com
 CHAPTER 1. ÉTUDE PHÉNOMÉNOLOGIQUE DES FLUIDES

2 Modélisation d'un uide

2.1 Dénition d'un uide
On dénit un uide comme étant un milieu matériel qui peut s'écouler. Les gaz et les liquides sont les uides
étudiés classiquement. On peut dire aussi qu'un uide est un milieu susceptible de se déformer s'il est soumis
à des contraintes tangentielles. Ces milieux sont constitués d'un nombre colossale de particules (de l'ordre
de NA ). Ce qui implique que le traitement de ces particules ne peut pas être de manière individuelle d'où la
nécessité d'introduire un autre modèle.

H
2.2 Modèle de milieu continu
A l'échelle macroscopique, un uide est caractérisé par les paramètres d'état utilisées en thermodynamique à

C
savoir la température, la pression, la quantité de matière...etc. Alors que à l'échelle microscopique, le uide
est un ensemble de particules qui peuvent être décrit par leurs positions et leurs quantité de mouvement.
Le modèle de milieu continu consiste à moyenner les caractéristiques des particules constituants le milieu
sur un échelle dit l'échelle mésoscopique. Dans ce modèle, on remplace les caractéristiques microscopiques

O
des particules par des champs que l'on notes v(M, t), ρ(M, t), T (M, t) ...etc.

Dans cette échelle (l'échelle mésoscopique), on dénit une particule uide comme étant un volume dτ

D
√
du uide de masse xe dont la longueur typique d ∼ dτ est très grande devant le libre parcours moyen l
3

et très faible devant la longueur caractéristique L de l'écoulement.

A
l << d << L

Dans la majorité des cas, l est de l'ordre de µm et L est de l'ordre de cm jusqu'au mètre alors la longueur d
caractéristique des mm.

3 Caractéristiques d'un uide

3.1 Pression d'un uide
Dans un uide, chaque particule applique une force sur les particules qui sont en contact avec cette particule
−→
c'est-ce-qu'on appelle la force de contrainte. La composante normale dF n de cette force permet de dénir la
pression p du uide. La force de pression pression est alors donnée par :

−→ −
→ −
→
dF n = −p dS dS
Le signe moins exprime que le sens de la force est l'opposé du sens
−
→
de dS . La pression s'exprime en pascal (P a ≡ N/m2 ), on trouve dS
aussi :
1 bar = 105 P a ;1 atm = 101 325 P a
1
1 mmHg = 133, 322 P a = atm −→
760 dF n

En mécanique des uides, on utilise des grandeurs volumiques ce qui revient à deux causes principales :
la première c'est que le volume du uide étudié est généralement dicile à préciser. Alors que la deuxième

H.ADOCH 4 / 35 adoch.hasaan@gmail.com
 CHAPTER 1. ÉTUDE PHÉNOMÉNOLOGIQUE DES FLUIDES

c'est que les grandeurs utilisés pour étudier un uide dépend généralement de la position et ils ne sont pas
uniforme que dans une particule uide. La forme volumique de la force de pression est alors :

→
− −−→
d F n = −grad(p) dτ

3.2 Notion de viscosité

La partie tangentielle de la force de contrainte donne une autre force appelée force de viscosité ou force
de cisaillement. Cette force apparaît entre deux particules uide (ou entre deux couches du uide) lors du

H
mouvement du uide.
Les expériences montrons que lorsque un uide est en écoulement

C
dans une conduite avec une vitesse par exemple → −
v = vx (y)→
−
e x. Y
On remarque l'apparition d'un gradient de vitesse. Pour expliquer →
−
v (y + dy)
ces observations, on considère deux particules uides en contact

O
sur une surface dS . La force exercée par la première particule −→
uide sur la deuxième est donnée par l'expression : dF n
dS

D
−→ ∂vx →
−
dF v = η dS ex →
−
∂y v (y)

A
Le coecient η est appelé viscosité dynamique du uide.
Remarques :
• Le coecient de viscosité s'exprime en poiseuille (P l) qui est équivalent à P a.s. Ce coecient dépend
de la température et de la contrainte appliquée. On donne un ordre de grandeur de coecient de
viscosité dynamique, pour quelques uides :
Fluide Air Eau Glycérine Miel Graisse
η à 25 C (en P l) 10−5 10−3
o
1,5 10 103

• Ce coecient de viscosité quantie les frottements entre les couches du uide lors de l'écoulement.
En eet, la couche rapide applique une force sur la couche la moins rapide pour augmenter sa vitesse
(on peut dire l'inverse aussi).
• Dans certains contextes, on remarque l'apparition d'une "viscosité massique" c'est ce qu'on appelle
la viscosité cinématique et qui est dénie par :
η
ν=
ρ

• Dans le cadre de notre programme, on considère que le coecient de viscosité ne dépend pas de la
contrainte appliquée sur le uide, on appelle ce type les uides newtoniens.
• La formulation "générale" (volumique) de la force de viscosité sur une particule uide nécessite
l'addition des forces appliquées sur chacune des surfaces. Ce qui donne à trois dimensions pour un
uide newtonien en écoulement incompressible :

H.ADOCH 5 / 35 adoch.hasaan@gmail.com
 CHAPTER 1. ÉTUDE PHÉNOMÉNOLOGIQUE DES FLUIDES

→
−
d F v = η ∆→
−
v dτ

Telle que un écoulement est incompressible si le champ de vitesse vérie :

div(→
−
v)=0

H
3.3 Condition aux limites cinématiques
Les écoulements étudiés généralement ont lieu en présence d'un frontière qui peuvent être constituées par des

C
obstacles rigides (mobiles ou non). Il est donc important de préciser quelles sont les conditions aux limites
qui résultent pour le champ de vitesse. Ces conditions aux limites c'est ce qu'on appelle les conditions aux
limites cinématiques. On a deux cas :

O
• Cas d'un parois xe : Lorsque un uide en écoulement dans une conduite, les parois de celle-ci impose
un débit massique nul à travers la paroi (ceci est valable aussi pour un obstacle plongé dans l'écoulement).
Ce qui s'écrire dans le référentiel lié au parois sous la forme :

D
→
− −
→
v f luide (M, t).dS = 0

A
−
→
dS est le vecteur surface élémentaire correspondant à la surface de la paroi étudié. Il est donc perpendiculaire
sur la surface considérée.
• Cas d'un Parois mobile : A cause de la viscosité, le uide au contact de la paroi possède la même vitesse
de celle-ci →
−
v p . Mathématiquement, on écrit :

→
−
v f luide (M, t).→
−
τ =→
−
v p .→
−
τ

Telle que →
−
v f luide (M, t) est la vitesse du uide en tout point de la surface de contact avec la paroi. →
−
τ est
un vecteur unitaire tangente à la paroi en M . La vitesse du uide au contact de l'obstacle par rapport à
l'obstacle est alors tangente à l'obstacle en tout point de contact.
→
−
vp
→
− →
−
vp
dS

−
→ →
−
vp
dS

→
−
dS

Obstacle xe dans un uide en mouvement Paroi mobile tangentiellement à un uide

H.ADOCH 6 / 35 adoch.hasaan@gmail.com
CHAPTER 2

CINÉMATIQUE DES FLUIDES

C H
1 Description du mouvement d'un uide

O
1.1 Description Lagrangienne

D
Pour décrire le mouvement d'un système en mécanique du point, on a utilisé la description de Lagrange.
Cette approche consiste à suivre la particule M lors de sont mouvement.

A
−−→
Pour ce faire, on utilise le vecteur position OM qui donne les coordonnées du système à chaque instant.
L'ensemble des positions de ce système c'est ce qu'on appelle sa trajectoire. En mécanique des uides,
cette approche consiste à suivre chaque particule uide au cours de temps, ce qui n'est pas très utile car on
s'intéresse à l'ensemble de l'écoulement.
Remarques :
• Cette description ne permet pas de prendre en compte les conditions aux limites.
• Les coordonnées dépendent du temps dans cette description (x(t), y(t) et z(t) par exemple).
d→
−v
• Dans cette description, l'accélération est dénie par : →
−
a = .
dt

1.2 Description Eulérienne

Une autre approche pour décrire un uide c'est s'intéresser au champ de vitesse → −
v (M, t) dans tout l'écoulement
à chaque instant. C'est à dire qu'on n'étudie pas l'évolution de la position des particules uide mais on regarde
l'évolution du champ des vitesses dans l'espace et au cours du temps.

Dans cette description, le concept de trajectoire n'a pas de sens. Alors que le concept de ligne de champ
est plus convenable. Une ligne de champ (ou ligne de courant) est l'ensemble des points dont le champ
vitesse est tangent à t xé. On utilise aussi le concept de ligne d'émission qui est dénie comme étant une
courbe décrivant la position à un instant t0 donné de l'ensemble des particule Mi d'un uide qui sont passées
antérieurement en un même point donnée E .

7
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

→
− M1 M2 M3 M4
v M5
→
−
v

E
→
− t0
v
Ligne de courant (ou de champ) Ligne d'émission
Remarques :

H
• Les coordonnées spatiales sont indépendantes du temps.
• Cette approche tient compte des conditions aux limites.

C
∂→
−v
• La dénition de l'accélération n'est plus donnée par .
∂t

O
1.3 Dérivée particulaire
Chaque description possède donc des avantages et ses inconvénients. Dans toute la suite, on adopte la

D
description eulérienne pour étudier les uides. Le problème qui se pose c'est l'expression de l'accélération
dans cette description. Pour obtenir cette expression, on revient à la dénition de l'accélération à savoir
l'accélération d'une particule uide c'est la variation de sa vitesse par rapport au temps. Le champ de vitesse

A
→
−v (M, t) dépend de quatre variables par exemple (x, y , z et t) alors la diérentielle s'écrit :
∂→
−
v ∂→
−
v ∂→
−
v ∂→
−v
d→
−
v = dx + dy + dz + dt
∂x ∂y ∂z ∂t
Ce qu'on peut écrire sous la forme :
d→
−v ∂→
−
v ∂→
−
v ∂→
−
v ∂→
−v
= vx + vy + vz +
dt ∂x ∂y ∂z ∂t
On dénit alors l'accélération comme étant la dérivée particulaire de la vitesse qu'on écrit dans le cas général
sous la forme :
→
− D→
−
v ∂→
−v  −−→
a = + + →
−
v .grad →
−
v
Dt ∂t
Cette expression contient deux termes :
∂→
−v
∗ est la dérivée locale de →−
v qui indique un caractère non permanent de la vitesse.
∂t −−→
∗ →−v .grad →−v c'est la dérivée convective de →−
v qui indique un caractère non uniforme de la vitesse.

Remarques :
• La dérivée particulaire tient compte alors du caractère non uniforme du champ des vitesses (c'est à
dire que la vitesse change avec la position dans le uide) et du caractère non permanent de ce champ
(le champ des vitesse dépend de l'instant considéré).
• Pour une grandeur X quelconque, la dérivée particulaire s'écrit :
DX ∂X → −−→
= + −
v .grad X
Dt ∂t

H.ADOCH 8 / 35 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

• Si la grandeur X est uniforme quelque soit l'instant t, la dérivée particulaire devient :

DX ∂X
=
Dt ∂t
• Si l'écoulement est permanent alors :

D X → −−→
= −
v .grad X
Dt

H
2 Conservation de la masse

C
2.1 Débit massique

O
Soit un uide de masse volumique en écoulement. On dénit le débit massique Dm à travers une surface
Σ comme étant la masse du uide qui traverse cette surface par unité de temps :

D
dm
(kg/s) Dm =
dt −
→
dS
On dénit aussi le vecteur courant massique

A
→
− →
−
j (M, t) par : j
¨
− −
→ →
Dm = j .dS
Σ

2.2 Équation de conservation

Soit un uide de masse volumique ρ(M, t) en écoulement. Le champ de vitesse associé est noté →
−
v (M, t). Pour
déterminer l'équation de conservation de la masse eectuant un bilan de masse. Pour simplier considérons
un écoulement unidimensionnel selon l'axe OX . On considère un volume de contrôle dτ de forme cylindrique
de d'axe OX et de longueur dx et de section S .

X
x x + dx

• Bilan temporel : Déterminant la masse du uide contenant dans une portion de longueur dx d'un cylindre
ctif de section S à deux instants t et t + dt :
m(t) = ρ(x, t)dτ = ρ(x, t)Sdx et m(t + dt) = ρ(x, t)dτ = ρ(x, t + dt)Sdx

H.ADOCH 9 / 35 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

Alors la variation de la masse pendant dt est :

∂ρ ∂m ∂ρ
d2 m = m(t + dt) − m(t) = (ρ(x, t + dt) − ρ(x, t)) Sdx = S dx dt ⇒ = S dx
∂t ∂t ∂t
Cette variation est due à un ux de masse c'est-à-dire un débit massique, c'est le bilan spatial :
¨ ¨
∂m − −
→ → − −
→ →
= De − D s = j e .dS − j s .dS
∂t Σe Σs

On a un écoulement unidimensionnel donc :

→
− →
−
j e = j(x, t)→
−
ex et j s = j(x + dx, t)→
−
ex

H
Alors que :
−
→ −
→
dS e = −dS →
−
ex et dS s = dS →
−
ex

C
Ce qui donne :
¨ ¨ ¨ ¨
− −
→ → − −
→ → ∂j
j e .dS − j s .dS = j(x, t) dS − j(x + dx, t) dS = (j(x, t) − j(x + dx, t)) S = − S dx
∂x

O
Σe Σs Σe Σs

On obtient alors :
∂ρ ∂j ∂ρ ∂j
S dx = − S dx ⇒ + =0
∂t ∂x ∂t ∂x

D
C'est l'équation locale de conservation de la masse dans le cas unidimensionnel. Alors que dans le cas
général on trouve :

A
∂ρ →
−
+ div( j ) = 0
∂t

Remarques :
• La notation d2 m signie que c'est une masse élémentaire d'ordre deux (produit de deux termes petits
dx et dt).
• L'équation de conservation obtenue est analogue à l'équation de conservation de charge électrique vue
en électromagnétique II . Ce qui revient au caractère conservatif de la masse et de la charge électrique.

• Par analogie à l'électromagnétisme, on peut montrer que pour un uide de masse volumique ρ(M, t)
en écoulement, le vecteur densité de courant est donné par :

→
−
j (M, t) = ρ(M, t)→
−
v (M, t)

→
−
• En utilisant l'expression du j et les propriétés de l'opérateur div on peut montrer que l'équation de
conservation s'écrit :
Dρ
+ ρ div(→
−
v)=0
Dt

• La forme intégrale de l'équation de conservation de la masse s'écrit :

˚   ‹
∂ρ →
− dm − −
→ →
+ div( j ) = 0 ⇒ =− j .dS
∂t dt Σ

H.ADOCH 10 / 35 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

• Il ne faut pas confondre le débit massique Dm avec ledébit volumique qu'on dénit comme étant
la variation du volume du uide par unité de temps à travers une surface Σ fermé :
¨
dV →
− −
→
Dv = = v .dS (en m3 /s)
dt Σ

3 Écoulements particuliers

H
3.1 Écoulement incompressible

C
On dénit un écoulement incompressible comme étant un écoulement qui vérie :
Dρ
= 0 ce qui est équivalent à div(→
−
v)=0

O
Dt
Considérons un uide en écoulement, un tube de courant est un volume délimité par un ensemble de ligne
de courant. Sur un telle tube, le débit volumique est conservé. C'est à dire :

D
‹
→
− −
→
v .dS = 0

A
Pour le tube ci-contre, cette conservation s'écrit :

ve .Se = vs Ss

Ss
Se
→
−
ve →
−
vs

Remarques :
• Un écoulement peut être incompressible même si :
∂ρ
6= 0
∂t
• Un écoulement peut être incompressible même si la masse volumique varie avec la position dans le
uide.
• Il ne faut pas confondre un écoulement incompressible avec un uide incompressible.

H.ADOCH 11 / 35 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

3.2 Écoulement irrotationnel

Considérons un uide en écoulement, le champ des vitesses correspondant est noté →
−
v (M, t). On dénit le
vecteur tourbillon (ou vecteur vorticité) par :

→
− 1 −→ −
Ω = rot(→
v) (en s−1 )
2

On dit qu'un écoulement est irrotationnel (ou non tourbillonnaire) si le vecteur tourbillon est nul en tout
point M de l'espace. Ce qui correspond alors à un champ de vitesse qui vérie :

H
˛
−→ → →
− →
−
rot(−
v)= 0 ⇐⇒ →
−
v . dl = 0

C
Remarques :

O
• Le caractère irrotationnel d'un écoulement dépend de référentiel d'étude.
• Pour un écoulement irrotationnel qu'on appelle aussi écoulement potentiel, on peut dénir une

D
fonction appelée potentiel de vitesse notée Φ(M, t). Cette fonction est dénie par analogie au
potentiel électrostatique par la relation :

A
→
− −−→
v = −grad(Φ)

• Pour un écoulement irrotationnel incompressible, le potentiel vérie l'équation de Laplace :

−−→
div(→
−
v ) = 0 = div(−grad(Φ)) = −∆Φ ⇒ ∆Φ = 0

• Les lignes de courant sont perpendiculaire aux surfaces équipotentielles dénies par :

Φ = cte

3.3 Écoulement stationnaire

Considérons un uide en écoulement dans une conduite par exemple. On dit que l'écoulement est stationnaire
si tous les champs eulériens (→
−
v (M, t); ρ(M, t) ...etc) sont indépendants du temps.

Analysant quelques cartes de champ d'écoulements stationnaires :

• Écoulement uniforme : Il s'agit de l'écoulement le plus simple dont les particule se déplacent sans se
déformer ni tourner. Donc c'est irrotationnel et incompressible.
• Écoulement de Couette plan : Il s'agit d'un écoulement cisaillant entre un plan mobile et un plan xe.
Une portion du uide est entraînée en subissant un fort cisaillement incompressible ainsi qu'une rotation
locale donc c'est rotationnel.
• Écoulement Maelström 1 : c'est un écoulement de type tourbillonnaire mais incompressible
• Écoulement Maelström 2 : c'est un écoulement tourbillonnaire et compressible.

H.ADOCH 12 / 35 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

H
Écoulement uniforme Écoulement de Couette plan

DO C
Écoulement de Maelström 1 Écoulement de Maelström 2

A
4 Exemple d'application : Tourbillon de Rankin
Considérons une tornade, dont on propose de déterminer le champ de vitesse → −
v (M, t). On modélise cette
tornade par un cylindre de rayon a et de hauteur L >> a. L'écoulement de l'air possède alors une symétrie
cylindrique, et on considère que l'écoulement est stationnaire. On donne pour cette tornade, le vecteur
tourbillon : →
−
 Ω0 e z si a ≥ r ≥ 0

→
−
Ω (M ) =
 →
−
0 si r > a
La hauteur de cylindre est très grande donc on peut le considéré comme inni alors par invariance la vitesse
dépend seulement de r. D'autre part, Le plan (M, → −
e r, →
−
e z ) est un plan d'anti-symétrie de la tornade donc le
champ de vitesse est perpendiculaire à ce plan alors :
→
−
v (M, t) = v(r)→
−
eϕ

En utilisant le théorème de Stockes, on écrit :

˛ ¨ ˛ ¨
→
− →
− −→ → −
→ →
− − −
→ →
v . dl = rot(−
v ).dS ⇒ →
−
v . dl = 2 Ω .dS
L Σ L Σ

Telle que L est un contour fermé et Σ est la surface délimitée par ce contour.

H.ADOCH 13 / 35 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

On choisit un contour circulaire de rayon r et centré sur l'axe OZ .

Alors :
Z
˛
→
− →
−
dl = rdϕ→
−
eθ ⇒ →
−
v . dl = 2πr v(r)
L

Alors que la deuxième partie de la relation dépend de la position

du point M : M
• M à l'extérieur de la tornade : alors : r

H
¨ ˆ 2π ˆ a
− −
→ → a
2 Ω .dS = 2 Ω0 rdrdϕ = 2Ω0 πa2
Σ 0 0

C
• M à l'intérieur de la tornade : alors :
¨ ˆ 2π ˆ r
− −
→ →
2 Ω .dS = 2 Ω0 rdrdϕ = 2Ω0 πr2

O
Σ 0 0

De ces relations, on déduire que :

D
Ω r→
− si a ≥ r ≥ 0

eϕ
 0

→
−
v (M ) =

A
2
 Ω0 a →
 −
e ϕ si r > a
r

→
−
La représentation du module du champ vitesse en fonction de r montre que même si le vecteur tourbillon Ω
n'est pas continue en r = a mais la vitesse est continue.

v(r)

r
a

H.ADOCH 14 / 35 adoch.hasaan@gmail.com
CHAPTER 3

DYNAMIQUE DES FLUIDES

1 Généralités
1.1 Rappels

O C H
D
Lors de l'étude d'un uide, on le modélise par un milieu continu. Dans ce modèle, on étudie une particule
uide en utilisant la description eulérienne qui consiste à associer à l'écoulement un ensemble de champs.
Principalement, un champ de vitesse → −
v (M, t), un champ de pression p(M, t) et un champ de masse volumique

A
µ(M, t).
Dans ce cadre l'accélération de la particule uide est donnée par la dérivée particulaire de la vitesse qui
s'écrit sous la forme : →
− →
−
→
− Dv ∂v −−→ −
+ →
−
v .grad →
 
a (M, t) = = v
Dt ∂t
Un uide en écoulement vérie aussi l'équation de conservation de masse qui s'écrit dans le cas général sous
la forme :
Dµ(M, t)) ∂µ(M, t))
+ µ(M, t) div(→
−
v (M, t)) = 0 ⇐⇒ + div (µ(M, t)→
−
v (M, t)) = 0
Dt ∂t
On dénit un écoulement homogène comme étant un écoulement dont la masse volumique µ(M, t) ne
dépend pas la position.

1.2 Les forces usuels

Dans un écoulement, la particule uide est soumis -dans le cas général- à plusieurs forces à savoir :
• Le poids : Le champ de pesanteur est considéré homogène, ce qui permet d'écrire que la forme volumique
du poids est :
→
− →
−
f m = µ→
−
g car P = µdτ →
−
g
• Le champ de pression : dont l'expression est démontrée en statique du uide :
→
− −−→
f p = −grad(p)

• Force de viscosité : pour un uide incompressible on a :

→
−
f v = η ∆→
−
v

15
 CHAPTER 3. DYNAMIQUE DES FLUIDES

• Force de Laplace : pour une particule chargée :

→
− →
− →
−
f L = ρ E + ρ→
−
v ∧B

avec ρ la densité de charge volumique.

• Force d'inertie : si le référentiel est non galiléen :
→
−
f i = −µ→
−
a e − µ→
−
ac

avec →
−
a e est l'accélération d'entraînement alors que →
−
a c est l'accélération de Coriolis (Cf. mécanique du point).

H
2 Équation d'Euler - Théorème de Bernoulli

C
2.1 Équation d'Euler
Considérons un uide de masse volumique µ en écoulement telle que :

O
• L'écoulement est parfait : c'est à dire que la force de viscosité est négligeable devant les autres forces.
• L'écoulement est homogène.

D
• On néglige tous les phénomènes diusifs (diusion thermique, de particules ...etc).

Appliquons le principe de la résultante cinétique sur une particule uide de volume dτ dans un référentiel

A
galiléen lié à la conduite, ce qui donne :
 →
∂−

v −−→ → −−→
µ dτ + ( v .grad) v = µ dτ →
→
− − −
g − grad(p) dτ
∂t

D'où l'équation d'Euler :

 →
∂−

v −−→ → −−→
µ + ( v .grad) v = µ→
→
− − −
g − grad(p)
∂t

On remarque bien que pour un uide au repos, on retrouve l'équation de l'hydrostatique.

2.2 Théorème de Bernoulli

Considérons un uide en écoulement telle que :
• Le uide homogène : c'est-à-dire que la masse volumique µ ne dépend pas du point M .
• L'écoulement est parfait et incompressible. : la force de viscosité est négligeable devant les autres forces.
• L'écoulement est stationnaire : les champs caractérisant le uide (principalement le champ de vitesse) ne
dépendent pas du temps.

Le principe de la résultante cinétique s'écrit alors :

−−→ − −−→
µ (→
−
v .grad)→
v = µ→
−
g − grad(p)

Or :
−−→ − 1 −−→ − 2
−→ →
(→
−
v .grad)→
v = grad →
v + rot(−v)∧→
−
v
2

H.ADOCH 16 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

et on choisit →
−
g = g→
−
e z alors :
−−→
µ→
−
g = grad (µgz)
Le principe de la résultante cinétique s'écrit donc :
 
1 −−→ →
−→ − −−→ −−→
µ v 2 + rot(→
grad − v)∧→
−
v = grad (µgz) − grad(p)
2
L'opérateur gradient est linéaire alors et le uide est homogène donc :
 2 
−−→ v −→ −
grad µ + p + µgz = −µ rot(→ v)∧→
−
v
2

H
→
−
Alors le long d'une ligne de courant (−→−
rot(→
v)∧→
−
v = 0 ) on a :

v2

C
µ + p + µgz = cte
2

C'est le théorème de Bernoulli qu'on écrit souvent entre deux points A et B d'une ligne de courant sous

O
la forme :
vA2 v2
µ + pA + µ g z A = µ B + pB + µ g z B

D
2 2

Remarques :

A
• Pour un écoulement irrotationnel, le théorème de Bernoulli est valable pour tout point du uide.
• Attention ici l'axe (Oz) est implicitement orienté vers la verticale ascendante.
• Lorsqu'on a une dépression PB < PA , le uide accélère vB > vA et vis-versa c'est l'eet Venturi.
• Le théorème de Bernoulli est en faite une loi de conservation de l'énergie du uide. En eet :
v2
∗ µ est l'énergie cinétique volumique du uide.
2
∗ µgz est l'énergie potentielle volumique associée au poids.
∗ p c'est l'énergie volumique associée au travail des forces de pression.

2.3 Exemple d'application : eet Venturi

Considérons un uide homogène en écoulement dans une conduite cylindrique de rayon variable. On considère
que l'écoulement est incompressible, parfait et stationnaire. on relie à cette conduite un manomètre diérentiel
remplie par un liquide incompressible de masse volumique ρ0 qui permet de mesurer la diérence de pression.

Ss
Se B
D
A h
C

H.ADOCH 17 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

La relation de Bernoulli appliquée entre les points A et B qui appartient au même ligne de courant s'écrit :
vA2 v2
µ + pA + µ g z A = µ B + pB + µ g z B
2 2
D'autre part, le uide est au repos entre A et C et entre B et D alors l'équation de l'hydrostatique s'écrit :
pA + µ g z A = p C + µ g z C

et : pB + µ g zB = pD + µ g zD
Or l'application de principe fondamentale de l'hydrostatique sur la uide du manomètre donne :

H
pC + ρ0 g zC = pD + ρ0 g zD ⇒ pC − pD = ρ 0 g h

C
Ces relation permet de trouver que :
2 2gh
(pA + µgzA ) − (pB + µgzB ) = ρ0 gh − µgh ⇒ vB2 − vA2 = (ρ0 gh − µgh) = (ρ0 − µ)

O
µ µ
Pour déduire la vitesse en A ou en B il sut d'utiliser le caractère incompressible de l'écoulement qui signie
que le débit volumique est conservé alors :

D
Ss
Dv = vA Se = vB Ss ⇒ vB = vA
Se

A
La vitesse en A par exemple est donnée par :
v
u 2gh(ρ − µ)
0
vA = u
u
 2 !
u Ss
tµ −1
Se

Dans cas ci-dessus, on a Ss < Se donc on a deux cas :

• ρ0 < µ est donc la vitesse en A est plus grande que la vitesse en B .
• ρ0 > µ donc la hauteur du uide dans le manomètre doit être plus faible en B c'est-à-dire que h < 0 et

3 Équation de Navier-Stockes
3.1 Démonstration
Pour un uide newtonienne en écoulement incompressible, le principe de la résultante cinétique s'écrit :
 →
∂−

v −−→ → −−→
µ dτ + ( v .grad) v = µ dτ →
→
− − −
g − grad(p) dτ + η∆→
−
v dτ
∂t

Ce qui donne l'équation de Navier-Stockes :

 →
∂−

v −−→ → −−→
µ + ( v .grad) v = µ →
→
− − −
g − grad(p) + η∆→
−
v
∂t

H.ADOCH 18 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

Remarques :
• C'est une équation diérentielle dont l'inconnu principal c'est le champ de vitesse →
−
v (M, t).
→
− −−→ →−
• Le terme ( v .grad) v est non linéaire alors cette équation est non linéaire.
• Dans le cas général, on montre (physique des milieux continues) que pour un uide newtonien,
l'équation devient :

∂µ→
−v −−→ − −−→ η −−→
+ (→
−
v .grad)(µ→
v)=µ→
−
g − grad(p) + η∆→
−
v + grad (div (→
−
v ))
∂t 3

C H
3.2 Nombre de Reynolds
Le nombre de Reynolds est un nombre sans dimension qui caractérise l'écoulement, plus précisément, il permet
de comparer le terme d'inertie (c-à-d le transfert de quantité de mouvement par convection) et le terme de

O
viscosité (c-à-d le transfert de quantité de mouvement par diusion) dans un écoulement. C'est-à-dire :

D
−−→ −
µ(→
−v .grad)→ v
Re =
||η ∆→−
v ||

A
En utilisant des grandeurs moyennes, on trouve que :
V2
µ µV L LV
Re = L ⇒ Re = =
V η ν
η 2
L
µ La masse volumique du uide.
η La viscosité dynamique du uide.
η
ν = La viscosité cinématique du uide.
µ
V La vitesse caractéristique du uide.
L La longueur caractéristique associé à l'écoulement.

Dans une canalisation L ce sera le diamètre de la canalisation, pour un écoulement autours d'un obstacle
la largeur de l'obstacle. Ce nombre permet de déterminer le type de l'écoulement :
• Écoulement est laminaire : c'est un écoulement dont les lignes de champs sont presque parallèles.
Ceci correspond à une dominance du transport diusive c'est-à-dire un petit nombre de Reynolds (pour une
conduite cylindrique Re < 2000).
• Écoulement turbulent : c'est qu'on les lignes de champ sont chaotiques ce qui correspond à une dominance
du transport convective (terme non linéaire) c'est-à-dire un grand nombre de Reynolds (Pour une conduite
cylindrique Re > 4000).

H.ADOCH 19 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

Turbulent Laminaire

O C H
D
3.3 Exemple d'application : écoulement de Poiseuille
Considérons l'écoulement incompressible et stationnaire d'un uide de viscosité dynamique η et de masse

A
volumique ρ dans une canalisation cylindrique de rayon R et de longueur L >> R. On note OZ l'axe de
révolution de la conduite. Pour que le uide s'écoule, on impose une pression d'entrée pe en z = 0 supérieur
à la pression de sortie ps en z = L. On néglige l'eet de la pesanteur.
La symétrie cylindrique de la conduite impose que le champ de vitesse s'écrit :
→
−
v (M, t) = v(r)→
−
ez

L'équation de Navier-Stockes s'écrit dans ce cas :

−−→ − −−→
µ (→
−
v .grad)→
v = −grad(p) + η ∆→
−
v

Or →
−
v (M, t) = v(r)→
−
e z donc :
−−→ − d →
−
(→
− v = v(r) v(r)→
v .grad)→ −
ez= 0
dz
Alors l'équation de Navier-Stockes devient :
−−→
grad(p) = η ∆→
−
v

La projection sur la base cylindrique donne :

 
dp →
− dp →
− dp − 1 ∂ ∂v →
er+ eϕ+ → ez=η r −
ez
dr rdϕ dz r ∂r ∂r
Alors :  
dp dp dp 1 ∂ ∂v
= = 0 et =η r
dr dϕ dz r ∂r ∂r
Donc la pression p ne dépend ni de r ni de θ alors les deux termes de l'égalité sont des constantes :
 
dp 1 ∂ ∂v
=η r =K
dz r ∂r ∂r

H.ADOCH 20 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

Commençant par la pression p(z) qu'on intègre facilement :

dp
=K ⇒ p(z) = Kz + cte
dz
Or p(z = 0) = pe et p(z = L) = ps alors :

ps − pe
p(z) = z + pe
L

H
Alors que le champ de vitesse vérie :
   
1 ∂ ∂v ps − pe ∂ ∂v ps − pe
⇒

C
η r =K= r = r
r ∂r ∂r L ∂r ∂r ηL

Donc :
∂v ps − pe 2

O
r = r + cte
∂r 2ηL
En remplaçant par r = 0 on trouve que la constante est nulle cte = 0. alors :

D
∂v ps − pe
= r
∂r 2ηL

A
Donc :
ps − pe 2
v(r) = r +A
4ηL
La conduite est au repos donc à cause de la viscosité v(r = R) = 0 donc :
ps − pe 2
A=− R
4ηL
Donc le champ de vitesse dans le cylindre est donnée par :

pe − p s
R2 − r 2


v(r) =
4ηL

A partir de cette expression, on peut déterminer le débit volumique à travers une section. Pour ce faire,on
utilise la dénition : ¨ ˆ R ˆ 2π
Dv = − −
→ →
v .dS =
pe − ps 2 2


R −r rdrdϕ
0 0 4ηL

H.ADOCH 21 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

Ce qui donne la loi de Hagen-Poiseuille:

πR4 πR4
Dm = (pe − ps ) = ∆p
8ηL 8ηL

Par analogie à l'électrocinétique, on peut dénit une résistance hydrodynamique Rh de la conduite par :

∆p 8ηL L
Rh = ⇒ Rh = 4
= 8ηR2
Dm

H
πR S

L
Électricité Courant Tension

C
Re = ρ ρ
S
L 8η
Hydraulique Débit volumique Diérence de pression Rh = 8ηR2
S R2

O
4 Pertes de charge
4.1 Présentation

D
A cause de la viscosité du uide, un gradient de pression doit être appliqué an que le uide s'écoule dans
la canalisation. On a donc une diminution de la pression de l'amont vers l'aval de la canalisation. Cette

A
diérence de pression est la perte de charge. On appelle alors la perte de charge entre deux points A et B
d'une canalisation, la diérence de pression entre ces deux points ∆P = |pA − pB |.
On trouve deux types de pertes de charge :
• Pertes de charge singulières : C'est une perte qui correspond à des changements de géométrie se
produisant localement. Il s'agit par exemple de modications du diamètre de la conduite (rétrécissement ou
élargissement), de changements de direction (coudes), de vannes, clapets etc.
• Pertes de charge régulières : elles correspond à une diminution de la pression le long d'une conduite à
cause des petites défauts de la surface de la conduite.
Dans les deux cas la puissance perdue lors de l'écoulement est donnée par :

P = ∆p S V

telle que S est la section de la conduite et V la vitesse moyenne du uide.

4.2 Pertes de charge singulières

Les pertes de charge singulières sont dues à un changement dans géométrie de la conduite. Ces pertes sont
données par l'expression :
1
∆p = K µ V 2
2

Telle que V la vitesse moyenne de l'écoulement, µ la masse volumique et K un coecient qui dépend de la
forme de singularité.

H.ADOCH 22 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

S1 S1

S1 S1

 2  2
S1 S2
K = 1− K = 1−
S2 S1

H
Elargissement brusque Rétrécissement brusque

O C
D α α

D
R

A
  72 ! α
α D
K= 0.131 + 1.847 K = sin2 (α) + 2 sin4
π R 2

Coude arrondi Coude brusque

K = 0, 5 K = 0, 04

Entrée brusque à une canalisation Entrée progressive à une canalisation

4.3 Pertes de charge régulières

Lors d'un écoulement dans une conduite monodirectionnelle, les particules du uide rencontrent les défauts
de la surface de la conduite car elle n'est pas parfaitement lisse. Des tels défauts sont caractérisées par une
grandeur appelée rugosité absolue qu'on note généralement ε.

H.ADOCH 23 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

H
Les pertes de charge dans ce cas s'écrit sous la forme :

1 L

C
∆p = ξ µ V 2
2 d

Telle que :

O
µ est la masse volumique du uide.
V est la vitesse moyenne du uide.

D
L La longueur de la canalisation.
d Le diamètre de la canalisation.
ξ (ou f ) coecient de perte de charge ou coecient de friction ou coecient de frottement. Ce coecient

A
dépend du type de l'écoulement, en eet :
• Pour un écoulement laminaire : le coecient est donné par :

64
ξ=
Re

• Pour un écoulement turbulent : le coecient de friction dépend du nombre de Reynolds Re et de la

rugosité relative dénie par :
ε
r=
d

La valeur de ce coecient est donnée par le diagramme de Moody dont on donne un exemple ci-dessous.
Ce diagramme représente le coecient ξ en fonction du nombre de Reynolds pour diérentes valeurs de la
rugosité relative r.

H.ADOCH 24 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

O C H
5.1 PrésentationA D
5 Force de traînée et de portance

Considérons un uide en écoulement dans une conduite. A cause de la viscosité le uide au contact de la
parois (ou d'un obstacle) possède une vitesse nulle. On aura donc un gradient de vitesse, ce qui forme une
couche de transition dont la vitesse passe de façon continue de 0 à la vitesse de l'écoulement loin de l'obstacle
c'est ce qu'on appelle une couche limite. Plus précisément, la couche limite est la région de l'écoulement
d'épaisseur δ au voisinage de la paroi, où la viscosité a une inuence notable sur l'écoulement.

→
−
V∞

H.ADOCH 25 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

L'écoulement dans cette couche peut être laminaire ou turbulent selon le nombre de Reynolds de la couche
qu'on appelle nombre de Reynolds local dont l'expression est :

µV∞ δ
ReL =
η

Alors que le nombre de Reynolds de la conduite est dit nombre de Reynolds globale :

µV∞ D
Reg =

H
η

Telle que D est la taille caractéristique de la conduite (Dans le cas d'un obstacle dans l'écoulement, D est la
taille caractéristique de l'obstacle).

C
Pour évaluer l'épaisseur de la couche limite, on utilise le faite que dans cette couche la viscosité domine

O
alors que en dehors de cette couche le convection domine. A la frontière de la couche limite, on peut écrire
que :
−−→ −
µ →
−
v .grad →
v ∼ η∆→
−
 
v

D
Alors en utilisant les ordres de grandeur :

A
V∞2 V∞ D
µ =η 2 ⇒ δ=p
D δ Reg

Donc l'épaisseur de la couche limite est généralement très faible devant la taille de la conduite (ou l'obstacle).
D'autre part, on trouve que :
p
ReL = Reg

Alors même si l'écoulement est turbulent en dehors de la couche limite, il peut être laminaire dans la couche
limite. Dans la gure ci-dessous on présente quelques formes de la couche limite au voisinage d'un obstacle
sphérique selon le nombre de Reynolds.

H.ADOCH 26 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

5.2 Force de traînée et de portance

Considérons un obstacle placé dans un uide de masse volumique µ et de viscosité dynamique η en écoulement
caractérisé par une vitesse V∞ loin de l'obstacle. On constate qu'après un régime transitoire la vitesse de
l'obstacle devient constante. Il y a donc une force qui s'oppose à l'action du poids.

O C H
Cette force peut se décomposer en deux parties, la première est dirigé selon la direction de l'écoulement
c'est la force de traînée Ft et l'autre est perpendiculaire à l'écoulement c'est ce qu'on appelle la force de
portance Fp . Les expressions de ces deux forces sont:

D
1 1
Ft = µ Cx S V∞2 et Fp = µ Cz S V∞2
2 2

A
telles que µ la masse volumique du uide, V∞ la vitesse de l'écoulement loin de l'objet, S c'est le maître-couple
qui est déni comme étant la projection de l'obstacle sur la direction de l'écoulement alors que Cx et Cz sont
des coecients dépendent principalement du nombre de Reynolds Re mais aussi de l'angle d'inclinaison α de
l'objet par rapport à la direction de l'écoulement. Cette dépendance est donnée par des courbes expérimentales
dont on présente un exemple ci-dessous.

Cx

Re

Evolution de coecient Cx en fonction de Re (échelle logarithmique)

H.ADOCH 27 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

O C H
A D
Evolution de coecient Cz en fonction de l'angle d'inclinaison α

Pour le coecient de traînée Cx , on remarque alors que :

• Pour les faibles nombre de Reynolds (Re < 1) le coecient de traînée s'écrit sous la forme :

24 6
Cx = alors : Ft = S η V∞
Re D

Pour une sphérique de rayon R on a D = R et S = πR2 , alors on trouve la formule de Stockes (force
frottement visqueux (c.f mécanique)) :
Ft = 6πR η V∞

• Pour un nombre de Reynolds qui vérie (103 < Re < 105 ) le coecient de traînée est constant Cx = 0, 4
alors la force de traînée s'écrit :
1
Ft = µ Cx S V∞2 = cte V∞2
2

C'est l'expression de la force de frottement visqueux pour des fortes vitesses.

• Pour un nombre de Reynolds de l'ordre de 105 on a une chute brusque de coecient de traînée c'est ce
qu'on appelle la crise de traînée. A ce stade la couche limite, on aura une transition laminaire turbulente
de la couche limite. On assiste alors à un décollement de la couche limite qui se déplace dans la zone aval.
Des irrégularités à la surface favorisent l'apparition d'une couche limite turbulente.

H.ADOCH 28 / 35 adoch.hasaan@gmail.com
 CHAPTER 3. DYNAMIQUE DES FLUIDES

Pour le coecient de portance Cz , qui possède une grande importance dans le décollage des avions. On
remarque qu'il varie rapidement avec l'angle d'inclinaison ou d'incidence de l'aile par rapport à la direction
de l'écoulement. En eet :
• Pour un angle nulle, le coecient de portance est nulle donc la force de portance est aussi nulle. D'où
l'absence d'une force de portance pour un obstacle sphère placé dans un écoulement.
• Pour un angle inférieur à 15o , le coecient de portance augment quasi-linéairement. La traînée augmente
légèrement mais reste faible. La couche limite laminaire se décolle sur l'extrados (face supérieure de l'aile) et
remonte au fur et à mesure que l'angle α augmente. L'augmentation de la portance est essentiellement liée à
la vitesse plus élevée du uide sur l'extrados que sur l'intrados (face inférieure de l'aile). A cette diérence
de vitesse, on associe une diérence de pression.

H
• Lorsque l'angle α atteint une quinzaine de degrés se produit alors le décrochage : le coecient de portance
chute et la force de portance devient alors plus faible.

DO C
A
H.ADOCH 29 / 35 adoch.hasaan@gmail.com
CHAPTER 4

BILANS DYNAMIQUES ET THERMODYNAMIQUES

1 Présentation
1.1 Cadre d'étude

O C H
D
Considérons un uide en écoulement stationnaire dans une conduite de forme quelconque. On note R le
référentiel galiléen lié cette conduite.
Dans un premier temps, on dénit un volume de contrôle

A
Sortie
noté Σ qui est une partie du système entouré par une P s , cs
surface appelée surface de contrôle. La vitesse d'entrée
du uide est notée ce et cs celle du sortie.
As
Le caractère stationnaire permet d'écrire que tout les
grandeurs associée au volume de contrôle ne dépend pas Volume de contrôle
Σ
explicitement du temps.
m(t) = m(t + dt) ; u(t) = u(t + dt) ; s(t) = s(t + dt)
→
−
p (t) = →
− Entrée
p (t + dt) ; ec (t) = ec (t + dt) Pe , ce , Ae
telle que m la masse, u l'énergie interne massique, s l'entropie massique, →
−
p la quantité de mouvement
massique et ec l'énergie cinétique massique
Pour eectuer un bilan, il faut dénir un système fermé car le volume Σ est ouvert. Pour ce faire, on
ajoute à Σ une petite masse du uide dm. on note le nouveau système Σ∗ , ce système est fermé.

→
−
Fs
dms
Σ∗ Σ∗
dme →
− Ae
Fe

Le système Σ∗ à l'instant t Le système Σ∗ à l'instant t + dt

30
 CHAPTER 4. BILANS DYNAMIQUES ET THERMODYNAMIQUES

1.2 Bilan de masse

le système Σ est fermé alors sa masse se conserve au cours du temps donc :
m∗ (t) = m∗ (t + dt) ⇒ m(t) + dme = m(t + dt) + dms

Or le régime est stationnaire (permanent) alors m(t) = m(t + dt) donc la conservation de la masse se traduit
par :
dme = dms

Dans tout la suit on note dm = dme = dms

H
2 Bilans dynamiques

C
2.1 Bilan de quantité de mouvement

O
Appliquant le principe de la résultante cinétique sur le système Σ∗ dans le référentiel R on écrit :
d→
−
p∗ → − →
− →
−

D
= F e + F s + F conduite/f luide
dt
Or :
d→
−
p∗ →
−
p ∗ (t + dt) − →
−
p ∗ (t) (→
−
p (t + dt) + dm→
−c ) − (→
−
p (t) + dm→
−c )

A
s e
= =
dt dt dt
Le régime est permanent alors :
d→
−
p∗ dm →
= (−c e − →
−c ) = D (→
s
− →
−
m c e − c s)
dt dt
D'autre part :
→
− →
−
F e = Pe Ae →
−u e et F s = −Ps As → −
us
avec Ae est la section d'entrée et →
−
u e le vecteur unitaire entrant perpendiculaire à cette surface. Alors que As
est la section de sortie et u le vecteur unitaire sortant perpendiculaire à la surface de sortie.
→
−
s

Pour déterminer l'expression de la force exercée par la conduite sur la uide, on applique le principe de la
résultante cinétique sur cette partie de la conduite (sans le uide) qu'on considère au repos. Ce qui donne :
→
− →
− → − →
− →
−
0 = P + F f luide/conduite + F f ixation + F pression

Alors : →
− →
− →
− → − →
−
F conduite/f luide = − F f luide/conduite = P + F f ixation + F pression
On note M la masse de la conduite entourant le système Σ, le poids de cette partie de la conduite est alors :
→
−
P = M→
−
g

La force de pression correspond à la pression atmosphérique P0 qu'on considère homogène (si la conduite est
dans l'air) alors : ¨
→
− −
→
F pression = − P0 dS
L

H.ADOCH 31 / 35 adoch.hasaan@gmail.com
 CHAPTER 4. BILANS DYNAMIQUES ET THERMODYNAMIQUES

Pour déterminer l'expression de cette force latérale peut être déduire en utilisant le théorème de gradient :
‹ ˚
−
→ −−→ →
− −−→ →
−
−P0 dS = − grad(P0 )dτ = 0 car grad(P0 ) = 0
S

On choisit comme surface fermée S la surface latérale fermée par la section d'entrée et la section de sortie.
Alors: ‹ ¨ ¨ ¨
−
→ −
→ −
→ −
→ → −
− P0 dS = − P0 dS − P0 dS − P0 dS = 0
S L Ae As

Donc la force de pression est :

H
¨ ¨ ¨
→
− −
→ −
→ −
→
F pression = − P0 dS = P0 dS + P0 dS = −P0 Ae →
−
u e + P0 As →
−
us
L Ae As

C
→
−
Alors que la force de xation (si elle existe) F f ixation dépend de type de xation de la conduite.
En remplaçant les résultats obtenus dans l'expression du principe de la résultante cinétique, on obtient :

O
→
−
Dm (→
−c − →
e
−c ) = P A →
s e
− →
− →
− →
− →
−
e u e − Ps As u s + M g − P0 Ae u e + P0 As u s + F f ixation

Donc le bilan de quantité de mouvement s'écrit :

D
Dm (→
−c − →
−c ) = (P − P ) A →
− →
− − →
→ −
e s e 0 e u e + (P0 − Ps ) As u s + M g + F f ixation

A
Dans le cas d'une conduite cylindrique d'axe OZ c'est-à-dire un écoulement unidimensionnel, la projection
sur l'axe OZ donne :

Dm (ce − cs ) = (Pe − P0 ) Ae + (P0 − Ps ) As ⇒ Dm (ce − cs ) = Pe Ae − Ps Ae

2.2 Bilan de moment cinétique

Le moment cinétique du système Σ∗ par rapport à un point O dans un référentiel R est donné par :
→
− →
− −−→
• A l'instant t : L ∗O (t) = L O (t) + OM e ∧ dm → −c
e
→
− →
− −−→
• A l'instant t + dt : L ∗O (t + dt) = L O (t + dt) + OM s ∧ dm →
−c
s

La variation du moment cinétique en régime permanent est alors :

−
→ →
− →
− −−→
−c − −−→ − 
ddL∗O = L ∗O (t + dt) − L ∗O (t) = dm OM s ∧ → s OM e ∧ →
ce

Avec Me est un point à l'entrée du système et Ms est un point à la sortie du système.

Le théorème de moment cinétique appliquée au système Σ∗ par rapport au point O s'écrit :
−
→
ddL∗O −→ →
− 
= MO F ext
dt
Ce qui donne :
−−→
−c − −
Dm OM s ∧ →
−→ −  −
OM e ∧ →
→ →
− 
c e = MO F ext
s

H.ADOCH 32 / 35 adoch.hasaan@gmail.com
 CHAPTER 4. BILANS DYNAMIQUES ET THERMODYNAMIQUES

2.3 Bilan d'énergie cinétique

Appliquons le théorème d'énergie cinétique sur le système fermé Σ∗ on écrit :
dEc∗ →
− →
−
= P( F ext ) + P( F int )
dt
Examinons chaque terme seul :
   
1 2 1 2
Ec (t + dt) + dmcs − Ec (t) + dmce
dEc∗ Ec∗ (t + dt) − Ec∗ (t) 2 2 1
= = = Dm (c2s − c2e )

H
dt dt dt 2
D'autre part :
→
−

C
P( F ext ) = Ppression + Ppoids + Pu
telle que la puissance des forces de pression :

O
 
→
− − →
− − P e Ps
Ppression = F e .→
c e − F s .→
c s = Pe Ae ce − Ps As cs = Dm −
µ µ

et la puissance du poids est :

D
δWpoids dEp∗
Ppoids = =− = −Dm g (zs − ze )
dt dt

A
D'autre part, la puissance intérieure comporte deux termes à savoir :
→
− →
− →
−
P( F int ) = P( F int,pression ) + P( F int,viscosit )

Ces deux puissances dépend de la nature de l'écoulement. Alors le bilan d'énergie cinétique s'écrit :
 
1 P e Ps →
− →
−
Dm (c2s − c2e ) = Dm − − Dm g (zs − ze ) + P( F int,pression ) + P( F int,viscosit ) + Pu
2 µ µ

Remarques :
• La puissance Pu est la puissance associée aux pièces mobiles si elles existent.
• Pour un écoulement est parfait la puissance des forces de viscosité est nulle :

→
−
P( F int,viscosit ) = 0

• Pour un écoulement incompressible la puissance des forces de pression est nulle :

→
−
P( F int,pression ) = 0

• Attention, la force de pression intérieure est diérente à la force de pression extérieurs qui dépend de
la pression à l'entrée du système et celle à la sortie.

H.ADOCH 33 / 35 adoch.hasaan@gmail.com
 CHAPTER 4. BILANS DYNAMIQUES ET THERMODYNAMIQUES

3 Bilans thermodynamiques
3.1 Bilan énergétique
Le premier principe appliqué à Σ∗ s'écrit :
dE ∗ = d U ∗ + Ec∗ + Ep,ext
∗


= δW + δQ

Pour tout partie X ∗ de l'énergie totale E ∗ , on a :

dX ∗ = X ∗ (t + dt) − X ∗ (t) = X(t + dt) + dmxs − (X(t) + dmxe )

H
X est la grandeur caractéristique du système Σ. Alors que xe et xs sont les grandeurs massiques associées à
la masse dm à l'entrée et la sortie. L'écoulement est permanent donc :

C
X(t + dt) = X(t) ⇐⇒ dX ∗ = dm(xs − xe )

O
Alors en appliquant ceci on trouve :
1
dU ∗ = dm(us − ue ) ; dEc∗ = dm(ec,s − ec,e ) = dm (c2s − c2e ) ; ∗
dEp,ext = dm(ep,s − ep,e )

D
2
Ce qui donne :  
∗ 1 2 1 2

A
dE = dm (us + cs + ep,s ) − (ue + ce + ep,e )
2 2
L'énergie potentielle associée aux force extérieur, correspond généralement au poids donc :
Ep = mgz ⇐⇒ ep = gz

Alors :  
∗ 1 2 1 2
dE = dm (us + cs + gzs ) − (ue + ce + gze )
2 2
Maintenant, on se propose d'étudier les échanges énergétiques :
• Le transfert thermique qui peut s'écrire sous la forme δQ = q dm avec q est le transfert thermique
massique.
• Le travail qui se présente dans le cas général sous plusieurs formes qu'on sépare on deux : le travail des
forces de pression δWp et le travail communiqué à des pièces mobiles (turbine par exemple) placées dans le
uide δWu = wu dm. Les forces de pression sont toujours dirigées vers l'intérieur du système. On écrit alors :
− −−−→
→ → −−−→
−
δWp = F e .dOM e + Fs .dOM s = Pe S.ce dt − Ps Scs dt = Pe Ve − Ps Vs = Pe ve dm − Ps vs dm

Ou encore :
δWp = −Pe (0 − Ve ) − Ps (Vs − 0) = Pe Ve − Ps Vs = Pe ve dm − Ps vs dm
Avec Se la section d'entrée, Ss la section de sortie, Ve = ve dm le volume occupé par la masse dm à l'entrée,
Vs = vs dm le volume occupé par la masse dm à la sortie alors que ve et vs sont les volumes massiques à
l'entrée et à la sortie.

En remplaçant dans le premier principe on obtient :

 
1 1
dm (us + c2s + gzs ) − (ue + c2e + gze ) = wu dm + qdm + (Pe ve − Ps vs )dm
2 2

H.ADOCH 34 / 35 adoch.hasaan@gmail.com
 CHAPTER 4. BILANS DYNAMIQUES ET THERMODYNAMIQUES

Ce qu'on peut écrire en utilisant l'enthalpie massique h = u + P v , d'où le premier industriel :

 
1 2 1 2
(hs + cs + gzs ) − (he + ce + gze ) = wu + q
2 2

D'autre forme de ce principe là s'écrivent :

1 1
∆(h + c2 + gz) = q + wu ou Dm ∆(h + c2 + gz) = Pth + Pu
2 2

H
Remarques :
• Pth est la puissance thermique échangé avec le milieu extérieur.

C
• Pu est la puissance mécanique utile.
• Les deux sont exprimé en J/s ou W .

O
• On peut déduire la relation obtenue pour le bilan d'énergie cinétique à partir du premier principe
industriel, il sut de supposer que la transformation est isentropique ce qui signie que :

D
∆p
Pth = 0 et ∆h =
µ

A
3.2 Bilan entropique
Réalisant un bilan entropique pour le système fermé Σ∗ entre l'instant t et l'instant t + dt toujours en régime
permanent et avec les mêmes conventions d'écriture :
• A l'instant t l'entropie du système est : S ∗ (t) = S(t) + dm sentr
• A l'instant t + dt l'entropie du système est : S ∗ (t + dt) = S(t + dt) + dm ssort
• La diérence en régime permanent (S(t) = S(t + dt)) donne :

dS ∗ = S ∗ (t + dt) − S ∗ (t) = dm (ssort − sentr )

Or :  
∗ δQ q
dS = + δSc = dm + sc
Tf ront Tf ront
Donc le bilan entropique s'écrit:

q
ssort − sentr = + sc = sech + sc
Tf ront

avec sech est l'entropie massique échangée avec le milieu extérieur lors de l'écoulement, q est le transfert
thermique massique échangé avec l'extérieur alors que Tf ront est la température à la frontière du système.

H.ADOCH 35 / 35 adoch.hasaan@gmail.com