Cours de Mécanique des fluides
Equations de la
mécanique des
fluides
Olivier LOUISNARD
Forces extérieures
Poids : g dV
V Déjà spécifiées pour l’hydrostatique
Pression : -pn dS
S
Mais dans un fluide en mouvement, il y a aussi
des frottements visqueux
• frottement fluide / fluide
• adhérence fluide aux parois solides
• dissipent de l’énergie
• origine microscopique : mouvement thermique + interactions dan
les liquides
Viscosité : expérience de Couette
y
t1 t2 > t1 t U0
0
v
h v v
x
Constatations expérimentales :
• u = U0 sur la plaque supérieure
•u=0 sur la plaque inférieure
• profil linéaire de u au bout d’un temps assez grand
• avec S surface mouillée
Le coefficient de proportionnalité ne dépend que du fluide
= viscosité dynamique
Viscosité
• homogène à kg.m-1.s-1 = Pa.s = Pl (Poiseuille)
on utilise le Poise (Po) et surtout le Centipoise (cPo)
Eau à 20°C : 10-3 Pa.s = 1 cPo
Air à 20°C : 1.85 10-5 Pa.s
• augmente avec T pour un gaz
indépendant de p pour un gaz
diminue avec T pour un liquide (cf. huile dans une poêle)
augmente avec p pour un liquide
Contrainte visqueuse
(Rappel MMC) : Contrainte = dF/dS
S
dFv=v dS
n dFp= -pn dS
dS
V n p= force de pression/u. de surface
dS
v= force visqueuse/u. de surface
Contrainte de pression p = -pn y
Exprimée facilement en fonction de n yx
Question : peut-on exprimer v en fonction de n ?
exv
Oui sous forme tensorielle xx
x
zx
v = = .n z
Pour les fluides
On montrera (Cours 7) : dits « newtonien »
Bilan de quantité de mouvement
S dFv=v. n dSMaintenant qu’on connaît toutes
dFp= -pn dS les forces extérieures,
n dS on peut écrire le bilan de QDM :
V dV dS n
vdV= - v (v.n) dS
dFp = gdV V S
vdV= v (v.n) dS + g dV -pn dS + v. n dS
- S V S S
V
Variation de QDM QDM Poids Pression Frottement
du fluide transportée visqueux
dans le volume V par le fluide
rentrant - sortant
A RETENIR
Bilan d’énergie
Le calcul de inclut les puissances de
toutes les forces extérieures
• Forces de pression : dFp = pn dS d p= pn.v dS
• Forces visqueuses : dFv = vn dS d v = ( vn).v dS
• Poids : dFg = g dV d g = g.v dV
S
dFv=vn dS
n dFp= -pn dS
dS
V n
dV dS
v
v dFp = gdV
Bilan d’énergie
u + v2/2)dV u + v2/2)(v.n) dS
V S
Variation d’énergie totale Energie totale transportée
du fluide par le fluide
dans le volume V rentrante - sortante
g.v dV -pv n dS vn) v dS
V S S
Puissance du Puissance des Puissance des Puissance
poids forces de frottements calorifique
pression visqueux
Equations locales
Objectif : remplacer le bilan sur un volume V
par des relations différentielles valables en chaque point du fluide
Moyens mathématiques :
• théorèmes analyse vectorielle
• puis passage à la limite V 0
Intérêt :
• calcul analytique
• calcul numérique
Un exemple
Conservation de la masse :
==>
Vrai quel que soit V donc :
Equations locales de la mécadef
Masse masse volumique
pression
énergie interne
QDM
vitesse
Energie
Système complet ? 1 équation vectorielle 1 inconnue vectorielle
2 équations scalaires 3 inconnues scalaires
+ 2 équations scalaires + 1 inconnue scalaire
Il manque une équation d’état :
Quelques équations d’état
Gaz parfait : (compresseurs, turbines à gaz)
BAROTROPES
GP isotherme : (rare) Equation de l’énergie
découplée de M et QDM
(acoustique, ondes de chocs,
GP isentropique : écoulements gazeux en général)
(explosions sous-marines,
Liquide écoulements liquides
compressible : supersoniques, rare)
(hydraulique, presque tous
Fluide les écoulements liquide +
incompressible : écoulements gaz faible Mach)
Autres écritures QDM
s’écrit aussi
ou encore
= a accélération du fluide
Des modèles pour simplifier
Ces équations sont des EDP très
complexes
On cherche donc des approximations à l’aide
d’hypothèses supplémentaires
• Modèle de fluide incompressible : masse volumique constante
• Modèle de fluide parfait : frottement visqueux négligés
Dans les TD à suivre, on utilisera en général les deux
On parlera ensuite de la validité ...
Fluide incompressible
(x,y,z,t) = M= dV = V Constante
V
Conservation de la masse :
S = Se + Ss
Se Ss
v
Général v.n dS = 0 n dS
S V n
v
veSe = vsSs Ce qui rentre = Ce qui sort
Tube de courant Accumulation de masse impossible
Ve Vs
V = vS
Equation locale débit volumique
(noté aussi Q)
Validité fluide incompressible
Correct si : Ma nombre de Mach
c vitesse du son dans le fluide
déduite de l’équation d’état
Exemple pour un gaz parfait: = 340 m/s à 298 K
Validité indépendante du caractère gazeux ou liquide
nutilisable si Ma > 0,3
nutilisable pour rendre compte de certains phénomènes (acoustique, choc
En pratique presque toujours valable dans les liquides
Modèle de fluide parfait
Permet de négliger les frottements visqueux
• mouvement non dissipatif
• conservation de l’énergie mécanique
• pas d’adhérence aux parois solides : le fluide
« glisse »
• ouvre de nombreuses simplifications mathématiques
Limitations évidentes. Ne rend pas compte :
• du freinage visqueux d’un corps ou d’un fluide (voiture économique !)
• de l’amortissement des ondes (vagues, acoustiques, ...)
• de la nécessité de pomper un fluide
Validité ?
Validité fluide parfait
Nombre de Reynolds
• Ecoulements externes :
Si Re >> 1, valable à l’extérieur de la couche limite (qui est petite)
(cf. cours 9)
Si Re << 1, totalement invalide,
à traiter par théorie écoulements rampants.
(cf. cours 7, 8)
• Ecoulements en conduite :
Fluide parfait applicable (Bernoulli) ...
avec correction pour pertes de charges
(cf. cours 6)
Conservation QDM en fluide parfait
Forme globale
vdV = v (v.n) dS + g dV -pn dS
V
- S V S
Variation de QDM QDM transportée Poids Pression
du fluide par le fluide
dans le volume V rentrant - sortant
Forme locale Equations d’Euler
Fluide parfait incompressible
Incompressible Parfait
Equations locales :
Masse
QDM
Une grande simplification
est possible :
Loi de Bernoulli
Loi de Bernoulli : démonstration
: // v => (rot v v) . dM = (v dM) . rot v = 0
On suit une ligne de courant dM
v 2
v
M
1 dM v v
On suppose régime permanent =>
On projette la conservation QDM sur la ligne de courant
( grad v2/2 + rot v v) . dM = (g grad p) . dM
De plus, on peut écrire g = grad (gz) si z orienté vers le haut
grad( v2/2 + p + gz) . dM = 0
v12/2 + p1 + gz1 = v22/2 + p2 + gz2
Loi de Bernoulli : énoncé
Sous les hypothèses :
• Fluide parfait
• Fluide incompressible Il existe une version en
• Régime permanent compressible
Peut être généralisé en
instationnaire dans quelques cas
rares (cf. TD)
La quantité p + v2/2 + gz est constante le long d’une ligne de courant
Energie Energie cinétique Energie
potentielle
potentielle de
de pesanteur
pression
Traduit la conservation de l’énergie mécanique
A retenir pour les TDs fluide parfait.
Fluide parfait
incompressible :
Formule de
Bernoulli :
p + v2/2 + gz = Cte le long d’une ligne de courant
Ecoulements (démontré ultérieurement)
unidirectionnels
Dans la direction transverse à un écoulement unidirectionnel,
la pression varie de façon hydrostatique.
Conditions aux limites en pression :
Aux points de contact entre un écoulement et l’atmosphère,
la pression vaut patm.