Abstract:In order that it be correctly characterized, irreversible turbulent mixing in stratified fluids must distinguish between adiabatic ‘stirring’ and diabatic ‘mixing’. Such a distinction has been formalized through the definition of a diapycnal diffusivity, $K_{{\it\rho}}$ (Winters & D’Asaro, J. Fluid Mech. , vol. 317, 1996, pp. 179–193) and an appropriate mixing efficiency, $\mathscr{E}$ (Caulfield & Peltier, J. Fluid Mech. , vol. 413, 2000, pp. 1–47). Equivalent attention has not been paid to the definitions of a corresponding momentum diffusivity $K_{m}$ and hence an appropriately defined turbulent Prandtl number $\mathit{Pr}_{t}=K_{m}/K_{{\it\rho}}$ . In this paper, the diascalar framework of Winters & D’Asaro (1996) is first reformulated to obtain an ‘Osborn-like’ formula in which the correct definition of irreversible mixing efficiency $\mathscr{E}$ is shown to replace the flux Richardson number which Osborn ( J. Phys. Oceanogr. , vol. 10, 1980, pp. 83–89) assumed to characterize this efficiency. We advocate the use of this revised representation for diapycnal diffusivity since the proposed reformulation effectively removes the simplifying assumptions on which the original Osborn formula was based. We similarly propose correspondingly reasonable definitions for $K_{m}$ and $\mathit{Pr}_{t}$ by eliminating the reversible component of the momentum production term. To explore implications of the reformulations for both diapycnal and momentum diffusivity we employ an extensive series of direct numerical simulations (DNS) to investigate the properties of the shear-induced density-stratified turbulence that is engendered through the breaking of a freely evolving Kelvin–Helmholtz wave. The DNS results based on the proposed reformulation of $K_{{\it\rho}}$ are compared with available estimations due to the mixing length model, as well as both the Osborn–Cox and the Osborn models. Estimates based upon the Osborn–Cox formulation are shown to provide the closest approximation to the diapycnal diffusivity delivered by the exact representation. Through compilation of the complete set of DNS results we explore the characteristic dependence of $K_{{\it\rho}}$ on the buoyancy Reynolds number $\mathit{Re}_{b}$ as originally investigated by Shih et al. ( J. Fluid Mech. , vol. 525, 2005, pp. 193–214) in their idealized study of homogeneous stratified and sheared turbulence, and show that the validity of their results is only further reinforced through analysis of the turbulence produced in the more geophysically relevant Kelvin–Helmholtz wave life-cycle ansatz. In contrast to the results described by Shih et al. (2005) however, we show that, besides $\mathit{Re}_{b}$ , a vertically averaged measure of the gradient Richardson number $\mathit{Ri}_{b}$ may equivalently characterize the turbulent mixing at high $\mathit{Re}_{b}$ . Based on the dominant driving processes involved in irreversible mixing, we categorize the intermediate (i.e. $\mathit{Re}_{b}=O(10^{1}{-}10^{2})$ ) and high (i.e. $\mathit{Re}_{b}>O(10^{2})$ ) range of $\mathit{Re}_{b}$ as ‘buoyancy-dominated’ and ‘shear-dominated’ mixing regimes, which together define a transition value of $\mathit{Ri}_{b}\sim 0.2$ . Mixing efficiency varies non-monotonically with both $\mathit{Re}_{b}$ and $\mathit{Ri}_{b}$ , with its maximum (on the order of 0.2–0.3) occurring in the ‘buoyancy-dominated’ regime. Unlike $K_{{\it\rho}}$ which is very sensitive to the correct choice of $\mathscr{E}$ (i.e. $K_{{\it\rho}}\propto \mathscr{E}/(1-\mathscr{E})$ ), we show that $K_{m}$ is almost insensitive to the choice of $\mathscr{E}$ (i.e. $K_{m}\propto 1/(1-\mathscr{E})$ ) so long as $\mathscr{E}$ is not close to unity, which implies $K_{m}\approx \mathit{Ri}_{b}\mathit{Re}_{b}$ for the entire range of $\mathit{Re}_{b}$ . The turbulent Prandtl number is consequently shown to decrease monotonically with $\mathit{Re}_{b}$ and may be (to first order) simply approximated by $\mathit{Re}_{b}$ itself. Assuming $\mathit{Pr}_{t}=1$ , or $\mathit{Pr}_{t}=10$ (as is common in large-scale numerical models of the ocean general circulation), is also suggested to be a questionable assumption.

Diapycnal diffusivity, turbulent Prandtl number and mixing efficiency in Boussinesq stratified turbulence

Self-similar mixing in stratified plane Couette flow for varying Prandtl number

Understanding mixing efficiency in the oceans: Do the nonlinearities of the equation of state matter?

Understanding the effect of Prandtl number on momentum and scalar mixing rates in neutral and stably stratified flows using gradient field dynamics

Prandtl number effects on extreme mixing events in forced stratified turbulence

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

Regimes of stratified turbulence at low Prandtl number

On the Mixing Models for Stratified Flows Subjected to Concomitant Stable and Unstable Stratifications

Ocean Turbulence. Part II: Vertical Diffusivities of Momentum, Heat, Salt, Mass, and Passive Scalars

Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parameterization

Ocean Turbulence. Part I: One-Point Closure Model—Momentum and Heat Vertical Diffusivities

Turbulent transport and mixing in transitional Rayleigh-Taylor unstable flow: A priori assessment of gradient-diffusion and similarity modeling

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models, and the ocean heat engine controversy

Mixing Efficiency in the Ocean

Fluid Mixing In Stratified Gravity Currents: The Prandtl Mixing Length

Data-driven algebraic models of the turbulent Prandtl number for buoyancy-affected flow near a vertical surface

Reynolds number dependence of turbulence induced by the Richtmyer-Meshkov instability using direct numerical simulations

Scaling the Mixing Efficiency of Sediment-Stratified Turbulence

Prandtl Number Dependence of Flow Topology in Quasi-Two-dimensional Turbulent Rayleigh-Benard Convection