### Abstract

Low-order finite element discretizations of the linear elasticity system suffer increasingly from locking effects and ill-conditioning, when the material approaches the incompressible limit, if only the displacement variables are used. Mixed finite elements using both displacement and pressure variables provide a well-known remedy, but they yield larger and indefinite discrete systems for which the design of scalable and efficient iterative solvers is challenging. Two-level overlapping Schwarz preconditioners for the almost incompressible system of linear elasticity, discretized by mixed finite elements with discontinuous pressures, are constructed and analyzed. The preconditioned systems are accelerated either by a GMRES (generalized minimum residual) method applied to the resulting discrete saddle point problem or by a PCG (preconditioned conjugate gradient) method applied to a positive definite, although extremely ill-conditioned, reformulation of the problem obtained by eliminating all pressure variables on the element level. A novel theoretical analysis of the algorithm for the positive definite reformulation is given by extending some earlier results by Dohrmann and Widlund. The main result of the paper is a bound on the condition number of the algorithm which is cubic in the relative overlap and grows logarithmically with the number of elements across individual subdomains but is otherwise independent of the number of subdomains, their diameters and mesh sizes, the incompressibility of the material, and possible discontinuities of the material parameters across the subdomain interfaces. Numerical results in the plane confirm the theory and also indicate that an analogous result should hold for the saddle point formulation, as well as for spectral element discretizations.

Original language | English (US) |
---|---|

Pages (from-to) | A811-A830 |

Journal | SIAM Journal on Scientific Computing |

Volume | 37 |

Issue number | 2 |

DOIs | |

State | Published - 2015 |

### Fingerprint

### Keywords

- Almost incompressible linear elasticity
- Domain decomposition methods
- Mixed finite and spectral elements
- Saddle point problems
- Two-level overlapping Schwarz preconditioners

### ASJC Scopus subject areas

- Applied Mathematics
- Computational Mathematics

### Cite this

*SIAM Journal on Scientific Computing*,

*37*(2), A811-A830. https://doi.org/10.1137/140981861

**Overlapping schwarz methods with a standard coarse space for almost incompressible linear elasticity.** / Cai, Mingchao; Pavarino, Luca F.; Widlund, Olof B.

Research output: Contribution to journal › Article

*SIAM Journal on Scientific Computing*, vol. 37, no. 2, pp. A811-A830. https://doi.org/10.1137/140981861

}

TY - JOUR

T1 - Overlapping schwarz methods with a standard coarse space for almost incompressible linear elasticity

AU - Cai, Mingchao

AU - Pavarino, Luca F.

AU - Widlund, Olof B.

PY - 2015

Y1 - 2015

N2 - Low-order finite element discretizations of the linear elasticity system suffer increasingly from locking effects and ill-conditioning, when the material approaches the incompressible limit, if only the displacement variables are used. Mixed finite elements using both displacement and pressure variables provide a well-known remedy, but they yield larger and indefinite discrete systems for which the design of scalable and efficient iterative solvers is challenging. Two-level overlapping Schwarz preconditioners for the almost incompressible system of linear elasticity, discretized by mixed finite elements with discontinuous pressures, are constructed and analyzed. The preconditioned systems are accelerated either by a GMRES (generalized minimum residual) method applied to the resulting discrete saddle point problem or by a PCG (preconditioned conjugate gradient) method applied to a positive definite, although extremely ill-conditioned, reformulation of the problem obtained by eliminating all pressure variables on the element level. A novel theoretical analysis of the algorithm for the positive definite reformulation is given by extending some earlier results by Dohrmann and Widlund. The main result of the paper is a bound on the condition number of the algorithm which is cubic in the relative overlap and grows logarithmically with the number of elements across individual subdomains but is otherwise independent of the number of subdomains, their diameters and mesh sizes, the incompressibility of the material, and possible discontinuities of the material parameters across the subdomain interfaces. Numerical results in the plane confirm the theory and also indicate that an analogous result should hold for the saddle point formulation, as well as for spectral element discretizations.

AB - Low-order finite element discretizations of the linear elasticity system suffer increasingly from locking effects and ill-conditioning, when the material approaches the incompressible limit, if only the displacement variables are used. Mixed finite elements using both displacement and pressure variables provide a well-known remedy, but they yield larger and indefinite discrete systems for which the design of scalable and efficient iterative solvers is challenging. Two-level overlapping Schwarz preconditioners for the almost incompressible system of linear elasticity, discretized by mixed finite elements with discontinuous pressures, are constructed and analyzed. The preconditioned systems are accelerated either by a GMRES (generalized minimum residual) method applied to the resulting discrete saddle point problem or by a PCG (preconditioned conjugate gradient) method applied to a positive definite, although extremely ill-conditioned, reformulation of the problem obtained by eliminating all pressure variables on the element level. A novel theoretical analysis of the algorithm for the positive definite reformulation is given by extending some earlier results by Dohrmann and Widlund. The main result of the paper is a bound on the condition number of the algorithm which is cubic in the relative overlap and grows logarithmically with the number of elements across individual subdomains but is otherwise independent of the number of subdomains, their diameters and mesh sizes, the incompressibility of the material, and possible discontinuities of the material parameters across the subdomain interfaces. Numerical results in the plane confirm the theory and also indicate that an analogous result should hold for the saddle point formulation, as well as for spectral element discretizations.

KW - Almost incompressible linear elasticity

KW - Domain decomposition methods

KW - Mixed finite and spectral elements

KW - Saddle point problems

KW - Two-level overlapping Schwarz preconditioners

UR - http://www.scopus.com/inward/record.url?scp=84928974761&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84928974761&partnerID=8YFLogxK

U2 - 10.1137/140981861

DO - 10.1137/140981861

M3 - Article

VL - 37

SP - A811-A830

JO - SIAM Journal of Scientific Computing

JF - SIAM Journal of Scientific Computing

SN - 1064-8275

IS - 2

ER -