Comparative analysis of four second-moment turbulence closure models for the oceanic mixed layer

In this comparative study, four different algebraic second-moment turbulenc e closure models are investigated in detail. These closure schemes differ i n the number of terms considered for the closure of the pressure-strain cor relations. These four turbulence closures result in the eddy-diffusivity pr inciple such that the closure assumptions are contained in dimensionless so -called stability functions. Their performance in terms of Prandtl number, Monin-Obukhov similarity theory, and length scale ratios are first tested a gainst data for simple flows. The turbulence closure is then completed by m eans of a k-epsilon two-equation model, but other models such as the two-eq uation model by Mellor and Yamada could also be used. The concept of the st eady-state Richardson number for homogeneous shear layers is exploited for calibrating the sensitivity of the four models to shear and stable stratifi cation. Idealized simulations of mixed layer entrainment into stably strati fied flow due to surface stress and due to free convection are carried out. For the latter experiment, comparison to recent large eddy simulation data is made. Finally, the well-known temperature profile data at OWS Papa are simulated for an annual cycle. The main result of this paper is that the ov erall performance of the new second-moment closure model by Canuto et al.-e xpressed as nondimensional stability functions-is superior compared to the others in terms of physical soundness, predictability, computational econom y, and numerical robustness.