Chaotic dynamos generated by a turbulent flow of liquid sodium

Multiple States in Turbulent Large-Aspect-Ratio Thermal Convection: What Determines the Number of Convection Rolls?

Wall-bounded turbulent flows can take different statistically stationary turbulent states, with different transport properties, even for the very same values of the control parameters. What state the system takes depends on the initial conditions. Here we analyze the multiple states in large-aspect ratio (Γ) two-dimensional turbulent Rayleigh-Bénard flow with no-slip plates and horizontally periodic boundary conditions as model system. We determine the number n of convection rolls, their mean aspect ratios Γ_{r}=Γ/n, and the corresponding transport properties of the flow (i.e., the Nusselt number Nu), as function of the control parameters Rayleigh (Ra) and Prandtl number. The effective scaling exponent β in Nu∼Ra^{β} is found to depend on the realized state and thus Γ_{r}, with a larger value for the smaller Γ_{r}. By making use of a generalized Friedrichs inequality, we show that the elliptical shape of the rolls and viscous damping determine the Γ_{r} window for the realizable turbulent states. The theoretical results are in excellent agreement with our numerical finding 2/3≤Γ_{r}≤4/3, where the lower threshold is approached for the larger Ra. Finally, we show that the theoretical approach to frame Γ_{r} also works for free-slip boundary conditions.

Dynamo Enhancement and Mode Selection Triggered by High Magnetic Permeability

We present results from consistent dynamo simulations, where the electrically conducting and incompressible flow inside a cylinder vessel is forced by moving impellers numerically implemented by a penalization method. The numerical scheme models jumps of magnetic permeability for the solid impellers, resembling various configurations tested experimentally in the von Kármán sodium experiment. The most striking experimental observations are reproduced in our set of simulations. In particular, we report on the existence of a time-averaged axisymmetric dynamo mode, self-consistently generated when the magnetic permeability of the impellers exceeds a threshold. We describe a possible scenario involving both the turbulent flow in the vicinity of the impellers and the high magnetic permeability of the impellers.

Stochastic reversal dynamics of two interacting magnetic dipoles: A simple model experiment

We report an experimental study of the dynamics of two coupled magnetic dipoles. The experiment consists in two coplanar permanent disk magnets separated by a distance d, each allowed to rotate on a fixed parallel axis-each magnet's axis being perpendicular to its dipolar moment vector. A torque of adjustable strength can be externally applied to one of the magnets, the other magnet being free. The driving torque may be time-independent or temporally fluctuating. We study the influence of the parameters of the driving torque on the dynamics of the coupled system, in particular the emergence of dynamical regimes such as stochastic reversals. We report transitions between stationary and stochastic reversal regimes. All the observed features can be understood by a simple mechanical dynamical model. The transition between statistically stationary regimes and reversals is explained introducing an effective potential energy incorporating both the coupling between magnets and the external driving. Relations between this simple experimental model with macroscopic models of magnetic spin coupling, as well as with chaotic reversals of turbulent dynamos, are discussed.

Recurrent bursts via linear processes in turbulent environments

Large-scale instabilities occurring in the presence of small-scale turbulent fluctuations are frequently observed in geophysical or astrophysical contexts but are difficult to reproduce in the laboratory. Using extensive numerical simulations, we report here on intense recurrent bursts of turbulence in plane Poiseuille flow rotating about a spanwise axis. A simple model based on the linear instability of the mean flow can predict the structure and time scale of the nearly periodic and self-sustained burst cycles. Poiseuille flow is suggested as a prototype for future studies of low-dimensional dynamics embedded in strongly turbulent environments.

Decay rates of magnetic modes below the threshold of a turbulent dynamo

We measure the decay rates of magnetic field modes in a turbulent flow of liquid sodium below the dynamo threshold. We observe that turbulent fluctuations induce energy transfers between modes with different symmetries (dipolar and quadrupolar). Using symmetry properties, we show how to measure the decay rate of each mode without being restricted to the one with the smallest damping rate. We observe that the respective values of the decay rates of these modes depend on the shape of the propellers driving the flow. Dynamical regimes, including field reversals, are observed only when the modes are both nearly marginal. This is in line with a recently proposed model.

Effect of soft-iron impellers on the von Kármán-sodium dynamo

The explanation for the observed axisymmetric magnetic field in the von Kármán-sodium (VKS) dynamo experiment is still an unresolved question. In this paper, the integral equation approach is extended to investigate the VKS dynamo action by taking into account the discontinuity of the magnetic permeability and electrical conductivity in the conducting region. When the relative magnetic permeability of the soft-iron impellers is set to 65, a steady toroidal field that is apparently axisymmetric is excited at the critical magnetic Reynolds number, Rmc≈27.23, which is close to the experimental result, Rmc≈30. Our results show that the critical magnetic Reynolds number declines as the relative magnetic permeability of the impellers increases. Furthermore, when the relative magnetic permeability is not greater than 37, an equatorial magnetic field with an azimuthal wave number m=1 is the dominant mode, otherwise a steady toroidal field with an azimuthal wave number m=0 predominates the magnetic field generated by the VKS dynamo action.

On the genesis of the Earth's magnetism

Few areas of geophysics are today progressing as rapidly as basic geomagnetism, which seeks to understand the origin of the Earth's magnetism. Data about the present geomagnetic field pours in from orbiting satellites, and supplements the ever growing body of information about the field in the remote past, derived from the magnetism of rocks. The first of the three parts of this review summarizes the available geomagnetic data and makes significant inferences about the large scale structure of the geomagnetic field at the surface of the Earth's electrically conducting fluid core, within which the field originates. In it, we recognize the first major obstacle to progress: because of the Earth's mantle, only the broad, slowly varying features of the magnetic field within the core can be directly observed. The second (and main) part of the review commences with the geodynamo hypothesis: the geomagnetic field is induced by core flow as a self-excited dynamo. Its electrodynamics define 'kinematic dynamo theory'. Key processes involving the motion of magnetic field lines, their diffusion through the conducting fluid, and their reconnection are described in detail. Four kinematic models are presented that are basic to a later section on successful dynamo experiments. The fluid dynamics of the core is considered next, the fluid being driven into motion by buoyancy created by the cooling of the Earth from its primordial state. The resulting flow is strongly affected by the rotation of the Earth and by the Lorentz force, which alters fluid motion by the interaction of the electric current and magnetic field. A section on 'magnetohydrodynamic (MHD) dynamo theory' is devoted to this rotating magnetoconvection. Theoretical treatment of the MHD responsible for geomagnetism culminates with numerical solutions of its governing equations. These simulations help overcome the first major obstacle to progress, but quickly meet the second: the dynamics of Earth's core are too complex, and operate across time and length scales too broad to be captured by any single laboratory experiment, or resolved on present-day computers. The geophysical relevance of the experiments and simulations is therefore called into question. Speculation about what may happen when computational power is eventually able to resolve core dynamics is given considerable attention. The final part of the review is a postscript to the earlier sections. It reflects on the problems that geodynamo theory will have to solve in the future, particularly those that core turbulence presents.

Dynamo threshold detection in the von Kármán sodium experiment

Predicting dynamo self-generation in liquid metal experiments has been an ongoing question for many years. In contrast to simple dynamical systems for which reliable techniques have been developed, the ability to predict the dynamo capacity of a flow and the estimate of the corresponding critical value of the magnetic Reynolds number (the control parameter of the instability) has been elusive, partly due to the high level of turbulent fluctuations of flows in such experiments (with kinetic Reynolds numbers in excess of 10(6)). We address these issues here, using the von Kármán sodium experiment and studying its response to an externally applied magnetic field. We first show that a dynamo threshold can be estimated from analysis related to critical slowing down and susceptibility divergence, in configurations for which dynamo action is indeed observed. These approaches are then applied to flow configurations that have failed to self-generate magnetic fields within operational limits, and we quantify the dynamo capacity of these configurations.

On the dynamic nature of azimuthal thermoacoustic modes in annular gas turbine combustion chambers

This paper deals with the dynamics of standing and rotating azimuthal thermoacoustic modes in annular combustion chambers. Simultaneous acoustic measurements have been made at multiple circumferential positions in an annular gas turbine combustion chamber. A detailed statistical analysis of the spatial Fourier amplitudes extracted from these data reveals that the acoustic modes are continuously switching between standing, clockwise and counter-clockwise travelling waves. A theoretical framework from which the modal dynamics can be explained is proposed and supported by real gas turbine data. The stochastic differential equations that govern these systems have been derived and used as a basis for system identification of the measured engine data. The model describes the probabilities of the two azimuthal wave components as a function of the random source intensity, the asymmetry in the system and the strength of the thermoacoustic interaction. The solution of the simplified system is in good agreement with experimental observations on a gas turbine combustion chamber.

Symmetry-breaking Hopf bifurcations to 1-, 2-, and 3-tori in small-aspect-ratio counterrotating Taylor-Couette flow

The nonlinear dynamics of Taylor-Couette flow in a small-aspect-ratio wide-gap annulus in the counterrotating regime is investigated by solving the full three-dimensional Navier-Stokes equations. The system is invariant under arbitrary rotations about the axis, reflection about the annulus midplane, and time translations. A systematic investigation is presented both in terms of the flow physics elucidated from the numerical simulations and from a dynamical system perspective provided by equivariant normal form theory. The dynamics are primarily associated with the behavior of the jet of angular momentum that emerges from the inner cylinder boundary layer at about the midplane. The sequence of bifurcations as the differential rotation is increased consists of an axisymmetric Hopf bifurcation breaking the reflection symmetry of the basic state leading to an axisymmetric limit cycle with a half-period-flip spatiotemporal symmetry. This undergoes a Hopf bifurcation breaking axisymmetry, leading to quasiperiodic solutions evolving on a 2-torus that is setwise symmetric. These undergo a further Hopf bifurcation, introducing a third incommensurate frequency leading to a 3-torus that is also setwise symmetric. On the 3-torus, as the differential rotation is further increased, a saddle-node-invariant-circle bifurcation takes place, destroying the 3-torus and leaving a pair of symmetrically related 2-tori states on which all symmetries of the system have been broken.

Role of large-scale velocity fluctuations in a two-vortex kinematic dynamo

This paper presents an analysis of the Dudley-James two-vortex flow, which inspired several laboratory-scale liquid-metal experiments, in order to better demonstrate its relation to astrophysical dynamos. A coordinate transformation splits the flow into components that are axisymmetric and nonaxisymmetric relative to the induced magnetic dipole moment. The reformulation gives the flow the same dynamo ingredients as are present in more complicated convection-driven dynamo simulations. These ingredients are currents driven by the mean flow and currents driven by correlations between fluctuations in the flow and fluctuations in the magnetic field. The simple model allows us to isolate the dynamics of the growing eigenvector and trace them back to individual three-wave couplings between the magnetic field and the flow. This simple model demonstrates the necessity of poloidal advection in sustaining the dynamo and points to the effect of large-scale flow fluctuations in exciting a dynamo magnetic field.

Experimental observation of spatially localized dynamo magnetic fields

We report the first experimental observation of a spatially localized dynamo magnetic field, a common feature of astrophysical dynamos and convective dynamo simulations. When the two propellers of the von Kármán sodium experiment are driven at frequencies that differ by 15%, the mean magnetic field's energy measured close to the slower disk is nearly 10 times larger than the one close to the faster one. This strong localization of the magnetic field when a symmetry of the forcing is broken is in good agreement with a prediction based on the interaction between a dipolar and a quadrupolar magnetic mode.

Connecting flow structures and heat flux in turbulent Rayleigh-Bénard convection

The aspect ratio (Γ) dependence of the heat transfer (Nusselt number Nu in dimensionless form) in turbulent (two-dimensional) Rayleigh-Bénard convection is numerically studied in the regime 0.4≤Γ≤1.25 for Rayleigh numbers 10(7)≤Ra≤Ra(9) and Prandtl numbers Pr=0.7 (gas) and 4.3 (water). Nu(Γ) shows a very rich structure with sudden jumps and sharp transitions. We connect these structures to the way the flow organizes itself in the sample and explain why the aspect ratio dependence of Nu is more pronounced for small Pr. Even for fixed Γ different turbulent states (with different resulting Nu) can exist, between which the flow can or cannot switch. In the latter case the heat transfer thus depends on the initial conditions.

Reducing global turbulent resistivity by eliminating large eddies in a spherical liquid-sodium experiment

Three-wave turbulent interactions and the role of eddy size on the turbulent electromotive force are studied in a spherical liquid-sodium dynamo experiment. A symmetric, equatorial baffle reduces the amplitude of the largest-scale turbulent eddies, which is inferred from the magnetic fluctuations spectrum (measured by a 2D array of surface probes). Differential rotation in the mean flow is >2 times more effective in generating mean toroidal magnetic fields from the applied poloidal field (via the Ω effect) when the largest-scale eddies are eliminated, thus demonstrating that the global turbulent resistivity (the β effect from the largest-scale eddies) is reduced by a similar amount.

Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh-Bénard convection

In turbulent thermal convection in cylindrical samples with an aspect ratio Γ≡D/L (D is the diameter and L the height), the Nusselt number Nu is enhanced when the sample is rotated about its vertical axis because of the formation of Ekman vortices that extract additional fluid out of thermal boundary layers at the top and bottom. We show from experiments and direct numerical simulations that the enhancement occurs only above a bifurcation point at a critical inverse Rossby number 1/Ro(c), with 1/Ro(c)∝1/Γ. We present a Ginzburg-Landau-like model that explains the existence of a bifurcation at finite 1/Ro(c) as a finite-size effect. The model yields the proportionality between 1/Ro(c) and 1/Γ and is consistent with several other measured or computed system properties.

Dipole-quadrupole dynamics during magnetic field reversals

The shape and the dynamics of reversals of the magnetic field in a turbulent dynamo experiment are investigated. We report the evolution of the dipolar and the quadrupolar parts of the magnetic field in the VKS experiment, and show that the experimental results are in good agreement with the predictions of a recent model of reversals: when the dipole reverses, part of the magnetic energy is transferred to the quadrupole, reversals begin with a slow decay of the dipole and are followed by a fast recovery, together with an overshoot of the dipole. Random reversals are observed at the borderline between stationary and oscillatory dynamos.

Statistical mechanics of Beltrami flows in axisymmetric geometry: theory reexamined

A simplified thermodynamic approach of the incompressible axisymmetric Euler equations is considered based on the conservation of helicity, angular momentum, and microscopic energy. Statistical equilibrium states are obtained by maximizing the Boltzmann entropy under these sole constraints. We assume that these constraints are selected by the properties of forcing and dissipation. The fluctuations are found to be Gaussian, while the mean flow is in a Beltrami state. Furthermore, we show that the maximization of entropy at fixed helicity, angular momentum, and microscopic energy is equivalent to the minimization of macroscopic energy at fixed helicity and angular momentum. This provides a justification of this selective decay principle from statistical mechanics. These theoretical predictions are in good agreement with experiments of a von Kármán turbulent flow and provide a way to measure the temperature of turbulence and check fluctuation-dissipation relations. Relaxation equations are derived that could provide an effective description of the dynamics toward the Beltrami state and the progressive emergence of a Gaussian distribution. They can also provide a numerical algorithm to determine maximum entropy states or minimum energy states.

Mechanisms for magnetic field reversals

We present a review of the different models that have been proposed to explain reversals of the magnetic field generated by a turbulent flow of an electrically conducting fluid (fluid dynamos). We then describe a simple mechanism that explains several features observed in palaeomagnetic records of the Earth's magnetic field, in numerical simulations and in a recent dynamo experiment. A similar model can also be used to understand reversals of large-scale flows that often develop on a turbulent background.

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<or=Re<or=120 and 0<or=Ra<or=30 000 .

Wave-driven dynamo action in spherical magnetohydrodynamic systems

Hydrodynamic and magnetohydrodynamic numerical studies of a mechanically forced two-vortex flow inside a sphere are reported. The simulations are performed in the intermediate regime between the laminar flow and developed turbulence, where a hydrodynamic instability is found to generate internal waves with a characteristic m=2 zonal wave number. It is shown that this time-periodic flow acts as a dynamo, although snapshots of the flow as well as the mean flow are not dynamos. The magnetic fields' growth rate exhibits resonance effects depending on the wave frequency. Furthermore, a cyclic self-killing and self-recovering dynamo based on the relative alignment of the velocity and magnetic fields is presented. The phenomena are explained in terms of a mixing of nonorthogonal eigenstates of the time-dependent linear operator of the magnetic induction equation. The potential relevance of this mechanism to dynamo experiments is discussed.

Transition by intermittency in granular matter: from discontinuous avalanches to continuous flow

We investigate, in the rotating drum configuration, the transition from the regime of discontinuous avalanches observed at low angular velocity to the regime of continuous flow observed at higher velocity. Instead of the hysteretic transition reported previously by Rajchenbach [Phys. Rev. Lett. 65, 2221 (1990)], with an apparent bistability of the two flow regimes in a range of drum velocities, we observe intermittency with spontaneous erratic switches from one regime to the other. Both scenarios of transition are recovered by a model dynamic equation for the avalanche flow with two sources of stochasticity: a Langevin noise during the avalanche flow and a distributed maximal stability angle at which avalanches start.

Transitions between turbulent states in rotating Rayleigh-Bénard convection

Weakly rotating turbulent Rayleigh-Bénard convection was studied experimentally and numerically. With increasing rotation and large enough Rayleigh number a supercritical bifurcation from a turbulent state with nearly rotation-independent heat transport to another with enhanced heat transfer is observed at a critical inverse Rossby number 1/Roc approximately 0.4. The strength of the large-scale convection roll is either enhanced or essentially unmodified depending on parameters for 1/Ro<1/Roc, but the strength increasingly diminishes beyond 1/Roc where it competes with Ekman vortices that cause vertical fluid transport and thus heat-transfer enhancement.