Sediment-laden fresh water above salt water: linear stability analysis

Abstract

AbstractWhen a layer of particle-laden fresh water is placed above clear, saline water, both Rayleigh–Taylor and double diffusive fingering instabilities may arise. For quasi-steady base profiles, we obtain linear stability results for such situations by means of a rational spectral approximation method with adaptively chosen grid points, which is able to resolve multiple steep gradients in the base state density profile. In the absence of salinity and for a step-like concentration profile, the dominant parameter is the ratio of the particle settling velocity to the viscous velocity scale. As long as this ratio is small, particle settling has a negligible influence on the instability growth. However, when the particles settle more rapidly than the instability grows, the growth rate decreases inversely proportional to the settling velocity. This damping effect is a result of the smearing of the vorticity field, which in turn is caused by the deposition of vorticity onto the fluid elements passing through the interface between clear and particle-laden fluid. In the presence of a stably stratified salinity field, this picture changes dramatically. An important new parameter is the ratio of the particle settling velocity to the diffusive spreading velocity of the salinity, or alternatively the ratio of the unstable layer thickness to the diffusive interface thickness of the salinity profile. As long as this quantity does not exceed unity, the instability of the system and the most amplified wavenumber are primarily determined by double diffusive effects. In contrast to situations without salinity, particle settling can have a destabilizing effect and significantly increase the growth rate. Scaling laws obtained from the linear stability results are seen to be largely consistent with earlier experimental observations and theoretical arguments put forward by other authors. For unstable layer thicknesses much larger than the salinity interface thickness, the particle and salinity interfaces become increasingly decoupled, and the dominant instability mode becomes Rayleigh–Taylor-like, centred at the lower boundary of the particle-laden flow region.