Block BDDC/FETI-DP Preconditioners for Three-Field Mixed Finite Element Discretizations of Biot’s Consolidation Model

Abstract. In this paper, we construct and analyze a block dual-primal preconditioner for Biot’s consolidation model approximated by three-field mixed finite elements based on a displacement, pressure, and total pressure formulation. The domain is decomposed into nonoverlapping subdomains, and the continuity of the displacement component across the subdomain interface is enforced by introducing a Lagrange multiplier. After eliminating all displacement variables and the independent subdomain interior components of pressure and total pressure, the problem is reduced to a symmetric positive definite linear system for the subdomain interface pressure, total pressure, and the Lagrange multiplier. This reduced system is solved by a preconditioned conjugate gradient method, with a block dual-primal preconditioner using a balancing domain decomposition by constraints preconditioner for both the interface total pressure block and the interface pressure blocks, as well as a finite element tearing and interconnecting–dual primal preconditioner for the Lagrange multiplier block. By analyzing the conditioning of the preconditioned subsystem associated with the interface pressure and total pressure components, we obtain a condition number bound of the preconditioned system, which is scalable in the number of subdomains, poly-logarithmic in the ratio of subdomain and mesh sizes, and robust with respect to the parameters of the model. Extensive numerical experiments confirm the theoretical result of the proposed algorithm.