BLOCK FETI-DP/BDDC PRECONDITIONERS FOR MIXED ISOGEOMETRIC DISCRETIZATIONS OF THREE-DIMENSIONAL ALMOST INCOMPRESSIBLE ELASTICITY

A block FETI-DP/BDDC preconditioner for mixed formulations of almost incompressible elasticity is constructed and analyzed; FETI-DP (Finite Element Tearing and Interconnecting Dual-Primal) and BDDC (Balancing Domain Decomposition by Constraints) are two very successful domain decomposition algorithms for a variety of elliptic problems. The saddle point problems of the mixed problems are discretized with mixed isogeometric analysis with continuous pressure fields. As in previous work by Tu and Li (2015) for finite element discretizations of the incompressible Stokes system, the proposed preconditioner is applied to a reduced positive definite system involving only the pressure interface variable and the Lagrange multiplier of the FETI-DP algorithm. In this work, we extend the theory to a wider class of saddle point problems and we propose a novel block-preconditioning strategy, which consists in using BDDC with deluxe scaling for the interface pressure block as well as deluxe scaling for the FETI-DP preconditioner for the Lagrange multiplier block. A convergence rate analysis is presented with a condition number bound for the preconditioned operator which depends on the inf-sup parameter of the fully assembled problem and the condition number of a closely related BDDC algorithm for compressible elasticity. This bound is scalable in the number of subdomains, poly-logarithmic in the ratio of subdomain and element sizes, and robust with respect to material incompressibility. Parallel numerical experiments validate the theory, demonstrate robustness in the presence of discontinuities of the Lame parameters, and indicate how the rate of convergence varies with respect to the spline polynomial degree and regularity and the deformation of the domain.

Automated summary

A block FETI-DP/BDDC preconditioner is constructed and analyzed for mixed formulations of almost incompressible elasticity. The strategy combines mixed isogeometric analysis with continuous pressure fields and deluxe scaling algorithms. Parallel numerical experiments validate the theory and demonstrate robustness in various scenarios.