Jump to content

Kelvin–Helmholtz instability

From Wikipedia, the free encyclopedia
Numerical simulation of a temporal Kelvin–Helmholtz instability

The Kelvin–Helmholtz instability (after Lord Kelvin and Hermann von Helmholtz) is a fluid instability that occurs when there is velocity shear in a single continuous fluid or a velocity difference across the interface between two fluids. Kelvin-Helmholtz instabilities are visible in the atmospheres of planets and moons, such as in cloud formations on Earth or the Great Red Spot on Jupiter, and the atmosphere of the Sun.[1][2]

Spatially developing 2D Kelvin-Helmholtz instability at low Reynolds number. Small perturbations, imposed at the inlet on the tangential velocity, evolve in the computational box. High Reynolds number would be marked with an increase of small scale motions.

Theory overview and mathematical concepts

[edit]
A KH instability rendered visible by clouds, known as fluctus,[3] over Mount Duval in Australia
A KH instability on the planet Saturn, formed at the interaction of two bands of the planet's atmosphere
Kelvin-Helmholtz billows 500m deep in the Atlantic Ocean
Animation of the KH instability, using a second order 2D finite volume scheme
Kelvin-Helmholtz cloud formation over Hartford, Connecticut
KH instability in the Sun's atmosphere, photographed on April 8, 2010.[1][4]

Fluid dynamics predicts the onset of instability and transition to turbulent flow within fluids of different densities moving at different speeds.[5] If surface tension is ignored, two fluids in parallel motion with different velocities and densities yield an interface that is unstable to short-wavelength perturbations for all speeds. However, surface tension is able to stabilize the short wavelength instability up to a threshold velocity.

Suppose we have two parallel plane streams of fluid of same density with speeds U1 and U2. Then the perturbed interface with the normal-mode decomposition ei(kx−ωt) obeys the dispersion relation (for ideal flow with no surface tension),

ω=kρ1U1+ρ2U2ρ1+ρ2+iiρ1ρ2ρ1+ρ2|U1−U2|.

If the density and velocity vary continuously in space (with the lighter layers uppermost, so that the fluid is RT-stable), the dynamics of the Kelvin-Helmholtz instability is described by the Taylor–Goldstein equation: (U−c)2(d2ϕ~dz2−k2ϕ~)+[N2−(U−c)d2Udz2]ϕ~=0, where N=g/Lρ denotes the Brunt–Väisälä frequency, U is the horizontal parallel velocity, k is the wave number, c is the eigenvalue parameter of the problem, ϕ~ is complex amplitude of the stream function. Its onset is given by the Richardson number Ri. Typically the layer is unstable for Ri<0.25. These effects are common in cloud layers. The study of this instability is applicable in plasma physics, for example in inertial confinement fusion and the plasma–beryllium interface. In situations where there is a state of static stability (where there is a continuous density gradient), the Rayleigh–Taylor instability is often insignificant compared to the magnitude of the Kelvin–Helmholtz instability.

Numerically, the Kelvin–Helmholtz instability is simulated in a temporal or a spatial approach. In the temporal approach, the flow is considered in a periodic (cyclic) box "moving" at mean speed (absolute instability). In the spatial approach, simulations mimic a lab experiment with natural inlet and outlet conditions (convective instability).


Discovery and history

[edit]

The existence of the Kelvin-Helmholtz instability was first discovered by German physiologist and physicist Hermann von Helmholtz in 1868. Helmholtz identified that "every perfect geometrically sharp edge by which a fluid flows must tear it asunder and establish a surface of separation".[6][5] Following that work, in 1871, collaborator William Thomson (later Lord Kelvin), developed a mathematical solution of linear instability whilst attempting to model the formation of ocean wind waves.[7]

Throughout the early 20th Century, the ideas of Kelvin-Helmholtz instabilities were applied to a range of stratified fluid applications. In the early 1920s, Lewis Fry Richardson developed the concept that such shear instability would only form where shear overcame static stability due to stratification, encapsulated in the Richardson Number.

Geophysical observations of the Kelvin-Helmholtz instability were made through the late 1960s/early 1970s, for clouds,[8] and later the ocean.[9]


Kelvin–Helmholtz instability in compressible flows: Landau stabilization

[edit]

The Kelvin–Helmholtz instability in compressible flows was first studied by Lev Landau in 1944,[10] and was later refined by Sergei Ivanovich Syrovatskii in 1954.[11] Suppose U2=0 and U1=U (which can always be achieved by a suitable choice of the frame of reference); then the dispersion relation for the normal mode ei(kxx+kyy−ωt) is given by

[1ω2−1(ω−Ukx)2][1c2k2−1ω2−1(ω−Ukx)2]=0,k=kx2+ky2,

where c is the sound speed. The zero of the first factor is ω=12Ukx, which is always real. The zeros of the second factor are given by solving the bi-quadratic form in terms of (ω−12Ukx):

ω=12Ukx±14U2kx2+c2k2±c2k21+U2kx2c2k2.

These roots are real only if Ukx/ck>22≈2.83, as shown elegantly by Landau. Let ϕ be the angle between the velocity flow U𝐞x and the wavevector 𝐤=kx𝐞x+ky𝐞y. Since kx=kcos⁡ϕ, the instability condition can be rewritten as

Mcos⁡ϕ<22,M=Uc;

here M is the Mach number. The above condition yields two distinct physical regimes based on the Mach number of the flow:[12]

  • If the flow velocity is relatively low, specifically M<22, the inequality holds for all angles, meaning the flow is unstable to perturbations propagating in any direction along the interface.
  • If the flow is highly supersonic such that M≥22, the interface becomes completely stabilized against directional perturbations traveling within a forward cone where cos⁡ϕ≥22/M; this is known as the Landau stabilization. In this supersonic regime, the instability only appears for oblique waves propagating at sufficiently large angles relative to the flow, satisfying cos⁡ϕ<22/M.

Landau in his original 1944 paper considered only the case where ϕ=0 (i.e., ky=0), leading to the prediction that the interface becomes completely stable at high supersonic speeds. However, Yakov Zeldovich subsequently conducted a laboratory experiment testing this phenomenon and found that the interface remains highly unstable and transitions into turbulence. Sergei I. Syrovatskii's 1954 generalization for arbitrary values of ϕ provided the missing mathematical framework, proving that oblique waves (ky≠0) will always grow, thereby explaining the experimental instability observed by Zeldovich.

See also

[edit]

Notes

[edit]
  1. 1 2 Fox, Karen C. (30 December 2014). "NASA's Solar Dynamics Observatory Catches "Surfer" Waves on the Sun". NASA-The Sun-Earth Connection: Heliophysics. NASA. Archived from the original on 20 November 2021. Retrieved 21 July 2011.
  2. ↑ Inyama, Fechi (October 2026). "The Sun's Swirlies". ScienceNews. 208 (10): 8–9. Archived from the original on 2026-09-25. Retrieved 2026-09-25.
  3. ↑ Sutherland, Scott (March 23, 2017). "Cloud Atlas joins the digital age while adding a dozen new cloud types". The Weather Network. Pelmorex Media. Retrieved 13 August 2025.
  4. ↑ Ofman, L.; Thompson, B. J. (2011-06-01). "SDO/AIA Observation of Kelvin-Helmholtz Instability in the Solar Corona". The Astrophysical Journal. 734 (1): L11. arXiv:1101.4249. Bibcode:2011ApJ...734L..11O. doi:10.1088/2041-8205/734/1/L11. ISSN 0004-637X. Archived from the original on 2023-07-14. Retrieved 2023-07-14.
  5. 1 2 Drazin, P. G. (2003). Encyclopedia of Atmospheric Sciences. Elsevier Ltd. pp. 1068–1072. doi:10.1016/B978-0-12-382225-3.00190-0.
  6. ↑ von Helmholtz, Hermann (1 November 1868). "XLIII. On discontinuous movements of fluids". The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science. 36 (244): 337–346. doi:10.1080/14786446808640073.
  7. ↑ Matsuoka, Chihiro (2014-03-31). "Kelvin-Helmholtz Instability and Roll-up". Scholarpedia. 9 (3) 11821. Bibcode:2014SchpJ...911821M. doi:10.4249/scholarpedia.11821. ISSN 1941-6016.
  8. ↑ Ludlam, F. H. (October 1967). "Characteristics of billow clouds and their relation to clear-air turbulence". Quarterly Journal of the Royal Meteorological Society. 93 (398): 419–435. Bibcode:1967QJRMS..93..419L. doi:10.1002/qj.49709339803.
  9. ↑ Woods, J. D. (18 June 1968). "Wave-induced shear instability in the summer thermocline". Journal of Fluid Mechanics. 32 (4): 791–800. Bibcode:1968JFM....32..791W. doi:10.1017/S0022112068001035. S2CID 67827521.
  10. ↑ Landau, L. (1944). Stability of tangential discontinuities in compressible fluid. In Dokl. Akad. Nauk SSSR (Vol. 44, No. 4, pp. 151-153).
  11. ↑ Syrovatskii, S. I. (1954). INSTABILITY OF TANGENTIAL DISCONTINUITIES IN A COMPRESSIBLE MEDIUM.
  12. ↑ Gerwin, R. A. (1968). Stability of the interface between two fluids in relative motion. Reviews of Modern Physics, 40(3), 652.

References

[edit]
[edit]