Extreme events in transitional turbulence.

Transitional localized turbulence in shear flows is known to either decay to an absorbing laminar state or to proliferate via splitting. The average passage times from one state to the other depend super-exponentially on the Reynolds number and lead to a crossing Reynolds number above which proliferation is more likely than decay. In this paper, we apply a rare-event algorithm, Adaptative Multilevel Splitting, to the deterministic Navier-Stokes equations to study transition paths and estimate large passage times in channel flow more efficiently than direct simulations. We establish a connection with extreme value distributions and show that transition between states is mediated by a regime that is self-similar with the Reynolds number. The super-exponential variation of the passage times is linked to the Reynolds number dependence of the parameters of the extreme value distribution. Finally, motivated by instantons from Large Deviation theory, we show that decay or splitting events approach a most-probable pathway. This article is part of the theme issue 'Mathematical problems in physical fluid dynamics (part 2)'.

A fluid mechanic's analysis of the teacup singularity.

The mechanism for singularity formation in an inviscid wall-bounded fluid flow is investigated. The incompressible Euler equations are numerically simulated in a cylindrical container. The flow is axisymmetric with the swirl. The simulations reproduce and corroborate aspects of prior studies reporting strong evidence for a finite-time singularity. The analysis here focuses on the interplay between inertia and pressure, rather than on vorticity. The linearity of the pressure Poisson equation is exploited to decompose the pressure field into independent contributions arising from the meridional flow and from the swirl, and enforcing incompressibility and enforcing flow confinement. The key pressure field driving the blowup of velocity gradients is that confining the fluid within the cylinder walls. A model is presented based on a primitive-variables formulation of the Euler equations on the cylinder wall, with closure coming from how pressure is determined from velocity. The model captures key features in the mechanics of the blowup scenario.

The rise of fully turbulent flow.

Over a century of research into the origin of turbulence in wall-bounded shear flows has resulted in a puzzling picture in which turbulence appears in a variety of different states competing with laminar background flow. At moderate flow speeds, turbulence is confined to localized patches; it is only at higher speeds that the entire flow becomes turbulent. The origin of the different states encountered during this transition, the front dynamics of the turbulent regions and the transformation to full turbulence have yet to be explained. By combining experiments, theory and computer simulations, here we uncover a bifurcation scenario that explains the transformation to fully turbulent pipe flow and describe the front dynamics of the different states encountered in the process. Key to resolving this problem is the interpretation of the flow as a bistable system with nonlinear propagation (advection) of turbulent fronts. These findings bridge the gap between our understanding of the onset of turbulence and fully turbulent flows.

Prediction of frequencies in thermosolutal convection from mean flows.

Motivated by studies of the cylinder wake, in which the vortex-shedding frequency can be obtained from the mean flow, we study thermosolutal convection driven by opposing thermal and solutal gradients. In the archetypal two-dimensional geometry with horizontally periodic and vertical no-slip boundary conditions, branches of traveling waves and standing waves are created simultaneously by a Hopf bifurcation. Consistent with similar analyses performed on the cylinder wake, we find that the traveling waves of thermosolutal convection have the RZIF property, meaning that linearization about the mean fields of the traveling waves yields an eigenvalue whose real part is almost zero and whose imaginary part corresponds very closely to the nonlinear frequency. In marked contrast, linearization about the mean field of the standing waves yields neither zero growth nor the nonlinear frequency. It is shown that this difference can be attributed to the fact that the temporal power spectrum for the traveling waves is peaked, while that of the standing waves is broad. We give a general demonstration that the frequency of any quasimonochromatic oscillation can be predicted from its temporal mean.

Asymptotic dynamics of reflecting spiral waves.

Resonantly forced spiral waves in excitable media drift in straight-line paths, their rotation centers behaving as pointlike objects moving along trajectories with a constant velocity. Interaction with medium boundaries alters this velocity and may often result in a reflection of the drift trajectory. Such reflections have diverse characteristics and are known to be highly nonspecular in general. In this context we apply the theory of response functions, which via numerically computable integrals, reduces the reaction-diffusion equations governing the whole excitable medium to the dynamics of just the rotation center and rotation phase of a spiral wave. Spiral reflection trajectories are computed by this method for both small- and large-core spiral waves in the Barkley model. Such calculations provide insight into the process of reflection as well as explanations for differences in trajectories across parameters, including the effects of incidence angle and forcing amplitude. Qualitative aspects of these results are preserved far beyond the asymptotic limit of weak boundary effects and slow resonant drift.

Non-specular reflections in a macroscopic system with wave-particle duality: spiral waves in bounded media.

Spiral waves in excitable media possess both wave-like and particle-like properties. When resonantly forced (forced at the spiral rotation frequency) spiral cores travel along straight trajectories, but may reflect from medium boundaries. Here, numerical simulations are used to study reflections from two types of boundaries. The first is a no-flux boundary which waves cannot cross, while the second is a step change in the medium excitability which waves do cross. Both small-core and large-core spirals are investigated. The predominant feature in all cases is that the reflected angle varies very little with incident angle for large ranges of incident angles. Comparisons are made to the theory of Biktashev and Holden. Large-core spirals exhibit other phenomena such as binding to boundaries. The dynamics of multiple reflections is briefly considered.

Simplifying the complexity of pipe flow.

Transitional pipe flow is modeled as a one-dimensional excitable and bistable medium. Models are presented in two variables, turbulence intensity and mean shear, that evolve according to established properties of transitional turbulence. A continuous model captures the essence of the puff-slug transition as a change from excitability to bistability. A discrete model, which additionally incorporates turbulence locally as a chaotic repeller, reproduces almost all large-scale features of transitional pipe flow. In particular, it captures metastable localized puffs, puff splitting, slugs, localized edge states, a continuous transition to sustained turbulence via spatiotemporal intermittency (directed percolation), and a subsequent increase in turbulence fraction toward uniform, featureless turbulence.

The onset of turbulence in pipe flow.

Shear flows undergo a sudden transition from laminar to turbulent motion as the velocity increases, and the onset of turbulence radically changes transport efficiency and mixing properties. Even for the well-studied case of pipe flow, it has not been possible to determine at what Reynolds number the motion will be either persistently turbulent or ultimately laminar. We show that in pipes, turbulence that is transient at low Reynolds numbers becomes sustained at a distinct critical point. Through extensive experiments and computer simulations, we were able to identify and characterize the processes ultimately responsible for sustaining turbulence. In contrast to the classical Landau-Ruelle-Takens view that turbulence arises from an increase in the temporal complexity of fluid motion, here, spatial proliferation of chaotic domains is the decisive process and intrinsic to the nature of fluid turbulence.

Influence of counter-rotating von Kármán flow on cylindrical Rayleigh-Bénard convection.

The axisymmetric flow in an aspect-ratio-one cylinder whose upper and lower bounding disks are maintained at different temperatures and rotate at equal and opposite velocities is investigated. In this combined Rayleigh-Bénard/von Kármán problem, the imposed temperature gradient is measured by the Rayleigh number Ra and the angular velocity by the Reynolds number Re. Although fluid motion is present as soon as Re not equal 0 , a symmetry-breaking transition analogous to the onset of convection takes place at a finite Rayleigh number higher than that for Re=0 . For Re<95 , the transition is a pitchfork bifurcation to a pair of steady states, while for Re>95 , it is a Hopf bifurcation to a limit cycle. The steady states and limit cycle are connected via a pair of saddle-node infinite-period bifurcations except very near the Takens-Bogdanov codimension-two point, where the scenario includes global bifurcations. Detailed phase portraits and bifurcation diagrams are presented, as well as the evolution of the leading part of the spectrum, over the parameter ranges 0

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.

Confined three-dimensional stability analysis of the cylinder wake.

A spatially confined stability analysis is reported for the cylinder wake at Reynolds numbers 190 and 260. The two three-dimensional instabilities at these Reynolds numbers are shown to be driven by the flow just behind the cylinder.

Computational study of turbulent laminar patterns in couette flow.

Turbulent-laminar patterns near transition are simulated in plane Couette flow using an extension of the minimal-flow-unit methodology. Computational domains are of minimal size in two directions but large in the third. The long direction can be tilted at any prescribed angle to the streamwise direction. Three types of patterned states are found and studied: periodic, localized, and intermittent. These correspond closely to observations in large-aspect-ratio experiments.

Large-excitability asymptotics for scroll waves in three-dimensional excitable media.

Three-dimensional scroll waves are considered in a reaction-diffusion model of excitable media in the large excitability limit. Coordinates based on the scroll filament are defined and shown to provide a natural extension of the coordinates used for two-dimensional spiral waves. The leading-order free-boundary equations for interface motion in three dimensions are explicitly derived in these coordinates. Three specific examples are considered: straight twisted scroll waves, axisymmetric scroll waves, and helical scroll waves. The equations for the fields at leading and first order in the core region are given.

Cookbook asymptotics for spiral and scroll waves in excitable media.

Algebraic formulas predicting the frequencies and shapes of waves in a reaction-diffusion model of excitable media are presented in the form of four recipes. The formulas themselves are based on a detailed asymptotic analysis (published elsewhere) of the model equations at leading order and first order in the asymptotic parameter. The importance of the first order contribution is stressed throughout, beginning with a discussion of the Fife limit, Fife scaling, and Fife regime. Recipes are given for spiral waves and detailed comparisons are presented between the asymptotic predictions and the solutions of the full reaction-diffusion equations. Recipes for twisted scroll waves with straight filaments are given and again comparisons are shown. The connection between the asymptotic results and filament dynamics is discussed, and one of the previously unknown coefficients in the theory of filament dynamics is evaluated in terms of its asymptotic expansion. (c) 2002 American Institute of Physics.

A dynamical systems approach to spiral wave dynamics.

A simple system of five nonlinear ordinary differential equations is shown to reproduce many dynamical features of spiral waves in two-dimensional excitable media.