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

Title: Multiscale two-stage solver for Biot’s poroelasticity equations in subsurface media

Abstract

In this paper, we propose a two-stage preconditioner for accelerating the iterative solution by a Krylov subspace method of Biot’s poroelasticity equations based on a displacement-pressure formulation. The spatial discretization combines a finite element method for mechanics and a finite volume approach for flow. The fully implicit backward Euler scheme is used for time integration. The result is a 2 × 2 block linear system for each timestep. The preconditioning operator is obtained by applying a two-stage scheme. The first stage is a global preconditioner that employs multiscale basis functions to construct coarse-scale coupled systems using a Galerkin projection. This global stage is effective at damping low-frequency error modes associated with long-range coupling of the unknowns. The second stage is a local block-triangular smoothing preconditioner, which is aimed at high-frequency error modes associated with short-range coupling of the variables. Lastly, various numerical experiments are used to demonstrate the robustness of the proposed solver.

Authors:
 [1];  [2];  [3]; ORCiD logo [2]
  1. Stanford Univ., CA (United States); Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Stanford Univ., CA (United States)
  3. Delft Univ. of Technology (Netherlands)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1548341
Report Number(s):
LLNL-JRNL-747875
Journal ID: ISSN 1420-0597; 932724
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
Computational Geosciences
Additional Journal Information:
Journal Volume: 23; Journal Issue: 2; Journal ID: ISSN 1420-0597
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; 97 MATHEMATICS AND COMPUTING; Poromechanics; Multiscale methods; Preconditioners; Iterative methods

Citation Formats

Castelletto, Nicola, Klevtsov, Sergey, Hajibeygi, Hadi, and Tchelepi, Hamdi A. Multiscale two-stage solver for Biot’s poroelasticity equations in subsurface media. United States: N. p., 2018. Web. doi:10.1007/s10596-018-9791-z.
Castelletto, Nicola, Klevtsov, Sergey, Hajibeygi, Hadi, & Tchelepi, Hamdi A. Multiscale two-stage solver for Biot’s poroelasticity equations in subsurface media. United States. https://doi.org/10.1007/s10596-018-9791-z
Castelletto, Nicola, Klevtsov, Sergey, Hajibeygi, Hadi, and Tchelepi, Hamdi A. Fri . "Multiscale two-stage solver for Biot’s poroelasticity equations in subsurface media". United States. https://doi.org/10.1007/s10596-018-9791-z. https://www.osti.gov/servlets/purl/1548341.
@article{osti_1548341,
title = {Multiscale two-stage solver for Biot’s poroelasticity equations in subsurface media},
author = {Castelletto, Nicola and Klevtsov, Sergey and Hajibeygi, Hadi and Tchelepi, Hamdi A.},
abstractNote = {In this paper, we propose a two-stage preconditioner for accelerating the iterative solution by a Krylov subspace method of Biot’s poroelasticity equations based on a displacement-pressure formulation. The spatial discretization combines a finite element method for mechanics and a finite volume approach for flow. The fully implicit backward Euler scheme is used for time integration. The result is a 2 × 2 block linear system for each timestep. The preconditioning operator is obtained by applying a two-stage scheme. The first stage is a global preconditioner that employs multiscale basis functions to construct coarse-scale coupled systems using a Galerkin projection. This global stage is effective at damping low-frequency error modes associated with long-range coupling of the unknowns. The second stage is a local block-triangular smoothing preconditioner, which is aimed at high-frequency error modes associated with short-range coupling of the variables. Lastly, various numerical experiments are used to demonstrate the robustness of the proposed solver.},
doi = {10.1007/s10596-018-9791-z},
journal = {Computational Geosciences},
number = 2,
volume = 23,
place = {United States},
year = {2018},
month = {11}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 18 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems
journal, January 1988


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

Multi-scale finite-volume method for elliptic problems in subsurface flow simulation
journal, May 2003


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

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

Adaptive fully implicit multi-scale finite-volume method for multi-phase flow and transport in heterogeneous porous media
journal, September 2006


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

Cell-centered finite volume discretizations for deformable porous media: CELL-CENTERED FINITE VOLUME DISCRETIZATIONS FOR ELASTICITY
journal, August 2014

  • Nordbotten, Jan Martin
  • International Journal for Numerical Methods in Engineering, Vol. 100, Issue 6
  • DOI: 10.1002/nme.4734

Cell-centered nonlinear finite-volume methods for the heterogeneous anisotropic diffusion problem
journal, February 2017

  • Terekhov, Kirill M.; Mallison, Bradley T.; Tchelepi, Hamdi A.
  • Journal of Computational Physics, Vol. 330
  • DOI: 10.1016/j.jcp.2016.11.010

An iterative solution method for linear systems of which the coefficient matrix is a symmetric $M$-matrix
journal, January 1977


Monotone nonlinear finite-volume method for challenging grids
journal, January 2018


Abstract robust coarse spaces for systems of PDEs via generalized eigenproblems in the overlaps
journal, August 2013


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

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


Use of Multiple Multiscale Operators To Accelerate Simulation of Complex Geomodels
journal, December 2017

  • Lie, Knut-Andreas; Møyner, Olav; Natvig, Jostein R.
  • SPE Journal, Vol. 22, Issue 06
  • DOI: 10.2118/182701-PA

Stability and convergence of sequential methods for coupled flow and geomechanics: Fixed-stress and fixed-strain splits
journal, March 2011

  • Kim, J.; Tchelepi, H. A.; Juanes, R.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 200, Issue 13-16
  • DOI: 10.1016/j.cma.2010.12.022

Convergence analysis of a new mixed finite element method for Biot's consolidation model
journal, February 2014

  • Yi, Son-Young
  • Numerical Methods for Partial Differential Equations, Vol. 30, Issue 4
  • DOI: 10.1002/num.21865

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

Convergence of iterative coupling for coupled flow and geomechanics
journal, October 2012


Thermo-hydro-mechanical modeling of fluid geological storage by Godunov-mixed methods: THERMO-HYDRO-MECHANICAL MODELING BY GODUNOV-MIXED METHODS
journal, March 2012

  • Castelletto, N.; Ferronato, M.; Gambolati, G.
  • International Journal for Numerical Methods in Engineering
  • DOI: 10.1002/nme.3352

Stability and monotonicity for some discretizations of the Biot’s consolidation model
journal, January 2016

  • Rodrigo, C.; Gaspar, F. J.; Hu, X.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 298
  • DOI: 10.1016/j.cma.2015.09.019

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

Constrained pressure residual multiscale (CPR-MS) method for fully implicit simulation of multiphase flow in porous media
journal, October 2015

  • Cusini, Matteo; Lukyanov, Alexander A.; Natvig, Jostein
  • Journal of Computational Physics, Vol. 299
  • DOI: 10.1016/j.jcp.2015.07.019

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


A systematic comparison of coupled and distributive smoothing in multigrid for the poroelasticity system
journal, March 2004

  • Gaspar, F. J.; Lisbona, F. J.; Oosterlee, C. W.
  • Numerical Linear Algebra with Applications, Vol. 11, Issue 23
  • DOI: 10.1002/nla.372

Finite volume methods for elasticity with weak symmetry: FINITE VOLUME METHODS FOR ELASTICITY WITH WEAK SYMMETRY
journal, April 2017

  • Keilegavlen, Eirik; Nordbotten, Jan Martin
  • International Journal for Numerical Methods in Engineering, Vol. 112, Issue 8
  • DOI: 10.1002/nme.5538

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

Algebraic dynamic multilevel method for compositional flow in heterogeneous porous media
journal, February 2018

  • Cusini, Matteo; Fryer, Barnaby; van Kruijsdijk, Cor
  • Journal of Computational Physics, Vol. 354
  • DOI: 10.1016/j.jcp.2017.10.052

Upscaling of Mixed Finite Element Discretization Problems by the Spectral AMGe Method
journal, January 2016

  • Kalchev, D. Z.; Lee, C. S.; Villa, U.
  • SIAM Journal on Scientific Computing, Vol. 38, Issue 5
  • DOI: 10.1137/15M1036683

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

On the relationship between the multiscale finite-volume method and domain decomposition preconditioners
journal, January 2008


Block-preconditioned Krylov Methods for Coupled Multiphase Reservoir Flow and Geomechanics
conference, August 2016

  • Klevtsov, S.; Castelletto, N.; White, J. A.
  • ECMOR XV - 15th European Conference on the Mathematics of Oil Recovery, Proceedings
  • DOI: 10.3997/2214-4609.201601900

Adaptive algebraic multiscale solver for compressible flow in heterogeneous porous media
journal, November 2015


Multiscale finite element coarse spaces for the application to linear elasticity
journal, January 2013


A multiscale restriction-smoothed basis method for high contrast porous media represented on unstructured grids
journal, January 2016


A Coupled Reservoir and Geomechanical Simulation System
journal, September 1998

  • Settari, Antonin; Mourits, F. M.
  • SPE Journal, Vol. 3, Issue 03
  • DOI: 10.2118/50939-PA

Iterative multiscale finite-volume method
journal, October 2008

  • Hajibeygi, Hadi; Bonfigli, Giuseppe; Hesse, Marc Andre
  • Journal of Computational Physics, Vol. 227, Issue 19
  • DOI: 10.1016/j.jcp.2008.06.013

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

A reduced order model-based preconditioner for the efficient solution of transient diffusion equations: A ROM-BASED PRECONDITIONER FOR TRANSIENT DIFFUSION PDES
journal, July 2016

  • Pasetto, Damiano; Ferronato, Massimiliano; Putti, Mario
  • International Journal for Numerical Methods in Engineering, Vol. 109, Issue 8
  • DOI: 10.1002/nme.5320

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


A Time-Stepping Scheme for Coupled Reservoir Flow and Geomechanics
conference, April 2013

  • Gai, Xiuli; Sun, Shuyu; Wheeler, Mary Fanett
  • SPE Annual Technical Conference and Exhibition
  • DOI: 10.2118/97054-MS

Monotone multiscale finite volume method
journal, August 2015

  • Wang, Yixuan; Hajibeygi, Hadi; Tchelepi, Hamdi A.
  • Computational Geosciences, Vol. 20, Issue 3
  • DOI: 10.1007/s10596-015-9506-7

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


Finite volume schemes for diffusion equations: Introduction to and review of modern methods
journal, May 2014


Two-Stage Algebraic Multiscale Linear Solver for Highly Heterogeneous Reservoir Models
journal, June 2012

  • Zhou, Hui; Tchelepi, Hamdi A.
  • SPE Journal, Vol. 17, Issue 02
  • DOI: 10.2118/141473-PA

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


Successful Application of Multiscale Methods in a Real Reservoir Simulator Environment
conference, August 2016

  • Lie, K. A.; Møyner, O.; Natvig, J. R.
  • ECMOR XV - 15th European Conference on the Mathematics of Oil Recovery, Proceedings
  • DOI: 10.3997/2214-4609.201601893

Stable Cell-Centered Finite Volume Discretization for Biot Equations
journal, January 2016

  • Nordbotten, Jan Martin
  • SIAM Journal on Numerical Analysis, Vol. 54, Issue 2
  • DOI: 10.1137/15M1014280

Novel preconditioners for the iterative solution to FE-discretized coupled consolidation equations
journal, May 2007

  • Bergamaschi, Luca; Ferronato, Massimiliano; Gambolati, Giuseppe
  • Computer Methods in Applied Mechanics and Engineering, Vol. 196, Issue 25-28
  • DOI: 10.1016/j.cma.2007.01.013

Numerical solution of saddle point problems
journal, April 2005


A generalized multiscale finite element method for poroelasticity problems II: Nonlinear coupling
journal, May 2016

  • Brown, Donald L.; Vasilyeva, Maria
  • Journal of Computational and Applied Mathematics, Vol. 297
  • DOI: 10.1016/j.cam.2015.11.007

A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media
journal, June 1997

  • Hou, Thomas Y.; Wu, Xiao-Hui
  • Journal of Computational Physics, Vol. 134, Issue 1
  • DOI: 10.1006/jcph.1997.5682

Algebraic multiscale solver for flow in heterogeneous porous media
journal, February 2014

  • Wang, Yixuan; Hajibeygi, Hadi; Tchelepi, Hamdi A.
  • Journal of Computational Physics, Vol. 259
  • DOI: 10.1016/j.jcp.2013.11.024

TP or not TP, that is the question
journal, January 2014


Multigrid method for nonlinear poroelasticity equations
journal, October 2015

  • Luo, P.; Rodrigo, C.; Gaspar, F. J.
  • Computing and Visualization in Science, Vol. 17, Issue 5
  • DOI: 10.1007/s00791-016-0260-8

A Real-Field Multiscale Black-Oil Reservoir Simulator
conference, February 2015

  • Kozlova, A.; Li, Z.; Natvig, J. R.
  • SPE Reservoir Simulation Symposium
  • DOI: 10.2118/173226-MS

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

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

Multiscale model reduction for shale gas transport in poroelastic fractured media
journal, January 2018

  • Akkutlu, I. Yucel; Efendiev, Yalchin; Vasilyeva, Maria
  • Journal of Computational Physics, Vol. 353
  • DOI: 10.1016/j.jcp.2017.10.023

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

RMCP: Relaxed Mixed Constraint Preconditioners for saddle point linear systems arising in geomechanics
journal, May 2012

  • Bergamaschi, Luca; Martínez, Ángeles
  • Computer Methods in Applied Mechanics and Engineering, Vol. 221-222
  • DOI: 10.1016/j.cma.2012.02.004

Algebraic multiscale method for flow in heterogeneous porous media with embedded discrete fractures (F-AMS)
journal, September 2016

  • Ţene, Matei; Al Kobaisi, Mohammed Saad; Hajibeygi, Hadi
  • Journal of Computational Physics, Vol. 321
  • DOI: 10.1016/j.jcp.2016.06.012

Multiscale finite-element method for linear elastic geomechanics
journal, February 2017

  • Castelletto, Nicola; Hajibeygi, Hadi; Tchelepi, Hamdi A.
  • Journal of Computational Physics, Vol. 331
  • DOI: 10.1016/j.jcp.2016.11.044

Stabilized low-order finite elements for coupled solid-deformation/fluid-diffusion and their application to fault zone transients
journal, September 2008

  • White, Joshua A.; Borja, Ronaldo I.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issue 49-50
  • DOI: 10.1016/j.cma.2008.05.015

On stability and convergence of finite element approximations of Biot's consolidation problem
journal, February 1994

  • Murad, Márcio A.; Loula, Abimael F. D.
  • International Journal for Numerical Methods in Engineering, Vol. 37, Issue 4
  • DOI: 10.1002/nme.1620370407

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

Hybrid Multiscale Formulation for Coupled Flow and Geomechanics
conference, August 2016

  • Castelletto, N.; Hajibeygi, H.; Tchelepi, H. A.
  • ECMOR XV - 15th European Conference on the Mathematics of Oil Recovery, Proceedings
  • DOI: 10.3997/2214-4609.201601888

Preconditioning Techniques for Large Linear Systems: A Survey
journal, November 2002


A Generalized Multiscale Finite Element Method for poroelasticity problems I: Linear problems
journal, March 2016

  • Brown, Donald L.; Vasilyeva, Maria
  • Journal of Computational and Applied Mathematics, Vol. 294
  • DOI: 10.1016/j.cam.2015.08.007

Stability, Accuracy, and Efficiency of Sequential Methods for Coupled Flow and Geomechanics
journal, June 2011

  • Kim, Jihoon; Tchelepi, Hamdi A.; Juanes, Ruben
  • SPE Journal, Vol. 16, Issue 02
  • DOI: 10.2118/119084-PA

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

Coupling multiscale finite element method for consolidation analysis of heterogeneous saturated porous media
journal, February 2009


Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques
conference, April 2013

  • Christie, M. A.; Blunt, M. J.
  • SPE Reservoir Simulation Symposium
  • DOI: 10.2118/66599-ms

Successful application of multiscale methods in a real reservoir simulator environment
journal, March 2017


A Locally Conservative Finite Element Framework for the Simulation of Coupled Flow and Reservoir Geomechanics
conference, September 2006

  • Juanes, Ruben; Jha, Birendra
  • SPE Annual Technical Conference and Exhibition, All Days
  • DOI: 10.2118/102812-ms

An Iterative Solution Method for Linear Systems of Which the Coefficient Matrix is a Symmetric M-Matrix
journal, January 1977

  • Meijerink, J. A.; van der Vorst, H. A.
  • Mathematics of Computation, Vol. 31, Issue 137
  • DOI: 10.2307/2005786

Two-Stage Algebraic Multiscale Linear Solver for Highly Heterogeneous Reservoir Models
conference, February 2011

  • Zhou, H. .; Tchelepi, H. A.
  • SPE Reservoir Simulation Symposium, All Days
  • DOI: 10.2118/141473-ms

Adaptive Algebraic Multiscale Solver for Compressible Flow in Heterogeneous Porous Media
text, January 2017


Algebraic multiscale method for flow in heterogeneous porous media with embedded discrete fractures (F-AMS)
text, January 2017


Anderson accelerated fixed-stress splitting schemes for consolidation of unsaturated porous media
text, January 2018


Works referencing / citing this record:

On the optimization of the fixed‐stress splitting for Biot's equations
journal, June 2019

  • Storvik, Erlend; Both, Jakub W.; Kumar, Kundan
  • International Journal for Numerical Methods in Engineering, Vol. 120, Issue 2
  • DOI: 10.1002/nme.6130