La mécanique des uides

H
PSI

DO C
A
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
0 Statique des uides

H
Force surfacique - force volumique . . . . .
Types de forces . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
1
1
1

D
1.2 Le poids . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.3 Force de pression . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
2 Équation fondamentale de l'hydrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

A
2.1 Démonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2 Fluide incompressible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.3 Fluide compressible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3 Action sur un corps solide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.1 Action sur les parois . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.2 Poussé d'Archimède . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1 Étude phénoménologique des uides 8
1 Modélisation d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.1 Dénition d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2 Modèle de milieu continu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 Caractéristiques d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1 Pression d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2 Notion de viscosité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.3 Condition aux limites cinématiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2 Cinématique des uides 12
1 Description du mouvement d'un uide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.1 Description Lagrangienne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.2 Description Eulérienne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3 Dérivée particulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2 Conservation de la masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.1 Débit massique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2 Équation de conservation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3 Écoulements particuliers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.1 Écoulement incompressible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.2 Écoulement irrotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

i
 TABLE DES MATIÈRES

3.3 Écoulement stationnaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

4 Exemple d'application : Tourbillon de Rankin . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3 Dynamique des uides 19
1 Généralités . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.1 Rappels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.2 Les forces usuels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2 Équation d'Euler - Théorème de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.1 Équation d'Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.2 Théorème de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20

H
2.3 Exemple d'application : eet Venturi . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3 Équation de Navier-Stockes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.1 Démonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

C
3.2 Nombre de Reynolds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.3 Exemple d'application : écoulement de Poiseuille . . . . . . . . . . . . . . . . . . . . . 23
4 Pertes de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

O
4.1 Présentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.2 Pertes de charge singulières . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
4.3 Pertes de charge régulières . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

D
5 Force de traînée et de portance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
5.1 Présentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
5.2 Force de traînée et de portance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

A
4 Bilans dynamiques et thermodynamiques 33
1 Présentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
1.1 Cadre d'étude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
1.2 Bilan de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2 Bilans dynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.1 Bilan de quantité de mouvement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.2 Bilan de moment cinétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.3 Bilan d'énergie cinétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3 Bilans thermodynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.1 Bilan énergétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.2 Bilan entropique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

H.ADOCH ii / 38 adoch.hasaan@gmail.com
CHAPTER 0

STATIQUE DES FLUIDES

C H
1 Force surfacique - force volumique
1.1 Types de forces

O
D
Lorsqu'un corps est soumis à une force, elle peut être :
• Ponctuelle: c'est-à-dire qu'elle est appliquée en un point M (Exemples : tension d'un l ou ressort).
• Surfacique : c'est-à-dire qu'elle est appliquée en une surface (Exemples : pression). On dénit la force

A
surfacique ou la densité surfacique de cette force par :
→
−
2 →
− dF
(N/m ) fs=
dS
• Volumique : c'est-à-dire qu'elle est appliquée en tout point de volume du corps (Exemples : coulomb,
gravitation). On dénit la force volumique ou la densité volumique de cette force par :
→
−
3 →
− dF
(N/m ) fv=
dτ
Dans un uide, on trouve plusieurs forces de nature diérentes, on se propose d'étudier quelques exemples.

1.2 Le poids
Le poids est dû principalement à l'interaction gravitationnelle entre la terre et le corps étudié. Cette interaction
est une force volumique qui s'applique sur tout point du corps. Considérons un uide (liquide ou gaz) de
masse dm et de volume dτ , son poids est donné par :
→
−
d P = dm →
−
g = ρdτ →
−
g

ρ est la masse volumique du uide. La force volumique associée au poids est donc :

→
−
fv=ρ→
−
g

1
 CHAPTER 0. STATIQUE DES FLUIDES

1.3 Force de pression

Pour la force de pression c'est une force surfacique donc l'expression est donnée par :
→
− −
→
d F = P dS
−
→ →
−
dS = dS →
−
n est le vecteur dont le module dS est égale à la surface soumise à la force d F et sa direction →
−
n est
perpendiculaire sur la surface et dans le sens de la force.
La force surfacique associée à cette force est alors :
→
−
→
− dF

H
fs= = P→
−
n
dS
Maintenant, déterminant l'équivalent volumique de cette force. Pour ce faire, on considère un petit volume

C
dτ = dxdydz d'un uide. Étudiant la force de pression appliquée par le reste du uide sur ce cube.
Y
La force de pression appliquée sur la face 1 est :

O
→
−
d F x,1 = P (x + dx, y, z)dydz(−→
−
e x)

Pour la face 2, on a : dy

D
2 1
→
−
d F x,2 = P (x, y, z)dydz(+→
−
e x)

A
X
La force totale selon OX est donc : dz
→
−
d F x = − (P (x + dx, y, z) − P (x, y, z)) dydz →
−
ex Z dx

Le développement limité à l'ordre 1 s'écrit :

∂P
P (x + dx) = P (x) + dx
∂x
Alors la force totale selon OX devient :
→
− ∂P ∂P →
dF x = − dxdydz →
−
ex=− dτ −
ex
∂x ∂x
De même on trouve dans les deux autres directions :
→
− ∂P → →
− ∂P →
dF y = − dτ −
ey df z = − dτ −
ez
∂y ∂z
La force totale appliquée par le reste du uide sur cette cube est la somme des trois forces ce qui donne :
 
→
− ∂P →
− ∂P →
− ∂P →
−
dF = − ex+ ey+ e z dτ
∂x ∂y ∂z

La force volumique associée à la force de pression est alors :

→
−  
→
− dF ∂P →
− ∂P →
− ∂P →
− →
− −−→
fv= =− ex+ ey+ ez ⇒ f v = −grad(P )
dτ ∂x ∂y ∂z
−−→
grad est l'opérateur vectoriel gradient dont l'expression dépend du système de coordonnées :

H.ADOCH 2 / 38 adoch.hasaan@gmail.com
 CHAPTER 0. STATIQUE DES FLUIDES

• En cartésienne :
−−→ ∂P →
− ∂P →
− ∂P →
−
grad(P ) = ex+ ey+ ez
∂x ∂y ∂z
• En cylindrique :
−−→ ∂P →
− 1 ∂P →
− ∂P →
−
grad(P ) = er+ eϕ+ ez
∂r r ∂ϕ ∂z
• En sphérique :
−−→ ∂P →
− 1 ∂P →
− 1 ∂P →
−
grad(P ) = er+ eθ+ eϕ
∂r r ∂θ rsin(θ) ∂ϕ

H
2 Équation fondamentale de l'hydrostatique

C
2.1 Démonstration
Considérons une particule M uide de masse dm et de volume dτ d'un uide au repos dans un référentiel
R supposé galiléen. Cette particule est soumis à deux forces à savoir le poids et la force de pression. Le

O
principe de la résultante cinétique (PFD pour un corps non ponctuel) appliqué à cette particule s'écrit :
→
− →
−
dm→
−

D
a (M/R) = d P + d F p
→
−
La particule est au repos donc →
−
a (M/R) = 0 alors :

A
→
− →
− →
−
dP + dF p = 0

D'après les paragraphes précédents :

−−→
ρdτ →
−
g − grad(p)dτ = 0
D'où la loi local de l'hydrostatique :
−−→
grad(p) = ρ→
−
g

Remarques :
• C'est un système de trois équations diérentielles de premier ordre. Donc la résolution donne le champ de
pression P (M ) dans le uide.
• Lors de la résolution, il faut faire attention à l'accélération de pesanteur qui dépend à priori de la position.
• La nature du uide, permet de déterminer la relation entre la masse volumique ρ et la pression P .

2.2 Fluide incompressible

Considérons un liquide (supposé incompressible) de masse volumique ρ placé dans une enceinte au repos.
Le champ de pesanteur est supposé uniforme. La loi de l'hydrostatique s'écrit :
−−→
grad(P ) = ρ→
−
g

La projection sur OZ donne :

dP
= ρg
dz

H.ADOCH 3 / 38 adoch.hasaan@gmail.com
 CHAPTER 0. STATIQUE DES FLUIDES

La masse volumique des liquides est pratiquement indépendante →

−
g
de la pression donc :
ˆ
P (z) = ρgdz

le liquide est homogène donc ρ ne dépend pas de z donc :

ˆ Liquide
P (z) = ρgdz = ρgz + cte

H
Ce qu'on peut écrire :

C
P (z) = P (z0 ) + ρg(z − z0 )
Z

O
La pression augmente alors avec la profondeur. Par exemple dans l'eau dont la masse volumique ρ =
1 kg/L = 103 kg/m3 la pression augmente par 1Bar = 105 P a chaque 10 mètre de profondeur.

D
2.3 Fluide compressible
Pour les gaz, la masse volumique dépend de la pression. Cette dépendance est selon la nature de gaz. On

A
se propose d'étudier un exemple classique à savoir le modèle de l'atmosphère isotherme dont lequel l'air est
modélisé par un gaz parfait de température homogène (T = T0 = cte). On choisit l'axe OZ vers le haut.
La loi de l'hydrostatique s'écrit :
−−→ Z →
−
g
grad(P ) = ρ→
−
g

La projection sur OZ donne :

dP
= −ρg
dz
La masse volumique est relié à la pression par l'équation des gaz Air
parfait :
m PM
PV = RT ⇒ P M = ρRT ⇒ ρ=
M RT
Donc l'équation de l'hydrostatique devient :
dP PM
=− g
dz RT
C'est une équation diérentielle de première ordre qu'on peut réécrire sous la forme :
dP Mg
=− dz
P RT
Une intégration de cette équation donne :
ˆ   ˆ
dP P Mg
= Ln = − dz
P P0 RT

H.ADOCH 4 / 38 adoch.hasaan@gmail.com
 CHAPTER 0. STATIQUE DES FLUIDES

La température est homogène (atmosphère isotherme) et l'accélération de pesanteur est considéré homogène
donc :  
P Mg
Ln =− z + cte
P0 RT0
Donc la pression est donnée par :
   
Mg z − z0
P (z) = P0 exp − (z − z0 ) = P0 exp −
RT0 δ

H
P0 est pression correspondante à l'altitude de référence z0 . Généralement on choisit z0 = 0 au niveau de la
mer, la pression est alors P0 = 1 atm.
Remarques :

C
• La grandeur δ est homogène à une longueur, en faite c'est la longueur caractéristique de l'atténuation de la
pression. Cette distance de l'ordre de 10km.
• Pour une augmentation d'altitude ∆z = δ la pression se divise par e.

O
• La fonction exponentielle contient en faite un rapport de l'énergie potentielle molaire associée au poids
(Mgz) et l'énergie cinétique associée à l'agitation thermique (RT). Un telle facteur est appelée facteur de
Boltzmann.

D
• A haute température l'énergie de l'agitation thermique est grande devant l'énergie donc la variation de
pression avec l'altitude est faible.

A
• Une autre interprétation de ce facteur, vient de l'expression de la densité des particules qui est donnée par :
 
N NA NA M gz
P V = nRT0 = RT0 ⇒ ∗
n (z) = P (z) = P0 exp −
NA RT0 RT0 RT0

Ce qu'on peut écrire sous la forme :

 
Mg
∗
n (z) = n∗0 exp − (z)
RT0

Alors le nombre de particules par unité de volume dépend aussi de l'altitude par un facteur de Boltzmann.
En eet, si la température est très basse la majorité de particules seront au sol (z ' 0) alors que lorsqu'on
augmente la température l'agitation thermique permet aux particules d'accéder à des altitudes plus grandes.

P n∗
P0 T2 > T1 n∗0 T2 > T1
T1 > T0 T1 > T0
T0 T0

z z

Evolution de la pression Evolution du nombre de particule

H.ADOCH 5 / 38 adoch.hasaan@gmail.com
 CHAPTER 0. STATIQUE DES FLUIDES

3 Action sur un corps solide

3.1 Action sur les parois
Lorsqu'on on met un uide dans une enceinte, celui la applique une force de pression sur chaque paroi de
récipient. Pour la calculer, on utilise la relation :
¨
→
− −
→
F = (P − P0 )dS

Par exemple pour la force appliquée sur la surface perpendiculaire M

H
sur OX on a :
¨
→
− H2 →
ρgzdydz →
− −

C
F = e x = ρgL ex X
2
La force de pression est une fonce surfacique mais on peut
la modélisée par une force ponctuelle. Pour ce faire, il faut

O
O
déterminer donner la force résultante et le point d'application C
appelée centre de poussé.
Le calcul de la force résultante ce fait par la relation précédente alors que le centre de poussé est déterminer

D
en utilisant le moment de la force :
ˆ ¨

A
−→ → − −−→ →
− −−→
OC ∧ F = OM ∧ d F = OM ∧ ρgz dS

O un point quelconque généralement le bas de la surface alors que M est un point quelconque de la surface.

3.2 Poussé d'Archimède

Lorsqu'on met un solide dans un uide, il est soumis à deux forces : son poids et la pression du liquide, la
résultante de la force de pression donne une force appelé poussé d'Archimède dont l'expression est admise

→
−
F a = −ρl Vemerge →
−
g

Avec :
ρl est la masse volumique du uide.
Vemerge est le volume de solide émergé dans le uide.

Exemple d'application :
Considérons un parallélépipède de longueur a, de largeur a, de hauteur h > a et de masse volumique ρs placée
verticalement dans un liquide de masse volumique ρl . Déterminons la hauteur du parallélépipède émergé
hemerg dans l'eau.

H.ADOCH 6 / 38 adoch.hasaan@gmail.com
 CHAPTER 0. STATIQUE DES FLUIDES

Le solide a son équilibre est soumis à deux forces dont les

expressions sont : →
−
g
→
− →
−
P = m→
−
g F a = −ρl Vemerge →
−
g
h
Le solide est à l'équilibre donc :
→
− → − →
−
P + Fa= 0 ⇒ m = ρl Vemerg

En utilisant la dénition de la masse on obtient :

H
m = ρs V = ρl Vemerg

Or V = a2 h et Vemerg = a2 hemerg donc :

C
ρs
hemerg = h
ρl

O
On a donc trois cas :
• ρs = ρl dans ce cas h = hemerg , le parallélépipède est émergé totalement dans l'eau.
• ρs ≤ ρl alors hemerg < h donc le parallélépipède est émergé partiellement dans l'eau. si hemerg << h une

D
partie très grande otte.
• ρs > ρl dans ce cas h < hemerg , ce qui signie que le parallélépipède est émergé totalement dans l'eau mais
n'est pas en équilibre. L'accélération du corps lorsqu'il est émergé totalement dans l'eau est donnée par le

A
principe fondamentale de la dynamique:
→
− → − ρs − ρl →
m→
−
a = P + F a = (ρs − ρl )V →
−
g ⇒ →
−
a = −
g
ρs

H.ADOCH 7 / 38 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 Modélisation d'un uide
1.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.

1.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 à
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
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τ
du uide de masse xe dont la longueur typique d ∼ 3 dτ est très grande devant le libre parcours moyen l
et très faible devant la longueur caractéristique L de l'écoulement.

l << d << L

Dans la majorité des cas, lp 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.

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

2 Caractéristiques d'un uide

2.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

H
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

C
aussi :
1 bar = 105 P a ; 1 atm = 101 325 P a
1

O
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 :

D
la première c'est que le volume du uide étudié est généralement dicile à préciser. Alors que la deuxième
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 :

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

2.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
mouvement du uide.
Les expériences montrons que lorsque un uide est en écoulement
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
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
−→ ∂vx →
−
dF v = η dS ex →
−
∂y v (y)

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.
• 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

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

l'inverse aussi).
• 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.
• 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 :
η
ν=
ρ

H
• On donne un ordre de grandeur de coecient de viscosité dynamique, pour quelques uides :
Fluide Air Eau Glycérine Miel Graisse

C
η à 25 C (en P l) 10−5 10−3
o
1,5 10 103

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

O
écoulement incompressible :
→
−
d F v = η ∆→
−
v dτ

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

A
div(→
−
v)=0

2.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
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 :
• 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 :

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

−
→
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 : la même contrainte précédente permet d'écrire que 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).dS = →
v p .dS

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

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

→
−
vp

−
→
dS

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

DO C
A
H.ADOCH 11 / 38 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. 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

A
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 .

→
− M1 M2 M3 M4
v M5
→
−
v

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

12
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

Remarques :
• Les coordonnées spatiales sont indépendantes du temps.
• Cette approche tient compte des conditions aux limites.
∂→
−v
• La dénition de l'accélération n'est plus donnée par .
∂t

1.3 Dérivée particulaire

Chaque description possède donc des avantages et ses inconvénients. Dans toute la suite, on adopte la
description eulérienne pour étudier les uides. Le problème qui se pose c'est l'expression de l'accélération

H
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
→
−v (M, t) dépend de quatre variables par exemple (x, y , z et t) alors la diérentielle s'écrit :

C
→
− ∂→
−
v ∂→
−
v ∂→
−
v ∂→
−v
dv = dx + dy + dz + dt
∂x ∂y ∂z ∂t

O
Ce qu'on peut écrire sous la forme :
d→
−v ∂→
−
v ∂→
−
v ∂→
−
v ∂→
−v

D
= 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

A
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
• 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.ADOCH 13 / 38 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

2 Conservation de la masse
2.1 Débit massique
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 :
dm
(kg/s) Dm =
dt −
→
dS
On dénit aussi le vecteur courant massique

H
→
− →
−
j (M, t) par : j
¨

C
− −
→ →
Dm = j .dS
Σ

O
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
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 .

A
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

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

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

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

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
Σe Σs Σe Σs ∂x
On obtient alors :
∂ρ ∂j ∂ρ ∂j
S dx = − S dx ⇒ + =0
∂t ∂x ∂t ∂x
C'est l'équation locale de conservation de la masse dans le cas unidimensionnel. Alors que dans le cas
général on trouve :

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

Remarques :

C
• La notation d2 m signie que c'est une masse élémentaire d'ordre deux (produit de deux termes petits dx et
dt).

O
• 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

D
écoulement, le vecteur densité de courant est donné par :

→
−

A
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 Σ

• Il ne faut pas confondre le débit massique Dm avec le dé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
3.1 Écoulement incompressible
On dénit un écoulement incompressible comme étant un écoulement qui vérie :
Dρ
= 0 ce qui est équivalent à div(→
−
v)=0
Dt
H.ADOCH 15 / 38 adoch.hasaan@gmail.com
 CHAPTER 2. CINÉMATIQUE DES FLUIDES

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 :
‹ Ss
→
− −
→
v .dS = 0 Se
→
−
ve →
−
vs
Pour le tube ci-contre, cette conservation s'écrit :

H
ve .Se = vs Ss

Remarques :

C
∂ρ
• 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.

O
3.2 Écoulement irrotationnel

D
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 :

A
→
− 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 :
˛
−→ → →
− →
−
rot(−
v)= 0 ⇐⇒ →
−
v . dl = 0

Remarques :
• 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 fonction
appelée potentiel de vitesse notée Φ(M, t). Cette fonction est dénie par analogie au potentiel électrostatique
par la relation :
→
− −−→
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

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

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

H
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.

DO C
A
Écoulement uniforme Écoulement de Couette plan

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

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

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

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.
On choisit un contour circulaire de rayon r et centré sur l'axe OZ .

H
Alors :
Z
˛
→
− →
−
dl = rdϕ→
− →
−

C
eθ ⇒ v . dl = 2πr v(r)
L

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

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

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

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

De ces relations, on déduire que :

Ω r→
− si a ≥ r ≥ 0

eϕ
 0

→
−
v (M ) = 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 18 / 38 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))
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

19
 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 20 / 38 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 −−→ →
−→ − −−→ −−→
µ grad −
v 2 + rot(→
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 :

C
v2
µ + 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

Remarque :

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.
La relation de Bernoulli appliquée entre les
points A et B qui appartient au même ligne
de courant s'écrit :
vA2 v2 Ss
µ + pA + µ g z A = µ B + pB + µ g z B
2 2 Se B
D'autre part, le uide est au repos entre D
A et C et entre B et D alors l'équation de
l'hydrostatique s'écrit : A h

pA + µ g z A = pC + µ g z C C

et : pB + µ g z B = pD + µ g z D

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

Or l'application de principe fondamentale de l'hydrostatique sur la uide du manomètre donne :

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

Ces relation permet de trouver que :

2 2gh
(pA + µgzA ) − (pB + µgzB ) = ρ0 gh − µgh ⇒ vB2 − vA2 = (ρ0 gh − µgh) = (ρ0 − µ)
µ µ
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 :

H
Ss
Dv = vA Se = vB Ss ⇒ vB = vA
Se

C
La vitesse en A par exemple est donnée par :

O
v
u 2gh(ρ − µ)
0
vA = u
u
 2 !
u Ss
tµ −1

D
Se

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

A
• ρ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

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.
• L'équation de Navier-Stockes n'a pas de solution analytique exacte.

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

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
viscosité (c-à-d le transfert de quantité de mouvement par diusion) dans un écoulement. C'est-à-dire :

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

H
En utilisant des grandeurs moyennes, on trouve que :
V2
µ

C
µV L LV
Re = L ⇒ Re = =
V η ν
η 2
L

O
µ La masse volumique du uide.
η La viscosité dynamique du uide.

D
ν La viscosité cinématique du uide.
V La vitesse caractéristique du uide
L La longueur caractéristique associé à l'écoulement.

A
Turbulent
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 Laminaire
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).

3.3 Exemple d'application : écoulement de Poiseuille

Considérons l'écoulement incompressible et stationnaire d'un uide de viscosité dynamique η et de masse
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

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

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 :

H
 
dp →
− dp →
− dp − 1 ∂ ∂v →
er+ eϕ+ → ez=η r −
ez
dr rdϕ dz r ∂r ∂r

C
Alors :  
dp dp dp 1 ∂ ∂v
= = 0 et =η r
dr dϕ dz r ∂r ∂r

O
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

D
dz r ∂r ∂r
Commençant par la pression p(z) qu'on intègre facilement :
dp

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

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

Alors que le champ de vitesse vérie :

   
1 ∂ ∂v ps − pe ∂ ∂v ps − pe
η r =K= ⇒ r = r
r ∂r ∂r L ∂r ∂r ηL
Donc :
∂v ps − pe 2
r = r + cte
∂r 2ηL
En remplaçant par r = 0 on trouve que la constante est nulle cte = 0. alors :
∂v ps − pe
= r
∂r 2ηL
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

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

A partir de cette expression, on peut déterminer le débit volumique à travers une section. Pour ce faire,on

H
utilise la dénition : ¨ ˆ R ˆ 2π
Dv = − −
→ →
v .dS =
pe − ps 2 2


R −r rdrdϕ
0 0 4ηL

C
Ce qui donne la loi de Hagen-Poiseuille:

πR4 πR4

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

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

∆p 8ηL 2L

A
Rh = ⇒ Rh = = 8ηR
Dm πR4 S

L
Électricité Courant Tension Re = ρ ρ
S
L 8η
Hydraulique Débit volumique Diérence de pression Rh = 8ηR2
S R2

4 Pertes de charge
4.1 Présentation
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
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.

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

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
S1 S1

S1 S1

O C
 2  2
S1 S2
K = 1− K = 1−
S2 S1

D
Elargissement brusque Rétrécissement brusque

A
D α α

  72 ! α
α D 2 4
K= 0.131 + 1.847 K = sin (α) + 2 sin
π R 2

Coude arrondi Coude brusque

K = 0, 5 K = 0, 04

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

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

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 ε.

C H
L

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

O
1 L
∆p = ξ µ V 2

D
2 d

Telle que :

A
µ est la masse volumique du uide.
V est la vitesse moyenne du uide.
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
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 de nombre de Reynolds et de 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 27 / 38 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 28 / 38 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 29 / 38 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 30 / 38 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 31 / 38 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 32 / 38 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

33
 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 34 / 38 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 35 / 38 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 36 / 38 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 37 / 38 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 .
• On peut déduire la relation obtenue pour le bilan d'énergie cinétique à partir du premier principe industriel,

O
il sut de supposer que la transformation est isentropique ce qui signie que :
∆p
Pth = 0 et ∆h =

D
µ

3.2 Bilan entropique

A
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 38 / 38 adoch.hasaan@gmail.com