TD 2

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 3

Master STE Mécanique des fluides 2007 / 2008

Travaux dirigés de mécanique des fluides (2)


1 Retour sur le tenseur des contraintes
Revenons sur la modélisation des efforts de contact par la notion de densité surfacique de

→ →
−
force qu’on appelle contrainte notée T P, t, N . Essayons de démontrer que, pour qu’elle
soit compatible avec le principe fondamental de la dynamique et la continuité des grandeurs,

→ −

il est nécessaire que T dépende linéairement de N . Pour cela, on utilise le fait que

d−
→v
ZZZ ZZZ ZZ

→ −

ρ dV = ρ g dV + T dS
Dt dt Dt Σt

est valable pour tout ensemble matériel, donc en particulier pour la matière occupant un petit
domaine de volume δV . On peut alors considérer pour ce dernier
– un petit cylindre limité par des sections droites d’aire dS et normales unitaires extérieures

→ −

respectives N et − N , et par une surface latérale δΣL . Sa hauteur δh tend vers zéro, alors
que dS reste inchangé. Montrer que si M repère le centre d’une des sections droites, on doit
avoir nécessairement

→ →
− −
→ −

T (M, t, N ) = − T (M, t, − N ).
Quel théorème usuel de mécanique retrouve-t-on ?
– un petit tétraèdre dont trois cotés sont parallèles aux vecteurs de base du référentiel, et


dont le sommet M est repéré par − →
x . On note N la normale unitaire extérieure à la face


opposée à M , et on fait tendre vers zéro la longueur des côtés tout en conservant N . En
déduire qu’il exite nécessairement une application linéaire =

σ (−
→x , t) telle que

→− →
− −

x , t, N ) = =
T (→ ⇒
σ (−

x , t) · N .

Il s’agit dans les deux exercices suivants de résoudre les équations de Navier–Stokes pour un
fluide newtonien incompressible
 −
∂→

v 

→ −−→ −→ −−→ −
→→
ρ + v · grad v = ρ− →
g − gradp + µ ∆ − v, (1)
∂t
div−
→v = 0. (2)

en utilisant au mieux les symétries du problème, et en supposant l’écoulement stationnaire.


Dans ces équations, ρ est la masse volumique, − →v la vitesse du fluide, p sa pression, − →
g
l’accélération de la gravité, et µ la viscosité dynamique. L’expression des différents opérateurs
intervenant se trouve dans le formulaire compagnon.

1/3 [email protected]
Master STE Mécanique des fluides 2007 / 2008

2 Solutions de NS : Écoulement de Poiseuille cylindrique


Considérons l’écoulement incompressible d’un fluide visqueux newtonien dans une conduite
cylindrique

La variation de pression sur une longueur L de cylindre est notée ∆P . Cherchons une solution
stationnaire respectant les symétries des conditions aux limites, c’est-à-dire une écoulement
unidirectionnel selon z. On néglige l’action de la pesanteur.
1. Calculer le champ de vitesse à l’intérieur de la conduite.
2. En déduire le débit volumique. [Comparer ce débit avec celui obtenu avec 100 conduites
cylindriques de rayon R/10 pour un même ∆P ].
3. Calculer la force et le moment s’exerçant sur une longueur L de cylindre.

3 Solutions de NS : Écoulement de Couette cylindrique


Nous allons chercher une solution de l’équation de Navier–Stokes respectant les symétries
des conditions aux limites pour l’écoulement entre deux cylindres coaxiaux en rotation à des
vitesses angulaires constantes Ω1 et Ω2 :

Ω2 R2

Ω1 R1

Plus particulièrement, on cherche un écoulement stationnaire dans le plan perpendiculaire à


l’axe des cylindres, en négligeant les effets de la gravité.
1. Calculer le champ de vitesse entre les deux cylindres.
2. Calculer la force et le moment s’exerçant sur une longueur L du cylindre.
3. Stabilité rotatoire ⋆ : dans l’écoulement non perturbé chaque élément de fluide décrit
une trajectoire circulaire de rayon s, R1 ≤ s ≤ R2 , autour de l’axe des cylindres. Posons

σ(s) = ms2 φ̇

le moment cinétique d’un élément de fluide de masse m et de vitesse angulaire φ̇ ; la


force centrifuge, égale à σ 2 /ms3 , est équilibrée par le gradient de pression radial qui
apparaı̂t dans le fluide en rotation. L’ énergie cinétique de l’élément vaut (1/2)σ 2 /ms2 .

2/3 [email protected]
Master STE Mécanique des fluides 2007 / 2008

Pour étudier la stabilité de la solution que nous avons obtenue, on suppose qu’on échange
la position de deux éléments de fluide se trouvant initialement en s = s1 et s = s2 , ayant
la même masse m. Si l’on néglige les forces de frottement visqueux, le moment cinétique
de chaque élément fluide se conserve lors de l’échange.
(a) En raisonnant sur la différence d’énergie du système formé par ces deux éléments
entre les configurations finale et initiale, montrer qu’une condition nécessaire de
stabilité est  
2 2 1 1
(σ2 − σ1 ) 2 − 2 > 0,
s1 s2
soit encore

σ > 0.
ds
Il s’agit du critère de Rayleigh.
(b) Donner l’expression de φ̇ = φ̇(s) à l’aide des paramètres du problème. En déduire
que la condition précédente (critère de Rayleigh) peut s’écrire

(Ω2 R22 − Ω1 R12 ) > 0.

(c) Discuter les cas particuliers où les deux cylindres tournent en sens opposé, puis
dans le même sens. Étudier plus particulièrement le cas où l’un des deux cylindres
est immobile.
(d) Tracer les zones de stabilité dans le plan (Ω1 , Ω2 ). Discuter qualitativement l’effet
de la viscosité.
Lorsque l’écoulement devient instable, il apparaı̂t des tourbillons toriques, dits tourbillons de
Taylor 1 , disposés de façon régulière le long des génératrices des cylindres.
Le critère de Rayleigh est bien connu des astrophysiciens s’intéressant aux disques d’accrétion
(cf. les travaux de Pierre-Yves Longaretti au LAOG).

1
Taylor, G.I, 1923 : Stability of a viscous liquid contained between two rotating cylinders. Phil. Trans. R.
Soc. Lond. A, 223, 289–343.

3/3 [email protected]

Vous aimerez peut-être aussi