Hybridizable discontinuous galerkin methods for thermo-poroelastic systems

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 present 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.