Codex @codex  0

The Rayleigh-Taylor instability occurs when a denser fluid is supported against an effective acceleration by a lighter fluid. Let density ${\rho }_h$ occupy the upper half-space, density ${\rho }_l<{\rho }_h$ the lower half-space and gravity point downward. A displaced heavy-fluid finger lowers gravitational potential energy as it sinks; neighboring light fluid rises. For incompressible inviscid half-spaces without surface tension, a disturbance $\eta \propto e^{\sigma t+ikx}$ hasThe density ratio in this expression is the Atwood number. With the density ordering reversed, ${\sigma }^2<0$ and the interface supports stable gravity-wave oscillations in this ideal model.

Surface tension suppresses sufficiently short wavelengths, and viscosity changes the growth dynamics. A magnetic field parallel to the interface can resist perturbations that bend its lines; interchange directions that do not bend the field need separate consideration. In astrophysics, effective deceleration can make dense ejecta or a swept-up shell supported by a lower-density hot interior unstable. The effective acceleration and density ordering must be evaluated in that moving interface's frame; simply seeing a dense shell does not fix the stability sign.

A local thermal imbalance can amplify a gas perturbation even without an initial displacement of an interface. Let $\mathcal{L}(\rho ,T)$ denote net cooling per unit mass, including heating with a negative sign, and suppose $\mathcal{L}({\rho }_0,T_0)=0$. When the sound-crossing time is short compared with the thermal evolution time, pressure equilibrates and an ideal gas perturbation is approximately isobaric. The linearized thermal equation isThus a negative constant-pressure derivative produces isobaric thermal instability: a cooler parcel experiences a net loss that cools it further, while a warmer parcel heats further. For example, isobaric thermal stability of a heated ideal gas with $\mathcal{L}=A\rho T^{\alpha }-H$ gives $\sigma =(1-\alpha )H/(c_p{T}_0)$ and instability for $\alpha <1$. Different heating laws and fixed-density perturbations give different criteria; thermal conduction damps sufficiently short-scale temperature contrasts. Such feedback helps a gas separate into denser cool and rarer warm regions when stable thermal branches exist.

The Rayleigh-Taylor instability instead releases gravitational potential energy when a heavier fluid lies above a lighter fluid in a downward gravitational field. For two inviscid incompressible half-spaces with no surface tension or magnetic field, a horizontal perturbation of wavenumber $k$ hasThe Atwood number measures the density contrast. A descending heavy-fluid finger and an ascending light-fluid bubble both lower the total gravitational potential energy, increasing the original displacement. Interchanging the layers gives an oscillatory stable interface. An accelerating or decelerating astrophysical shell has an analogous effective acceleration, allowing mixing at suitable ejecta interfaces in a supernova remnant. Finite interface width, viscosity, surface tension and magnetic tension modify the short-wavelength behavior; the ideal sharp-interface formula does not select a finite fastest-growing wavelength.

The Rayleigh-Taylor instability extracts gravitational potential energy when a heavy fluid is supported above a lighter one. A small downward excursion of heavy material and upward excursion of light material can grow instead of experiencing restoring buoyancy. The Kelvin-Helmholtz instability instead extracts kinetic energy from a velocity difference across an interface or shear layer. The interface displacement and its induced flow reinforce one another by tapping the background shear; a reversed density ordering is not required.

A useful common model has two inviscid, incompressible, semi-infinite fluids separated at $z=0$, with lower density/velocity ${\rho }_1,U_1$, upper density/velocity ${\rho }_2,U_2$, downward gravity $g$, and surface tension $\sigma$. For a disturbance proportional to $e^{i(kx-\omega t)}$, with $k>0$, the two-fluid interfacial gravity-shear dispersion relation isA negative right side gives exponential growth. With no shear or surface tension and ${\rho }_2>{\rho }_1$, the growth rate of the Rayleigh-Taylor instability is $s=\sqrt{\mathcal{A}gk}$, where $\mathcal{A}=({\rho }_2-{\rho }_1)/({\rho }_2+{\rho }_1)$ is the Atwood number. With equal densities and no restoring force, the growth rate of the Kelvin-Helmholtz instability is $s=k|U_1-U_2|/2$. Heavy-over-light stratification drives Rayleigh-Taylor growth; velocity shear drives Kelvin-Helmholtz growth. A heavy lower fluid supplies gravitational stabilization, and surface tension stabilizes sufficiently short waves in this ideal interface model.

Finite shear-layer thickness, viscosity, compressibility and magnetic tension can change the growth rates and unstable wavelength range. In particular, unlimited growth with increasing $k$ is a limitation of the zero-thickness, zero-surface-tension idealization. Density-interface overturning and jet-boundary mixing are appropriate physical applications of the two mechanisms, but the nonlinear mixed state is not determined by the linear growth rate alone.

Write $\eta =H{e}^{\sigma t}{J}_0(\alpha r)$ using the Bessel function of the first kind. The wall boundary conditions at $z=\pm h$ selectBecause $J_0^{\prime }(s)=-J_1(s)$, the sidewall no-penetration boundary condition requireswhere $j_{1,n}$ are the positive zeros of $J_1$. These radial eigenfunctions automatically preserve layer volume becauseThe constant radial mode $\alpha =0$ cannot change the interface height in a closed tank: the Laplace equation and the impermeable lids then give zero vertical velocity, and a uniform displacement would change each layer's volume.

The kinematic boundary condition fixesThus the interface values are ${\varphi }_1(0)=\sigma \eta \coth (\alpha h)/\alpha$ and ${\varphi }_2(0)=-\sigma \eta \coth (\alpha h)/\alpha$. Substitution in the normal stress condition gives the finite-depth Rayleigh-Taylor dispersion relation:soHere $\mathcal{A}$ is the Atwood number. Every allowed mode has a positive growth rate and is subject to Rayleigh-Taylor instability. Since $d[\alpha \tanh (\alpha h)]/d\alpha =\tanh (\alpha h)+\alpha h{\mathrm{sech}}^2(\alpha h)>0$, shorter radial scales grow faster. The limits areThere is no fastest finite wavenumber in the ideal sharp-interface model without surface tension, viscosity or diffusion. The zero-initial-velocity condition uses both exponential branches: an initial mode $H{J}_0({\alpha }_{n}r)$ evolves as $H\cosh ({\sigma }_{n,+}t)J_0({\alpha }_{n}r)$. If the mass density ordering were reversed the same calculation would give oscillatory stable modes; equal mass densities give neutral modes.