A Matrix Dependent/Algebraic Multigrid Approach for Extruded Meshes with Applications to Ice Sheet Modeling

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 le...

[1]  Jinchao Xu,et al.  Convergence estimates for multigrid algorithms without regularity assumptions , 1991 .

[2]  Philippe Huybrechts,et al.  An ice-shelf model test based on the Ross Ice Shelf, Antarctica , 1996, Annals of Glaciology.

[3]  P. M. De Zeeuw,et al.  Matrix-dependent prolongations and restrictions in a blackbox multigrid solver , 1990 .

[4]  L. Morland.  Unconfined Ice-Shelf Flow , 1987 .

[5]  C. Schoof.  COULOMB FRICTION AND OTHER SLIDING LAWS IN A HIGHER-ORDER GLACIER FLOW MODEL , 2010 .

[6]  Douglas J. Brinkerhoff,et al.  Data assimilation and prognostic whole ice sheet modelling with the variationally derived, higher order, open source, and fully parallel ice sheet model VarGlaS , 2013 .

[7]  Georg Stadler,et al.  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 , 2014, J. Comput. Phys..

[8]  Frank Pattyn,et al.  Antarctic subglacial conditions inferred from a hybrid ice sheet/ice stream model , 2010 .

[9]  G. Stadler,et al.  Optimal initial conditions for coupling ice sheet models to Earth system models , 2014 .

[10]  H. Blatter.  Velocity and stress fields in grounded glaciers: a simple algorithm for including deviatoric stress gradients , 1995 .

[11]  Todd D. Ringler,et al.  A parallel high‐order accurate finite element nonlinear Stokes ice sheet model and benchmark experiments , 2012 .

[12]  D. Goldberg,et al.  Data assimilation using a hybrid ice flow model , 2010 .

[13]  Georg Stadler,et al.  Solution of Nonlinear Stokes Equations Discretized By High-Order Finite Elements on Nonconforming and Anisotropic Meshes, with Application to Ice Sheet Dynamics , 2014, SIAM J. Sci. Comput..

[14]  Jacob B. Schroder,et al.  A new perspective on strength measures in algebraic multigrid , 2010, Numer. Linear Algebra Appl..

[15]  Thomas A. Manteuffel,et al.  An energy‐based AMG coarsening strategy , 2006, Numer. Linear Algebra Appl..

[16]  Mika Malinen,et al.  Capabilities and performance of Elmer/Ice, a new-generation ice sheet model , 2013 .

[17]  William H. Lipscomb,et al.  A Community Ice Sheet Model for Sea Level Prediction , 2009 .

[18]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[19]  Yair Shapira.  Matrix-Based Multigrid: Theory and Applications , 2008 .

[20]  A. Brandt,et al.  The Multi-Grid Method for the Diffusion Equation with Strongly Discontinuous Coefficients , 1981 .

[21]  John Frederick Nye,et al.  The distribution of stress and velocity in glaciers and ice-sheets , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[22]  William H. Lipscomb,et al.  Toward a new generation of ice sheet models , 2007 .

[23]  Mauro Perego,et al.  Enhanced basal lubrication and the contribution of the Greenland ice sheet to future sea-level rise , 2013, Proceedings of the National Academy of Sciences.

[24]  Daniel F. Martin,et al.  Adaptive mesh, finite volume modeling of marine ice sheets , 2013, J. Comput. Phys..

[25]  Irina K. Tezaur,et al.  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 , 2015, ICCS.

[26]  David Pollard,et al.  Modelling West Antarctic ice sheet growth and collapse through the past five million years , 2009, Nature.

[27]  William H. Lipscomb,et al.  Consistent approximations and boundary conditions for ice-sheet dynamics from a principle of least action , 2010, Journal of Glaciology.

[28]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[29]  P. Wesseling.  An Introduction to Multigrid Methods , 1992 .

[30]  Ralf Hiptmair,et al.  Analysis of tensor product multigrid , 2001, Numerical Algorithms.

[31]  K. Hutter.  Theoretical Glaciology: Material Science of Ice and the Mechanics of Glaciers and Ice Sheets , 1983 .

[32]  Robert D. Falgout,et al.  Semicoarsening Multigrid on Distributed Memory Machines , 1999, SIAM J. Sci. Comput..

[33]  W. Paterson,et al.  The Physics of Glaciers Ed. 4 , 2010 .

[34]  Ed Bueler,et al.  Shallow shelf approximation as a “sliding law” in a thermomechanically coupled ice sheet model , 2008, 0810.3449.

[35]  M. P. Ida,et al.  A Semicoarsening Multigrid Algorithm for SIMD Machines , 1992, SIAM J. Sci. Comput..

[36]  J. W. Ruge,et al.  4. Algebraic Multigrid , 1987 .

[37]  J. E. Dendy.  Black box multigrid for nonsymmetric problems , 1983 .

[38]  Eric Rignot,et al.  Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM) , 2012 .

[39]  Daniel F. Martin,et al.  Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate , 2015 .

[40]  STEVE SCHAFFER,et al.  A Semicoarsening Multigrid Method for Elliptic Partial Differential Equations with Highly Discontinuous and Anisotropic Coefficients , 1998, SIAM J. Sci. Comput..

[41]  F. Pattyn.  A new three-dimensional higher-order thermomechanical ice sheet model: Basic sensitivity, ice stream development, and ice flow across subglacial lakes , 2003 .

[42]  Christian Schoof,et al.  Thin-Film Flows with Wall Slip: An Asymptotic Analysis of Higher Order Glacier Flow Models , 2010 .

[43]  Marian Brezina,et al.  Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems , 2005, Computing.

[44]  Yuri A. Melnikov,et al.  Influence Functions and Matrices , 1998 .