Efficient parameter-robust preconditioners for linear poroelasticity and elasticity in the primal formulation

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Huang, Weizhang, Wang, Zhuoran
Format: Preprint
Publié: 2025
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866912451647242240
author Huang, Weizhang
Wang, Zhuoran
author_facet Huang, Weizhang
Wang, Zhuoran
contents 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/λ$, where $λ$ 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 systems. 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 locking parameters, 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 require 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.
format Preprint
id arxiv_https___arxiv_org_abs_2506_21361
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Efficient parameter-robust preconditioners for linear poroelasticity and elasticity in the primal formulation
Huang, Weizhang
Wang, Zhuoran
Numerical Analysis
65M60, 65F08, 65F10, 74F10
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/λ$, where $λ$ 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 systems. 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 locking parameters, 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 require 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.
title Efficient parameter-robust preconditioners for linear poroelasticity and elasticity in the primal formulation
topic Numerical Analysis
65M60, 65F08, 65F10, 74F10
url https://arxiv.org/abs/2506.21361