A matrix dependent/algebraic multigrid approach for extruded meshes with applications to ice sheet modeling
Abstract
A multigrid method is proposed that combines ideas from matrix dependent multigrid for structured grids and algebraic multigrid for unstructured grids. It targets problems where a three-dimensional mesh can be viewed as an extrusion of a two-dimensional, unstructured mesh in a third dimension. Our motivation comes from the modeling of thin structures via finite elements and, more specifically, the modeling of ice sheets. Extruded meshes are relatively common for thin structures and often give rise to anisotropic problems when the thin direction mesh spacing is much smaller than the broad direction mesh spacing. Within our approach, the first few multigrid hierarchy levels are obtained by applying matrix dependent multigrid to semicoarsen in a structured thin direction fashion. After sufficient structured coarsening, the resulting mesh contains only a single layer corresponding to a two-dimensional, unstructured mesh. Algebraic multigrid can then be employed in a standard manner to create further coarse levels, as the anisotropic phenomena is no longer present in the single layer problem. The overall approach remains fully algebraic, with the minor exception that some additional information is needed to determine the extruded direction. Furthermore, this facilitates integration of the solver with a variety of different extruded mesh applications.
- Authors:
-
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sandia National Lab. (SNL-CA), Livermore, CA (United States)
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Publication Date:
- Research Org.:
- Sandia National Lab. (SNL-CA), Livermore, CA (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Biological and Environmental Research (BER)
- OSTI Identifier:
- 1334163
- Report Number(s):
- SAND-2015-8309J
Journal ID: ISSN 1064-8275; 606326
- Grant/Contract Number:
- AC04-94AL85000
- Resource Type:
- Accepted Manuscript
- Journal Name:
- SIAM Journal on Scientific Computing
- Additional Journal Information:
- Journal Volume: 38; Journal Issue: 5; Journal ID: ISSN 1064-8275
- Publisher:
- SIAM
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; ice sheets; iterative solvers; algebraic method
Citation Formats
Tuminaro, Raymond S., Perego, Mauro, Tezaur, Irina Kalashnikova, Salinger, Andrew G., and Price, Stephen. A matrix dependent/algebraic multigrid approach for extruded meshes with applications to ice sheet modeling. United States: N. p., 2016.
Web. doi:10.1137/15m1040839.
Tuminaro, Raymond S., Perego, Mauro, Tezaur, Irina Kalashnikova, Salinger, Andrew G., & Price, Stephen. A matrix dependent/algebraic multigrid approach for extruded meshes with applications to ice sheet modeling. United States. https://doi.org/10.1137/15m1040839
Tuminaro, Raymond S., Perego, Mauro, Tezaur, Irina Kalashnikova, Salinger, Andrew G., and Price, Stephen. Thu .
"A matrix dependent/algebraic multigrid approach for extruded meshes with applications to ice sheet modeling". United States. https://doi.org/10.1137/15m1040839. https://www.osti.gov/servlets/purl/1334163.
@article{osti_1334163,
title = {A matrix dependent/algebraic multigrid approach for extruded meshes with applications to ice sheet modeling},
author = {Tuminaro, Raymond S. and Perego, Mauro and Tezaur, Irina Kalashnikova and Salinger, Andrew G. and Price, Stephen},
abstractNote = {A multigrid method is proposed that combines ideas from matrix dependent multigrid for structured grids and algebraic multigrid for unstructured grids. It targets problems where a three-dimensional mesh can be viewed as an extrusion of a two-dimensional, unstructured mesh in a third dimension. Our motivation comes from the modeling of thin structures via finite elements and, more specifically, the modeling of ice sheets. Extruded meshes are relatively common for thin structures and often give rise to anisotropic problems when the thin direction mesh spacing is much smaller than the broad direction mesh spacing. Within our approach, the first few multigrid hierarchy levels are obtained by applying matrix dependent multigrid to semicoarsen in a structured thin direction fashion. After sufficient structured coarsening, the resulting mesh contains only a single layer corresponding to a two-dimensional, unstructured mesh. Algebraic multigrid can then be employed in a standard manner to create further coarse levels, as the anisotropic phenomena is no longer present in the single layer problem. The overall approach remains fully algebraic, with the minor exception that some additional information is needed to determine the extruded direction. Furthermore, this facilitates integration of the solver with a variety of different extruded mesh applications.},
doi = {10.1137/15m1040839},
journal = {SIAM Journal on Scientific Computing},
number = 5,
volume = 38,
place = {United States},
year = {Thu Oct 06 00:00:00 EDT 2016},
month = {Thu Oct 06 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
The Multi-Grid Method for the Diffusion Equation with Strongly Discontinuous Coefficients
journal, December 1981
- Alcouffe, R. E.; Brandt, Achi; Dendy, Jr., J. E.
- SIAM Journal on Scientific and Statistical Computing, Vol. 2, Issue 4
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
Convergence estimates for multigrid algorithms without regularity assumptions
journal, September 1991
- Bramble, James H.; Pasciak, Joseph E.; Wang, Jun Ping
- Mathematics of Computation, Vol. 57, Issue 195
An energy-based AMG coarsening strategy
journal, January 2006
- Brannick, J.; Brezina, M.; MacLachlan, S.
- Numerical Linear Algebra with Applications, Vol. 13, Issue 2-3
Data assimilation and prognostic whole ice sheet modelling with the variationally derived, higher order, open source, and fully parallel ice sheet model VarGlaS
journal, January 2013
- Brinkerhoff, D. J.; Johnson, J. V.
- The Cryosphere, Vol. 7, Issue 4
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
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
Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate
journal, January 2015
- Cornford, S. L.; Martin, D. F.; Payne, A. J.
- The Cryosphere, Vol. 9, Issue 4
Black box multigrid for nonsymmetric problems
journal, January 1983
- Dendy, J. E.
- Applied Mathematics and Computation, Vol. 13, Issue 3-4
A Semicoarsening Multigrid Algorithm for SIMD Machines
journal, November 1992
- Dendy, Jr., J. E.; Ida, M. P.; Rutledge, J. M.
- SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 6
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
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
Data assimilation using a hybrid ice flow model
journal, January 2011
- Goldberg, D. N.; Sergienko, O. V.
- The Cryosphere, Vol. 5, Issue 2
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
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet
journal, September 2015
- Isaac, Tobin; Petra, Noemi; Stadler, Georg
- Journal of Computational Physics, Vol. 296
Solution of Nonlinear Stokes Equations Discretized By High-Order Finite Elements on Nonconforming and Anisotropic Meshes, with Application to Ice Sheet Dynamics
journal, January 2015
- Isaac, Tobin; Stadler, Georg; Ghattas, Omar
- SIAM Journal on Scientific Computing, Vol. 37, Issue 6
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
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
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
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
An ice-shelf model test based on the Ross Ice Shelf, Antarctica
journal, January 1996
- MacAyeal, D. R.; Rommelaere, V.; Huybrechts, P.
- Annals of Glaciology, Vol. 23
A new three-dimensional higher-order thermomechanical ice sheet model: Basic sensitivity, ice stream development, and ice flow across subglacial lakes
journal, January 2003
- Pattyn, Frank
- Journal of Geophysical Research, Vol. 108, Issue B8
Antarctic subglacial conditions inferred from a hybrid ice sheet/ice stream model
journal, July 2010
- Pattyn, Frank
- Earth and Planetary Science Letters, Vol. 295, Issue 3-4
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
Modelling West Antarctic ice sheet growth and collapse through the past five million years
journal, March 2009
- Pollard, David; DeConto, Robert M.
- Nature, Vol. 458, Issue 7236
A Semicoarsening Multigrid Method for Elliptic Partial Differential Equations with Highly Discontinuous and Anisotropic Coefficients
journal, January 1998
- Schaffer, Steve
- SIAM Journal on Scientific Computing, Vol. 20, Issue 1
Coulomb Friction and Other Sliding laws in a Higher-Order Glacier flow Model
journal, January 2010
- Schoof, Christian
- Mathematical Models and Methods in Applied Sciences, Vol. 20, Issue 01
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
A new perspective on strength measures in algebraic multigrid: STRENGTH MEASURES IN ALGEBRAIC MULTIGRID
journal, November 2009
- Olson, Luke N.; Schroder, Jacob; Tuminaro, Raymond S.
- Numerical Linear Algebra with Applications, Vol. 17, Issue 4
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
Albany/FELIX : a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis
journal, January 2015
- Tezaur, I. K.; Perego, M.; Salinger, A. G.
- Geoscientific Model Development, Vol. 8, Issue 4
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
journal, January 2015
- Tezaur, Irina K.; Tuminaro, Raymond S.; Perego, Mauro
- Procedia Computer Science, Vol. 51
Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems
journal, September 1996
- Vaněk, P.; Mandel, J.; Brezina, M.
- Computing, Vol. 56, Issue 3
Matrix-dependent prolongations and restrictions in a blackbox multigrid solver
journal, December 1990
- De Zeeuw, P. M.
- Journal of Computational and Applied Mathematics, Vol. 33, Issue 1
Works referencing / citing this record:
Exploring basal sliding with a fluidity-based, ice-sheet model using FOSLS: Fluidity-Based Ice-Sheet Models
journal, March 2018
- Allen, Jeffery; Manteuffel, Tom; Rajaram, Harihar
- Numerical Linear Algebra with Applications, Vol. 25, Issue 3
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