Search

Stably-stratified layers may be present at the top of the electrically-conducting fluid layers of many planets either because the temperature gradient is locally subadiabatic or because a stable composition gradient is maintained by the segregation of chemical elements. Here we study the double-diffusive processes taking place in such a stable layer, considering the case of Mercury's core where the temperature gradient is stable but the composition gradient is unstable over a 800km-thick layer. The large difference in the molecular diffusivities leads to the development of buoyancy-driven instabilities that drive radial flows known as fingering convection. We model fingering convection using hydrodynamical simulations in a rotating spherical shell and varying the rotation rate and the stratification strength. For small Rayleigh numbers (i.e. weak background temperature and composition gradients), fingering convection takes the form of columnar flows aligned with the rotation axis and with an azimuthal size comparable with the layer thickness. For larger Rayleigh numbers, the flows retain a columnar structure but the azimuthal size is drastically reduced leading to thin sheet-like structures that are elongated in the meridional direction. The azimuthal length decreases when the thermal stratification increases, following closely the scaling law expected from the linear non-rotating planar theory (Stern, 1960). We find that the radial flows always remain laminar with local Reynolds number of order 1-10. Equatorially-symmetric zonal flows form due to latitudinal variations of the axisymmetric composition. The zonal velocity exceeds the non-axisymmetric velocities at the largest Rayleigh numbers. We discuss plausible implications for planetary magnetic fields., Comment: accepted for publication in JGR: Planets

We study dynamo action using numerical simulations of planar Boussinesq convection at rapid rotation (low Ekman numbers, Ek), focusing on subcritical dynamo action in which the dynamo is sustained for Rayleigh numbers, Ra, below the critical Rayleigh number for the onset of nonmagnetic convection, Ra$_c$. These solutions are found by first investigating the supercritical regime, in which the dynamo is able to generate a large-scale magnetic field that significantly influences the convective motions, with an associated Elsasser number of order Ek$^{1/3}$. Subcritical solutions are then found by tracking this solution branch into the subcritical regime, taking a supercritical solution and then gradually lowering the corresponding Rayleigh number. We show that decreasing the Ekman number leads to an extension of the subcritical range of Ra/Ra$_c$, down to an optimal value of Ek$=10^{-5}$. For magnetic Prandtl numbers of order unity, subcriticality is then hampered by the emergence of a large-scale mode at lower Ekman numbers when the dynamo driven by the smaller scale convection generates relatively stronger large-scale magnetic field. The inability of the large-scale mode to sustain dynamo action leads to an intermittent behaviour that appears to inhibit subcriticality. The subcritical solutions are also sensitive to the value of the magnetic Reynolds number (or equivalently, the magnetic Prandtl number, Pm), as values of the magnetic Reynolds number greater than 70 are required to produce dynamo action, but large values lead to fluctuations that are able to push the system too far from the subcritical branch and towards the trivial conducting state., Comment: Accepted for publication in Phys. Rev. Fluids

Motivated by the dynamics in the deep interiors of many stars, we study the interaction between overshooting convection and the large-scale poloidal fields residing in radiative zones. We have run a suite of 3D Boussinesq numerical calculations in a spherical shell that consists of a convection zone with an underlying stable region that initially compactly contains a dipole field. By varying the strength of the convective driving, we find that, in the less turbulent regime, convection acts as turbulent diffusion that removes the field faster than solely molecular diffusion would do. However, in the more turbulent regime, turbulent pumping becomes more efficient and partially counteracts turbulent diffusion, leading to a local accumulation of the field below the overshoot region. These simulations suggest that dipole fields might be confined in underlying stable regions by highly turbulent convective motions at stellar parameters. The confinement is of large-scale field in an average sense and we show that it is reasonably modeled by mean-field ideas. Our findings are particularly interesting for certain models of the Sun, which require a large-scale, poloidal magnetic field to be confined in the solar radiative zone in order to explain simultaneously the uniform rotation of the latter and the thinness of the solar tachocline., Comment: Accepted to MNRAS, 14 figures

The possibility of subcritical behaviour in the geostrophic turbulence regime of rapidly rotating thermally driven convection is explored. In this regime a non-local inverse energy transfer may compete with the more traditional and local direct cascade. We show that, even for control parameters for which no inverse cascade has previously been observed, a subcritical transition towards a large-scale vortex state can occur when the system is initialized with a vortex dipole of finite amplitude. This new example of bistability in a turbulent flow, which may not be specific to rotating convection, opens up new avenues for studying energy transfer in strongly anisotropic three-dimensional flows., Comment: 12 pages, 6 figures

Convection is a fundamental physical process in the fluid cores of planets because it is the primary transport mechanism for heat and chemical species and the primary energy source for planetary magnetic fields. Key properties of convection, such as the characteristic flow velocity and lengthscale, are poorly quantified in planetary cores due to their strong dependence on planetary rotation, buoyancy driving and magnetic fields, which are all difficult to model under realistic conditions. In the absence of strong magnetic fields, the core convective flows are expected to be in a regime of rapidly-rotating turbulence, which remains largely unexplored to date. Here we use a combination of numerical models designed to explore this low-viscosity regime to show that the convective lengthscale becomes independent of the viscosity and is entirely determined by the flow velocity and planetary rotation. For the Earth's core, we find that the characteristic con-vective lengthscale is approximately 30km and below this scale, motions are very weak. The 30-km cutoff scale rules out small-scale dynamo action and supports large-eddy simulations of core dynamics. Furthermore, it implies that our understanding of magnetic reversals from numerical geodynamo models does not relate to the Earth, because they require too intense flows. Our results also indicate that the liquid core of the Moon might still be in an active convective state despite the absence of a present-day dynamo.

One of the most prominent dynamical features of turbulent, rapidly-rotating convection is the formation of large-scale coherent structures, driven by Reynolds stresses resulting from the small-scale convective flows. In spherical geometry, such structures consist of intense zonal flows that are invariant along the rotation axis. In planar geometry, long-lived, depth-invariant structures also form at large scales, but, in the absence of horizontal anisotropy, they consist of vortices that grow to the domain size. In this work, through the introduction of horizontal anisotropy into a numerical model of planar rotating convection by the adoption of unequal horizontal box sizes (i.e. $L_x \le L_y$, where the $xy$-plane is horizontal), we investigate whether unidirectional flows and large-scale vortices can coexist. We find that only a small degree of anisotropy is required to bring about a transition from dynamics dominated by persistent large-scale vortices to dynamics dominated by persistent unidirectional flows parallel to the shortest horizontal direction. When the anisotropy is sufficiently large, the unidirectional flow consists of multiple jets, generated on a timescale smaller than a global viscous timescale, thus signifying that the upscale energy transfer does not spontaneously feed the largest available mode in the system. That said, the multiple jets merge on much longer timescales. Large-scale vortices of size comparable with $L_x$ systematically form in the flanks of the jets and can be persistent or intermittent. This indicates that large-scale vortices, either coexisting with jets or not, are a robust dynamical feature of planar rotating convection., Comment: submitted to Physical Review Fluids

We study rapidly-rotating Boussinesq convection driven by internal heating in a full sphere. We use a numerical model based on the quasi-geostrophic approximation for the velocity field, whereas the temperature field is three-dimensional. This approximation allows us to perform simulations for Ekman numbers down to 1e-8, Prandtl numbers relevant for liquid metals (~0.1) and Reynolds numbers up to 3e4. Persistent zonal flows composed of multiple jets form as a result of the mixing of potential vorticity. For the largest Rayleigh numbers computed, the zonal velocity is larger than the convective velocity despite the presence of boundary friction. The convective structures and the zonal jets widen when the thermal forcing increases. Prograde and retrograde zonal jets are dynamically different: in the prograde jets (which correspond to weak potential vorticity gradients) the convection transports heat efficiently and the mean temperature tends to be homogenised; by contrast, in the cores of the retrograde jets (which correspond to steep gradients of potential vorticity) the dynamics is dominated by the propagation of Rossby waves, resulting in the formation of steep mean temperature gradients and the dominance of conduction in the heat transfer process. Consequently, in quasi-geostrophic systems, the width of the retrograde zonal jets controls the efficiency of the heat transfer., Comment: 25 pages, 13 figues, published in Geophys. J. Int

Several recent studies have demonstrated how large-scale vortices may arise spontaneously in rotating planar convection. Here we examine the dynamo properties of such flows in rotating Boussinesq convection. For moderate values of the magnetic Reynolds number ($100 \lesssim Rm \lesssim 550$, with $Rm$ based on the box depth and the convective velocity), a large-scale (i.e. system-size) magnetic field is generated. The amplitude of the magnetic energy oscillates in time, nearly out of phase with the oscillating amplitude of the large-scale vortex. The large-scale vortex is disrupted once the magnetic field reaches a critical strength, showing that these oscillations are of magnetic origin. The dynamo mechanism relies on those components of the flow that have length scales lying between that of the large-scale vortex and the typical convective cell size; smaller-scale flows are not required. The large-scale vortex plays a crucial role in the magnetic induction despite being essentially two-dimensional; we thus refer to this dynamo as a large-scale-vortex dynamo. For larger magnetic Reynolds numbers, the dynamo is small scale, with a magnetic energy spectrum that peaks at the scale of the convective cells. In this case, the small-scale magnetic field continuously suppresses the large-scale vortex by disrupting the correlations between the convective velocities that allow it to form. The suppression of the large-scale vortex at high $Rm$ therefore probably limits the relevance of the large-scale-vortex dynamo to astrophysical objects with moderate values of $Rm$, such as planets. In this context, the ability of the large-scale-vortex dynamo to operate at low magnetic Prandtl numbers is of great interest., Comment: 28 pages, 17 figures, accepted for publication in J. Fluid Mech

We study nonlinear convection in a rapidly rotating sphere with internal heating for values of the Prandtl number relevant for liquid metals ($Pr\in[10^{-2},10^{-1}]$). We use a numerical model based on the quasi-geostrophic approximation, in which variations of the axial vorticity along the rotation axis are neglected, whereas the temperature field is fully three-dimensional. We identify two separate branches of convection close to onset: (i) a well-known weak branch for Ekman numbers greater than $10^{-6}$, which is continuous at the onset (supercritical bifurcation) and consists of thermal Rossby waves, and (ii) a novel strong branch at lower Ekman numbers, which is discontinuous at the onset. The strong branch becomes subcritical for Ekman numbers of the order of $10^{-8}$. On the strong branch, the Reynolds number of the flow is greater than $10^3$, and a strong zonal flow with multiple jets develops, even close to the nonlinear onset of convection. We find that the subcriticality is amplified by decreasing the Prandtl number. The two branches can co-exist for intermediate Ekman numbers, leading to hysteresis ($Ek=10^{-6}$, $Pr=10^{-2}$). Nonlinear oscillations are observed near the onset of convection for $Ek=10^{-7}$ and $Pr=10^{-1}$., Comment: 30 pages, 16 figures, published in JFM

We propose a new self-consistent dynamo mechanism for the generation of large-scale magnetic fields in natural objects. Recent computational studies have described the formation of large-scale vortices (LSVs) in rotating turbulent convection. Here we demonstrate that for magnetic Reynolds numbers below the threshold for small-scale dynamo action, such turbulent flows can sustain large-scale magnetic fields --- i.e. fields with a significant component on the scale of the system., Comment: 8 pages, 6 figures, accepted for publication in Physical Review E

Using numerical simulations of rapidly rotating Boussinesq convection in a Cartesian box, we study the formation of long-lived, large-scale, depth-invariant coherent structures. These structures, which consist of concentrated cyclones, grow to the horizontal size of the box, with velocities significantly larger than the convective motions. We vary the rotation rate, the thermal driving and the aspect ratio in order to determine the domain of existence of these large-scale vortices (LSV). We find that two conditions are required for their formation. First, the Rayleigh number, a meaure of the thermal driving, must be several times its value at the linear onset of convection; this corresponds to Reynolds numbers, based on the convective velocity and the box depth, $\gtrsim 100$. Second, the rotational constraint on the convective structures must be strong. This requires that the local Rossby number, based on the convective velocity and the horizontal convective scale, $\lesssim 0.15$. Simulations in which certain wavenumbers are artificially suppressed in spectral space suggest that the LSV are produced by the interactions of small-scale, depth-dependent convective motions. The presence of LSV significantly reduces the efficiency of the convective heat transport., Comment: 31 pages, 17 figures, published in J. Fluid Mech

We present a numerical study of dynamo action in a conducting fluid encased in a metallic spherical shell. Motions in the fluid are driven by differential rotation of the outer metallic shell, which we refer to as "the wall". The two hemispheres of the wall are held in counter-rotation, producing a steady, axisymmetric interior flow consisting of differential rotation and a two-cell meridional circulation with radial inflow in the equatorial plane. From previous studies, this type of flow is known to maintain a stationary equatorial dipole by dynamo action if the magnetic Reynolds number is larger than about 300 and if the outer boundary is electrically insulating. We vary independently the thickness, electrical conductivity, and magnetic permeability of the wall to determine their effect on the dynamo action. The main results are: (a) Increasing the conductivity of the wall hinders the dynamo by allowing eddy currents within the wall, which are induced by the relative motion of the equatorial dipole field and the wall. This processes can be viewed as a skin effect or, equivalently, as the tearing apart of the dipole by the differential rotation of the wall, to which the field lines are anchored by high conductivity. (b) Increasing the magnetic permeability of the wall favors dynamo action by constraining the magnetic field lines in the fluid to be normal to the wall, thereby decoupling the fluid from any induction in the wall. (c) Decreasing the wall thickness limits the amplitude of the eddy currents, and is therefore favorable for dynamo action, provided that the wall is thinner than the skin depth. We explicitly demonstrate these effects of the wall properties on the dynamo field by deriving an effective boundary condition in the limit of vanishing wall thickness., Comment: accepted for publication in Physical Review E

We have performed numerical simulations of boundary-driven dynamos using a three-dimensional non-linear magnetohydrodynamical model in a spherical shell geometry. A conducting fluid of magnetic Prandtl number Pm=0.01 is driven into motion by the counter-rotation of the two hemispheric walls. The resulting flow is of von K\'arm\'an type, consisting of a layer of zonal velocity close to the outer wall and a secondary meridional circulation. Above a certain forcing threshold, the mean flow is unstable to non-axisymmetric motions within an equatorial belt. For fixed forcing above this threshold, we have studied the dynamo properties of this flow. The presence of a conducting outer wall is essential to the existence of a dynamo at these parameters. We have therefore studied the effect of changing the material parameters of the wall (magnetic permeability, electrical conductivity, and thickness) on the dynamo. In common with previous studies, we find that dynamos are obtained only when either the conductivity or the permeability is sufficiently large. However, we find that the effect of these two parameters on the dynamo process are different and can even compete to the detriment of the dynamo. Our self-consistent approach allow us to analyze in detail the dynamo feedback loop. The dynamos we obtain are typically dominated by an axisymmetric toroidal magnetic field and an axial dipole component. We show that the ability of the outer shear layer to produce a strong toroidal field depends critically on the presence of a conducting outer wall, which shields the fluid from the vacuum outside. The generation of the axisymmetric poloidal field, on the other hand, occurs in the equatorial belt and does not depend on the wall properties., Comment: accepted for publication in Physical Review E

We present a dynamo mechanism arising from the presence of barotropically unstable zonal jet currents in a rotating spherical shell. The shear instability of the zonal flow develops in the form of a global Rossby mode, whose azimuthal wavenumber depends on the width of the zonal jets. We obtain self-sustained magnetic fields at magnetic Reynolds numbers greater than 1000. We show that the propagation of the Rossby waves is crucial for dynamo action. The amplitude of the axisymmetric poloidal magnetic field depends on the wavenumber of the Rossby mode, and hence on the width of the zonal jets. We discuss the plausibility of this dynamo mechanism for generating the magnetic field of the giant planets. Our results suggest a possible link between the topology of the magnetic field and the profile of the zonal winds observed at the surface of the giant planets. For narrow Jupiter-like jets, the poloidal magnetic field is dominated by an axial dipole whereas for wide Neptune-like jets, the axisymmetric poloidal field is weak., Comment: published in Icarus

We numerically investigate the efficiency of a spherical Couette flow at generating a self-sustained magnetic field. No dynamo action occurs for axisymmetric flow while we always found a dynamo when non-axisymmetric hydrodynamical instabilities are excited. Without rotation of the outer sphere, typical critical magnetic Reynolds numbers $Rm_c$ are of the order of a few thousands. They increase as the mechanical forcing imposed by the inner core on the flow increases (Reynolds number $Re$). Namely, no dynamo is found if the magnetic Prandtl number $Pm=Rm/Re$ is less than a critical value $Pm_c\sim 1$. Oscillating quadrupolar dynamos are present in the vicinity of the dynamo onset. Saturated magnetic fields obtained in supercritical regimes (either $Re>2 Re_c$ or $Pm>2Pm_c$) correspond to the equipartition between magnetic and kinetic energies. A global rotation of the system (Ekman numbers $E=10^{-3}, 10^{-4}$) yields to a slight decrease (factor 2) of the critical magnetic Prandtl number, but we find a peculiar regime where dynamo action may be obtained for relatively low magnetic Reynolds numbers ($Rm_c\sim 300$). In this dynamical regime (Rossby number $Ro\sim -1$, spheres in opposite direction) at a moderate Ekman number ($E=10^{-3}$), a enhanced shear layer around the inner core might explain the decrease of the dynamo threshold. For lower $E$ ($E=10^{-4}$) this internal shear layer becomes unstable, leading to small scales fluctuations, and the favorable dynamo regime is lost. We also model the effect of ferromagnetic boundary conditions. Their presence have only a small impact on the dynamo onset but clearly enhance the saturated magnetic field in the ferromagnetic parts. Implications for experimental studies are discussed.