FSAI-based parallel Mixed Constraint Preconditioners for saddle point problems arising in geomechanics

  • Authors:
  • Luca Bergamaschi;ÁNgeles MartíNez

  • Affiliations:
  • -;-

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2011

Quantified Score

Hi-index 7.30

Visualization

Abstract

In this paper we propose and describe a parallel implementation of a block preconditioner for the solution of saddle point linear systems arising from Finite Element (FE) discretization of 3D coupled consolidation problems. The Mixed Constraint Preconditioner developed in [L. Bergamaschi, M. Ferronato, G. Gambolati, Mixed constraint preconditioners for the solution to FE coupled consolidation equations, J. Comput. Phys., 227(23) (2008), 9885-9897] is combined with the parallel FSAI preconditioner which is used here to approximate the inverses of both the structural (1, 1) block and an appropriate Schur complement matrix. The resulting preconditioner proves effective in the acceleration of the BiCGSTAB iterative solver. Numerical results on a number of test cases of size up to 2x10^6 unknowns and 1.2x10^8 nonzeros show the perfect scalability of the overall code up to 256 processors.