Instabilité de Rayleigh-Taylor

L’instabilité de Rayleigh–Taylor, nommée en hommage aux physiciens britanniques Lord Rayleigh et G. I. Taylor, est une instabilité de l’interface séparant deux fluides de densités différentes, qui résulte de la poussée du fluide le plus lourd sur le fluide le plus léger (l'accélération dans le cas d'un système dynamique ou la gravité pour un système initialement statique est dirigée vers la phase légère)[1],[2]. Ce phénomène est produit par exemple par l'onde de choc à l'origine des nuages interstellaires. Dans ce cas particulier où le choc est à l'origine de la mise en vitesse du système, on parlera d'instabilité de Richtmyer-Meshkov. Il se produit une situation analogue lorsque la gravité affecte deux fluides de densités différentes (le fluide le plus dense se trouvant au-dessus du fluide le moins dense) comme de l'huile minérale à la surface de l'eau[2].

Manifestation de l'instabilité de Rayleigh–Taylor dans la nébuleuse du Crabe.

Considérons deux couches de fluides immiscibles superposées dans deux plans parallèles, la plus lourde surplombant la plus légère et toutes deux soumises à la pesanteur terrestre. L’équilibre est instable à la moindre perturbation : toute perturbation va s'amplifier et libérer de l’énergie potentielle, le fluide le plus lourd gagnant progressivement la moitié inférieure sous l'effet du champ de gravitation, et le fluide léger passe au-dessus. C'est cette configuration qu'a étudiée Lord Rayleigh[2]. La découverte importante de G. I. Taylor a consisté à montrer que cette situation est équivalente à celle qui se produit lorsque les fluides (hors de toute gravité) sont accélérés, le fluide léger étant propulsé à l'intérieur du fluide le plus lourd[2]. Cela se produit notamment lorsque l'on projette un verre à terre avec une accélération supérieure à la pesanteur terrestre[2] g.

Lorsque l’instabilité développe ses effets, des irrégularités (« fossettes ») se propagent vers le bas en polypes de Rayleigh–Taylor qui finissent même par se mélanger. C'est pourquoi on qualifie parfois l’instabilité de Rayleigh–Taylor d’instabilité à traines (fingering instability)[3]. Le fluide le plus léger s'expand vers le haut comme un champignon nucléaire[4],[5].

On observe ce phénomène dans plusieurs situations courantes, non seulement dans les dômes salins ou les couches d’inversion, mais aussi en astrophysique et en électrocinétique. Les polypes de Rayleigh-Taylor sont particulièrement visibles dans la Nébuleuse du Crabe, où le plérion engendré par le pulsar du Crabe déborde les projections issues de l’explosion de la supernova il y a 1 000 ans[6].

Il ne faut pas confondre l’instabilité de Rayleigh–Taylor avec l’instabilité de Plateau-Rayleigh (parfois appelée « instabilité du tuyau d'arrosage ») : cette dernière, qui se produit dans les jets de liquide, est due à la tension superficielle, qui tend à disperser un jet cylindrique en une projection de gouttelettes de même volume mais de surface spécifique moindre.

Analyse linéaire de la stabilité

modifier

L’instabilité bidimensionnelle non-visqueuse de Rayleigh–Taylor constitue un excellent banc d'essai pour l'étude mathématique de la stabilité du fait de la nature extrêmement simple de la configuration initiale[7], décrite par un champ de vitesse moyenne tel que   où le champ gravitationnel est   Une interface en   sépare les fluides de densités   dans la zone supérieure, et   dans la zone inférieure. On montre que dans cette section, lorsque le fluide le plus lourd se trouve au-dessus, la moindre perturbation de l’interface s’amplifie exponentiellement, avec le taux[2]

 

  est le taux de croissance,   est le nombre d'onde spatial et   est le Nombre d'Atwood.

 
Simulation hydrodynamique d’un « polype » isolé de l’instabilité de Rayleigh–Taylor[9]. Observez la formation d'une instabilité de Kelvin-Helmholtz, à partir de la deuxième image (se développant à partir du niveau  ), ainsi que la formation d’un chapeau de champignon à un stade évolué dans les troisième et quatrième images de cette séquence.

Le temps caractéristique de croissance de la surface libre   initialement en   est donné par :

 

qui croît exponentiellement avec le temps. Ici B désigne l’amplitude de la perturbation initiale, et   est la partie réelle de l’expression complexe entre parenthèses.

En général, la condition pour que l’instabilité soit linéaire est que la partie imaginaire de la célérité complexe c soit positive. Finalement, le rétablissement de la tension superficielle diminue c2 en module et a donc un effet stabilisant. En effet, il existe un domaine d'ondes courtes pour lesquelles la tension superficielle stabilise le système et empêche l’instabilité.

Comportement à long terme

modifier

L’analyse qui précède n'est plus valable quand on a affaire à une perturbation de grande amplitude : dans ce cas, la croissance est non-linéaire, les polypes et les bulles s'entrelacent et s'enroulent en tourbillons. Comme l'illustre la figure ci-contre, il faut recourir à la simulation numérique[10], pour décrire mathématiquement le système.

Notes et références

modifier
  1. D.H. Sharp, « An Overview of Rayleigh-Taylor Instability », Physica D, vol. 12,‎ , p. 3–18 (DOI 10.1016/0167-2789(84)90510-4)
  2. a b c d e et f Drazin (2002) p. 50–51.
  3. H. B. Che, B. Hilko et E. Panarella, « The Rayleigh–Taylor instability in the spherical pinch », Journal of Fusion Energy, vol. 13, no 4,‎ , p. 275–280 (DOI 10.1007/BF02215847)
  4. (en) C.-Y. Wang et R. A. Chevalier, « Instabilities and Clumping in Type Ia Supernova Remnants », version v1, .
  5. (en) R. J. Tayler (dir.), W. Hillebrandt et P. Höflich, Stellar Astrophysics, Supernova 1987a in the Large Magellanic Cloud, Bristol/Philadelphia, CRC Press, , 356 p. (ISBN 0-7503-0200-3), p. 249–302 : cf. page 274.
  6. J. Jeff Hester, « The Crab Nebula: an Astrophysical Chimera », Annual Review of Astronomy and Astrophysics, vol. 46,‎ , p. 127–155 (DOI 10.1146/annurev.astro.45.051806.110608)
  7. a et b Drazin (2002) p. 48–52.
  8. On trouve un calcul similaire dans l'ouvrage de Chandrasekhar (1981), §92, pp. 433–435.
  9. Shengtai Li et Hui Li, « Parallel AMR Code for Compressible MHD or HD Equations », Los Alamos National Laboratory (consulté le )
  10. IUSTI, « Simulation numérique des instabilités de Richtmyer-Meshkov », (consulté le )

Voir aussi

modifier

Bibliographie

modifier

Sources historiques

modifier

Bibliographie récente

modifier

Articles connexes

modifier

Liens externes

modifier