We study a finite element computational model for solving the coupled problem arising in the interaction between a free fluid and a fluid in a poroelastic medium. The free fluid is governed by the Stokes equations, while the flow in the poroelastic medium is modeled using the Biot poroelasticity system. Equilibrium and kinematic conditions are imposed on the interface. A mixed Darcy formulation is employed, resulting in continuity of flux condition of essential type. A Lagrange multiplier method is employed to impose weakly this condition. A stability and error analysis is performed for the semi-discrete continuous-in-time and the fully discrete formulations. A series of numerical experiments is presented to confirm the theoretical convergence rates and to study the applicability of the method to modeling physical phenomena and the sensitivity of the model with respect to its parameters.

We develop a multipoint stress mixed finite element method for linear elasticity with weak stress symmetry on quadrilateral grids, which can be reduced to a symmetric and positive definite cell centered system. The method is developed on simplicial grids in [4]. The method utilizes the lowest order Brezzi-Douglas-Marini finite element spaces for the stress and the trapezoidal quadrature rule in order to localize the interaction of degrees of freedom, which allows for local stress elimination around each vertex. We develop two variants of the method. The first uses a piecewise constant rotation and results in a cell-centered system for displacement and rotation. The second uses a continuous piecewise bilinear rotation and trapezoidal quadrature rule for the asymmetry bilinear form. This allows for further elimination of the rotation, resulting in a cell-centered system for the displacement only. Stability and error analysis is performed for both methods. First-order convergence is established for all variables in their natural norms. A duality argument is employed to prove second order superconvergence of the displacement at the cell centers. Numerical results are presented in confirmation of the theory.

No abstract

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.