Albany/FELIX : A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis
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 libra...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , , , , |
Language: | unknown |
Published: |
2023
|
Subjects: | |
Online Access: | http://www.osti.gov/servlets/purl/1214225 https://www.osti.gov/biblio/1214225 https://doi.org/10.5194/gmd-8-1197-2015 |
id |
ftosti:oai:osti.gov:1214225 |
---|---|
record_format |
openpolar |
spelling |
ftosti:oai:osti.gov:1214225 2023-07-30T04:03:54+02:00 Albany/FELIX : A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis Tezaur, I. K. Perego, M. Salinger, A. G. Tuminaro, R. S. Price, S. F. 2023-06-26 application/pdf http://www.osti.gov/servlets/purl/1214225 https://www.osti.gov/biblio/1214225 https://doi.org/10.5194/gmd-8-1197-2015 unknown http://www.osti.gov/servlets/purl/1214225 https://www.osti.gov/biblio/1214225 https://doi.org/10.5194/gmd-8-1197-2015 doi:10.5194/gmd-8-1197-2015 97 MATHEMATICS AND COMPUTING 54 ENVIRONMENTAL SCIENCES 2023 ftosti https://doi.org/10.5194/gmd-8-1197-2015 2023-07-11T09:02:50Z 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. Other/Unknown Material Greenland Ice Sheet SciTec Connect (Office of Scientific and Technical Information - OSTI, U.S. Department of Energy) Greenland Geoscientific Model Development 8 4 1197 1220 |
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 |
97 MATHEMATICS AND COMPUTING 54 ENVIRONMENTAL SCIENCES |
spellingShingle |
97 MATHEMATICS AND COMPUTING 54 ENVIRONMENTAL SCIENCES 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 |
topic_facet |
97 MATHEMATICS AND COMPUTING 54 ENVIRONMENTAL SCIENCES |
description |
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. |
author |
Tezaur, I. K. Perego, M. Salinger, A. G. Tuminaro, R. S. Price, S. F. |
author_facet |
Tezaur, I. K. Perego, M. Salinger, A. G. Tuminaro, R. S. Price, S. F. |
author_sort |
Tezaur, I. K. |
title |
Albany/FELIX : A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis |
title_short |
Albany/FELIX : A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis |
title_full |
Albany/FELIX : A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis |
title_fullStr |
Albany/FELIX : A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis |
title_full_unstemmed |
Albany/FELIX : A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis |
title_sort |
albany/felix : a parallel, scalable and robust, finite element, first-order stokes approximation ice sheet solver built for advanced analysis |
publishDate |
2023 |
url |
http://www.osti.gov/servlets/purl/1214225 https://www.osti.gov/biblio/1214225 https://doi.org/10.5194/gmd-8-1197-2015 |
geographic |
Greenland |
geographic_facet |
Greenland |
genre |
Greenland Ice Sheet |
genre_facet |
Greenland Ice Sheet |
op_relation |
http://www.osti.gov/servlets/purl/1214225 https://www.osti.gov/biblio/1214225 https://doi.org/10.5194/gmd-8-1197-2015 doi:10.5194/gmd-8-1197-2015 |
op_doi |
https://doi.org/10.5194/gmd-8-1197-2015 |
container_title |
Geoscientific Model Development |
container_volume |
8 |
container_issue |
4 |
container_start_page |
1197 |
op_container_end_page |
1220 |
_version_ |
1772815031881170944 |