We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling. This poroelasticity problem suffers from rapidly oscillating material parameters, which calls for a thorough numerical treatment. In this paper, we propose a method based on the local orthogonal decomposition technique and motivated by a similar approach used for linear thermoelasticity. Therein, local corrector problems are constructed in line with the static equations, whereas we propose to consider the full system. This allows to benefit from the given saddle point structure and results in two decoupled corrector problems for the displacement and the pressure. We prove the optimal first-order convergence of this method and verify the result by numerical experiments.

In this work, we employ the Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM) to solve the problem of linear heterogeneous poroelasticity with coefficients of high contrast. The proposed method makes use of the idea of energy minimization with suitable constraints in order to generate efficient basis functions for the displacement and the pressure. These basis functions are constructed by solving a class of local auxiliary optimization problems based on eigenfunctions containing local information on the heterogeneity. Techniques of oversampling are adapted to enhance the computational performance. Convergence of first order is shown and illustrated by a number of numerical tests.

No abstract

In this paper, we consider a pressure-velocity formulation of the heterogeneous wave equation and employ the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) to solve this problem. The proposed method provides a flexible framework to construct crucial multiscale basis functions for approximating the pressure and velocity. These basis functions are constructed by solving a class of local auxiliary optimization problems over the eigenspaces that contain local information on the heterogeneity. Techniques of oversampling are adapted to enhance the computational performance. The first-order convergence of the proposed method is proved and illustrated by several numerical tests.

No abstract

In this research, an online basis enrichment strategy for the constraint energy minimizing generalized multiscale finite element method in mixed formulation is proposed. The online approach is based on the technique of oversampling. One makes use of the information of residual and the data in the partial differential equation such as the source function. The analysis presented shows that the proposed online enrichment leads to a fast convergence from multiscale approximation to the fine-scale solution. The error reduction can be made sufficiently large by suitably selecting oversampling regions and the number of oversampling layers. Also, the convergence rate of the enrichment can be tuned by a user-defined parameter. Numerical results are provided to illustrate the efficiency of the proposed method.

No abstract

No abstract