Magnetic boundary layers in numerical dynamos with heterogeneous outer boundary heat flux

Terra-Nova, Filipe; Amit, Hagay

doi:10.1016/j.pepi.2020.106589

Cited by 5 publications

(1 citation statement)

References 94 publications

(158 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…For instance, MagIC was used to study the "top-heavy" regime of double-diffusive convection, when thermal and compositional background gradients are destabilizing (Tassin et al, 2021). In a recent paper, MagIC has been used to study magnetic boundary layers with heterogeneous outer boundary heat flux, with results suggesting a significant deviation from classical estimates of the diffusion time and the magnetic Reynold number (Terra-Nova and Amit, 2020). In another paper, the dynamics of a possible stable stratification layer atop Earth's core was explored (e.g.…”

Section: Introductionmentioning

confidence: 99%

MagIC v5.10: a two-dimensional MPI distribution for pseudo-spectral magneto hydrodynamics simulations in spherical geometry

Lago

Gastine

Dannert

et al. 2021

Preprint

View full text Add to dashboard Cite

Abstract. We discuss two parallelization schemes for MagIC, an open-source, high-performance, pseudo-spectral code for the numerical solution of the magneto hydrodynamics equations in a rotating spherical shell. MagIC calculates the non-linear terms on a numerical grid in spherical coordinates while the time step updates are performed on radial grid points with a spherical harmonic representation of the lateral directions. Several transforms are required to switch between the different representations. The established hybrid implementation of MagIC uses MPI-parallelization in radius and relies on existing fast spherical transforms using OpenMP. Our new two-dimensional MPI decomposition implementation also distributes the latitudes or the azimuthal wavenumbers across the available MPI tasks/compute cores. We discuss several non-trivial algorithmic optimizations and the different data distribution layouts employed by our scheme. In particular, the two-dimensional distribution data layout yields a code that strongly scales well beyond the limit of the current one-dimensional distribution. We also show that the two-dimensional distribution implementation, although not yet fully optimized, can already be faster than the existing finely optimized hybrid implementation when using many thousands of CPU cores. Our analysis indicates that the two-dimensional distribution variant can be further optimized to also surpass the performance of the one-dimensional distribution for a few thousand cores.

show abstract

Section: Introductionmentioning

confidence: 99%