Subcritical transition to turbulence in wallbounded flows: the case of plane Poiseuille flow
Abstract
In wallbounded flows, the laminar regime remain linearly stable up to large values of the Reynolds number while competing with nonlinear turbulent solutions issued from finite amplitude perturbations. The transition to turbulence of plane channel flow (plane Poiseuille flow) is more specifically considered via numerical simulations. Previous conflicting observations are reconciled by noting that the twodimensional directed percolation scenario expected for the decay of turbulence may be interrupted by a symmetrybreaking bifurcation favoring localized turbulent bands. At the other end of the transitional range, a preliminary study suggests that the laminarturbulent pattern leaves room to a featureless regime beyond a well defined threshold to be determined with precision.
1 Context
Closed systems such as Rayleigh–Bénard convection become turbulent according to a globally supercritical scenario starting from an intuitively simple linear instability. As the control parameter is increased, this primary instability is then followed by a small number of successive, secondary, tertiary,… instabilities, continuously replacing a bifurcating regime by a nearby bifurcated regime, rendering the whole process, at least in principle, accessible to weakly nonlinear perturbation theory [1]. Open flows with inflectional velocity profiles are prone to an inertial instability (Kelvin–Helmholtz) and, as observed in mixing layers, jets, or wakes, similarly follow a globally supercritical scenario [2].
This is in strong contrast with the case of wallbounded shear flows depicted in Fig. 1. With velocity profiles deprived from inflection points, unidirectional flows along walls, in tubes or plane channels, boundary layers, etc. remain inertially stable and may experience instability against Tollmien–Schichting waves, but only beyond some high linear threshold , where is a suitably defined Reynolds number [2]. can even be pushed to infinity, e.g. for flows in ducts of circular or square sections, or for plane Couette flow (PCF). Whether this instability is supercritical (a) or subcritical (b) [case of plane channel (Poiseuille) flow] is irrelevant since, in practice, much below , nontrivial states may exist in addition to the trivial laminar state owing to nonlinearities inherent in the advection term of Navier–Stokes equation (NSE). The flow can therefore “jump” to the nontrivial branch populated with chaotic solutions as a result of the amplification of finite amplitude perturbations, hence a direct, discontinuous, globally subcritical transition to turbulence, and the definition of a global stability threshold below which laminar flow is recovered from any initial condition after a possibly long transient. A general review can be found in [3].
Relying on the identification of essential features of the processes sustaining nontrivial states, a cyclical sequence of streamwise vortices inducing streaks by liftup and subsequent breakdown – Waleffe’s SSP [4], a large body of research on this transition has been focussing on the elucidation of the structure of the phase space in the framework of lowdimensional dynamical systems. The concept of minimal flow unit (MFU) introduced in [5] –minimal periodized domain in which nontrivial states can persist– was instrumental in such studies, ending in a picture where unstable coherent structures and manifolds attached to them are key elements [6]. Interesting outcomes of this approach are the interpretation of turbulence’s lifetime distributions in terms of chaotic transients and the discovery of edge states sitting on the boundary of the attraction basin of the laminar solution, from the neighborhood of which trajectories can be found to land on the turbulent attractor or else are visited during decay to laminar flow.
By nature unable to deal with spatial modulations in physical space the MFU approach had to be extended to treat the passage from attractor coexistence in phase space to what is actually observed: the coexistence in physical space of domains filled with one of the possible states, either trivial or nontrivial, and interfaces between them, which –by the way– is generic in extended systems experiencing a subcritical instability. Accordingly, localized edge states were found, e.g. [7], that can serve as germ for the expansion of turbulence above or show up at the latest stages of decay. This directly leads to the interpretation of the transition in terms of spatiotemporal intermittency promoted by Pomeau [8] who conjectured that this regime should decay following a scenario interpreted within the phasetransition framework as a transition in the universality class of a stochastic contamination process called directed percolation (DP). This nonequilibrium process describe the invasion of a uniformly quiescent absorbing regime, here obviously laminar flow, by an active state, here chaotic or turbulent, as the contamination probability increases. In this framework, the global stability threshold corresponds to the DP threshold and scaling properties are expected in its vicinity. This statisticalphysics viewpoint has received some support from experiments in a quasionedimensional cylindrical Couette flow (CCF) configuration [9] and by numerical simulations of a model of shear flow without walls in a quasitwodimensional geometry [10], among a few other cases.
Channel flow is one such additional case for which agreement has been found with 2DDP universality in a decay experiment from uniform turbulence [11]. This result is however contradicted by the observation of oblique localized turbulent bands (LTBs) in numerical [12] or laboratory experiments [13] below reported in [11]. Our numerical study, to be developed below [14] reveals how these two contradictory results can be reconciled by identifying a transition that preempts the DP scenario just above its supposed critical point.
As further increases, the spatiotemporally intermittent laminarturbulent distribution progressively disappears leaving a state of essentially uniform turbulence called featureless [15]. In PCF [16] and CCF at large circumferential aspect ratio [15, 16], as well as in channel flow [17], and a few other cases, the laminarturbulent alternation displays a wellorganized oblique turbulent band (or spiral) pattern. A threshold can sometimes be defined for the establishment of the featureless regime, e.g. in PCF [16, 18], as also suggested by our preliminary results for channel flow (see below).
Summarizing this long introduction, wallbounded flows generically experience a direct and wild transition transition to turbulence with different regimes as illustrated in Fig. 1: (i) nontrivial solution branch disconnected from laminar flow, to be studied using statisticalphysics tools in use for critical phenomena, with decay after long chaotic transients below ; (ii) above , regular pattern observable up to an upper threshold and featureless turbulent flow above ; (iii) at the moderate of interest, viscous effects strong enough to ensure the persistence of coherent structures, making tools from the theory of deterministic, dissipative, lowdimensional dynamical systems still relevant.
Remark also that laminar flow is still a possible linearly stable solution to the NSE all along the transitional range and finite amplitude germs are required to leave the trivial branch, which ordinarily is sufficient to establish the subcritical character of the transition. This character is however often expressed trough an expected discontinuous behavior of observables at the corresponding threshold, here . This qualitative expectation may be misleading since, quantitatively, the “distance” between the trivial and nontrivial branch may tend to zero in terms of turbulent fraction, as observed in several cases [9, 10, 11].
2 Plane channel flow
We now turn to our own results for channel flow obtained by direct numerical simulations of NSE using a program developed by one of us. Details can be found in [14]. The flow is driven by a constant body force . The centerline velocity of the laminar flow induced by this force and the halfdistance between the plates are used to turn the equations dimensionless, hence , with the kinematic viscosity of the fluid, and . A domain is considered with periodic boundary conditions in the streamwise and spanwise directions () and the usual noslip condition in the wallnormal direction . A spectral Fourier–Chebyshev–Fourier scheme is developed in a velocityvorticity formulation involving the wallnormal velocity and vorticity [19]. Most simulations have been performed for and and complementary studies not reported here in four times larger or smaller domains. The working resolution is with Fourier wavenumbers in each inplane direction and Chebyshev polynomials along . Aliasing is fully removed by evaluating nonlinear terms using a total of mode amplitudes. The simulations appear to be wellresolved in the range of of interest to the transition.
Overview of the transitional range in physical space.
At given below [20], besides the laminar solution, NSE possess solutions triggered by finite perturbations and belonging to the nontrivial branch mentioned in §1. Such solutions are associated to a mean streamwise velocity smaller than the corresponding laminar value (in units of ). Due to turbulence, speed , often called bulk velocity, is an empirical quantity that fluctuates and has to be averaged in time. Accordingly, we define to facilitate the comparison with other works that use the bulk velocity in the definition of the Reynolds number ( for laminar flow). The different flow regimes we have observed are depicted in Fig. 2. A preliminary simulation at , starting from a finite amplitude initial random condition, has provided us with a pattern that has next been evolved upon decreasing or increasing regularly (values up to have been considered). Solutions presented have reached a statistically steady state.
LTBs described in [12, 13] are here shown for , in the onesided propagation regime where all the bands go in the same direction. They are seen to recede and disappear for . At they propagate in both directions, forming a strongly intermittent loose discontinuous laminarturbulent network, with turbulent arms broken by laminar gaps. At , laminar gaps have disappeared and the network is now continuous but still loose and intermittent. When is further increased this crisscrossed pattern gets tight and more regular with conspicuous domains of one or the other orientation separated by grain boundaries, here at . This pattern next fades away, being hardly visible at . A close parallel can be drawn with the case of CCF (and PCF) examined by the Saclay group [16].
The lower transitional range.
Below , the laminarturbulent alternation is sufficiently marked for the definition of a turbulent fraction to make sense. This turbulent fraction was determined using a twolevel momentpreserving thresholding method [20] that automatically places the cutoff between laminar and turbulent local states so as to reduce the distribution of the graylevels to an optimal blackandwhite distribution.
Figure 3 (left) displays the variation of as a function of .
The line corresponds to a fit of the data against the function
At decreasing , the symmetrybreaking transition is due to the rapidly decreasing probability of lateral branching compared to splitting along the LTBs. The two orientations appear with essentially equal weights and the rate of lateral branching is comparable to that of parallel splitting in the symmetrical regime ( in Fig. 2). But this rate rapidly decreases in the green zones appearing in Fig. 3 and becomes negligible for in the yellow zones. An observable measuring the lack of symmetry has been designed as plotted in Fig. 3 (right). A simple phenomenological model has been developed for this symmetry restoration [14], supporting a linear variation of the asymmetry amplitude that helped us extrapolate it to zero and define the corresponding threshold at .
The upper transitional range.
For , the spatiotemporally intermittent laminarturbulent alternation is less marked rendering the analysis in terms of turbulence fraction less appropriate while the pattern becomes both more steady and more regular, suggesting the recourse to Fourier analysis. The transverse perturbation energy field is considered here and series of Fourier spectra are recorded at statistically steady state. Means and standard deviations of time series of the intensity of each mode with , i.e. , are computed for 50 snapshots separated by , hence essentially independent.
Figure 4 (top) display the means of these spectra . At large inside the supposedly featureless regime, wide bumps are found, corresponding to highly fluctuating, low amplitude, modulations, less and less pronounced as increases (). In contrast, when the pattern is well established the spectra are dominated by sharp peaks (). The variation of the maximum as a function of is displayed using linlog scale in Fig. 4 (bottom, left). It shows that the inflection at () corresponds to a rapid increase, interpreted as a sudden condensation of the weak turbulenceintensity modulations into a genuine pattern. Figure 4 (bottom, left) displays the position of these maxima in the wavevector plane as a function of . It is seen that for , the twocomponents of the wave vectors increase roughly linearly with , i.e. a regular wavelength decrease from at to at . They however maintain a nearly constant ratio , hence turbulent bands essentially making an angle with the streamwise direction all along this range. For , the wavelength of the turbulenceintensity fluctuations continues to decrease slightly as they become weaker and weaker with increasing .
3 Discussion/Conclusion
The transition to turbulence in wallbounded flows, briefly reviewed in §1, remains a good testing ground for the theory of nonequilibrium phase transitions. In §2, we illustrated the case of plane channel flow driven by a constant bulk force using our own simulations. Results presented in the recent literature, among others [11, 12, 13, 17], were recovered. In particular, the relevance of the 2DDP scenario for turbulence decay [11] was confirmed but reconciled with other results showing the presence of localized turbulent states below the DPthreshold [12, 13]. Within our protocol, a symmetrybreaking bifurcation due to the rapid decrease of the rate of lateral branching in the loose spatiotemporally intermittent network regime forces the system to a new regime preempting the critical behavior inherent in the DP framework [14]. Ongoing work is devoted to improve the statistical study of the flow in the upper transitional range. The physical mechanisms underlying the branching processes at increasing near and those responsible for the patterning observed at decreasing at the still putative threshold remain to be elucidated.
Acknowledgments.
This work is supported by JSPS KAKENHI Grant Number JP17K14588. Simulations are performed using NIFS’s FUJITSU FX100 “Plasma Simulator” (contract NIFS16KNSS083). We would like to thank Y. Duguet (LIMSI) for interesting discussions about the problem.
Footnotes
 The control parameter is rather than , yielding an excellent fit from up to .
References
 P. Manneville, Instabilities, Chaos and Turbulence, 2nd edition, Imperial College Press (2010).
 P. Huerre & M. Rossi, Hydrodynamic instabilities in open flows, in C. Godrèche & P. Manneville, Hydrodynamics and Nonlinear Instabilities, Cambridge University Press (1998).

P. Manneville, Mechanical Engineering Reviews 3, 15 (2016)
[https://doi.org/10.1299/mer.1500684].  F. Waleffe, Physics of Fluids 9, 883 (1997).
 J. Jiménez & P. Moin, J. Fluid Mech. 225, 213 (1991).
 G. Kawahara et al., Annu. Rev. Fluid Mech. 44, 203 (2011).
 Y. Duguet et al., Phys. Fluids 21, 111701 (2009).
 Y. Pomeau, Physica D 23, 3 (1986).
 G. Lemoult et al., Nature Physics, 12, 254 (2016).
 M. Chantry et al., J. Fluid Mech. 824, R1 (2017).
 M. Sano & K. Tamai, Nature Physics 12, 249 (2016).
 X. Xiong et al., Phys. Fluids 27, 041702 (2015).

C.S. Paranjape et al.,
http://www.delegia.com/app/data/8684/Abstract/29766/ETC16_Paranjape.pdf  M. Shimizu & P. Manneville, arXiv:1808.06479 [physics.fludyn].
 C.D. Andereck et al., J. Fluid Mech. 164, 155 (1986).
 A. Prigent et al., Physica D 174, 100 (2003).
 T. Tsukahara et al., 4th TSFP Symposium, Williamsburg, Proceedings p. 935ff (2005).
 L.S. Tuckerman & D. Barkley, Phys. Fluids 23, 041301 (2011).
 L. Quartapelle, Numerical Solution of the Incompressible Navier–Stokes Equations, ISNM vol. 113, Birkhäuser (1993).
 S. Orszag, J. Fluid Mech. 50, 689 (1971).
 W.H. Tsai, Computer vision, graphics, and image processing 29, 377 (1985).