Enhanced Relaxed Physical Factorization preconditioner for coupled poromechanics

07/29/2020
by   Matteo Frigo, et al.
0

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.

READ FULL TEXT
research
07/11/2021

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-...
research
01/11/2021

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...
research
04/04/2022

A Vanka-based parameter-robust multigrid relaxation for the Stokes-Darcy Brinkman problems

We propose a block-structured multigrid relaxation scheme for solving th...
research
03/12/2023

hp-Multigrid preconditioner for a divergence-conforming HDG scheme for the incompressible flow problems

In this study, we present an hp-multigrid preconditioner for a divergenc...
research
12/01/2020

Bifurcation Analysis of the Eigenstructure of the Discrete Single-curl Operator in Three-dimensional Maxwell's Equations with Pasteur Media

This paper focuses on studying the bifurcation analysis of the eigenstru...
research
06/18/2021

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

We introduce a software package called HIFIR for preconditioning sparse,...

Please sign up or login with your details

Forgot password? Click here to reset