A computational macroscale model for the time fractional poroelasticity problem in fractured and heterogeneous media

Abstract

In this paper, we consider the poroelasticity problem with a time memory formalism that couples the pressure and displacement, and we assume this multiphysics process occurs in multicontinuum media. A coupled system of equations for pressures in each continuum and elasticity equations for displacements of the medium are included in the mathematical model. Based on the Caputo's time fractional derivative, we derive an implicit finite difference approximation for time discretization. Also, a Discrete Fracture Model (DFM) is used to model fluid flow through fractures and treat the complex network of fractures. Further, we develop a coarse grid approximation based on the Generalized Multiscale Finite Element Method (GMsFEM), where we solve local spectral problems for construction of the multiscale basis functions. The main idea of the proposed method is to reduce the dimensionality of the problem because our model equation has multiple fractional powers, there multiple unknowns with memory effects. Consequently, the solution is on a coarse grid, which saves some computational time. We present numerical results for the two-dimensional model problems in fractured heterogeneous porous media. After, we investigate error analysis between reference (fine-scale) solution and multiscale solution with different numbers of multiscale basis functions. The results show that on a coarse grid, the proposed approach can achieve good accuracy.