Sequential symmetric interior penalty discontinuous Galerkin method for fully coupled quasi-static thermo-poroelasticity problems

In this paper, we investigate a sequentially decoupled numerical method for solving the fully coupled quasi-static thermo-poroelasticity problems with nonlinear convective transport. The symmetric interior penalty discontinuous Galerkin method is employed for spatial discretization and the backward Euler method for temporal discretization. Unlike other splitting algorithms, this type of sequential method does not require any internal iterations and the computational efficiency is higher than that of the fully implicit nonlinear numerical scheme. In the theoretical analysis, a cut-off operator is introduced to prove the existence and uniqueness of numerical solution and the stability analysis of numerical scheme is conducted. Then, we derive the optimal convergence order estimates in space and time. Finally, several numerical examples are presented to illustrate the accuracy and efficiency of our proposed method.