ABSTRACT In this article, we propose and analyze a parallel finite element method for the closed‐loop geothermal system based on a fully overlapping domain decomposition approach. The fundamental idea is to decompose the spatial domain using a fully overlapping method and solve the equations independently in parallel over time on each processor via the backward Euler scheme. Within this framework, local solutions are obtained from a global grid that is computed with local refinement by each processor. By employing a series of parallel solvers, this algorithm ensures low communication costs and ease of implementation, resulting in significant computational savings and improved efficiency for the coupled system. This method exhibits its effectiveness in addressing problems that involve complex geometries and large scales. Some local a priori estimates, which are crucial for the theoretical analysis, are obtained. Finally, some numerical experiments are constructed to illustrate the effectiveness of our algorithms.
Li et al. (Thu,) studied this question.