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

Applied and Computational Mathematics Seminar
Monday, April 4, 2011 - 2:00pm for 1 hour (actually 50 minutes)
Skiles 005
Lili Ju – Department of Mathematics, University of South Carolina – http://www.math.sc.edu/~ju/
Haomin Zhou
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.