skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A matrix dependent/algebraic multigrid approach for extruded meshes with applications to ice sheet modeling

Journal Article · · SIAM Journal on Scientific Computing
DOI:https://doi.org/10.1137/15m1040839· OSTI ID:1334163
 [1];  [1];  [2];  [1];  [3]
  1. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
  2. Sandia National Lab. (SNL-CA), Livermore, CA (United States)
  3. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
+ Show Author Affiliations

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.

Research Organization:
Sandia National Lab. (SNL-CA), Livermore, CA (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Biological and Environmental Research (BER)
Grant/Contract Number:
AC04-94AL85000
OSTI ID:
1334163
Report Number(s):
SAND-2015-8309J; 606326
Journal Information:
SIAM Journal on Scientific Computing, Vol. 38, Issue 5; ISSN 1064-8275
Publisher:
SIAMCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 8 works
Citation information provided by
Web of Science

References (34)

The Multi-Grid Method for the Diffusion Equation with Strongly Discontinuous Coefficients journal December 1981
Velocity and stress fields in grounded glaciers: a simple algorithm for including deviatoric stress gradients journal January 1995
Convergence estimates for multigrid algorithms without regularity assumptions journal September 1991
An energy-based AMG coarsening strategy journal January 2006
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
Semicoarsening Multigrid on Distributed Memory Machines journal January 2000
Adaptive mesh, finite volume modeling of marine ice sheets journal January 2013
Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate journal January 2015
Black box multigrid for nonsymmetric problems journal January 1983
A Semicoarsening Multigrid Algorithm for SIMD Machines journal November 1992
Consistent approximations and boundary conditions for ice-sheet dynamics from a principle of least action journal January 2010
Capabilities and performance of Elmer/Ice, a new-generation ice sheet model journal January 2013
Data assimilation using a hybrid ice flow model journal January 2011
An overview of the Trilinos project journal September 2005
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
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
Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM): ICE SHEET SYSTEM MODEL journal March 2012
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
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
Toward a New Generation of Ice Sheet Models journal January 2007
An ice-shelf model test based on the Ross Ice Shelf, Antarctica journal January 1996
A new three-dimensional higher-order thermomechanical ice sheet model: Basic sensitivity, ice stream development, and ice flow across subglacial lakes journal January 2003
Antarctic subglacial conditions inferred from a hybrid ice sheet/ice stream model journal July 2010
Optimal initial conditions for coupling ice sheet models to Earth system models: PEREGO ET AL. journal September 2014
Modelling West Antarctic ice sheet growth and collapse through the past five million years journal March 2009
A Semicoarsening Multigrid Method for Elliptic Partial Differential Equations with Highly Discontinuous and Anisotropic Coefficients journal January 1998
Coulomb Friction and Other Sliding laws in a Higher-Order Glacier flow Model journal January 2010
Thin-Film Flows with Wall Slip: An Asymptotic Analysis of Higher Order Glacier Flow Models journal January 2010
A new perspective on strength measures in algebraic multigrid: STRENGTH MEASURES IN ALGEBRAIC MULTIGRID journal November 2009
Enhanced basal lubrication and the contribution of the Greenland ice sheet to future sea-level rise journal August 2013
Albany/FELIX : a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis journal January 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 journal January 2015
Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems journal September 1996
Matrix-dependent prolongations and restrictions in a blackbox multigrid solver journal December 1990

Cited By (3)

Exploring basal sliding with a fluidity-based, ice-sheet model using FOSLS: Fluidity-Based Ice-Sheet Models journal March 2018
MPAS-Albany Land Ice (MALI): a variable-resolution ice sheet model for Earth system modeling using Voronoi grids journal January 2018
A Study on the Performance Portability of the Finite Element Assembly Process Within the Albany Land Ice Solver book February 2020

Similar Records

Nonlinear multigrid based on local spectral coarsening for heterogeneous diffusion problems
Journal Article · Tue Dec 01 00:00:00 EST 2020 · Computer Methods in Applied Mechanics and Engineering · OSTI ID:1334163
+2 more
Comparison of multigrid algorithms for high‐order continuous finite element discretizations
Journal Article · Thu Apr 02 00:00:00 EDT 2015 · Numerical Linear Algebra with Applications · OSTI ID:1334163
Efficient Operator-Coarsening Multigrid Schemes for Local Discontinuous Galerkin Methods
Journal Article · Tue Dec 17 00:00:00 EST 2019 · SIAM Journal on Scientific Computing · OSTI ID:1334163

Related Subjects

97 MATHEMATICS AND COMPUTING
ice sheets
iterative solvers
algebraic method