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

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