Auto-Stabilized Weak Galerkin Finite Element Methods for Biot's consolidation model on Non-Convex Polytopal Meshes

Abstract

This paper presents an auto-stabilized weak Galerkin (WG) finite element method for the Biot's consolidation model within the classical displacement-pressure two-field formulation. Unlike traditional WG approaches, the proposed scheme achieves numerical stability without the requirement of traditional stabilizers. Spatial discretization is performed using weak Galerkin finite elements for both displacement and pressure approximations, while a backward Euler scheme is employed for temporal discretization to ensure a fully implicit and stable formulation. We establish the well-posedness of the resulting linear system at each time step and provide a rigorous error analysis, deriving optimal-order convergence. A significant merit of this WG scheme is its flexibility on general shape-regular polytopal meshes, including those with non-convex geometries. By utilizing bubble functions as a primary analytical tool, the method produces stable, oscillation-free pressure approximations without specialized treatment. Numerical experiments are presented to validate the theoretical convergence rates and demonstrate the computational efficiency and robustness of the auto-stabilized formulation.