Analysis of FETI methods for multiscale PDEs

In this paper, we study a variant of the finite element tearing and interconnecting (FETI) method which is suitable for elliptic PDEs with highly heterogeneous (multiscale) coefficients α(x); in particular, coefficients with strong variation within subdomains and/or jumps that are not aligned with the subdomain interfaces. Using energy minimisation and cut-off arguments we can show rigorously that for an arbitrary (positive) coefficient function the condition number of the preconditioned FETI system can be bounded by C(α) (1 + log(H/h))² where H is the subdomain diameter and h is the mesh size, and where the function C(α) depends only on the coefficient variation in the vicinity of subdomain interfaces. In particular, if varies only mildly in a layer Ω _i,η of width η near the boundary of each of the subdomains Ω _i , then , independent of the variation of α in the remainder Ω _i \Ω _i,η of each subdomain and independent of any jumps of α across subdomain interfaces. The quadratic dependence of C(α) on H/η can be relaxed to a linear dependence under stronger assumptions on the behaviour of α in the interior of the subdomains. Our theoretical findings are confirmed in numerical tests. C. Pechstein was supported by the Austrian Science Funds (FWF) under grant F1306.