Source author record

Andrew J. Wathen

Andrew J. Wathen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Computational Engineering, Finance, and Science physics.geo-ph Numerical Analysis

Catalog footprint

What is connected

3works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2015arXiv

Three-field block-preconditioners for models of coupled magma/mantle dynamics

For a prescribed porosity, the coupled magma/mantle flow equations can be formulated as a two-field system of equations with velocity and pressure as unknowns. Previous work has shown that while optimal preconditioners for the two-field formulation can be obtained, the construction of preconditioners that are uniform with respect to model parameters is difficult. This limits the applicability of two-field preconditioners in certain regimes of practical interest. We address this issue by reformulating the governing equations as a three-field problem, which removes a term that was problematic in the two-field formulation in favour of an additional equation for a pressure-like field. For the three-field problem, we develop and analyse new preconditioners and we show numerically that they are optimal in terms of problem size and less sensitive to model parameters, compared to the two-field preconditioner. This extends the applicability of optimal preconditioners for coupled mantle/magma dynamics into parameter regimes of physical interest.

preprint2014arXiv

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

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.

preprint2010arXiv

Stability and preconditioning for a hybrid approximation on the sphere

This paper proposes a new preconditioning scheme for a linear system with a saddle-point structure arising from a hybrid approximation scheme on the sphere, an approximation scheme that combines (local) spherical radial basis functions and (global) spherical polynomials. Making use of a recently derived inf-sup condition [13] and the Brezzi stability and convergence theorem for this approximation scheme, we show that the linear system can be optimally preconditioned with a suitable block-diagonal preconditioner. Numerical experiments with a non-uniform distribution of data points support the theoretical conclusions.