On the use of rigid body modes in the deflated preconditioned conjugate gradient method

T. B. Jönsthövel, M. B. Van Gijzen, C. Vuik, A. Scarpas

Research output: Contribution to journalArticlepeer-review

13 Scopus citations

Abstract

Large discontinuities in material properties, such as those encountered in composite materials, lead to ill-conditioned systems of linear equations. These discontinuities give rise to small eigenvalues that may negatively affect the convergence of iterative solution methods such as the preconditioned conjugate gradient method. This paper considers the deflated preconditioned conjugate gradient method for solving such systems. Our deflation technique uses as the deflation space the rigid body modes of sets of elements with homogeneous material properties. We show that in the deflated spectrum the small eigenvalues are mapped to zero and no longer negatively affect the convergence. We justify our approach through mathematical analysis and show with numerical experiments on both academic and realistic test problems that the convergence of our DPCG method is independent of discontinuities in the material properties.

Original languageBritish English
Pages (from-to)B207-B225
JournalSIAM Journal on Scientific Computing
Volume35
Issue number1
DOIs
StatePublished - 2013

Keywords

  • Conjugate gradients
  • CT scan
  • Deflation
  • Preconditioners
  • Rigid body modes
  • Structural mechanics

Fingerprint

Dive into the research topics of 'On the use of rigid body modes in the deflated preconditioned conjugate gradient method'. Together they form a unique fingerprint.

Cite this