On the scalability of the Albany/FELIX first-order Stokes approximation ice sheet solver for large-scale simulations of the Greenland and Antarctic ice sheets
Abstract
We examine the scalability of the recently developed Albany/FELIX finite-element based code for the first-order Stokes momentum balance equations for ice flow. We focus our analysis on the performance of two possible preconditioners for the iterative solution of the sparse linear systems that arise from the discretization of the governing equations: (1) a preconditioner based on the incomplete LU (ILU) factorization, and (2) a recently-developed algebraic multigrid (AMG) preconditioner, constructed using the idea of semi-coarsening. A strong scalability study on a realistic, high resolution Greenland ice sheet problem reveals that, for a given number of processor cores, the AMG preconditioner results in faster linear solve times but the ILU preconditioner exhibits better scalability. In addition, a weak scalability study is performed on a realistic, moderate resolution Antarctic ice sheet problem, a substantial fraction of which contains floating ice shelves, making it fundamentally different from the Greenland ice sheet problem. We show that as the problem size increases, the performance of the ILU preconditioner deteriorates whereas the AMG preconditioner maintains scalability. This is because the linear systems are extremely ill-conditioned in the presence of floating ice shelves, and the ill-conditioning has a greater negative effect on the ILU preconditioner than onmore »
- Authors:
-
- Sandia National Lab. (SNL-CA), Livermore, CA (United States)
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Publication Date:
- Research Org.:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States); Sandia National Laboratories (SNL), Albuquerque, NM, and Livermore, CA (United States)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1201743
- Alternate Identifier(s):
- OSTI ID: 1214669; OSTI ID: 1214670
- Grant/Contract Number:
- AC04-94AL85000; AC52-06NA25396
- Resource Type:
- Journal Article: Accepted Manuscript
- Journal Name:
- Procedia Computer Science
- Additional Journal Information:
- Journal Volume: 51; Journal Issue: C; Journal ID: ISSN 1877-0509
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 54 ENVIRONMENTAL SCIENCES; 97 MATHEMATICS AND COMPUTING; ice sheet model; first-order Stokes approximation; finite element method; scalability; ILU preconditioner; algebraic multigrid (AMG) preconditioner; semi-coarsening; Greenland; Antarctica; Algebraic multigrid (AMG) preconditioner
Citation Formats
Tezaur, Irina K., Tuminaro, Raymond S., Perego, Mauro, Salinger, Andrew G., and Price, Stephen F. On the scalability of the Albany/FELIX first-order Stokes approximation ice sheet solver for large-scale simulations of the Greenland and Antarctic ice sheets. United States: N. p., 2015.
Web. doi:10.1016/j.procs.2015.05.467.
Tezaur, Irina K., Tuminaro, Raymond S., Perego, Mauro, Salinger, Andrew G., & Price, Stephen F. On the scalability of the Albany/FELIX first-order Stokes approximation ice sheet solver for large-scale simulations of the Greenland and Antarctic ice sheets. United States. https://doi.org/10.1016/j.procs.2015.05.467
Tezaur, Irina K., Tuminaro, Raymond S., Perego, Mauro, Salinger, Andrew G., and Price, Stephen F. 2015.
"On the scalability of the Albany/FELIX first-order Stokes approximation ice sheet solver for large-scale simulations of the Greenland and Antarctic ice sheets". United States. https://doi.org/10.1016/j.procs.2015.05.467. https://www.osti.gov/servlets/purl/1201743.
@article{osti_1201743,
title = {On the scalability of the Albany/FELIX first-order Stokes approximation ice sheet solver for large-scale simulations of the Greenland and Antarctic ice sheets},
author = {Tezaur, Irina K. and Tuminaro, Raymond S. and Perego, Mauro and Salinger, Andrew G. and Price, Stephen F.},
abstractNote = {We examine the scalability of the recently developed Albany/FELIX finite-element based code for the first-order Stokes momentum balance equations for ice flow. We focus our analysis on the performance of two possible preconditioners for the iterative solution of the sparse linear systems that arise from the discretization of the governing equations: (1) a preconditioner based on the incomplete LU (ILU) factorization, and (2) a recently-developed algebraic multigrid (AMG) preconditioner, constructed using the idea of semi-coarsening. A strong scalability study on a realistic, high resolution Greenland ice sheet problem reveals that, for a given number of processor cores, the AMG preconditioner results in faster linear solve times but the ILU preconditioner exhibits better scalability. In addition, a weak scalability study is performed on a realistic, moderate resolution Antarctic ice sheet problem, a substantial fraction of which contains floating ice shelves, making it fundamentally different from the Greenland ice sheet problem. We show that as the problem size increases, the performance of the ILU preconditioner deteriorates whereas the AMG preconditioner maintains scalability. This is because the linear systems are extremely ill-conditioned in the presence of floating ice shelves, and the ill-conditioning has a greater negative effect on the ILU preconditioner than on the AMG preconditioner.},
doi = {10.1016/j.procs.2015.05.467},
url = {https://www.osti.gov/biblio/1201743},
journal = {Procedia Computer Science},
issn = {1877-0509},
number = C,
volume = 51,
place = {United States},
year = {Thu Jan 01 00:00:00 EST 2015},
month = {Thu Jan 01 00:00:00 EST 2015}
}
Web of Science
Works referenced in this record:
Velocity and stress fields in grounded glaciers: a simple algorithm for including deviatoric stress gradients
journal, January 1995
- Blatter, Heinz
- Journal of Glaciology, Vol. 41, Issue 138
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
Bedmap2: improved ice bed, surface and thickness datasets for Antarctica
journal, January 2013
- Fretwell, P.; Pritchard, H. D.; Vaughan, D. G.
- The Cryosphere, Vol. 7, Issue 1
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
Albany/FELIX: a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis
posted_content, November 2014
- Kalashnikova, I.; Perego, M.; Salinger, A. G.
- Geoscientific Model Development
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
Ice Flow of the Antarctic Ice Sheet
journal, August 2011
- Rignot, E.; Mouginot, J.; Scheuchl, B.
- Science, Vol. 333, Issue 6048
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
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
Works referencing / citing this record:
MPAS-Albany Land Ice (MALI): a variable-resolution ice sheet model for Earth system modeling using Voronoi grids
journal, January 2018
- Hoffman, Matthew J.; Perego, Mauro; Price, Stephen F.
- Geoscientific Model Development, Vol. 11, Issue 9
A Study on the Performance Portability of the Finite Element Assembly Process Within the Albany Land Ice Solver
book, February 2020
- Watkins, Jerry; Tezaur, Irina; Demeshko, Irina
- Numerical Methods for Flows: FEF 2017 Selected Contributions, p. 177-188
LIVVkit 2.1: automated and extensible ice sheet model validation
journal, January 2019
- Evans, Katherine J.; Kennedy, Joseph H.; Lu, Dan
- Geoscientific Model Development, Vol. 12, Issue 3
Progress in Fast, Accurate Multi-scale Climate Simulations
journal, January 2015
- Collins, W. D.; Johansen, H.; Evans, K. J.
- Procedia Computer Science, Vol. 51
Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11)
journal, January 2017
- Habbal, Feras; Larour, Eric; Morlighem, Mathieu
- Geoscientific Model Development, Vol. 10, Issue 1
An ice sheet model validation framework for the Greenland ice sheet
journal, January 2017
- Price, Stephen F.; Hoffman, Matthew J.; Bonin, Jennifer A.
- Geoscientific Model Development, Vol. 10, Issue 1
Interactive comment on “An ice sheet model validation framework for the Greenland ice sheet” by Stephen F. Price et al.
journal, August 2016
- Price, Stephen F.
- Geoscientific Model Development Discussions