Enhanced balancing Neumann-Neumann preconditioning in computational fluid and solid mechanics

Research output: Contribution to journalArticleResearchpeer-review

12 Citations (Scopus)

Abstract

In this work, we propose an enhanced implementation of balancing Neumann-Neumann (BNN) preconditioning together with a detailed numerical comparison against the balancing domain decomposition by constraints (BDDC) preconditioner. As model problems, we consider the Poisson and linear elasticity problems. On one hand, we propose a novel way to deal with singular matrices and pseudo-inverses appearing in local solvers. It is based on a kernel identification strategy that allows us to efficiently compute the action of the pseudo-inverse via local indefinite solvers. We further show how, identifying a minimum set of degrees of freedom to be fixed, an equivalent definite system can be solved instead, even in the elastic case. On the other hand, we propose a simple implementation of the algorithm that reduces the number of Dirichlet solvers to only one per iteration, leading to similar computational cost as additive methods. After these improvements of the BNN preconditioned conjugate gradient algorithm, we compare its performance against that of the BDDC preconditioners on a pair of large-scale distributed-memory platforms. The enhanced BNN method is a competitive preconditioner for three-dimensional Poisson and elasticity problems and outperforms the BDDC method in many cases.

Original languageEnglish
Pages (from-to)203-230
Number of pages28
JournalInternational Journal for Numerical Methods in Engineering
Volume96
Issue number4
DOIs
Publication statusPublished - 24 Jul 2013
Externally publishedYes

Keywords

  • Balancing domain decomposition
  • BDDC
  • BNN
  • Coarse-grid correction
  • domain decomposition
  • elasticity
  • parallelization
  • scalability

Cite this