Multiscale mortar mixed finite element methods for the Biot system of poroelasticity

TL;DR

This paper introduces a multiscale mortar mixed finite element method for the Biot system of poroelasticity, enabling efficient and accurate simulations on non-matching grids with reduced computational complexity.

Contribution

It develops a novel multiscale mortar domain decomposition approach with theoretical analysis and efficient algorithms for the Biot system of poroelasticity.

Findings

The method achieves stability and error estimates under certain mortar space conditions.

Numerical experiments confirm the multiscale method's accuracy and efficiency.

The approach reduces computational cost by decoupling subdomain solves from iteration counts.

Abstract

We develop a mixed finite element domain decomposition method on non-matching grids for the Biot system of poroelasticity. A displacement-pressure vector mortar function is introduced on the interfaces and utilized as a Lagrange multiplier to impose weakly continuity of normal stress and normal velocity. The mortar space can be on a coarse scale, resulting in a multiscale approximation. We establish existence, uniqueness, stability, and error estimates for the semidiscrete continuous-in-time formulation under a suitable condition on the richness of the mortar space. We further consider a fully-discrete method based on the backward Euler time discretization and show that the solution of the algebraic system at each time step can be reduced to solving a positive definite interface problem for the composite mortar variable. A multiscale stress-flux basis is constructed, which makes the number of subdomain solves independent of the number of iterations required for the interface problem, as well as the number of time steps. We present numerical experiments verifying the theoretical results and illustrating the multiscale capabilities of the method for a heterogeneous benchmark problem.