Articles | Volume 6, issue 6
Geosci. Model Dev., 6, 2099–2119, 2013
https://doi.org/10.5194/gmd-6-2099-2013

Special issue: Isaac Newton Institute programme on multiscale numerics for...

Geosci. Model Dev., 6, 2099–2119, 2013
https://doi.org/10.5194/gmd-6-2099-2013

Development and technical paper 17 Dec 2013

Development and technical paper | 17 Dec 2013

Automating the solution of PDEs on the sphere and other manifolds in FEniCS 1.2

M. E. Rognes1, D. A. Ham2,3, C. J. Cotter2, and A. T. T. McRae2,4 M. E. Rognes et al.
  • 1Center for Biomedical Computing, Simula Research Laboratory, P.O. Box 134, 1325 Lysaker, Norway
  • 2Department of Mathematics, Imperial College London, London SW7 2AZ, UK
  • 3Department of Computing, Imperial College London, London SW7 2AZ, UK
  • 4The Grantham Institute for Climate Change, Imperial College London, London SW7 2AZ, UK

Abstract. Differential equations posed over immersed manifolds are of particular importance in studying geophysical flows; for instance, ocean and atmosphere simulations crucially rely on the capability to solve equations over the sphere. This paper presents the extension of the FEniCS software components to the automated solution of finite element formulations of differential equations defined over general, immersed manifolds. We describe the implementation and, in particular detail, how the required extensions essentially reduce to the extension of the FEniCS form compiler to cover this case. The resulting implementation has all the properties of the FEniCS pipeline and we demonstrate its flexibility by an extensive range of numerical examples covering a number of geophysical benchmark examples and test cases. The results are all in agreement with the expected values. The description here relates to DOLFIN/FEniCS 1.2.