Mécanique des fluides Partie 2
Partie 2
- Équations de Navier-Stokes :
- Rappels (tenseur des contraintes) ;
- Équations de Navier-Stokes ;
- Problématique concernant l’accélération ;
- Loi de transport (théorème de Reynolds).
- Couche limite & écoulements autour d’obstacles
- Notion de couche limite ;
- Écoulements autour d’obstacles ;
- Écoulement autour d’un cylindre ;
- Écoulement autour d’un profil d’aile d’avion.
EIGSI – Casablanca
Mécanique des fluides Partie 2
Rappels : tenseur des contraintes
Forces appliquées sur Ω :
S
- Force à distance : sur le volume Ω ; S
- Forces de contact : sur la surface SΩ.
Forces de contact :
- Forces de contact normales à SΩ ;
- Forces de contact tangentielles à SΩ.
T (M ,n) = (M )n
Ti = ij n j
11 12 13
= 21 22 23
31 32 33
12 = 21 ; 13 = 31 ; 23 = 32
EIGSI – Casablanca 2
1
Mécanique des fluides Partie 2
Équations de Navier-Stokes :
- Équations aux dérivées partielles non linéaires.
- Description du mouvement des fluides :
- Modélisation des écoulement autour d’obstacles : bâtiments, avions, voitures, trains… ;
- Modélisation des masses d’aire de l’atmosphère (météorologie) ;
- Modélisation des courants océaniques ;
-…
Principe fondamentale de la dynamique :
Fext = m a
Comment
Fsurfacique + Fvolumique = m a calculer
l’accélération??
Peuvent être exprimées
Gravité dans la
on utilisant le tenseur de
majorité des cas
contraintes
EIGSI – Casablanca 3
Mécanique des fluides Partie 2
Problématique concernant l’accélération :
- La description Eulérienne est la plus adaptée à l’étude du comportement des fluides.
- L’accélération est une grandeur Lagrangienne par définition.
Écoulement permanent:
A B
U B = cste t
U A = cste t
Ne représente pas la variation de la
U
=0 vitesse d’une particule mais la variation
t de la vitesse en un point fixe.
EIGSI – Casablanca 4
2
Mécanique des fluides Partie 2
Dérivée particulaire :
dU U U U U
= +Ux +U y +Uz
dt t x y z
Dérivée temporelle Dérivée convective
(variation locale) (variation spatiale)
dU U U
dt
=
t
(
+ U .grad U =
t
) (
+ U . U )
U x U x U x U x
t + U x x + U y y + U z z
dU U y U y U y U y
En coordonnées cartésiennes : = +Ux +U y +Uz
dt t x y z
U z U z U z U z
t + U x x + U y y + U z z
EIGSI – Casablanca 5
Mécanique des fluides Partie 2
Principe fondamentale de la dynamique :
m a = Fvolumique + Fsurfacique
dU
dt
dV = g dV + n dS
S
Théorème de Green-Ostrogradski
(théorème de la divergence)
dU
dt dV = g dV + div dV
U
+ (U . )U = g + .
t
EIGSI – Casablanca 6
3
Mécanique des fluides Partie 2
Tenseur des contraintes :
- Le tenseur des contraintes peut être décomposé en deux parties :
- partie relative aux forces de pression ;
- partie relative aux effets visqueux.
= p + v = − pI + v
Pression hydrostatique : existe même
Effets visqueux
s’il n’y a pas de mouvement
v = 2 D
− p 0 0
p = 0 − p 0 Loi de comportement de Navier-Stokes
0 0 − p pour un fluide newtonien en
écoulement incompressible
1 U U j
Dij = i +
2 x j xi
Tenseur de vitesse de déformation
EIGSI – Casablanca 7
Mécanique des fluides Partie 2
Équations de Navier-Stokes :
U
+ (U . )U = g − p + U
2
t
En coordonnées cartésiennes :
U x U x U x U x p 2U x 2U x 2U x
+Ux +U y +Uz = g x − + 2 + 2 + 2
t x y z x x y z
U y U y U y U y p U y U y 2U y
2 2
+ U + U + U = g − + 2 + 2 + 2
x z
t x y z y
x y z y
y
U z + U U z + U U z + U U z = g − p + U z + U z + U z
2 2 2
t
x 2 2
x y z z z
x y z z
y
2
EIGSI – Casablanca 8
4
Mécanique des fluides Partie 2
Loi de transport : (théorème de Reynolds)
La loi de transport consiste à écrire le taux de variation à l’instant t d’une quantité scalaire ou vectorielle
définie sur un volume matériel quelconque ( t ) délimité par sa surface ( t ) d’un milieu fluide en
écoulement.
V2 ( t )
V1 ( t )
Soit f ( x, t ) une fonction définie sur le volume matériel (t ) de surface (t ) .
Le taux de variation de l’intégrale de f ( x , t ) sur ( t ) est:
d
f ( x , t ) dV
dt (t )
EIGSI – Casablanca 9
Mécanique des fluides Partie 2
Par définition de la dérivée :
d 1
f ( x , t ) dV = lim f ( x , t + t ) dV − f ( x , t ) dV
t →0 t
dt (t ) (t +t ) ( t )
On a :
V ( t ) volume balayé par S1
( t + t ) = ( t ) − V1 ( t ) + V2 ( t ) avec 1
V2 ( t ) volume balayé par S2
Donc :
f ( x , t + t ) dV +
d
1 ( t )
f ( x , t + t ) dV
V2 ( t )
f ( x , t ) dV = lim
dt (t )
t →0 t
− f ( x , t + t ) dV − f ( x , t ) dV
V1 ( t ) ( t )
On a :
1 f ( x , t )
lim f ( x , t + t ) dV − f ( x , t ) dV = dV
t → 0 t (t ) ( t ) (t ) t
EIGSI – Casablanca 10
5
Mécanique des fluides Partie 2
On a :
1
lim f ( x , t + t ) dV − f ( x , t + t ) dV
t → 0 t V2 ( t ) V1 ( t ) ( )
dV = U . n t dS
1
= lim
t → 0 t S2
( ) ( )
f ( x , t + t ) U . n t dS + f ( x , t + t ) U . n t dS
S1
( )
= lim f ( x , t + t ) U . n dS + f ( x , t + t ) U . n dS
t → 0
( )
S2 S1
= f ( x , t ) (U . n ) dS = ( f ( x , t )U ) dV
( t ) ( t )
Théorème de Reynolds :
f ( x , t )
d
f ( x , t ) dV =
t
( )
+ . f ( x , t ) U dV
( t )
dt (t )
Taux de variation de f sur Variation de f par rapport à t Flux de f à travers la surface
EIGSI – Casablanca 11
Mécanique des fluides Partie 2
Notion de couche limite :
Couche de faible épaisseur qui se développe à l’interface d’un corps solide.
a
a
Le bord
d’attaque
Zone 1: écoulement laminaire. Zone laminaire: l’épaisseur de la Ux
Rex =
Zone 2: écoulement transitoire. couche limite augmente en fonction de x.
Zone 3: écoulement turbulent, Zone turbulente: l’épaisseur de la
avec une sous-couche visqueuse couche limite augmente plus vite que Nombre de Reynolds à
(laminaire). dans la zone laminaire. une abscisse donnée
EIGSI – Casablanca 12
6
Mécanique des fluides Partie 2
Épaisseur de la couche limite :
Écoulement sur une plaque plane de longueur l.
a
a
Couche limite laminaire: Couche limite turbulente:
5, 48 0,37
= x = x
Rex Re1x 5
EIGSI – Casablanca 13
Mécanique des fluides Partie 2
Écoulement d’un fluide à l’entrée d’une conduite :
Le
Longueur d’entrée : (longueur d’établissement)
Régime laminaire: (Re<2000) Régime transitoire et turbulent: (Re>2000)
Le Le
0, 0065 Re 4, 4 Re1 6
D D
EIGSI – Casablanca 14
7
Mécanique des fluides Partie 2
Écoulements autour d’obstacles :
Bernoulli entre a et b :
U 2
p0 = p + la pression d ' arrêt
2
EIGSI – Casablanca 15
Mécanique des fluides Partie 2
Forces agissant sur un obstacle : (profil d’aile)
Force de l’écoulement du fluide agissant sur l’aile : Trainée :
Tp = p dS cos
- Effort de pression (normal) ;
S
- Effort de frottement (tangentiel). T = Tp + T f
T f = dS sin
F =T +L
S
Portance :
Lp = − p dS sin
T : Traînée.
S
L L = Lp + L f
L f = dS cos
: Portance.
S
EIGSI – Casablanca 16
8
Mécanique des fluides Partie 2
Exemples de traînée de pression et de frottement :
EIGSI – Casablanca 17
Mécanique des fluides Partie 2
Écoulement autour d’un cylindre : (fluide visqueux)
- La couche limite se sépare en deux parties
symétriques.
- Deux tourbillons se forment derrière le cylindre.
Re 47
- Oscillation périodique du sillage.
- Les tourbillons se détachent alternativement du
bas et du haut du cylindre.
Allées de Von Karman
47 Re 3,8 105
EIGSI – Casablanca 18
9
Mécanique des fluides Partie 2
- Oscillation de plus en plus intense.
- Point de décollement se déplace en amont du
cylindre.
Allées de Von Karman
47 Re 3,8 105
La couche limite devient turbulente
Re 3,8 105
EIGSI – Casablanca 19
Mécanique des fluides Partie 2
Écoulement autour d’un profil d’aile d’avion :
- La traînée et la portance sont difficiles à obtenir théoriquement sauf pour les objets de forme simple.
- Des formules empiriques (expérimentales) sont utilisées.
1 CD : coefficient de traînée.
T= U 2 CD S S : surface projetée sur le plan vertical à la direction de
2
l’écoulement (surface frontale).
1 CL : coefficient de portance.
L= U 2 CL S S : surface projetée sur le plan horizontal à la direction
2
de l’écoulement (surface portante ou surface alaire).
EIGSI – Casablanca 20
10
Mécanique des fluides Partie 2
Exemple : (A380-800)
A380-800
Longueur : 72 m
Envergure : 80 m
Hauteur : 24 m
Poids max : 560 t
Surface alaire : 845 m2
EIGSI – Casablanca 21
Mécanique des fluides Partie 2
Coefficient de portance en fonction de l’angle d’attaque :
Décrochage
(perte de portance)
EIGSI – Casablanca 22
11
Mécanique des fluides Partie 2
Angle de décrochage :
- Début de décollement de la couche limite sur l’extrados (face supérieure de l’aile) à partir du bord
de fuite.
- Élargissement et déplacement du décollement vers le bord d’attaque.
- Couche limite turbulente.
EIGSI – Casablanca 23
Mécanique des fluides Partie 2
Application : (écoulement de Couette plan)
On considère l’écoulement plan entre deux plaques distantes de h dont l’une est mobile à la vitesse U et
l’autre est fixe. La largeur des parois est grande devant h ainsi que la longueur notée L.
Le fluide a un comportement newtonien et l’écoulement est incompressible et permanent. On néglige le
poids du fluide devant les forces de viscosité.
1- Montrer que le champ de vitesse est de la forme : U = (U x ( y ) , 0, 0 ) .
2- A partir des équations de Navier-Stokes, montrer que le champ de pression est uniforme en tout point
et égal à la pression atmosphérique p0 et que Ux ( y) = U y h .
3- Calculer le débit massique de l’écoulement.
EIGSI – Casablanca 24
12
Mécanique des fluides Partie 2
Bibliographie :
Mécanique expérimentale des fluides – Tome 1 : Statique et dynamique des fluides non visqueux –
Raymond Comolet – Dunod.
Mécanique expérimentale des fluides – Tome 2 : Dynamique des fluides réels, turbomachines –
Raymond Comolet – Dunod.
Mécanique des fluides appliquée – Roger Ouziaux & Jean Perrier – Dunod.
EIGSI – Casablanca 25
13
