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
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...
Published in: | Procedia Computer Science |
---|---|
Main Authors: | , , , , |
Language: | unknown |
Published: |
2023
|
Subjects: | |
Online Access: | http://www.osti.gov/servlets/purl/1201743 https://www.osti.gov/biblio/1201743 https://doi.org/10.1016/j.procs.2015.05.467 |
id |
ftosti:oai:osti.gov:1201743 |
---|---|
record_format |
openpolar |
spelling |
ftosti:oai:osti.gov:1201743 2023-07-30T03:56:40+02:00 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 Tezaur, Irina K. Tuminaro, Raymond S. Perego, Mauro Salinger, Andrew G. Price, Stephen F. 2023-06-26 application/pdf http://www.osti.gov/servlets/purl/1201743 https://www.osti.gov/biblio/1201743 https://doi.org/10.1016/j.procs.2015.05.467 unknown http://www.osti.gov/servlets/purl/1201743 https://www.osti.gov/biblio/1201743 https://doi.org/10.1016/j.procs.2015.05.467 doi:10.1016/j.procs.2015.05.467 54 ENVIRONMENTAL SCIENCES 97 MATHEMATICS AND COMPUTING 2023 ftosti https://doi.org/10.1016/j.procs.2015.05.467 2023-07-11T09:01:32Z 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. 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. Here, 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. Other/Unknown Material Antarc* Antarctic Greenland Ice Sheet Ice Shelves SciTec Connect (Office of Scientific and Technical Information - OSTI, U.S. Department of Energy) Antarctic Greenland Procedia Computer Science 51 2026 2035 |
institution |
Open Polar |
collection |
SciTec Connect (Office of Scientific and Technical Information - OSTI, U.S. Department of Energy) |
op_collection_id |
ftosti |
language |
unknown |
topic |
54 ENVIRONMENTAL SCIENCES 97 MATHEMATICS AND COMPUTING |
spellingShingle |
54 ENVIRONMENTAL SCIENCES 97 MATHEMATICS AND COMPUTING 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 |
topic_facet |
54 ENVIRONMENTAL SCIENCES 97 MATHEMATICS AND COMPUTING |
description |
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. 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. Here, 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. |
author |
Tezaur, Irina K. Tuminaro, Raymond S. Perego, Mauro Salinger, Andrew G. Price, Stephen F. |
author_facet |
Tezaur, Irina K. Tuminaro, Raymond S. Perego, Mauro Salinger, Andrew G. Price, Stephen F. |
author_sort |
Tezaur, Irina K. |
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 |
title_short |
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 |
title_full |
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 |
title_fullStr |
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 |
title_full_unstemmed |
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 |
title_sort |
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 |
publishDate |
2023 |
url |
http://www.osti.gov/servlets/purl/1201743 https://www.osti.gov/biblio/1201743 https://doi.org/10.1016/j.procs.2015.05.467 |
geographic |
Antarctic Greenland |
geographic_facet |
Antarctic Greenland |
genre |
Antarc* Antarctic Greenland Ice Sheet Ice Shelves |
genre_facet |
Antarc* Antarctic Greenland Ice Sheet Ice Shelves |
op_relation |
http://www.osti.gov/servlets/purl/1201743 https://www.osti.gov/biblio/1201743 https://doi.org/10.1016/j.procs.2015.05.467 doi:10.1016/j.procs.2015.05.467 |
op_doi |
https://doi.org/10.1016/j.procs.2015.05.467 |
container_title |
Procedia Computer Science |
container_volume |
51 |
container_start_page |
2026 |
op_container_end_page |
2035 |
_version_ |
1772814085419696128 |