P. N. Vabishchevich

A. E. Kolesov, P. N. Vabishchevich

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.

P. N. Vabishchevich

The work is devoted to the development and computational implementation of the homogenization method for modeling unsteady flows of a viscous incompressible fluid in periodic porous media taking into account memory effects. At the macrolevel, the flow is described by an integro-differential Darcy law with a tensor memory kernel determined by solving unsteady problems on the periodicity cell. The developed approach to computational homogenization is based on finding the steady-state and unsteady components of the conductivity tensor from solving auxiliary boundary value and spectral problems on the periodicity cell. The nonlocal macroscopic problem is transformed into a local system of differential equations by approximating the memory kernel as a sum of exponentials. Issues of spatial finite element approximation are discussed, and stable two-level schemes in time are constructed. The results of applying the developed computational homogenization technology for unsteady filtration problems in porous media to a two-dimensional test problem are presented.

P. N. Vabishchevich

Schemes with the second-order approximation in time are considered for numerical solving the Cauchy problem for an evolutionary equation of first order with a self-adjoint operator. The implicit two-level scheme based on the Padé polynomial approximation is unconditionally stable. It demonstrates good asymptotic properties in time and provides an adequate evolution in time for individual harmonics of the solution (has spectral mimetic stability). In fact, the only drawback of this scheme is the necessity to solve an equation with an operator polynomial of second degree at each time level. We consider modifications of these schemes, which are based on solving equations with operator polynomials of first degree. Such computational implementations occur, for example, if we apply the fully implicit two-level scheme (the backward Euler scheme). A three-level modification of the SM-stable scheme is proposed. Its unconditional stability is established in the corresponding norms. The emphasis is on the scheme, where the numerical algorithm involves two stages, namely, the backward Euler scheme of first order at the first (prediction) stage and the following correction of the approximate solution using a factorized operator. The SM-stability is established for the proposed scheme. To illustrate the theoretical results of the work, a model problem is solved numerically.