Mécanique des fluides

Correction TD

A. El Khalfi 1
 Correction Exercice 1
On étudie l'écoulement de fluide visqueux
incompressible (de caractéristiques r et m)
qui se produit lorsque l'on rapproche à
vitesse constante 𝑢 = −𝑈 𝑒𝑧 un disque
mobile d'un disque fixe de même rayon R et
de même axe (cf. figure).
On suppose que le fluide à l'extérieur des
disques est à une pression p=p0 uniforme et
qu’on vérifie : h/R<<1 .
D'autre part, on néglige l'influence des
forces de masse.
On choisit le repère d'étude le référentiel
cylindrique.

A. El Khalfi 2
 Hypothèses du problème

Le repère de travail est le repère relatif : 𝑒𝑟 , 𝑒𝜃 , 𝑒𝑧

Variables du problème : 𝑟, 𝜃, 𝑧

Symétrie cylindrique :
𝜕
𝜕𝜃

Ecoulement plan (r,z): 𝑉𝜃 = 0

Inconnues du problème : 𝑝 𝑟, 𝑧, 𝑡 , vr (𝑟, 𝑧, 𝑡), vz (𝑟, 𝑧, 𝑡)

A. El Khalfi 3
 Equation du problème

Equation de continuité : ∂vr + ∂vz + vr = 0

∂r ∂z r
Equations de Navier-Stokes (en coordonnées cylindriques :
∂vr ∂vr ∂vr 1 ∂p μ ∂2 vr ∂2 vr 1 ∂vr vr
+ vr + vz =− + + 2 + − 2
∂t ∂r ∂z ρ ∂r ρ ∂r 2 ∂z r ∂r r

∂vz ∂vz ∂vz 1 ∂p μ ∂2 vz ∂2 vz 1 ∂vz

+ vr + vz =− + + 2 +
∂t ∂r ∂z ρ ∂z ρ ∂r 2 ∂z r ∂r
Conditions aux limites :
Adhérence au disque inférieur immobile :
vr r, z = 0, t = 0 vz r, z = 0, t = 0
Adhérence au disque supérieur mobile :
vr r, z = h, t = 0 vz r, z = h, t = −U
A l’extérieur des disques, la pression est : 𝑝(𝑟 = 𝑅, 𝑧, 𝑡) = 𝑝
0
A. El Khalfi 4
 Analyse de grandeurs
Echelles des longueurs : 𝑜 𝑟 = 𝑅 ,𝑜 𝑧 = ℎ
𝜕 1 𝜕 1
Echelles des opérateurs : 𝑜 = ,𝑜 =
𝜕𝑟 𝑅 𝜕𝑧 ℎ
Echelles des vitesses : 𝑜 vz = 𝑈 , 𝑜 vr =?
∂vr ∂vz vr
L’équation de continuité : ∂r + ∂z + r = 0
permet d’écrire : ∂v
∂vz vr 𝑈𝑅
r
o =o =o 𝑜 vr =
∂r ∂z r ℎ
ℎ
Echelles de temps : 𝑜 𝑡 = C’est le temps caractéristique
𝑈 d’une particule fluide

Le nombre de Reynolds caractéristique 𝑈ℎ 𝑈ℎ

de l’écoulement est défini par : 𝑅𝑒 = 𝜌 =
𝜇 𝜈

A. El Khalfi 5
 On suppose que : Re<<1
Analyse de grandeur de la 1ere équations de Navier Stokes
∂vr ∂vr ∂vr 1 ∂p μ ∂2 vr ∂2 vr 1 ∂vr vr
+ vr + vz =− + + 2 + − 2
∂t ∂r ∂z ρ ∂r ρ ∂r 2 ∂z r ∂r r
𝑅𝑈 2 𝑅𝑈 2 𝑅𝑈 2 ?? 𝑈 2 1 𝑅𝑈 2 1 𝑈 2 1 𝑈2 1
ℎ2 ℎ2 ℎ2 𝑅 𝑅𝑒 ℎ2 𝑅𝑒 𝑅 𝑅𝑒 𝑅 𝑅𝑒
1 2 3 4 5 6 7
2
On compare : 𝑜 1 𝑅𝑈 2 ℎ2 𝑅𝑒 𝑜 4 ℎ
= 2 x 2
= 𝑅𝑒 ≪ 1 𝑒𝑡 = ≪1
𝑜(5) ℎ 𝑅𝑈 𝑜(5) 𝑅
On peut donc négliger les termes (1), (2), (3), (4), (6) et (7) devant le
terme (5).
La 1ere équation de Navier-Stokes se réduit à :
∂p ∂2 vr
− +μ 2 =0
∂r ∂z
A. El Khalfi 6
 Analyse de grandeur de la 2nd équation de Navier Stokes
∂vz ∂vz ∂vz 1 ∂p μ ∂2 vz 1 ∂vz ∂2 vz
+ vr + vz =− + + + 2
∂t ∂r ∂z ρ ∂z ρ ∂r 2 r ∂r ∂z

𝑈2 𝑈2 𝑈2 ?? 1 𝑈2 ℎ 1 𝑈2 ℎ 1 𝑈2
ℎ ℎ ℎ 𝑅𝑒 𝑅2 𝑅𝑒 𝑅2 𝑅𝑒 ℎ
1 2 3 4 5 6

On compare :
𝑜 1 𝑈2 ℎ 𝑜 4 1 𝑈2 ℎ ℎ ℎ 2
= 𝑥 2 𝑅𝑒 = 𝑅𝑒 ≪ 1 = 𝑥 2 𝑅𝑒 = ≪1
𝑜(6) ℎ 𝑈 𝑜(6) 𝑅𝑒 𝑅 2 𝑈 𝑅
On peut donc négliger les termes (1), (2), (3), (4) et (5) devant le terme
(6).
La 2nd équation de Navier-Stokes se réduit à :
∂p ∂2 vz
− +μ 2 =0
∂z ∂z
A. El Khalfi 7
 Finalement, on a :
∂p ∂2 vr
− +μ 2 =0
∂r ∂z
∂p ∂2 vr
− +μ 2 =0
∂z ∂z
Comparons les gradients de pression selon r et z, on aura donc :

∂p ∂2 vr UR
o =o μ 2 =μ 3
∂r ∂z h
∂p ∂2 vr U
o =o μ 2 =μ 2
∂z ∂z h
𝜕𝑝 𝑈
𝑜 𝜇 ℎ
𝜕𝑧 ℎ 2
On compare : = = ≪1
𝜕𝑝 𝑈𝑅 𝑅
𝑜 𝜇 3
𝜕𝑟 ℎ
𝜕𝑝 𝜕𝑝
On déduit : 𝑜 ≪𝑜
𝜕𝑧 𝜕𝑟
A. El Khalfi 8
 Dans l’hypothèse où on néglige la variation ∂p ∂2 vr
de pression selon z, les équations de − +μ 2 =0
∂r ∂z
Navier-Stokes se réduisent à : ∂2 vz
2
=0
∂z
Pour la suite du problème, on pose : dp(r, t) = K = Cte
dr
Calcul du champ de vitesse :
∂2 vz
2
= 0 ⇒ vz r, z, t = A z + B
∂z
Conditions aux limites du problème:
vz r, z = 0, t = 0 ⇒ B = 0 z
−U vz r, z, t = −U
vz r, z = h, t = −U ⇒ A = h
h

A. El Khalfi 9
 ∂2 vr dp K 2
μ 2 = = K ⇒ vr r, z, t = z +Cz+D
∂z dr 2μ

Conditions aux limites du problème:

vr r, z = 0, t = 0 ⇒ D = 0
h
vr r, z = h, t = 0 ⇒ C = −K
2μ

K
vr r, z, t = z(z − h)
2μ

A. El Khalfi 10
 Bilan intégral sur le volume de contrôle de la conservation de la masse

Le volume de control W est de frontières S=S1US2US3

A. El Khalfi 11
 Le fluide étant incompressible : 𝑑𝑖𝑣 𝑉 = 0

Aussi, on peut écrire : 𝑑𝑖𝑣 𝑉 𝑑𝑣 = 𝑉 . 𝑛 𝑑Σ = 0

Ω Σ
Donc

𝑉 . 𝑒𝑧 𝑑Σ + 𝑉 . 𝑒𝑟 𝑑Σ + 𝑉. (−𝑒𝑧 ) 𝑑Σ = 0
Σ1 Σ2 Σ3

A. El Khalfi 12
 2𝜋 𝑅
𝑉 . 𝑒𝑧 𝑑Σ = vz 𝑟, 𝑧 = ℎ, 𝑡 𝑟𝑑𝑟𝑑𝜃
Σ1 0 0
2𝜋 𝑅
z
= −U 𝑟𝑑𝑟𝑑𝜃 = −𝑈𝜋𝑅2
0 0 h

𝑉 . −𝑧 𝑑Σ = 0 𝑐𝑎𝑟 𝑙𝑒 𝑑𝑖𝑠𝑞𝑢𝑒 𝑒𝑠𝑡 𝑖𝑚𝑚𝑜𝑏𝑖𝑙𝑒

Σ3
A. El Khalfi 13
 ℎ 2𝜋
𝑉 . 𝑒𝑟 𝑑𝛴 = v𝑟 𝑟 = 𝑅, 𝑧, 𝑡 𝑅𝑑𝜃 𝑑𝑧
𝛴2 0 0
ℎ 2𝜋
𝐾
= 𝑧 𝑧 − ℎ 𝑅𝑑𝜃 𝑑𝑧
0 0 2𝜇
2𝜋 ℎ 3
𝐾𝑅 𝐾𝑅𝜋 ℎ
= 𝑑𝜃 𝑧 2 − ℎ𝑧 𝑑𝑧 = −
2𝜇 0 0 𝜇 6
A. El Khalfi 14
 On a donc : 𝐾𝑅𝜋 ℎ3
− − 𝑈𝜋𝑅2 = 0
𝜇 6

D’où finalement : 6𝜇𝑈𝑅

𝐾=− 3
ℎ

Enfin, on calcule :
𝑑𝑝(𝑟, 𝑡) 6𝜇𝑈𝑅 6𝜇𝑈𝑅 2 𝑟
= 𝐾 = − 3 ⇒ 𝑝 𝑟, 𝑡 = 3
1− + 𝑝0
𝑑𝑟 ℎ ℎ 𝑅
3𝑈𝑅
vr 𝑟, 𝑧, 𝑡 = 3 𝑧(ℎ − 𝑧)
ℎ

z
vz r, z, t = −U
h

A. El Khalfi 15
 Calcul de la force qu’il faut appliquer sur le disque supérieur
Le disque supérieur est soumis à la force de pressions suivantes :

𝐹= 𝑇𝑑Σ = 𝜎. 𝑒𝑧 𝑑Σ
Σ1 Σ1
𝜕𝑣𝑧 𝜕𝑣𝑧 3𝜇𝑈𝑅
𝜎𝑟𝑧 =𝜇 + =
𝜕𝑟 𝜕𝑟 ℎ2
𝜎. 𝑒𝑧 = 𝜎𝜃𝑧 =0
𝜕v𝑧
𝜎𝑧𝑧 = −𝑝 𝑟, 𝑡 + 𝜇
𝜕𝑧

𝐹= 𝜎𝑟𝑧 𝑧=ℎ 𝑑Σ 𝑒𝑟 + 𝜎𝑧𝑧 𝑧=ℎ 𝑑Σ 𝑒𝑧

Σ1 Σ1
Force tangentielle due aux Force normale exercée dans le
effets de frottement sens du mouvement du disque
A. El Khalfi 16
 𝑅 2𝜋
3𝜇𝑈𝑅 3𝜇𝜋𝑈𝑅3
𝜎𝑟𝑧 𝑧=ℎ 𝑑Σ = − 2 𝑟𝑑𝑟 𝑑𝜃 = −
Σ1 0 0 ℎ ℎ2

𝑅 2𝜋
2𝜇𝑈 6𝜇𝑅2 𝑈 6𝜇𝑅𝑈
𝜎𝑧𝑧 𝑧=ℎ 𝑑Σ = − 𝑝0 + + 3
− 3
𝑟 𝑟𝑑𝑟𝑑𝜃
Σ1 0 0 ℎ ℎ ℎ
2
2𝜇𝑈 𝑅
= −𝑝0 𝜋𝑅 2 − 𝜋𝑅2 (1 + 2 )
ℎ ℎ

A. El Khalfi 17