Scalable algorithms for three-field mixed finite element coupled poromechanics

JOURNAL OF COMPUTATIONAL PHYSICS

卷:327

Article

Castelletto, Nicola¹  White, Joshua A.²  Ferronatoc, Massimiliano³ 

[1] Stanford Univ, Energy Resources Engn, Stanford, CA 94305 USA

[2] Lawrence Livermore Natl Lab, Atmospher Earth & Energy Div, Livermore, CA USA

[3] Univ Padua, Dept Civil Environm & Architectural Engn, I-35100 Padua, Italy

关键词: Poromechanics;    Preconditioners;    Iterative methods;    Mixed formulation;    Algebraic multigrid;   

DOI  :  10.1016/j.jcp.2016.09.063

来源: Elsevier

【 摘 要 】

We introduce a class of block preconditioners for accelerating the iterative solution of coupled poromechanics equations based on a three-field formulation. The use of a displacement/velocity/pressure mixed finite-element method combined with a first order backward difference formula for the approximation of time derivatives produces a sequence of linear systems with a 3 x 3 unsymmetric and indefinite block matrix. The preconditioners are obtained by approximating the two-level Schur complement with the aid of physically-based arguments that can be also generalized in a purely algebraic approach. A theoretical and experimental analysis is presented that provides evidence of the robustness, efficiency and scalability of the proposed algorithm. The performance is also assessed for a real-world challenging consolidation experiment of a shallow formation. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free