Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 5.154 IF 5.154
  • IF 5-year value: 5.697 IF 5-year
    5.697
  • CiteScore value: 5.56 CiteScore
    5.56
  • SNIP value: 1.761 SNIP 1.761
  • IPP value: 5.30 IPP 5.30
  • SJR value: 3.164 SJR 3.164
  • Scimago H <br class='hide-on-tablet hide-on-mobile'>index value: 59 Scimago H
    index 59
  • h5-index value: 49 h5-index 49
Volume 6, issue 6
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.

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
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.

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.

Publications Copernicus
Download
Citation