It is demonstrated that the azimuthal magnetorotational instability (AMRI) also works with radially increasing rotation rates contrary to the standard magnetorotational instability for axial fields which requires negative shear. The stability against non-axisymmetric perturbations of a conducting Taylor–Couette flow with positive shear under the influence of a toroidal magnetic field is considered if the background field between the cylinders is current free. For small magnetic Prandtl number $Pm\rightarrow 0$ the curves of neutral stability converge in the (Hartmann number,Reynolds number) plane approximating the stability curve obtained in the inductionless limit $Pm=0$. The numerical solutions for $Pm=0$ indicate the existence of a lower limit of the shear rate. For large $Pm$ the curves scale with the magnetic Reynolds number of the outer cylinder but the flow is always stable for magnetic Prandtl number unity as is typical for double-diffusive instabilities. We are particularly interested to know the minimum Hartmann number for neutral stability. For models with resting or almost resting inner cylinder and with perfectly conducting cylinder material the minimum Hartmann number occurs for a radius ratio of $r_{\text{in}}=0.9$. The corresponding critical Reynolds numbers are smaller than $10^{4}$.

A mean-field theory of differential rotation in a density stratified turbulent convection has been developed. This theory is based on the combined effects of the turbulent heat flux and anisotropy of turbulent convection on the Reynolds stress. A coupled system of dynamical budget equations consisting in the equations for the Reynolds stress, the entropy fluctuations and the turbulent heat flux has been solved. To close the system of these equations, the spectral $\unicode[STIX]{x1D70F}$ approach, which is valid for large Reynolds and Péclet numbers, has been applied. The adopted model of the background turbulent convection takes into account an increase of the turbulence anisotropy and a decrease of the turbulent correlation time with the rotation rate. This theory yields the radial profile of the differential rotation which is in agreement with that for the solar differential rotation.

The homogeneous dynamo effect is at the root of magnetic field generation in cosmic bodies, including planets, stars and galaxies. While the underlying theory had increasingly flourished since the middle of the 20th century, hydromagnetic dynamos were not realized in the laboratory until 1999. On 11 November 1999, this situation changed with the first observation of a kinematic dynamo in the Riga experiment. Since that time, a series of experimental campaigns has provided a wealth of data on the kinematic and the saturated regime. This paper is intended to give a comprehensive survey about these experiments, to summarize their main results and to compare them with numerical simulations.

Mean field electrodynamics (MFE) facilitates practical modelling of secular, large scale properties of astrophysical or laboratory systems with fluctuations. Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages. Often however, real systems do not exhibit such scale separation. This raises two questions: (I) What are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) How precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel. Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise. The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work.

We discuss a mean-field theory of the generation of large-scale vorticity in a rotating density stratified developed turbulence with inhomogeneous kinetic helicity. We show that the large-scale non-uniform flow is produced due to either a combined action of a density stratified rotating turbulence and uniform kinetic helicity or a combined effect of a rotating incompressible turbulence and inhomogeneous kinetic helicity. These effects result in the formation of a large-scale shear, and in turn its interaction with the small-scale turbulence causes an excitation of the large-scale instability (known as a vorticity dynamo) due to a combined effect of the large-scale shear and Reynolds stress-induced generation of the mean vorticity. The latter is due to the effect of large-scale shear on the Reynolds stress. A fast rotation suppresses this large-scale instability.

In most numerical simulations of the Earth’s core the dynamo is located outside the tangent cylinder and, in a zero-order sense, takes the form of a classical $\unicode[STIX]{x1D6FC}^{2}$ dynamo. Such a dynamo usually requires a distribution of helicity, $h$, which is asymmetric about the equator and in the simulations it is observed that, outside the tangent cylinder, the helicity is predominantly negative in the north and positive in the south. If we are to extrapolate the results of these simulations to the planets, we must understand how this asymmetry in helicity is established and ask if the same mechanism is likely to operate in a planet. In some of the early numerical dynamos, which were too viscous by a factor of at least $10^{9}$, as measured by the Ekman number, the asymmetric helicity distribution was attributed to Ekman pumping. However, Ekman pumping plays much less of a role in more recent, and less viscous, numerical dynamos, and almost certainly plays no significant role in the core of a planet. So the question remains: what establishes the asymmetric helicity distribution in the simulations and is this mechanism likely to carry over to planetary cores? In this paper we review the evidence that planetary dynamos, and their numerical analogues, might be maintained by helical waves, especially inertial waves, excited in and around the equatorial regions. This cartoon arises from the observation that there tends to be a statistical bias in the buoyancy flux towards the equatorial regions, and so waves are preferentially excited there. Moreover, upward (downward) propagating inertial waves carry negative (positive) helicity, which leads naturally to a segregation in $h$.

Using a variational optimization method we find the smallest flow perturbations that can trigger kinematic dynamo action in Kolmogorov flow. In comparison to previous work, a second-order mean field dynamo model is used to track down the optimal dynamos in the high magnetic Reynolds number limit ($Rm$). The magnitude of minimal perturbation flows decays inversely proportional to the magnetic Reynolds number. We reveal the asymptotic high-$Rm$ structure of the optimal flow perturbation and the magnetic eigenmode. We identify the optimal dynamo as of $\unicode[STIX]{x1D6FC}{-}\unicode[STIX]{x1D6FA}$ type, with magnetic fluctuations that localize on a critical layer.

We apply a nonlinear mean-field dynamo model which includes a budget equation for the dynamics of Wolf numbers to predict solar activity. This dynamo model takes into account the algebraic and dynamic nonlinearities of the $\unicode[STIX]{x1D6FC}$ effect, where the equation for the dynamic nonlinearity is derived from the conservation law for the magnetic helicity. The budget equation for the evolution of the Wolf number is based on a formation mechanism of sunspots related to the negative effective magnetic pressure instability. This instability redistributes the magnetic flux produced by the mean-field dynamo. To predict solar activity on the time scale of one month we use a method based on a combination of the numerical solution of the nonlinear mean-field dynamo equations and the artificial neural network. A comparison of the results of the prediction of the solar activity with the observed Wolf numbers demonstrates a good agreement between the forecast and observations.

Mean-field dynamo equations are addressed with the aid of the path integral method. The evolution of magnetic field is treated as a three-dimensional Wiener random process, and the mean magnetic-field equations are obtained with the Wiener integrals taken over all the trajectories of the fluid particles. The form of the equations is just the same as the conventional mean-field equations, but here the equations are derived with the velocity field realisation affected by the force exerted by the magnetic field. In this sense, we derive nonlinear dynamo equations.

Stellar radiative zones are typically assumed to be motionless in standard models of stellar structure but there is sound theoretical and observational evidence that this cannot be the case. We investigate by direct numerical simulations a three-dimensional and time-dependent model of stellar radiation zones consisting of an electrically conductive and stably stratified anelastic fluid confined to a rotating spherical shell and driven by a baroclinic torque. As the baroclinic driving is gradually increased a sequence of transitions from an axisymmetric and equatorially symmetric time-independent flow to flows with a strong poloidal component and lesser symmetry are found. It is shown that all flow regimes characterised by significant non-axisymmetric components are capable of generating a self-sustained magnetic field. As the value of the Prandtl number is decreased and the value of the Ekman number is decreased, flows become strongly time dependent with progressively complex spatial structure and dynamos can be generated at lower values of the magnetic Prandtl number.

We discuss the effect of different types of fluctuations on dynamos generated in the limit of scale separation. We first recall that the magnetic field observed in the VKS (von Karman flow of liquid sodium) experiment is not the one that would be generated by the mean flow alone and that smaller scale turbulent fluctuations therefore play an important role. We then consider how velocity fluctuations affect the dynamo threshold in the framework of mean-field magnetohydrodynamics. We show that the detrimental effect of turbulent fluctuations observed with many flows disappears in the case of helical flows with scale separation. We also find that fluctuations of the electrical conductivity of the fluid, for instance related to temperature fluctuations in convective flows, provide an efficient mechanism for dynamo action. Finally, we conclude by describing an experimental configuration that could be used to test the validity of mean-field magnetohydrodynamics in strongly turbulent flows.

A dynamo driven by motions unaffected by viscous forces is termed ‘magnetostrophic’, but cannot be found through today’s numerical simulations, which require substantial viscosity to stabilize solutions of the full magnetohydrodynamic (MHD) dynamo equations. By using an alternative numerical technique, proposed by Taylor (Proc. R. Soc. Lond. A, vol. 274, 1963, pp. 274–283), we recently obtained the first magnetostrophic dynamo solutions ever derived (Wu & Roberts, Geophys. Astrophys. Fluid Dyn., vol. 109, 2014, pp. 84–110). These were axisymmetric and of mean-field type. In an earlier paper (Roberts & Wu, Geophys. Astrophys. Fluid Dyn., vol. 108, 2014, pp. 696–715), we proposed an extension of Taylor’s method. Here we explore its numerical implications, comparing them to the consequences of Taylor’s original proposal. One of the differences between the two approaches is that our modification retains torsional waves but Taylor’s theory does not. A more important difference is that our extension of Taylor’s method is, for reasons presented here, the most general possible that does not suffer from the limitations imposed by viscosity on today’s numerical simulations.

The predictions of mean-field electrodynamics can now be probed using direct numerical simulations of random flows and magnetic fields. When modelling astrophysical magnetohydrodynamics, it is important to verify that such simulations are in agreement with observations. One of the main challenges in this area is to identify robust quantitative measures to compare structures found in simulations with those inferred from astrophysical observations. A similar challenge is to compare quantitatively results from different simulations. Topological data analysis offers a range of techniques, including the Betti numbers and persistence diagrams, that can be used to facilitate such a comparison. After describing these tools, we first apply them to synthetic random fields and demonstrate that, when the data are standardized in a straightforward manner, some topological measures are insensitive to either large-scale trends or the resolution of the data. Focusing upon one particular astrophysical example, we apply topological data analysis to H i observations of the turbulent interstellar medium (ISM) in the Milky Way and to recent magnetohydrodynamic simulations of the random, strongly compressible ISM. We stress that these topological techniques are generic and could be applied to any complex, multi-dimensional random field.

Recent advances in mean-field theory are reviewed and applications to the Sun, late-type stars, accretion disks, galaxies and the early Universe are discussed. We focus particularly on aspects of spatio-temporal non-locality, which provided some of the main new qualitative and quantitative insights that emerged from applying the test-field method to magnetic fields of different length and time scales. We also review the status of nonlinear quenching and the relation to magnetic helicity, which is an important observational diagnostic of modern solar dynamo theory. Both solar and some stellar dynamos seem to operate in an intermediate regime that has not yet been possible to model successfully. This regime is bracketed by antisolar-like differential rotation on one end and stellar activity cycles belonging to the superactive stars on the other. The difficulty in modelling this regime may be related to shortcomings in simulating solar/stellar convection. On galactic and extragalactic length scales, the observational constraints on dynamo theory are still less stringent and more uncertain, but recent advances both in theory and observations suggest that more conclusive comparisons may soon be possible also here. The possibility of inversely cascading magnetic helicity in the early Universe is particularly exciting in explaining the recently observed lower limits of magnetic fields on cosmological length scales. Such magnetic fields may be helical with the same sign of magnetic helicity throughout the entire Universe. This would be a manifestation of parity breaking.

The aim of this short paper is simply to demonstrate the effect of interactions between distinct waves induced by an oscillatory forcing on generation of the large-scale electromotive force (EMF) in a highly electrically conducting plasma in the dynamical regime, i.e. when the back reaction of the Lorentz force on the flow is included. The mean EMF is induced by waves with distinct but close phase speeds, which are called ‘beating waves’ by analogy with the acoustic effect of ‘beat’, in the presence of a locally uniform seed magnetic field. These waves may have either positive or negative helicity, and it is supposed that waves of a single sign of helicity are preferentially excited by a symmetry-breaking mechanism. The formula for the mean EMF in a highly conducting plasma is derived, conditions relevant to those of the early universe before and during galaxy formation. Within the scope of the weak magnetohydrodynamic turbulence dominated by the linear waves an $\unicode[STIX]{x1D6FC}$-effect, that is generation of the large-scale EMF is studied, which leads to amplification of the seed mean magnetic field. Possible subsequent equilibration due to quenching resulting from the presence of the Lorentz force is briefly discussed.

The quasi-geostrophic dynamo model (QGDM) is a multiscale, fully nonlinear Cartesian dynamo model that is valid in the asymptotic limit of low Rossby number. In the additional limit of small magnetic Prandtl number investigated here, the QGDM is a self-consistent, asymptotically exact form of an $\unicode[STIX]{x1D6FC}^{2}$ large-scale dynamo. This article explores methods for simulating the multiscale QGDM and investigates how convection is altered by the magnetic field in the planetary regime of small Rossby number and small magnetic Prandtl number. At present, this combination is beyond the reach of direct numerical simulations. We use a simplified class of solutions whose horizontal structure is restricted to a periodic hexagonal lattice characterized by a single horizontal wavenumber (single mode). In contrast with previous kinematic investigations of the QGDM, the Lorentz force is included to study saturated, self-consistent dynamos. Two methodologies are used to assess handling of the multiple time scales of the QGDM: a stiff, common-in-time approach where all time scales are converted to a single time variable and a heterogeneous multiscale modelling approach employing fast time averaging on the Reynolds, magnetic and buoyancy eddy fluxes that feed back onto the slow scales. These strategies produce consistent results and each illustrates self-similar dynamics as the time-averaging window is increased. The properties of the convection are significantly altered by the dynamo-generated magnetic field. All solutions show a decrease in the overall heat transfer efficiency as compared to non-magnetic convection, suggesting that a change in length scale or flow planform plays a critical role in the enhanced heat transfer efficiency observed in previous dynamo studies. All dynamo solutions show a trend of increasing ohmic dissipation relative to viscous dissipation as the buoyancy forcing is increased.

The theory of mean field electrodynamics, now celebrating its fiftieth birthday, has had a profound influence on our modelling of cosmical dynamos, greatly enhancing our understanding of how such dynamos may operate. Here I discuss some of its undoubted triumphs, but also some of the problems that can arise in a mean field approach to dynamos in fluids (or plasmas) that are both highly turbulent and also extremely good electrical conductors, as found in all astrophysical settings.

Fully compressible magnetohydrodynamic (MHD) turbulence is investigated in the framework of the multiple-scale direct-interaction approximation. With the aid of the propagators (correlation and Green’s functions), fluctuating fields are solved, and turbulent correlations are estimated in highly compressible turbulence. We focus on the expression of the turbulent electromotive force (EMF). Obliqueness between the mean magnetic field and the mean-density gradient, the mean internal density gradient and the non-equilibrium mean velocity contributes to the EMF in the presence of the density variance, which is ubiquitous in turbulence in strongly variable density flows such as the shock-front region. This density-variance effect is expected to locally enhance the turbulence intensity across the shock front, leading to a fast reconnection.

We develop a mean-field theory of compressibility effects in turbulent magnetohydrodynamics and passive scalar transport using the quasi-linear approximation and the spectral $\unicode[STIX]{x1D70F}$-approach. We find that compressibility decreases the $\unicode[STIX]{x1D6FC}$ effect and the turbulent magnetic diffusivity both at small and large magnetic Reynolds numbers, $Rm$. Similarly, compressibility decreases the turbulent diffusivity for passive scalars both at small and large Péclet numbers, $Pe$. On the other hand, compressibility does not affect the effective pumping velocity of the magnetic field for large $Rm$, but it decreases it for small $Rm$. Density stratification causes turbulent pumping of passive scalars, but it is found to become weaker with increasing compressibility. No such pumping effect exists for magnetic fields. However, compressibility results in a new passive scalar pumping effect from regions of low to high turbulent intensity both for small and large Péclet numbers. It can be interpreted as compressible turbophoresis of non-inertial particles and gaseous admixtures, while the classical turbophoresis effect exists only for inertial particles and causes them to be pumped to regions with lower turbulent intensity.

We explore the growth of large-scale magnetic fields in a shear flow, due to helicity fluctuations with a finite correlation time, through a study of the Kraichnan–Moffatt model of zero-mean stochastic fluctuations of the $\unicode[STIX]{x1D6FC}$ parameter of dynamo theory. We derive a linear integro-differential equation for the evolution of the large-scale magnetic field, using the first-order smoothing approximation and the Galilean invariance of the $\unicode[STIX]{x1D6FC}$-statistics. This enables construction of a model that is non-perturbative in the shearing rate $S$ and the $\unicode[STIX]{x1D6FC}$-correlation time $\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D6FC}}$. After a brief review of the salient features of the exactly solvable white-noise limit, we consider the case of small but non-zero $\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D6FC}}$. When the large-scale magnetic field varies slowly, the evolution is governed by a partial differential equation. We present modal solutions and conditions for the exponential growth rate of the large-scale magnetic field, whose drivers are the Kraichnan diffusivity, Moffatt drift, shear and a non-zero correlation time. Of particular interest is dynamo action when the $\unicode[STIX]{x1D6FC}$-fluctuations are weak; i.e. when the Kraichnan diffusivity is positive. We show that in the absence of Moffatt drift, shear does not give rise to growing solutions. But shear and Moffatt drift acting together can drive large-scale dynamo action with growth rate $\unicode[STIX]{x1D6FE}\propto |S|$.