---
abstract: 'The transition from Tayor vortex flow to wavy-vortex flow is revisited. The Self-Sustaining Process (SSP) of Waleffe \[Phys. Fluids [**9**]{}, 883–900 (1997)\] proposes that a key ingredient in transition to turbulence in wall-bounded shear flows is a three-step process involving rolls advecting streamwise velocity, leading to streaks which become unstable to a wavy perturbation whose nonlinear interaction with itself feeds the rolls. We investigate this process in Taylor-Couette flow. The instability of Taylor-vortex flow to wavy-vortex flow, a process which is the inspiration for the second phase of the SSP, is shown to be caused by the streaks, with the rolls playing a negligible role, as predicted by Jones \[J. Fluid Mech. [**157**]{}, 135–162 (1985)\] and demonstrated by Martinand et al. \[Phys. Fluids [**26**]{}, 094102 (2014)\]. In the third phase of the SSP, the nonlinear interaction of the waves with themselves reinforces the rolls. We show this both quantitatively and qualitatively, identifying physical regions in which this reinforcement is strongest, and also demonstrate that this nonlinear interaction depletes the streaks.'
author:
- Tommy Dessup
- 'Laurette S. Tuckerman'
- José Eduardo Wesfreid
- Dwight Barkley
- 'Ashley P. Willis'
title: 'The Self-Sustaining Process in Taylor-Couette Flow'
---

[**Keywords:**]{} Taylor-Couette instability, transition to turbulence .5cm

Introduction {#sec:Introduction}
============

The first successful linear stability analysis for a viscous fluid was carried out in 1923 by G.I. Taylor [@taylor1923stability] for the flow between two concentric differentially rotating cylinders. What then became known as Taylor-Couette flow has played a central role in hydrodynamic stability theory ever since. In the standard configuration of a stationary outer cylinder, as the inner cylinder rotation rate is increased, laminar flow is succeeded by axisymmetric Taylor vortices via the centrifugal instability first explained by Rayleigh [@rayleigh1916]. The Taylor vortices subsequently develop azimuthal waves, seen in experiments by researchers such as Coles [@coles1965transition], Swinney and co-workers [@gorman1982spatial; @king1984wave; @andereck1986flow] and others [@hegseth1996bifurcations; @wereley1998spatio]. Wavy-vortex flow was studied computationally when this became possible in the 1980s by authors such as Jones [@jones1981nonlinear; @jones1985transition], Marcus [@marcus1984a; @marcus1984b] and others [@edwards1991onset; @antonijoan2002stable].

Taylor-Couette flow has also been studied as a way of approaching plane Couette flow, which undergoes transition to three-dimensional turbulence despite being linearly stable at all Reynolds numbers. The azimuthal, radial, and axial directions of Taylor-Couette flow play the role of the streamwise, cross-channel, and spanwise direction, respectively. As the ratio between the cylinder radii approaches one, the correspondence between the two flows becomes exact. The possibility of approaching plane Couette flow via Taylor-Couette flow has been used by many authors for many different purposes. Nagata [@nagata1990three] used homotopy to calculate otherwise inaccessible unstable steady states of plane Couette flow. Hristova et al. [@hristova2002transient] and Meseguer et al. [@meseguer2002energy] compared transient growth rates between the two flows. Prigent et al. [@prigent2002large] extended the observation of coexisting turbulent and laminar regions seen in Taylor-Couette by Coles [@coles1965transition] to plane Couette flow. Faisst and Eckhardt [@faisst2000transition] used Taylor-Couette flow to approach the turbulent lifetimes and intermittency of plane Couette flow. A very narrow gap Taylor-Couette geometry was used as a proxy for plane Couette flow by Shi et al. [@shi2013scale] to calculate the statistical threshold of sustained turbulence and by Lemoult et al. [@lemoult2016directed] to establish that this transition was manifested as a directed percolation phase transition.

Here, we take the analogy in the opposite direction: extending an idea developed for plane Couette flow to Taylor-Couette flow. Waleffe [@waleffe1995hydrodynamic; @hamilton1995regeneration; @waleffe1997self] has proposed a now widely-accepted three-part mechanism, by which streamwise rolls (damped in the plane Couette case), cause streamwise streaks (by simple advection of the streamwise velocity contours), which become wavy (through instability), acquiring streamwise dependence. The nonlinear self-interaction of the wavy streaks drives the streamwise rolls, thus closing the cycle. The mechanism is similar to that proposed by Hall and co-workers [@hall1991strongly; @hall2010streamwise; @blackburn2013lower] and by Beaume and co-workers [@beaume2014exact; @beaume2015reduced; @beaume2016modulated]. Experimental evidence for the SSP in plane boundary layer and channel flow has been reported by Wesfreid and colleagues in [@duriez2009self; @klotz2017experiments]. These experiments show a strong correlation between the growth of rolls and the presence of waves: both phenomena occur above the same Reynolds-number threshold.

Although the SSP was influenced by these phenomena in Taylor-Couette flow, it has not actually been applied to Taylor-Couette flow itself. The main purpose of this paper is to see how the Self-Sustaining Process plays out in Taylor-Couette flow, where the analogous structures, i.e. axisymmetric and wavy Taylor vortices, are actually stable equilibrium states.

Equations, Methods and Parameters
=================================

The equations governing Taylor-Couette flow and the methods for computing it are sufficiently well known as to warrant only a very brief exposition. The inner and outer cylinders have radii and angular velocities $R_j$ and $\Omega_j$. From these, along with the kinematic viscosity $\nu$, can be constructed the length scale ${d}\equiv R_2-R_1$, the time scale ${d}^2/\nu$, the two Reynolds numbers $Re_j\equiv R_j\Omega_j{d}/\nu$, and the radius ratio $\eta\equiv R_1/R_2$. The non-dimensionalized governing equations and boundary conditions are then

$$\begin{aligned}
{\partial}_t {\mathbf{U}}&= {\mathbf{U}}\times{\nabla\times}{\mathbf{U}}-\nabla P + {\nabla^2}{\mathbf{U}}\\
  {\nabla\cdot}{\mathbf{U}}&= 0 \\
  {\mathbf{U}}&= Re_j {\mathbf{e}}_\theta\mbox{~at~} r=r_j \equiv R_j/d,\,\, j=1,2\end{aligned}$$

\[eq:NS\]

We will restrict our consideration to the classic inner-cylinder-rotation case with $\Omega_2=0$ so that $Re_2=0$ and hence we use $Re$ to denote the inner Reynolds number $Re_1$. Nonlinear Taylor-vortex and wavy-vortex flows, denoted by TVF and WVF or ${\mathbf{U}}_{\rm TVF}$ and ${\mathbf{U}}_{\rm WVF}$, are calculated by solving the evolution equations numerically. For linear stability analysis, the nonlinear code has been adapted to solve the linearized equations

$$\begin{aligned}
{\partial}_t {\mathbf{u}}&= {\mathbf{U}}\times{\nabla\times}{\mathbf{u}}+ {\mathbf{u}}\times{\nabla\times}{\mathbf{U}}-\nabla p + {\nabla^2}{\mathbf{u}}\\
  {\nabla\cdot}{\mathbf{u}}&= 0 \\
  {\mathbf{u}}&= 0 \quad\mbox{~at~} r=r_j,\,\, j=1,2\end{aligned}$$

\[eq:linstab\]

where ${\mathbf{U}}$ is the flow whose stability is sought. Temporal integration of effectively carries out the power method, converging to the eigenvector whose eigenvalue has largest real part. Most commonly, we take ${\mathbf{U}}$ to be Taylor-vortex flow, ${\mathbf{U}}_{\rm TVF}$, and the power method returns the wavy vortex eigenvector ${\mathbf{u}}_{\rm wvf}$ and corresponding eigenvalue.

The code we use represents functions on a spatial Chebyshev grid in the radial direction $r$ and on equally spaced points in the azimuthal $\theta$ and axial $z$ directions, with spatial derivatives taken via finite differences in $r$ and by differentiation of Fourier series in $\theta$, $z$. Multiplications are carried out in the grid space representation by Fourier transforming in $\theta,z$. Taylor-vortex flow is calculated in an axisymmetric domain with $N_r=33$ radial points and $N_z=16$ points over the axial domain $[0,{L}_z]$ or, equivalently, multiples of the wavenumber $2\pi/{L}_z$. Computations of wavy-vortex flow eigenvectors use a single azimuthal mode ${M_0}$. Nonlinear wavy-vortex flow is calculated using $N_\theta=16$ points in the azimuthal sector $[0,2\pi/{M_0}]$ or, equivalently, multiples of the wavenumber ${M_0}$.

One difficulty is deciding which of the many TVFs or WVFs to study. Each TVF is characterized by an axial wavenumber, and each WVF has an axial and an azimuthal wavenumber. States with different wavenumbers can be simultaneously stable, as emphasized by Coles [@coles1965transition] and by many subsequent researchers [@gorman1982spatial; @king1984wave; @hegseth1996bifurcations]. Jones [@jones1985transition] and Antonijoan & Sanchez [@antonijoan2002stable] have shown the complexity of the bifurcations and ranges of existence of wavy-vortex states with different azimuthal wavenumbers as the radius ratio and the axial wavelength are varied. We select the radius ratio to be $\eta=0.92$, corresponding to $r_1=11.5$ and $r_2=12.5$. To make connection with the SSP in plane Couette flow, we take the axial wavelength to be ${L}_z=2$, corresponding to a spanwise wavelength of 4 half-gaps, near the length considered by Waleffe [@waleffe1995hydrodynamic; @hamilton1995regeneration; @waleffe1997self]. (Note that the length scale in the Taylor-Couette problem is the full gap.) We use the term circumferential wavelength to denote a length at the midgap $r=\bar{r}$, in contrast with an azimuthal wavelength, which is expressed in radians and necessarily a fraction of $2\pi$. To approximate the streamwise wavelength of 10 half-gaps studied by Waleffe, we first express the circumferential wavelength $L_\theta$ of a wavy-vortex state with azimuthal wavenumber ${M_0}$ in units of the gap $${L}_\theta = \frac{2\pi \bar{r}}{{M_0}}
  =\frac{2\pi}{{M_0}}\;\frac{({r_2}+{r_1})/2}{{r_2}-{r_1}}
  =\frac{\pi}{{M_0}}\;\frac{1+\frac{{r_1}}{{r_2}}}{1-\frac{{r_1}}{{r_2}}}
  =\frac{\pi}{{M_0}}\;\frac{1+\eta}{1-\eta}$$ Setting $\eta=0.92$ and ${L}_\theta=5$, corresponding to 10 in half-gaps, leads to $${M_0}=\frac{\pi}{{L}_\theta}\;\frac{1+\eta}{1-\eta}
 =\frac{\pi}{5}\;\frac{1.92}{0.08}\approx 15$$ The critical Reynolds number for onset of Taylor-vortex flow in which only the inner cylinder rotates is approximately $$Re_{\rm TVF}\approx\sqrt{\frac{1708}{\eta(1-\eta)}}\frac{1+\eta}{2}$$ which diverges as the narrow-gap (or plane Couette) limit $\eta\rightarrow 1$ is approached [@faisst2000transition]. The dependence of the critical Reynolds on $\eta$ is shown in Fig. \[fig:detpar\](a), together with the relationship between ${M_0}$ and $\eta$ for ${L}_\theta=5$. For $\eta=0.92$, Taylor vortices appear above $Re_{\rm TVF}\approx 146$. For these values of $\eta$, ${L}_z$, and ${M_0}$, Taylor-vortex flow remains stable until $Re_{\rm WVF} \approx 201$, above which the flow becomes unstable to wavy Taylor vortices. Figure \[fig:vis1\](a) shows the Taylor-vortex flow while Fig. \[fig:vis2\] shows the wavy-vortex flow, both computed at $Re=300$.

![(a) Dependence of critical Reynolds number $Re^{\rm TVF}$ on the radius ratio $\eta$ for transition to Taylor-vortex flow (right scale). Also shown is the relationship between azimuthal wavenumber ${M_0}$ and $\eta$ for circumferential wavelength $L_\theta=5$ (left scale). Our chosen parameter values are $\eta=0.92$ and ${M_0}=15$. (b) Schematic decomposition of flow into ${\mathbf{U}}_{\rm Cou}$, ${\mathbf{U}}_{\rm mean}$, ${\mathbf{U}}_{\rm roll}$, ${\mathbf{U}}_{\rm streak}$, ${\mathbf{U}}_{\rm wave}$ according to axial and azimuthal Fourier modes $k$ and $m$.[]{data-label="fig:detpar"}](figure1){width="0.85\columnwidth"}

Analysis in terms of Self-Sustaining Process
============================================

We begin our analysis by introducing notation. Flow fields ${\mathbf{U}}$ can be decomposed as follows; see Fig. \[fig:detpar\](b).

$$\begin{aligned}
  {\mathbf{U}}&=& \sum_k\sum_m \left(\hat{U}^{k,m}_r(r) {\mathbf{e}}_r+ \hat{U}^{k,m}_\theta(r) {\mathbf{e}}_\theta + \hat{U}^{k,m}_z(r) {\mathbf{e}}_z\right) e^{i(kz/{L}_z+m{M_0}\theta)}\qquad
  \label{eq:fourdecomp}\\
  &=& {\mathbf{U}}_{\rm Cou}+{\mathbf{U}}_{\rm mean} + {\mathbf{U}}_{\rm roll} + {\mathbf{U}}_{\rm streak} + {\mathbf{U}}_{\rm wave}\end{aligned}$$

\[eq:decomp1\]

where

$$\begin{aligned}
{\mathbf{U}}_{\rm Cou} &\equiv \left(A r + \frac{B}{r}\right){\mathbf{e}}_\theta \label{eq:Cou}\\
{\mathbf{U}}_{\rm mean} &\equiv \hat{U}^{0,0}_\theta(r) {\mathbf{e}}_\theta-{\mathbf{U}}_{\rm Cou}\label{eq:mean}\\
{\mathbf{U}}_{\rm roll} &\equiv \sum_{k\neq 0}
\left(\hat{U}_r^{k,0}(r){\mathbf{e}}_r + \hat{U}_z^{k,0}(r){\mathbf{e}}_z\right)\:e^{ikz/{L}_z}\\
  {\mathbf{U}}_{\rm streak} &\equiv \sum_{k\neq 0}
  \hat{U}_\theta^{k,0}(r){\mathbf{e}}_\theta \:e^{ikz/{L}_z}\\
  {\mathbf{U}}_{\rm wave} &\equiv
\sum_k\sum_{m\neq 0}
  \hat{{\mathbf{U}}}^{km}(r)
  e^{i(kz/{L}_z+m{M_0}\theta)}
  \label{eq:Ewave}\end{aligned}$$

\[eq:decomp2\]

Note that defines ${\mathbf{U}}_{\rm mean}$ to be the $(\theta,z)$-independent deviation from laminar Couette flow ${\mathbf{U}}_{\rm Cou}$, in contrast to Waleffe [@waleffe1995hydrodynamic; @hamilton1995regeneration; @waleffe1997self] whose mode ${\rm M}$ includes the laminar Couette solution . In terms of this decomposition, Taylor-vortex flow and wavy-vortex flow take the form

$$\begin{aligned}
  {\mathbf{U}}_{\rm TVF}&={\mathbf{U}}_{\rm Cou}+{\mathbf{U}}_{\rm mean}+{\mathbf{U}}_{\rm roll}+{\mathbf{U}}_{\rm
    streak} \label{eq:TVF} 
  \\
  {\mathbf{U}}_{\rm WVF}&={\mathbf{U}}_{\rm Cou}+{\mathbf{U}}_{\rm mean}+{\mathbf{U}}_{\rm roll}+{\mathbf{U}}_{\rm
    streak}+{\mathbf{U}}_{\rm wave}  \end{aligned}$$

\[eq:decomp3\]

Waleffe’s Self-Sustaining Process [@waleffe1995hydrodynamic; @hamilton1995regeneration; @waleffe1997self] describes three steps involving the components ${\mathbf{U}}_{\rm roll}$, ${\mathbf{U}}_{\rm streak}$, and ${\mathbf{U}}_{\rm wave}$:\
A) ${\mathbf{U}}_{\rm roll} \Longrightarrow {\mathbf{U}}_{\rm streak}$. This is a statement of kinematic advection of the azimuthal velocity.\
B) ${\mathbf{U}}_{\rm streak} \Longrightarrow {\mathbf{U}}_{\rm wave}$. This is described by Waleffe as a linear instability.\
C) ${\mathbf{U}}_{\rm wave}\Longrightarrow {\mathbf{U}}_{\rm roll}$. The nonlinear interaction of the wave with itself reinforces the rolls.

Rolls to streaks
----------------

The SSP begins with streamwise invariant rolls ${\mathbf{U}}_{\rm roll}$ and considers the development of streaks from these rolls. Rolls transport fluid with high azimuthal velocity from the inner cylinder towards the outer cylinder and vice versa, causing the azimuthal velocity profile to vary along $z$ with the axial periodicity of the rolls. In plane Couette flow, or Waleffe’s free-slip version [@waleffe1997self] now sometimes called Waleffe flow [@beaume2014exact; @beaume2015reduced; @beaume2016modulated; @chantry2016turbulent; @chantry2017universal], rolls are not themselves an equilibrium state. Hence in the planar case it is necessary to initiate the SSP by inserting rolls into the flow and observing the resulting streak development. Permanent rolls and streaks have been produced in variants of plane Couette flow by including a spanwise-oriented wire or ribbon experimentally [@dauchot1995finite; @bottin1997intermittency; @bottin1998experimental] or numerically [@barkley1999stability]. For the Taylor-Couette problem, however, this phase is straightforward. The rolls and the streaks that they generate are contained in Taylor-vortex flow, which bifurcates supercritically and exists as a stable nonlinear equilibrium. In Fig. \[fig:vis1\](a), calculated at $Re=300$, the rolls are the meridional-plane flow indicated by arrows. The streaks are the axial variation in the azimuthal flow driven by the rolls and are seen as the colored patches.

![Visualizations in the meridional $(r,z)$ plane of (a) Taylor-vortex flow (without laminar Couette flow), (b) the ${M_0}=15$ eigenvector leading to wavy-vortex flow, and (c) nonlinear interaction of this eigenvector with itself. The parameters are $Re=300$ and $\eta=0.92$. The inner cylinder is on the left and the outer cylinder on the right. In each case, the meridional velocity within the plane is indicated by arrows and the azimuthal velocity perpendicular to it is indicated by colors. Red indicates a positive deviation of the azimuthal velocity from laminar Couette flow, blue a negative deviation, and green no deviation. Thus, in (a) the arrows show the rolls and the colors show the streaks of Taylor-vortex flow. The white dashed boxes in (c) and (a) highlight the alignment between the axial components (arrows) of $\langle{\mathbf{u}}_{\rm
    wvf}\times\nabla\times{\mathbf{u}}_{\rm wvf}\rangle$ and of the rolls of ${\mathbf{U}}_{\rm TVF}$, which comprise the third step of the SSP.[]{data-label="fig:vis1"}](figure2){width="0.75\columnwidth"}

Streaks to waves
----------------

![Wavy-vortex flow (including Taylor-vortex flow but not laminar Couette flow) at $Re=300$. Above: four meridional planes over azimuthal period $[0,2\pi/15]$. Azimuthal velocity indicated by colors, meridional velocity by arrows. Below: one azimuthal period $[0,2\pi/15]$ at mid-gap $\bar{r}=12$. Radial velocity indicated by colors. The dashed lines indicate the positions of the four meridional planes shown above. []{data-label="fig:vis2"}](figure3){width="0.75\columnwidth"}

![Comparison of linearization about ${\mathbf{U}}_{\text{TVF}}$ and about ${\mathbf{U}}_{\text{TVF}}-{\mathbf{U}}_{\text{roll}}$. Eigenvectors ${\mathbf{u}}_{\text{wvf}}$ resulting from linearization about (a) only ${\mathbf{U}}_{\text{TVF}}-{\mathbf{U}}_{\text{roll}}$ and (b) the full ${\mathbf{U}}_{\text{TVF}}$ for $Re=300$. Azimuthal velocity designated by color (red for positive, blue for negative, green for zero) and radial and axial velocity by arrows. (c) Growth rate (real part of eigenvalue) for linearization about ${\mathbf{U}}_{\text{TVF}}$ (black, solid), ${\mathbf{U}}_{\text{TVF}}-{\mathbf{U}}_{\text{roll}}$ (blue, dashed), as a function of $Re$. Since omitting ${\mathbf{U}}_{\rm roll}$ from the base flow barely changes the eigenvector or eigenvalue, it is clear that it plays no role in the instability.[]{data-label="fig:eigcompare"}](figure4){width="0.85\columnwidth"}

We now turn to the second stage of the SSP in which the streaks become unstable to waviness. Once again, the situation in the Tayor-Couette problem is much more clear cut than in the planar case. The onset of waviness is a distinct supercritical instability – the transition from Taylor-vortex flow ${\mathbf{U}}_{\rm TVF}$ to wavy-vortex flow ${\mathbf{U}}_{\rm WVF}$. In the ${\mathbf{U}}_{\rm WVF}$ state shown in Fig. \[fig:vis2\], the flow has azimuthal variation (waviness) and is an azimuthally travelling wave. In 1985, Jones [@jones1985transition] suggested that the instability arose from the streaks, i.e. the axial variation of the azimuthal flow, which he called azimuthal jets. Thirty years later, Martinand, Serre and Lueptow [@martinand2014mechanisms] confirmed this idea by constructing the linear operator governing the wavy instability and showing that the eigenvalues of the portion of the operator arising from the azimuthal shear, i.e. the streaks, best matched the eigenvalues of the entire operator. They also demonstrated a number of common features between the transition to wavy vortex flow and the Kelvin-Helmholtz instability, notably a phase speed intermediate between that of the two cylinders and the multiplicity of possible azimuthal wavenumbers.

We show this by a different procedure, carrying out linearization about ${\mathbf{U}}_{\rm TVF}$ and about ${\mathbf{U}}_{\rm TVF}-{\mathbf{U}}_{\rm roll}$, i.e. the Taylor vortex flow without its radial or axial components; see Eqs. . Fig. \[fig:eigcompare\] compares the eigenvectors and growth rates resulting from these two linearizations. Since omitting ${\mathbf{U}}_{\rm roll}$ from the base flow barely changes the eigenvector or eigenvalue, it is clear that it plays no role in the instability. In contrast, linearization about ${\mathbf{U}}_{\rm
  TVF}-{\mathbf{U}}_{\rm streak}$, i.e. omitting the axial dependence of the azimuthal flow, leads to eigenvalues with very small growth rate and eigenvectors with no resemblance to those of ${\mathbf{U}}_{\rm TVF}$. (These results are not displayed.) These numerical experiments confirm that the instability mechanism responsible for the transition of ${\mathbf{U}}_{\rm TVF}$ to ${\mathbf{U}}_{\rm WVF}$ is the axial variation of the azimuthal velocity.

In addition to linearization, we examine the energy content in the flow components of the nonlinear states. We decompose both Taylor-vortex flow ${\mathbf{U}}_{\rm TVF}$ and the wavy-vortex flow ${\mathbf{U}}_{\rm WVF}$ into components given in Eqs. - and compute the energy of each. Fig. \[fig:bifdiag\](a) shows the variation of the energy components as a function of Reynolds number. (The much larger energy of ${\mathbf{U}}_{\rm Cou}$ and a contribution combining ${\bf U}_{\rm Cou}$ and ${\bf U}_{\rm mean}$ are not shown.) ${\mathbf{U}}_{\rm TVF}$ appears at $Re=146$ and ${\mathbf{U}}_{\rm WVF}$ appears at $Re=201$. It can be seen that ${\rm E}^{\rm WVF}_{\rm streak}$, the energy of the streaks in ${{\mathbf{U}}}_{\rm WVF}$, is substantially decreased from the analogous quantity ${\rm
  E}^{\rm TVF}_{\rm streak}$ in ${\mathbf{U}}_{\rm TVF}$. This decrease is almost exactly counterbalanced by the energy in the waviness, ${\rm E}^{\rm WVF}_{\rm wave}$, suggesting that the energy in the waviness is extracted from the streaks. The energy in the rolls is small and is almost the same in the two states. Thus, in addition to the linear instability mechanism, the comparison between the energy content of the saturated nonlinear states with and without waves shows that streaks feed the waves. As stated by Waleffe [@waleffe1997self], it is not the rolls but the streaks whose energy is drained by the waves.

![(a) Energy decomposition for ${\mathbf{U}}_{\rm TVF}$ and ${\mathbf{U}}_{\rm WVF}$. (See Fig. \[fig:detpar\](b) for definitions of this decomposition.) Curves marked with crosses correspond to the components of ${\mathbf{U}}_{\rm TVF}$ originating at $Re=146$. Curves marked with circles correspond to the energy components of ${\mathbf{U}}_{\rm WVF}$ which bifurcates at $Re=201$. The streak energy is lower for WVF than it is for TVF; the difference between the two is close to the energy in the waves (which is necessarily zero for TVF). The energy in the deviation of the mean from Couette flow is also lower for WVF than for TVF. The energy in the rolls is approximately the same for the two flows. (b) Normalized inner product of nonlinear self-interaction $\langle {\bf NL},{{\mathbf{U}}}_{\rm roll}\rangle/||{{\mathbf{U}}}_{\rm roll}||$, $\langle {\bf NL},{{\mathbf{U}}}_{\rm streak}\rangle/||{{\mathbf{U}}}_{\rm streak}||$ and $\langle {\bf NL},{{\mathbf{U}}}_{\rm mean}\rangle/||{{\mathbf{U}}}_{\rm mean}||$ for rolls, streaks, and deviation of the mean from Couette flow. The nonlinear term ${\bf NL}$ feeds the rolls and mean but drains the streaks.[]{data-label="fig:bifdiag"}](figure5){width="0.85\columnwidth"}

Waves to rolls
--------------

The key novelty of the SSP is the positive feedback of the waviness on the rolls. To study this in Taylor-Couette flow, we calculate the eigenvector ${\mathbf{u}}_{\rm wvf}$ responsible for the bifurcation to wavy vortices, shown in Fig. \[fig:vis1\](b). (This complex eigenvector is shown here at one spatial or temporal phase.) We then compute the nonlinear interaction of ${\mathbf{u}}_{\rm wvf}$ with itself, in the form ${\mathbf{u}}_{\rm wvf} \times \nabla \times
{\mathbf{u}}_{\rm wvf}$. Since ${\mathbf{u}}_{\rm wvf} \sim e^{\pm i {M_0}\theta}$, this quadratic term leads to azimuthal dependence of the form $e^{\pm 2 i {M_0}\theta}$ (second harmonic) and $1$ (constant). We are interested in the constant contribution which has the form: $${\bf NL} \equiv \langle {\mathbf{u}}_{\rm wvf} \times \nabla \times {\mathbf{u}}_{\rm wvf}\rangle
  \equiv {\mathbf{u}}_{\rm wvf}^R \times \nabla \times {\mathbf{u}}_{\rm wvf}^R
  + {\mathbf{u}}_{\rm wvf}^I \times \nabla \times {\mathbf{u}}_{\rm wvf}^I
  \label{eq:quad}$$ This term feeds back on the $\theta$-independent contributions ${\mathbf{U}}_{\rm
  roll}$, ${\mathbf{U}}_{\rm streak}$ and ${\mathbf{U}}_{\rm mean}$. A visualization of this vector quantity is shown in Fig. \[fig:vis1\](c). On a qualitative level, by comparing the arrows of Fig. \[fig:vis1\](c) with those of Fig. \[fig:vis1\](a), one can see the feedback of this term on ${\mathbf{U}}_{\rm roll}$. The white-dashed boxes highlight regions in which the axial component of the Taylor-vortex flow is strong and aligned with the axial component of ${\bf \rm NL}$. The resemblance is especially strong on near-axial curves in ${\bf NL}$ converging towards saddles above and below regions with high azimuthal component shown in red.

A more quantitative picture of the feedback is presented in Fig. \[fig:bifdiag\](b). Shown is the normalized inner product between ${\rm
  NL} $ and each of ${\mathbf{U}}_{\rm roll}$, ${\mathbf{U}}_{\rm streak}$ and ${\mathbf{U}}_{\rm mean}$ defined by: $$\begin{aligned}
\langle {\bf NL},{\bf U}_{---}\rangle &=  \int_0^{L_z} dz \int_{r_1}^{r_2} r \:dr\: {\bf NL}(r,z)\cdot{\bf U}_{---}(r,z)
  \end{aligned}$$ where ${\bf U}_{---}$ is any of ${\mathbf{U}}_{\rm roll}$, ${\mathbf{U}}_{\rm streak}$ and ${\mathbf{U}}_{\rm mean}$. It can be seen that ${\bf NL}$ has a positive overlap with ${\bf U}_{\rm roll}$, meaning that indeed, the nonlinear interaction of ${\mathbf{u}}_{\rm wvf}$ with itself acts as a driving mechanism for rolls. ${\bf NL}$ also drives ${\mathbf{U}}_{\rm mean}$. In contrast, ${\bf NL}$ has a negative overlap with ${\bf U}_{\rm streak}$ and hence this term tends to suppress the streaks.

Conclusion
==========

According to the self-sustaining process (SSP) of [@waleffe1997self], the building block of transition to turbulence in plane Couette flow and other wall-bounded shear flows, rolls induce streaks, which in turn undergo an instability to waviness, whose nonlinear interaction feeds the rolls. In plane Couette flow, laminar flow (the analogue of ${\bf U}_{\rm Cou}$) is stable for all Reynolds numbers; there is no equivalent of the steady Taylor-vortex flow. For Taylor-vortex flow, however, most of the steps of the SSP are already in place. Vortices (rolls) induce streaks (axially periodic variation of the azimuthal flow) kinematically via advection, as in plane Couette flow. We have confirmed that the instability to wavy-vortex flow is due to this variation [@martinand2014mechanisms]. In addition, we have shown that the energy of the waves in nonlinear wavy-vortex flow compensates almost exactly for the decreased energy in the streaks, as compared to the energy in the streaks of nonlinear Taylor-vortex flow. The third step is the feedback of the waves on the rolls, which is crucial for the SSP since in plane Couette flow the rolls do not arise from a linear instability leading to a nonlinear equilibrium. We have shown that this feedback mechanism exists in Taylor-Couette flow and that it is the rolls that are fed and not the streaks. The nonlinear self-interaction of the waves generates localized regions with strong axial forcing: this is the nature of the feedback on the Taylor vortices which closes the SSP.

This research was partly supported by the grant TRANSFLOW, provided by the Agence Nationale de la Recherche (ANR).

[39]{}ifxundefined \[1\][ ifx[\#1]{} ]{}ifnum \[1\][ \#1firstoftwo secondoftwo ]{}ifx \[1\][ \#1firstoftwo secondoftwo ]{}““\#1””@noop \[0\][secondoftwo]{}sanitize@url \[0\][‘\
12‘\$12 ‘&12‘\#12‘12‘\_12‘%12]{}@startlink\[1\]@endlink\[0\]@bib@innerbibempty @noop [****,  ()]{} @noop [****,  ()]{} @noop [****, ()]{} @noop [****,  ()]{} @noop [****, ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****, ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****, ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****, ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} [****,  ()](\doibase 10.1103/PhysRevLett.103.144502) @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****,  ()]{} @noop [****, ()]{} @noop [****,  ()]{} @noop [****,  ()]{}
