In this paper we investigate the effect of stable stratification on plane Couette flow when gravity is oriented in the spanwise direction. When the flow is turbulent we observe near-wall layering and associated new mean flows in the form of large-scale spanwise-flattened streamwise rolls. The layers exhibit the expected buoyancy scaling $l_{z}\sim U/N$ where $U$ is a typical horizontal velocity scale and $N$ the buoyancy frequency. We associate the new coherent structures with a stratified modification of the well-known large-scale secondary circulation in plane Couette flow. We find that the possibility of the transition to sustained turbulence is controlled by the relative size of this buoyancy scale to the spanwise spacing of the streaks. In parts of parameter space transition can also be initiated by a newly discovered linear instability in this system (Facchini et al., J. Fluid Mech., vol. 853, 2018, pp. 205–234). When wall turbulence can be sustained the linear instability opens up new routes in phase space for infinitesimal disturbances to initiate turbulence. When the buoyancy scale suppresses turbulence the linear instability leads to more ordered nonlinear states, with the possibility for intermittent bursts of secondary shear instability.

Recent experimental observations (Kühnen et al., Nat. Phys., vol. 14, 2018b, pp. 386–390) have shown that flattening a turbulent streamwise velocity profile in pipe flow destabilises the turbulence so that the flow relaminarises. We show that a similar phenomenon exists for laminar pipe flow profiles in the sense that the nonlinear stability of the laminar state is enhanced as the profile becomes more flattened. The flattening of the laminar base profile is produced by an artificial localised body force designed to mimic an obstacle used in the experiments of Kühnen et al. (Flow Turbul. Combust., vol. 100, 2018a, pp. 919–943) and the nonlinear stability measured by the size of the energy of the initial perturbations needed to trigger transition. Significant drag reduction is also observed for the turbulent flow when triggered by sufficiently large disturbances. In order to make the nonlinear stability computations more efficient, we examine how indicative the minimal seed – the disturbance of smallest energy for transition – is in measuring transition thresholds. We first show that the minimal seed is relatively robust to base profile changes and spectral filtering. We then compare the (unforced) transition behaviour of the minimal seed with several forms of randomised initial conditions in the range of Reynolds numbers $Re=2400$ – $10\,000$ and find that the energy of the minimal seed after the Orr and oblique phases of its evolution is close to that of a critical localised random disturbance. In this sense, the minimal seed at the end of the oblique phase can be regarded as a good proxy for typical disturbances (here taken to be the localised random ones) and is thus used as initial condition in the simulations with the body force. The enhanced nonlinear stability and drag reduction predicted in the present study are an encouraging first step in modelling the experiments of Kühnen et al. and should motivate future developments to fully exploit the benefits of this promising direction for flow control.

We consider turbulence in a stratified ‘Kolmogorov’ flow, driven by horizontal shear in the form of sinusoidal body forcing in the presence of an imposed background linear stable stratification in the third direction. This flow configuration allows the controlled investigation of the formation of coherent structures, which here organise the flow into horizontal layers by inclining the background shear as the strength of the stratification is increased. By numerically converging exact steady states from direct numerical simulations of chaotic flow, we show, for the first time, a robust connection between linear theory predicting instabilities from infinitesimal perturbations to the robust finite-amplitude nonlinear layered state observed in the turbulence. We investigate how the observed vertical length scales are related to the primary linear instabilities and compare to previously considered examples of shear instability leading to layer formation in other horizontally sheared flows.

Efforts to model accretion disks in the laboratory using Taylor–Couette flow apparatus are plagued with problems due to the substantial impact the end plates have on the flow. We explore the possibility of mitigating the influence of these end plates by imposing stable stratification in their vicinity. Numerical computations and experiments confirm the effectiveness of this strategy for restoring the axially homogeneous quasi-Keplerian solution in the unstratified equatorial part of the flow for sufficiently strong stratification and moderate layer thickness. If the rotation ratio is too large, however (e.g. ${\it\Omega}_{o}/{\it\Omega}_{i}=(r_{i}/r_{o})^{3/2}$ , where ${\it\Omega}_{o}/{\it\Omega}_{i}$ is the angular velocity at the outer/inner boundary and $r_{i}/r_{o}$ is the inner/outer radius), the presence of stratification can make the quasi-Keplerian flow susceptible to the stratorotational instability. Otherwise (e.g. for ${\it\Omega}_{o}/{\it\Omega}_{i}=(r_{i}/r_{o})^{1/2}$ ), our control strategy is successful in reinstating a linearly stable quasi-Keplerian flow away from the end plates. Experiments probing the nonlinear stability of this flow show only decay of initial finite-amplitude disturbances at a Reynolds number $Re=O(10^{4})$ . This observation is consistent with most recent computational (Ostilla-Mónico, et al. J. Fluid Mech., vol. 748, 2014, R3) and experimental results (Edlund & Ji, Phys. Rev. E, vol. 89, 2014, 021004) at high $Re$ , and reinforces the growing consensus that turbulence in cold accretion disks must rely on additional physics beyond that of incompressible hydrodynamics.

We consider long-time simulations of two-dimensional turbulence body forced by $\sin 4y\hat {\boldsymbol{x}} $ on the torus $(x, y)\in \mathop{[0, 2\mathrm{\pi} ] }\nolimits ^{2} $ with the purpose of extracting simple invariant sets or ‘exact recurrent flows’ embedded in this turbulence. Each recurrent flow represents a sustained closed cycle of dynamical processes which underpins the turbulence. These are used to reconstruct the turbulence statistics using periodic orbit theory. The approach is found to be reasonably successful at a low value of the forcing where the flow is close to but not fully in its asymptotic (strongly) turbulent regime. Here, a total of 50 recurrent flows are found with the majority buried in the part of phase space most populated by the turbulence giving rise to a good reproduction of the energy and dissipation p.d.f. However, at higher forcing amplitudes now in the asymptotic turbulent regime, the generated turbulence data set proves insufficiently long to yield enough recurrent flows to make viable predictions. Despite this, the general approach seems promising providing enough simulation data is available since it is open to extensive automation and naturally generates dynamically important exact solutions for the flow.

We propose a general strategy for determining the minimal finite amplitude disturbance that triggers transition to turbulence in shear flows. This involves constructing a variational problem that searches over all disturbances of fixed initial amplitude which respect the boundary conditions, incompressibility and the Navier–Stokes equations, to maximize a chosen functional over an asymptotically long time period. The functional must be selected such that it identifies turbulent velocity fields by taking significantly enhanced values compared to those for laminar fields. We illustrate this approach using the ratio of the final to initial perturbation kinetic energies (energy growth) as the functional and the energy norm to measure amplitudes in the context of pipe flow. Our results indicate that the variational problem yields a smooth converged solution provided that the initial amplitude is below the threshold for transition. This optimal is the nonlinear analogue of the well-studied (linear) transient growth optimal. At the critical threshold, the optimization seeks out a disturbance that is on the ‘edge’ of turbulence during the period. Above this threshold, when disturbances trigger turbulence by the end of the period, convergence is then practically impossible. The first disturbance found to trigger turbulence as the amplitude is increased identifies the ‘minimal seed’ for the given geometry and forcing (Reynolds number). We conjecture that it may be possible to select a functional such that the converged optimal below threshold smoothly converges to the minimal seed at threshold. Our choice of the energy growth functional is shown to come close to this for the pipe flow geometry investigated here.

Fully three-dimensional computations of flow through a long pipe demand a huge number of degrees of freedom, making it very expensive to explore parameter space and difficult to isolate the structure of the underlying dynamics. We therefore introduce a ‘2+ε-dimensional’ model of pipe flow, which is a minimal three-dimensionalization of the axisymmetric case: only sinusoidal variation in azimuth plus azimuthal shifts are retained; yet the same dynamics familiar from experiments are found. In particular the model retains the subcritical dynamics of fully resolved pipe flow, capturing realistic localized ‘puff-like’ structures which can decay abruptly after long times, as well as global ‘slug’ turbulence. Relaminarization statistics of puffs reproduce the memoryless feature of pipe flow and indicate the existence of a Reynolds number about which lifetimes diverge rapidly, provided that the pipe is sufficiently long. Exponential divergence of the lifetime is prevalent in shorter periodic domains. In a short pipe, exact travelling-wave solutions are found near flow trajectories on the boundary between laminar and turbulent flow. In a long pipe, the attracting state on the laminar–turbulent boundary is a localized structure which resembles a smoothened puff. This ‘edge’ state remains localized even for Reynolds numbers at which the turbulent state is global.