Instability of stationary unbounded stratified fluid

Abstract

Suppose that the density of stationary unbounded viscous fluid is a sinusoidal function of the vertical position coordinate z. Is this body of fluid gravitationally unstable to small disturbances, and, if so, under what conditions, and to what type of disturbance? These questions are considered herein, and the answers are that the fluid is indeed unstable, for any non-zero value of the amplitude of the sine wave, to disturbances with large horizontal wavelength. These disturbances have approximately vertical velocity everywhere and tilt the alternate layers of heavier and of lighter fluid, causing the fluid in the former to slide down and that in the latter to slide up, leading to a sinusoidal variation of the vertically averaged density and thereby to reinforcement of the vertical motion. The identification of this novel and efficient global instability mechanism prompts a consideration of the stability of other cases of unbounded fluid stratified in layers. Two other types of undisturbed density distribution, the first an isolated central layer of heavier or lighter fluid, with density varying say as a Gaussian function, and the second an isolated layer of fluid in which the density varies as the derivative of a Gaussian function, are found to be unstable, at all values of the magnitude of the density variation, to disturbances having the same global character. For the first of these two types of density distribution, the behaviour of a disturbance with long horizontal wavelength depends only on the net excess mass of unit area of the central layer, and for the second it depends only on the first moment of the density in the central layer. For the second type there arises another global instability mechanism in which light fluid is stripped away from one side of the layer and heavy fluid from the other without any tilting. In all cases the properties of a neutral disturbance are determined numerically, and the growth rate is found as a function of the Rayleigh number, the Prandtl number, and the horizontal wavenumber of the disturbance. An energy argument gives results easily for the inviscid non-diffusive limit, when all disturbances grow, and reveals the tilting-sliding mechanism of the instability of a disturbance with large horizontal wavelength in its simplest form.