Dynamique des fluides

Analyse de la quantité de
mouvement

Bovic Kilundu 2
 Introduction: 3 équations de conservation en
mécanique des fluides
1. Conservation de masse
𝑑 𝑚𝑠𝑦𝑠𝑡
▪ Pour un système fermé subissant un changement: 𝑚𝑠𝑦𝑠𝑡 = 𝐶𝑡𝑒 ou =0
𝑑𝑡
𝑑𝑚𝐸𝑉
▪ Pour un élément de volume: 𝑚𝑖𝑛
ሶ − 𝑚ሶ 𝑜𝑢𝑡 =
𝑑𝑡
2. Conservation de la quantité de mouvement
▪ Pour un corps rigide: 𝑝Ԧ = 𝑚𝑐Ԧ
𝑑 𝑝Ԧ
▪ Seconde loi de Newton: 𝐹 = 𝑚𝑎Ԧ ou ∑𝐹 = donc 𝑝Ԧ = 𝑐𝑡𝑒 si ∑𝐹Ԧ = 0
Ԧ Ԧ
𝑑𝑡
3. Conservation de l’énergie
▪ Transfert d’énergie pour un système fermé: sous forme de chaleur et de travail
▪ Conservation:
′
𝐸𝑛𝑒𝑟𝑔𝑖𝑒 𝑡𝑜𝑡𝑎𝑙𝑒 𝑒𝑛𝑡𝑟𝑎𝑛𝑡𝑒 − 𝐸𝑛𝑒𝑟𝑔𝑖𝑒 𝑡𝑜𝑡𝑎𝑙𝑒 𝑠𝑜𝑟𝑡𝑎𝑛𝑡𝑒 =
𝑉𝑎𝑟𝑖𝑎𝑡𝑖𝑜𝑛 𝑑𝑒 𝑙 é𝑛𝑒𝑟𝑔𝑖𝑒 𝑡𝑜𝑡𝑎𝑙𝑒 𝑑𝑢 𝑠𝑦𝑠𝑡è𝑚𝑒
𝑑𝐸𝑠𝑦𝑠𝑡
➢ Soit 𝐸𝑖𝑛 − 𝐸𝑜𝑢𝑡 = Δ𝐸𝑠𝑦𝑠𝑡 ሶ ሶ
ou 𝐸𝑖𝑛 − 𝐸𝑜𝑢𝑡 = pour un élément de
𝑑𝑡
volume
Prof Bovic Kilundu 3
 Lois de Newton
𝑑 𝑚𝑐Ԧ
Ԧ
▪ Seconde loi: 𝐹 = 𝑚𝑎Ԧ = . En mécanique des fluides, la seconde
𝑑𝑡
loi de Newton est appelée Équation linéaire de mouvement
▪ Pour les solides en rotation:
𝑑𝜔 𝑑 𝐼𝜔 d𝐻
𝑀 = 𝐼 𝛼Ԧ = 𝐼 = = , 𝐻 est le moment angulaire
𝑑𝑡 dt 𝑑𝑡
▪ Pour un corps solide en rotation autour d’un axe 𝑥:
𝑑𝜔𝑥 𝑑𝐻𝑥
𝑀𝑥 = 𝐼𝑥 =
𝑑𝑡 𝑑𝑡

Bovic Kilundu 4
 Forces agissant sur un élément de volume
▪ Les forces qui agissent sur un EV comprennent:
1. Les forces de volume (gravité ou électromagnétique): d𝐹Ԧ𝑔𝑟𝑎𝑣𝑖𝑡é = 𝜌𝑔𝑑𝑉Ԧ
o Vecteur gravité: 𝑔Ԧ = −𝑔𝑘
o Force de volume totale: ∑𝑣𝑜𝑙 𝐹Ԧ𝑣𝑜𝑙 = ‫𝑉𝑑𝑔 𝜌 𝑉𝐸׬‬
Ԧ = 𝑚𝐸𝑉 𝑔Ԧ
2. Les forces de surface
o Pas si faciles à analyser  existence des composantes normales et tangentielles
o Tenseur des contraintes σ𝑖𝑗 pour décrire les contraintes de surface en un point d’écoulement

𝜎𝑥𝑥 𝜎𝑥𝑦 𝜎𝑥𝑧

𝜎𝑖𝑗 = 𝜎𝑦𝑥 𝜎𝑦𝑦 𝜎𝑦𝑧
𝜎𝑧𝑥 𝜎𝑧𝑦 𝜎𝑧𝑧

o Composantes diagonale= contraintes normales, dépendent de la pression et des contraintes visqueuses

o Les composantes hors diagonales sont appelées contraintes de cisaillement , essentiellement constituées
des contraintes de viscosité

Bovic Kilundu 5
▪ 𝜎𝑖𝑗 : contrainte dans la direction j dans un plan
perpendiculaire à la direction i

▪ Orientation positive des contraintes:

▪ Le produit scalaire d’un tenseur et d’un vecteur donne

un vecteur → produit contracté ou produit interne
▪ Le produit interne du tenseur des contraintes 𝜎𝑖𝑗 avec
le vecteur normal 𝑛 d’un petit élément de surface
donne la force par unité de surface:
𝑑 𝐹Ԧ𝑠𝑢𝑟𝑓 = 𝜎𝑖𝑗 . 𝑛 𝑑𝐴 ⇒ ∑𝐹Ԧ𝑠𝑢𝑟𝑓 = න 𝜎𝑖𝑗 . 𝑛 𝑑𝐴
𝐸𝑆

Bovic Kilundu 6
▪ Force totale sur l’EV:
∑𝐹Ԧ = ∑𝑣𝑜𝑙𝑢𝑚𝑒 𝐹Ԧ + ∑𝑠𝑢𝑟𝑓𝑎𝑐𝑒 𝐹Ԧ = ∑𝐹Ԧ𝑔𝑟𝑎𝑣𝑖𝑡é + ∑ 𝐹Ԧ𝑠𝑢𝑟𝑓𝑎𝑐𝑒 = ‫𝑔𝜌 𝑉𝐸׬‬Ԧ 𝑑𝑉 + ‫𝐴𝑑𝑛 𝑗𝑖𝜎 𝑆𝐸׬‬

Cette équation est très utile pour l’analyse différentielle de la quantité de

mouvement. Mais pour l’analyse pratique, on peut utiliser l’expression
suivante:

∑𝐹Ԧ = ∑𝐹Ԧ𝑔𝑟𝑎𝑣𝑖𝑡é + ∑𝐹Ԧ𝑝𝑟𝑒𝑠𝑠𝑖𝑜𝑛 + ∑𝐹Ԧ𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑡é + ∑𝐹Ԧ𝑎𝑢𝑡𝑟𝑒𝑠

Force totale Forces de Forces

volume surfaciques

Bovic Kilundu 7
▪ ∑𝐹Ԧ𝑎𝑢𝑡𝑟𝑒𝑠 : forces de réaction nécessaires pour faire tourner
l’écoulement; forces sur les boulons, les câbles, les supports ou les
parois que rencontre l’écoulement.

Clé de robinet

Bovic Kilundu 8
 Equation linéaire de la quantité de
mouvement
▪ Comme la vitesse et la densité peuvent changer d’un point à un autre
de l’écoulement, la seconde loi de Newton s’écrit:
𝑑
∑𝐹Ԧ = න 𝜌𝑐Ԧ 𝑑𝑉
𝑑𝑡 𝑠𝑦𝑠𝑡
▪ Cette équation est écrite pour un élément de masse alors qu’en
mécanique des fluides on travaille avec des éléments de volume.
▪ La formulation pour un EV est:

𝑑 𝑚𝑐Ԧ 𝑠𝑦𝑠𝑡 𝑑
= න 𝜌 𝑐Ԧ 𝑑𝑉 + න 𝜌𝑐Ԧ 𝑐Ԧ𝑟 . 𝑛 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝐸𝑉 𝐸𝑆

Bovic Kilundu 9
Soit
𝑑
∑𝐹Ԧ = න 𝜌 𝑐Ԧ 𝑑𝑉 + න 𝜌𝑐Ԧ 𝑐Ԧ𝑟 . 𝑛 𝑑𝐴
𝑑𝑡 𝐸𝑣 𝐸𝑆

Somme de Variation dans le temps du Vitesse d’écoulement de la

toutes les forces changement de la quantité quantité de mouvement en
extérieures de mouvement du contenu sortie de l’élément de surface
agissant sur l’EV de l’EV par écoulement de masse

𝑐Ԧ𝑟 = 𝑐Ԧ − 𝑐Ԧ𝐸𝑆 : vitesse relative du fluide par rapport à l’ES

Pour un EV fixe, sans mouvement ni déformation), 𝑐Ԧ𝑟 = 𝑐Ԧ

Bovic Kilundu 10
▪ Pour un EV fixe:
𝑑
∑𝐹Ԧ = න 𝜌 𝑐Ԧ 𝑑𝑉 + න 𝜌𝑐Ԧ 𝑐. Ԧ 𝑛 𝑑𝐴
𝑑𝑡 𝐸𝑉 𝐸𝑆
▪ L’équation de la quantité de mouvement est une
équation vectorielle → on considèrera ses
composantes suivant x, y et z en coordonnées
cartésiennes
▪ Dans la plupart des cas ∑𝐹=Ԧ somme du poids, des
forces de pression et des forces de réaction.

▪ Cette équation est souvent utilisée pour déterminer

les forces induites par un écoulement.
Bovic Kilundu 11
 Cas particuliers

▪ Ecoulement permanent:
∑𝐹Ԧ = න 𝜌𝑐Ԧ 𝑐Ԧ𝑟 . 𝑛 𝑑𝐴
𝐸𝑆
▪ On peut utiliser les vitesses moyennes pour éviter
l’utilisation des intégrales

න 𝜌𝑐Ԧ 𝑐Ԧ𝑟 . 𝑛 𝑑𝐴𝑐 = 𝜌𝑐𝑚𝑜𝑦 𝐴𝑐 . 𝑐Ԧ𝑚𝑜𝑦 = 𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦

𝐴𝑐

Bovic Kilundu 12
 Facteur de correction du flux de quantité
de mouvement: 𝜷
▪ Le profil des vitesses à travers les entrées et les sorties n’est pas
uniforme → on corrige la quantité de mouvement calculée avec la
vitesse moyenne en utilisant le facteur 𝛽 ( coefficient de Boussinesq)

𝑑
∑𝐹Ԧ = න 𝜌 𝑐Ԧ 𝑑𝑉 + ෍ 𝛽𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦 − ෍ 𝛽 𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦
𝑑𝑡 𝐸𝑉
𝑠𝑜𝑟𝑡𝑖𝑒 𝑒𝑛𝑡𝑟é𝑒

Le coefficient de Boussinesq est défini tel que

න 𝜌𝑐Ԧ 𝑐.
Ԧ 𝑛 𝑑𝐴𝑐 = 𝛽𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦
𝐴𝑐
Bovic Kilundu 13
▪ Si la masse volumique est uniforme sur toute l’entrée ou la sortie et
𝑐Ԧ est dans la même direction que 𝑐Ԧ𝑚𝑜𝑦 , on résout pour avoir 𝛽:

‫𝑐𝜌 𝐴׬‬Ԧ 𝑐.
Ԧ 𝑛 𝑑𝐴𝑐 ‫𝑐𝜌 𝐴׬‬Ԧ 𝑐.
Ԧ 𝑛 𝑑𝐴𝑐
𝑐 𝑐
𝛽= =
𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦 𝜌𝑐𝑚𝑜𝑦 𝐴𝑐 𝑐𝑚𝑜𝑦

On peut simplifier pour avoir

2
1 𝑐
𝛽= න 𝑑𝐴𝑐
𝐴𝑐 𝐴𝑐 𝑐𝑚𝑜𝑦
On peut montrer que 𝛽 est toujours supérieur à 1
Bovic Kilundu 14
 Ecoulement permanent
𝑑
∑𝐹Ԧ = න 𝜌 𝑐Ԧ 𝑑𝑉 + ෍ 𝛽𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦 − ෍ 𝛽 𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦
𝑑𝑡 𝐸𝑉
𝑠𝑜𝑟𝑡𝑖𝑒 𝑒𝑛𝑡𝑟é𝑒

▪ Si l’écoulement est permanent, on a

∑𝐹Ԧ = ෍ 𝛽 𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦 − ෍ 𝛽 𝑚ሶ 𝑐Ԧ𝑚𝑜𝑦
𝑠𝑜𝑟𝑡𝑖𝑒 𝑒𝑛𝑡𝑟é𝑒
▪ Ecoulement permanent à une seule entrée et une seule sortie:
∑𝐹Ԧ = 𝑚ሶ 𝛽2 𝑐Ԧ2 − 𝛽1 𝑐Ԧ1

Bovic Kilundu 15
 La relation ∑𝐹Ԧ = 𝑚ሶ 𝛽2 𝑐Ԧ2 − 𝛽1 𝑐Ԧ1 est
vectorielle

Bovic Kilundu 16
 Exemples
▪ Travailler les exemples 6-1 à 6-7 dans Çengel , p. 248-257

Bovic Kilundu 17
 Exemple
Un coude réducteur est utilisé pour dévier de
30° un écoulement d'eau à un débit de 14 kg/s
dans un tuyau horizontal vers le haut tout en
l'accélérant. Le coude rejette l'eau dans
l'atmosphère. La section transversale du
coude est de 113 cm2 à l'entrée et de 7 cm2 à
la sortie. La différence d'altitude entre les
centres de la sortie et de l'entrée est de 30 cm.
Le poids du coude et de l'eau qu'il contient est
considéré comme négligeable. Déterminez (a)
la pression manométrique au centre de
l'entrée du coude et (b) la force d'ancrage
nécessaire pour maintenir le coude en place.

Bovic Kilundu 18
 Travail étudiants
▪ Lecture de la section 6.5

Bovic Kilundu 19
Approche différentielle

Bovic Kilundu 20
 Introduction
▪ Forces agissant sur une particule fluide
1. Forces de volume: 𝒇𝒃 comme le poids (𝜌𝐠)
𝐃𝒄
2. Forces d’inertie: 𝜌
𝐷𝑡
3. Forces de surface ou pressions 𝒇𝒔
▪ Seconde loi de Newton:
𝐷𝒄
𝜌 =𝒇
𝐷𝑡

𝒇 = 𝒇𝒃 + 𝒇𝒔 : force par unité de volume

Prof Bovic Kilundu 21

▪ La force de volume est dite conservative si:
𝒇𝒃 = 𝛁G

G est le potentiel scalaire

▪ Si la gravité est la seule force de volume, on a
𝐟𝐛 = ρ𝐠
𝐠 = 𝒈: force de gravité par unité de masse

Bovic Kilundu 22
 Forces de surface
▪ Considérons la force dans la direction 𝑥1
▪ Chacune des 6 face a une composante dans cette
direction
+
𝜕𝜎11 𝑑𝑥1
𝜎11 = 𝜎11 + .
𝜕𝑥1 2

−
𝜕𝜎11 𝑑𝑥1
𝜎11 = 𝜎11 − .
𝜕𝑥1 2

Les faces normales à 𝑥1 contribuent pour une force:

𝜕𝜎11
. 𝑑𝑥1 𝑑𝑥2 𝑑𝑥3
𝜕𝑥1
Bovic Kilundu 23
▪ En incluant la contribution des autres faces, la force de surface totale
par unité de volume suivant 𝑥1 est:

𝜕𝜎11 𝜕𝜎21 𝜕𝜎31

𝑓𝑠1 = + +
𝜕𝑥1 𝜕𝑥2 𝜕𝑥3
▪ L’équation de la quantité de mouvement dans la direction 𝑥1 :
𝐷𝑢1 𝜕𝜎11 𝜕𝜎21 𝜕𝜎31
𝜌 = + + + 𝜌𝑓1
𝐷𝑡 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3
▪ On peut généraliser cette expression à d’autres direction, en notation
indicielle on obtient
𝐷𝑢𝑗 𝜕𝜎𝑖𝑗
𝜌 = + 𝜌𝑓𝑗
𝐷𝑡 𝜕𝑥𝑖
Bovic Kilundu 24
▪ Pour d’autres systèmes de coordonnées:
D𝐜
= 𝛁. 𝛔 + ρ𝐟
Dt
▪ En coordonnées cartésiennes, les trois composantes de cette équations
sont:
𝐷𝑢1 𝜕𝜎11 𝜕𝜎21 𝜕𝜎31
𝜌 = + + + 𝜌𝑓1
𝐷𝑡 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3

𝐷𝑢2 𝜕𝜎12 𝜕𝜎22 𝜕𝜎32

𝜌 = + + + 𝜌𝑓2
𝐷𝑡 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3

𝐷𝑢3 𝜕𝜎13 𝜕𝜎23 𝜕𝜎33

𝜌 = + + + 𝜌𝑓3
𝐷𝑡 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3
Bovic Kilundu 25
 𝐷𝑢𝑗 𝜕𝜎𝑖𝑗
𝜌 = + 𝜌𝑓𝑗
𝐷𝑡 𝜕𝑥𝑖
𝜕𝜌 𝜕 𝜌𝑢𝑖
▪ En multipliant l’équation de continuité + = 0 par 𝑢𝑖 et en
𝜕𝑡 𝜕𝑥𝑖
additionnant avec l’équation de la quantité de mouvement ci-dessus:

𝜕 𝜌𝑢𝑗 𝜕 𝜌𝑢𝑖 𝑢𝑗 𝜕𝜎𝑖𝑗

+ = + 𝜌𝑓𝑗
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑖
Ou
𝜕 𝜌𝒄
+ 𝛁. 𝜌𝒄𝒄 = 𝛁. 𝝈 + 𝜌𝒇
𝜕𝑡
C’est une autre forme de l’équation de la quantité de mouvement appelée
équation de Cauchy
▪ Le terme 𝒄𝒄 est le produit dyadique, il représente un tenseur de deuxième
ordre

Bovic Kilundu 26
 Equation de Cauchy en coordonnées
cartésiennes
𝜕 𝜌𝒄
+ 𝛁. 𝜌𝒄𝒄 = 𝛁. 𝝈 + 𝜌𝒇
𝜕𝑡
𝜕 𝜌𝑢1 𝜕 𝜌𝑢1 𝑢1 𝜕 𝜌𝑢2 𝑢1 𝜕 𝜌𝑢3 𝑢1 𝜕𝜎11 𝜕𝜎21 𝜕𝜎31
+ + + = 𝜌𝑓1 + + +
𝜕𝑡 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3

𝜕 𝜌𝑢2 𝜕 𝜌𝑢1 𝑢2 𝜕 𝜌𝑢2 𝑢2 𝜕 𝜌𝑢3 𝑢2 𝜕𝜎12 𝜕𝜎22 𝜕𝜎32

+ + + = 𝜌𝑓2 + + +
𝜕𝑡 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3

𝜕 𝜌𝑢3 𝜕 𝜌𝑢1 𝑢3 𝜕 𝜌𝑢2 𝑢3 𝜕 𝜌𝑢3 𝑢3 𝜕𝜎13 𝜕𝜎23 𝜕𝜎33

+ + + = 𝜌𝑓3 + + +
𝜕𝑡 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3 𝜕𝑥1 𝜕𝑥2 𝜕𝑥3
Bovic Kilundu 27
Equations de Navier-Stokes

Bovic Kilundu 28