Towards automatic finite-element methods for geodynamics via Firedrake

• Published in: Geoscientific Model Development Vol. 15; no. 13; pp. 5127 - 5166
• Main Authors: Davies, D. Rhodri, Kramer, Stephan C, Ghelichkhan, Sia, Gibson, Angus
• Format: Journal Article
• Language: English
• Published: Katlenburg-Lindau Copernicus GmbH 05.07.2022; Copernicus Publications
• Subjects: Algorithms; Analysis; Approximation; Automation; Boundary conditions; C plus plus; Cartesian coordinates; Collaboration; Complexity; Compressibility; Convection; Cores; Differential equations; Earth; Finite element analysis; Finite element method; Flexibility; Fluid dynamics; Geodynamics; Geophysical fluids; Hydrodynamics; Libraries; Mantle; Mantle convection; Mathematical models; Methods; Partial differential equations; Physics; Plate motion; Plate tectonics; Rheological properties; Simulation; Software; Tectonophysics
• ISSN: 1991-9603; 1991-959X; 1991-962X; 1991-9603; 1991-962X
• DOI: 10.5194/gmd-15-5127-2022

Firedrake is an automated system for solving partial differential equations using the finite-element method. By applying sophisticated performance optimisations through automatic code-generation techniques, it provides a means of creating accurate, efficient, flexible, easily extensible, scalable, transparent and reproducible research software that is ideally suited to simulating a wide range of problems in geophysical fluid dynamics. Here, we demonstrate the applicability of Firedrake for geodynamical simulation, with a focus on mantle dynamics. The accuracy and efficiency of the approach are confirmed via comparisons against a suite of analytical and benchmark cases of systematically increasing complexity, whilst parallel scalability is demonstrated up to 12 288 compute cores, where the problem size and the number of processing cores are simultaneously increased. In addition, Firedrake's flexibility is highlighted via straightforward application to different physical (e.g. complex non-linear rheologies, compressibility) and geometrical (2-D and 3-D Cartesian and spherical domains) scenarios. Finally, a representative simulation of global mantle convection is examined, which incorporates 230 Myr of plate motion history as a kinematic surface boundary condition, confirming Firedrake's suitability for addressing research problems at the frontiers of global mantle dynamics research.