Robust BDDC algorithms for finite volume element methods

Abstract

The balancing domain decomposition by constraints (BDDC) method is applied to the linear system arising from the finite volume element method (FVEM) discretization of a scalar elliptic equation. The FVEMs share nice features of both finite element and finite volume methods and are flexible for complicated geometries with good conservation properties. However, the resulting linear system usually is asymmetric. The generalized minimal residual (GMRES) method is used to accelerate convergence. The proposed BDDC methods allow for jumps of the coefficient across subdomain interfaces. When jumps of the coefficient appear inside subdomains, the BDDC algorithms adaptively choose the primal variables deriving from the eigenvectors of some local generalized eigenvalue problems. The adaptive BDDC algorithms with advanced deluxe scaling can ensure good performance with highly discontinuous coefficients. A convergence analysis of the BDDC method with a preconditioned GMRES iteration is provided, and several numerical experiments confirm the theoretical estimate.