Balancing Neumann-Neumann preconditioners for mixed approximations of heterogeneous problems in linear elasticity

Summary. Balancing Neumann-Neumann methods are extented to mixed formulations of the linear elasticity system with discontinuous coefficients, discretized with mixed finite or spectral elements with discontinuous pressures. These domain decomposition methods implicitly eliminate the degrees of freedom associated with the interior of each subdomain and solve iteratively the resulting saddle point Schur complement using a hybrid preconditioner based on a coarse mixed elasticity problem and local mixed elasticity problems with natural and essential boundary conditions. A polylogarithmic bound in the local number of degrees of freedom is proven for the condition number of the preconditioned operator in the constant coefficient case. Parallel and serial numerical experiments confirm the theoretical results, indicate that they still hold for systems with discontinuous coefficients, and show that our algorithm is scalable, parallel, and robust with respect to material heterogeneities. The results on heterogeneous general problems are also supported in part by our theory.Mathematics Subject Classification (1991):65N55, 65N30, 65N35, 65F10, 65Y05This work was supported by a scholarship of CNPq, of the Ministry for Science and Technology of Brazil, under process 201205/97-1. The work was developed in part at MCS/ANL-DOE, under a Givens Research Associate appointment in Summer 2001.This work was supported in part by the National Science Foundation under Grant NSF-CCR-9732208 and in part by MIUR.This work was supported in part by the National Science Foundation under Grants qNSF-CCR-9732208, and in part by the U.S. Department of Energy under contracts DE-FC02-01ER25482 and DE-FG02-92ER25127.