Chute d’une bille dans un fluide très visqueux

On s’intéresse à une bille d’acier sphérique de rayon 𝑅, de masse 𝑚, de masse volumique 𝜌acier=7850 kg m−3, de vitesse 𝑣⃗=−𝑣𝑒⃗𝑧 lâchée sans vitesse initiale dans une éprouvette remplie de glycérine (de masse volumique 𝜌glycérine=1260 kg m−3 et de viscosité de l’ordre de 𝜂∼0⋅100 Pa ⁠s).

Le champ de pesanteur 𝑔⃗ est uniforme. On pose 𝑔′=𝑔(1−𝜌glycérine𝜌acier). L’axe (𝑂𝑧) est vertical ascendant.

L’éprouvette est de diamètre très supérieur à celui de la bille. La force de frottements visqueux exercée par le fluide sur la bille est 𝐹⃗=−6𝜋𝜂𝑅𝑣⃗.

1/ Déterminer l’équation différentielle vérifiée par la vitesse 𝑣 de la bille.

Coup de pouce 1
Quelles sont les 3 forces qui s’exercent sur la bille ?
Coup de pouce 2
Appliquer le théorème de la résultante cinétique à la bille.
Corrigé

Le théorème de la résultante cinétique appliqué à la bille s’écrit

𝑚d𝑣⃗d𝑡=𝑚𝑔⃗+𝜌glycérine𝑉𝑔𝑒⃗𝑧−6𝜋𝜂𝑅𝑣⃗

Comme 𝑣⃗=−𝑣𝑒⃗𝑧, la projection sur 𝑒⃗𝑧 donne

−𝑚d𝑣d𝑡=−𝑚𝑔+𝜌glycérine𝑉𝑔+6𝜋𝜂𝑅𝑣=−𝑚𝑔+𝜌glycérine𝑚𝜌acier𝑔+6𝜋𝜂𝑅𝑣=−𝑚𝑔′+6𝜋𝜂𝑅𝑣

Dans sa forme canonique, cette équation s’écrit

d𝑣d𝑡+6𝜋𝜂𝑅𝑚𝑣=𝑔′

2/ Calculer la durée caractéristique 𝜏 associée à l’équation différentielle pour 𝑅=1,5 mm.

Corrigé

La durée caractéristique de l’équation différentielle est

𝜏=𝑚6𝜋𝜂𝑅=43𝜋𝑅3𝜌acier6𝜋𝜂𝑅=2𝑅2𝜌acier9𝜂=3⋅10−3 s

3/ Exprimer la vitesse limite en fonction de 𝑅, 𝑔′, 𝜌acier et 𝜂.

Coup de pouce 1
Quelle est la solution particulière de l’équation différentielle ?
Coup de pouce 2
Exprimer la masse de la bille en fonction de 𝜌 et de 𝑅.
Corrigé
𝑣lim=𝑔′𝑚6𝜋𝜂𝑅=𝑔′43𝜋𝑅3𝜌acier6𝜋𝜂𝑅=2𝑔′𝑅2𝜌acier9𝜂

4/ On mesure 𝑣1 s, la norme de la vitesse une seconde après avoir lâché la bille pour différentes tailles de bille. La courbe ci-dessous montre l’évolution de 𝑣1 s, en fonction de 𝑅2. En déduire la valeur de la viscosité 𝜂 de la glycérine.

Coup de pouce 1
En comparant 1⁠ ⁠s et 𝜏, dans quel régime se trouve-t-on ? Quelle est l’expression de la vitesse dans ce régime ?
Corrigé

1 s≫𝜏, la vitesse mesurée au bout de 1⁠ ⁠s est donc la vitesse limite 𝑣lim=2𝑔′𝑅2𝜌acier9𝜂.

Le point le plus à droite correspond à une bille de 1,5⁠ ⁠mm. Pour les autres, le temps caractéristique est plus court encore que celui précédemment calculé.

Les points expérimentaux sont alignés. La pente correspond à 2𝑔′𝜌acier9𝜂.

Une lecture graphique donne une pente de 𝑝=0,029 m s−12,4 mm2=1,2⋅104 s−1 ⁠m−1.

On a donc

𝜂=2𝑔′𝜌acier9𝑝=1,2 Pl

5/ Calculer le nombre de Reynolds pour la plus grosse sphère. Est-il légitime de considérer des frottements fluides linéaires ?

Coup de pouce 1
Le nombre de Reynolds croît-il ou décroit-il avec 𝑅 ?
Corrigé

La plus grosse bille (𝑅=1,5 mm) atteint 𝑣lim=2,6⋅10−2 m s−1, d’où

𝑅𝑒=𝑣𝐷𝜈=2𝑣𝑅𝜌glycérine𝜂=7,7⋅10−2

𝑅𝑒 est d’autant plus grand que 𝑅 est grand. Il est très inférieur à 1 pour la plus grande des billes, donc pour toutes les autres aussi. Il est donc légitime de considérer un modèle linéaire des frottements fluides.