Search

We present $\texttt{cunusht}$, a general-purpose Python package that wraps a highly efficient CUDA implementation of the nonuniform spin-$0$ spherical harmonic transform. The method is applicable to arbitrary pixelization schemes, including schemes constructed from equally-spaced iso-latitude rings as well as completely nonuniform ones. The algorithm has an asymptotic scaling of $\mathrm{O}{(\ell_{\rm max}^3)}$ for maximum multipole $\ell_{\rm max}$ and achieves machine precision accuracy. While $\texttt{cunusht}$ is developed for applications in cosmology in mind, it is applicable to various other interpolation problems on the sphere. We outperform the fastest available CPU algorithm by a factor of up to 5 for problems with a nonuniform pixelization and $\ell_{\rm max}>4\cdot10^3$ when comparing a single modern GPU to a modern 32-core CPU. This performance is achieved by utilizing the double Fourier sphere method in combination with the nonuniform fast Fourier transform and by avoiding transfers between the host and device. For scenarios without GPU availability, $\texttt{cunusht}$ wraps existing CPU libraries. $\texttt{cunusht}$ is publicly available and includes tests, documentation, and demonstrations.

Planetary cores are the seat of rich and complex fluid dynamics, in which the effects of rotation and magnetic field combine. The equilibria governing the strength of the magnetic field produced by the dynamo effect, the organisation and amplitude of the flow, and those of the density field, remain debated despite remarkable progress made in their numerical simulation. This paper proposes a new approach based on the explicit consideration of the variation of time scales {\tau} with spatial scales {\ell} for the different physical phenomena involved. The {\tau}-{\ell} diagrams thus constructed constitute a very complete graphic summary of the dynamics of the object under study. We highlight the role of the available convective power in controlling this dynamics, together with the relevant force balance, for which we derive a very telling {\tau}-{\ell} translation. Several scenarios are constructed and discussed for the Earth's core, shedding new light on the width of convective columns and on the force equilibria to be considered. A QG-MAC scenario adapted from Aubert [2019] gives a good account of the observations. A diversion to Venus reveals the subtlety and relativity of the notion of 'fast rotator'. We discuss scaling laws and their validity domain, and illustrate 'path strategies'. A complete toolbox is provided, allowing everyone to construct a {\tau}-{\ell} diagram of a numerical simulation, a laboratory experiment, a theory, or a natural object. Supplementary material for this article is supplied as a separate archive Nataf_Schaeffer_SupMat.zip, the related data is displayed in document Nataf_Schaeffer_SupMat.pdf., Comment: Supplementary Material available with python programs at: https://gricad-gitlab.univ-grenoble-alpes.fr/natafh/shell_tau-ell_programs

The generation of mean flows is a long-standing issue in rotating fluids. Motivated by planetary objects, we consider here a rapidly rotating fluid-filled spheroid, which is subject to weak perturbations of either the boundary (e.g. tides) or the rotation vector (e.g. in direction by precession, or in magnitude by longitudinal librations). Using boundary-layer theory, we determine the mean zonal flows generated by nonlinear interactions within the viscous Ekman layer. These flows are of interest because they survive in the relevant planetary regime of both vanishing forcings and viscous effects. We extend the theory to take into account (i) the combination of spatial and temporal perturbations, providing new mechanically driven zonal flows (e.g. driven by latitudinal librations), and (ii) the spheroidal geometry relevant for planetary bodies. Wherever possible, our analytical predictions are validated with direct numerical simulations. The theoretical solutions are in good quantitative agreement with the simulations, with expected discrepancies (zonal jets) in the presence of inertial waves generated at the critical latitudes (as for precession). Moreover, we find that the mean zonal flows can be strongly affected in spheroids. Guided by planetary applications, we also revisit the scaling laws for the geostrophic shear layers at the critical latitudes, and the influence of a solid inner core., Comment: 45 pages, 20 figures. Published online 14 April 2021 in J. Fluid Mech. Correction of a minor typo in the caption of figure 17 of the published version

Convection in a spherical shell is widely used to model fluid layers of planets and stars. The choice of thermal boundary conditions in such models is not always straightforward. To understand the implications of this choice, we report on the effects of the thermal boundary condition on thermal convection, in terms of instability onset, fully developed transport properties and flow structure. We use the Boussinesq approximation, and impose a Robin boundary condition at the top. This enforces the temperature anomaly and its radial derivative to be linearly coupled with a proportionality factor $\beta$. Using the height H of the fluid layer, we introduce the non-dimensional Biot number Bi* = $\beta$H. Varying Bi* allows us to transition from fixed temperature for Bi* = +$\infty$, to fixed thermal flux for Bi* = 0. The bottom boundary of the shell is kept isothermal. We find that the onset of convection is only affected by Bi* in the non-rotating case. Far from onset, considering an effective Rayleigh number and a generalized Nusselt number, we show that the Nusselt and P{\'e}clet numbers follow standard universal scaling laws, independent of Bi* in all cases considered. However, the large-scale flow structure keeps the signature of the boundary condition with more vigorous large scales for smaller Bi* , even though the global heat transfer and kinetic energy are the same. Finally, for all practical purposes, the Robin condition can be safely replaced by a fixed flux when Bi* < 0.03 and by a fixed temperature for Bi* > 30.

In planetary fluid cores, the density depends on temperature and chemical composition, which diffuse at very different rates. This leads to various instabilities, bearing the name of double-diffusive convection. We investigate rotating double-diffusive convection (RDDC) in fluid spheres. We use the Boussinesq approximation with homogeneous internal thermal and compositional source terms. We focus on the finger regime, in which the thermal gradient is stabilising whereas the compositional one is destabilising. First, we perform a global linear stability analysis in spheres. The critical Rayleigh numbers drastically drop for stably stratified fluids, yielding large-scale convective motions where local analyses predict stability. We evidence the inviscid nature of this large-scale double-diffusive instability, enabling the determination of the marginal stability curve at realistic planetary regimes. In particular, we show that in stably stratified spheres, the Rayleigh numbers $Ra$ at the onset evolve like $Ra \sim Ek^{-1}$, where $Ek$ is the Ekman number. This differs from rotating convection in unstably stratified spheres, for which $Ra \sim Ek^{-4/3}$. The domain of existence of inviscid convection thus increases as $Ek^{-1/3}$. Second, we perform nonlinear simulations. We find a transition between two regimes of RDDC, controlled by the strength of the stratification. Furthermore, far from the RDDC onset, we find a dominating equatorially anti-symmetric, large-scale zonal flow slightly above the associated linear onset. Unexpectedly, a purely linear mechanism can explain this phenomenon, even far from the instability onset, yielding a symmetry breaking of the nonlinear flow at saturation. For even stronger stable stratification, the flow becomes mainly equatorially-symmetric and intense zonal jets develop. Finally, we apply our results to the early Earth core. Double diffusion can reduce the critical Rayleigh number by four decades for realistic core conditions. We suggest that the early Earth core was prone to turbulent RDDC, with large-scale zonal flows.

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.

Precession of planets or moons affects internal liquid layers by driving flows, instabilities and possibly dynamos.The energy dissipated by these phenomena can influence orbital parameters such as the planet's spin rate.However, there is no systematic study of these flows in the spherical shell geometry relevant for planets, and the lack of scaling law prevents convincing extrapolation to celestial bodies.We have run more than 900 simulations of fluid spherical shells affected by precession, to systematically study basic flows, instabilities, turbulence, and magnetic field generation.We observe no significant effects of the inner core on the onset of the instabilities.We obtain an analytical estimate of the viscous dissipation, mostly due to boundary layer friction in our simulations.We propose theoretical onsets for hydrodynamic instabilities, and document the intensity of turbulent fluctuations.We extend previous precession dynamo studies towards lower viscosities, at the limits of today's computers.In the low viscosity regime, precession dynamos rely on the presence of large-scale vortices, and the surface magnetic fields are dominated by small scales.Interestingly, intermittent and self-killing dynamos are observed.Our results suggest that large-scale planetary magnetic fields are unlikely to be produced by a precession-driven dynamo in a spherical core.But this question remains open as planetary cores are not exactly spherical, and thus the coupling between the fluid and the boundary does not vanish in the relevant limit of small viscosity.Moreover, the fully turbulent dissipation regime has not yet been reached in simulations.Our results suggest that the melted lunar core has been in a turbulent state throughout its history.Furthermore, in the view of recent experimental results, we propose updated formulas predicting the fluid mean rotation vector and the associated dissipation in both the laminar and the turbulent regimes.

High-performance computing trends towards many-core systems are expected to continue over the next decade. As a result, parallel-in-time methods, mathematical formulations which exploit additional degrees of parallelism in the time dimension, have gained increasing interest in recent years. In this work we study a massively parallel rational approximation of exponential integrators (REXI). This method replaces a time integration of stiff linear oscillatory and diffusive systems by the sum of the solutions of many decoupled systems, which can be solved in parallel. Previous numerical studies showed that this reformulation allows taking arbitrarily long time steps for the linear oscillatory parts. The present work studies the non-linear shallow-water equations on the rotating sphere, a simplified system of equations used to study properties of space and time discretization methods in the context of atmospheric simulations. After introducing time integrators, we first compare the time step sizes to the errors in the simulation, discussing pros and cons of different formulations of REXI. Here, REXI already shows superior properties compared to explicit and implicit time stepping methods. Additionally, we present wallclock-time-to-error results revealing the sweet spots of REXI obtaining either an over 6x higher accuracy within the same time frame or an about 3x reduced time-to-solution for a similar error threshold. Our results motivate further explorations of REXI for operational weather/climate systems.

Stellar magnetism plays an important role in stellar evolution theory. Approximatively 10% of observed main sequence (MS) and pre-main-sequence (PMS) radiative stars exhibit surface magnetic fields above the detection limit, raising the question of their origin. These stars host outer radiative envelopes, which are stably stratified. Therefore, they are assumed to be motionless in standard models of stellar structure and evolution. We focus on rapidly rotating, radiative stars which may be prone to the tidal instability, due to an orbital companion. Using direct numerical simulations in a sphere, we study the interplay between a stable stratification and the tidal instability, and assess its dynamo capability. We show that the tidal instability is triggered regardless of the strength of the stratification (Brunt-V{\"a}is{\"a}l{\"a} frequency). Furthermore, the tidal instability can lead to both mixing and self-induced magnetic fields in stably stratified layers (provided that the Brunt-V{\"a}is{\"a}l{\"a} frequency does not exceed the stellar spin rate in the simulations too much). The application to stars suggests that the resulting magnetic fields could be observable at the stellar surfaces. Indeed, we expect magnetic field strengths up to several Gauss.Consequently, tidally driven dynamos should be considered as a (complementary) dynamo mechanism, possibly operating in radiative MS and PMS stars hosting orbital companions. In particular, tidally driven dynamos may explain the observed magnetism of tidally deformed and rapidly rotating Vega-like stars.

We investigate free hydromagnetic eigenmodes of an incompressible, inviscid and ideal electrically conducting fluid in rotating triaxial ellipsoids. The container rotates with an angular velocity tilted from its figure. The magnetic base state is a uniform current density also tilted. Three-dimensional perturbations upon the base state are expanded onto a finite-dimensional polynomial basis. By combining symbolic and numerical computations, we are able to get the eigenmodes of high spatial complexity. Hydromagnetic modes of the sphere still exist in triaxial geometry. A plane-wave analysis is also carried on, explaining the dispersion relation observed in our model. Without magnetic field, the modes reduce to the inertial modes of the ellipsoids, which form a complete basis. We propose to use these modes to study the weakly nonlinear saturation of inertial instabilities, especially the elliptical one.

We present an attempt to reach realistic turbulent regime in direct numerical simulations of the geodynamo. We rely on a sequence of three convection-driven simulations in a rapidly rotating spherical shell. The most extreme case reaches towards the Earth's core regime by lowering viscosity (magnetic Prandtl number Pm=0.1) while maintaining vigorous convection (magnetic Reynolds number Rm>500) and rapid rotation (Ekman number E=1e-7), at the limit of what is feasible on today's supercomputers. A detailed and comprehensive analysis highlights several key features matching geomagnetic observations or dynamo theory predictions -- all present together in the same simulation -- but it also unveils interesting insights relevant for Earth's core dynamics.In this strong-field, dipole-dominated dynamo simulation, the magnetic energy is one order of magnitude larger than the kinetic energy. The spatial distribution of magnetic intensity is highly heterogeneous, and a stark dynamical contrast exists between the interior and the exterior of the tangent cylinder (the cylinder parallel to the axis of rotation that circumscribes the inner core).In the interior, the magnetic field is strongest, and is associated with a vigorous twisted polar vortex, whose dynamics may occasionally lead to the formation of a reverse polar flux patch at the surface of the shell. Furthermore, the strong magnetic field also allows accumulation of light material within the tangent cylinder, leading to stable stratification there. Torsional Alfv{\'e}n waves are frequently triggered in the vicinity of the tangent cylinder and propagate towards the equator.Outside the tangent cylinder, the magnetic field inhibits the growth of zonal winds and the kinetic energy is mostly non-zonal. Spatio-temporal analysis indicates that the low-frequency, non-zonal flow is quite geostrophic (columnar) and predominantly large-scale: an m=1 eddy spontaneously emerges in our most extreme simulations, without any heterogeneous boundary forcing.Our spatio-temporal analysis further reveals that (i) the low-frequency, large-scale flow is governed by a balance between Coriolis and buoyancy forces -- magnetic field and flow tend to align, minimizing the Lorentz force; (ii) the high-frequency flow obeys a balance between magnetic and Coriolis forces; (iii) the convective plumes mostly live at an intermediate scale, whose dynamics is driven by a 3-term 1 MAC balance -- involving Coriolis, Lorentz and buoyancy forces. However, small-scale (E^{1/3}) quasi-geostrophic convection is still observed in the regions of low magnetic intensity., Comment: 48 pages, 28 figures, accepted for publication in GJI

The Derviche Tourneur Sodium experiment, a spherical Couette magnetohydrodynamics experiment with liquid sodium as the medium and a dipole magnetic field imposed from the inner sphere, recently underwent upgrades to its diagnostics to better characterize the flow and induced magnetic fields with global rotation. In tandem with the upgrades, a set of direct numerical simulations were run with the xshells code to give a more complete view of the fluid and magnetic dynamics at various rotation rates of the inner and outer spheres. These simulations reveal several dynamic regimes, determined by the Rossby number. At positive differential rotation there is a regime of quasigeostrophic flow, with low levels of fluctuations near the outer sphere. Negative differential rotation shows a regime of what appear to be saturated hydrodynamic instabilities at low negative differential rotation, followed by a regime where filamentary structures develop at low latitudes and persist over five to ten differential rotation periods as they drift poleward. We emphasize that all these coherent structures emerge from turbulent flows. At least some of them seem to be related to linear instabilities of the mean flow. The simulated flows can produce the same measurements as those that the physical experiment can take, with signatures akin to those found in the experiment. This paper discusses the relation between the internal velocity structures of the flow and their magnetic signatures at the surface., Comment: 15 pages, 11 figures

Torsional Alfv{\'e}n waves propagating in the Earth's core have been inferred by inversion techniques applied to geomagnetic models. They appear to propagate across the core but vanish at the equator, exchanging angular momentum between core and mantle. Assuming axial symmetry, we find that an electrically conducting layer at the bottom of the mantle can lead to total absorption of torsional waves that reach the equator. We show that the reflection coefficient depends on G Br , where Br is the strength of the radial magnetic field at the equator, and G the conductance of the lower mantle there. With Br = 7e-4 T., torsional waves are completely absorbed when they hit the equator if G = 1.3e8 S. For larger or smaller G, reflection occurs. As G is increased above this critical value, there is less attenuation and more angular momentum exchange. Our finding dissociates efficient core-mantle coupling from strong ohmic dissipation in the mantle.

We test the ability of large scale velocity fields inferred from geomagnetic secular variation data to produce the global magnetic field of the Earth.Our kinematic dynamo calculations use quasi-geostrophic (QG) flows inverted from geomagnetic field models which, as such, incorporate flow structures that are Earth-like and may be important for the geodynamo.Furthermore, the QG hypothesis allows straightforward prolongation of the flow from the core surface to the bulk.As expected from previous studies, we check that a simple quasi-geostrophic flow is not able to sustain the magnetic field against ohmic decay.Additional complexity is then introduced in the flow, inspired by the action of the Lorentz force.Indeed, on centenial time-scales, the Lorentz force can balance the Coriolis force and strict quasi-geostrophy may not be the best ansatz.When the columnar flow is modified to account for the action of the Lorentz force, magnetic field is generated for Elsasser numbers larger than 0.25 and magnetic Reynolds numbers larger than 100.This suggests that our large scale flow captures the relevant features for the generation of the Earth's magnetic field and that the invisible small scale flow may not be directly involved in this process.Near the threshold, the resulting magnetic field is dominated by an axial dipole, with some reversed flux patches.Time-dependence is also considered, derived from principal component analysis applied to the inverted flows.We find that time periods from 120 to 50 years do not affect the mean growth rate of the kinematic dynamos.Finally we notice the footprint of the inner-core in the magnetic field generated deep in the bulk of the shell, although we did not include one in our computations.

Seismic waves sensitive to the outermost part of the Earth's liquid core seem to be affected by a stably stratified layer at the core-mantle boundary. Such a layer could have an observable signature in both long-term and short-term variations of the magnetic field of the Earth, which are used to probe the flow at the top of the core. Indeed, with the recent SWARM mission, it seems reasonable to be able to identify waves propagating in the core with period of several months, which may play an important role in the large-scale dynamics. In this paper, we characterize the influence of a stratified layer at the top of the core on deep quasi-geostrophic (Rossby) waves. We compute numerically the quasi-geostrophic eigenmodes of a rapidly rotating spherical shell, with a stably stratified layer near the outer boundary. Two simple models of stratification are taken into account, which are scaled with commonly adopted values of the Brunt-V{\"a}is{\"a}l{\"a} frequency in the Earth's core. In the absence of magnetic field, we find that both azimuthal wavelength and frequency of the eigenmodes control their penetration into the stratified layer: the higher the phase speed, the higher the permeability of the stratified layer to the wave motion. We also show that the theory developed by Takehiro and Lister [2001] for thermal convection extends to the whole family of Rossby waves in the core. Adding a magnetic field, the penetrative behaviour of the quasi-geostrophic modes (the so-called fast branch) is insensitive to the imposed magnetic field and only weakly sensitive to the precise shape of the stratification. Based on these results, the large-scale and high frequency modes (1 to 2 month periods) may be detectable in the geomagnetic data measured at the Earth's surface, especially in the equatorial area where the modes can be trapped.

The contribution of small scale turbulent fluctuations to the induction of mean magnetic field is investigated in our liquid sodium spherical Couette experiment with an imposed magnetic field.An inversion technique is applied to a large number of measurements at $Rm \approx 100$ to obtain radial profiles of the $\alpha$ and $\beta$ effects and maps of the mean flow.It appears that the small scale turbulent fluctuations can be modeled as a strong contribution to the magnetic diffusivity that is negative in the interior region and positive close to the outer shell.Direct numerical simulations of our experiment support these results.The lowering of the effective magnetic diffusivity by small scale fluctuations implies that turbulence can actually help to achieve self-generation of large scale magnetic fields., Comment: Rajout d'un erratum

We present a reconstruction of the mean axisymmetric azimuthal and meridional flows in the DTS liquid sodium experiment. The experimental device sets a spherical Couette flow enclosed between two concentric spherical shells where the inner sphere holds a strong dipolar magnet, which acts as a magnetic propeller when rotated. Measurements of the mean velocity, mean induced magnetic field and mean electric potentials have been acquired inside and outside the fluid for an inner sphere rotation rate of 9 Hz (Rm 28). Using the induction equation to relate all measured quantities to the mean flow, we develop a nonlinear least square inversion procedure to reconstruct a fully coherent solution of the mean velocity field. We also include in our inversion the response of the fluid layer to the non-axisymmetric time-dependent magnetic field that results from deviations of the imposed magnetic field from an axial dipole. The mean azimuthal velocity field we obtain shows super-rotation in an inner region close to the inner sphere where the Lorentz force dominates, which contrasts with an outer geostrophic region governed by the Coriolis force, but where the magnetic torque remains the driver. The meridional circulation is strongly hindered by the presence of both the Lorentz and the Coriolis forces. Nevertheless, it contributes to a significant part of the induced magnetic energy. Our approach sets the scene for evaluating the contribution of velocity and magnetic fluctuations to the mean magnetic field, a key question for dynamo mechanisms.

The flow of liquid metal inside the Earth's core produces the geomagnetic field and its time variations. Understanding the variability of those deep currents is crucial to improve the forecast of geomagnetic field variations, which affect human spacial and aeronautic activities. Moreover, it may provide relevant information on the core dynamics. The main goal of this study is to extract and characterize the leading variability modes of core flows over centennial periods, and to assess their statistical robustness. To this end, we use flows that we invert from two geomagnetic field models (gufm1 and COV-OBS), and apply Principal Component Analysis and Singular Value Decomposition of coupled fields. The quasi geostrophic (QG) flows inverted from both geomagnetic field models show similar features. However, COV-OBS has a less energetic mean and larger time variability. The statistical significance of flow components is tested from analyses performed on subareas of the whole domain. Bootstrapping methods are also used to extract significant flow features required by both gufm1 and COV-OBS. Three main empirical circulation modes emerge, simultaneously constrained by both geomagnetic field models and expected to be robust against the particular a priori used to build them (large scale QG dynamics). Mode 1 exhibits three large vortices at medium/high latitudes, with opposite circulation under the Atlantic and the Pacific hemispheres. Mode 2 interestingly accounts for most of the variations of the Earth's core angular momentum. In this mode, the regions close to the tangent cylinder and to the equator are correlated, and oscillate with a period between 80 and 90 years. Each of these two modes is energetic enough to alter the mean flow, sometimes reinforcing the eccentric gyre, and other times breaking it up into smaller circulations. The three main circulation modes added together to the mean flow account for about 70% of the flows variability, 90% of the root mean square total velocities, and 95% of the secular variation induced by the total flows. Direct physical interpretation of the computed modes is not straightforward. Nonethe-less, similarities found between the two first modes and time/spatial features identified in different studies of core dynamics, suggest that our approach can help to pinpoint the relevant physical processes inside the core on centennial timescales., Comment: 33 pages

Several teams have reported peculiar frequency spectra for flows in a spherical shell. To address their origin, we perform numerical simulations of the spherical Couette flow in a dipolar magnetic field, in the configuration of the DTS experiment. The frequency spectra computed from time-series of the induced magnetic field display similar bumpy spectra, where each bump corresponds to a given azimuthal mode number m. The bumps show up at moderate Reynolds number (2 600) if the time-series are long enough (>300 rotations of the inner sphere). We present a new method that permits to retrieve the dominant frequencies for individual mode numbers m, and to extract the modal structure of the full non-linear flow. The maps of the energy of the fluctuations and the spatio-temporal evolution of the velocity field suggest that fluctuations originate in the outer boundary layer. The threshold of instability if found at Re_c = 1 860. The fluctuations result from two coupled instabilities: high latitude B\"odewadt-type boundary layer instability, and secondary non-axisymmetric instability of a centripetal jet forming at the equator of the outer sphere. We explore the variation of the magnetic and kinetic energies with the input parameters, and show that a modified Elsasser number controls their evolution. We can thus compare with experimental determinations of these energies and find a good agreement. Because of the dipolar nature of the imposed magnetic field, the energy of magnetic fluctuations is much larger near the inner sphere, but their origin lies in velocity fluctuations that initiate in the outer boundary layer., Comment: 23 pages

In this paper, we report on very efficient algorithms for the spherical harmonic transform (SHT). Explicitly vectorized variations of the algorithm based on the Gauss-Legendre quadrature are discussed and implemented in the SHTns library which includes scalar and vector transforms. The main breakthrough is to achieve very efficient on-the-fly computations of the Legendre associated functions, even for very high resolutions, by taking advantage of the specific properties of the SHT and the advanced capabilities of current and future computers. This allows us to simultaneously and significantly reduce memory usage and computation time of the SHT. We measure the performance and accuracy of our algorithms. Even though the complexity of the algorithms implemented in SHTns are in $O(N^3)$ (where N is the maximum harmonic degree of the transform), they perform much better than any third party implementation, including lower complexity algorithms, even for truncations as high as N=1023. SHTns is available at https://bitbucket.org/nschaeff/shtns as open source software., Comment: 8 pages

Alfv\'en waves propagate in electrically conducting fluids in the presence of a magnetic field. Their reflection properties depend on the ratio between the kinematic viscosity and the magnetic diffusivity of the fluid, also known as the magnetic Prandtl number Pm. In the special case Pm=1, there is no reflection on an insulating, no-slip boundary, and the wave energy is entirely dissipated in the boundary layer. We investigate the consequences of this remarkable behaviour for the numerical modeling of torsional Alfv\'en waves (also known as torsional oscillations), which represent a special class of Alfv\'en waves, in rapidly rotating spherical shells. They consist of geostrophic motions and are thought to exist in the fluid cores of planets with internal magnetic field. In the geophysical limit Pm << 1, these waves are reflected at the core equator, where they are entirely absorbed for Pm=1. Our numerical calculations show that the reflection coefficient at the equator of these waves remains below 0.2 for Pm > 0.3, which is the range of values for which geodynamo numerical models operate. As a result, geodynamo models with no-slip boundary conditions cannot exhibit torsional oscillation normal modes., Comment: 12 pages

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

Rapidly rotating spherical kinematic dynamos are computed using the combination of a quasi geostrophic (QG) model for the velocity field and a classical spectral 3D code for the magnetic field. On one hand, the QG flow is computed in the equatorial plane of a sphere and corresponds to Rossby wave instabilities of a geostrophic internal shear layer produced by differential rotation. On the other hand, the induction equation is computed in the full sphere after a continuation of the QG flow along the rotation axis. Differential rotation and Rossby-wave propagation are the key ingredients of the dynamo process which can be interpreted in terms of $\alpha\Omega$ dynamo. Taking into account the quasi geostrophy of the velocity field to increase its time and space resolution enables us to exhibit numerical dynamos with very low Ekman (rapidly rotating) and Prandtl numbers (liquid metals) which are asymptotically relevant to model planetary core dynamos.

We study the destabilization of a shear layer, produced by differential rotation of a rotating axisymmetric container. For small forcing, this produces a shear layer, which has been studied by Stewartson and is almost invariant along the rotation axis. When the forcing increases, instabilities develop. To study the asymptotic regime (very low Ekman number $E$), we develop a quasi-geostrophic two-dimensional model, whose main original feature is to handle the mass conservation correctly, resulting in a divergent two-dimensional flow, and valid for any container provided that the top and bottom have finite slopes. We use it to derive scalings and asymptotic laws by a simple linear theory, extending the previous analyses to large slopes (as in a sphere), for which we find different scaling laws. For a flat container, the critical Rossby number for the onset of instability evolves as $E^{3/4}$ and may be understood as a Kelvin-Helmoltz shear instability. For a sloping container, the instability is a Rossby wave with a critical Rossby number proportional to $\beta E^{1/2}$, where $\beta$ is related to the slope. We also investigate the asymmetry between positive and negative differential rotation and propose corrections for finite Ekman and Rossby numbers. Implemented in a numerical code, our model allows us to study the onset over a broad range of parameters, determining the threshold but also other features such as the spatial structure. We also present a few experimental results, validating our model and showing its limits., Comment: 28 pages, 15 figures. * New version including discussion of the recent work of Hollerbach, and much more. * A sign error in the Ekman pumping has been corrected. This has almost no influence on the results presented in the previous version