Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations

Laboratory experiments on stably stratified grid turbulence have suggested that turbulent diffusivity κ_ρ can be expressed in terms of a turbulence activity parameter ε/ν N² , with different power-law relations appropriate for different levels of ε/ν N² . To further examine the applicability of these findings to both a wider range of the turbulence intensity parameter ε/ν N² and different forcing mechanisms, DNS data of homogeneous sheared stratified turbulence generated by Shih et al. (2000) and Venayagamoorthy et al. (2003) are analysed in this study. Both scalar eddy diffusivity κ_ρ and eddy viscosity κ_ν are found to be well-correlated with ε/ν N² , and three distinct regimes of behaviour depending on the value of ε/ν N² are apparent. In the diffusive regime D , corresponding to low values of ε/ν N² and decaying turbulence, the total diffusivity reverts to the molecular value; in the intermediate regime I , corresponding to 7 \,< ε/ν N² \,<\, 100 and stationary turbulence, diffusivity exhibits a linear relationship with ε/ν N² , as predicted by Osborn (1980); finally, in the energetic regime E , corresponding to higher values of ε/ν N² and growing turbulence, the diffusivity scales with (ε/ν N²)1/2 . The dependence of the flux Richardson number Rf on explains the shift in power law between regimes I and E . Estimates for the overturning length scale and velocity scales are found for the various ε/ν N² regimes. It is noted that ε/ν N² \,~\, Re/ Ri\,~\, ReFr² , suggesting that such Reynolds–Richardson number or Reynolds–Froude number aggregates are more descriptive of stratified turbulent flow conditions than the conventional reliance on Richardson number alone.