Instabilité de Kelvin-Helmholtz — Écoulement cisaillé biphasique

Deux couches de fluide glissant l'une contre l'autre se déstabilisent à leur interface en formant les billons caractéristiques de KH. Réglez le rapport de densité, le cisaillement, la tension superficielle et activez la stabilisation magnétique (MHD).

Évolution de l'interface et champ de vitesse

Le fluide inférieur lourd (bleu) glisse contre le fluide supérieur léger (orange). Quand le cisaillement l'emporte sur la gravité, la tension et la raideur magnétique, l'interface s'enroule en billons KH. Les flèches montrent le champ de vitesse perturbé.

Taux de croissance γ(k)

Partie imaginaire de ω. Pic au nombre d'onde le plus instable k* ; l'interface amplifie préférentiellement cette longueur d'onde.

Dispersion ωᵣ(k)

Partie réelle de la fréquence. Là où γ=0 c'est une onde stable gravité-capillarité-Alfvén ; là où γ>0 l'onde est remplacée par une croissance exponentielle.

Le mécanisme de Kelvin-Helmholtz

Quand deux couches de fluide glissent l'une contre l'autre, le saut de vitesse à l'interface est une nappe de vorticité. Toute petite corrugation de l'interface abaisse la pression à la crête (Bernoulli) et soulève le creux, donc la corrugation croît : la classique instabilité KH. Pour des densités égales sans gravité ni tension, toutes les longueurs d'onde sont instables et le taux de croissance augmente avec le nombre d'onde k, donc les ondes les plus courtes gagnent jusqu'à ce que la tension superficielle les coupe. La théorie linéaire donne la relation de dispersion (ρ₁+ρ₂)ω² − 2(ρ₁U₁+ρ₂U₂)kω + (ρ₁U₁²+ρ₂U₂² − B²/μ₀)k² − (ρ₁−ρ₂)gk − σk³ = 0, dont les racines complexes ω = ωᵣ − iγ donnent un taux de croissance positif γ là où le cisaillement domine.

Ce qui stabilise l'interface

Trois effets combattent l'instabilité. (1) Stratification de densité stable — quand le fluide le plus lourd est en dessous (ρ₁>ρ₂), la gravité (ρ₁−ρ₂)gk agit comme force de rappel sur les ondes longues ; c'est pourquoi la surface de la mer ne se brise en vagues qu'au-dessus d'une vitesse de vent critique. (2) La tension superficielle σk³ tue les longueurs d'onde les plus courtes, fixant une coupure à haut k ; c'est pourquoi les petites rides sont lisses. (3) Un champ magnétique horizontal ajoute la tension magnétique B²k²/μ₀ (effet Alfvén), qui peut supprimer entièrement KH — c'est pourquoi la magnétopause du vent solaire peut rester nette. Augmenter ΔU ou le rapport de densité ramène vers l'instabilité.

Où apparaissent les billons KH

L'instabilité KH est omniprésente. Le vent sur l'eau sème les vagues océaniques et les embruns. L'atmosphère striée de Jupiter est une galerie de tourbillons enroulés par KH, y compris la frontière de cisaillement de la Grande Tache Rouge. La magnétopause terrestre — frontière entre le vent solaire et la magnétosphère — est une couche de cisaillement magnétisée dont les ondes KH laissent le plasma du vent solaire s'infiltrer. KH pilote aussi le mélange dans les débordements océaniques, les couches nuageuses cisaillées (le classique « nuage en billon »), les jets d'échappement et l'interface combustible-air dans la combustion scramjet. En astrophysique il façonne les bords de restes de supernova et la turbulence des disques d'accrétion.

Mode d'emploi

Commencez par Cisaillement classique : densités égales, cisaillement pur — toutes les longueurs d'onde sont instables et les billons s'enroulent vite. Ajoutez maintenant la Gravité (ρ₂<ρ₁) pour mimer une stratification stable et regardez les ondes longues se supprimer. Montez la Tension superficielle pour rogner les ondes les plus courtes et voyez la courbe de croissance gagner une coupure à haut k. Activez le Mode MHD et augmentez la tension magnétique : le taux de croissance chute sur tout k, et un champ assez fort stabilise entièrement l'interface (cas magnétopause). Le panneau Taux de croissance met en évidence le mode le plus instable k* ; le panneau Dispersion montre la fréquence d'onde réelle qui survit quand γ=0.