Hybrid domain decomposition algorithms for compressible and almost incompressible elasticity

Clark R. Dohrmann, Olof B. Widlund

Research output: Contribution to journalArticle

Abstract

Overlapping Schwarz methods are considered for mixed finite element approximations of linear elasticity, with discontinuous pressure spaces, as well as for compressible elasticity approximated by standard conforming finite elements. The coarse components of the preconditioners are based on spaces, with a number of degrees of freedom per subdomain which are uniformly bounded, which are similar to those previously developed for scalar elliptic problems and domain decomposition methods of iterative substructuring type, i.e. methods based on nonoverlapping decompositions of the domain. The local components of the new preconditioners are based on solvers on a set of overlapping subdomains. In the current study, the dimension of the coarse spaces is smaller than in recently developed algorithms; in the compressible case all independent face degrees of freedom have been eliminated while in the almost incompressible case five out of six are not needed. In many cases, this will result in a reduction of the dimension of the coarse space by about one half compared with that of the algorithm previously considered. In addition, in spite of using overlapping subdomains to define the local components of the preconditioner, values of the residual and the approximate solution need only to be retained on the interface between the subdomains in the iteration of the new hybrid Schwarz algorithm. The use of discontinuous pressures makes it possible to work exclusively with symmetric, positive-definite problems and the standard preconditioned conjugate gradient method. Bounds are established for the condition number of the preconditioned operators. The bound for the almost incompressible case grows in proportion to the square of the logarithm of the number of degrees of freedom of individual subdomains and the third power of the relative overlap between the overlapping subdomains, and it is independent of the Poisson ratio as well as jumps in the Lamé parameters across the interface between the subdomains. Numerical results illustrate the findings.

Original languageEnglish (US)
Pages (from-to)157-183
Number of pages27
JournalInternational Journal for Numerical Methods in Engineering
Volume82
Issue number2
DOIs
StatePublished - Apr 9 2010

Fingerprint

Decomposition Algorithm
Domain Decomposition
Overlapping
Elasticity
Preconditioner
Decomposition
Degree of freedom
Domain decomposition methods
Conjugate gradient method
Iterative Substructuring
Poisson ratio
Schwarz Methods
Problem Decomposition
Preconditioned Conjugate Gradient Method
Poisson's Ratio
Mixed Finite Elements
Linear Elasticity
Domain Decomposition Method
Condition number
Finite Element Approximation

Keywords

  • Almost incompressible elasticity
  • Domain decomposition
  • Iterative methods
  • Mixed finite element methods
  • Overlapping schwarz
  • Preconditioners

ASJC Scopus subject areas

  • Engineering(all)
  • Applied Mathematics
  • Numerical Analysis

Cite this

Hybrid domain decomposition algorithms for compressible and almost incompressible elasticity. / Dohrmann, Clark R.; Widlund, Olof B.

In: International Journal for Numerical Methods in Engineering, Vol. 82, No. 2, 09.04.2010, p. 157-183.

Research output: Contribution to journalArticle

@article{90d262bfb50249ffabe355183f5a2644,
title = "Hybrid domain decomposition algorithms for compressible and almost incompressible elasticity",
abstract = "Overlapping Schwarz methods are considered for mixed finite element approximations of linear elasticity, with discontinuous pressure spaces, as well as for compressible elasticity approximated by standard conforming finite elements. The coarse components of the preconditioners are based on spaces, with a number of degrees of freedom per subdomain which are uniformly bounded, which are similar to those previously developed for scalar elliptic problems and domain decomposition methods of iterative substructuring type, i.e. methods based on nonoverlapping decompositions of the domain. The local components of the new preconditioners are based on solvers on a set of overlapping subdomains. In the current study, the dimension of the coarse spaces is smaller than in recently developed algorithms; in the compressible case all independent face degrees of freedom have been eliminated while in the almost incompressible case five out of six are not needed. In many cases, this will result in a reduction of the dimension of the coarse space by about one half compared with that of the algorithm previously considered. In addition, in spite of using overlapping subdomains to define the local components of the preconditioner, values of the residual and the approximate solution need only to be retained on the interface between the subdomains in the iteration of the new hybrid Schwarz algorithm. The use of discontinuous pressures makes it possible to work exclusively with symmetric, positive-definite problems and the standard preconditioned conjugate gradient method. Bounds are established for the condition number of the preconditioned operators. The bound for the almost incompressible case grows in proportion to the square of the logarithm of the number of degrees of freedom of individual subdomains and the third power of the relative overlap between the overlapping subdomains, and it is independent of the Poisson ratio as well as jumps in the Lam{\'e} parameters across the interface between the subdomains. Numerical results illustrate the findings.",
keywords = "Almost incompressible elasticity, Domain decomposition, Iterative methods, Mixed finite element methods, Overlapping schwarz, Preconditioners",
author = "Dohrmann, {Clark R.} and Widlund, {Olof B.}",
year = "2010",
month = "4",
day = "9",
doi = "10.1002/nme.2761",
language = "English (US)",
volume = "82",
pages = "157--183",
journal = "International Journal for Numerical Methods in Engineering",
issn = "0029-5981",
publisher = "John Wiley and Sons Ltd",
number = "2",

}

TY - JOUR

T1 - Hybrid domain decomposition algorithms for compressible and almost incompressible elasticity

AU - Dohrmann, Clark R.

AU - Widlund, Olof B.

PY - 2010/4/9

Y1 - 2010/4/9

N2 - Overlapping Schwarz methods are considered for mixed finite element approximations of linear elasticity, with discontinuous pressure spaces, as well as for compressible elasticity approximated by standard conforming finite elements. The coarse components of the preconditioners are based on spaces, with a number of degrees of freedom per subdomain which are uniformly bounded, which are similar to those previously developed for scalar elliptic problems and domain decomposition methods of iterative substructuring type, i.e. methods based on nonoverlapping decompositions of the domain. The local components of the new preconditioners are based on solvers on a set of overlapping subdomains. In the current study, the dimension of the coarse spaces is smaller than in recently developed algorithms; in the compressible case all independent face degrees of freedom have been eliminated while in the almost incompressible case five out of six are not needed. In many cases, this will result in a reduction of the dimension of the coarse space by about one half compared with that of the algorithm previously considered. In addition, in spite of using overlapping subdomains to define the local components of the preconditioner, values of the residual and the approximate solution need only to be retained on the interface between the subdomains in the iteration of the new hybrid Schwarz algorithm. The use of discontinuous pressures makes it possible to work exclusively with symmetric, positive-definite problems and the standard preconditioned conjugate gradient method. Bounds are established for the condition number of the preconditioned operators. The bound for the almost incompressible case grows in proportion to the square of the logarithm of the number of degrees of freedom of individual subdomains and the third power of the relative overlap between the overlapping subdomains, and it is independent of the Poisson ratio as well as jumps in the Lamé parameters across the interface between the subdomains. Numerical results illustrate the findings.

AB - Overlapping Schwarz methods are considered for mixed finite element approximations of linear elasticity, with discontinuous pressure spaces, as well as for compressible elasticity approximated by standard conforming finite elements. The coarse components of the preconditioners are based on spaces, with a number of degrees of freedom per subdomain which are uniformly bounded, which are similar to those previously developed for scalar elliptic problems and domain decomposition methods of iterative substructuring type, i.e. methods based on nonoverlapping decompositions of the domain. The local components of the new preconditioners are based on solvers on a set of overlapping subdomains. In the current study, the dimension of the coarse spaces is smaller than in recently developed algorithms; in the compressible case all independent face degrees of freedom have been eliminated while in the almost incompressible case five out of six are not needed. In many cases, this will result in a reduction of the dimension of the coarse space by about one half compared with that of the algorithm previously considered. In addition, in spite of using overlapping subdomains to define the local components of the preconditioner, values of the residual and the approximate solution need only to be retained on the interface between the subdomains in the iteration of the new hybrid Schwarz algorithm. The use of discontinuous pressures makes it possible to work exclusively with symmetric, positive-definite problems and the standard preconditioned conjugate gradient method. Bounds are established for the condition number of the preconditioned operators. The bound for the almost incompressible case grows in proportion to the square of the logarithm of the number of degrees of freedom of individual subdomains and the third power of the relative overlap between the overlapping subdomains, and it is independent of the Poisson ratio as well as jumps in the Lamé parameters across the interface between the subdomains. Numerical results illustrate the findings.

KW - Almost incompressible elasticity

KW - Domain decomposition

KW - Iterative methods

KW - Mixed finite element methods

KW - Overlapping schwarz

KW - Preconditioners

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

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

U2 - 10.1002/nme.2761

DO - 10.1002/nme.2761

M3 - Article

VL - 82

SP - 157

EP - 183

JO - International Journal for Numerical Methods in Engineering

JF - International Journal for Numerical Methods in Engineering

SN - 0029-5981

IS - 2

ER -