Anisotropy of emergent large-scale dynamics in forced stratified shear flows

Abstract: Stably stratified shear flows, where the base velocity shear is quasi-continuously forced externally, can arise in many geophysically relevant circumstances. It is important to determine the emergent dynamics of the ensuing statistically steady stratified turbulence. We investigate this phenomenon in a series of three-dimensional direct numerical simulations using spectral element methods. We force the flow to relax back towards vertical hyperbolic tangent profiles of streamwise velocity and buoyancy, with characteristic half-depth $\dO$, half-velocity jump $\UO$, and half-buoyancy jump $\BO$, with a relaxation time $t_r=100 \tauadv$ where $\tauadv:=\dO/\UO$. We consider computational domains of vertical extent $\Gz=48$ with a range of horizontal extents $16\leq\Gh\leq512$. We simulate a fluid with Prandtl number $\Pr:=\nu/\kappa=1$, and set the initial bulk Reynolds number $\ReO:=\UO\dO/\nu=50$ and Richardson number $\RiO:=\BO\dO/\UO^{2}=1/80$. At these parameters, the flow is initially unstable to a primary Kelvin-Helmholtz instability. We simulate the continuously forced flow over about $5000\tauadv$, and investigate the dynamically emergent length scales and turbulence properties of the statistically stationary flow, in particular the local turbulent flux coefficient. We find that the shear layer half depth converges to $d\approx8$ with markedly increased $\Re\approx400$ and associated convergent vertical mixing properties only for $\Gh\gtrsim\Ghcrit=96$. However, emergent dominant yet large-scale spanwise or streamwise flow structures appear to extend up to $\Lambda_{y}\approx50$ or even $\Lambda_{x}\approx115$, respectively. Our observations demonstrate the marked anisotropy of characteristic emergent length scales, and are consistent with the possibility that an `imprint' of the primary instability continues to survive in such turbulent flows.