A Scalable Auxiliary Space Preconditioner for High-Order Finite Element Methods

Authors Young-Ju Lee, Wei Leng, Chensong Zhang

Publication date 11 April 2012

Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55)

Abstract: In this paper, we revisit an auxiliary space preconditioning method proposed by Xu [Computing 56, 1996], in which low-order finite element spaces are employed as auxiliary spaces for solving linear algebraic systems arising from high-order finite element discretizations. We provide a new convergence rate estimate and parallel implementation of the proposed algorithm. We show that this method is user-friendly and can play an important role in a variety of Poisson-based solvers for more challenging problems such as the Navier--Stokes equation. We investigate the performance of the proposed algorithm using the Poisson equation and the Stokes equation on 3D unstructured grids. Numerical results demonstrate the advantages of the proposed algorithm in terms of efficiency, robustness, and parallel scalability.

This page was built for publication: A Scalable Auxiliary Space Preconditioner for High-Order Finite Element Methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6232274)