We will describe some recent work in the numerical simulation of problems of Geodynamics. The relevant partial differential equations share some of the features of the well-known Stokes equations, but there are significant differences. Our work has been to create rapid solvers for the large systems of equations arising from finite element approximation. We will briefly describe the relevant models and our preconditioned Krylov subspace iterative solvers which enable some of the first computations on these models.
This is joint work with Sander Rhebergen, Richard Katz, Garth Wells, John Rudge and Laura Alisic.
Additional Information
Location: ESB 4133 Andrew Wathen, Mathematics Department, University of Oxford, UK