Dynamique des fluides visqueux PSI 2021-2022

TD n◦8 Dynamique des fluides visqueux

ENCPB - Pierre-Gilles de Gennes

Résumé
? Exercice niveau CCP
• Exercice niveau Centrale/Mines
 Exercice nécessitant un sens physique particulier.

1. Ecoulements de type Couette 3. On utilise le viscosimètre pour mesurer la viscosité de l’huile d’olive. Un
opérateur exerce un couple C pour maintenir le cylindre extérieur fixe, le
1.1 Ecoulement sur un plan incliné cylindre intérieur tournant à la vitesse constante ω1 = 3 tours par seconde. On
On reprend l’écoulement de Couette vu en cours mais cette fois, on considère un mesure C = 9, 4.10−2 N·m. On donne R1 = 4, 8 cm, R2 = 5, 0 cm, h = 15 cm.
écoulement libre d’huile sur un plan incliné d’un angle α = 10◦ . L’écoulement est Déterminer la viscosité de l’huile d’olive.
supposé permanent, homogène et incompressible.
1. Faire un schéma et introduire un système de coordonnées adaptés à la des- 1.3 Plaque tectonique (difficile) •
cription du problème. De quel(les) variable(s) dépend le champ des vitesses à On considère une plaque tectonique aussi appelée lithosphère, de longueur L,
l’intérieur du fluide ? d’épaisseur hL , qui se déplace suivant l’axe Ox à une vitesse constante et uniforme
2. En appliquant le principe fondamental de la dynamique à un système bien v0 . Cette plaque repose sur l’asthénosphère qui a une épaisseur hA . Les densités de
choisi, déterminer une équation différentielle vérifiée par le champ des vitesses. la lithosphère et de l’asthénosphère sont ici supposées égales. On prendra l’origine
des ordonnées sous la plaque, à la surface de l’asthénosphère. On considèrera que
3. Quelles sont les conditions aux limites imposées par l’écoulement ? L >> hL et hA et on admettra que la vitesse dans l’asthénosphère vx (z) est parallèle
4. En déduire l’expression du champ des vitesses. à l’axe Ox. Au bas de l’asthénosphère (z = −hA ), la vitesse est nulle. On note p(x, z)
la pression dans l’asthénosphère.
1.2 Viscosimètre de Couette
On considère l’écoulement d’un fluide visqueux de viscosité η entre deux cy-
lindres concentriques de rayons R1 et R2 > R1 , tournant autour de leur axe commun
aux vitesses angulaires respectives ω1 et ω2 . On propose un champ des vitesses du
fluide de la forme  
B
~v = Ar + ~uθ
r
1. Déterminer complètement le champ des vitesses.
2. Par analogie avec le cas de l’écoulement plan, exprimer la force de cisaillement
exercée par le fluide intérieur sur le fluide extérieur à travers une surface
élémentaire d’aire dS normale à ~ur . Déterminer le couple exercé par le fluide
sur le cylindre extérieur.
 TD n◦ 8 Dynamique des fluides visqueux — 2/??

1. En supposant le gradient horizontal de pression nul, exprimer la pression en

fonction de z. On négligera la pression à la surface de la plaque.
2. En effectuant un bilan des forces sur une particule fluide, calculer le profil
des vitesses dans l’asthénosphère. Quelle est la contrainte (force par unité de
surface) tangentielle visqueuse qui s’exerce sous la plaque en mouvement ?
3. Exprimer en fonction de v0 , hL et hA le débit volumique Dv induit par le
déplacement de la lithosphère et de l’asthénosphère par unité d’épaisseur
suivant Oy.
4. Le déplacement global de matière en surface est maintenant compensé par un
flux de retour qui assure un débit volumique globalement nul. Montrer que
cela impose l’existence d’un gradient de pression horizontal ∂∂ Px . Calculer ∂∂ Px
en fonction de hL , hA , v0 et   le profil de vitesse vx (z) résultant. On
η ainsi que
hL 1
pourra ensuite poser A = 6 hA + 2 pour simplifer les expressions.
5. A quelle profondeur la vitesse s’annulle-t-elle dans l’asthénosphère ? Quelle
est la vitesse maximale de retour ? Tracer vx (z)/v0 en fonction de z/hA en On s’intéresse à une portion de liquide contenue dans un cylindre de rayon r
prenant A = 4. (r < R), compris entre les plans de pressions P1 et P2 distants de L.

6. Exprimer la contrainte tangentielle visqueuse qui s’exerce sous la plaque en 1. En coordonnées cylindriques, montrer que la vitesse ne dépend ni de z ni de θ .
mouvement en fonction de η, v0 , hL et hA ? Comment se compare-t-elle à la 2. Que vaut la vitesse du fluide en contact avec la paroi ? Quel phénomène permet
question du 1). d’attribuer cette valeur ? Quel est le signe de dv
dr ?
7. Montrer que l’existence d’un gradient de pression horizontal doit conduire 3. Donner l’expression des forces de viscosité ~Fv s’exerçant sur la portion de
à une inclinaison de la lithosphère (on admettra que l’inclinaison est suffi- liquide contenue dans le cylindre de rayon r. Quelles sont les autres forces
sament faiblepour que la vitesse de l’écoulement dans l’asthénosphère reste présentes s’exerçant sur la portion de liquide contenue dans le cylindre de
horizontale). Exprimer la différence de hauteur entre les deux extrémités de la rayon r ? Donner leur expression.
plaque. 4. Montrer que l’accélération de la portion de liquide est nulle.En déduire que le
champ des vitesse au sein du fluide vérifie l’équation :
2. Ecoulements de type Poiseuille dv
= kr
2.1 Turbulence? dr
À partir de quelle vitesse U l’écoulement de l’eau dans un tube de diamètre où k est une constante qu’on exprimera en fonction des données.
D = 10 cm devient-il turbulent ? Qu’en est-il pour l’air dans les conditions normales 5. Déterminer l’expression complète de v(r) et tracer le profil des vitesses.
de température et de pression ? 6. Déterminer l’expression du débit volumique dans l’ensemble du capillaire en
fonction des données.
2.2 Viscosimètre?
On considère un capillaire constitué d’un cylindre de rayon R, de longueur L, 2.3 Château d’eau?
rempli d’un liquide incompressible de masse volumique µ, de viscosité dynamique On s’intéresse à un réseau de distribution d’eau domestique. Ce type de circuit
η. Le liquide s’écoule longitudinalement dans le capillaire, depuis le haut vers le est alimenté par des châteaux d’eau qui assurent la mise en pression du réseau. L’eau
bas, c’est-à-dire depuis la pression P1 vers la pression P2 . On suppose l’écoulement est supposée incompressible, de masse volumique µ = 1kg.L−1 .
permanent (stationnaire) incompressible et homogène.
 TD n◦ 8 Dynamique des fluides visqueux — 3/??

génération p = 1, les bronches p = 2, et ainsi de suite. On se place en régime

stationnaire et l’air est assimilé à un fluide de viscosité η = 1, 8.10−5 Pl. Il y a
23 générations de voies aériennes dont les 16 premières sont conductrices. Une
bronchiole de génération p est assimilée à un cylindre de rayon r p et de longueur ` p .

1. Quel est l’ordre de grandeur de la pression qui peut être attendue dans une
canalisation au pied d’un château d’eau de 25 m de eau ? On supposera que le
débit dans la canalisation est suffisament faible pour ne pas perturber le champ
de pression dans le château d’eau. On prendra g = 10 m.s−2 .
Dans la conduite, l’écoulement est supposé laminaire et permanent. La pesanteur
est négligée. La conduite est cylindrique, de section S = 1, 0 cm2 et de longueur
L = 100m. La canalisation part du pied du château d’eau et son autre extrémité est On admet que la loi de Hagen-Poiseuille est valable pour p < 16. La figure
à l’air libre. La pression en entrée est notée P1 et celle en sortie P2 . Le fluide est ci-dessous représente les quatre premières générations d’un arbre bronchique. À
newtonien, de viscosité dynamique η = 10−3 P`. chaque génération chaque dimension longueur et rayon est multipliée par h, constante
inférieure à un, identique, pour les deux dimensions.
2. Calculer le débit volumique DV en sortie de la canalisation ainsi que la vitesse 1. Déterminer le nombre N(p) de bronchioles à la pème génération en fonction
moyenne de l’écoulement. de p ainsi que le rayon r p et la longueur ` p de la bronchiole de génération p
3. Rappeler la loi de Poiseuille ainsi que ses conditions d’application. L’exprimer en fonction de p, h, r1 et `1 , valeurs pour p = 1.
en fonction des données de l’énoncé. Donner l’expression de la résistance 2. Exprimer le volume Vp d’une bronchiole de génération p en fonction de V1 , h
hydraulique RH de la canalisation et faire l’application numérique. et p. En déduire le volume total Vpt de génération p. On posera X = 2h3 .
4. Par analogie avec l’effet Joule, exprimer la puissance dissipée par viscosité en Exprimer le volume Vt de l’arbre composé de n générations.
fonction de RH et DV . Sachant que la capacité thermique massique de l’eau 3. Calculer la résistance hydraulique R p d’une bronchiole de génération p en fonc-
vaut ce = 4, 18 J.K−1 .kg−1 , calculer l’élévation de température en 1h si l’eau tion de R1 (résistance hydraulique pour p = 1) et p. En déduire la résistance
récupère toute l’énergie dissipée. hydraulique totale de la génération p puis la résistance hydraulique totale de
l’arbre qui contient n générations.
5. Calculer le nombre de Reynolds pour cet écoulement et conclure.
4. Montrer que le volume total diverge quand n → ∞ pour h supérieur à une
2.4 Arbre bronchique valeur critique hc dont on précisera la valeur numérique. À quelle condition
Dans un arbre bronchique, les voies respiratoires se divisent par dichotomie avec sur h la résistance hydraulique diverge-t-elle ?
une réduction systématique de la longueur et du diamètre. Dans le problème on 5. Pour l’homme, h a été mesuré à h = 0, 85. Estimer la variation de pression
suppose que la trachée se divise en deux bronches. Chacune d’elles se divise à son entre l’entrée et la sortie des 16 premières générations correspondant à une
tour en deux autres, et ainsi de suite. Nous notons ” générations ” les différentes inspiration normale.
subdivisions qui seront indicées par les nombres successifs, p : la trachée est la 6. La modélisation par la loi de Hagen-Poiseuille est-elle pertinente ici ?
 TD n◦ 8 Dynamique des fluides visqueux — 4/??

2.5 Champ des vitesses et de pression dans un écoulement de Poiseuille•

On considère comme dans le cours un écoulement de Poiseuille cylindrique dans
une conduite horizontale de rayon R et de largeur L. La viscosité dynamique du
fluide est notée η et sa masse volumique µ. L’écoulement est supposé homogène,
stationnaire et incompressible.

1. En effectuant un bilan des forces sur une particule de fluide située à la distance
r de l’axe de la conduite, déterminer le champ de pression et le champ des
vitesses dans l’écoulement. La géométrie de la particule de fluide est rappelée
ci-dessous.

1. Question préliminaire : exprimer le volume dV d’une petite cerne cylindrique

(cf. figure ci-dessus) et montrer qu’en sommant le volume de toutes les cernes
comprises dans le cylindre, on retrouve le volume total du cylindre.
2. Exprimer la résultante d~F des forces de pression s’exerçant en amont et en
aval de la cerne.
3. On rappelle l’expression du champ des vitesses d’un écoulement de Poiseuille
cylindrique (on pourra utiliser directement le résultat) :
R2 ∆P r2
 
~v(r) = 1 − 2 ~ur
4η L R
Exprimer la puissance dPdiss des forces de pression s’exerçant sur la cerne en
amont et en aval puis en déduire l’expression de dPdiss .
2.6 Puissance dissipée dans un écoulement de Poiseuille•
On a vu dans le cours qu’un écoulement de Poiseuille vérifie la loi d’Ohm 3. Ecoulement autour d’un objet
hydraulique : 3.1 Mesure de viscosité?
On souhaite mesurer précisément la viscosité de la glycérine, dont on sait qu’elle
8ηL est de l’ordre de 1 Pl. La masse volumique de la glycérine est µg = 1, 260.103 kg·m−3 .
∆P = Pe − Ps = Rth Dv avec Rth = Pour cela, on fait chuter une bille de silice de masse volumique µs = 2, 5.103 kg·m−3
πR4
et de rayon R = 0, 75 ± 0, 01mm dans un récipient cylindrique rempli de glycérine.
c’est à dire que l’extérieur (une pompe par exemple) doit imposer une différence On souhaite mesurer la vitesse en régime permanent pour en déduire η.
de pression ∆P pour assurer un débit Dv . On cherche dans cet exercice à démontrer 1. Déterminer un ordre de grandeur de la vitesse limite vlim suivant que le régime
que la puissance à fournir par la pompe vaut : Pdiss = Rth Dv2 . Pour une conduite est laminaire ou turbulent. Conclure.
horizontale, cette puissance est intégralement dissipée par les forces de viscosité puis 2. Comment faut-il choisir les dimensions du récipient (rayon et hauteur) ?
convertie en chaleur. L’expérience nécessite-t-elle l’utilisation d’une caméra rapide ? D’une caméra ?
Pour calculer Pdiss , on commence par s’intéresser à une petite cerne cylindrique 3. Une fois la vitesse limite atteinte, on mesure une distance parcourue h = 10, 0±
de rayon interne r et d’épaisseur dr (figure ci-dessous). 0, 2 cm en une durée T = 98, 0 ± 0, 2 s. En déduire η avec son incertitude.
 TD n◦ 8 Dynamique des fluides visqueux — 5/??

3.2 Aviron• est atteinte, le pilote actionne la gouverne de profondeur, ce qui provoque la
Un skiff est une embarcation de course avec avirons, dans laquelle prend place rotation de l’avion (le nez s’élève). Identifier cette gouverne sur le schéma
un unique rameur. Il existe aussi des embarcations regroupant jusqu’à 8 rameurs. ci-dessous. À quoi servent les autres parties mobiles ?
Le record du monde de vitesse en skiff est de 17,4 km.h−1 . Avec 8 rameurs, il est
de 21,4 km.h−1 . L’objet de cet exercice est de comprendre pourquoi la différence
n’est pas plus importante, et de dégager des ordres de grandeurs grâce à des facteurs
d’échelle. On donne la viscosité dynamique de l’eau η = 1, 0.10−3 Pl. On suppose
que les embarcations se déplacent en ligne droite avec vitesse de norme v constante.
1. Évaluer le nombre de Reynolds pour un skiff. En déduire la façon dont la force
de traı̂née dépend de la norme v de la vitesse.
2. Si on double la puissance développée par le (ou les) rameur(s), par combien v
est-elle multipliée ?
3. Une personne développant une puissance P0 fait avancer son skiff à 14 km/h.
Quelle est en fonction de P0 la puissance développée par le recordman du
monde ?
4. En supposant, pour simplifier, que la partie immergée d’une embarcation pour 3.4 Portance et traı̂née d’un avion de tourisme léger?
8 rameurs est homothétique de celle d’une embarcation pour 1 seul, et en Les courbes ci-dessous, obtenues par des essais en souffleries à l’aide d’une
faisant les hypothèses qui vous semblent légitimes, établir un lien entre la balance aérodynamique, représentent l’évolution des coefficients de portance Cz et
vitesse d’un skiff et celle d’une embarcation à 8 rameurs. Commenter, compte de traı̂née Cx en fonction de l’angle d’incidence α pour un profil d’aile.
tenu des records de vitesse pour les deux types de bateaux.

3.3 A 380?

1. Calculer la vitesse d’un A380 au décollage au niveau de la mer, à une température

de 20 ◦ C, pour une masse de 421 tonnes, une surface portante de 845 m2 et un
coefficient de portance Cz = 1,38. Calculer la variation relative de cette vitesse
due à une variation de température de 20 ◦ C. Commenter.
2. Vaut-il mieux décoller avec le vent de face, de derrière ou de côté ?
3. On distingue 3 axes de rotation pour l’avion : l’axe de tangage, dans la di- 1. (a) Commenter l’allure des courbes ainsi que les ordres de grandeur des
rection des ailes, l’axe de roulis dans la direction de l’avion et l’axe de lacet valeurs respectives de Cz et Cx pour une incidence donnée. Discuter
perpendiculaire aux deux autres. Représenter ces 3 axes sur un schéma. l’intérêt de ce profil d’aile.
4. Différentes parties mobiles (cf. schéma ci dessous) permettent de jouer sur (b) À partir de quelle incidence l’aile est-elle potentiellement en situation
les rotations autour de ces axes. Par exemple, lorsque la vitesse de décollage de décrochage ? Proposez une interprétation en vous appuyant sur vos
 TD n◦ 8 Dynamique des fluides visqueux — 6/??

connaissance sur le couche limite. Dans la comparaison des courbes de

Cz et de Cx était-il possible de prévoir l’augmentation de traı̂née lorsque
la portance chute à forte incidence ?
2. Dans toute la suite, on assimilera pour simplifier les coefficients de portance et
de traı̂née de l’aile à ceux de l’avion et on supposera que le pilote ne fait pas
usage des volets pour décoller. La corde fait avec l’axe longitudinal de l’avion
un angle β de 4◦ (angle de calage). La force de traction ~T qu’exerce le moteur
de l’avion sera prise de même support que son axe longitudinal.

(a) Dans la phase de roulage précédent l’instant du décollage, le pilote tire

momentanément sur le manche de sorte que la queue de l’avion s’enfonce.
L’avion étant à charge maximale, à quel angle δ le pilote doit-il choisir de
placer l’axe longitudinal de l’avion par rapport au sol ? Faire un schéma
en représentant les forces s’exerçant sur l’avion et le vent relatif à l’infini.
Quelle est la vitesse minimum à l’instant où les roues quittent le sol ? On
donne la surface alaire S = 12 m2 et la masse maximale m = 750 kg.

(b) Juste après le décollage, le pilote remet le manche au neutre et fait ainsi
monter l’avion sur une trajectoire confondue avec le support de son axe
longitudinal incliné du même angle δ par rapport au sol qu’à l’instant du
décollage. Comment est modifié le schéma précédent ? Sachant que la
puissance maximale du moteur équipant l’appareil est de 75 kW, à quelle
pourcentage de la puissance maximale faut-il alors se placer pour assurer
cette phase ascensionnelle stabilisée ?