Preconditioning by Projectors in the Solution of Contact Problems: A Parallel Implementation |
| |
Authors: | Z. Dostál A. Friedlander F.A.M. Gomes S.A. Santos |
| |
Affiliation: | (1) Department of Applied Mathematics, VB-Technical University Ostrava, Ostrava, Czech Republic;(2) Department of Applied Mathematics, IMECC, University of Campinas, Campinas, Brazil |
| |
Abstract: | A non-overlapping domain decomposition algorithm of the Neumann–Neumann type for solving contact problems of elasticity is presented. Using the duality theory of convex programming, the discretized problem turns into a quadratic one with equality and bound constraints. The dual problem is modified by orthogonal projectors to the natural coarse space. The resulting problem is solved by an augmented Lagrangian algorithm. The projectors ensure an optimal convergence rate for the solution of the auxiliary linear problems by the preconditioned conjugate gradient method. Relevant aspects on the numerical linear algebra of these problems are presented, together with an efficient parallel implementation of the method. |
| |
Keywords: | contact problems parallel computation |
本文献已被 SpringerLink 等数据库收录! |
|