A new numerically stable sequential algorithm for coupled finite-strain elastoplastic geomechanics and flow
In this work, we propose a new numerically stable sequential implicit method for coupled flow and finite-strain multiplicative elastoplastic geomechanics. We find from stability analysis that the sequential method that solves the flow problem first by fixing the first Piola-Kirchhoff total stress, solving the geomechanics (solid deformation) problem at the next step, is unconditionally stable (contractive and B-stable). In this sense, this method named the fixed first Piola-Kirchhoff stress method is an extension of the fixed stress method in coupled flow and infinitesimal geomechanics. We also study the fixed second Piola-Kirchhoff stress method, comparing it with the fixed first Piola-Kirchhoff stress method, because the constitutive relations are formulated by the second Piola-Kirchhoff total stress, although fixing the second Piola-Kirchhoff total stress field does not provide theoretical unconditional stability. In space discretization, we use the finite element method for the geomechanics problem with the total Lagrangian approach, while employing the finite volume method for the flow problem. Geometrical nonlinearity from the total Lagrangian approach results in full-tensor permeability even if the initial permeability is isotropic. To consider the full-tensor permeability accurately, we employ the multipoint flux approximation in solving the flow problem. In time discretization, the backward Euler method is used. Then, we perform numerical experiments of the two fixed-stress sequential methods with various scenarios, and find superiority of their numerical stability. Specifically, the two fixed-stress sequential methods provide almost identical results for elasticity. For J2 plasticity, the two sequential methods yield slightly different numerical results although the trends of the results are similar. For the Drucker-Prager plasticity, we find that the results from the two methods are almost identical. For all the test cases, the two sequential methods are numerically stable.
- Research Organization:
- Texas A & M Univ., College Station, TX (United States). Texas A & M Engineering Experiment Station
- Sponsoring Organization:
- USDOE Office of Fossil Energy (FE)
- Grant/Contract Number:
- FE0028973
- OSTI ID:
- 1845977
- Alternate ID(s):
- OSTI ID: 1548549
- Report Number(s):
- DOE-TAMU-FE0028973-11
- Journal Information:
- Computer Methods in Applied Mechanics and Engineering, Vol. 335; ISSN 0045-7825
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Large deformation poromechanics with local mass conservation: An enriched Galerkin finite element framework: Large deformation poromechanics with local mass conservation
|
journal | August 2018 |
On the optimization of the fixed‐stress splitting for Biot's equations
|
journal | June 2019 |
Similar Records
Generalized finite strains, generalized stresses, and a principle of virtual work for finite-element computer programs
Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO2 in Fractured Aquifers