Repository logo
 

Analysis of block preconditioners for models of coupled magma/mantle dynamics


Type

Article

Change log

Authors

Rhebergen, S 
Wells, GN 
Katz, RF 
Wathen, AJ 

Abstract

This article considers the iterative solution of a finite element discretisation of the magma dynamics equations. In simplified form, the magma dynamics equations share some features of the Stokes equations. We therefore formulate, analyse and numerically test a Elman, Silvester and Wathen-type block preconditioner for magma dynamics. We prove analytically and demonstrate numerically the optimality of the preconditioner. The presented analysis highlights the dependence of the preconditioner on parameters in the magma dynamics equations that can affect convergence of iterative linear solvers. The analysis is verified through a range of two- and three-dimensional numerical examples on unstructured grids, from simple illustrative problems through to large problems on subduction zone-like geometries. The computer code to reproduce all numerical examples is freely available as supporting material.

Description

Keywords

magma dynamics, mantle dynamics, finite element method, preconditioners

Journal Title

SIAM Journal on Scientific Computing

Conference Name

Journal ISSN

1064-8275
1095-7197

Volume Title

36

Publisher

Society for Industrial & Applied Mathematics (SIAM)
Sponsorship
Natural Environment Research Council (NE/I023929/1)