A Parallel High-Order Accurate Finite Element Nonlinear Stokes Ice-Sheet Model and Benchmark Experiments

Applied and Computational Mathematics Seminar
Monday, April 4, 2011 - 14:00
1 hour (actually 50 minutes)
Skiles 005
Department of Mathematics, University of South Carolina
In this talk, we present a parallel finite element implementation ontetrahedral  grids of the nonlinear three-dimensional nonlinear Stokes model for thedynamics and evolution of ice-sheets. Discretization is based on a high-orderaccurate  scheme using the Taylor-Hood element pair. Both no-slip and sliding boundary conditions at the ice-bedrock boundary are studied. In addition, effective solvers using preconditioning techniques for the saddle-point system resulting fromthe  discretization are discussed and implemented. We demonstrate throughestablished ice-sheet benchmark experiments that our finite element nonlinear Stokesmodel  performs at least as well as other published and established Stokes modelsin the  field, and the parallel solver is shown to be efficient, robust, and scalable.