Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Geosci. Model Dev., 6, 2099-2119, 2013
https://doi.org/10.5194/gmd-6-2099-2013
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.
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 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.

Citation: Rognes, M. E., Ham, D. A., Cotter, C. J., and McRae, A. T. T.: Automating the solution of PDEs on the sphere and other manifolds in FEniCS 1.2, Geosci. Model Dev., 6, 2099-2119, https://doi.org/10.5194/gmd-6-2099-2013, 2013.
Publications Copernicus
Download
Share