DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Enhanced relaxed physical factorization preconditioner for coupled poromechanics

Abstract

The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled poromechanics. For its application, however, it is necessary to invert blocks with the algebraic form C^ = (C + βFFT), where C is a symmetric positive definite matrix, FFT a rank-deficient term, and β a real non-negative coefficient. The inversion of C^, performed in an inexact way, can become unstable for large values of β, as it usually occurs at some stages of a full poromechanical simulation. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with C^. This strategy can prove useful in other problems as well where such an issue might arise, such as augmented Lagrangian preconditioning techniques for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix C^. Second, we develop a technique based on the use of a proper projection operator annihilating the near-kernel modes of C^. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Furthermore, effectiveness and robustness of the proposed algorithms are demonstratedmore » in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.« less

Authors:
 [1];  [2]; ORCiD logo [1]
  1. Univ. of Padova (Italy)
  2. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1841859
Report Number(s):
LLNL-JRNL-798397
Journal ID: ISSN 0898-1221; 1000384
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
Computers and Mathematics with Applications (Oxford)
Additional Journal Information:
Journal Name: Computers and Mathematics with Applications (Oxford); Journal Volume: 106; Journal ID: ISSN 0898-1221
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Preconditioning; Krylov subspace methods; Poromechanics

Citation Formats

Frigo, Matteo, Castelletto, Nicola, and Ferronato, Massimiliano. Enhanced relaxed physical factorization preconditioner for coupled poromechanics. United States: N. p., 2021. Web. doi:10.1016/j.camwa.2021.11.015.
Frigo, Matteo, Castelletto, Nicola, & Ferronato, Massimiliano. Enhanced relaxed physical factorization preconditioner for coupled poromechanics. United States. https://doi.org/10.1016/j.camwa.2021.11.015
Frigo, Matteo, Castelletto, Nicola, and Ferronato, Massimiliano. Wed . "Enhanced relaxed physical factorization preconditioner for coupled poromechanics". United States. https://doi.org/10.1016/j.camwa.2021.11.015. https://www.osti.gov/servlets/purl/1841859.
@article{osti_1841859,
title = {Enhanced relaxed physical factorization preconditioner for coupled poromechanics},
author = {Frigo, Matteo and Castelletto, Nicola and Ferronato, Massimiliano},
abstractNote = {The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled poromechanics. For its application, however, it is necessary to invert blocks with the algebraic form C^ = (C + βFFT), where C is a symmetric positive definite matrix, FFT a rank-deficient term, and β a real non-negative coefficient. The inversion of C^, performed in an inexact way, can become unstable for large values of β, as it usually occurs at some stages of a full poromechanical simulation. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with C^. This strategy can prove useful in other problems as well where such an issue might arise, such as augmented Lagrangian preconditioning techniques for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix C^. Second, we develop a technique based on the use of a proper projection operator annihilating the near-kernel modes of C^. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Furthermore, effectiveness and robustness of the proposed algorithms are demonstrated in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.},
doi = {10.1016/j.camwa.2021.11.015},
journal = {Computers and Mathematics with Applications (Oxford)},
number = ,
volume = 106,
place = {United States},
year = {Wed Dec 15 00:00:00 EST 2021},
month = {Wed Dec 15 00:00:00 EST 2021}
}

Works referenced in this record:

Statistic inversion of multi-zone transition probability models for aquifer characterization in alluvial fans
journal, June 2015

  • Zhu, Lin; Dai, Zhenxue; Gong, Huili
  • Stochastic Environmental Research and Risk Assessment, Vol. 30, Issue 3
  • DOI: 10.1007/s00477-015-1089-2

A stabilized element-based finite volume method for poroelastic problems
journal, July 2018

  • Honório, Hermínio T.; Maliska, Clovis R.; Ferronato, Massimiliano
  • Journal of Computational Physics, Vol. 364
  • DOI: 10.1016/j.jcp.2018.03.010

A nonconforming finite element method for the Biot’s consolidation model in poroelasticity
journal, January 2017

  • Hu, Xiaozhe; Rodrigo, Carmen; Gaspar, Francisco J.
  • Journal of Computational and Applied Mathematics, Vol. 310
  • DOI: 10.1016/j.cam.2016.06.003

On an Uzawa smoother in multigrid for poroelasticity equations: On an Uzawa smoother in multigrid for poroelasticity equations
journal, November 2016

  • Luo, P.; Rodrigo, C.; Gaspar, F. J.
  • Numerical Linear Algebra with Applications, Vol. 24, Issue 1
  • DOI: 10.1002/nla.2074

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media
journal, January 2020

  • Bui, Quan M.; Osei-Kuffuor, Daniel; Castelletto, Nicola
  • SIAM Journal on Scientific Computing, Vol. 42, Issue 2
  • DOI: 10.1137/19M1256117

Parameter-Robust Discretization and Preconditioning of Biot's Consolidation Model
journal, January 2017

  • Lee, Jeonghun J.; Mardal, Kent-Andre; Winther, Ragnar
  • SIAM Journal on Scientific Computing, Vol. 39, Issue 1
  • DOI: 10.1137/15M1029473

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: the discrete-in-time case
journal, March 2007


Anderson accelerated fixed-stress splitting schemes for consolidation of unsaturated porous media
journal, March 2019

  • Both, Jakub Wiktor; Kumar, Kundan; Nordbotten, Jan Martin
  • Computers & Mathematics with Applications, Vol. 77, Issue 6
  • DOI: 10.1016/j.camwa.2018.07.033

Convergence of iterative coupling of geomechanics with flow in a fractured poroelastic medium
journal, June 2016

  • Girault, Vivette; Kumar, Kundan; Wheeler, Mary F.
  • Computational Geosciences, Vol. 20, Issue 5
  • DOI: 10.1007/s10596-016-9573-4

Efficient solvers for hybridized three-field mixed finite element coupled poromechanics
journal, June 2021

  • Frigo, Matteo; Castelletto, Nicola; Ferronato, Massimiliano
  • Computers & Mathematics with Applications, Vol. 91
  • DOI: 10.1016/j.camwa.2020.07.010

Robust fixed stress splitting for Biot’s equations in heterogeneous media
journal, June 2017

  • Both, Jakub Wiktor; Borregales, Manuel; Nordbotten, Jan Martin
  • Applied Mathematics Letters, Vol. 68
  • DOI: 10.1016/j.aml.2016.12.019

Convergence analysis of multirate fixed-stress split iterative schemes for coupling flow with geomechanics
journal, November 2016

  • Almani, T.; Kumar, K.; Dogru, A.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 311
  • DOI: 10.1016/j.cma.2016.07.036

A locally conservative finite element framework for the simulation of coupled flow and reservoir geomechanics
journal, July 2007


Parameter estimates for the Relaxed Dimensional Factorization preconditioner and application to hemodynamics
journal, March 2016

  • Benzi, Michele; Deparis, Simone; Grandperrin, Gwenol
  • Computer Methods in Applied Mechanics and Engineering, Vol. 300
  • DOI: 10.1016/j.cma.2015.11.016

A fully coupled 3-D mixed finite element model of Biot consolidation
journal, June 2010

  • Ferronato, Massimiliano; Castelletto, Nicola; Gambolati, Giuseppe
  • Journal of Computational Physics, Vol. 229, Issue 12
  • DOI: 10.1016/j.jcp.2010.03.018

GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986

  • Saad, Youcef; Schultz, Martin H.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
  • DOI: 10.1137/0907058

A multiscale fixed stress split iterative scheme for coupled flow and poromechanics in deep subsurface reservoirs
journal, January 2018


A Relaxed Physical Factorization Preconditioner for Mixed Finite Element Coupled Poromechanics
journal, January 2019

  • Frigo, Matteo; Castelletto, Nicola; Ferronato, Massimiliano
  • SIAM Journal on Scientific Computing, Vol. 41, Issue 4
  • DOI: 10.1137/18M120645X

New stabilized discretizations for poroelasticity and the Stokes’ equations
journal, November 2018

  • Rodrigo, C.; Hu, X.; Ohm, P.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 341
  • DOI: 10.1016/j.cma.2018.07.003

A macroelement stabilization for mixed finite element/finite volume discretizations of multiphase poromechanics
journal, June 2020

  • Camargo, Julia T.; White, Joshua A.; Borja, Ronaldo I.
  • Computational Geosciences, Vol. 25, Issue 2
  • DOI: 10.1007/s10596-020-09964-3

Scalable algorithms for three-field mixed finite element coupled poromechanics
journal, December 2016

  • Castelletto, Nicola; White, Joshua A.; Ferronato, Massimiliano
  • Journal of Computational Physics, Vol. 327
  • DOI: 10.1016/j.jcp.2016.09.063

On the causes of pressure oscillations in low-permeable and low-compressible porous media: PRESSURE OSCILLATIONS IN LOW-PERMEABLE POROUS MEDIA
journal, July 2011

  • Haga, Joachim Berdal; Osnes, Harald; Langtangen, Hans Petter
  • International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 36, Issue 12
  • DOI: 10.1002/nag.1062

An Augmented Lagrangian Approach to Linearized Problems in Hydrodynamic Stability
journal, January 2008

  • Olshanskii, Maxim A.; Benzi, Michele
  • SIAM Journal on Scientific Computing, Vol. 30, Issue 3
  • DOI: 10.1137/070691851

Parameter-Robust Convergence Analysis of Fixed-Stress Split Iterative Method for Multiple-Permeability Poroelasticity Systems
journal, January 2020

  • Hong, Qingguo; Kraus, Johannes; Lymbery, Maria
  • Multiscale Modeling & Simulation, Vol. 18, Issue 2
  • DOI: 10.1137/19M1253988

A priori error estimates for a discretized poro-elastic–elastic system solved by a fixed-stress algorithm
journal, January 2019

  • Girault, Vivette; Wheeler, Mary F.; Almani, Tameem
  • Oil & Gas Science and Technology – Revue d’IFP Energies nouvelles, Vol. 74
  • DOI: 10.2516/ogst/2018071

Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques
journal, August 2001

  • Christie, M. A.; Blunt, M. J.
  • SPE Reservoir Evaluation & Engineering, Vol. 4, Issue 04
  • DOI: 10.2118/72469-PA

Large scale micro finite element analysis of 3D bone poroelasticity
journal, July 2014


Block-preconditioned Newton–Krylov solvers for fully coupled flow and geomechanics
journal, April 2011


A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case
journal, March 2007


A stabilized hybrid mixed finite element method for poroelasticity
journal, July 2020


Convergence analysis of two-grid fixed stress split iterative scheme for coupled flow and deformation in heterogeneous poroelastic media
journal, November 2018

  • Dana, Saumik; Wheeler, Mary F.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 341
  • DOI: 10.1016/j.cma.2018.07.018

Stable discretization of poroelasticity problems and efficient preconditioners for arising saddle point type matrices
journal, August 2012

  • Axelsson, Owe; Blaheta, Radim; Byczanski, Petr
  • Computing and Visualization in Science, Vol. 15, Issue 4
  • DOI: 10.1007/s00791-013-0209-0

An Augmented Lagrangian Preconditioner for the 3D Stationary Incompressible Navier--Stokes Equations at High Reynolds Number
journal, January 2019

  • Farrell, Patrick E.; Mitchell, Lawrence; Wechsung, Florian
  • SIAM Journal on Scientific Computing, Vol. 41, Issue 5
  • DOI: 10.1137/18M1219370

Modified augmented Lagrangian preconditioners for the incompressible Navier-Stokes equations
journal, May 2011

  • Benzi, Michele; Olshanskii, Maxim A.; Wang, Zhen
  • International Journal for Numerical Methods in Fluids, Vol. 66, Issue 4
  • DOI: 10.1002/fld.2267

Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems
journal, March 1992

  • van der Vorst, H. A.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 2
  • DOI: 10.1137/0913035

Robust iterative schemes for non-linear poromechanics
journal, April 2018

  • Borregales, Manuel; Radu, Florin A.; Kumar, Kundan
  • Computational Geosciences, Vol. 22, Issue 4
  • DOI: 10.1007/s10596-018-9736-6

General Theory of Three‐Dimensional Consolidation
journal, February 1941

  • Biot, Maurice A.
  • Journal of Applied Physics, Vol. 12, Issue 2
  • DOI: 10.1063/1.1712886

A general preconditioning framework for coupled multiphysics problems with application to contact- and poro-mechanics
journal, December 2019

  • Ferronato, Massimiliano; Franceschini, Andrea; Janna, Carlo
  • Journal of Computational Physics, Vol. 398
  • DOI: 10.1016/j.jcp.2019.108887

A parallel block preconditioner for large-scale poroelasticity with highly heterogeneous material parameters
journal, March 2012

  • Haga, Joachim Berdal; Osnes, Harald; Langtangen, Hans Petter
  • Computational Geosciences, Vol. 16, Issue 3
  • DOI: 10.1007/s10596-012-9284-4

Accuracy and convergence properties of the fixed-stress iterative solution of two-way coupled poromechanics: Accuracy and convergence properties of the fixed-stress iterative solution of two-way coupled poromechanics
journal, June 2015

  • Castelletto, N.; White, J. A.; Tchelepi, H. A.
  • International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 39, Issue 14
  • DOI: 10.1002/nag.2400

Robust Preconditioners for a New Stabilized Discretization of the Poroelastic Equations
journal, January 2020

  • Adler, J. H.; Gaspar, F. J.; Hu, X.
  • SIAM Journal on Scientific Computing, Vol. 42, Issue 3
  • DOI: 10.1137/19M1261250

A Mixed Finite Element--Finite Volume Formulation of the Black-Oil Model
journal, January 1998

  • Bergamaschi, Luca; Mantica, Stefano; Manzini, Gianmarco
  • SIAM Journal on Scientific Computing, Vol. 20, Issue 3
  • DOI: 10.1137/S1064827595289303

On the fixed-stress split scheme as smoother in multigrid methods for coupling flow and geomechanics
journal, November 2017

  • Gaspar, Francisco J.; Rodrigo, Carmen
  • Computer Methods in Applied Mechanics and Engineering, Vol. 326
  • DOI: 10.1016/j.cma.2017.08.025

Field-of-Values Convergence Analysis of Augmented Lagrangian Preconditioners for the Linearized Navier–Stokes Problem
journal, January 2011

  • Benzi, Michele; Olshanskii, Maxim A.
  • SIAM Journal on Numerical Analysis, Vol. 49, Issue 2
  • DOI: 10.1137/100806485

Consolidation Des Sols (Étude Mathématique)
journal, September 1953


The 3‐D Facies and Geomechanical Modeling of Land Subsidence in the Chaobai Plain, Beijing
journal, March 2020

  • Zhu, Lin; Franceschini, Andrea; Gong, Huili
  • Water Resources Research, Vol. 56, Issue 3
  • DOI: 10.1029/2019WR027026

A Relaxed Dimensional Factorization preconditioner for the incompressible Navier–Stokes equations
journal, July 2011

  • Benzi, Michele; Ng, Michael; Niu, Qiang
  • Journal of Computational Physics, Vol. 230, Issue 16
  • DOI: 10.1016/j.jcp.2011.04.001

ILUT: A dual threshold incomplete LU factorization
journal, July 1994


Parameter-robust stability of classical three-field formulation of Biot's consolidation model
journal, January 2018

  • Hong, Qingguo; Kraus, Johannes
  • ETNA - Electronic Transactions on Numerical Analysis, Vol. 48
  • DOI: 10.1553/etna_vol48s202

On some mixed finite element methods for incompressible and nearly incompressible finite elasticity
journal, November 1996


HSL_MI20 : An efficient AMG preconditioner for finite element problems in 3D: HSL_MI20 : AN EFFICIENT AMG PRECONDITIONER
journal, October 2009

  • Boyle, Jonathan; Mihajlović, Milan; Scott, Jennifer
  • International Journal for Numerical Methods in Engineering, Vol. 82, Issue 1
  • DOI: 10.1002/nme.2758

Space–time finite element approximation of the Biot poroelasticity system with iterative coupling
journal, June 2017

  • Bause, M.; Radu, F. A.; Köcher, U.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 320
  • DOI: 10.1016/j.cma.2017.03.017

Incomplete Cholesky Factorizations with Limited Memory
journal, January 1999


Block-partitioned solvers for coupled poromechanics: A unified framework
journal, May 2016

  • White, Joshua A.; Castelletto, Nicola; Tchelepi, Hamdi A.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 303
  • DOI: 10.1016/j.cma.2016.01.008

Enriched Galerkin finite elements for coupled poromechanics with local mass conservation
journal, November 2018

  • Choo, Jinhyun; Lee, Sanghyun
  • Computer Methods in Applied Mechanics and Engineering, Vol. 341
  • DOI: 10.1016/j.cma.2018.06.022

Scalable algorithms for three-field mixed finite element coupled poromechanics
journal, December 2016

  • Castelletto, Nicola; White, Joshua A.; Ferronato, Massimiliano
  • Journal of Computational Physics, Vol. 327
  • DOI: 10.1016/j.jcp.2016.09.063

New stabilized discretizations for poroelasticity and the Stokes' equations
text, January 2017