The dynamo effect is the most popular candidate to explain the non-primordial magnetic fields of astrophysical objects. Although many systematic studies of parameters have already been made to determine the different dynamical regimes explored by direct numerical geodynamo simulations, it is only recently that the regime corresponding to the outer core of the Earth characterized by a balance of forces between the Coriolis and Lorentz forces is accessible numerically. In most previous studies, the Lorentz force played a relatively minor role. For example, they have shown that a purely hydrodynamic parameter (the local Rossby number $Ro_\ell$ determines the stability domain of dynamos dominated by the axial dipole (dipolar dynamos). In this study, we show that this result cannot hold when the Lorentz force becomes dominant. We model turbulent geodynamo simulations with a strong Lorentz force by varying the important parameters over several orders of magnitude. This method enables us to question previous results and to argue on the applications of numerical dynamos in order to better understand the geodynamo problem. Strong dipolar fields considerably affect the kinetic energy distribution of convective motions which enables the maintenance of this field configuration. The relative importance of each force depends on the spatial length scale, whereas $Ro_\ell$ is a global output parameter which ignores the spatial dependency. We show that inertia does not induce a dipole collapse as long as the Lorentz and the Coriolis forces remain dominant at large length scales.

Stably stratified layers are present in stellar interiors (radiative zones) as well as planetary interiors - recent observations and theoretical studies of the Earth's magnetic field seem to indicate the presence of a thin, stably stratified layer at the top of the liquid outer core. We present direct numerical simulations of this region, which is modelled as an axisymmetric spherical Couette flow for a stably stratified fluid embedded in a dipolar magnetic field. For strong magnetic fields, a super-rotating shear layer, rotating nearly 30% faster than the imposed rotation rate difference between the inner convective dynamo region and the outer boundary, is generated in the stably stratified region. In the Earth context, and contrary to what was previously believed, we show that this super-rotation may extend towards the Earth magnetostrophic regime if the density stratification is sufficiently large. The corresponding differential rotation triggers magnetohydrodynamic instabilities and waves in the stratified region, which feature growth rates comparable to the observed timescale for geomagnetic secular variations and jerks. In the stellar context, we perform a linear analysis which shows that similar instabilities are likely to arise, and we argue that it may play a role in explaining the observed magnetic dichotomy among intermediate-mass stars.

We investigate the impact of a solid-body rotation $\Omega_0$ on the large-scale dynamics of an incompressible magnetohydrodynamic turbulent flow in presence of a background magnetic field $\bf B_0$ and at low Rossby number. Three-dimensional direct numerical simulations are performed in a periodic box, at unit magnetic Prandtl number and with a forcing at intermediate wavenumber $k_f=20$. When $\Omega_0$ is aligned with $\bf B_0$ (i.e. $\theta \equiv \widehat{\left(\Omega_{0}, \bf B_0 \right)} = 0$), inverse transfer is found for the magnetic spectrum at $k<k_f$. This transfer is stronger when the forcing excites preferentially right-handed (rather than left-handed) fluctuations; it is smaller when $\theta>0$ and becomes weak when $\theta \ge 35^o$. These properties are understood as the consequence of an inverse cascade of hybrid helicity which is an inviscid/ideal invariant of this system when $\theta=0$. Hybrid helicity emerges, therefore, as a key element for understanding rotating dynamos. Implication of these findings on the origin of the alignment of the magnetic dipole with the rotation axis in planets and stars is discussed.

During recent decades, data from space missions have provided strong evidence of deep liquid oceans underneath a thin outer icy crust on several moons of Jupiter, particularly Europa. But these observations have also raised many unanswered questions regarding the oceanic motions generated under the ice, or the mechanisms leading to the geological features observed on Europa. By means of direct numerical simulations of Europa's interior, we show here that Jupiter's magnetic field generates a retrograde oceanic jet at the equator, which may influence the global dynamics of Europa's ocean and contribute to the formation of some of its surface features by applying a unidirectional torque on Europa's ice shell.

We investigate the nature of the dynamo bifurcation in a configuration applicable to the Earth's liquid outer core. Numerical studies on the stability domain of dipolar magnetic fields found a dichotomy between non-reversing dipole-dominated dynamos and the reversing non-dipole-dominated multipolar solutions. We show that, by considering weak initial fields, the above transition disappears and is replaced by a region of bistability. Such a result was also observed in models with free-slip boundaries in which the geostrophic zonal flow can develop and participate to the dynamo mechanism for non-dipolar fields. We show that a similar process develops in no-slip models when viscous effects are reduced sufficiently. The following three regimes are distinguished: (i) Close to the onset of convection ($Ra_c$) with only the most critical convective mode (wave number) being present, dynamos set in supercritically and are dipole-dominated. (ii) in the range $3<Ra/Ra_c<Ra_c$, the bifurcations are subcritical and only dipole-dominated dynamos exist. (iii) $(Ra/Ra_c>10)$, the relative importance of zonal flows increases with $Ra$ in non-magnetic models. The field topology depends on the magnitude of the initial magnetic field. The dipolar branch has a subcritical behavior whereas the multipolar branch has a supercritical behavior. By approaching more realistic parameters, the extension of this bistable regime increases. A hysteretic behavior questions the common interpretation for geomagnetic reversals. Far above the dynamo threshold (by increasing the magnetic Prandtl number), Lorentz forces contribute to the first order force balance, as predicted for planetary dynamos. When $Ra$ is sufficiently high, dipolar fields affect significantly the flow speed, the flow structure and heat transfer which is reduced by the Lorentz force regardless of the field strength.

Recent interferometric data have been used to constrain the brightness distribution at the surface of nearby stars, in particular the so-called gravity darkening that makes fast rotating stars brighter at their poles than at their equator. However, good models of gravity darkening are missing for stars that posses a convective envelope. In order to better understand how rotation affects the heat transfer in stellar convective envelopes, we focus on the heat flux distribution in latitude at the outer surface of numerical models. We carry out a systematic parameter study of three-dimensional, direct numerical simulations of anelastic convection in rotating spherical shells. Restricting our investigations to hydrodynamical models with a thermal Prandtl number fixed to unity, we consider both thick and thin (solar-like) shells, and vary the stratification over three orders of magnitude. We measure the heat transfer efficiency in terms of the Nusselt number, defined as the output luminosity normalised by the conductive state luminosity. We report diverse Nusselt number profiles in latitude, ranging from brighter (usually at the onset of convection) to darker equator and uniform profiles. We find that the variations of the surface brightness are mainly controlled by the surface value of the local Rossby number: when the Coriolis force dominates the dynamics, the heat flux is weakened in the equatorial region by the zonal wind and enhanced at the poles by convective motions inside the tangent cylinder. In the presence of a strong background density stratification however, as expected in real stars, the increase of the local Rossby number in the outer layers leads to uniformisation of the surface heat flux distribution.

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.

Numerical modelling of convection driven dynamos in the Boussinesq approximation revealed fundamental characteristics of the dynamo-generated magnetic fields and the fluid flow. Because these results were obtained for an incompressible fluid, their validity for gas planets and stars remains to be assessed. A common approach is to take some density stratification into account with the so-called anelastic approximation. The validity of previous results obtained in the Boussinesq approximation is tested for anelastic models. We point out and explain specific differences between both types of models, in particular with respect to the field geometry and the field strength, but we also compare scaling laws for the velocity amplitude, the magnetic dissipation time, and the convective heat flux. Our investigation is based on a systematic parameter study of spherical dynamo models in the anelastic approximation. We make use of a recently developed numerical solver and provide results for the test cases of the anelastic dynamo benchmark. The dichotomy of dipolar and multipolar dynamos identified in Boussinesq simulations is also present in our sample of anelastic models. Dipolar models require that the typical length scale of convection is an order of magnitude larger than the Rossby radius. However, the distinction between both classes of models is somewhat less explicit than in previous studies. This is mainly due to two reasons: we found a number of models with a considerable equatorial dipole contribution and an intermediate overall dipole field strength. Furthermore, a large density stratification may hamper the generation of dipole dominated magnetic fields. Previously proposed scaling laws, such as those for the field strength, are similarly applicable to anelastic models. It is not clear, however, if this consistency necessarily implies similar dynamo processes in both settings.

Context. Three-dimensional spherical dynamo simulations carried out within the framework of the anelastic approximation have revealed that the established distinction between dipolar and multipolar dynamos tends to be less clear than it was in Boussinesq studies. This result was first interpreted as a direct consequence of the existence of a larger number of models with a high equatorial dipole contribution, together with an intermediate dipole field strength. However, this finding has not been clearly related to specific changes that would have been introduced by the use of the anelastic approximation. Aims. In this paper, we primarily focus on the effects of choosing of a different mass distribution. Indeed, it is likely to have as large consequences as would taking a stratified reference state into account, especially when comparing our results to previous Boussinesq studies. Methods. Our investigation is based on a systematic parameter study of weakly stratified anelastic dynamo models. Results. We show that the tendencies highlighted in previous anelastic dynamo simulations are already present in the Boussinesq limit. Thus they cannot be systematically related to anelastic effects. Actually, a central mass distribution can result in changes in the magnetic field topology that are mainly due to the concentration of convective cells close to the inner sphere.

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.

Magnetic fields of low-mass stars and planets are thought to originate from self-excited dynamo action in their convective interiors. Observations reveal a variety of field topologies ranging from large-scale, axial dipole to more structured magnetic fields. In this article, we investigate more than 70 three-dimensional, self-consistent dynamo models obtained by direct numerical simulations. The control parameters, the aspect ratio and the mechanical boundary conditions have been varied to build up this sample of models. Both, strongly dipolar and multipolar models have been obtained. We show that these dynamo regimes can in general be distinguished by the ratio of a typical convective length scale to the Rossby radius. Models with a predominantly dipolar magnetic field were obtained, if the convective length scale is at least an order of magnitude larger than the Rossby radius. Moreover, we highlight the role of the strong shear associated with the geostrophic zonal flow for models with stress-free boundary conditions. In this case, the above transition disappears and is replaced by a region of bistability for which dipolar and multipolar dynamos co-exist. We interpret our results in terms of dynamo eigenmodes using the so-called test-field method. We can thus show that models in the dipolar regime are characterized by an isolated 'single mode'. Competing overtones become significant as the boundary to multipolar dynamos is approached. We discuss how these findings relate to previous models and to observations.

Context: Large-scale magnetic fields resulting from hydromagnetic dynamo action may differ substantially in their time dependence. Cyclic field variations, characteristic for the solar magnetic field, are often explained by an important omega-effect, i.e. by the stretching of field lines due to strong differential rotation. Aims: The dynamo mechanism of a convective, oscillatory dynamo model is investigated. Methods: We solve the MHD-equations for a conducting Boussinesq fluid in a rotating spherical shell. For a resulting oscillatory model, dynamo coefficients have been computed with the help of the so-called test-field method. Subsequently, these coefficients have been used in a mean-field calculation in order to explore the underlying dynamo mechanism. Results: Although the rather strong differential rotation present in this model influences the magnetic field, the omega-effect alone is not responsible for its cyclic time variation. If the omega-effect is suppressed, the resulting alpha^2-dynamo remains oscillatory. Surprisingly, the corresponding alpha-omega dynamo leads to a non-oscillatory magnetic field. Conclusions: The assumption of an alpha-omega mechanism does not explain the occurrence of magnetic cycles satisfactorily.

We show that a simple, modified version of the Magnetorotational Instability (MRI) can develop in the outer liquid core of the Earth, in the presence of a background shear. It requires either thermal wind, or a primary instability, such as convection, to drive a weak differential rotation within the core. The force balance in the Earth's core is very unlike classical astrophysical applications of the MRI (such as gaseous disks around stars). Here, the weak differential rotation in the Earth core yields an instability by its constructive interaction with the planet's much larger rotation rate. The resulting destabilising mechanism is just strong enough to counteract stabilizing resistive effects, and produce growth on geophysically interesting timescales. We give a simple physical explanation of the instability, and show that it relies on a force balance appropriate to the Earth's core, known as magnetostrophic balance.