LP n°8 : Notion de viscosité d'un uide. Ecoulements visqueux.

(LP n°3 du Concours

Externe Spécial)
Ugo Hincelin

5 Décembre 2016

Niveau 2.1.1 Origine . . . . . . . . . . . . . . . . . . . . . . . . . 3

2.1.2 Fluide newtonien . . . . . . . . . . . . . . . . . . . 3
Objectifs 2.1.3 Fluide non-newtonien . . . . . . . . . . . . . . . . . 3
2.2 Equation de Navier-Stokes . . . . . . . . . . . . . . . . . . 3
 Comprendre la notion de viscosité et mesurer sa valeur 2.2.1 Terme advectif . . . . . . . . . . . . . . . . . . . . 4
 Etablir une équation de la dynamique des uides qui relie forces et
2.2.2 Terme diusif . . . . . . . . . . . . . . . . . . . . . 4
champ de vitesse
 Dénir le nombre de Reynolds et ses implications 2.2.3 Remarque sur la diérence entre l'advection et la
convection . . . . . . . . . . . . . . . . . . . . . . . 4

Prérequis 3 Ecoulement 4
3.1 Nombre de Reynolds . . . . . . . . . . . . . . . . . . . . . 4
 Statique des uides (équation locale, poussée d'Archimède) 3.2 Chute d'une boule dans un uide (expérience quantitative) 4
 Cinématique des uides (écoulement stationnaire et incompressible) 3.3 Ecoulement de Couette . . . . . . . . . . . . . . . . . . . . 5
3.3.1 Champ de vitesse . . . . . . . . . . . . . . . . . . . 5
Table des matières 3.3.2 Dissipation d'énergie . . . . . . . . . . . . . . . . . 5
3.4 Ecoulement de Poiseuille . . . . . . . . . . . . . . . . . . . 6
1 Introduction 2
4 Conclusion 6
2 Viscosité 2
2.1 Nécessité d'une contrainte tangentielle . . . . . . . . . . . 2 5 Questions et commentaires 7
1
1 Introduction
La viscosité est une caractéristique d'un milieu liée à ses propriétés
d'écoulements et aux interactions des éléments du milieu lorsque celui-
ci est mis en mouvement [2]. La viscosité inue la façon dont les uides
s'écoulent. Cela se voit facilement pour des objets de la vie quotidienne.
Des liquides comme le miel et l'eau s'écoulent diéremment car ils ont des
viscosités diérentes (Figure 1). Lorsque l'eau s'écoule d'un robinet, la
vitesse d'écoulement inuence l'apparition d'une instabilité plus ou moins
rapidement (Figure 2). La viscosité d'un uide inuence le mouvement de Figure 2  Régime d'écoulement de l'eau en sortie d'un robinet
particules au sein du uide : des particules en suspension vont sédimenter
plus ou moins vite selon le uide dans lequel elles se trouvent (Figure 3).

Figure 3  Sédimentation
Figure 1  Ecoulement d'un miel et de l'eau
le couteau met une couche de miel en mouvement, que cette couche met
en mouvement la couche située sous elle, etc. Les forces de pression à
2 Viscosité −−→
l'interface entre deux couches dFn = ±p(M, t) dS → −
uy , forces exercées par
[1] p.299 une couche sur une autre, sont verticales donc incapables de transmettre
un mouvement horizontal. Il existe donc des forces de contact tangentielles
telles que :
2.1 Nécessité d'une contrainte tangentielle −→
dFt = σ(M, t) dS →
−
ux (1)
On étale du miel sur une tartine. La tartine est supposée plane et xe,
et la lame du couteau est plane et parallèle à la tartine. Le couteau a
→
−
une vitesse U = U →−
ux selon l'axe (Ox) (voir Figure 6). On constate que où σ(M, t) est la contrainte tangentielle au point M à l'instant t.

2
2.1.1 Origine 2.1.3 Fluide non-newtonien
A l'échelle microscopique, le transfert de la quantité de mouvement Pour un tel uide, la viscosité dépend du gradient de vitesse. Exemple :
entre les couches se fait par l'agitation moléculaire et les interactions entre le yaourt. Sa viscosité diminue quand on le bat rapidement. Le sang est
les molécules de type  liaison H  et  van der Waals . Les molécules aussi un exemple.
les plus rapides augmentent la vitesse des molécules les plus lentes et
réciproquement.
2.2 Equation de Navier-Stokes
2.1.2 Fluide newtonien On suppose le uide newtonien et incompressible.
Détails du calcul donnés dans [1] p.302 (remarque : utiliser le Laplacien
Pour un uide newtonien, la contrainte tangentielle est proportionnelle tel que 4 et non →−
4 ).
à l'inhomogénéité du champ des vitesses. Dans le cas de notre exemple, la −→ −−→ −−−−→
dF = dFy + dFy+dy = −η ∂v∂y x (y)
dS →
−
ux + η ∂vx (y+dy) dS →
−
ux
contrainte tangentielle s'écrit : ∂y
Avec Taylor ordre 1 : f (a + h) = f (a) + f (a)h, on obtient :
0
−→
dS dy →
− dx dy dz →
−
2 v (y) 2 v (y)
dF = η ∂ ∂yx
2 ux = η ∂ ∂yx
2 ux , ce qui donne (on admet la
∂vx (M, t)
σ(M, t) = η (2) généralisation à 3D) l'expression de la force volumique de viscosité pour
∂y
un uide newtonien incompressible :
où η est la viscosité dynamique. [η] = M L−1 T −1 ; unité : Poiseuille 1 −→
dF
1PI = 1Pa.s = η 4→
−
v (3)
La viscosité dépend de la température et de la pression. Pour un liquide, dτ
elle diminue avec la température car on  casse  les liaisons. Pour un √
gaz,
On écrit le principe fondamental de la dynamique pour une particule
elle augmente avec la température (théorie cinétique des gaz, η ∝ T ,
de uide de masse dm = ρ dτ :
[5] p.69) Le tableau suivant donne des valeurs dans les conditions usuelles →
− −−→ −→ −→
(CATP - conditions ambiantes température pression - 25°C et 1 bar) : ρ dτ ddtv = −gradp dτ + η 4→ −
v dτ + dF 0 où dF 0 désigne la résultante des
−−−→ −−→
Corps Air Eau Huile d'olive Glycérine Miel autres forces appliquées à la particule (comme le poids). dFp = − gradp dτ
Viscosité (PI) 1, 8.10−5 1, 0.10−3 0,1 1,4 ∼ 10 est la résultante de forces de pression. −
→ →
−
Avec l'expression de la dérivée particulaire ([3] p.4) dv = ∂∂tv dt +
1. Jean-Léonard-Marie Poiseuille1 (1797 - 1869) est un physicien et médecin fran- →
−
∂ v
→
− →
− −
→
dx + ∂∂yv dy + ∂∂zv dz ⇔ dv
→
−
= ∂∂tv + (→
− −−→ −
v .grad)→
v , on obtient l'équation
çais. On lui doit diérents mémoires sur le c÷ur et la circulation du sang dans les ∂x
de Navier-Stokes :
dt

vaisseaux (l'hémodynamique), qui lui permirent d'établir en 1844  à travers son ou-
vrage Le Mouvement des liquides dans les tubes de petits diamètres  les lois de
l'écoulement laminaire des uides visqueux dans les tuyaux cylindriques. Source : wi- d→
− −−→0
 →
∂−

v v −−→ → −−→ dF
kipédia. ρ
dt
=ρ
∂t
→
− − →
−
+ ( v .grad) v = −grad p + η 4 v +
dτ
(4)

3
2.2.1 Terme advectif  Re  1 : transport advectif plus important. Ecoulement turbulent.
−−→  →
− −−→ 
→
−  Re  1 : transport diusif dominant. Ecoulement laminaire.
ρ (→−
v .grad)→
−
v est le terme d'advection. Si ρ ∂∂tv + (→
−
v .grad)→
−v = 0
Voir la Figure 4.
→
−
dv →
−
alors dt
= 0 . Conservation de la quantité de mouvement lors de l'écou-
lement, i.e. qu'il y a transport de la quantité de mouvement.

2.2.2 Terme diusif

→
−
η 4→−v est le terme de diusion (due à la viscosité). Si ρ ∂∂tv = η 4→
−
v
→
− →
−
alors ∂ v − 4→
∂t ρ
η −v = 0 . C'est une équation de diusion. ν = est appelée
ρ
η

 viscosité cinématique .

2.2.3 Remarque sur la diérence entre l'advection et la convec-

tion
L'advection est le transport d'une quantité dû à son entraînement
par le mouvement d'un milieu, quelque soit l'origine de ce mouvement
(gravitaire ou entraînement forcé). Exemple : goutte d'encre dans un uide
en écoulement (il y a advection + diusion).
La convection est une advection causée par une instabilité locale du
milieu. Exemples : eau chaude sous eau froide dans une casserole, ou masse
d'air chaud dans l'atmosphère.
Source : [2]
Figure 4  Ecoulements laminaire et turbulent (extrait de [4])
3 Ecoulement
3.2 Chute d'une boule dans un uide (expérience
3.1 Nombre de Reynolds quantitative)
Le nombre de Reynolds évalue l'importance relative du terme advectif On lme la chute d'une bille d'acier dans la glycérine (Figure 5). Acqui-
sur le terme diusif : sition vidéo et traitement avec les logiciels Atelier Scientique et Regressi.
−−→ − On observe un régime transitoire et un régime permanent sur le graphe
ρ (→
advection −
v .grad)→ v
(5) vitesse = fonction(temps). La vitesse limite nous donne la valeur de la
v ρv L ρvL
Re = = ∼ =
diusion η 4→−
v η Lv2 η viscosité dynamique de la glycérine.

4
 3.3.1 Champ de vitesse
→
− →
−
Ecoulement stationnaire ( ∂∂tv = 0 ) et incompressible (div → −v = ∇→
−v =
z 0, ici cela donne ∂x = 0) (voir Figure 6).
∂vx
Bille d'acier
Conditions aux limites (interface uide solide) : uide réel en admettant
que les composantes tangentielles des vitesses du uide et du solide sont
égales (condition de non glissement). Remarque : pour un uide parfait,
Glycérine
la viscosité est nulle.
On obtient le prol de vitesse suivant :
Uz
Eprouvette
vx (z) = (7)
e
Figure 5  Chute d'une bille d'acier dans la glycérine
→
−
f = −6πη R →
−
v est la force de traînée à faible Reynolds, solution ana-
lytique par Stokes de l'équation de Navier-Stokes.
→
−
P = m→ −g = 43 πR3 ρa →
−
g est le poids de la bille
−
→ 3 → −
FA = − 3 πR ρg g est la poussée d'Archimède
4

On écrit le PFD et on obtient

t
v = C(1 − exp(− )) (6)
τ

avec C = (ρa6πηR et τ = 2R9ηρa

−ρg )gV 2

L'expérience donne vlimite = (26, 0 ± 0, 9)mm.s−1 et τ = (81 ± 36)ms.

2
La viscosité est alors η = 2R (ρval−ρg )g = (1, 25 ± 0, 04)Pa.s (ou PI).
Remarque : justication de la hauteur de l'éprouvette telle que la vitesse Figure 6  Ecoulement de Couette. Prol de vitesse en vert.
limite est atteinte bien avant de toucher le fond (exercice 4 p.345 de [1]).

3.3.2 Dissipation d'énergie

3.3 Ecoulement de Couette
Voir [1] p.305.
voir p.304-306 de [1] On calcule la puissance reçue par la particule uide via les forces de
Le miel sur la tartine est un écoulement de Couette. viscosité situées en haut et en bas de la particule. On obtient la puissance

5
volumique dissipée dans l'écoulement :
dP η U2
= 2 (8)
dτ e
En eet, la vitesse de la particule est constante (écoulement stationnaire)
donc son énergie cinétique est indépendante du temps. La puissance reçue
est donc dissipée au sein de la particule, son énergie interne augmente, sa
température augmente.
Un écoulement de uide visqueux s'accompagne d'une dissipation
d'énergie au sein du uide.
Figure 7  Ecoulement de Poiseuille plan. Prol de vitesse en vert.
3.4 Ecoulement de Poiseuille
couche limite où l'écoulement est visqueux, et au delà de cette couche
(Suggestion d'ajout de cette partie par le correcteur) limite, l'écoulement peut être décrit comme parfait.
Traitée dans [5] en p.162 (très succinct), mieux décrit dans [1] en p.334 Notons quelques applications relatives à la viscosité et aux écoulements
(bien aussi en p.558 dans le Pérez Méca). visqueux :
Ecoulement stationnaire, laminaire et incompressible dans une conduite
cylindrique d'axe (Oz), de rayon R et de longueur L (voir Figure 7). La
pesanteur est négligée. On applique une diérence de pression selon z ce Application : la vitesse d'un cycliste La viscosité est source de dis-
qui met en mouvement le uide. En coordonnées cylindrique on obtient le sipation d'énergie au sein d'un uide. Un cycliste, qui cherche à atteindre
prol de vitesse suivant : une vitesse maximale, va donc développer une puissance qui va en partie
permettre de lutter contre les forces de traînée exercées par l'air. Cette
2 2
(p(0) − p(L))(R − r ) énergie fournie sera dissipée dans l'air, d'autant plus que les forces de
vz (r) = (9) traînée seront importantes. Notons l'intérêt qu'a le cycliste à utiliser des
4ηL
vêtements de texture optimisée pour réduire les actions de contact entre
l'air et ses vêtements, à pratiquer en altitude pour diminuer la masse vo-
Application Important pour les écoulements dans les canalisations. lumique de l'air, ou encore à prendre une position aérodynamique pour
diminuer la taille de la couche limite en diminuant les zones où le gradient
de vitesse est grand et ainsi diminuer les forces de viscosité.
4 Conclusion
Couche limite Pour faire le lien entre les uides parfaits et mon miel qui Application : le temps de séjour dans une installation d'épura-
accroche, parler de la couche limite (section 1.6 p.309 de [1]). A l'échelle tion L'expérience de la bille dans la glycérine peut être assimilée au
moléculaire, une couche de liquide se xe sur le solide. Il y a alors une processus de sédimentation dans les stations d'épuration. An de traiter

6
plus d'eau, on cherche à diminuer le temps de sédimentation, et donc à Références
augmenter la vitesse limite de chute des particules en suspension dans
l'eau. Pour cela, on traite l'eau an d'agréger les particules an d'aug- [1] S. Olivier, K. Lewis, 2014, Tec & Doc Lavoisier, Compétences prépas
menter leur rayon. Physique PC PC*
[2] R. Taillet, L. Villain, P. Febvre, 2008, Dictionnaire de physique
[3] M. Rieutord, 2014, Une introduction à la dynamique des uides
5 Questions et commentaires [4] E. Guyon, J.-P. Hulin, L. Petit, 2005, Ce que disent les uides, Belin
Pour la science
 Calculer le Reynolds de l'eau du robinet. Reynolds important (∼ [5] E. Guyon, J.-P. Hulin, L. Petit, 2012, Hydrodynamique physique, EDP
106 ), donc on peut avoir un régime laminaire à haut Reynolds. En Sciences, CNRS Editions
revanche, au bout d'un moment, on va voir se développer une tur- [6] [1] essentiellement et [5] pour compléter. Les autres ponctuellement.
bulence, d'autant plus rapidement que la vitesse d'écoulement de
l'eau est grande. Les écoulements aux échelles macro sont toujours à
grand Reynolds. La question est : quelle distance va parcourir l'eau
en écoulement stable ? Petit Re : capillaire (écoulement sanguin)
 Pourquoi des bulles se forment à faible débit du robinet ? Surface
optimale permettant de minimiser l'énergie de tension supercielle.
Voir section 4.1 p.326 de [1].
 Exemple de uide non newtonien ? Peinture. Existence de particules
plus grandes que l'échelle atomique mais assez petite à l'échelle du
uide. Voir section 4.4 p.140 dans [5]
 Qu'est-ce qu'un uide parfait ? v diérent de 0 aux limites.
 Et si le uide est compressible ? apparition d'un deuxième coecient
de viscosité
 Attention à l'opérateur Laplacien : 4→ v = ( ∂∂xv2x + ∂∂yv2x + ∂∂zv2x ) →
− −
2 2 2
ux +
∂ 2 vy →
)−
2 ∂ 2 vy
( ∂∂xv2y + ∂y 2
+ ∂z 2
uy + ...
 Est-ce que la formule de Stokes (expression de la force de traînée
de la bille dans la glycérine) est valable durant le régime transitoire
(car c'est un régime non stationnaire) ? Il faut comparer l'échelle de
temps hydro avec l'échelle de temps méca. Si τhydro  τméca alors
OK. τméca est le temps de transition mesuré τ = 92 R ηρa . τhydro = ρ ∂v
2

∂t
décroît rapidement par rapport à τméca .