Computational geodynamics to model mantle convection and magma oceans with the Thermal Lattice Boltzmann Method and strong scaling on exascale clusters

Summary The Thermal Lattice Boltzmann Method (TLBM) offers an alternative to classical PDE based numerical methods for modelling planetary dynamics. It solves the Lattice Boltzmann Equation on a unitary discrete lattice which requires no matrix operations, and hence scales extremely well on HPC clusters. We describe the TLBM equations which are capable of modelling solid state mantle convection and turbulent magma ocean dynamics, and we present how to obtain the physical (dimensional) time scale from the lattice time scale (number of time steps), convert quantities in physical (dimensional) units (η, κ, σ…) to unitary (dimensionless) lattice units, and strong scaling performance in 3D to 300K cores. We present examples to illustrate use of the TLBM to model plate tectonics, 2D and 3D mantle convection, turbulent magma ocean dynamics, and planetary accretion. We present a table of run times which demonstrate the TLBM’s high throughput performance. The TLBM’s performance enables 3D global convection modelling of the whole mantle on a lattice with up to ∼2 × 1011 grid points corresponding to a resolution of 3 km, and magma ocean dynamics modelling into the ultimate turbulent regime of Ra ≳ 1014, to at least Ra = 1015 (demonstrated to date in 2D for an aspect ratio of 1), and projected to be possible up to Ra = 1018 in 2D or Ra = 1015 in 3D for low aspect ratio models. The TLBM is fast and makes it possible to run huge models which opens exascale computing to planetary dynamics research. We believe that the TLBM provides a valuable new means to advance geodynamical research into the future by enabling fast and high-resolution planetary dynamics simulations including of global mantle convection, magma ocean dynamics into the ultimate turbulent regime, and planetary accretion.