A parallel block preconditioner for large-scale poroelasticity with highly heterogeneous material parameters

作者：Haga Joachim Berdal^*; Osnes Harald; Langtangen Hans Petter

来源：Computational Geosciences, 2012, 16(3): 723-734.

DOI：10.1007/s10596-012-9284-4

• Parallel computing
• Block preconditioning
• Poroelasticity
• Low-permeable media
• Finite elements

Large-scale simulations of coupled flow in deformable porous media require iterative methods for solving the systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in realistic geological applications, such as basin evolution at regional scales. The success of iterative methods for such problems depends strongly on finding effective preconditioners with good parallel scaling properties, which is the topic of the present paper. We present a parallel preconditioner for Biot's equations of coupled elasticity and fluid flow in porous media. The preconditioner is based on an approximation of the exact inverse of the two-by-two block system arising from a finite element discretisation. The approximation relies on a highly scalable approximation of the global Schur complement of the coefficient matrix, combined with generally available state-of-the-art multilevel preconditioners for the individual blocks. This preconditioner is shown to be robust on problems with highly heterogeneous material parameters. We investigate the weak and strong parallel scaling of this preconditioner on up to 512 processors and demonstrate its ability on a realistic basin-scale problem in poroelasticity with over eight million tetrahedral elements.

• 出版日期2012-6