High-definition turbulent geodynamo
Using the xshells code on thousands of computer cores allowed to simulate the generation of the Earth's magnetic field with unprecedented resolution.
Schaeffer et al. 2017 doi:10.1093/gji/ggx265
Powered by SHTns, the blazingly fast spherical harmonic transform library.
Xshells was designed for speed from the ground up.
Easy to use, with a nice user manual and example parameter files for geodynamo, double-diffusive convection, precession, spherical Couette, ...
Runs from your laptop to massively parallel supercomputers, using both MPI and OpenMP.
Distributed under the CeCILL License (compatible with GNU GPL). Everybody is free to use, modify, contribute.
Tools for post-processing and visualization in python and paraview are provided.
At high thermal diffusivity (like a liquid metal) and at high rotation rate (like a planetary core), the onset of convection is subcritial,
the flow is always turbulent, while retaining the columnar structure due to the Coriolis force.
Kaplan et al. 2017 doi:10.1103/PhysRevLett.119.094501
Guervilly et al. 2019 doi:10.1038/s41586-019-1301-5