Onset of three-dimensionality in rapidly rotating turbulent flows

Abstract: Turbulent flows driven by a vertically invariant body force were proven to become exactly two-dimensional above a critical rotation rate, using upper bound theory. This transition in dimensionality of a turbulent flow has key consequences for the energy dissipation rate. However, its location in parameter space is not provided by the bounding procedure. To determine this precise threshold between exactly 2D and partially 3D flows, we perform a linear stability analysis over a fully turbulent 2D base state. This requires integrating numerically a quasi-2D set of equations over thousands of turnover times, to accurately average the growth rate of the 3D perturbations over the statistics of the turbulent 2Dbase flow. We leverage the capabilities of modern GPUs to achieve this task, which allows us to investigate the parameter space up to Re = 10^5. At Reynolds numbers typical of 3D DNS and laboratory experiments, Re in [10^2, 5x10^3], the turbulent 2D flow becomes unstable to 3D motion through a centrifugal-type instability. However, at even higher Reynolds number another instability takes over. A candidate mechanism for the latter instability is the parametric excitation of inertial waves by the modulated 2D flow, a phenomenon that we illustrate with an oscillatory 2D Kolmogorov flow.