TY - JOUR
T1 - Coupled Models and Parallel Simulations for Three-Dimensional Full-Stokes Ice Sheet Modeling
AU - Huai Zhang, Lili Ju, Max Gunzburger, Todd Ringler & Stephen Price
JO - Numerical Mathematics: Theory, Methods and Applications
VL - 3
SP - 396
EP - 418
PY - 2011
DA - 2011/04
SN - 4
DO - http://doi.org/10.4208/nmtma.2011.m1031
UR - https://global-sci.org/intro/article_detail/nmtma/5975.html
KW - Ice sheet modeling, nonlinear Stokes equation, finite element method, parallel implementation, centroial Voronoi Delaunay meshes.
AB - <p style="text-align: justify;">A three-dimensional full-Stokes computational model is considered
for determining the dynamics, temperature, and thickness of ice
sheets. The governing thermo-mechanical equations consist of the
three-dimensional full-Stokes system with nonlinear rheology for the
momentum, an advective-diffusion energy equation for temperature
evolution, and a mass conservation equation for ice-thickness
changes. Here, we discuss the variable resolution meshes, the finite
element discretizations, and the parallel algorithms employed by the
model components. The solvers are integrated through a well-designed
coupler for the exchange of parametric data between components. The
discretization utilizes high-quality, variable-resolution centroidal
Voronoi Delaunay triangulation meshing and existing parallel
solvers. We demonstrate the gridding technology, discretization
schemes, and the efficiency and scalability of the parallel solvers
through computational experiments using both simplified geometries
arising from benchmark test problems and a realistic Greenland ice
sheet geometry.</p>