Operator-Splitting Methods Respecting Eigenvalue Problems For Shallow Shelf Equations With Basal Drag

Jürgen Geiser , Reinhard Calov, Thomas Recknagel

Preprint series: Institut für Mathematik, Humboldt-Universität zu Berlin (ISSN 0863-0976), 31 pp.

MSC 2000

80A20 Heat and mass transfer, heat flow

80M25 Other numerical methods

Abstract
We discuss different numerical methods for solving the shallow shelf equations with basal drag. The coupled equations are decomposed into operators for membranes stresses, basal shear stress and driving stress. Applying reasonable parameter values, we demonstrate that the operator of the membrane stresses is much stiffer than operator of the basal shear stress. Therefore, we propose a new splitting method, which alternates between the iteration on the membrane-stress operator and the basal-shear operator, with a stronger iteration on the operator of the membrane stress. We show that this splitting improves the computational performance of the numerical method, although the choice of the (standard) method to solve for all operators in one step speeds up the scheme too. (Based on the delicate and coupled equation we propose a new decomposition method to decouple into simpler solvable sub-equations. After a number of approximations we consider the error of the method and proposed a choice of the operators.)