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