A two-stage preconditioner for multiphase poromechanics in reservoir simulation

White, Joshua A.; Castelletto, Nicola; Klevtsov, Sergey; Bui, Quan M.; Osei-Kuffuor, Daniel; Tchelepi, Hamdi A.

doi:10.1016/j.cma.2019.112575

Title: A two-stage preconditioner for multiphase poromechanics in reservoir simulation

Abstract

Many applications involving porous media–notably reservoir engineering and geologic applications–involve tight coupling between multiphase fluid flow, transport, and poromechanical deformation. While numerical models for these processes have become commonplace in research and industry, the poor scalability of existing solution algorithms has limited the size and resolution of models that may be practically solved. In this work, we propose a two-stage Newton–Krylov solution algorithm to address this shortfall. The proposed solver exhibits rapid convergence, good parallel scalability, and is robust in the presence of highly heterogeneous material properties. Here, the key to success of the solver is a block-preconditioning strategy that breaks the fully-coupled system of mass and momentum balance equations into simpler sub-problems that may be readily addressed using targeted algebraic methods. Numerical results are presented to illustrate the performance of the solver on challenging benchmark problems.

Authors:

; Castelletto, Nicola; Klevtsov, Sergey; Bui, Quan M.; Osei-Kuffuor, Daniel;

Publication Date:: Sun Dec 01 00:00:00 EST 2019

Research Org.:: Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)

Sponsoring Org.:: USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1564406

Alternate Identifier(s):: OSTI ID: 1569188

Report Number(s):: LLNL-JRNL-763612
Journal ID: ISSN 0045-7825; S0045782519304402; 112575; PII: S0045782519304402

Grant/Contract Number:: AC52-07NA27344

Resource Type:: Published Article

Journal Name:: Computer Methods in Applied Mechanics and Engineering

Additional Journal Information:: Journal Name: Computer Methods in Applied Mechanics and Engineering Journal Volume: 357 Journal Issue: C; Journal ID: ISSN 0045-7825

Publisher:: Elsevier

Country of Publication:: Netherlands

Language:: English

Subject:: 58 GEOSCIENCES; 98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL PROTECTION

Citation Formats


                    White, Joshua A., Castelletto, Nicola, Klevtsov, Sergey, Bui, Quan M., Osei-Kuffuor, Daniel, and Tchelepi, Hamdi A. A two-stage preconditioner for multiphase poromechanics in reservoir simulation.  Netherlands: N. p., 2019. 
Web.  doi:10.1016/j.cma.2019.112575.

Copy to clipboard


                    White, Joshua A., Castelletto, Nicola, Klevtsov, Sergey, Bui, Quan M., Osei-Kuffuor, Daniel, & Tchelepi, Hamdi A. A two-stage preconditioner for multiphase poromechanics in reservoir simulation.  Netherlands.  https://doi.org/10.1016/j.cma.2019.112575

Copy to clipboard


                    White, Joshua A., Castelletto, Nicola, Klevtsov, Sergey, Bui, Quan M., Osei-Kuffuor, Daniel, and Tchelepi, Hamdi A. Sun .  
"A two-stage preconditioner for multiphase poromechanics in reservoir simulation".  Netherlands.  https://doi.org/10.1016/j.cma.2019.112575.

Copy to clipboard


                    
@article{osti_1564406,

  title        = {A two-stage preconditioner for multiphase poromechanics in reservoir simulation},

  author       = {White, Joshua A. and Castelletto, Nicola and Klevtsov, Sergey and Bui, Quan M. and Osei-Kuffuor, Daniel and Tchelepi, Hamdi A.},

  abstractNote = {Many applications involving porous media–notably reservoir engineering and geologic applications–involve tight coupling between multiphase fluid flow, transport, and poromechanical deformation. While numerical models for these processes have become commonplace in research and industry, the poor scalability of existing solution algorithms has limited the size and resolution of models that may be practically solved. In this work, we propose a two-stage Newton–Krylov solution algorithm to address this shortfall. The proposed solver exhibits rapid convergence, good parallel scalability, and is robust in the presence of highly heterogeneous material properties. Here, the key to success of the solver is a block-preconditioning strategy that breaks the fully-coupled system of mass and momentum balance equations into simpler sub-problems that may be readily addressed using targeted algebraic methods. Numerical results are presented to illustrate the performance of the solver on challenging benchmark problems.},

  doi          = {10.1016/j.cma.2019.112575},

  journal      = {Computer Methods in Applied Mechanics and Engineering},

  number       = C,

  volume       = 357,

  place        = {Netherlands},

  year         = {Sun Dec 01 00:00:00 EST 2019},

  month        = {Sun Dec 01 00:00:00 EST 2019}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Publisher's Version of Record
https://doi.org/10.1016/j.cma.2019.112575

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 28 works

Citation information provided by
Web of Science

Figures / Tables:

Algorithm 1: Application of the preconditioner M^-1

All figures and tables (7 total)

Save / Share:

Export Metadata

Save to My Library

Works referencing / citing this record:

Aspects of Solvers for Large-Scale Coupled Problems in Porous Media
journal, September 2019

Nägel, Arne; Logashenko, Dmitry; Schroder, Jacob B.
Transport in Porous Media, Vol. 130, Issue 1
DOI: 10.1007/s11242-019-01323-w

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

Similar Records in DOE PAGES and OSTI.GOV collections:

Preconditioners for multiphase poromechanics with strong capillarity

Journal Article Camargo, Julia T. ; White, Joshua A. ; Castelletto, Nicola ; ... - International Journal for Numerical and Analytical Methods in Geomechanics

Abstract This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two‐phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection‐dominated to diffusion‐dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.
https://doi.org/10.1002/nag.3192

Full Text Available
A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Journal Article Bui, Quan M. ; Osei-Kuffuor, Daniel ; Castelletto, Nicola ; ... - SIAM Journal on Scientific Computing

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-artmore »« less
Cited by 19
https://doi.org/10.1137/19M1256117

Full Text Available
Linear and nonlinear solvers for simulating multiphase flow within large-scale engineered subsurface systems

Journal Article Park, Heeho D. ; Hammond, Glenn E. ; Valocchi, Albert J. ; ... - Advances in Water Resources

Multiphase flow simulation is well-known to be computationally demanding, and modeling large-scale engineered subsurface systems entails significant additional numerical challenges. These challenges arise from: (a) the presence of small-scale discrete features like shafts, tunnels, waste packages, and barriers; (b) the need to accurately represent both the waste form processes at the small spatial scale of the repository and the large-scale transport processes throughout heterogeneous geological formations; (c) the strong contrast in material properties such as porosity and permeability, as well as the nonlinear constitutive relations for multiphase flow. Numerical solution entails discretization of the coupled system of nonlinear governing equationsmore »« less
https://doi.org/10.1016/j.advwatres.2021.104029

Full Text Available
Initial performance of fully-coupled AMG and approximate block factorization preconditioners for solution of implicit FE resistive MHD.

Conference Shadid, John Nicolas ; Lin, Paul Tinphone ; Pawlowski, Roger Patrick ; ...

This brief paper explores the development of scalable, nonlinear, fully-implicit solution methods for a stabilized unstructured finite element (FE) discretization of the 2D incompressible (reduced) resistive MHD system. The discussion considers the stabilized FE formulation in context of a fully-implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton-Krylov methods, which are preconditioned using fully-coupled algebraic multilevel (AMG) techniques and a new approximate block factorization (ABF) preconditioner. The intent of these preconditioners is to enable robust, scalable and efficient solution approaches for the large-scale sparse linear systems generated by the Newton linearization. We present results for themore »« less
Enhanced relaxed physical factorization preconditioner for coupled poromechanics

Journal Article Frigo, Matteo ; Castelletto, Nicola ; Ferronato, Massimiliano - Computers and Mathematics with Applications (Oxford)

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 + βFF^T), where C is a symmetric positive definite matrix, FF^T 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,more »« less
https://doi.org/10.1016/j.camwa.2021.11.015

Full Text Available

Similar Records

Title: A two-stage preconditioner for multiphase poromechanics in reservoir simulation

Abstract

Citation Formats

Figures / Tables:

Aspects of Solvers for Large-Scale Coupled Problems in Porous Media journal, September 2019

Aspects of Solvers for Large-Scale Coupled Problems in Porous Media
journal, September 2019