Parallel application of a novel domain decomposition preconditioner for the adaptive finite-element solution of three-dimensional convection-dominated PDEs

We describe and analyse the parallel implementation of a novel domain decomposition preconditioner for the fast iterative solution of linear systems of algebraic equations arising from the discretization of elliptic partial differential equations (PDEs) in three dimensions. In previous theoretical work, this preconditioner has been proved to be optimal for symmetric positive-definite (SPD) linear systems.

In this paper, we provide details of our three-dimensional parallel implementation and demonstrate that the technique may be generalized to the solution of non-symmetric algebraic systems, such as those arising when convection-diffusion problems are discretized using either Galerkin or stabilized finite-element methods (FEMs). Furthermore, we illustrate the potential of the preconditioner for use within an adaptive finite-element framework by successfully solving convection-dominated problems on locally, rather than globally, refined meshes.