Voir la source de Forme d'une goutte de pluie

[[File:Raindrops sizes.svg|thumb|upright=1.2|Dynamiques de la formation des gouttes de pluies.]]
L'article [[Goutte_(physique)|goutte]] montre que si la plupart du temps, les gouttes d'eau dans l'air, par exemple, sont sphériques (brouillard, bruine et petites pluies), elles peuvent également prendre des formes très diverses :
sphériques, disque incurvé comme un globule rouge, 
parachute avec bourrelet.

== Analyse physique ==

La forme de la goutte est déterminée par un équilibre entre la [[tension superficielle]] qui tend à conserver une forme sphérique et la résistance de l'air (ou traînée aérodynamique). Cette dernière dépend de la vitesse de chute de la goutte. Prédire la forme de la goutte consiste à déterminer sa vitesse et la  traînée aérodynamique puis à en déduire sa déformation.

La goutte, initialement sphérique de rayon <math>R</math>,
chute sous l'effet de son poids <math>\rho_{\rm liq}\,\frac{4}{3} \pi R^3g</math>
et est ralentie par la traînée aérodynamique.

Sous l'effet de la  traînée aérodynamique, la goutte se déforme.
Une variation de courbure <math>C</math> 
d'une extrémité à l'autre de la goutte,
du fait de la [[tension superficielle]]
correspond à un écart de pression 
<math>\Delta p=\gamma\,C</math> (cf [[Pression de Laplace]]).

La [[Grandeur sans dimension]] pertinente dans ce problème est le [[nombre de Weber]] qui compare les effets inertiels aux effets de tension de surface.

== Accélération et vitesse terminale ==

Dans le référentiel de l'air,
l'accélération de la goutte 
résulte d'une part de son poids (la poussée d'Archimède due à l'air est ici négligée),
d'autre part de la traînée aérodynamique.
[[File:CX_SPHERE.png|thumb|Valeur du <math>C_x</math> de la sphère à tous les Reynolds possibles.]]

Pour cette traînée aérodynamique, on peut utiliser la formulation inertielle, propre aux hauts [[nombre de Reynolds]]<ref>Cette formulation n'a pas de signification physique au bas Reynolds, mais elle est numériquement exacte ; on peut donc en faire usage.</ref>. Cette formulation inertielle utilise le <math>C_x</math> bien connu dont la valeur dans toute la plage des Reynolds possibles dessine la courbe rouge du graphe ci-contre.

On peut donc écrire :

<math>\frac{4}{3} \pi R^3 \rho_{\rm liq}\,\frac{{\rm d}V}{{\rm d}t}
= \frac{4}{3} \pi R^3\rho_{\rm liq}\,g
-\frac{1}{2}\rho_{\rm air}\,V^2\, \pi R^2 C_x</math>

Autrement dit :

<math>\rho_{\rm liq}\,R^3\,\left(
\frac{4}{3} \pi g - \frac{4}{3} \pi\frac{{\rm d}V}{{\rm d}t}\right)
= \frac {1}{2}\rho_{\rm air}\,V^2\, \pi R^2 C_x </math>


Le <math>C_x</math> présent dans cette équation peut être pris, '''en ordre de grandeur''', comme unitaire (voir la courbe rouge ci-contre).<br>

Après simplifications, on gagne à présenter la précédente égalité sous la forme suivante :

:<math> \frac { {\rm d}V} {{\rm d}t} +

\bigg[ \frac { \rho_{\rm air}} {\rho_{\rm liq}} 

 \ \frac {3}{8R} \ C_x \bigg] \ V^2 = g </math>
En effet, on peut y reconnaître l'équation différentielle classique :

:<math>V' + a\ V^2 = g</math>

:la variable de <math>V</math> et <math>V'</math>étant bien sûr le temps <math>t</math>, <math>a</math> représentant le contenu des grands crochets et <math>g</math> toujours l'accélération de la pesanteur.<br>

Pour intégrer cette équation différentielle, il faut s'assurer que le contenu <math>a</math> des grands crochets est une constante. Ce sera le cas si le <math>C_x</math> est pris comme constant<ref>En réalité le <math>C_x</math> de la sphère est très variable avec sa vitesse de chute, à diamètre et fluide donnés. Cependant, les aérodynamiciens considèrent qu'en dessous du Reynolds diamétral {{formatnum:300000}}, le <math>C_x</math> peut être pris comme constant et valant <math>0,5</math></ref>. Au demeurant, il est facile de comprendre que, quel que soit le <math>C_x</math> (variable) éprouvé par la goutte lors de son accélération dans l'air, sa vitesse terminale finira par être fonction de son <math>C_x</math> terminal (c.-à-d. le <math>C_x</math> correspondant à cette vitesse terminale, ou plus exactement au Reynolds à cette vitesse terminale). Il n'y a donc pas d'objection à l'intégration.<br>

La solution de l'intégration de l'équation différentielle ci-dessus est, classiquement :
: <math>V =\frac { \sqrt g \; Tanh (a \ \sqrt g \ t )} {\sqrt a}</math>

:du moins si l'on fait usage de la condition initiale <math>V(t=0) = 0 </math>

Cette intrusion de la fonction [[tangente hyperbolique]] est habituelle dans les calculs de balistique.<br>

Par chance, pour le temps très grand auquel la goutte va connaître sa vitesse de chute stabilisée <math>V_\infty</math>, la fonction <math>Tanh</math> tend vers l'unité. On dégage donc :

:<math>V_\infty = \sqrt \frac {  g } { a}</math>

:soit :

:<math>V_\infty = \sqrt \frac {  g } {\bigg[ \frac { \rho_{\rm air}} {\rho_{\rm liq}} 

 \ \frac {3}{4D} \ C_x \bigg]}</math> <br>

:du moins si l'on utilise préférentiellement le diamètre D de la sphère. Ce qui peut s'écrire encore :

:<math>V_\infty = \sqrt{ \frac{\rho_{\rm liq}} {\rho_{\rm air}} \frac {  4 \ g \ D } {  3 \ C_x} }</math>

=== Observation de ce résultat ===
[[File:Vitesse et Reynolds gouttes d'eau selon diamètre, Beard & Pruppacher.png|thumb|left|Vitesse et Reynolds gouttes d'eau dans l'air selon leur diamètre, Beard & Pruppacher]]
[[File:Vitesse terminale des gouttes d'eau de diamètre inférieur à 1mm.png|thumb|Vitesse des gouttes d'eau de diamètre < {{unité|1|mm}}.]]
[[File:Anim déformation goutte pluie.gif|thumb|left|Animation de la déformation des gouttes de pluie selon leur diamètre]]

* Dans leur ouvrage essentiel, Clift, Grace et Weber<ref>P. 113, Bubbles, Drops, and Particles, Clift, Grace and Weber, Academic Press, 1978,[https://books.google.com/books?hl=fr&lr=&id=UUrOmD8niUQC&oi=fnd&pg=PP1&ots=pHxRYFdPbU&sig=KLSSMVQ8kqFmA_JnP8wrjZSkjuM]</ref> proposent le même libellé pour la vitesse terminale des sphères rigides (si l'on néglige également la Poussée d'Archimède de l'air sur la goutte).
:Les mêmes auteurs notent qu'entre les Reynolds 750 et {{formatnum:350000}} le <math>Cx</math> de sphères rigides se montre à peu près constant et égal à <math> 0,445</math> mais cette plage de Reynolds est justement celle où les gouttes perdent progressivement leur forme sphérique (ce que l'on observe ci-contre à gauche sur le graphe  du <math>C_x</math> de la sphère en fonction du diamètre au fait que la courbe du <math>C_x</math> diverge de la courbe fuchsia réservée aux sphères lisses).<br>
Voir aussi à gauche l'animation montrant la déformation des gouttes de pluie (les cercles rouges représentent les ''sphères équivalentes'' aux gouttes).


*À l'observation de ce résultat, on peut être tenté de penser que la vitesse terminale <math>V_\infty </math> est fonction de la racine carrée du diamètre, alors que les graphes de Beard et Pruppacher (ci-contre à droite) indiquent une vitesse terminale évoluant linéairement avec le diamètre (au moins pour les diamètres inférieurs à 1 mm). La raison de ce paradoxe apparent est que le <math>C_x</math> n'est pas constant mais, au contraire, variable avec le diamètre et le Reynolds<ref>Cette variabilité aurait pu interdire l'intégration en Tangente hyperbolique, sauf à fixer, par itération, le <math>C_x</math> à sa valeur maximale (que la goutte finira toujours pas atteindre) pour l'intégration.</ref>, comme l'indique le graphe ci-dessus à gauche ainsi que le grand graphe du Cx de la sphère selon le Reynolds.

*Tous les calculs ci-dessus ont été faits en utilisant comme diamètre des gouttes '''le diamètre de la sphère équivalente''', c.-à-d. le diamètre de la sphère de même volume que la goutte. C'est une habitude chez les physiciens-météorologues puisqu'ils sont confrontés à des gouttes dont la forme n'est pas forcément sphérique (au moins pour les plus grosses)(c'est d'ailleurs ce qui explique que sur le graphe du Cx de la sphère aux différents Reynolds, la courbe "brouillard et pluie" diverge de la courbe rouge des sphères rigides). Il faut d'ailleurs ajouter que le <math>C_x</math> utilisé dans ces mêmes calculs et dans les différents graphes est établi en référence à la surface frontale de la sphère de même volume que la goutte ; de même, le Reynolds est construit en référence au diamètre de cette sphère de même volume.

*Avec ce libellé de la vitesse de chute stabilisée (ou vitesse terminale), on peut calculer la Traînée <math> \textstyle \frac{1}{2}  \rho_{\rm air}  V_\infty^2 \ (\pi R^2) \ C_x </math> : on trouve alors comme libellé de cette traînée le poids de la goutte.
[[File:Eclatement de goutte de pluie.gif|thumb|Éclatement d'une grosse goutte de pluie.]]
*Quand le rayon <math>R</math> est fort, la vitesse <math>V_\infty</math> est aussi plus forte. Inversement, cela explique que dans la nature, les très petites gouttes de liquide (comme celles créant le brouillard) chutent avec une vitesse très faible alors que les grosses gouttes de pluie chutent avec une forte vitesse : Lorsque les gouttes de pluie atteignent un diamètre de 5,5 à 6 mm, leur vitesse est si forte que les efforts aérodynamiques les font éclater ; ceci explique qu'il n'existe pas de goutte de pluie de diamètre supérieur à 5,5 ou 6 mm (animation ci-contre).

*Pour une masse volumique <math>\rho_{\rm air}</math> plus forte, la vitesse terminale <math>V_\infty</math> est plus faible. Inversement, dans le vide (pour <math>\rho_{\rm air} = 0</math>) , une goutte acquerrait une vitesse quasi infinie (du moins si la goutte ne s'évaporait pas immédiatement).
*Quand le <math>C_x</math> de la goutte est plus fort (goutte déformée, par exemple), la <math>V_\infty</math> est plus faible : le <math>C_x</math> a donc bien un rôle de freinage...

== Gravité subjective et forme ==

Dans le référentiel de la goutte,
l'accélération de la gravité est diminuée
de l'accélération propre de la goutte.
L'écart correspondant de pression hydrostatique
d'une extrémité à l'autre de la goutte vaut donc :

<math>\Delta p=\rho_{\rm liq}\,\left(
g-\frac{{\rm d}V}{{\rm d}t}\right)\,R</math>

À cause de la [[pression de Laplace]], cet écart de pression correspond 
à un écart de courbure, comme il a été dit 
(<math>\Delta p=\gamma\,C</math>).


== Forme et vitesse ==

Des équations ci-dessus, il résulte
que tout au cours du temps :

<math>\gamma\,C
=\rho_{\rm air}\,V^2</math>

Cette équation donne la déformation de la goutte par compétition entre les effets de tension superficielle et les effets de frottement dans l'air. Ainsi, une forte déformation de la goutte est atteinte
lorsque <math>C\simeq 1/R</math>.
En réalité, cela est vrai uniquement 
en l'absence d'oscillations de la goutte,
comme nous verrons plus loin.

Dans le régime de chute libre, lorsque la goutte a atteint son régime de vitesse terminale, cette équation se réécrit :

<math>\gamma\,C
=\rho_{\rm liquide}\,R\,g</math>

Si on considère que la goutte est déformée à partir du moment où la courbure est de l'ordre de <math> 1/R</math>, cela donne une taille caractéristique de la goutte en dessous de laquelle elle reste sphérique :

<math>R_c=\sqrt{\frac{\gamma}{\rho g}} </math>
Cette longueur est, un peu fortuitement, la [[longueur capillaire]].<br>

[[File:Eclatement de goutte de pluie.gif|thumb|Éclatement des grosses gouttes de pluie.]]

Par conséquent, une goutte de pluie qui fait moins de {{Unité|3|mm}} de diamètre (longueur capillaire de l'eau) aura une forme complètement sphérique. Une goutte plus grosse pourra prendre une forme plus étrange (cacahuète, parachute ou sac à bourrelet...)(animation ci-contre). Dans ce dernier cas, l'air qui s'engouffre à l'intérieur du sac fait bientôt éclater la goutte en une multitude de gouttelettes de taille inférieure à {{Unité|3|mm}}<ref>Les auteurs de MICROPHYSICS OF CLOUDS AND PRECIPITATION 
[https://vdoc.pub/download/microphysics-of-clouds-and-precipitation-7hiplhj3pdp0] commentent ainsi les expérimentations : "Tombant dans de l'air calme, les gouttes peuvent être aussi grosse que 4,5 mm (en diamètre équivalent) avant d'exploser". D'autres auteurs indiquent qu’il n’existe pas de gouttes de pluie de diamètre supérieur à 6 mm.</ref>.{{clr}}

== Collision avec le sol ou un objet ==
Quand la goutte touche le sol ou une surface horizontale, elle s’aplatit et éventuellement éclate ou rebondit en créant de minuscules gouttes d'[[aérosol]]s, facilement inhalées et qui seraient responsables de l'[[odeur]] développée par certaines pluies et éventuellement de la libération de virus et bactéries à partir du sol ou de l'air<ref>[https://www.youtube.com/watch?v=Waqmq_GTyjA#t=47  ''Rainfall can release aerosols, high-speed video shows'' ] vidéo postée par le MIT (Massachusetts Institute of Technology), à la suite d'un travail présenté dans la revue ''Nature Communication''</ref>.

== Notes et références ==
{{Références}}

== Voir aussi ==
=== Bibliographie ===
* É. Reyssat, F. Chevy, A.-L. Biance, L. Petitjean and D. Quéré, "Shape and instability of free-falling liquid globules",'' EPL'', '''80''', 34005 (2007)

=== Articles connexes ===
* [[Goutte (physique)]]
* [[Distribution des gouttes de pluie]]
* [[Capteur de gouttelettes]]
* [[Pluie]]
* [[Molécule d'eau]]

{{Portail|Physique|Eau}}

[[Catégorie:Dynamique des fluides]]
[[Catégorie:Eau]]