Enhanced Relaxed Physical Factorization preconditioner for coupled poromechanics

by   Matteo Frigo, et al.

In this work, we focus on the relaxed physical factorization (RPF) preconditioner for the block linear systems arising from the three-field formulation (displacement/velocity/pressure) of coupled poromechanics. Inspired by the relaxed dimensional factorization developed for the Navier-Stokes equations by Benzi et al. [J. Comput. Phys., 230 (2011), pp. 6185–6202; Comput. Methods Appl. Mech. Engrg., 300 (2016), pp. 129–145], the RPF preconditioner was recently advanced by Frigo et al. [SIAM J. Sci. Comp., 41 (2019), pp. B694–B720] and relies on: (i) combining proper physics-based splittings of the block matrix by field type, and (ii) introducing an optimal relaxation parameter α. However, a possible drawback arises from the need of inverting blocks in the form Ĉ = ( C + β F F^T) for large values of the real coefficient β, where C is a regular square matrix and FF^T is a rank-deficient term. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with Ĉ, which can also prove useful in other problems where such an issue might arise, such as augmented Lagrangian preconditioning strategies for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix Ĉ. Second, we develop a technique based on the use of a proper projection operator onto the range of F. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Effectiveness and robustness of the proposed algorithms are demonstrated in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.



There are no comments yet.


page 17


A projection-based, semi-implicit time-stepping approach for the Cahn-Hilliard Navier-Stokes equations on adaptive octree meshes

We present a projection-based framework for solving a thermodynamically-...

Exact pressure elimination for the Crouzeix-Raviart scheme applied to the Stokes and Navier-Stokes problems

We show that, using the Crouzeix-Raviart scheme, a cheap algebraic trans...

An Energy-Stable Scheme for Incompressible Navier-Stokes Equations with Periodically Updated Coefficient Matrix

We present an energy-stable scheme for simulating the incompressible Nav...

Preconditioners for computing multiple solutions in three-dimensional fluid topology optimization

Topology optimization problems generally support multiple local minima, ...

HIFIR: Hybrid Incomplete Factorization with Iterative Refinement for Preconditioning Ill-conditioned and Singular Systems

We introduce a software package called HIFIR for preconditioning sparse,...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.