Search

We present dynamos computed using a hybrid QG-3D numerical scheme in a thick spherical shell geometry. Our model is based on a quasi-geostrophic convection code extended with a 3-D treatment of heat transport and magnetic induction. We find a collection of self-sustained, multipolar, weak field dynamos with magnetic energy one or two orders of magnitude lower than the kinetic energy. The poloidal magnetic energy is weak and, by construction, there is a lack of equatorially antisymmetric components in the Buoyancy and Lorentz forces. This leads to configurations where the velocity field is only weakly impacted by the magnetic field, similar to dynamos found in 3-D simulations where zonal flows and the Ω-effect dominate. The time-dependence of these dynamos is characterized by quasi-periodic oscillations that we attribute to dynamo waves. The QG-3D dynamos found so far are not Earth-like. The inability of our setup to produce strong, dipole-dominated, magnetic fields likely points to a missing ingredient in our QG flows, and a related lack of helicity and α -effect. The models presented here may be more relevant for studying stellar dynamos where zonal flows are known to dominate. This study was carried out at modest control parameters, however moving to lower Ekman numbers, when smaller values of both the magnetic and hydrodynamic Prandtl numbers can be of interest, our approach will be able to gain in efficiency by using relatively coarse grids for the 3-D magnetic and temperature fields and a finer grid for the QG velocity field. [ABSTRACT FROM AUTHOR]

In fluid dynamics, the scaling behaviour of flow length scales is commonly used to infer the governing force balance of a system. The key to a successful approach is to measure length scales that are simultaneously representative of the energy contained in the solution (energetically relevant) and also indicative of the established force balance (dynamically relevant). In the case of numerical simulations of rotating convection and magnetohydrodynamic dynamos in spherical shells, it has remained difficult to measure length scales that are both energetically and dynamically relevant, a situation that has led to conflicting interpretations, and sometimes misrepresentations of the underlying force balance. By analysing an extensive set of magnetic and non-magnetic models, we focus on two length scales that achieve both energetic and dynamical relevance. The first one is the peak of the poloidal kinetic energy spectrum, which we successfully compare to crossover points on spectral representations of the force balance. In most dynamo models, this result confirms that the dominant length scale of the system is controlled by a previously proposed quasi-geostrophic (QG-) MAC (Magneto-Archimedean-Coriolis) balance. In non-magnetic convection models, the analysis generally favours a QG-CIA (Coriolis-Inertia-Archimedean) balance. Viscosity, which is typically a minor contributor to the force balance, does not control the dominant length scale at high convective supercriticalities in the non-magnetic case, and in the dynamo case, once the generated magnetic energy largely exceeds the kinetic energy. In dynamo models, we introduce a second energetically relevant length scale associated with the loss of axial invariance in the flow. We again relate this length scale to another crossover point in scale-dependent force balance diagrams, which marks the transition between large-scale geostrophy (the equilibrium of Coriolis and pressure forces) and small-scale magnetostrophy, where the Lorentz force overtakes the Coriolis force. Scaling analysis of these two energetically and dynamically relevant length scales suggests that the Earth's dynamo is controlled by a QG-MAC balance at a dominant scale of about |$200 \, \mathrm{km}$|⁠ , while magnetostrophic effects are deferred to scales smaller than |$50 \, \mathrm{km}$|⁠. [ABSTRACT FROM AUTHOR]

We study rotating thermal convection in spherical shells as prototype for flow in the cores of terrestrial planets, gas planets or in stars. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from 10-3 to 10-5. The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh-Bénard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Simulations with convection zones of different thickness show that the zonal flows become less energetic in thicker shells, and, therefore, their effect on heat-transfer efficiency decreases. Furthermore, we explore the influence of the magnetic field on the nonzonal flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at E = 10-5 magnetic field significantly affects the convection and a maximum gain of about 30 per cent (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at E = 10-5 is probably in a 'magnetostrophic' regime. [ABSTRACT FROM AUTHOR]