Double-diffusive sedimentation

Eckart Heinz Meiburg and Peter Burns, University of California Santa Barbara, Santa Barbara, CA, United States
Abstract:
When a layer of particle-laden fresh water is placed above clear, saline water, both double-diffusive and Rayleigh-Taylor instabilities may arise. The present investigation extends the linear stability analysis of Burns and Meiburg (2012) into the nonlinear regime, by means of two- and three-dimensional direct numerical simulations (DNS). The vigorous growth of individual fingers gives rise to a secondary instability, and eventually to the formation of intense plumes that become detached from the interfacial region. The simulations show that the presence of particles with a Stokes settling velocity modifies the traditional double-diffusive fingering by creating an unstable `nose region' in the horizontally averaged profiles, located between the upward moving salinity and the downward moving sediment interface. The effective thickness ls (lc) of the salinity (sediment) interface grows diffusively, as does the height H of the nose region. The ratio H/ls initially grows and then plateaus, at a value that is determined by the balance between the flux of sediment into the rose region from above, the double-diffusive/Rayleigh-Taylor flux out of the nose region below, and the rate of sediment accumulation within the nose region. For small values of H/ls<O(0.1), double-diffusive fingering dominates, while for larger values H/ls>O(0.1) the sediment and salinity interfaces become increasingly separated in space and the dominant instability mode becomes Rayleigh-Taylor-like. A scaling analysis based on the results of a parametric study indicates that H/ls is a linear function of a single dimensionless grouping that can be interpreted as the ratio of in- and outflow of sediment into the nose region. The simulation results furthermore indicate that double-diffusive and Rayleigh-Taylor instability mechanisms cause the effective settling velocity of the sediment to scale with the overall buoyancy velocity of the system, which can be orders of magnitude larger than the Stokes settling velocity.