Global bifurcations to subcritical magnetorotational dynamo action in Keplerian shear flow

Magnetorotational dynamo action in Keplerian shear flow is a three-dimensional, non-linear magnetohydrodynamic process whose study is relevant to the understanding of accretion processes and magnetic field generation in astrophysics. Transition to this form of dynamo action is subcritical and shares many characteristics of transition to turbulence in non-rotating hydrodynamic shear flows. This suggests that these different fluid systems become active through similar generic bifurcation mechanisms, which in both cases have eluded detailed understanding so far. In this paper, we build on recent work on the two problems to investigate numerically the bifurcation mechanisms at work in the incompressible Keplerian magnetorotational dynamo problem in the shearing box framework. Using numerical techniques imported from dynamical systems research, we show that the onset of chaotic dynamo action at magnetic Prandtl numbers larger than unity is primarily associated with global homoclinic and heteroclinic bifurcations of nonlinear magnetorotational dynamo cycles. These global bifurcations are found to be supplemented by local bifurcations of cycles marking the beginning of period-doubling cascades. The results suggest that nonlinear magnetorotational dynamo cycles provide the pathway to turbulent injection of both kinetic and magnetic energy in incompressible magnetohydrodynamic Keplerian shear flow in the absence of an externally imposed magnetic field. Studying the nonlinear physics and bifurcations of these cycles in different regimes and configurations may subsequently help to better understand the physical conditions of excitation of magnetohydrodynamic turbulence and instability-driven dynamos in a variety of astrophysical systems and laboratory experiments. The detailed characterization of global bifurcations provided for this three-dimensional subcritical fluid dynamics problem may also prove useful for the problem of transition to turbulence in hydrodynamic shear flows

Visualizing the geometry of state space in plane Couette flow

Motivated by recent experimental and numerical studies of coherent structures in wall-bounded shear flows, we initiate a systematic exploration of the hierarchy of unstable invariant solutions of the Navier-Stokes equations. We construct a dynamical, 10^5-dimensional state-space representation of plane Couette flow at Re = 400 in a small, periodic cell and offer a new method of visualizing invariant manifolds embedded in such high dimensions. We compute a new equilibrium solution of plane Couette flow and the leading eigenvalues and eigenfunctions of known equilibria at this Reynolds number and cell size. What emerges from global continuations of their unstable manifolds is a surprisingly elegant dynamical-systems visualization of moderate-Reynolds turbulence. The invariant manifolds tessellate the region of state space explored by transiently turbulent dynamics with a rigid web of continuous and discrete symmetry-induced heteroclinic connections.Comment: 32 pages, 13 figures submitted to Journal of Fluid Mechanic

On the state space geometry of the Kuramoto-Sivashinsky flow in a periodic domain

The continuous and discrete symmetries of the Kuramoto-Sivashinsky system restricted to a spatially periodic domain play a prominent role in shaping the invariant sets of its chaotic dynamics. The continuous spatial translation symmetry leads to relative equilibrium (traveling wave) and relative periodic orbit (modulated traveling wave) solutions. The discrete symmetries lead to existence of equilibrium and periodic orbit solutions, induce decomposition of state space into invariant subspaces, and enforce certain structurally stable heteroclinic connections between equilibria. We show, on the example of a particular small-cell Kuramoto-Sivashinsky system, how the geometry of its dynamical state space is organized by a rigid `cage' built by heteroclinic connections between equilibria, and demonstrate the preponderance of unstable relative periodic orbits and their likely role as the skeleton underpinning spatiotemporal turbulence in systems with continuous symmetries. We also offer novel visualizations of the high-dimensional Kuramoto-Sivashinsky state space flow through projections onto low-dimensional, PDE representation independent, dynamically invariant intrinsic coordinate frames, as well as in terms of the physical, symmetry invariant energy transfer rates.Comment: 31 pages, 17 figures; added references, corrected typos. Due to file size restrictions some figures in this preprint are of low quality. A high quality copy may be obtained from http://www.cns.gatech.edu/~predrag/papers/preprints.html#rp

Low-dimensional models for turbulent plane Couette flow in a minimal flow unit

We model turbulent plane Couette flow in the minimal flow unit (MFU) – a domain whose spanwise and streamwise extent is just sufficient to maintain turbulence – by expanding the velocity field as a sum of optimal modes calculated via proper orthogonal decomposition from numerical data. Ordinary differential equations are obtained by Galerkin projection of the Navier–Stokes equations onto these modes. We first consider a 6-mode (11-dimensional) model and study the effects of including losses to neglected modes. Ignoring these, the model reproduces turbulent statistics acceptably, but fails to reproduce dynamics; including them, we find a stable periodic orbit that captures the regeneration cycle dynamics and agrees well with direct numerical simulations. However, restriction to as few as six modes artificially constrains the relative magnitudes of streamwise vortices and streaks and so cannot reproduce stability of the laminar state or properly account for bifurcations to turbulence as Reynolds number increases. To address this issue, we develop a second class of models based on ‘uncoupled’ eigenfunctions that allow independence among streamwise and cross-stream velocity components. A 9-mode (31-dimensional) model produces bifurcation diagrams for steady and periodic states in qualitative agreement with numerical Navier–Stokes solutions, while preserving the regeneration cycle dynamics. Together, the models provide empirical evidence that the ‘backbone’ for MFU turbulence is a periodic orbit, and support the roll–streak–breakdown–roll reformation picture of shear-driven turbulence

Diffusion of a passive scalar from a no-slip boundary into a two-dimensional chaotic advection field

Using a time-periodic perturbation of a two-dimensional steady separation bubble on a plane no-slip boundary to generate chaotic particle trajectories in a localized region of an unbounded boundary layer flow, we study the impact of various geometrical structures that arise naturally in chaotic advection fields on the transport of a passive scalar from a local 'hot spot' on the no-slip boundary. We consider here the full advection-diffusion problem, though attention is restricted to the case of small scalar diffusion, or large Peclet number. In this regime, a certain one-dimensional unstable manifold is shown to be the dominant organizing structure in the distribution of the passive scalar. In general, it is found that the chaotic structures in the flow strongly influence the scalar distribution while, in contrast, the flux of passive scalar from the localized active no-slip surface is, to dominant order, independent of the overlying chaotic advection. Increasing the intensity of the chaotic advection by perturbing the velocity held further away from integrability results in more non-uniform scalar distributions, unlike the case in bounded flows where the chaotic advection leads to rapid homogenization of diffusive tracer. In the region of chaotic particle motion the scalar distribution attains an asymptotic state which is time-periodic, with the period the same as that of the time-dependent advection field. Some of these results are understood by using the shadowing property from dynamical systems theory. The shadowing property allows us to relate the advection-diffusion solution at large Peclet numbers to a fictitious zero-diffusivity or frozen-field solution, corresponding to infinitely large Peclet number. The zero-diffusivity solution is an unphysical quantity, but is found to be a powerful heuristic tool in understanding the role of small scalar diffusion. A novel feature in this problem is that the chaotic advection field is adjacent to a no-slip boundary. The interaction between the necessarily non-hyperbolic particle dynamics in a thin near-wall region and the strongly hyperbolic dynamics in the overlying chaotic advection field is found to have important consequences on the scalar distribution; that this is indeed the case is shown using shadowing. Comparisons are made throughout with the flux and the distributions of the passive scalar for the advection-diffusion problem corresponding to the steady, unperturbed, integrable advection field