 |
| Titel |
Albany/FELIX: a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis |
| VerfasserIn |
I. K. Tezaur, M. Perego, A. G. Salinger, R. S. Tuminaro, S. F. Price |
| Medientyp |
Artikel
|
| Sprache |
Englisch
|
| ISSN |
1991-959X
|
| Digitales Dokument |
URL |
| Erschienen |
In: Geoscientific Model Development ; 8, no. 4 ; Nr. 8, no. 4 (2015-04-27), S.1197-1220 |
| Datensatznummer |
250116285
|
| Publikation (Nr.) |
 copernicus.org/gmd-8-1197-2015.pdf |
|
|
|
|
|
| Zusammenfassung |
| 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. |
| |
|
| Teil von |
|
|
|
|
|
|
|
|