Splitting strategies for the fully-coupled nonlinear thermo-hydro-mechanical problem

Abstract

We propose novel semi-decoupled and fully-decoupled iterative algorithms for efficiently solving the fully-coupled nonlinear four-field thermo-poroelastic model discretized in space by discontinuous Galerkin method on polytopal grids. We present the model problem, its four-field formulation, and the arbitrary-order weighted symmetric interior penalty scheme exploited for its spatial discretization. Such a scheme is robust with respect to strong heterogeneities in the model coefficients. Then, we present the two solution strategies and prove that under suitable conditions both schemes are convergent. A wide set of numerical simulations is presented to assess the convergence and robustness properties of the proposed method. Moreover, we test the scheme with literature and physically sound test cases for proof-of-concept applications in the geophysical context.