A priori and a posteriori error bounds for the fully mixed FEM formulation of poroelasticity with stress-dependent permeability

Arbaz Khan,Bishnu P. Lamichhane,Ricardo Ruiz-Baier,Segundo Villa-Fuentes

2024-09-05

Abstract:We develop a family of mixed finite element methods for a model of nonlinear poroelasticity where, thanks to a rewriting of the constitutive equations, the permeability depends on the total poroelastic stress and on the fluid pressure and therefore we can use the Hellinger--Reissner principle with weakly imposed stress symmetry for Biot's equations. The problem is adequately structured into a coupled system consisting of one saddle-point formulation, one linearised perturbed saddle-point formulation, and two off-diagonal perturbations. This system's unique solvability requires assumptions on regularity and Lipschitz continuity of the inverse permeability, and the analysis follows fixed-point arguments and the Babuška--Brezzi theory. The discrete problem is shown uniquely solvable by applying similar fixed-point and saddle-point techniques as for the continuous case. The method is based on the classical PEERS$_k$ elements, it is exactly momentum and mass conservative, and it is robust with respect to the nearly incompressible as well as vanishing storativity limits. We derive a priori error estimates, we also propose fully computable residual-based a posteriori error indicators, and show that they are reliable and efficient with respect to the natural norms, and robust in the limit of near incompressibility. These a posteriori error estimates are used to drive adaptive mesh refinement. The theoretical analysis is supported and illustrated by several numerical examples in 2D and 3D.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to model the interaction between fluid flow and the mechanical response of materials in nonlinear poroelastic materials. Specifically, the research focuses on developing a Mixed Finite Element Method (MFEM) for solving a nonlinear poroelastic model in which the permeability depends on the total poroelastic stress and fluid pressure. Since this dependence of permeability is non - monotonic, it is difficult to directly apply traditional tools for monotonic saddle - point problems. Therefore, the main objectives of this paper are:

1. **Propose a new Mixed Finite Element Method**: By reformulating the constitutive equation of permeability so that it depends on the total poroelastic stress and fluid pressure, the Hellinger - Reissner principle can be used to avoid explicitly solving for strain. This method is robust when dealing with near - incompressibility and when the storage coefficient is close to zero.

2. **Analyze the well - posedness of the model**: Through fixed - point theory and Babuška - Brezzi theory, the unique solvability of the continuous and discrete problems is proven. In particular, the article considers the Lipschitz continuity and positive - definiteness conditions of the inverse of permeability in the \(L^2\) space, which are crucial for ensuring the well - posedness of the problem.

3. **Derive a priori error estimates**: Through theoretical analysis, a priori error estimates are obtained, and these estimates are still valid in the case of near - incompressibility and when the storage coefficient is close to zero.

4. **Propose a posteriori error estimators**: Fully computable residual - based a posteriori error estimators are developed, and it is proven that they are reliable and efficient in the natural norm and remain robust in the near - incompressibility limit. These a posteriori error estimators can be used to drive adaptive mesh refinement.

5. **Numerical verification**: Through multiple two - dimensional and three - dimensional numerical experiments, the effectiveness of the theoretical analysis is verified, and the performance of the proposed Mixed Finite Element Method in practical applications is demonstrated.

In summary, this paper aims to solve the numerical simulation challenges of fluid - structure interaction problems in nonlinear poroelastic materials by developing and analyzing a new Mixed Finite Element Method, especially when the permeability depends on stress and pressure.