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 8, issue 4
Geosci. Model Dev., 8, 1197-1220, 2015
https://doi.org/10.5194/gmd-8-1197-2015
© Author(s) 2015. This work is distributed under
the Creative Commons Attribution 3.0 License.
Geosci. Model Dev., 8, 1197-1220, 2015
https://doi.org/10.5194/gmd-8-1197-2015
© Author(s) 2015. This work is distributed under
the Creative Commons Attribution 3.0 License.

Model description paper 27 Apr 2015

Model description paper | 27 Apr 2015

Albany/FELIX: a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis

I. K. Tezaur1, M. Perego2, A. G. Salinger2, R. S. Tuminaro2, and S. F. Price3 I. K. Tezaur et al.
  • 1Quantitative Modeling and Analysis Department, Sandia National Laboratories, P.O. Box 969, MS 9159, Livermore, CA 94551, USA
  • 2Computational Mathematics Department, Sandia National Laboratories, P.O. Box 5800, MS 1320, Albuquerque, NM 87185, USA
  • 3Fluid Dynamics and Solid Mechanics Group, Los Alamos National Laboratory, P.O. Box 1663, MS B216, Los Alamos, NM 87545, USA

Abstract. This paper describes a new parallel, scalable and robust finite element based solver for the first-order Stokes momentum balance equations for ice flow. The solver, known as Albany/FELIX, is constructed using the component-based approach to building application codes, in which mature, modular libraries developed as a part of the Trilinos project are combined using abstract interfaces and template-based generic programming, resulting in a final code with access to dozens of algorithmic and advanced analysis capabilities. Following an overview of the relevant partial differential equations and boundary conditions, the numerical methods chosen to discretize the ice flow equations are described, along with their implementation. The results of several verification studies of the model accuracy are presented using (1) new test cases for simplified two-dimensional (2-D) versions of the governing equations derived using the method of manufactured solutions, and (2) canonical ice sheet modeling benchmarks. Model accuracy and convergence with respect to mesh resolution are then studied on problems involving a realistic Greenland ice sheet geometry discretized using hexahedral and tetrahedral meshes. Also explored as a part of this study is the effect of vertical mesh resolution on the solution accuracy and solver performance. The robustness and scalability of our solver on these problems is demonstrated. Lastly, we show that good scalability can be achieved by preconditioning the iterative linear solver using a new algebraic multilevel preconditioner, constructed based on the idea of semi-coarsening.

Publications Copernicus
Download
Short summary
In this manuscript, we discuss the development and validation of a new momentum balance solver for modeling the flow of glaciers and ice sheets based on the 1st-order Stokes equations. We demonstrate the numerical convergence of our solver (with respect to computational mesh spacing), its flexibility (with respect to both the choice of mesh and finite element type), and its computational performance (robustness and scalability when applied to both idealized and realistic ice sheet simulations).
In this manuscript, we discuss the development and validation of a new momentum balance solver...
Citation
Share