Hybridizable Discontinuous Galerkin Methods for Thermo-Poroelastic Systems

Abstract: We propose a high-order hybridizable discontinuous Galerkin (HDG) formulation for the fully dynamic, linear thermo-poroelasticity problem. The governing equations are formulated as a first-order hyperbolic system incorporating solid and fluid velocities, heat flux, effective stress, pore pressure, and temperature as state variables. We establish well-posedness of the continuous problem using semigroup theory and develop an energy-consistent HDG discretization. The method exploits computational advantages of HDG-including locality and static condensation-while maintaining energy conservation for the coupled system. We establish an $hp$-convergence analysis and support it with comprehensive numerical experiments, confirming the theoretical rates and showcasing the method's effectiveness for thermo-poroelastic wave propagation in heterogeneous media.