skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

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 onmore » 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.« less

Authors:
 [1];  [2];  [2];  [2];  [3]
  1. Sandia National Lab. (SNL-CA), Livermore, CA (United States)
  2. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
  3. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Sandia National Laboratories (SNL), Albuquerque, NM, and Livermore, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1208879
Alternate Identifier(s):
OSTI ID: 1214225
Grant/Contract Number:  
AC52-06NA25396; AC04-94AL85000
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Geoscientific Model Development (Online)
Additional Journal Information:
Journal Volume: 8; Journal Issue: 4; Journal ID: ISSN 1991-9603
Publisher:
European Geosciences Union
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; 97 MATHEMATICS AND COMPUTING; 54 ENVIRONMENTAL SCIENCES

Citation Formats

Tezaur, I. K., Perego, M., Salinger, A. G., Tuminaro, R. S., and Price, S. F. Albany/FELIX: A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis. United States: N. p., 2015. Web. doi:10.5194/gmd-8-1197-2015.
Tezaur, I. K., Perego, M., Salinger, A. G., Tuminaro, R. S., & Price, S. F. Albany/FELIX: A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis. United States. doi:10.5194/gmd-8-1197-2015.
Tezaur, I. K., Perego, M., Salinger, A. G., Tuminaro, R. S., and Price, S. F. Mon . "Albany/FELIX: A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis". United States. doi:10.5194/gmd-8-1197-2015. https://www.osti.gov/servlets/purl/1208879.
@article{osti_1208879,
title = {Albany/FELIX: A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis},
author = {Tezaur, I. K. and Perego, M. and Salinger, A. G. and Tuminaro, R. S. and Price, S. F.},
abstractNote = {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.},
doi = {10.5194/gmd-8-1197-2015},
journal = {Geoscientific Model Development (Online)},
issn = {1991-9603},
number = 4,
volume = 8,
place = {United States},
year = {2015},
month = {4}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 8 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Exact solutions to the thermomechanically coupled shallow-ice approximation: effective tools for verification
journal, January 2007


Consistent approximations and boundary conditions for ice-sheet dynamics from a principle of least action
journal, January 2010

  • Dukowicz, John K.; Price, Stephen F.; Lipscomb, William H.
  • Journal of Glaciology, Vol. 56, Issue 197
  • DOI: 10.3189/002214310792447851

Modelling West Antarctic ice sheet growth and collapse through the past five million years
journal, March 2009


Greenland flow variability from ice-sheet-wide velocity mapping
journal, January 2010


Implementation and Initial Evaluation of the Glimmer Community Ice Sheet Model in the Community Earth System Model
journal, October 2013

  • Lipscomb, William H.; Fyke, Jeremy G.; Vizcaíno, Miren
  • Journal of Climate, Vol. 26, Issue 19
  • DOI: 10.1175/JCLI-D-12-00557.1

Semicoarsening Multigrid on Distributed Memory Machines
journal, January 2000

  • Brown, Peter N.; Falgout, Robert D.; Jones, Jim E.
  • SIAM Journal on Scientific Computing, Vol. 21, Issue 5
  • DOI: 10.1137/S1064827598339141

Adaptive mesh, finite volume modeling of marine ice sheets
journal, January 2013

  • Cornford, Stephen L.; Martin, Daniel F.; Graves, Daniel T.
  • Journal of Computational Physics, Vol. 232, Issue 1
  • DOI: 10.1016/j.jcp.2012.08.037

Amesos2 and Belos: Direct and Iterative Solvers for Large Sparse Linear Systems
journal, January 2012

  • Bavier, Eric; Hoemmen, Mark; Rajamanickam, Sivasankaran
  • Scientific Programming, Vol. 20, Issue 3
  • DOI: 10.1155/2012/243875

A parallel high-order accurate finite element nonlinear Stokes ice sheet model and benchmark experiments: A PARALLEL FEM STOKES ICE SHEET MODEL
journal, January 2012

  • Leng, Wei; Ju, Lili; Gunzburger, Max
  • Journal of Geophysical Research: Earth Surface, Vol. 117, Issue F1
  • DOI: 10.1029/2011JF001962

Manufactured solutions and the verification of three-dimensional Stokes ice-sheet models
journal, January 2013


Enhanced basal lubrication and the contribution of the Greenland ice sheet to future sea-level rise
journal, August 2013

  • Shannon, S. R.; Payne, A. J.; Bartholomew, I. D.
  • Proceedings of the National Academy of Sciences, Vol. 110, Issue 35
  • DOI: 10.1073/pnas.1212647110

Thin-Film Flows with Wall Slip: An Asymptotic Analysis of Higher Order Glacier Flow Models
journal, January 2010

  • Schoof, C.; Hindmarsh, R. C. A.
  • The Quarterly Journal of Mechanics and Applied Mathematics, Vol. 63, Issue 1
  • DOI: 10.1093/qjmam/hbp025

Coupling ice flow models of varying orders of complexity with the Tiling method
journal, January 2012

  • Seroussi, Helene; Ben Dhia, Hachmi; Morlighem, Mathieu
  • Journal of Glaciology, Vol. 58, Issue 210
  • DOI: 10.3189/2012JoG11J195

Greenland ice sheet contribution to sea-level rise from a new-generation ice-sheet model
journal, January 2012


Toward a New Generation of Ice Sheet Models
journal, January 2007

  • Little, Christopher M.; Oppenheimer, Michael; Alley, Richard B.
  • Eos, Transactions American Geophysical Union, Vol. 88, Issue 52
  • DOI: 10.1029/2007EO520002

Achieving Textbook Multigrid Efficiency for Hydrostatic Ice Sheet Flow
journal, January 2013

  • Brown, Jed; Smith, Barry; Ahmadia, Aron
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 2
  • DOI: 10.1137/110834512

Solving PDEs with Intrepid
journal, January 2012

  • Bochev, P.; Edwards, H. C.; Kirby, R. C.
  • Scientific Programming, Vol. 20, Issue 2
  • DOI: 10.1155/2012/403902

Bifurcation Tracking Algorithms and Software for Large Scale Applications
journal, March 2005

  • Salinger, A. G.; Burroughs, E. A.; Pawlowski, R. P.
  • International Journal of Bifurcation and Chaos, Vol. 15, Issue 03
  • DOI: 10.1142/S0218127405012508

An overview of the Trilinos project
journal, September 2005

  • Heroux, Michael A.; Phipps, Eric T.; Salinger, Andrew G.
  • ACM Transactions on Mathematical Software, Vol. 31, Issue 3
  • DOI: 10.1145/1089014.1089021

An ice-shelf model test based on the Ross Ice Shelf, Antarctica
journal, January 1996


Data assimilation using a hybrid ice flow model
journal, January 2011


Transient glacier response with a higher-order numerical ice-flow model
journal, January 2002


Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM): ICE SHEET SYSTEM MODEL
journal, March 2012

  • Larour, E.; Seroussi, H.; Morlighem, M.
  • Journal of Geophysical Research: Earth Surface, Vol. 117, Issue F1
  • DOI: 10.1029/2011JF002140

Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model
journal, January 2009

  • Bueler, Ed; Brown, Jed
  • Journal of Geophysical Research, Vol. 114, Issue F3
  • DOI: 10.1029/2008JF001179

Globalization Techniques for Newton–Krylov Methods and Applications to the Fully Coupled Solution of the Navier–Stokes Equations
journal, January 2006


The Glimmer community ice sheet model
journal, January 2009

  • Rutt, I. C.; Hagdorn, M.; Hulton, N. R. J.
  • Journal of Geophysical Research, Vol. 114, Issue F2
  • DOI: 10.1029/2008JF001015

Parameter and state estimation with a time-dependent adjoint marine ice sheet model
journal, January 2013


Automating Embedded Analysis Capabilities and Managing Software Complexity in Multiphysics Simulation, Part I: Template-Based Generic Programming
journal, January 2012

  • Pawlowski, Roger P.; Phipps, Eric T.; Salinger, Andrew G.
  • Scientific Programming, Vol. 20, Issue 2
  • DOI: 10.1155/2012/202071

Optimal initial conditions for coupling ice sheet models to Earth system models: PEREGO ET AL.
journal, September 2014

  • Perego, Mauro; Price, Stephen; Stadler, Georg
  • Journal of Geophysical Research: Earth Surface, Vol. 119, Issue 9
  • DOI: 10.1002/2014JF003181

Parallel finite-element implementation for higher-order ice-sheet models
journal, January 2012

  • Perego, Mauro; Gunzburger, Max; Burkardt, John
  • Journal of Glaciology, Vol. 58, Issue 207
  • DOI: 10.3189/2012JoG11J063

Automating Embedded Analysis Capabilities and Managing Software Complexity in Multiphysics Simulation, Part II: Application to Partial Differential Equations
journal, January 2012

  • Pawlowski, Roger P.; Phipps, Eric T.; Salinger, Andrew G.
  • Scientific Programming, Vol. 20, Issue 3
  • DOI: 10.1155/2012/818262

Implementation of the Jacobian-free Newton–Krylov method for solving the first-order ice sheet momentum balance
journal, July 2011

  • Lemieux, Jean-François; Price, Stephen F.; Evans, Katherine J.
  • Journal of Computational Physics, Vol. 230, Issue 17
  • DOI: 10.1016/j.jcp.2011.04.037

Well-Posedness Results for a Nonlinear Stokes Problem Arising in Glaciology
journal, January 2013

  • Chen, Qingshan; Gunzburger, Max; Perego, Mauro
  • SIAM Journal on Mathematical Analysis, Vol. 45, Issue 5
  • DOI: 10.1137/110848694

Quantum computer aided design simulation and optimization of semiconductor quantum dots
journal, October 2013

  • Gao, X.; Nielsen, E.; Muller, R. P.
  • Journal of Applied Physics, Vol. 114, Issue 16
  • DOI: 10.1063/1.4825209

A Community Ice Sheet Model for Sea Level Prediction: Building a Next-Generation Community Ice Sheet Model; Los Alamos, New Mexico, 18–20 August 2008
journal, January 2009

  • Lipscomb, William; Bindschadler, Robert; Bueler, Ed
  • Eos, Transactions American Geophysical Union, Vol. 90, Issue 3
  • DOI: 10.1029/2009EO030004

A stabilized assumed deformation gradient finite element formulation for strongly coupled poromechanical simulations at finite strain: STABILIZED F-BAR FINITE ELEMENT FORMULATION FOR POROMECHANICS
journal, January 2013

  • Sun, WaiChing; Ostien, Jakob T.; Salinger, Andrew G.
  • International Journal for Numerical and Analytical Methods in Geomechanics
  • DOI: 10.1002/nag.2161

A new bed elevation dataset for Greenland
journal, January 2013

  • Bamber, J. L.; Griggs, J. A.; Hurkmans, R. T. W. L.
  • The Cryosphere, Vol. 7, Issue 2
  • DOI: 10.5194/tc-7-499-2013

Roles of marine ice, rheology, and fracture in the flow and stability of the Brunt/Stancomb-Wills Ice Shelf
journal, January 2009

  • Khazendar, A.; Rignot, E.; Larour, E.
  • Journal of Geophysical Research, Vol. 114, Issue F4
  • DOI: 10.1029/2008JF001124

Benchmark experiments for higher-order and full-Stokes ice sheet models (ISMIP–HOM)
journal, January 2008

  • Pattyn, F.; Perichon, L.; Aschwanden, A.
  • The Cryosphere, Vol. 2, Issue 2
  • DOI: 10.5194/tc-2-95-2008

Coulomb Friction and Other Sliding laws in a Higher-Order Glacier flow Model
journal, January 2010


Capabilities and performance of Elmer/Ice, a new-generation ice sheet model
journal, January 2013

  • Gagliardini, O.; Zwinger, T.; Gillet-Chaulet, F.
  • Geoscientific Model Development, Vol. 6, Issue 4
  • DOI: 10.5194/gmd-6-1299-2013

An adaptive Newton multigrid method for a model of marine ice sheets
journal, November 2013


An ice-shelf model test based on the Ross Ice Shelf, Antarctica
journal, January 1996


    Works referencing / citing this record: