In certain geophysical contexts such as lava lakes and mantle convection, a cold, viscous boundary layer forms over a deep pool. The following model problem investigates the buoyant instability of the layer. Beneath a shear-free horizontal boundary, a thin layer (thickness d 1 ) of very viscous fluid overlies a deep layer of less dense, much less viscous fluid; inertia and surface tension are negligible. After the initial unstable equilibrium is perturbed, a long-wave analysis describes the growth of the disturbance, including the nonlinear effects of large amplitude. The results show that nonlinear effects greatly enhance growth, so that initial local maxima in the thickness of the viscous film grow to infinite thickness in finite time, with a timescale 8µ/∆ρ gd 1 . In the final catastrophic growth the peak thickness is inversely proportional to the remaining time. (A parallel analysis for fluids with power-law rheology shows similar catastrophic growth.) While the small-slope approximation must fail before this singular time, the failure is only local, and a similarity solution describes how the peaks become downwelling plumes as the viscous film drains away.
BUOYANT INSTABILITY OF A VISCOUS FILM OVER A PASSIVE FLUID
INTRODUCTIONThis work examines the strongly nonlinear effects of finite amplitude in the RayleighTaylor instability of a horizontal viscous film under a shear-free boundary and over a much less viscous fluid. Inertia and surface tension are neglected, and, in the parameter range considered, the motion is limited by normal stresses in the more viscous fluid. The analysis exploits the fact that the most unstable wavelengths are long compared to the thickness of the film. The results show how the growth of disturbances to the interface becomes greatly enhanced when the disturbance amplitude becomes large, leading to the formation of downwelling sheets or plumes in a finite time.The motivation for this problem comes from certain geophysical situations, particularly the stability of the Earth's lithospheric plates. In simplified terms, the oceanic lithosphere (tectonic plates) can be considered a cold, stiff thermal boundary layer above the convecting mantle. Where two plates come together, one subducts under the other and flows downward due to its negative buoyancy. The question of how a new subduction zone is formed, how one large plate may break into two and thus allow some of the dense material to flow back down into the mantle, is not yet resolved. Other closely related geophysical situations include the surfaces of lava lakes, thermal convection in the mantles of other planets, and possibly convection in the Earth's solid core.This work examines a simple model of one possible mechanism for the initiation of subduction: the Rayleigh-Taylor instability. In this model, the lithosphere and the mantle are treated as distinct, highly viscous fluids, the lithosphere being denser (and much more viscous) than the mantle. In this unstable configuration, any variations in the lithosphere thickness t...