Some semi-decoupled algorithms with optimal convergence for a four-field linear thermo-poroelastic model

We propose three semi-decoupled algorithms for efficiently solving a four-field thermoporoelastic model. The first two algorithms adopt a sequential strategy: at the initial time step, all variables are computed simultaneously using a monolithic solver; thereafter, the system is split into a mixed linear elasticity subproblem and a coupled pressure-temperature reaction-diffusion subproblem. The two variants differ in the order in which these subproblems are solved. To further improve computational efficiency, we introduce a parallel semidecoupled algorithm. In this approach, the four-field system is solved monolithically only at the first time step, and the two subproblems are then solved in parallel at subsequent time levels. None of the three algorithms requires iterative procedures at each time step, and are free from stabilization. Rigorous analysis confirms their unconditional stability, optimal convergence rates, and robustness under a wide range of physical parameter settings. These theoretical results are further validated by numerical experiments.