Turbulent Transport by Diffusive Stratified Shear Flows: From Local to Global Models. II. Limitations of Local Models

Abstract

This paper continues the systematic investigation of diffusive shear instabilities initiated in Part I of this series. In this work, we primarily focus on quantifying the impact of nonlocal mixing, which is not taken into account in Zahn’s mixing model. We present the results of direct numerical simulations in a new model setup designed to contain coexisting laminar and turbulent shear layers. As in Part I, we use the low Péclet number approximation of Lignières to model the evolution of the perturbations. Our main findings are twofold. First, turbulence is not necessarily generated whenever Zahn’s nonlinear criterion J Pr < ( J Pr) c can also take place in a manner analogous to convective overshoot. We found that the turbulent kinetic energy decays roughly exponentially beyond the edge of the shear-unstable region, on a lengthscale δ that is directly proportional to the scale of the turbulent eddies, which are themselves of the order of the Zahn scale (see Part I). Our results suggest that mixing by diffusive shear instabilities should be modeled with more care than is currently standard in stellar evolution codes.