Splitting schemes with respect to physical processes for double-porosity poroelasticity problems

Abstract We consider unsteady poroelasticity problem in fractured porous medium within the classical Barenblatt double-porosity model. For numerical solution of double-porosity poroelasticity problems we construct splitting schemes with respect to physical processes, where transition to a new time level is associated with solving separate problem for the displacements and fluid pressures in pores and fractures. The stability of schemes is achieved by switching to three-level explicit-implicit difference scheme with some of the terms in the system of equations taken from the lower time level and by choosing a weight parameter used as a regularization parameter. The computational algorithm is based on the finite element approximation in space. The investigation of stability of splitting schemes is based on the general stability (well-posedness) theory of operator-difference schemes. A priori estimates for proposed splitting schemes and the standard two-level scheme are provided. The accuracy and stability of considered schemes are demonstrated by numerical experiments.

[1]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[2]  Francisco Armero,et al.  A new unconditionally stable fractional step method for non‐linear coupled thermomechanical problems , 1992 .

[3]  Mary F. Wheeler,et al.  Coupling multipoint flux mixed finite element methodswith continuous Galerkin methods for poroelasticity , 2013, Computational Geosciences.

[4]  D. Elsworth,et al.  Multiporosity/multipermeability approach to the simulation of naturally fractured reservoirs , 1993 .

[5]  W. Marsden.  I and J , 2012 .

[6]  A. Samarskii.  The Theory of Difference Schemes , 2001 .

[7]  Andro Mikelić,et al.  Convergence of iterative coupling for coupled flow and geomechanics , 2013, Computational Geosciences.

[8]  Ruben Juanes,et al.  A locally conservative finite element framework for the simulation of coupled flow and reservoir geomechanics , 2007 .

[9]  Dimitri E. Beskos,et al.  On the theory of consolidation with double porosity—II , 1982 .

[10]  S. G. Mikhlin,et al.  THE SPECTRUM OF A FAMILY OF OPERATORS IN THE THEORY OF ELASTICITY , 1973 .

[11]  Petr N. Vabishchevich,et al.  Splitting scheme for poroelasticity and thermoelasticity problems , 2013, FDM.

[12]  Anders Logg,et al.  DOLFIN: Automated finite element computing , 2010, TOMS.

[13]  Francisco Javier Lisbona,et al.  Finite-difference analysis of fully dynamic problems for saturated porous media , 2011, J. Comput. Appl. Math..

[14]  Petr N. Vabishchevich,et al.  Finite-Difference Analysis for the Linear Thermoporoelasticity Problem and Its Numerical Resolution by Multigrid Methods , 2012 .

[15]  Mary F. Wheeler,et al.  Overcoming the problem of locking in linear elasticity and poroelasticity: an heuristic approach , 2009 .

[16]  Petr N. Vabishchevich,et al.  Finite difference analysis of a double‐porosity consolidation model , 2012 .

[17]  Petr N. Vabishchevich,et al.  Splitting Scheme for Poroelasticity and Thermoelasticity Problems , 2014, FDM.

[18]  Ivan Yotov,et al.  Partitioning strategies for the interaction of a fluid with a poroelastic material based on a Nitsche’s coupling approach , 2014, 1403.5707.

[19]  Cornelis W. Oosterlee,et al.  An efficient multigrid solver for a reformulated version of the poroelasticity system , 2007 .

[20]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[21]  Vicente Hernández,et al.  SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems , 2005, TOMS.

[22]  G. Marchuk.  Splitting and alternating direction methods , 1990 .

[23]  Mary F. Wheeler,et al.  Iteratively coupled mixed and Galerkin finite element methods for poro‐elasticity , 2007 .

[24]  Jean E. Roberts,et al.  Mixed and hybrid finite element methods , 1987 .

[25]  Petr N. Vabishchevich,et al.  EXPLICIT-IMPLICIT SPLITTING SCHEMES FOR SOME SYSTEMS OF EVOLUTIONARY EQUATIONS , 2014 .

[26]  Ruben Juanes,et al.  Stability and convergence of sequential methods for coupled flow and geomechanics: Drained and undrained splits , 2011 .

[27]  Owe Axelsson,et al.  Stable discretization of poroelasticity problems and efficient preconditioners for arising saddle point type matrices , 2012, Comput. Vis. Sci..

[28]  Petr N. Vabishchevich,et al.  Additive Operator-Difference Schemes: Splitting Schemes , 2013 .

[29]  Petr N. Vabishchevich,et al.  A stabilized method for a secondary consolidation Biot's model , 2008 .

[30]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[31]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[32]  Anders Logg,et al.  The FEniCS Project Version 1.5 , 2015 .

[33]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[34]  Petr N. Vabishchevich,et al.  Difference Schemes with Operator Factors , 2002 .

[35]  P. Knabner,et al.  Numerical Methods for Elliptic and Parabolic Partial Differential Equations , 2003, Texts in Applied Mathematics.

[36]  H. Tchelepi,et al.  Stability and convergence of sequential methods for coupled flow and geomechanics: Fixed-stress and fixed-strain splits , 2011 .

[37]  P. Cochat,et al.  Et al , 2008, Archives de pediatrie : organe officiel de la Societe francaise de pediatrie.

[38]  Dimitri E. Beskos,et al.  On the theory of consolidation with double porosity—III A finite element formulation , 1984 .

[39]  Ivan Yotov,et al.  An operator splitting approach for the interaction between a fluid and a multilayered poroelastic structure , 2015 .

[40]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .