Traveling waves in pipe flow

Following marginal stability manifolds in quasilinear dynamical reductions of multiscale flows in two space dimensions

We derive a two-dimensional (2D) extension of a recently developed formalism for slow-fast quasilinear (QL) systems subject to fast instabilities. The emergent dynamics of these systems is characterized by a slow evolution of (suitably defined) mean fields coupled to marginally stable, fast fluctuation fields. By exploiting this scale separation, an efficient hybrid fast-eigenvalue/slow-initial-value solution algorithm can be developed in which the amplitude of the fast fluctuations is slaved to the slowly evolving mean fields to ensure marginal stability-and temporal scale separation-is maintained. For 2D systems, the fluctuation eigenfunctions are labeled by their Fourier wave numbers characterizing spatial variability in that extended spatial direction, and the marginal mode(s) must coincide with the fastest-growing mode(s) over all admissible Fourier wave numbers. Here we derive an ordinary differential equation governing the slow evolution of the wave number of the fastest-growing fluctuation mode that simultaneously must be slaved to the mean dynamics to ensure the mode has zero growth rate. We illustrate the procedure in the context of a 2D model partial differential equation that shares certain attributes with the equations governing strongly stratified shear flows and other strongly constrained forms of geophysical turbulence in extreme parameter regimes. The slaved evolution follows one or more marginal stability manifolds, which constitute select state-space structures that are not invariant under the full flow dynamics yet capture quasicoherent structures in physical space in a manner analogous to invariant solutions identified in, e.g., transitionally turbulent shear flows. Accordingly, we propose that marginal stability manifolds are central organizing structures in a dynamical systems description of certain classes of multiscale flows in which scale separation justifies a QL approximation of the dynamics.

High-Speed Schlieren Imaging of Vapor Formation in Electrospray Plume

Previous mechanistic descriptions of electrosprays mostly focused on the dynamics of Taylor cones, initial droplets, and progeny droplets. However, vapor formation during droplet desolvation in an electrospray plume has not been discussed to a great extent. Here, we implement a double-pass on-axis schlieren high-speed imaging system to observe generation and propagation of vapors in an offline electrospray source under different conditions. Switching between turbulent and laminar vapor flow was observed for all of the scanned conditions, which may be attributed to randomly occurring disturbances in the sample flow inside the electrospray emitter. Calculation of mean vapor flow velocity and analysis of vapor flow patterns were performed using in-house developed image processing programs. Experiments performed at different electrospray voltages (0-6 kV), solvent flow rates (100-600 μL min-1), and methanol concentrations (50-100%), indicate only a weak dependency between electrospray voltage and mean vapor velocity, implying that the vapor is mostly neutral; thus, the vapor is not accelerated by electric field. On the other hand, electrospraying solutions of analytes (with mass 151 Da or 12 kDa) did not remarkably increase the overall vapor flow velocity. The source of vapor's velocity is attributed to the inertia of the electrospray droplets. Although there are some differences between a modern electrospray ionization (ESI) setup and the setup used in our experiment (e.g., using a higher flow rate and larger emitter), we believe the findings of our study can be projected to a modern ESI setup.

Exploring the free-energy landscape of a rotating superfluid

The equilibrium state of a superfluid in a rotating cylindrical vessel is a vortex crystal-an array of vortex lines, which is stationary in the rotating frame. Experimental realizations of this behavior typically show a sequence of transient states before the free-energy-minimizing configuration is reached. Motivated by these observations, we construct a new method for a systematic exploration of the free-energy landscape via gradient-based optimization of a scalar loss function. Our approach is inspired by the pioneering numerical work of Campbell and Ziff [Phys. Rev. B. 20, 1886 (1979)] and makes use of automatic differentiation, which crucially allows us to include entire solution trajectories in the loss. We first use the method to converge thousands of low free-energy relative equilibria in the unbounded domain for vortex numbers in the range 10≤N≤30, which reveals an extremely dense set of mostly saddle-like solutions. As part of this search, we discover new continuous families of relative equilibria, which are often global minimizers of free energy. These continuous families all consist of crystals arranged in a double-ring configuration, and we assess which state from the family is most likely to be observed experimentally by computing energy-minimizing pathways from nearby local minima-identifying a common entry point into the family. The continuous families become discrete sets of equal-energy solutions when the wall is introduced in the problem. Finally, we develop an approach to compute homoclinic orbits and use it to examine the dynamics in the vicinity of the minimizing state by converging connections for low-energy saddles.

Observing a dynamical skeleton of turbulence in Taylor-Couette flow experiments

Recent work shows that recurrent solutions of the equations governing fluid flow play an important role in structuring the dynamics of turbulence. Here, an improved version of an earlier method (Krygier et al. 2021 J. Fluid. Mech. 923, A7 and Crowley et al. 2022 Proc. Natl Acad. Sci. USA 119, e2120665119) is used for detecting and analyzing intervals of time when turbulence 'shadows' (spatially and temporally mimics) recurrent solutions in both numerical simulations and laboratory experiments. We find that all the recurrent solutions shadowed in numerics are also shadowed in experiment, and the corresponding statistics of shadowing agree. Our results set the stage for experimentally grounded dynamical descriptions of turbulence in a variety of wall-bounded shear flows, enabling applications to forecasting and control. This article is part of the theme issue 'Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper (part 1)'.

Turbulence tracks recurrent solutions

Despite a long and rich history of scientific investigation, fluid turbulence remains one of the most challenging problems in science and engineering. One of the key outstanding questions concerns the role of coherent structures that describe frequently observed patterns embedded in turbulence. It has been suggested, but not proved, that coherent structures correspond to unstable, recurrent solutions of the governing equation of fluid dynamics. Here, we present experimental and numerical evidence that three-dimensional turbulent flow tracks, episodically but repeatedly, the spatial and temporal structure of multiple such solutions. Our results provide compelling evidence that coherent structures, grounded in the governing equations, can be harnessed to predict how turbulent flows evolve.

Continuing invariant solutions towards the turbulent flow

A new mathematical framework is proposed for characterizing the coherent motion of fluctuations around a mean turbulent channel flow. We search for statistically invariant coherent solutions of the unsteady Reynolds-averaged Navier-Stokes equations written in a perturbative form with respect to the turbulent mean flow, using a suitable approximation of the Reynolds stress tensor. This is achieved by setting up a continuation procedure of known solutions of the perturbative Navier-Stokes equations, based on the continuous increase of the turbulent eddy viscosity towards its turbulent value. The recovered solutions, being sustained only in the presence of the Reynolds stress tensor, are representative of the statistically coherent motion of turbulent flows. For small friction Reynolds number and/or domain size, the statistically invariant motion is almost identical to the corresponding invariant solution of the Navier-Stokes equations. Whereas, for sufficiently large friction number and/or domain size, it considerably departs from the starting invariant solution of the Navier-Stokes equations, presenting spatial structures, main wavelengths and scaling very close to those characterizing both large- and small-scale motion of turbulent channel flows. This article is part of the theme issue 'Mathematical problems in physical fluid dynamics (part 2)'.

Dynamical landscape of transitional pipe flow

The transition to turbulence in pipes is characterized by a coexistence of laminar and turbulent states. At the lower end of the transition, localized turbulent pulses, called puffs, can be excited. Puffs can decay when rare fluctuations drive them close to an edge state lying at the phase-space boundary with laminar flow. At higher Reynolds numbers, homogeneous turbulence can be sustained, and dominates over laminar flow. Here we complete this landscape of localized states, placing it within a unified bifurcation picture. We demonstrate our claims within the Barkley model, and motivate them generally. Specifically, we suggest the existence of an antipuff and a gap-edge-states which mirror the puff and related edge state. Previously observed laminar gaps forming within homogeneous turbulence are then naturally identified as antipuffs nucleating and decaying through the gap edge. We also discuss alternatives to the suggested bifurcation diagram, which could be relevant for wall-bounded flows other than straight pipes.

Coarse Graining the State Space of a Turbulent Flow Using Periodic Orbits

We show that turbulent dynamics that arise in simulations of the three-dimensional Navier-Stokes equations in a triply periodic domain under sinusoidal forcing can be described as transient visits to the neighborhoods of unstable time-periodic solutions. Based on this description, we reduce the original system with more than 10^{5} degrees of freedom to a 17-node Markov chain where each node corresponds to the neighborhood of a periodic orbit. The model accurately reproduces long-term averages of the system's observables as weighted sums over the periodic orbits.

Susceptibility of transient chimera states

Chaotic dynamics of a dynamical system is not necessarily persistent. If there is (without any active intervention from outside) a transition towards a (possibly nonchaotic) attractor, this phenomenon is called transient chaos, which can be observed in a variety of systems, e.g., in chemical reactions, population dynamics, neuronal activity, or cardiac dynamics. Also, chimera states, which show coherent and incoherent dynamics in spatially distinct regions of the system, are often chaotic transients. In many practical cases, the control of the chaotic dynamics (either the termination or the preservation of the chaotic dynamics) is desired. Although the self-termination typically occurs quite abruptly and can so far in general not be properly predicted, previous studies showed that in many systems a 'terminal transient phase" (TTP) prior to the self-termination existed, where the system was less susceptible against small but finite perturbations in different directions in state space. In this study, we show that, in the specific case of chimera states, these susceptible directions can be related to the structure of the chimera, which we divide into the coherent part, the incoherent part and the boundary in between. That means, in practice, if self-termination is close we can identify the direction of perturbation which is likely to maintain the chaotic dynamics (the chimera state). This finding improves the general understanding of the state space structure during the TTP, and could contribute also to practical applications like future control strategies of epileptic seizures which have been recently related to the collapse of chimera states.

Terminating transient chaos in spatially extended systems

In many real-life systems, transient chaotic dynamics plays a major role. For instance, the chaotic spiral or scroll wave dynamics of electrical excitation waves during life-threatening cardiac arrhythmias can terminate by itself. Epileptic seizures have recently been related to the collapse of transient chimera states. Controlling chaotic transients, either by maintaining the chaotic dynamics or by terminating it as quickly as possible, is often desired and sometimes even vital (as in the case of cardiac arrhythmias). We discuss in this study that the difference of the underlying structures in state space between a chaotic attractor (persistent chaos) and a chaotic saddle (transient chaos) may have significant implications for efficient control strategies in real life systems. In particular, we demonstrate that in the latter case, chaotic dynamics in spatially extended systems can be terminated via a relatively low number of (spatially and temporally) localized perturbations. We demonstrate as a proof of principle that control and targeting of high-dimensional systems exhibiting transient chaos can be achieved with exceptionally small interactions with the system. This insight may impact future control strategies in real-life systems like cardiac arrhythmias.

Basis for finding exact coherent states

One of the outstanding problems in the dynamical systems approach to turbulence is to find a sufficient number of invariant solutions to characterize the underlying dynamics of turbulence [Annu. Rev. Fluid Mech. 44, 203 (2012)10.1146/annurev-fluid-120710-101228]. As a practical matter, the solutions can be difficult to find. To improve this situation, we show how to find periodic orbits and equilibria in plane Couette flow by projecting pseudorecurrent segments of turbulent trajectories onto the left-singular vectors of the Navier-Stokes equations linearized about the relevant mean flow (resolvent modes). The projections are, subsequently, used to initiate Newton-Krylov-hookstep searches, and new (relative) periodic orbits and equilibria are discovered. We call the process project-then-search and validate the process by first applying it to previously known fixed point and periodic solutions. Along the way, we find new branches of equilibria, which include bifurcations from previously known branches, and new periodic orbits that closely shadow turbulent trajectories in state space.

Heteroclinic and homoclinic connections in a Kolmogorov-like flow

Recent studies suggest that unstable recurrent solutions of the Navier-Stokes equation provide new insights into dynamics of turbulent flows. In this study, we compute an extensive network of dynamical connections between such solutions in a weakly turbulent quasi-two-dimensional Kolmogorov flow that lies in the inversion-symmetric subspace. In particular, we find numerous isolated heteroclinic connections between different types of solutions-equilibria, periodic, and quasiperiodic orbits-as well as continua of connections forming higher-dimensional connecting manifolds. We also compute a homoclinic connection of a periodic orbit and provide strong evidence that the associated homoclinic tangle forms the chaotic repeller that underpins transient turbulence in the symmetric subspace.

Spontaneous termination of chaotic spiral wave dynamics in human cardiac ion channel models

Chaotic spiral or scroll wave dynamics can be found in diverse systems. In cardiac dynamics, spiral or scroll waves of electrical excitation determine the dynamics during life-threatening arrhythmias like ventricular fibrillation. In numerical studies it was found that chaotic episodes of spiral and scroll waves can be transient, thus they terminate spontaneously. We show in this study that this behavior can also be observed using models which describe the ion channel dynamics of human cardiomyocytes (Bueno-Orovio-Cherry-Fenton model and the Ten Tusscher-Noble-Noble-Panfilov model). For both models we find that the average lifetime of the chaotic transients grows exponentially with the system size. With this behavior, we classify the systems into the group of type-II supertransients. We observe a significant difference of the breakup behavior between the models, which results in a distinct dynamics during the final phase just before the termination. The observation of a (temporally) stable single-spiral state affects the prevailing description of the dynamics of type-II supertransients as being “quasi-stationary” and also the feasibility of predicting the spontaneous termination of the spiral wave dynamics. In the long term, the relation between the breakup behavior of spiral waves and properties of chaotic transients like predictability or average transient lifetime may contribute to an improved understanding and classification of cardiac arrhythmias.

Turbulence Appearance and Nonappearance in Thin Fluid Layers

Flows in fluid layers are ubiquitous in industry, geophysics, and astrophysics. Large-scale flows in thin layers can be considered two dimensional with bottom friction added. Here we find that the properties of such flows depend dramatically on the way they are driven. We argue that a wall-driven (Couette) flow cannot sustain turbulence, no matter how small the viscosity and friction. Direct numerical simulations (DNSs) up to the Reynolds number Re=10^{6} confirm that all perturbations die in a plane Couette flow. On the contrary, for sufficiently small viscosity and friction, perturbations destroy the pressure-driven laminar (Poiseuille) flow. What appears instead is a traveling wave in the form of a jet slithering between wall vortices. For 5×10^{3}<Re<3×10^{4}, the mean flow in most cases has remarkably simple structure: the jet is sinusoidal with a parabolic velocity profile, and vorticity is constant inside vortices, while the fluctuations are small. At higher Re, strong fluctuations appear, yet the mean traveling wave survives. Considering the momentum flux barrier in such a flow, we derive a new scaling law for the Re dependence of the friction factor and confirm it by DNS.

Unstable equilibria and invariant manifolds in quasi-two-dimensional Kolmogorov-like flow

Recent studies suggest that unstable, nonchaotic solutions of the Navier-Stokes equation may provide deep insights into fluid turbulence. In this article, we present a combined experimental and numerical study exploring the dynamical role of unstable equilibrium solutions and their invariant manifolds in a weakly turbulent, electromagnetically driven, shallow fluid layer. Identifying instants when turbulent evolution slows down, we compute 31 unstable equilibria of a realistic two-dimensional model of the flow. We establish the dynamical relevance of these unstable equilibria by showing that they are closely visited by the turbulent flow. We also establish the dynamical relevance of unstable manifolds by verifying that they are shadowed by turbulent trajectories departing from the neighborhoods of unstable equilibria over large distances in state space.

Terminal Transient Phase of Chaotic Transients

Transient chaos in spatially extended systems can be characterized by the length of the transient phase, which typically grows quickly with the system size (supertransients). For a large class of these systems, the chaotic phase terminates abruptly, without any obvious precursors in commonly used observables. Here we investigate transient spatiotemporal chaos in two different models of this class. By probing the state space using perturbed trajectories we show the existence of a "terminal transient phase," which occurs prior to the abrupt collapse of chaotic dynamics. During this phase the impact of perturbations is significantly different from the earlier transient and particular patterns of (non)susceptible regions in state space occur close to the chaotic trajectories. We therefore hypothesize that even without perturbations proper precursors for the collapse of chaotic transients exist, which might be highly relevant for coping with spatiotemporal chaos in cardiac arrhythmias or brain functionality, for example.

Forecasting Fluid Flows Using the Geometry of Turbulence

The existence and dynamical role of particular unstable solutions (exact coherent structures) of the Navier-Stokes equation is revealed in laboratory studies of weak turbulence in a thin, electromagnetically driven fluid layer. We find that the dynamics exhibit clear signatures of numerous unstable equilibrium solutions, which are computed using a combination of flow measurements from the experiment and fully resolved numerical simulations. We demonstrate the dynamical importance of these solutions by showing that turbulent flows visit their state space neighborhoods repeatedly. Furthermore, we find that the unstable manifold associated with one such unstable equilibrium predicts the evolution of turbulent flow in both experiment and simulation for a considerable period of time.

The relationship between free-stream coherent structures and near-wall streaks at high Reynolds numbers

The present work is based on our recent discovery of a new class of exact coherent structures generated near the edge of quite general boundary layer flows. The structures are referred to as free-stream coherent structures and were found using a large Reynolds number asymptotic approach to describe equilibrium solutions of the Navier-Stokes equations. In this paper, first we present results for a new family of free-stream coherent structures existing at relatively large wavenumbers. The new results are consistent with our earlier theoretical result that such structures can generate larger amplitude wall streaks if and only if the local spanwise wavenumber is sufficiently small. In a Blasius boundary layer, the local wavenumber increases in the streamwise direction so the wall streaks can typically exist only over a finite interval. However, here it is shown that they can interact with wall curvature to produce exponentially growing Görtler vortices through the receptivity process by a novel nonparallel mechanism. The theoretical predictions found are confirmed by a hybrid numerical approach. In contrast with previous receptivity investigations, it is shown that the amplitude of the induced vortex is larger than the structures in the free-stream which generate it.This article is part of the themed issue 'Toward the development of high-fidelity models of wall turbulence at large Reynolds number'.

Marginal turbulent state of viscoelastic fluids: A polymer drag reduction perspective

The laminar-turbulent (LT) transition of dilute polymer solutions is of great interest not only for the complex transition dynamics itself, but also for its potential link to the maximum drag reduction (MDR) phenomenon. We present an in-depth investigation of the edge state (ES), an asymptotic solution on the LT boundary, in viscoelastic channel flow. For given Re and simulation domain size, mean flow statistics of the ES do not vary with the introduction of polymers, proving that there is a region of turbulent states not susceptible to polymer drag reduction effects. The dynamics of the ES features low-frequency fluctuations and in the longer domains we studied it is nearly periodic with regular bursts of turbulent activities separated by extended quiescent periods. Its flow field is dominated by elongated vortices and streaks, with very weak extensional and rotational flow motions. Polymer stretching is almost exclusively contributed by the mean shear and polymer-turbulence interaction is minimal. Flow structures and the kinematics of the ES match hibernating turbulence, an MDR-like phase intermittently occurring in turbulent dynamics. Its observation now seems to result from recurrent visits to certain parts of the ES. The ES offers explanations for the existence and universality of MDR, the quantitative magnitude of which, however, still remains unsolved.

Low-dimensional representations of exact coherent states of the Navier-Stokes equations from the resolvent model of wall turbulence

We report that many exact invariant solutions of the Navier-Stokes equations for both pipe and channel flows are well represented by just a few modes of the model of McKeon and Sharma [J. Fluid Mech. 658, 336 (2010)]. This model provides modes that act as a basis to decompose the velocity field, ordered by their amplitude of response to forcing arising from the interaction between scales. The model was originally derived from the Navier-Stokes equations to represent turbulent flows and has been used to explain coherent structure and to predict turbulent statistics. This establishes a surprising new link between the two distinct approaches to understanding turbulence.

Sudden Relaminarization and Lifetimes in Forced Isotropic Turbulence

We demonstrate an unexpected connection between isotropic turbulence and wall-bounded shear flows. We perform direct numerical simulations of isotropic turbulence forced at large scales at moderate Reynolds numbers and observe sudden transitions from a chaotic dynamics to a spatially simple flow, analogous to the laminar state in wall bounded shear flows. We find that the survival probabilities of turbulence are exponential and the typical lifetimes increase superexponentially with the Reynolds number. Our results suggest that both isotropic turbulence and wall-bounded shear flows qualitatively share the same phase-space dynamics.

Effect of plate permeability on nonlinear stability of the asymptotic suction boundary layer

The nonlinear stability of the asymptotic suction boundary layer is studied numerically, searching for finite-amplitude solutions that bifurcate from the laminar flow state. By changing the boundary conditions for disturbances at the plate from the classical no-slip condition to more physically sound ones, the stability characteristics of the flow may change radically, both for the linearized as well as the nonlinear problem. The wall boundary condition takes into account the permeability K̂ of the plate; for very low permeability, it is acceptable to impose the classical boundary condition (K̂=0). This leads to a Reynolds number of approximately Re(c)=54400 for the onset of linearly unstable waves, and close to Re(g)=3200 for the emergence of nonlinear solutions [F. A. Milinazzo and P. G. Saffman, J. Fluid Mech. 160, 281 (1985); J. H. M. Fransson, Ph.D. thesis, Royal Institute of Technology, KTH, Sweden, 2003]. However, for larger values of the plate's permeability, the lower limit for the existence of linear and nonlinear solutions shifts to significantly lower Reynolds numbers. For the largest permeability studied here, the limit values of the Reynolds numbers reduce down to Re(c)=796 and Re(g)=294. For all cases studied, the solutions bifurcate subcritically toward lower Re, and this leads to the conjecture that they may be involved in the very first stages of a transition scenario similar to the classical route of the Blasius boundary layer initiated by Tollmien-Schlichting (TS) waves. The stability of these nonlinear solutions is also investigated, showing a low-frequency main unstable mode whose growth rate decreases with increasing permeability and with the Reynolds number, following a power law Re(-ρ), where the value of ρ depends on the permeability coefficient K̂. The nonlinear dynamics of the flow in the vicinity of the computed finite-amplitude solutions is finally investigated by direct numerical simulations, providing a viable scenario for subcritical transition due to TS waves.

Osborne Reynolds pipe flow: Direct simulation from laminar through gradual transition to fully developed turbulence

The precise dynamics of breakdown in pipe transition is a century-old unresolved problem in fluid mechanics. We demonstrate that the abruptness and mysteriousness attributed to the Osborne Reynolds pipe transition can be partially resolved with a spatially developing direct simulation that carries weakly but finitely perturbed laminar inflow through gradual rather than abrupt transition arriving at the fully developed turbulent state. Our results with this approach show during transition the energy norms of such inlet perturbations grow exponentially rather than algebraically with axial distance. When inlet disturbance is located in the core region, helical vortex filaments evolve into large-scale reverse hairpin vortices. The interaction of these reverse hairpins among themselves or with the near-wall flow when they descend to the surface from the core produces small-scale hairpin packets, which leads to breakdown. When inlet disturbance is near the wall, certain quasi-spanwise structure is stretched into a Lambda vortex, and develops into a large-scale hairpin vortex. Small-scale hairpin packets emerge near the tip region of the large-scale hairpin vortex, and subsequently grow into a turbulent spot, which is itself a local concentration of small-scale hairpin vortices. This vortex dynamics is broadly analogous to that in the boundary layer bypass transition and in the secondary instability and breakdown stage of natural transition, suggesting the possibility of a partial unification. Under parabolic base flow the friction factor overshoots Moody's correlation. Plug base flow requires stronger inlet disturbance for transition. Accuracy of the results is demonstrated by comparing with analytical solutions before breakdown, and with fully developed turbulence measurements after the completion of transition.

Coherent structures in wall-bounded turbulence

The inherent difficulty of understanding turbulence has led to researchers attacking the topic in many different ways over the years of turbulence research. Some approaches have been more successful than others, but most only deal with part of the problem. One approach that has seen reasonable success (or at least popularity) is that of attempting to deconstruct the complex and disorganised turbulent flow field into to a set of motions that are in some way organised. These motions are generally called "coherent structures". There are several strands to this approach, from identifying the coherent structures within the flow, defining their characteristics, explaining how they are created, sustained and destroyed, to utilising their features to model the turbulent flow. This review considers research on coherent structures in wall-bounded turbulent flows: a class of flow which is extremely interesting to many scientists (mainly, but not exclusively, physicists and engineers) due to their prevalence in nature, industry and everyday life. This area has seen a lot of activity, particularly in recent years, much of which has been driven by advances in experimental and computational techniques. However, several ideas, developed many years ago based on flow visualisation and intuition, are still both informative and relevant. Indeed, much of the more recent research is firmly indebted to some of the early pioneers of the coherent structures approach. Therefore, in this review, selected historical research is discussed along with the more contemporary advances in an attempt to provide the reader with a good overview of how the field has developed and to highlight the perspicacity of some of the early researchers, as well as providing an overview of our current understanding of the role of coherent structures in wall-bounded turbulent flows.

Reduced description of exact coherent states in parallel shear flows

A reduced description of exact coherent structures in the transition regime of plane parallel shear flows is developed, based on the Reynolds number scaling of streamwise-averaged (mean) and streamwise-varying (fluctuation) velocities observed in numerical simulations. The resulting system is characterized by an effective unit Reynolds number mean equation coupled to linear equations for the fluctuations, regularized by formally higher-order diffusion. Stationary coherent states are computed by solving the resulting equations simultaneously using a robust numerical algorithm developed for this purpose. The algorithm determines self-consistently the amplitude of the fluctuations for which the associated mean flow is just such that the fluctuations neither grow nor decay. The procedure is used to compute exact coherent states of a flow introduced by Drazin and Reid [Hydrodynamic Stability (Cambridge University Press, Cambridge, UK, 1981)] and studied by Waleffe [Phys. Fluids 9, 883 (1997)]: a linearly stable, plane parallel shear flow confined between stationary stress-free walls and driven by a sinusoidal body force. Numerical continuation of the lower-branch states to lower Reynolds numbers reveals the presence of a saddle node; the saddle node allows access to upper-branch states that are, like the lower-branch states, self-consistently described by the reduced equations. Both lower- and upper-branch states are characterized in detail.

Localization in a spanwise-extended model of plane Couette flow

We consider a nine-partial-differential-equation (1-space and 1-time) model of plane Couette flow in which the degrees of freedom are severely restricted in the streamwise and cross-stream directions to study spanwise localization in detail. Of the many steady Eckhaus (spanwise modulational) instabilities identified of global steady states, none lead to a localized state. Spatially localized, time-periodic solutions were found instead, which arise in saddle node bifurcations in the Reynolds number. These solutions appear global (domain filling) in narrow (small spanwise) domains yet can be smoothly continued out to fully spanwise-localized states in very wide domains. This smooth localization behavior, which has also been seen in fully resolved duct flow (S. Okino, Ph.D. thesis, Kyoto University, Kyoto, 2011), indicates that an apparently global flow structure does not have to suffer a modulational instability to localize in wide domains.

Distinct organizational States of fully developed turbulent pipe flow

Organizational states of turbulence are identified through novel analysis of large scale pipe flow experiments at a Reynolds number of 35 000. The distinct states are revealed by an azimuthal decomposition of the two-point spatial correlation of the streamwise velocity fluctuation. States with dominant azimuthal wave numbers corresponding to k_{θ}=2,3,4,5,6 are discovered and their structure revealed as a series of alternately rotating quasistreamwise vortices. Such organizational states are highly reminiscent of the nonlinear traveling wave solutions previously identified at Reynolds numbers an order of magnitude lower.

Canonical exact coherent structures embedded in high Reynolds number flows

The applications and implications of two recently addressed asymptotic descriptions of exact coherent structures in shear flows are discussed. The first type of asymptotic framework to be discussed was introduced in a series of papers by Hall & Smith in the 1990s and was referred to as vortex-wave interaction theory (VWI). New results are given here for the canonical VWI problem in an infinite region; the results confirm and extend the results for the infinite problem inferred the recent VWI computation of plane Couette flow. The results given define for the first time exact coherent structures in unbounded flows. The second type of canonical structure described here is that recently found for asymptomatic suction boundary layer and corresponds to freestream coherent structures (FCS), in boundary layer flows. Here, it is shown that the FCS can also occur in flows such as Burgers vortex sheet. It is concluded that both canonical problems can be locally embedded in general shear flows and thus have widespread applicability.

Genesis of streamwise-localized solutions from globally periodic traveling waves in pipe flow

The aim in the dynamical systems approach to transitional turbulence is to construct a scaffold in phase space for the dynamics using simple invariant sets (exact solutions) and their stable and unstable manifolds. In large (realistic) domains where turbulence can coexist with laminar flow, this requires identifying exact localized solutions. In wall-bounded shear flows, the first of these has recently been found in pipe flow, but questions remain as to how they are connected to the many known streamwise-periodic solutions. Here we demonstrate that the origin of the first localized solution is in a modulational symmetry-breaking Hopf bifurcation from a known global traveling wave that has twofold rotational symmetry about the pipe axis. Similar behavior is found for a global wave of threefold rotational symmetry, this time leading to two localized relative periodic orbits. The clear implication is that many global solutions should be expected to lead to more realistic localized counterparts through such bifurcations, which provides a constructive route for their generation.

Long-wavelength instability of coherent structures in plane Couette flow

We study the stability of coherent structures in plane Couette flow against long-wavelength perturbations in wide domains that cover several pairs of coherent structures. For one and two pairs of vortices, the states retain the stability properties of the small domains, but for three pairs new unstable modes are found. They are shown to be connected to bifurcations that break the translational symmetry and drive the coherent structures from the spanwise extended state to a modulated one that is a precursor to spanwise localized states. Tracking the stability of the orbits as functions of the spanwise wave length reveals a rich variety of additional bifurcations.

Unstable flow structures in the Blasius boundary layer

Finite amplitude coherent structures with a reflection symmetry in the spanwise direction of a parallel boundary layer flow are reported together with a preliminary analysis of their stability. The search for the solutions is based on the self-sustaining process originally described by Waleffe (Phys. Fluids 9, 883 (1997)). This requires adding a body force to the Navier-Stokes equations; to locate a relevant nonlinear solution it is necessary to perform a continuation in the nonlinear regime and parameter space in order to render the body force of vanishing amplitude. Some states computed display a spanwise spacing between streaks of the same length scale as turbulence flow structures observed in experiments (S.K. Robinson, Ann. Rev. Fluid Mech. 23, 601 (1991)), and are found to be situated within the buffer layer. The exact coherent structures are unstable to small amplitude perturbations and thus may be part of a set of unstable nonlinear states of possible use to describe the turbulent transition. The nonlinear solutions survive down to a displacement thickness Reynolds number Re * = 496 , displaying a 4-vortex structure and an amplitude of the streamwise root-mean-square velocity of 6% scaled with the free-stream velocity. At this Re* the exact coherent structure bifurcates supercritically and this is the point where the laminar Blasius flow starts to cohabit the phase space with alternative simple exact solutions of the Navier-Stokes equations.

Complexity of localised coherent structures in a boundary-layer flow

We study numerically transitional coherent structures in a boundary-layer flow with homogeneous suction at the wall (the so-called asymptotic suction boundary layer ASBL). The dynamics restricted to the laminar-turbulent separatrix is investigated in a spanwise-extended domain that allows for robust localisation of all edge states. We work at fixed Reynolds number and study the edge states as a function of the streamwise period. We demonstrate the complex spatio-temporal dynamics of these localised states, which exhibits multistability and undergoes complex bifurcations leading from periodic to chaotic regimes. It is argued that in all regimes the dynamics restricted to the edge is essentially low-dimensional and non-extensive.

Instability of streaks in pipe flow of shear-thinning fluids

This study is motivated by recent experimental results dealing with the transition to turbulence in a pipe flow of shear-thinning fluids, where a streaky flow with an azimuthal wave number n=1 is observed in the transitional regime. Here, a linear stability analysis of pipe flow of shear-thinning fluids modulated azimuthally by finite amplitude streaks is performed. The shear-thinning behavior of the fluid is described by the Carreau model. The streaky base flows considered are obtained from two-dimensional direct numerical simulation using finite amplitude longitudinal rolls as the initial condition and by extracting the velocity field at time t(max), where the amplitude of the streaks reaches its maximum, denoted by A(max). It is found that the amplitude A(max) increases with increasing Reynolds number as well as with increasing amplitude E(0) of the initial longitudinal rolls. For sufficiently large streaks amplitude, streamwise velocity profiles develop inflection points, leading to instabilities. Depending on the threshold amplitude A(c), two different modes may trigger the instability of the streaks. If A(c) exceeds approximately 41.5% of the centerline velocity, the instability mode is located near the axis of the pipe, i.e., it is a "center mode." For weaker amplitude A(c), the instability mode is located near the pipe wall, in the region of highest wall normal shear, i.e., it is a "wall mode." The threshold amplitude A(c) decreases with increasing shear-thinning effects. The energy equation analysis indicates that (i) wall modes are driven mainly by the work of the Reynolds stress against the wall normal shear and (ii) for center modes, the contribution of the normal wall shear remains dominant; however, it is noted that the contribution of the Reynolds stress against the azimuthal shear increases with increasing shear-thinning effects.

Fold-pitchfork bifurcation for maps with Z(2) symmetry in pipe flow

This study aims to provide a better understanding of recently identified transition scenarios exhibited by traveling wave solutions in pipe flow. This particular family of solutions are invariant under certain reflectional symmetry transformations and they emerge from saddle-node bifurcations within a two-dimensional parameter space characterized by the length of the pipe and the Reynolds number. The present work precisely provides a detailed analysis of a codimension-two saddle-node bifurcation arising in discrete dynamical systems (maps) with Z(2) symmetry. Normal form standard techniques are applied in order to obtain the reduced map up to cubic order. All possible bifurcation scenarios exhibited by this normal form are analyzed in detail. Finally, a qualitative comparison of these scenarios with the ones observed in the aforementioned hydrodynamic problem is provided.

Streamwise-localized solutions at the onset of turbulence in pipe flow

Although the equations governing fluid flow are well known, there are no analytical expressions that describe the complexity of turbulent motion. A recent proposition is that in analogy to low dimensional chaotic systems, turbulence is organized around unstable solutions of the governing equations which provide the building blocks of the disordered dynamics. We report the discovery of periodic solutions which just like intermittent turbulence are spatially localized and show that turbulent transients arise from one such solution branch.

Cartography of high-dimensional flows: a visual guide to sections and slices

Symmetry reduction by the method of slices quotients the continuous symmetries of chaotic flows by replacing the original state space by a set of charts, each covering a neighborhood of a dynamically important class of solutions, qualitatively captured by a "template." Together these charts provide an atlas of the symmetry-reduced "slice" of state space, charting the regions of the manifold explored by the trajectories of interest. Within the slice, relative equilibria reduce to equilibria and relative periodic orbits reduce to periodic orbits. Visualizations of these solutions and their unstable manifolds reveal their interrelations and the role they play in organizing turbulence/chaos.

Periodic orbits near onset of chaos in plane Couette flow

We track the secondary bifurcations of coherent states in plane Couette flow and show that they undergo a periodic doubling cascade that ends with a crisis bifurcation. We introduce a symbolic dynamics for the orbits and show that the ones that exist fall into the universal sequence described by Metropolis, Stein and Stein for unimodal maps. The periodic orbits cover much of the turbulent dynamics in that their temporal evolution overlaps with turbulent motions when projected onto a plane spanned by energy production and dissipation.

Transient chaotic rotating waves in a ring of unidirectionally coupled symmetric Bonhoeffer-van der Pol oscillators near a codimension-two bifurcation point

Propagating waves in a ring of unidirectionally coupled symmetric Bonhoeffer-van der Pol (BVP) oscillators were studied. The parameter values of the BVP oscillators were near a codimension-two bifurcation point around which oscillatory, monostable, and bistable states coexist. Bifurcations of periodic, quasiperiodic, and chaotic rotating waves were found in a ring of three oscillators. In rings of large numbers of oscillators with small coupling strength, transient chaotic waves were found and their duration increased exponentially with the number of oscillators. These exponential chaotic transients could be described by a coupled map model derived from the Poincaré map of a ring of three oscillators. The quasiperiodic rotating waves due to the mode interaction near the codimension-two bifurcation point were evidently responsible for the emergence of the transient chaotic rotating waves.

Edge state in pipe flow experiments

Recent numerical studies suggest that in pipe and related shear flows, the region of phase space separating laminar from turbulent motion is organized by a chaotic attractor, called an edge state, which mediates the transition process. We here confirm the existence of the edge state in laboratory experiments. We observe that it governs the dynamics during the decay of turbulence underlining its potential relevance for turbulence control. In addition we unveil two unstable traveling wave solutions underlying the experimental flow fields. This observation corroborates earlier suggestions that unstable solutions organize turbulence and its stability border.

Experimental scaling law for the subcritical transition to turbulence in plane Poiseuille flow

We present an experimental study of the transition to turbulence in a plane Poiseuille flow. Using a well-controlled perturbation, we analyze the flow by using extensive particle image velocimetry and flow visualization (using laser-induced fluorescence) measurements, and use the deformation of the mean velocity profile as a criterion to characterize the state of the flow. From a large parametric study, four different states are defined, depending on the values of the Reynolds number and the amplitude of the perturbation. We discuss the role of coherent structures, such as hairpin vortices, in the transition. We find that the minimal amplitude of the perturbation triggering transition scales asymptotically as Re(-1).

Periodic magnetorotational dynamo action as a prototype of nonlinear magnetic-field generation in shear flows

The nature of dynamo action in shear flows prone to magnetohydrodynamc instabilities is investigated using the magnetorotational dynamo in Keplerian shear flow as a prototype problem. Using direct numerical simulations and Newton's method, we compute an exact time-periodic magnetorotational dynamo solution to three-dimensional dissipative incompressible magnetohydrodynamic equations with rotation and shear. We discuss the physical mechanism behind the cycle and show that it results from a combination of linear and nonlinear interactions between a large-scale axisymmetric toroidal magnetic field and nonaxisymmetric perturbations amplified by the magnetorotational instability. We demonstrate that this large-scale dynamo mechanism is overall intrinsically nonlinear and not reducible to the standard mean-field dynamo formalism. Our results therefore provide clear evidence for a generic nonlinear generation mechanism of time-dependent coherent large-scale magnetic fields in shear flows and call for new theoretical dynamo models. These findings may offer important clues to understanding the transitional and statistical properties of subcritical magnetorotational turbulence.

Traveling hairpin-shaped fluid vortices in plane Couette flow

Traveling-wave solutions are discovered in plane Couette flow. They are obtained when the so-called steady hairpin vortex state found recently by Gibson [J. Fluid Mech. 638, 243 (2009)] and Itano and Generalis [Phys. Rev. Lett. 102, 114501 (2009)] is continued to sliding Couette flow geometry between two concentric cylinders by using the radius ratio as a homotopy parameter. It turns out that in the plane Couette flow geometry two traveling waves having the phase velocities with opposite signs are associated with their appearance from the steady hairpin vortex state, where the amplitude of the phase velocities increases gradually from zero as the Reynolds number is increased. The solutions obviously inherit the streaky structure of the hairpin vortex state, but shape preserving flow patterns propagate in the streamwise direction. Other striking features of the solution are asymmetric mean flow profiles and strong quasistreamwise vortices which occupy the vicinity of only the top or bottom moving boundary, depending on the sign of the phase velocity. Furthermore, we find that the pitchfork bifurcation associated with the appearance of the solution becomes imperfect when the flow is perturbed by a Poiseuille flow component.

Rapid path to transition via nonlinear localized optimal perturbations in a boundary-layer flow

Recent studies have suggested that in some cases transition can be triggered by some purely nonlinear mechanisms. Here we aim at verifying such an hypothesis, looking for a localized perturbation able to lead a boundary-layer flow to a chaotic state, following a nonlinear route. Nonlinear optimal localized perturbations have been computed by means of an energy optimization which includes the nonlinear terms of the Navier-Stokes equations. Such perturbations lie on the turbulent side of the laminar-turbulent boundary, whereas, for the same value of the initial energy, their linear counterparts do not. The evolution of these perturbations toward a turbulent flow involves the presence of streamwise-inclined vortices at short times and of hairpin structures prior to breakdown.

Optimally amplified large-scale streaks and drag reduction in turbulent pipe flow

The optimal amplifications of small coherent perturbations within turbulent pipe flow are computed for Reynolds numbers up to one million. Three standard frameworks are considered: the optimal growth of an initial condition, the response to harmonic forcing and the Karhunen-Loève (proper orthogonal decomposition) analysis of the response to stochastic forcing. Similar to analyses of the turbulent plane channel flow and boundary layer, it is found that streaks elongated in the streamwise direction can be greatly amplified from quasistreamwise vortices, despite linear stability of the mean flow profile. The most responsive perturbations are streamwise uniform and, for sufficiently large Reynolds number, the most responsive azimuthal mode is of wave number m=1 . The response of this mode increases with the Reynolds number. A secondary peak, where m corresponds to azimuthal wavelengths λ_{θ}^{+}≈70-90 in wall units, also exists in the amplification of initial conditions and in premultiplied response curves for the forced problems. Direct numerical simulations at Re=5300 confirm that the forcing of m=1,2 and m=4 optimal structures results in the large response of coherent large-scale streaks. For moderate amplitudes of the forcing, low-speed streaks become narrower and more energetic, whereas high-speed streaks become more spread. It is further shown that drag reduction can be achieved by forcing steady large-scale structures, as anticipated from earlier investigations. Here the energy balance is calculated. At Re=5300 it is shown that, due to the small power required by the forcing of optimal structures, a net power saving of the order of 10% can be achieved following this approach, which could be relevant for practical applications.

Towards minimal perturbations in transitional plane Couette flow

For parallel shear flows, transition to turbulence occurs only for perturbations of sufficiently large amplitude. It is therefore relevant to study the shape, amplitude, and dynamics of the least energetic initial disturbances leading to transition. We suggest a numerical approach to find such minimal perturbations, applied here to the case of plane Couette flow. The optimization method seeks such perturbations at initial time as a linear combination of a finite number of linear optimal modes. The energy threshold of the minimal perturbation for a Reynolds number Re=400 is only 2% less than for a pair of symmetric oblique waves. The associated transition scenario shows a long transient approach to a steady state solution with special symmetries. Modal analysis shows how the oblique-wave mechanism can be optimized by the addition of other oblique modes breaking the flow symmetry and whose nonlinear interaction generates spectral components of the edge state. The Re dependence of energy thresholds is revisited, with evidence for a O(Re(-2)) -scaling for both oblique waves and streamwise vortices scenarios.

Linear and nonlinear stability of a thermally stratified magnetically driven rotating flow in a cylinder

The stability of a thermally stratified liquid metal flow is considered numerically. The flow is driven by a rotating magnetic field in a cylinder heated from above and cooled from below. The stable thermal stratification turns out to destabilize the flow. This is explained by the fact that a stable stratification suppresses the secondary meridional flow, thus indirectly enhancing the primary rotation. The instability in the form of Taylor-Görtler rolls is consequently promoted. These rolls can only be excited by finite disturbances in the isothermal flow. A sufficiently strong thermal stratification transforms this nonlinear bypass instability into a linear one reducing, thus, the critical value of the magnetic driving force. A weaker temperature gradient delays the linear instability but makes the bypass transition more likely. We quantify the non-normal and nonlinear components of this transition by direct numerical simulation of the flow response to noise. It is observed that the flow sensitivity to finite disturbances increases considerably under the action of a stable thermal stratification. The capabilities of the random forcing approach to identify disconnected coherent states in a general case are discussed.

Dynamics and Control of Global Instabilities in Open-Flows: A Linearized Approach

This review article addresses the dynamics and control of low-frequency unsteadiness, as observed in some aerodynamic applications. It presents a coherent and rigorous linearized approach, which enables both to describe the dynamics of commonly encountered open-flows and to design open-loop and closed-loop control strategies, in view of suppressing or delaying instabilities. The approach is global in the sense that both cross-stream and streamwise directions are discretized in the evolution operator. New light will therefore be shed on the streamwise properties of open-flows. In the case of oscillator flows, the unsteadiness is due to the existence of unstable global modes, i.e., unstable eigenfunctions of the linearized Navier–Stokes operator. The influence of nonlinearities on the dynamics is studied by deriving nonlinear amplitude equations, which accurately describe the dynamics of the flow in the vicinity of the bifurcation threshold. These equations also enable us to analyze the mean flow induced by the nonlinearities as well as the stability properties of this flow. The open-loop control of unsteadiness is then studied by a sensitivity analysis of the eigenvalues with respect to base-flow modifications. With this approach, we manage to a priori identify regions of the flow where a small control cylinder suppresses unsteadiness. Then, a closed-loop control approach was implemented for the case of an unstable open-cavity flow. We have combined model reduction techniques and optimal control theory to stabilize the unstable eigenvalues. Various reduced-order-models based on global modes, proper orthogonal decomposition modes, and balanced modes were tested and evaluated according to their ability to reproduce the input-output behavior between the actuator and the sensor. Finally, we consider the case of noise-amplifiers, such as boundary-layer flows and jets, which are stable when viewed in a global framework. The importance of the singular value decomposition of the global resolvent will be highlighted in order to understand the frequency selection process in such flows.

Distinct large-scale turbulent-laminar states in transitional pipe flow

When fluid flows through a channel, pipe, or duct, there are two basic forms of motion: smooth laminar motion and complex turbulent motion. The discontinuous transition between these states is a fundamental problem that has been studied for more than 100 yr. What has received far less attention is the large-scale nature of the turbulent flows near transition once they are established. We have carried out extensive numerical computations in pipes of variable lengths up to 125 diameters to investigate the nature of transitional turbulence in pipe flow. We show the existence of three fundamentally different turbulent states separated by two distinct Reynolds numbers. Below Re (1) approximately equal 2,300, turbulence takes the form of familiar equilibrium (or longtime transient) puffs that are spatially localized and keep their size independent of pipe length. At Re (1) the flow makes a striking transition to a spatio-temporally intermittent flow that fills the pipe. Irregular alternation of turbulent and laminar regions is inherent and does not result from random disturbances. The fraction of turbulence increases with Re until Re (2) approximately equal 2,600 where there is a continuous transition to a state of uniform turbulence along the pipe. We relate these observations to directed percolation and argue that Re (1) marks the onset of infinite-lifetime turbulence.