<figure><img src="image_1.png" /></figure>

HAL
open science

# Intermittency and Critical Scaling in Annular Couette Flow

Kazuki Takeda, Yohann Duguet, Takahiro Tsukahara

► To cite this version:

Kazuki Takeda, Yohann Duguet, Takahiro Tsukahara. Intermittency and Critical Scaling in Annular Couette Flow. Entropy, 2020, 22 (9), pp.988. 10.3390/e22090988. hal-03032083

HAL Id: hal-03032083

https://hal.science/hal-03032083

Submitted on 30 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
Article

# Intermittency and Critical Scaling in Annular Couette Flow

Kazuki Takeda<fnref n="1" />, Yohann Duguet<fnref n="2" /> ID and Takahiro Tsukahara <fnref n="1" />,* ID

Received: 14 August 2020; Accepted: 1 September 2020; Published: 4 September 2020

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** The onset of turbulence in subcritical shear flows is one of the most puzzling manifestations of critical phenomena in fluid dynamics. The present study focuses on the Couette flow inside an infinitely long annular geometry where the inner rod moves with constant velocity and entrains fluid, by means of direct numerical simulation. Although for a radius ratio close to unity the system is similar to plane Couette flow, a qualitatively novel regime is identified for small radius ratio, featuring no oblique bands. An analysis of finite-size effects is carried out based on an artificial increase of the perimeter. Statistics of the turbulent fraction and of the laminar gap distributions are shown both with and without such confinement effects. For the wider domains, they display a cross-over from exponential to algebraic scaling. The data suggest that the onset of the original regime is consistent with the dynamics of one-dimensional directed percolation at onset, yet with additional frustration due to azimuthal confinement effects.

**Keywords:** subcritical phenomenon; transition to turbulence; direct numerical simulation

## 1. Introduction

The dynamics at the onset of turbulent fluid flow, as the parameters are varied, is one of the most puzzling issues of hydrodynamics. Subcritical flows are known to feature two regimes in competition, namely a laminar and a turbulent one. As the Reynolds number (their main control parameter) is varied, this competition takes the form of laminar-turbulent coexistence featuring some interesting analogies with phase transitions in thermodynamics. The onset of this coexistence in wall-bounded shear flows has been speculated to follow a statistical scenario called directed percolation (DP). It involves a critical point (a critical Reynolds number) in the vicinity of which fluctuations diverge algebraically [1,2]. The directed percolation scenario has gained theoretical importance because it appears as the usual rule for a one-dimensional systems obeying a set of specific properties, notably a unique absorbing state and short-range interactions [3,4]. However, it quickly proved difficult to isolate similar phenomena experimentally [5]. The main limitations happen to be finite-size effects, as well as the presence of defects [6–8] or issues revolving around nucleation rates [9,10]. The first experimental evidence for directed percolation in a two-dimensional physical system, with a complete set of critical exponents, occurred in electroconvection in nematic liquid crystals [11]. More recent experiments and numerical simulations with inert liquids were aimed at establishing the critical exponents relevant for the laminar-turbulent transition. The only meaningful experimental results are to be found in Ref. [12] for the flow inside an annulus driven by the revolutions of the outer wall, where all critical exponents match those of (1 + 1)-D DP. All other experimental attempts in effectively two-dimensional geometries have so far lead to ambiguous results [13,14]. A few numerical studies based on other geometries have also confirmed the DP hypothesis in one dimension, among them [15]. The most

<footnote n="1">¹ Department of Mechanical Engineering, Tokyo University of Science, Chiba 278-8510, Japan; tairyuuuu88@gmail.com</footnote>

<footnote n="2">² LIMSI-CNRS, Université Paris-Saclay, F-91400 Orsay, France; duguet@limsi.fr</footnote>

<footnote>* Correspondence: tsuka@rs.tus.ac.jp; Tel.: +81-4-7122-9352</footnote>
notorious system displaying one-dimensional spatiotemporal intermittency (STI) is cylindrical pipe flow. Although (1 + 1)-D DP has been widely speculated and is found in the most recent modelling approaches [16–18], clean experimental evidence seems to require facilities of a size beyond anything engineerable [19]. The only convincing two-dimensional study to date based on the (underresolved) Navier–Stokes equations and supporting the DP hypothesis is found in Ref. [20]. There again, a cost compromise was necessary between accuracy of the Navier–Stokes solutions and size effects. There the set of critical exponents differs from their unidimensional counterpart and corresponds to (2 + 1)-D DP. The status of the application of (2 + 1)-D DP to other planar flows is still open: for plane Couette flow (pCf), finite-size effects wrongly predict to discontinuous scenarios [21], whereas plane Poiseuille flow (pPf) seems to display a two-stage behavior so far poorly understood [22–24]. At a finite distance from the critical point, these two planar flows feature more structured arrays of turbulent stripes, all oblique to the mean flow direction (see, e.g., [21,25–30] for recent reviews).

Given the current status of DP affairs in shear flows, new flow candidates where to probe
the DP hypothesis are encouraged, irrespective of the effective dimension considered (one or two).
In the present article, we revisit transition in annular Couette flow (aCf) in the light of critical scaling.
This flow has a geometry similar to cylindrical pipe flow, however, with a solid cylinder at its centerline.
The geometry is determined by the radius ratio $\eta$ between the radius of the outer pipe and that of the
inner one. This flow supports both turbulence [31] as well as a linearly stable base flow for all Reynolds
number of interest, hence transition has to be of the subcritical type. Unlike annular pipe flow [32–34],
no pressure gradient is applied, instead the fluid is entrained by the translating motion of the inner
cylinder [35]. Earlier work by some of us [36] on this flow have lead to surprising results: although the
transitional flow reported for $\eta \ge 0.5$ consists of helical bands of turbulence wrapping around the inner
rod, for lower values of $\eta$, a new regime of laminar-turbulent alternations was reported. This regime is
characterized by slightly shorter streamwise correlations and non-oblique structures, explained by the
azimuthal confinement and by the impossibility to host azimuthal large-scale flows [37]. The aim of the
present article is to give a more detailed characterization of the novel low-$\eta$ intermittent regime and of
its onset. In particular, the azimuthal extension of aCf is investigated in a range of parameters beyond
that used by Kunii et al. [36]. As will be seen, this new choice of geometrical parameters leads to new
conclusions regarding the critical exponents. This new parametric study allows one to rationalize once
and for all the quantitative comparison between original geometry and the extended one.

The plan of this article unfolds as follows: the geometry and the numerical methods are explained
in Section 2, and the statistics of STI are reported in Section 3 and discussed in Section 4.

## 2. Set-Up and Methodology

### 2.1. *Geometry of aCf*

Annular Couette flow is the flow in the interstice between two coaxial cylinders of formally infinite length, driven by the motion at velocity $U_w > 0$ of the inner cylinder in the x-direction. The annular geometry of this flow is common to both Taylor–Couette flow and annular Pipe flow; however, the forcing is different and no spin of the walls is considered. A sketch of that geometry is displayed in Figure 1 with the usual notations for the cylindrical coordinates $(x, r, \theta)$. Assuming that the inner and outer cylinder have respective dimensional radii $r_{in}$ and $r_{out}$, the main geometrical parameter of this study is the radius ratio $\eta = r_{in}/r_{out}$, which varies in the open interval $(0, 1)$. We also introduce the gap $h$ between the two cylinders $h = r_{out} - r_{in}$.

Computationally, the pipes require to have either finite length or to be spatially periodic. The use
of a spectral Fourier-based method to solve the pressure Poisson equation requires axial and azimuthal
periodicity. This introduces the two wavelengths $L_x$ and $L_\theta$, respectively, as the domain length and
the angular periodicity. While $L_x$ is a free parameter, the natural value for $L_\theta$ is $2\pi$ because of the
cylindrical geometry. However, there is no computational obstruction to choosing other values for $L_\theta$,
for instance $L_\theta = 8\pi$ or $16\pi$ as in Ref. [36]. In what follows, we keep the generic notation $L_\theta$.
<figure><img src="image_3.png" /><figcaption>Figure 1. Sketch of annular Couette flow in the cylindrical coordinate system.</figcaption></figure>

Like in other wall-bounded shear flows, the main lengthscale ruling out the transitional dynamics at onset is the gap $h$ between the two solid walls, which here depends directly on the value $\eta$ via the relation $h = r_{out}(1 - \eta)$. The perimeter on the internal cylinder, at mid-gap or on the external cylinder, now expressed in units of $h$, is shown in Figure 2 when the original dimensional value of $L_{\theta}$ is $2\pi$ (Figure 2a). The inner perimeter is also displayed when $L_{\theta}$ is a multiple of $2\pi$ (Figure 2b), with $L_{\theta} = 2\pi n$. The theory developed in Refs. [33,36] shows that azimuthal large-scale flows cannot be accommodated by the geometry unless $L_{\theta}r/h \gg 1$ everywhere in the domain. The data for the inner cylinder play the role of a lower bound. For $L_{\theta} = 2\pi$, it is clear from Figure 2a that, for the lowest values of $\eta$, no azimuthal large-scale flow is possible. However, increasing $n$ leads to azimuthal large-scale flows being possible for smaller and smaller values of $\eta$. This leads to the possibility to artificially restore large-scale flows otherwise ruled out by geometrical confinement.

## 2.2. Governing Equations and Computational Methods

Whereas $\eta$ is a geometrical parameter only, we also introduce the Reynolds number $Re_w = U_w h / 4\nu$, based on the half velocity of the cylinder sliding $U_w / 2$, the half gap width $h / 2$, and the kinematic viscosity $\nu$ of the fluid. The reason why half-gap and half-velocities are considered to non-dimensionalize the equations is a simple way to reconnect with the standard conventions for pCf as $\eta$ goes towards unity. By choosing this convention for all values of $\eta$, the non-dimensional incompressible equations ruling the flow dynamics without any turbulence model read

$$ \nabla^* \cdot \mathbf{u}^* = 0, \qquad (1) $$

$$ \frac{\partial \mathbf{u}^*}{\partial t^*} + (\mathbf{u}^* \cdot \nabla^*) \mathbf{u}^* = -\nabla^* p^* + \frac{1}{4Re_w} \Delta^* \mathbf{u}^*, \qquad (2) $$

where superscripts $*$ indicate quantities non-dimensionalized with $U_w$ and $h$, and where $\mathbf{u} = (u_x, u_r, u_{\theta})$ and $p$ represent the velocity field and the pressure field, respectively.

Equation (2) is discretized in space using finite differences and with fine enough grid resolutions according to the standard criteria of direct numerical simulation (DNS) [26]. The time discretization is carried out using a second-order Crank–Nicolson scheme, and an Adams–Bashforth scheme for the wall-normal viscous term and the other terms, respectively. Further details about the numerical methods used here can be found in Ref. [38]. Table 1 lists the parameters used in this computational study.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) circumference of original annular pipe system at the outer cylinder, at mid-gap, and at the inner cylinder; (b) circumference at the inner cylinder for $L_\theta \ge 2\pi$.</figcaption></figure>

<table><caption>**Table 1.** Computational conditions. $L_i^*$: length of the computational domain in the direction $i$, non-dimensionalized by the gap width $h = (r_{out} - r_{in})$; $L_{out}^*$ (resp. $L_{in}^*$) the circumference of the outer (resp. inner) cylinder surfaces, normalized by $h$; $N_i$: the number of grids.</caption><thead><tr><th>η = r<sub>in</sub>/r<sub>out</sub></th><th>0.1</th><th>0.15</th><th>0.2</th><th>0.3</th></tr></thead><tbody><tr><td>L<sub>x</sub><sup>*</sup> × L<sub>r</sub><sup>*</sup></td><td colspan="4">512 × 1</td></tr><tr><td>L<sub>θ</sub></td><td>2π</td><td>32π</td><td>128π</td><td>409.6 × 1</td></tr><tr><td>L<sub>out</sub><sup>*</sup> (= L<sub>θ</sub> r<sub>out</sub><sup>*</sup>)</td><td>7.0</td><td>111.7</td><td>446.8</td><td>437.1</td></tr><tr><td>L<sub>in</sub><sup>*</sup> (= L<sub>θ</sub> r<sub>in</sub><sup>*</sup>)</td><td>0.7</td><td>11.2</td><td>44.7</td><td>71.0</td></tr><tr><td>N<sub>x</sub> × N<sub>r</sub></td><td colspan="4">2048 × 64</td></tr><tr><td>N<sub>θ</sub></td><td>32</td><td>512</td><td>2048</td><td>2048</td></tr></tbody></table>

### 3. Statistics at the Onset of Transition

#### 3.1. Global Stability and Coherent Structures Close to Onset

In the present subsection, we recall some key results of Ref. [36] together with some updated predictions. The investigation of the onset of turbulence starts with the determination of the global Reynolds number $Re_g$, defined as the highest Reynolds number below which no turbulence can survive (at least in the thermodynamic limit, i.e., over infinite observation times in unbounded domains). Since the flow is subcritical, using a given type of initial condition for this task can lead to overestimates of $Re_g$. The commonly adopted strategy, both in experiments and numerics, is that of an adiabatic descent [39] initiated from a turbulent state at sufficiently high Reynolds number. In the limit where the waiting time between successive diminutions of $Re$ is sufficient long, the value at which turbulence gets extinct is a good approximation of $Re_g$. Figure 3 displays information about $Re_g$ depending on the radius ratio $\eta$. For $L_\theta = 2\pi (n = 1)$, $Re_g$ increases monotonically with decreasing $\eta$. For larger $L_\theta$, $Re_g$ is always smaller than for the case with $L_\theta = 2\pi$ and the same value of $\eta$, with a now decreasing trend for $Re_g(\eta)$ which is even more marked once $\eta \le 0.3$. The values of $L_\theta$ needed to obtain this curve robustly are all listed in Table 1. As for the case of artificially extended aCf at $\eta = 0.1$, the result for $L_\theta = 128\pi$ is plotted in the figure. The parameter range strictly below $\eta = 0.1$ has not been investigated.
<figure><img src="image_5.png" /><figcaption>Figure 3. Radius ratio η dependency of the global critical Reynolds number Re<sub>g</sub>. The plot includes the pCf limit η → 1 from Ref. [21] (labeled "*1"), as well as DNS data from Ref. [36] for η = 0.5 and 0.8 (labeled "*2"). Triangles: original aCf with L<sub>θ</sub> = 2π is plotted using triangles; circles: artificially extended aCf (L<sub>θ</sub> > 2π).</figcaption></figure>

The fact that artificially extended systems display a lower threshold in Re indicates that some specific spatiotemporal regimes, specific to large L<sub>θ</sub> and not allowed for in narrow domains, are able to maintain themselves against relaminarization. As in Ref. [36], we can compare typical snapshots of the velocity fields in the corresponding regime in order to highlight the qualitative differences. Figures 4 and 5 display instantaneous snapshots of the radial velocity at mid-gap (i.e., r = (r<sub>in</sub> + r<sub>out</sub>)/2) at respectively η = 0.3 and 0.1, one very close to Re<sub>g</sub> (left column) and the other slightly above it (right column). Each row corresponds to a different value of the integer n (n = 1, 16, 48, and 64), i.e., another value of L<sub>θ</sub>. When n = 1, the one-dimensional intermittency is reminiscent of the dynamics in cylindrical pipe flow [40]. The differences between different values of η emerge only for higher n. For η = 0.3, the stripe patterns exhibit an obliqueness typical of most laminar-turbulent patterns [25,26,37,41]. However, it is visually clear that the situation is different for η = 0.1, with shorter structures and less pronounced obliqueness. It is not immediately clear whether the effective dimensionality of the proliferation process is rather one or two. These issues can be addressed using the determination of critical exponents, as will be done in the next subsection.

## 3.2. Data Binarization

Velocity fluctuations with respect to the mean flow are defined as $u' = u - \bar{u}$, where $\bar{u}$ is the space-averaged time-dependent velocity averaged along x and θ, as defined in Equation (3). Here, y denotes the (dimensional) distance from the inner cylinder to the outer cylinder as y = r - r<sub>in</sub>, instead of using r.

$$ \bar{u}(y,t) = \frac{1}{L_x L_\theta} \int_0^{L_x} \int_0^{L_\theta} u(x, y, \theta, t) \, dx \, d\theta. \quad (3) $$

The flow is separated into its laminar and turbulent components by postulating a threshold independently of the Reynolds number. The local criterion chosen is $|u'_r/U_w| \ge 0.01$ for turbulence and $|u'_r/U_w| < 0.01$ for laminar flow, with $u'_r$ the radial velocity component, which vanishes everywhere for strictly laminar flow. As in Figures 4 and 5, localized turbulent regions are visualized by contours of $u'_r^*$ in steps of ±0.01. The turbulent fraction $F_t$ is evaluated at mid-gap (y = h/2) by estimating the percentage of grid points for which the turbulent criterion above is fulfilled.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Contours of radial velocity fluctuations $u_r^*$ at mid-gap for $\eta = 0.3$ around $Re_w = Re_g$. Typical snapshots of instantaneous flow fields obtained after reaching each equilibrium state are shown here. The main flow is from left to right. (a,b) original aCf with $L_\theta = 2\pi$, and (c,d) artificially extended with $L_\theta = 96\pi$.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** The same as Figure 4, but for $\eta = 0.1$. (a,b) $L_\theta = 2\pi$; (c,d) $L_\theta = 32\pi$; and (e,f) $L_\theta = 128\pi$.</figcaption></figure>

The dynamics of the proliferation process for $\eta = 0.1$ and $0.3$ is illustrated in Figure 6 using space-time diagrams and compared one to another in the case $n = 1$. The spatial variable is $x - U_{ft}$,
i.e., the streamwise coordinate in a frame moving with constant velocity $U_f$, which is close to the average velocity $\overline{u_m}$. The space-time diagram is based on the binarized radial velocity $u'_r$. The absolute value of the radial velocity evaluated at mid-gap is first averaged azimuthally according to

$$ \langle u'_{r\text{rms}} \rangle_{\theta}(x, t) = \sqrt{\frac{1}{2\pi} \int_{0}^{2\pi} u'^2(x, h/2, \theta, t) \, d\theta}. \quad (4) $$

and the binarization criterion is $\langle u'_{r\text{rms}} \rangle_{\theta}/U_w \ge 0.01$. The frame velocity $U_f$ for $\eta = 0.3$ is chosen to be same with $\overline{u_m}$, which is estimated in two steps. First, a spatially average velocity is evaluated at every time $t$

$$ u_m(t) = \frac{1}{L_x(r_{\text{out}}^2 - r_{\text{in}}^2)L_{\theta}} \int_0^{L_x} \int_{r_{\text{in}}}^{r_{\text{out}}} \int_0^{L_{\theta}} u_x(x, y, \theta, t) r \, dx dr d\theta, \quad (5) $$

then it is time-averaged using a classical moving average technique over a time interval $\Delta T$ (with $\Delta T > 10^4 h/U_w$ after reaching equilibrium).

$$ \overline{u_m} = \frac{1}{\Delta T} \int_T^{T+\Delta T} u_m(t) \, dt. \quad (6) $$

We found that, for $\eta = 0.1$, an optimal value of $U_f$ for the frame to move with puffs was slightly slower than $\overline{u_m}$. For each value of $\eta$, three space-time diagrams are displayed, respectively below, close to and above the corresponding critical point $Re_g(\eta)$. The shorter aspect of the coherent structures for $\eta = 0.1$ is striking compared to $\eta = 0.3$. Many more splitting and decay events, qualitatively similar to the pipe flow case [40,42,43], occur for $\eta = 0.1$ despite equal pipe lengths. This suggests that the status of the present simulations for $\eta = 0.1$ is qualitatively much closer to the thermodynamic limit than it is for $\eta = 0.3$. As a by-product, the critical scaling is expected to converge at a lower price than at higher $\eta$. Given the cost obstacles induced by the diverging lengthscales/timescales in most critical phenomena, the above conclusion is positive news.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Space-time $(x - U_f t)/h$ diagram of original aCf ($L_\theta = 2\pi$) for $\eta = 0.1$ (three leftmost columns) and 0.3 (three rightmost columns). Black: turbulence according to the criterion $\langle u'_{r\text{rms}} \rangle_{\theta}/U_w \ge 0.01$. The values of the frame velocity $U_f$ for $\eta = 0.1$ are 0.288$U_w$ at $Re_w = 407$, 0.2875$U_w$ at $Re_w = 407.5$, 0.2815$U_w$ at $Re_w = 415$, and those for $\eta = 0.3$ are approximately equal to $\overline{u_m}$.</figcaption></figure>
### 3.3. Intermittency Statistics

The statistical post-processing protocol for STI is vastly similar to that used by other authors: the first step is to monitor the decay in the time of the turbulence fraction $F_t(t)$ when the system is initiated with turbulence everywhere. By dichotomy, this yields a good approximation of $Re_g$ and allows one to define the reduced control parameter $\epsilon = (Re_w - Re_g)/Re_g$. This decay is expected to be algebraic exactly at onset, i.e., of the form $F_t(t) = O(t^{-\alpha})$. This yields as well the so-called dynamic exponent $\alpha$. In a second phase, the equilibrium turbulent fraction (i.e., its time average) is monitored as a function of $\epsilon$. For $\epsilon > 0$, the data versus the expected scaling $F_t(t) = O(\epsilon^\beta)$ yield the exponent $\beta$. Eventually, the mean correlation length $\xi(Re_w)$ (either $\xi_x$ in the streamwise direction or $\xi_\theta$ in the azimuthal one) can be estimated at equilibrium by monitoring the cumulative distribution function (CDF) of the laminar gaps $P_{lam}(l_x > L)$, where $l_x$ stands for the length of a laminar trough and $L$ is a dummy variable. A critical exponent $\mu_\perp$ can be evaluated from fits as the algebraic decay exponent of the CDF.

We begin by describing the results from the critical quench experiments of Figure 7 for $\eta = 0.1$ and $n = 64$. The initial condition corresponds to a turbulent velocity field from a long simulation well above $Re_g$, here taken as $Re_w = 280$. The same initial condition is used for new simulations at another target value of $Re_w$, in principle such that $Re_w$ is “close” to $Re_g$. As expected, the flow relaminarizes (attested by the monotonic decrease of $F_t(t)$) for sufficiently low values of $Re_w$, whereas it stays turbulent for the higher values. In the latter case, the turbulent fraction reaches a non-zero mean value $\bar{F}_t$, which will be reported in the next figure. The set of colored curves in Figure 7a straddle the decay curve corresponding to the critical value $Re_w = Re_g$, whose best approximation in the figure is the red curve associated with $Re_w = 262.5$. For continuous phase transitions, the corresponding decay is expected to be of power-law type, i.e., $F_t = O(t^{-\alpha})$. This fact of $260 < Re_g < 262.5$ yields an approximation of $Re_g = 261.7$, which allows for defining $\epsilon$ as before. The present approach rests on the hypothesis of a critical scaling in the vicinity of the critical point. If that hypothesis is correct then, by rescaling time and turbulent fraction, the curves of Figure 7a should collapse onto two master curves, one for the relaminarization process and the other for the saturation process. This is tested in Figure 7b by plotting $t^\alpha F_t(t)$ as a function of the rescaled time $t|\epsilon|^{\nu\|}$. As for $\alpha$ and $\nu_\|$ the approximate values from (1 + 1)-D DP theory, respectively 0.451 and 1.733, have been used for the rescaling. The match is satisfying, which confirms that a critical range has been identified in this system.

As a by-product of Figure 7, the values of the mean turbulent fraction $\bar{F}_t$, obtained after reaching equilibrium, are reported in Figure 8 as functions of $Re_w$. Critical theories all predict a scaling $\bar{F}_t = O(\epsilon^\beta)$ close enough to the critical point. The algebraic scaling revealed in the previous plots of critical quench suggests that, for instance, $Re = 262.5$ belongs to the range where algebraic fits apply for $\eta = 0.1$ and $L_\theta = 128\pi$. Consequently, if, for these parameters, $\epsilon$ is defined using the approximated $Re_g = 261.7$, the dependence of $\bar{F}_t$ versus $\epsilon$ is also expected to be algebraic in the same range of values of $Re$. In that case, the power-law exponents can be classically estimated using log-log plots and compared to those from DP theories. Algebraic fits of $\bar{F}_t$ are shown in Figure 8 both for $\eta = 0.1$ (left) and 0.3 (right). For each case, the main plot of $\bar{F}_t$ versus $Re_w$ is displayed in linear coordinates, while the inset displays $\bar{F}_t$ versus $\epsilon$ in log-log coordinates, in order to highlight the quality of the estimation of the power-law exponent.
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Critical quenches from $Re_{w} = 280$ to each Reynolds number. Temporal variation of turbulent fraction $F_{t}$ for $\eta = 0.1$ and $L_{\theta} = 128\pi$ (log-log scale). In (a), the black dashed-dotted line and dashed line each indicate possible algebraic fits with the dynamic exponent $\alpha$ from (2 + 1)-D and (1 + 1)-D directed percolation (respectively $\alpha = 0.451$ and 0.159). See also Supplementary; (b) test of the 1D scaling hypothesis by plotting $t^{\alpha}F_{t}$ vs. $t\epsilon^{v_{\parallel}}$ (log-log scale), with $v_{\parallel}=1.733$ for (1 + 1)-D DP.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Reynolds-number dependence of the time-averaged turbulent fraction $\bar{F}_{t}$ vs $Re_{w}$ for the different radius ratios in the original domain ($L_{\theta} = 2\pi$) and in artificially extended domains ($L_{\theta} \gg 2\pi$). Vertical error bars: standard deviations of $F_{t}$ during the averaging period. Dashed/dashed-dotted line: algebraic fits $\bar{F}_{t} = O(\epsilon^{\beta})$, with exponent $\beta$ obtained either as best fit $\beta_{\text{fit}}$ or from the (1 + 1)-D DP universality class $\beta_{1D} = 0.276$. In each figure, the insets are plotted in log-log coordinates versus $\epsilon$ that is determined with $Re_{g}$ presented in Table 2.</figcaption></figure>

The details of the fitting procedure for the various parameters used are given in Table 2. It includes the values of the best fitted exponents as well as the approximate fitting range. As could already be deduced graphically from the insets in Figure 8a, for $\eta = 0.1$, the compatibility of the exponent $\beta$ with the theoretical value of $\beta_{1D} = 0.276$ from (1 + 1)-D DP is good (to the second digit). This is confirmed
for both $\eta = 32\pi$ and $\eta = 128\pi$, which suggests that the thermodynamic limit is already reached, at least as far as the determination of the exponent $\beta$ is concerned. For $L_\theta = 2\pi$ the approximated exponent is 0.31 which constitutes a less accurate, but still consistent approximation of the theoretical exponent. For $L_\theta = 2\pi$, the range of validity of the algebraic fits extends up to $\approx 5\%$, whereas it exceeds 10% for $L_\theta \ge 32\pi$. For $\eta = 0.3$, the situation is slightly different: for a large azimuthal extent $L_\theta = 96\pi$, there is a very good match with the 1D theoretical exponent all the way up to $\epsilon \approx 20\%$. For $L_\theta = 2\pi$, however, although an algebraic fit seems consistent with the data below $\epsilon < 1\%$ the measured exponent is closer to 0.12 than to 0.276: none of these values matches any of the percolation theories.

<table><caption>**Table 2.** Critical Reynolds number $Re_g$ and critical exponent $\beta$ depending on geometrical parameters $\eta$ (radius ratio) and $L_\theta$ (azimuthal extension). In addition, shown is the fitting range to estimate $Re_g$ and $\beta$. $^\dagger$: not measured.</caption><thead><tr><th>η = r<sub>in</sub>/r<sub>out</sub></th><th>L<sub>θ</sub></th><th>Fitting Range</th><th>Re<sub>g</sub></th><th>β</th></tr></thead><tbody><tr><td>0.10</td><td>2π</td><td>407.5–460.0</td><td>406.9</td><td>0.31(3)</td></tr><tr><td>0.10</td><td>32π</td><td>277.5–300.0</td><td>269.0</td><td>0.26(2)</td></tr><tr><td>0.10</td><td>128π</td><td>263.0–270.0</td><td>261.7</td><td>0.28(2)</td></tr><tr><td>0.15</td><td>128π</td><td>—<sup>†</sup></td><td>290.5</td><td>—<sup>†</sup></td></tr><tr><td>0.20</td><td>112π</td><td>—<sup>†</sup></td><td>303.5</td><td>—<sup>†</sup></td></tr><tr><td>0.30</td><td>2π</td><td>389.0–395.0</td><td>388.7</td><td>0.12(2)</td></tr><tr><td>0.30</td><td>96π</td><td>320.5–375.0</td><td>319.0</td><td>0.28(1)</td></tr><tr><td>(1 + 1)-D DP model</td><td></td><td>—</td><td>—</td><td>0.276</td></tr></tbody></table>

The interpretation is delicate. On one hand, algebraic fits seem always verified as soon as $\epsilon$ is small enough; on the other hand, (1 + 1)-D percolation exponents are well approximated only for sufficient azimuthal extension of the order of $100\pi$ or more. The original system with $L_\theta = 2\pi$ hence needs to be interpreted as a system with the DP property that experiences a geometrical frustration due to lateral confinement. The present data support the hypothesis that the frustration effect is stronger for $\eta = 0.3$ than for $\eta = 0.1$, and thus that the quality of the DP fit will be correspondingly worse. Conversely, the convergence towards the thermodynamic limit seems slower for larger $\eta$.

Importantly, we emphasize the main difference between the present conclusion and that by Kunii et al. [36], where the azimuthal extension for $\eta = 0.1$ was limited to $L_\theta = 16\pi$ (to be compared to the present values of $32\pi$ and $128\pi$). The fits reported in Figure 16 of that article suggested a fit compatible with the (2 + 1)-D exponent $\beta_{2D} = 0.583$. This former result, in the light of the present computations, is re-interpreted now as a finite-size effect.

A power-law dependence of $\bar{F}_t$ alone does not warrant the proximity to the critical point, as pointed out by Shimizu and Manneville [23] for pPf. Although the critical quenches reported earlier also suggest power-law statistics near the picked up values for $Re_g$, the classical determination relies on, at least, three independent algebraic exponents. In order to lift this ambiguity, we chose to report in Figure 9 statistics of laminar gap size for different values of $Re_w$ near the suspected critical point. Expecting possible anisotropy when the domain is artificially extended in $\theta$, two kinds of statistics have been monitored, similarly to the study of Chantry et al. [20]. The axial extent of the gaps for $\eta = 0.1$ and $L_\theta = 128\pi$ is shown in Figure 9a in log-log coordinates (and Figure 9c in lin-log representation). The azimuthal extent of the laminar gaps is shown in Figure 9b in log-log coordinates (and Figure 9d in lin-log representation). All four figures support a cross-over from exponential to power-law statistics as $Re_w$ approaches the value of 262.5, with a decay exponent graphically compatible with the decay exponent $\mu_\perp$ of (1 + 1)-D DP. The cross-over appears, however, more clearly in the azimuthal where the match with the theoretical value of $\mu_\perp$ is valid over a full decade. In the streamwise direction, the trend is not clear enough to extract a critical exponent with full accuracy. This confirms, however, that the present statistics are indeed gathered in a relevant neighborhood
of the critical point and that, for these parameters, $Re_w = 262.5$ is a decent working approximation of $Re_g$.

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Time-averaged distributions of laminar gap in (a,c) the streamwise direction and (b,d) the azimuthal direction, evaluated at mid-gap. (a,b) log-log plots vs. (c,d) lin-log plots. $L_\theta = 128\pi$, $\eta = 0.1$ as in Figure 5e,f. In both figures, black dashed-dotted line (· · ·) and dashed lien (---) indicate theoretical distributions $P(\Delta L^*) \sim \Delta L^{*-\mu_\perp}$ with exponents $\mu_\perp$ from the universality classes of (2 + 1)-D DP and (1 + 1)-D DP, respectively, i.e., $\mu_{\perp 2D} = 1.84$, and $\mu_{\perp 1D} = 1.748$.</figcaption></figure>

### 3.4. Dynamics of Localized Turbulent Patches

In this last subsection, we address the issue of the influence of azimuthal confinement/extension on the lower transition threshold $Re_g$, as the estimations from Figure 3 suggest. In Ref. [36], a similar trend was noted (from measurements in shorter and narrower domains). The mechanism suggested in this former work addressed the presence of oblique stripes rather than their influence on the value of $Re_g$. It was thereafter realized that the phenomenon governing the value of $Re_g$, and by extension all statistics of the turbulent fraction, is the way different coherent structures interact together dynamically rather than the shape of such individual structures (although that shape certainly influences the interactions). In analogy with pipe [16,42,43] and channel [44,45], the finite turbulent fraction is the result of a dynamical competition between the proliferation of coherent structures and their tendency to decay in number. The transitional range where $\bar{F}_t > 0$ is dominated by the splitting of coherent structures, whereas instantaneous relaminarizations become rare. We hence focus on the
dynamics of splitting events in two different computational domains, namely those with $L_\theta = 32\pi$ and $128\pi$. Figure 10 contains zooms on the radial velocity plotted for different values of $y = \text{cst}$ surfaces (a different value for each row) and for different times (different columns). In Figure 10, the value of $L_\theta$ is fixed to $32\pi$, but the circumference in terms of $r\theta/h$ varies according to $r$. The global dynamics of these flows can also been scrutinized in the videos made available as Supplementary Materials. The comparison of different values of $y$ is useful to confirm that, for all parameters, the spots remain coherent over the gap even during splitting events.

Lateral splitting events are considered in each of these figures and videos. Because of the different advection velocities in the azimuthal direction, spanwise collisions can occur. During spanwise collisions, usually one of the two spots disappears (see also Ref. [21] for similar observations in pCf). This tends to reduce the turbulent fraction while the other surviving spot is still active. In the presence of a short enough spanwise periodicity, a spot collides with itself rather than with a different neighbor. In such periodic domains, the local relaminarization of one spot is equivalent to the extinction of an infinity of identical spots. Hence, the turbulent fraction decreases more than in large domains where individual spots behave more like independent entities. We thus expect more turbulence to proliferate more for larger $L_\theta$. As a consequence, the critical Reynolds number $Re_g$, for which the rate of proliferation balances the probability to relaminarize locally, is lowered when $L_\theta$ is increased, consistently with the thresholds reported in Figure 3 and Table 2. This effect is more marked at lower $\eta$.

<figure><img src="image_12.png" /><figcaption>**Figure 10.** Snapshots of splitting and self-colliding events in aCf for $Re_w = 262.5$ with $L_\theta = 32\pi$ and $\eta = 0.1$. Radial velocity in a frame moving with bulk velocity $\bar{u}_m$. Here, $t = 0$ is an arbitrary time instant after reaching equilibrium. Top row, $y^* = y/h \approx 0.9$; center row, $y^* \approx 0.5$; lower row, $y^* \approx 0.1$.</figcaption></figure>

4. Conclusions

The present DNS study deals with the statistical aspects of the intermittent transitional regime of aCf, with an emphasis on the low values of the radius ratio $\eta$ close to 0.1. It is an extension of
the simulations reported recently by Ref. [36]. The paper compares two computational situations, respectively the case of a realistic geometry and the one where the azimuthal extent is larger than the original value of $2\pi$. In Ref. [36], this parametric trick was introduced in an explicit attempt to decouple the effects of wall curvature effects from the effects of azimuthal confinement induced by the geometry. The main conclusion for large $\eta$ was that the reported absence of oblique laminar-turbulent patterns was due to azimuthal confinement, since they could re-appear for $L_\theta > 2\pi$. In the present article, the same trick is introduced for $\eta = 0.1$; however, larger values of $L_\theta$ have been tried up to $128\pi$ (i.e., 64 times the original value). The oblique patterns do not reappear and a new percolating regime takes place with shorter spatial correlations. The statistical analysis of the STI is convergent as $L_\theta$ grows, and is consistent with (1 + 1)-D DP. This updates the results of Ref. [36] where (2 + 1)-D DP was suggested from fits with $L_\theta = 16\pi$. The present results suggest now that the $L_\theta = 16\pi$ algebraic statistics was still far from the true thermodynamic limit, while $L_\theta = 128\pi$ seems to yield more decent results.

To our knowledge, there has been only poor evidence for the cross-over from exponential to algebraic scaling in the shear flow literature, as far as well-resolved simulations of the Navier–Stokes equations are concerned [2]. An exception is the work by Shi et al. [46] in a tilted periodic domain of pCf, which again is not a fully realistic numerical domain. It is interesting to speculate how much the present results can teach us something about a fully realistic system such as cylindrical pipe flow. Naive homotopy of the turbulent regimes is ruled out because of the singularity near the centerline. Instead, we can compare the rate at which these two effectively one-dimensional percolating systems tend towards their own thermodynamic limit. This issue was raised recently in the experimental study by Mukund and Hof [19]. There, despite pipes as long as 3000 diameters, no critical regime (with power-law statistics) was identified, only classical STI as reported in Refs. [47,48]. This issue was attributed to the narrowness of the critical range, and to a clustering property of puffs which delays the convergence to the thermodynamic limit. Here, in aCf with $\eta = 0.1$, the situation is different but depends on this artificial parameter $L_\theta$. To our surprise, power-law statistics of the turbulent fraction as well as of the laminar gap distributions do appear in our simulations as $Re_w$ is reduced. All cases shown in Figure 8 suggest a cross-over from turbulent to power-law behavior as $Re_w$ is within ≈ 1% of the critical point. For $L_\theta = 2\pi$ or around, the turbulent fraction curve still suggests an unconverged power-law. For $L_\theta = 32\pi$ or $128\pi$, power-law statistics of $\bar{F}_t$ are fully consistent with one-dimensional DP appear. This occurs despite a value of $L_x$ of only 512h, i.e., much less than the pipe flow case and even less if one counts in outer pipe diameters. A possible interpretation is that azimuthal extension, by modifying the interaction with neighboring spots, can suppress the tendency to form clusters, and hence converge faster towards the thermodynamic limit. This is consistent with lower transition thresholds in $Re_w$ as well. One is left wondering if a similar approach to cylindrical pipe flow could also easily yield the percolation exponents from simulation measurements.

We conclude by noting that artificially modifying both the shape of turbulent patches and their interaction, as done here using azimuthal extension, is more than an esoteric thought experiment or an exotic parameter study. It is used here as a legitimate strategy in order to untangle complex phenomena, e.g., to decouple confinement from curvature effects. As demonstrated in our recent work using a simple modeling approach [49,50], wall roughness can have similar effects on transitional flows and change the way turbulence invades laminar flows. We expect similar strategies of artificial domain extension to be relevant to such cases too.

**Supplementary Materials:** Video S1: Time evolution of turbulent fraction $F_t(t)$ and of fluctuating velocity fields visualized at mid-gap, for $\eta = 0.1$ with an artificially extended azimuthal domain size of $L_\theta = 128\pi$. On the right column, contours show $x-\theta$ distributions of the radial velocity fluctuation $u'_r$ normalized by the inner-cylinder velocity $U_w$. Top (orange box and curve in the graph): above the global critical Reynolds number $Re_g$. Middle (red): near $Re_g$. Bottom (black): below $Re_g$. A supporting video article is available at https://doi.org/10.5281/zenodo.3985963.
**Author Contributions:** Conceptualization and methodology, T.T. and Y.D.; simulations and acquisition, K.T.; post-processing, K.T. and T.T.; writing Y.D. and T.T.; funding acquisition, T.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was funded by Grant-in-Aid for JSPS (Japan Society for the Promotion of Science) Fellowship JP16H06066 and JP19H02071.

**Acknowledgments:** Numerical simulations were performed on SX-ACE supercomputers at the Cybermedia Center of Osaka University and the Cyberscience Center of Tohoku University.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## Abbreviations

The following abbreviations are used in this manuscript:

aCf annular Couette flow

CDF cumulative distribution function

DNS direct numerical simulation

DP direct percolation

pCf plane Couette flow

pPf plane Poiseuille flow

rms root-mean-square value

STI spatiotemporal intermittency

## References

1. Pomeau, Y. Front motion, metastability and subcritical bifurcations in hydrodynamics. *Phys. D Nonlinear Phenom.* **1986**, *23*, 3–11. [CrossRef]

2. Eckhardt, B. Transition to turbulence in shear flows. *Phys. A Stat. Mech. Appl.* **2018**, *504*, 121–129. [CrossRef]

3. Grassberger, P. On phase transitions in Schlögl’s second model. *Zeitschrift für Physik B Condensed Matter* **1982**, *47*, 365–374. [CrossRef]

4. Janssen, H.K. On the nonequilibrium phase transition in reaction-diffusion systems with an absorbing stationary state. *Zeitschrift für Physik B Condensed Matter* **1981**, *42*, 151–154. [CrossRef]

5. Daviaud, F.; Dubois, M.; Bergé, P. Spatio-temporal intermittency in quasi one-dimensional Rayleigh-Bénard convection. *EPL (Europhys. Lett.)* **1989**, *9*, 441. [CrossRef]

6. Chaté, H.; Manneville, P. Spatio-temporal intermittency in coupled map lattices. *Phys. D Nonlinear Phenom.* **1988**, *32*, 409–422. [CrossRef]

7. Bohr, T.; van Hecke, M.; Mikkelsen, R.; Ipsen, M. Breakdown of universality in transitions to spatiotemporal chaos. *Phys. Rev. Lett.* **2001**, *86*, 5482. [CrossRef]

8. Duguet, Y.; Maistrenko, Y.L. Loss of coherence among coupled oscillators: From defect states to phase turbulence. *Chaos Interdiscip. J. Nonlinear Sci.* **2019**, *29*, 121103. [CrossRef]

9. Daviaud, F.; Bonetti, M.; Dubois, M. Transition to turbulence via spatiotemporal intermittency in one-dimensional Rayleigh-Bénard convection. *Phys. Rev. A* **1990**, *42*, 3388. [CrossRef]

10. Kreilos, T.; Khapko, T.; Schlatter, P.; Duguet, Y.; Henningson, D.S.; Eckhardt, B. Bypass transition and spot nucleation in boundary layers. *Phys. Rev. Fluids* **2016**, *1*, 043602. [CrossRef]

11. Takeuchi, K.A.; Kuroda, M.; Chaté, H.; Sano, M. Directed percolation criticality in turbulent liquid crystals. *Phys. Rev. Lett.* **2007**, *99*, 234503. [CrossRef] [PubMed]

12. Lemoult, G.; Shi, L.; Avila, K.; Jalikop, S.V.; Avila, M.; Hof, B. Directed percolation phase transition to sustained turbulence in Couette flow. *Nat. Phys.* **2016**, *12*, 254–258. [CrossRef]

13. Sano, M.; Tamai, K. A universal transition to turbulence in channel flow. *Nat. Phys.* **2016**, *12*, 249–253. [CrossRef]

14. Avila, K. Shear Flow Experiments: Characterizing the Onset of Turbulence as a Phase Transition. Ph.D. Thesis, Niedersächsische Staats-und Universitätsbibliothek Göttingen, Göttingen, Germany, 2013.

15. Hiruta, Y.; Toh, S. Subcritical laminar–turbulent transition as nonequilibrium phase transition in two-dimensional Kolmogorov flow. *J. Phys. Soc. Jpn.* **2020**, *89*, 044402. [CrossRef]
16. Barkley, D. Theoretical perspective on the route to turbulence in a pipe. *J. Fluid Mech.* **2016**, *803*, P1.
[CrossRef]

17. Linga, G. Fluid Flows with Complex Interfaces: Modelling and Simulation from Pore to Pipe. Ph.D. Thesis,
The Niels Bohr Institute, Faculty of Science, University of Copenhagen, København, Danmark, 2018.

18. Shih, H.-Y.; Lemoult, G.; Vasudevan, M.; Lopez, J.; Shih, H.Y.; Linga, G.; Hof, B.; Mathiesen, J.; Goldenfeld, N.; Hof, B. Statistical model and universality class for interacting puffs in transitional turbulence. *Bull. Am. Phys. Soc.* **2020**, *65*, W25.00007.

19. Mukund, V.; Hof, B. The critical point of the transition to turbulence in pipe flow. *J. Fluid Mech.* **2018**, *839*, 76-94. [CrossRef]

20. Chantry, M.; Tuckerman, L.S.; Barkley, D. Universal continuous transition to turbulence in a planar shear flow. *J. Fluid Mech.* **2017**, *824*, R1. [CrossRef]

21. Duguet, Y.; Schlatter, P.; Henningson, D.S. Formation of turbulent patterns near the onset of transition in plane Couette flow. *J. Fluid Mech.* **2010**, *650*, 119-129. [CrossRef]

22. Tao, J.; Eckhardt, B.; Xiong, X. Extended localized structures and the onset of turbulence in channel flow.
*Phys. Rev. Fluids* **2018**, *3*, 011902. [CrossRef]

23. Shimizu, M.; Manneville, P. Bifurcations to turbulence in transitional channel flow. *Phys. Rev. Fluids* **2019**, *4*, 113903. [CrossRef]

24. Kashyap, P.V.; Duguet, Y.; Dauchot, O. Flow statistics in the patterning regime of plane channel flow. *Entropy* **2020**, submitted.

25. Prigent, A.; Grégoire, G.; Chaté, H.; Dauchot, O.; van Saarloos, W. Large-scale finite-wavelength modulation within turbulent shear flows. *Phys. Rev. Lett.* **2002**, *89*, 014501. [CrossRef]

26. Tsukahara, T.; Seki, Y.; Kawamura, H.; Tochio, D. DNS of turbulent channel flow at very low Reynolds numbers. In *TSFP Digital Library Online*; Begel House Inc.: Danbury, CT, USA, 2005; pp. 935-940,

27. Brethouwer, G.; Duguet, Y.; Schlatter, P. Turbulent-laminar coexistence in wall flows with Coriolis, buoyancy or Lorentz forces. *J. Fluid Mech.* **2012**, *704*, 137-172. [CrossRef]

28. Manneville, P. Transition to turbulence in wall-bounded flows: Where do we stand? *Mech. Eng. Rev.* **2016**, *3*.
[CrossRef]

29. Manneville, P. Laminar-turbulent patterning in transitional flows. *Entropy* **2017**, *19*, 316. [CrossRef]

30. Tuckerman, L.S.; Chantry, M.; Barkley, D. Patterns in Wall-Bounded Shear Flows. *Annu. Rev. Fluid Mech.* **2020**, *52*, 343-367. [CrossRef]

31. Chung, S.Y.; Rhee, G.H.; Sung, H.J. Direct numerical simulation of turbulent concentric annular pipe flow: Part 1: Flow field. *Int. J. Heat Fluid Flow* **2002**, *23*, 426-440. [CrossRef]

32. Ishida, T.; Duguet, Y.; Tsukahara, T. Transitional structures in annular Poiseuille flow depending on radius ratio. *J. Fluid Mech.* **2016**, *794*, R2. [CrossRef]

33. Ishida, T.; Duguet, Y.; Tsukahara, T. Turbulent bifurcations in intermittent shear flows: From puffs to oblique stripes. *Phys. Rev. Fluids* **2017**, *2*, 073902. [CrossRef]

34. Fukuda, T.; Tsukahara, T. Heat transfer of transitional regime with helical turbulence in annular flow. *Int. J. Heat Fluid Flow* **2020**, *82*, 108555. [CrossRef]

35. Frei, C.; Lüscher, P.; Wintermantel, E. Thread-annular flow in vertical pipes. *J. Fluid Mech.* **2000**, *410*, 185-210.
[CrossRef]

36. Kunii, K.; Ishida, T.; Duguet, Y.; Tsukahara, T. Laminar-turbulent coexistence in annular Couette flow. *J. Fluid Mech.* **2019**, *879*, 579-603. [CrossRef]

37. Duguet, Y.; Schlatter, P. Oblique laminar-turbulent interfaces in plane shear flows. *Phys. Rev. Lett.* **2013**, *110*, 034502. [CrossRef]

38. Abe, H.; Kawamura, H.; Matsuo, Y. Direct numerical simulation of a fully developed turbulent channel flow with respect to the Reynolds number dependence. *J. Fluids Eng.* **2001**, *123*, 382-393. [CrossRef]

39. Prigent, A.; Dauchot, O. Transition to versus from turbulence in subcritical Couette flows. In *IUTAM Symposium on Laminar-Turbulent Transition and Finite Amplitude Solutions*; Springer: Dordrecht, The Netherlands, 2005; pp. 195-219.

40. Moxey, D.; Barkley, D. Distinct large-scale turbulent-laminar states in transitional pipe flow. *Proc. Natl. Acad. Sci. USA* **2010**, *107*, 8091-8096. [CrossRef]

41. Fukudome, K.; Tsukahara, T.; Ogami, Y. Heat and momentum transfer of turbulent stripe in transitional-regime plane Couette flow. *Int. J. Adv. Eng. Sci. Appl. Math.* **2018**, *10*, 291-298. [CrossRef]
42. Avila, K.; Moxey, D.; de Lozar, A.; Avila, M.; Barkley, D.; Hof, B. The onset of turbulence in pipe flow. *Science* **2011**, *333*, 192–196. [CrossRef]

43. Shimizu, M.; Manneville, P.; Duguet, Y.; Kawahara, G. Splitting of a turbulent puff in pipe flow. *Fluid Dyn. Res.* **2014**, *46*, 061403. [CrossRef]

44. Shimizu, M.; Kanazawa, T.; Kawahara, G. Exponential growth of lifetime of localized turbulence with its extent in channel flow. *Fluid Dyn. Res.* **2019**, *51*, 011404. [CrossRef]

45. Paranjape, C.S.; Duguet, Y.; Hof, B. Oblique stripe solutions of channel flow. *J. Fluid Mech.* **2020**, *897*, [CrossRef]

46. Shi, L.; Avila, M.; Hof, B. Scale invariance at the onset of turbulence in Couette flow. *Phys. Rev. Lett.* **2013**, *110*, 204502. [CrossRef] [PubMed]

47. Sreenivasan, K.; Ramshankar, R. Transition intermittency in open flows, and intermittency routes to chaos. *Phys. D Nonlinear Phenom.* **1986**, *23*, 246–258. [CrossRef]

48. Avila, M.; Hof, B. Nature of laminar-turbulence intermittency in shear flows. *Phys. Rev. E* **2013**, *87*, 063012. [CrossRef] [PubMed]

49. Ishida, T.; Brethouwer, G.; Duguet, Y.; Tsukahara, T. Laminar-turbulent patterns with rough walls. *Phys. Rev. Fluids* **2017**, *2*, 073901. [CrossRef]

50. Tsukahara, T.; Tomioka, T.; Ishida, T.; Duguet, Y.; Brethouwer, G. Transverse turbulent bands in rough plane Couette flow. *J. Fluid Sci. Technol.* **2018**, *13*, JFST0019. [CrossRef]

<figure><img src="image_13.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).