A partitioned solution procedure for geothermal reservoir analysis

The monolithic augmentation approach, when applied to the analysis of geothermal reservoirs, results in an unsymmetric coefficient matrix of the combined equations. The computational difficulties experienced with this approach suggest the use of a partitioned solution procedure. For this purpose, the fully coupled first-order three-field system of equations is integrated in the time domain by a two-point recurrence scheme. The resulting time-discrete algebraic system is then partitioned in such a way that the coefficient matrix remains symmetric. A stability analysis of the partitioned system is also carried out.