2 Overview

In the original Couette experiment of Ji et al. (Reference Ji, Burin, Schartman and Goodman2006), the endcaps were split into two independently rotating rings, in order to minimise the viscous stress and ensuing Ekman circulation. Ji et al. found no evidence that qK flow broke down into turbulence at a shear Reynolds number (hereafter  $R_{s}$ ) up to $10^{6}$ , by far the best measured constraint. This was based on the absence of any increase in the Reynolds stress, a signature of turbulent flow. Another experimental approach, without split endcaps but of greater vertical extent than Ji et al., tried to isolate and study the midsection of the apparatus. Paoletti & Lathrop (Reference Paoletti and Lathrop2011) found an increase in the outer torque (presumably via a heightened Reynolds stress) for similar $R_{s}$ values. Both results could not be correct. Or could they?

Avila (Reference Avila2012) undertook direct numerical simulations of these set-ups, and found that both were prone to secondary turbulent Ekman circulations for $R_{s}\sim 10^{3}$ . This accounted for the Paoletti & Lathrop (Reference Paoletti and Lathrop2011) findings, but seemed inconsistent with the null results of Ji et al. (Reference Ji, Burin, Schartman and Goodman2006).

Edlund & Ji (Reference Edlund and Ji2014) revisited the original Ji et al. (Reference Ji, Burin, Schartman and Goodman2006) investigation in a modified experimental set-up, and found that qK flow remained stable, even when strongly perturbed. They argued that the three orders-of-magnitude difference between the $R_{s}$ of Avila’s numerical simulations ( ${\sim}10^{3}$ ) and that of the laboratory set-up ( ${\sim}10^{6}$ ) was important: more than surpassing a threshold value is involved. Edlund & Ji (Reference Edlund and Ji2015) followed with a laboratory study and model of how departures from ideal Couette flow change with  $R_{s}$ .

The recent study of Lopez & Avila (Reference Lopez and Avila2017) has succeeded in carrying out numerical simulations of the laboratory configurations of Edlund & Ji (Reference Edlund and Ji2014) for $R_{s}$ approaching $10^{5}$ , an order of magnitude below those of the laboratory, but beyond previous numerical studies of laboratory experiments. The role of viscosity is double edged. A large $R_{s}$ lowers small-scale dissipation, making it easier to sustain turbulent fluctuations. But heightened viscous stresses at smaller $R_{s}$ are a source of correlations in the radial and azimuthal velocity fluctuations, which allow the flow to extract free energy from the shear, making it easier to nourish turbulence. Lopez & Avila find that at lower $R_{s}$ values, secondary Ekman circulation forms extensively, sometimes over the entire flow region. The secondary flows are prone to their own instabilities, unrelated to any hypothetical flow instability of the underlying qK profile. But at larger $R_{s}$ , the turbulent flow recedes from the bulk of the zone of interest, confining itself to thin layers near the endcaps of the apparatus. The presumption then is that at yet larger $R_{s}$ , the turbulence would be confined to very narrow boundary layers, leaving the bulk of the interior laminar. In brief, experiments and numerical simulations now seem to be broadly consistent with an absence of any nonlinear shear flow instability in qK flow.

Edlund & Ji (Reference Edlund and Ji2014) examined three configurations of a modified Couette apparatus. In the ‘HTX’ (hydrodynamic turbulence experiment) configuration, our focus here, the end cap is split into three annuli. Two of these rings are adjacent to the inner and outer cylinders, and spin at the same rate. The middle ring is independent. An ‘optimal boundary condition’ is defined by the requirement that the torques on the inner and outer cylinder be equal and opposite, while the endplate torques vanish. This corresponds to the ideal limit of axially infinite cylinders, and is obtained for a particular ratio of inner to outer cylinder rotation rate. The simulations of Lopez & Avila (Reference Lopez and Avila2017) were carried out under this condition.

Figure 1. Recession of radial velocity isosurfaces with increasing $R_{s}$ toward axial boundaries. See text for further details.

The take-home message is summarised in figure 1. In the HTX configuration, as $R_{s}$ increases, the level of turbulent fluctuations decreases and recedes toward the endcaps. This is clear evidence that turbulence in Couette experiments is very sensitive to $R_{s}$ : low values, far from suppressing turbulence, are likely to trigger it throughout a broad region of the apparatus via secondary Ekman flows. As $R_{s}$ increases, the region of turbulence is confined to a diminishing boundary layer and the root mean square fluctuations decrease. In short, there is no compelling reason to believe that qK shear flow, in itself, is intrinsically unstable.

3 The future

The findings of Lopez & Avila (Reference Lopez and Avila2017) are broadly consistent with what is known about Ekman layers: their thickness scales inversely with a power of $R_{s}$ , in classical applications an inverse square root. We see a clear demonstration of this trend in a numerical simulation of a realistic laboratory apparatus. This work heralds an expanded role for high quality numerical investigations of laboratory Couette flow, perhaps under a variety of different boundary conditions. At present, the laboratory is approximately an order of magnitude ahead in $R_{s}$ ( $10^{6}$ versus $10^{5}$ ), and we may soon anticipate convergence.

Astrophysicists, note the focus here is reduced to the nature of Ekman circulation, not shear instability in qK flow. Whatever putative instabilities may be lurking in the bulk of the main flow (and there is scant evidence there are any), this flow is hardly in the same vein as simple shear, hypersensitive to small perturbations. Dynamicists who remain convinced that high Reynolds number shear flow of any sort is unstable must be prepared to dig below the level of the unavoidable secondary Ekman circulation. (Ever more powerful numerical simulations of free-boundary discs are an alternative route.)

Other dynamics, notably baroclinic rotation, may be a source of flow destabilisation, but the goal here is not instability per se. It is not even the onset of turbulence. It is a mode of enhanced internal angular momentum transport. This is of course what turbulence sometimes does – but not always. The key is an enhanced correlation in the radial and azimuthal velocity fluctuations. If the cause of the breakdown of laminar flow is the radial shear itself, this correlation is more likely to appear (as in the magneto-rotational instability), than if the cause is, say, an adverse entropy gradient or related vertical velocity gradient.

Fluid dynamicists, note there is a grand challenge problem here. All indicators seem to imply that large Reynolds number qK flow is truly nonlinearly stable. Can this be proven analytically? What valuable insights might emerge by even a partially successful attempt to do so?