Theoretical results on a block preconditioner used in ice-sheet modeling: eigenvalue bounds for singular power-law fluids

• Main Authors: Helanow, Christian, Ahlkrona, Josefin
• Format: Text
• Language: unknown
• Published: 2023
• Subjects: Mathematics - Numerical Analysis; Physics - Computational Physics; 86A40; 76A05; 65F08; 86-08; Ice Sheet
• Online Access: http://arxiv.org/abs/2307.14688

The properties of a block preconditioner that has been successfully used in finite element simulations of large scale ice-sheet flow is examined. The type of preconditioner, based on approximating the Schur complement with the mass matrix scaled by the variable viscosity, is well-known in the context of Stokes flow and has previously been analyzed for other types of non-Newtonian fluids. We adapt the theory to hold for the regularized constitutive (power-law) equation for ice and derive eigenvalue bounds of the preconditioned system for both Picard and Newton linearization using \emph{inf-sup} stable finite elements. The eigenvalue bounds show that viscosity-scaled preconditioning clusters the eigenvalues well with only a weak dependence on the regularization parameter, while the eigenvalue bounds for the traditional non-viscosity-scaled mass-matrix preconditioner are very sensitive to the same regularization parameter. The results are verified numerically in two experiments using a manufactured solution with low regularity and a simulation of glacier flow. The numerical results further show that the computed eigenvalue bounds for the viscosity-scaled preconditioner are nearly independent of the regularization parameter. Experiments are performed using both Taylor-Hood and MINI elements, which are the common choices for \emph{inf-sup} stable elements in ice-sheet models. Both elements conform well to the theoretical eigenvalue bounds, with MINI elements being more sensitive to the quality of the meshes used in glacier simulations. Comment: 27 pages, 6 figures, 1 table