A Primal Hybrid Finite Element Method for Poroelasticity

  • Correa, Maicon (University of São Paulo)
  • Taraschi, Giovanni (Unicamp)

Please login to view abstract download link

Poroelasticity problems combine the deformation of a porous matrix with the flow of a fluid through its porous. The most classical linear model for such a problem is the Biot system of partial differential equations, where the Darcy law is used to model the flow, linear elasticity is assumed to the poromechanics, and the Terzaghi principle is assumed to be valid [1]. In the context of Finite Element approximations for Biot, $H^1$-conforming and mixed strategies are widely discussed in the literature [2]. Hybrid methods, on the other hand, have not received the same level of attention. In this work, we discuss a Primal hybrid method for the Biot system, where Lagrange multipliers are added to weakly enforce the continuity of the pressure and displacement fields [3] and recently applied to the elasticity. The main advantage of the proposed strategy, in comparison with $H^1$-conforming methods, is the better mass/equilibrium conservation properties of the discrete solutions, which allows for an element-wise post-processing strategy to recover the stress and the velocity fields, leading to global $H$(div)-conforming approximations. The recovered stress can be either the effective or the total stress, depending on the choices made in the hybridization process. In this work, we present stable approximation spaces for the proposed method in quadrilateral meshes, allowing different approximation orders for the elasticity and the flow subproblems. We also discuss the factors that influence the $H$(div) convergence for the recovered velocity and stress. This work is supported by National Council for Scientific and Technological Development - CNPq (grants 403673/2025-9 and 307679/2023-3) and the São Paulo Research Foundation - FAPESP (grant 2013/07375-0). [1] Correa M.R., Murad M.A., Fixed-stress sequential schemes for a black-oil model in poroelastic media, J. Comput. Phys., vol .543, Art. 114406, 2025. [2] M. A. Murad and A. F. Loula. Improved accuracy in finite element analysis of Biot's consolidation problem. Comput. Methods Appl. Mech. Eng. , vol. 95, n. 3, pp. 359–382, 1992. [3] P-A. Raviart and J-M. Thomas. Primal hybrid finite element methods for 2nd order elliptic equations. Math. Comput, Vol. 31, pp. 391-413, 1977.