Turbulent shear layers in a uniformly stratified background: DNS at high Reynolds number

Abstract

Direct numerical simulations are performed to investigate a stratified shear layer at high Reynolds number in a study where the Richardson number is varied among cases. Unlike previous work on a two-layer configuration in which the shear layer resides between two layers with constant density, an unbounded fluid with uniform stratification is considered here. The evolution of the shear layer includes a primary Kelvin-Helmholtz shear instability followed by a wide range of secondary shear and convective instabilities, similar to the two-layer configuration. During transition to turbulence, the shear layers at low exhibit a period of thickness contraction (not observed at lower) when the momentum and buoyancy fluxes are counter-gradient. The behaviour in the turbulent regime is significantly different from the case with a two-layer density profile. The transition layers, which are zones with elevated shear and stratification that form at the shear-layer edges, are stronger and also able to support a significant internal wave flux. After the shear layer becomes turbulent, mixing in the transition layers is shown to be more efficient than that which develops in the centre of the shear layer. Overall, the cumulative mixing efficiency is larger than the often assumed value of 1/6. Also, is found to be smaller than that in the two-layer configuration at moderate Ri. It is relatively less sensitive to background stratification, exhibiting little variation for. The dependence of mixing efficiency on buoyancy Reynolds number during the turbulence phase is qualitatively similar to homogeneous sheared turbulence.