The intermittency boundary in stratified plane Couette flow

We study stratified turbulence in plane Couette flow using direct numerical simulations. Two external dimensionless parameters control the dynamics, the Reynolds number Re=Uh/ and the bulk Richardson number Ri=gₕTh/U^2, where U and T are half the velocity and temperature difference between the two walls respectively, h is the half channel depth, is the kinematic viscosity and gₕ is the buoyancy parameter. We focus on spatio-temporal intermittency due to stratification and we explore the boundary between fully developed turbulence and intermittent flow in the Re-Ri plane. The structures populating the intermittent flow regime show coexistence between laminar and turbulent patches, and we demonstrate that there are qualitative differences between the previously studied low- Re low- Ri intermittent regime and the high- Re high- Ri intermittent regime. At low- Re low- Ri, turbulent regions span the entire gap, whereas at high- Re high- Ri, turbulence is confined vertically with complex dynamics arising from interacting turbulent layers. Consistent with a previous investigation of Flores & Riley (Boundary-Layer Meteorol. , vol. 129 (2), 2010, pp. 241–259), we present evidence suggesting that intermittency in the asymptotic regime of high- Re Couette flows appears for L^+<200, where L^+=Lu/, with L being the Monin–Obukhov length scale, L=u^3/Cqₖ, qₖ the wall heat flux, C the von Kármán constant and u=ₖ/{₀} the friction velocity determined from the wall shear stress ₖ, where ₀ is the constant background density. We also consider the mixing as quantified by various versions of the flux Richardson number Ri₅, defined as the ratio of the conversion rate from kinetic to potential energy to the turbulent kinetic energy injection rate due to shear. We investigate how laminar and turbulent regions separately contribute to the overall mixing. Remarkably, we find that although fluctuations are greatly suppressed in the laminar regions, Ri₅ does not change significantly compared with its value in turbulent regions. As we observe a tight coupling between the mean temperature and velocity fields, we demonstrate that both Monin–Obukhov self-similarity theory (Monin & Obukhov, Contrib. Geophys. Inst. Acad. Sci. USSR, vol. 151, 1954, pp. 163–187) and the explicit algebraic model of Lazeroms et al. (J. Fluid