A Relaxed Physical Factorization Preconditioner for Mixed Finite Element Coupled Poromechanics
In this work, we introduce a relaxed physical factorization (RPF) preconditioner for the efficient iterative solution of the linearized algebraic system arising from the mixed finite element discretization of coupled poromechanics equations. The preconditioner is obtained by using a proper factorization of the $$3\times3$$ block matrix and setting a relaxation parameter $$\alpha$$. The preconditioner is inspired by the relaxed dimensional factorization introduced by Benzi et al. [J. Comput. Phys., 230 (2011), pp. 6185--6202; Comput. Methods Appl. Mech. Engrg., 300 (2016), pp. 129--145]. A stable algorithm is advanced to compute the optimal value of $$\alpha$$, along with a lower bound to control the possible ill-conditioning of the $$\alpha$$ dependent inner blocks. Lastly, numerical experiments in both theoretical benchmarks and real-world applications are presented and discussed to investigate the RPF properties, performance, and robustness.
- Research Organization:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC52-07NA27344
- OSTI ID:
- 1548353
- Report Number(s):
- LLNL-JRNL-756041; 943354
- Journal Information:
- SIAM Journal on Scientific Computing, Vol. 41, Issue 4; ISSN 1064-8275
- Publisher:
- SIAMCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Monolithic Multigrid for a Reduced-Quadrature Discretization of Poroelasticity
A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media