Past exam of the astrophysical fluid dynamics course of the University of Cambridge / 2014 / Paper 1 / 2 / c

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.