No abstract

We develop non-overlapping domain decomposition methods for the Biot system of poroelasticity in a mixed form. The solid deformation is modeled with a mixed three-field formulation with weak stress symmetry. The fluid flow is modeled with a mixed Darcy formulation. We introduce displacement and pressure Lagrange multipliers on the subdomain interfaces to impose weakly continuity of normal stress and normal velocity, respectively. The global problem is reduced to an interface problem for the Lagrange multipliers, which is solved by a Krylov space iterative method. We study both monolithic and split methods. In the monolithic method, a coupled displacement-pressure interface problem is solved, with each iteration requiring the solution of local Biot problems. We show that the resulting interface operator is positive definite and analyze the convergence of the iteration. We further study drained split and fixed stress Biot splittings, in which case we solve separate interface problems requiring elasticity and Darcy solves. We analyze the stability of the split formulations. Numerical experiments are presented to illustrate the convergence of the domain decomposition methods and compare their accuracy and efficiency.

We develop a space-time mortar mixed finite element method for parabolic problems. The domain is decomposed into a union of subdomains discretized with non-matching spatial grids and asynchronous time steps. The method is based on a space-time variational formulation that couples mixed finite elements in space with discontinuous Galerkin in time. Continuity of flux (mass conservation) across space-time interfaces is imposed via a coarse-scale space-time mortar variable that approximates the primary variable. Uniqueness, existence, and stability, as well as a priori error estimates for the spatial and temporal errors are established. A space-time non-overlapping domain decomposition method is developed that reduces the global problem to a space-time coarse-scale mortar interface problem. Each interface iteration involves solving in parallel space-time subdomain problems. The spectral properties of the interface operator and the convergence of the interface iteration are analyzed. Numerical experiments are provided that illustrate the theoretical results and the flexibility of the method for modeling problems with features that are localized in space and time.

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-intime 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.

We develop a space-time mortar mixed finite element method for parabolic problems. The domain is decomposed into a union of subdomains discretized with non-matching spatial grids and asynchronous time steps. The method is based on a space-time variational formulation that couples mixed finite elements in space with discontinuous Galerkin in time. Continuity of flux (mass conservation) across space-time interfaces is imposed via a coarse-scale space-time mortar variable that approximates the primary variable. Uniqueness, existence, and stability, as well as a priori error estimates for the spatial and temporal errors are established. A space-time non-overlapping domain decomposition method is developed that reduces the global problem to a space-time coarse-scale mortar interface problem. Each interface iteration involves solving in parallel space-time subdomain problems. The spectral properties of the interface operator and the convergence of the interface iteration are analyzed. Numerical experiments are provided that illustrate the theoretical results and the flexibility of the method for modeling problems with features that are localized in space and time.