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

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.

[1]  Ray S. Tuminaro,et al.  Parallel Smoothed Aggregation Multigrid : Aggregation Strategies on Massively Parallel Machines , 2000, ACM/IEEE SC 2000 Conference (SC'00).

[2]  George Karypis,et al.  Parmetis parallel graph partitioning and sparse matrix ordering library , 1997 .

[3]  Ümit V. Çatalyürek,et al.  Hypergraph-based Dynamic Load Balancing for Adaptive Scientific Computations , 2007, 2007 IEEE International Parallel and Distributed Processing Symposium.

[4]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[5]  Hans Petter Langtangen,et al.  Efficient block preconditioners for the coupled equations of pressure and deformation in highly discontinuous media , 2011 .

[6]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[7]  R. Tuminaro,et al.  A parallel block multi-level preconditioner for the 3D incompressible Navier--Stokes equations , 2003 .

[8]  R. Fletcher.  Conjugate gradient methods for indefinite systems , 1976 .

[9]  Padma Raghavan,et al.  Parallel Processing for Scientific Computing , 2006, Software, Environments, Tools.

[10]  Mark F. Adams,et al.  Ultrascalable Implicit Finite Element Analyses in Solid Mechanics with over a Half a Billion Degrees of Freedom , 2004, Proceedings of the ACM/IEEE SC2004 Conference.

[11]  W. Hackbusch.  Iterative Solution of Large Sparse Systems of Equations , 1993 .

[12]  Jonathan J. Hu,et al.  ML 5.0 Smoothed Aggregation Users's Guide , 2006 .

[13]  Rajeev Thakur,et al.  Improving the Performance of Collective Operations in MPICH , 2003, PVM/MPI.

[14]  Edmond Chow,et al.  A Survey of Parallelization Techniques for Multigrid Solvers , 2006, Parallel Processing for Scientific Computing.

[15]  Ulrike Meier Yang,et al.  Parallel Algebraic Multigrid Methods — High Performance Preconditioners , 2006 .

[16]  Thomas A. Manteuffel,et al.  Adaptive Smoothed Aggregation (αSA) , 2004, SIAM J. Sci. Comput..

[17]  John N. Shadid,et al.  A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier-Stokes equations , 2008, J. Comput. Phys..

[18]  Hans Petter Langtangen,et al.  Computational Partial Differential Equations - Numerical Methods and Diffpack Programming , 1999, Lecture Notes in Computational Science and Engineering.

[19]  Aslak Tveito,et al.  Numerical solution of partial differential equations on parallel computers , 2006 .

[20]  SCALABLE ALGEBRAIC MULTIGRID ON 3500 PROCESSORS , .

[21]  Joachim Berdal Haga,et al.  On the causes of pressure oscillations in low‐permeable and low‐compressible porous media , 2012 .

[22]  Jonathan J. Hu,et al.  ML 3.1 smoothed aggregation user's guide. , 2004 .

[23]  A. George,et al.  On the Complexity of Sparse $QR$ and $LU$ Factorization of Finite-Element Matrices , 1988 .

[24]  K. Toh,et al.  An efficient diagonal preconditioner for finite element solution of Biot's consolidation equations , 2002 .

[25]  K. Toh,et al.  Block preconditioners for symmetric indefinite linear systems , 2004 .

[26]  Mark F. Adams,et al.  Evaluation of three unstructured multigrid methods on 3D finite element problems in solid mechanics , 2000 .

[27]  Takumi Washio,et al.  Ordering strategies and related techniques to overcome the trade-off between parallelism and convergence in incomplete factorizations , 1999, Parallel Comput..