Predictive Scaling Laws for Spherical Rotating Dynamos

State of the art numerical models of the Geodynamo are still performed in a parameter regime extremely remote from the values relevant to the physics of the Earth's core. In order to establish a connection between dynamo modeling and the geophysical motivation, {it is necessary to use} scaling laws. Such scaling laws establish the dependence of essential quantities (such as the magnetic field strength) on measured or controlled quantities. They allow for a direct confrontation of advanced models with geophysical {constraints}. (...) We show that previous empirical scaling laws for the magnetic field strength essentially reflect the statistical balance between energy production and dissipation for saturated dynamos. Such power based scaling laws are thus necessarily valid for any dynamo in statistical equilibrium and applicable to any numerical model, irrespectively of the dynamo mechanism. We show that direct numerical fits can provide contradictory results owing to biases in the parameters space covered in the numerics and to the role of a priori hypothesis on the fraction of ohmic dissipation. We introduce predictive scaling laws, i.e. relations involving input parameters of the governing equations only. We guide our reasoning on physical considerations. We show that our predictive scaling laws can properly describe the numerical database and reflect the dominant forces balance at work in these numerical simulations. We highlight the dependence of the magnetic field strength on the rotation rate. Finally, our results stress that available numerical models operate in a viscous dynamical regime, which is not relevant to the Earth's core

The dynamo bifurcation in rotating spherical shells

We investigate the nature of the dynamo bifurcation in a configuration applicable to the Earth's liquid outer core, i.e. in a rotating spherical shell with thermally driven motions. We show that the nature of the bifurcation, which can be either supercritical or subcritical or even take the form of isola (or detached lobes) strongly depends on the parameters. This dependence is described in a range of parameters numerically accessible (which unfortunately remains remote from geophysical application), and we show how the magnetic Prandtl number and the Ekman number control these transitions.Comment: 16 pages, 14 figure

Toward an asymptotic behaviour of the ABC dynamo

The ABC flow was originally introduced by Arnol'd to investigate Lagrangian chaos. It soon became the prototype example to illustrate magnetic-field amplification via fast dynamo action, i.e. dynamo action exhibiting magnetic-field amplification on a typical timescale independent of the electrical resistivity of the medium. Even though this flow is the most classical example for this important class of dynamos (with application to large-scale astrophysical objects), it was recently pointed out (Bouya Isma\"el and Dormy Emmanuel, Phys. Fluids, 25 (2013) 037103) that the fast dynamo nature of this flow was unclear, as the growth rate still depended on the magnetic Reynolds number at the largest values available so far $(\text{Rm} = 25000)$ . Using state-of-the-art high-performance computing, we present high-resolution simulations (up to 40963) and extend the value of $\text{Rm}$ up to $5\cdot10^5$ . Interestingly, even at these huge values, the growth rate of the leading eigenmode still depends on the controlling parameter and an asymptotic regime is not reached yet. We show that the maximum growth rate is a decreasing function of $\text{Rm}$ for the largest values of $\text{Rm}$ we could achieve (as anticipated in the above-mentioned paper). Slowly damped oscillations might indicate either a new mode crossing or that the system is approaching the limit of an essential spectrum

Transition between viscous dipolar and inertial multipolar dynamos

We investigate the transition from steady dipolar to reversing multipolar dynamos. The Earth has been argued to lie close to this transition, which could offer a scenario for geomagnetic reversals. We show that the transition between dipolar and multipolar dynamos is characterized by a three terms balance (as opposed to the usually assumed two terms balance), which involves the non-gradient parts of inertial, viscous and Coriolis forces. We introduce from this equilibrium the sole parameter ${{\rm Ro}}\,{{\rm E}}^{-1/3} \equiv {{\rm Re}}\,{{\rm E}}^{2/3}$, which accurately describes the transition for a wide database of 132 fully three dimensional direct numerical simulations of spherical rotating dynamos (courtesy of U. Christensen). This resolves earlier contradictions in the literature on the relevant two,terms balance at the transition. Considering only a two terms balance between the non-gradient part of the Coriolis force and of inertial forces, provides the classical ${{\rm Ro}}/{\ell_u}$ (Christensen and Aubert, 2006). This transition can be equivalently described by ${{\rm Re}} \, {\ell^{2}_u}$, which corresponds to the two terms balance between the non-gradient part of inertial forces and viscous forces (Soderlund {\it et al.}, 2012).Comment: 14 pages, 4 figure

Mechanisms of Planetary and Stellar Dynamos

We review some of the recent progress on modeling planetary and stellar dynamos. Particular attention is given to the dynamo mechanisms and the resulting properties of the field. We present direct numerical simulations using a simple Boussinesq model. These simulations are interpreted using the classical mean-field formalism. We investigate the transition from steady dipolar to multipolar dynamo waves solutions varying different control parameters, and discuss the relevance to stellar magnetic fields. We show that owing to the role of the strong zonal flow, this transition is hysteretic. In the presence of stress-free boundary conditions, the bistability extends over a wide range of parameters.Comment: Proceedings of IAUS 294 "Solar and Astrophysical Dynamos and Magnetic Activity" Editors A.G. Kosovichev, E.M. de Gouveia Dal Pino, & Y.Yan, Cambridge University Press, to appear (2013

Dipolar dynamos in stratified systems

Observations of low-mass stars reveal a variety of magnetic field topologies ranging from large-scale, axial dipoles to more complex magnetic fields. At the same time, three-dimensional spherical simulations of convectively driven dynamos reproduce a similar diversity, which is commonly obtained either with Boussinesq models or with more realistic models based on the anelastic approximation, which take into account the variation of the density with depth throughout the convection zone. Nevertheless, a conclusion from different anelastic studies is that dipolar solutions seem more difficult to obtain as soon as substantial stratifications are considered. In this paper, we aim at clarifying this point by investigating in more detail the influence of the density stratification on dipolar dynamos. To that end, we rely on a systematic parameter study that allows us to clearly follow the evolution of the stability domain of the dipolar branch as the density stratification is increased. The impact of the density stratification both on the dynamo onset and the dipole collapse is discussed and compared to previous Boussinesq results. Furthermore, our study indicates that the loss of the dipolar branch does not ensue from a specific modification of the dynamo mechanisms related to the background stratification, but could instead result from a bias as our observations naturally favour a certain domain in the parameter space characterized by moderate values of the Ekman number, owing to current computational limitations. Moreover, we also show that the critical magnetic Reynolds number of the dipolar branch is scarcely modified by the increase of the density stratification, which provides an important insight into the global understanding of the impact of the density stratification on the stability domain of the dipolar dynamo branch

Spin-down in a rapidly rotating cylinder container with mixed rigid and stress-free boundary conditions

A comprehensive study of the classical linear spin-down of a constant density viscous fluid (kinematic viscosity \nu) rotating rapidly (angular velocity \Omega) inside an axisymmetric cylindrical container (radius L, height H) with rigid boundaries, that follows the instantaneous small change in the boundary angular velocity at small Ekman number $E=\nu/H^2\Omega \ll 1$, was provided by Greenspan & Howard (1963). $E^{1/2}$-Ekman layers form quickly triggering inertial waves together with the dominant spin-down of the quasi-geostrophic (QG) interior flow on the $O(E^{-1/2}\Omega^{-1})$ time-scale. On the longer lateral viscous diffusion time-scale $O(L^2/\nu)$, the QG-flow responds to the $E^{1/3}$-side-wall shear-layers. In our variant the side-wall and top boundaries are stress-free; a setup motivated by the study of isolated atmospheric structures, such as tropical cyclones, or tornadoes. Relative to the unbounded plane layer case, spin-down is reduced (enhanced) by the presence of a slippery (rigid) side-wall. This is evinced by the QG-angular velocity, \omega*, evolution on the O(L^2/\nu) time-scale: Spatially, \omega* increases (decreases) outwards from the axis for a slippery (rigid) side-wall; temporally, the long-time ($\gg L^2/\nu)$ behaviour is dominated by an eigensolution with a decay rate slightly slower (faster) than that for an unbounded layer. In our slippery side-wall case, the $E^{1/2} \times E^{1/2}$ corner region that forms at the side-wall intersection with the rigid base is responsible for a $\ln E$ singularity within the $E^{1/3}$-layer causing our asymptotics to apply only at values of E far smaller than can be reached by our Direct Numerical Simulation (DNS) of the entire spin-down process. Instead, we solve the $E^{1/3}$-boundary-layer equations for given E numerically. Our hybrid asymptotic-numerical approach yields results in excellent agreement with our DNS.Comment: 33 pages, 10 figure