In this paper, an adapted model is developed for borehole heat exchangers (BHEs) to simulate geothermal applications such as heat storage on a large scale efficiently and with high accuracy. The adapted numerical model represents all BHE components, allowing for a detailed representation of the governing processes. The approach is calibrated and validated for a single U-tube BHE using a high-resolution experimental data set from a laboratory thermal response test. It is found that the computational effort can be reduced by factors of *50, *50 and *25 for single U-tube, double U-tube and coaxial BHEs, respectively, if an absolute deviation of less than 1 % compared to a conventional fully discretised model is allowed. Computation times can be reduced further by accepting higher deviations. The adapted modelling approach allows for a detailed and correct representation of the temporal and spatial temperature distribution under highly transient conditions by applying it to a high-temperature heat storage scenario using multiple BHEs. The model is especially suited to represent coupled flow and heat transport processes, to account for groundwater flow in the BHE region as well as geological heterogeneities and especially interaction between a large number of BHEs.Keywords Borehole thermal energy storage Á Borehole heat exchanger Á Numerical simulation Á Fully discretised models Á OpenGeoSys
List of symbols aSolid compressibility (Pa -1 ) bFluid compressibility (Pa -1 ) cq Volumetric heat capacity (J m -3 K -1 ) dThickness (m) D Heat diffusion dispersion tensor g Gravitational acceleration (m s -1 ) k Intrinsic permeability (m) L Length (m) m Side length (m) n Porosity (-) p Pressure (Pa) Q Sources and sinks (kg m -3 s -1 ) Q H Heat sources and sinks (W m -3 ) rRadius (m) R th Thermal resistance (K W -1 )Transport velocity (m s -1 ) z Depth (m) a Dispersivity (m) kThermal conductivity (W m -1 K -1 ) l Fluid dynamic viscosity (N s m -2 ) q Density (kg m -3 )