Search

In order to characterize the global circulation of the subsurface ocean of Jovian and Saturnian moons, we analyze the properties of 21 three-dimensional simulations of Boussinesq thermal convection in a rapidly rotating spherical shell. Flow is driven by an adverse temperature contrast imposed across the domain, and is subjected to no-slip boundary conditions. We cover a region of parameter space previously unexplored by global simulations, both in terms of rapid rotation and vigor of convective forcing, closer to, yet still admittedly far from, the conditions appropriate for the subsurface ocean of Ganymede, Europa, Enceladus, and Titan. Our most extreme simulations exhibit a dynamic global circulation that combines powerful east-west zonal jets, planetary waves, and vortices. A spectral analysis of the kinetic energy distribution performed in cylindrical geometry reveals a high degree of anisotropy of the simulated flows. Specifically, the axisymmetric zonal energy spectra follow a steep $-5$ slope in wavenumber space, with the energy amplitude exclusively controlled by the rotation rate. In contrast, the non-axisymmetric residual spectra display a gentle $-5/3$ slope, with the energy amplitude controlled by the thermal buoyancy input power. This spectral behavior conforms with the theory of zonostrophic turbulence and allows us to propose tentative extrapolations of these findings to the more extreme conditions of icy satellites. By assuming that kinetic energy dissipates via Ekman friction we predict an upper bound for the zonal velocity ranging from a few centimeters per second for Enceladus to about one meter per second for Ganymede, with residual velocities smaller than the zonal velocity by an order of magnitude on each moon. These predictions yield typical jets size approaching the ocean depth of Titan, Ganymede and Europa and $10$ to $40\%$ of the ocean depth on Enceladus., Comment: 24 pages, 12 figures, 3 tables, accepted for publication in Icarus

We use 123 three dimensional direct numerical simulations to study fingering convection in non-rotating spherical shells. We investigate the scaling behaviour of the flow lengthscale, the non-dimensional heat and compositional fluxes $Nu$ and $Sh$ and the mean convective velocity over the fingering convection instability domain defined by $1 \leq R_\rho < Le$, $R_\rho$ being the ratio of density perturbations of thermal and compositional origins and $Le$ the Lewis number. We show that the chemical boundary layers are marginally unstable and adhere to the laminar Prandtl-Blasius model, hence explaining the asymmetry between the inner and outer spherical shell boundary layers. We develop scaling laws for two asymptotic regimes close to the two edges of the instability domain, namely $R_\rho \lesssim Le$ and $R_\rho \gtrsim 1$. For the former, we develop novel power laws of a small parameter $\epsilon$ measuring the distance to onset, which differ from theoretical laws published to date in Cartesian geometry. For the latter, we find that the Sherwood number $Sh$ gradually approaches a scaling $Sh\sim Ra_\xi^{1/3}$ when $Ra_\xi \gg 1$; and that the P\'eclet number accordingly follows $Pe \sim Ra_\xi^{2/3} |Ra_T|^{-1/4}$, $Ra_\xi$ being the chemical Rayleigh number. When the Reynolds number exceeds a few tens, we report on a secondary instability which takes the form of large-scale toroidal jets which span the entire spherical domain. Jets distort the fingers resulting in Reynolds stress correlations, which in turn feed the jet growth until saturation. This nonlinear phenomenon can yield relaxation oscillation cycles., Comment: 45 pages, 21 figures, 3 tables, accepted for publication in JFM

Convection occurs ubiquitously on and in rotating geophysical and astrophysical bodies. Prior spherical shell studies have shown that the convection dynamics in polar regions can differ significantly from the lower latitude, equatorial dynamics. Yet most spherical shell convective scaling laws use globally-averaged quantities that erase latitudinal differences in the physics. Here we quantify those latitudinal differences by analyzing spherical shell simulations in terms of their regionalized convective heat transfer properties. This is done by measuring local Nusselt numbers in two specific, latitudinally separate, portions of the shell, the polar and the equatorial regions, $Nu_p$ and $Nu_e$, respectively. In rotating spherical shells, convection first sets in outside the tangent cylinder such that equatorial heat transfer dominates at small and moderate supercriticalities. We show that the buoyancy forcing, parameterized by the Rayleigh number $Ra$, must exceed the critical equatorial forcing by a factor of $\approx 20$ to trigger polar convection within the tangent cylinder. Once triggered, $Nu_p$ increases with $Ra$ much faster than does $Nu_e$. The equatorial and polar heat fluxes then tend to become comparable at sufficiently high $Ra$. Comparisons between the polar convection data and Cartesian numerical simulations reveal quantitative agreement between the two geometries in terms of heat transfer and averaged bulk temperature gradient. This agreement indicates that spherical shell rotating convection dynamics are accessible both through spherical simulations and via reduced investigatory pathways, be they theoretical, numerical or experimental., Comment: 11 pages, 6 figures, accepted for publication in JFM

Observations of surface magnetic fields of cool stars reveal a large diversity of configurations. Although there is now a consensus that these fields are generated through dynamo processes occurring within the convective zone, the physical mechanism driving such a variety of field topologies is still debated. This paper discusses the possible origins of dipole and multipole-dominated morphologies using three-dimensional numerical simulations of stratified systems where the magnetic feedback on the fluid motion is significant. Our main result is that dipolar solutions are found at Rossby numbers up to 0.4 in strongly stratified simulations, where previous works suggested that only multipolar fields should exist. We argue that these simulations are reminiscent of the outlier stars observed at Rossby numbers larger than 0.1, whose large-scale magnetic field is dominated by their axisymmetric poloidal component. As suggested in previous Boussinesq calculations, the relative importance of inertial over Lorentz forces is again controlling the dipolar to multipolar transition. Alternatively, we find that the ratio of kinetic to magnetic energies can equally well capture the transition in the field morphology. We test the ability of this new proxy to predict the magnetic morphology of a few M-dwarf stars whose internal structure matches that of our simulations and for which homogeneous magnetic field characterization is available. Finally, the magnitude of the differential rotation obtained in our simulations is compared to actual measurements reported in the literature for M-dwarfs. In our simulations, we find a clear relationship between anti-solar differential rotation and the emergence of dipolar fields., Comment: 15 pages, 15 figures, accepted for publication in MNRAS

We analyze the behaviour of an ensemble of time integrators applied to the semi-discrete problem resulting from the spectral discretization of the equations describing Boussinesq convection in a cylindrical annulus. The equations are cast in a vorticity-streamfunction formulation that yields a differential algebraic equation (DAE). The ensemble comprises 28 members: 4 implicit-explicit multistep schemes, 22 implicit-explicit Runge-Kutta (IMEX-RK) schemes, and 2 fully explicit schemes used for reference. The schemes are assessed for 11 different physical setups that cover laminar and turbulent regimes. Multistep and order 2 IMEX-RK methods exhibit their expected order of convergence under all circumstances. IMEX-RK methods of higher-order show occasional order reduction that impacts both algebraic and differential field variables. We ascribe the order reduction to the stiffness of the problem and, to a larger extent, the presence of the DAE. Using the popular Crank-Nicolson Adams-Bashforth of order 2 (CNAB2) integrator as reference, performance is defined by the ratio of maximum admissible time step to the cost of performing one iteration; the maximum admissible time step is determined by inspection of the time series of viscous dissipation within the system, which guarantees a physically acceptable solution. Relative performance is bounded between 0.5 and 1.5 across all studied configurations. Considering accuracy jointly with performance, we find that 6 schemes consistently outperform CNAB2, meaning that in addition to allowing for a more efficient calculation, the accuracy that they achieve at their operational limit of stability yields a lower error. In our most turbulent setup, where the behaviour of the methods is almost entirely dictated by their explicit component, 13 IMEX-RK integrators outperform CNAB2 in terms of accuracy and efficiency., Comment: 39 pages, 15 figures, 5 tables, accepted for publication in JCP

The ongoing NASA's Juno mission puts new constraints on the internal dynamics of Jupiter. Data gathered by its onboard magnetometer reveal a dipole-dominated surface magnetic field accompanied by strong localised magnetic flux patches. The gravity measurements indicate that the fierce surface zonal jets extend several thousands of kilometers below the cloud level before rapidly decaying below $0.94-0.96\,R_J$, $R_J$ being the mean Jovian radius at the one bar level. Several internal models suggest an intricate internal structure with a thin intermediate region in which helium would segregate from hydrogen, forming a compositionally-stratified layer. Here, we develop the first global Jovian dynamo which incorporates an intermediate stably-stratified layer between $0.82\,R_J$ and $0.86\,R_J$. Analysing the energy balance reveals that the magnetic energy is almost one order of magnitude larger than kinetic energy in the metallic region, while most of the kinetic energy is pumped into zonal motions in the molecular envelope. Those result from the different underlying force hierarchy with a triple balance between Lorentz, Archimedean and ageostrophic Coriolis forces in the metallic core and inertia, buoyancy and ageostrophic Coriolis forces controlling the external layers. The simulation presented here is the first to demonstrate that multiple zonal jets and dipole-dominated dynamo action can be consolidated in a global simulation. The inclusion of an stable layer is a necessary ingredient that allows zonal jets to develop in the outer envelope without contributing to the dynamo action in the deeper metallic region. Stable stratification however also smooths out the small-scale features of the magnetic field by skin effect. These constraints suggest that possible stable layers in Jupiter should be located much closer to the surface ($0.9-0.95\,R_J$)., Comment: 23 pages, 16 figures, 2 tables, accepted for publication in Icarus

Dynamo action in the Earth's outer core is expected to be controlled by a balance between pressure, Coriolis, buoyancy and Lorentz forces, with marginal contributions from inertia and viscous forces. Current numerical simulations of the geodynamo, however, operate at much larger inertia and viscosity because of computational limitations. This casts some doubt on the physical relevance of these models. Our work aims at finding dynamo models in a moderate computational regime which reproduce the leading-order force balance of the Earth. By performing a systematic parameter space survey with Ekman numbers in the range $10^{-6} \leq E \leq 10^{-4}$, we study the variations of the force balance when changing the forcing (Rayleigh number, $Ra$) and the ratio between viscous and magnetic diffusivities (magnetic Prandtl number, $Pm$). For dipole-dominated dynamos, we observe that the force balance is structurally robust throughout the investigated parameter space, exhibiting a quasi-geostrophic (QG) balance (balance between Coriolis and pressure forces) at zeroth order, followed by a first-order MAC balance between the ageostrophic Coriolis, buoyancy and Lorentz forces. At second order this balance is disturbed by contributions from inertia and viscous forces. Dynamos with a different sequence of the forces, where inertia and/or viscosity replace the Lorentz force in the first-order force balance, can only be found close to the onset of dynamo action and in the multipolar regime. Our study illustrates that most classical numerical dynamos are controlled by a QG-MAC balance, while cases where viscosity and inertia play a dominant role are the exception rather than the norm.

We present a new pseudo-spectral open-source code nicknamed pizza. It is dedicated to the study of rapidly-rotating Boussinesq convection under the 2-D spherical quasi-geostrophic approximation, a physical hypothesis that is appropriate to model the turbulent convection that develops in planetary interiors. The code uses a Fourier decomposition in the azimuthal direction and supports both a Chebyshev collocation method and a sparse Chebyshev integration formulation in the cylindrically-radial direction. It supports several temporal discretisation schemes encompassing multi-step time steppers as well as diagonally-implicit Runge-Kutta schemes. The code has been tested and validated by comparing weakly-nonlinear convection with the eigenmodes from a linear solver. The comparison of the two radial discretisation schemes has revealed the superiority of the Chebyshev integration method over the classical collocation approach both in terms of memory requirements and operation counts. The good parallelisation efficiency enables the computation of large problem sizes with $\mathcal{O}(10^4\times 10^4)$ grid points using several thousands of ranks. This allows the computation of numerical models in the turbulent regime of quasi-geostrophic convection characterised by large Reynolds $Re$ and yet small Rossby numbers $Ro$. A preliminary result obtained for a strongly supercritical numerical model with a small Ekman number of $10^{-9}$ and a Prandtl number of unity yields $Re\simeq 10^5$ and $Ro \simeq 10^{-4}$. pizza is hence an efficient tool to study spherical quasi-geostrophic convection in a parameter regime inaccessible to current global 3-D spherical shell models., Comment: 20 pages, 10 figures, 4 tables, accepted for publication in GJI

Rayleigh-B\'enard convection in rotating spherical shells can be considered as a simplified analogue of many astrophysical and geophysical fluid flows. Here, we use three-dimensional direct numerical simulations to study this physical process. We construct a dataset of more than 200 numerical models that cover a broad parameter range with Ekman numbers spanning $3\times 10^{-7} \leq E \leq 10^{-1}$, Rayleigh numbers within the range $10^3 < Ra < 2\times 10^{10}$ and a Prandtl number unity. We investigate the scaling behaviours of both local (length scales, boundary layers) and global (Nusselt and Reynolds numbers) properties across various physical regimes from onset of rotating convection to weakly-rotating convection. Close to critical, the convective flow is dominated by a triple force balance between viscosity, Coriolis force and buoyancy. For larger supercriticalities, a subset of our numerical data approaches the asymptotic diffusivity-free scaling of rotating convection $Nu\sim Ra^{3/2}E^{2}$ in a narrow fraction of the parameter space delimited by $6\,Ra_c \leq Ra \leq 0.4\,E^{-8/5}$. Using a decomposition of the viscous dissipation rate into bulk and boundary layer contributions, we establish a theoretical scaling of the flow velocity that accurately describes the numerical data. In rapidly-rotating turbulent convection, the fluid bulk is controlled by a triple force balance between Coriolis, inertia and buoyancy, while the remaining fraction of the dissipation can be attributed to the viscous friction in the Ekman layers. Beyond $Ra \simeq E^{-8/5}$, the rotational constraint on the convective flow is gradually lost and the flow properties vary to match the regime changes between rotation-dominated and non-rotating convection. The quantity $Ra E^{12/7}$ provides an accurate transition parameter to separate rotating and non-rotating convection., Comment: 42 pages, 20 figures, 3 tables, accepted for publication in JFM

We simulate numerically Boussinesq convection in non-rotating spherical shells for a fluid with a unity Prandtl number and Rayleigh numbers up to $10^9$. In this geometry, curvature and radial variations of the gravitationnal acceleration yield asymmetric boundary layers. A systematic parameter study for various radius ratios (from $\eta=r_i/r_o=0.2$ to $\eta=0.95$) and gravity profiles allows us to explore the dependence of the asymmetry on these parameters. We find that the average plume spacing is comparable between the spherical inner and outer bounding surfaces. An estimate of the average plume separation allows us to accurately predict the boundary layer asymmetry for the various spherical shell configurations explored here. The mean temperature and horizontal velocity profiles are in good agreement with classical Prandtl-Blasius laminar boundary layer profiles, provided the boundary layers are analysed in a dynamical frame, that fluctuates with the local and instantaneous boundary layer thicknesses. The scaling properties of the Nusselt and Reynolds numbers are investigated by separating the bulk and boundary layer contributions to the thermal and viscous dissipation rates using numerical models with $\eta=0.6$ and a gravity proportional to $1/r^2$. We show that our spherical models are consistent with the predictions of Grossmann \& Lohse's (2000) theory and that $Nu(Ra)$ and $Re(Ra)$ scalings are in good agreement with plane layer results., Comment: 43 pages, 25 figures, 2 tables, accepted for publication in JFM