View in EDS

Efficient Parameter-Robust Preconditioners for Linear Poroelasticity and Elasticity in the Primal Formulation.

Saved in:
Bibliographic Details
Title: Efficient Parameter-Robust Preconditioners for Linear Poroelasticity and Elasticity in the Primal Formulation.
Authors: Huang, Weizhang1 (AUTHOR) whuang@ku.edu, Wang, Zhuoran1 (AUTHOR) wangzr@ku.edu
Source: SIAM Journal on Scientific Computing. 2026, Vol. 48 Issue 3, pA1536-A1563. 28p.
Subjects: Poroelasticity, Schur complement, Numerical analysis, Elasticity, Finite element method, Iterative methods (Mathematics), Krylov subspace
Abstract: Poroelasticity problems play an important role in various engineering, geophysical, and biological applications. Their full discretization results in a large-scale saddle-point system at each time step that is becoming singular for locking cases and needs effective preconditioners for its fast iterative solution. Instead of constructing spectrally equivalent ones, we develop nonsingular preconditioners so that the eigenvalues of the preconditioned system consist of a cluster around 1 and an outlier in the order of \(1/\lambda\) , where \(\lambda\) is a Lamé constant that is large for locking cases. It is known that the convergence factor of GMRES is bounded by the radius of the cluster for this type of system. Both two- and three-field block triangular Schur complement preconditioners are studied. Upper bounds of the radius of the eigenvalue cluster for those systems are obtained and shown to be related to the inf-sup condition but independent of mesh size, time step, and the locking parameter \(\lambda\) , which reflects the robustness of the preconditioners with respect to parameter variations. Moreover, the developed preconditioners do not need to compute the Schur complement, and neither requires exact inversion of diagonal blocks except the leading one. A locking-free weak Galerkin finite element method and the implicit Euler scheme are used for the discretization of the governing equation. Both two- and three-dimensional numerical results are presented to confirm the effectiveness and parameter robustness of the developed preconditioners. [ABSTRACT FROM AUTHOR]
Copyright of SIAM Journal on Scientific Computing is the property of Society for Industrial & Applied Mathematics and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Engineering Source
Description
Abstract:Poroelasticity problems play an important role in various engineering, geophysical, and biological applications. Their full discretization results in a large-scale saddle-point system at each time step that is becoming singular for locking cases and needs effective preconditioners for its fast iterative solution. Instead of constructing spectrally equivalent ones, we develop nonsingular preconditioners so that the eigenvalues of the preconditioned system consist of a cluster around 1 and an outlier in the order of \(1/\lambda\) , where \(\lambda\) is a Lamé constant that is large for locking cases. It is known that the convergence factor of GMRES is bounded by the radius of the cluster for this type of system. Both two- and three-field block triangular Schur complement preconditioners are studied. Upper bounds of the radius of the eigenvalue cluster for those systems are obtained and shown to be related to the inf-sup condition but independent of mesh size, time step, and the locking parameter \(\lambda\) , which reflects the robustness of the preconditioners with respect to parameter variations. Moreover, the developed preconditioners do not need to compute the Schur complement, and neither requires exact inversion of diagonal blocks except the leading one. A locking-free weak Galerkin finite element method and the implicit Euler scheme are used for the discretization of the governing equation. Both two- and three-dimensional numerical results are presented to confirm the effectiveness and parameter robustness of the developed preconditioners. [ABSTRACT FROM AUTHOR]
ISSN:10648275
DOI:10.1137/24M1672882