Applied and Computational Mathematics Seminar
Monday, April 4, 2011 - 14:00
1 hour (actually 50 minutes)
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.