Energy pile groups transmit geothermal energy and have attracted widespread attention as one of new building energy-saving technologies. Accurately predicting the time-dependent behaviors of energy pile groups is a challenge, given the complex thermal and mechanical interactions between piles, surrounding soils and the pile cap. This study presents a semi-analytical solution for analyzing energy pile groups within heat exchangers. Utilizing the transformed differential quadrature method, a flexible coefficient matrix for the saturated surrounding soils is acquired, which accounts for both consolidation and heat transfer. The piles are segmented, and the discrete solving equations considering thermal stresses and expansion are formulated. To accurately reflect the interactions among piles-to-piles, piles-to-soils and piles-to-pile cap, the coupled matrix equations are constructed with involving both the displacement coordination and the force equilibrium at the pile-soil interface as well as the pile cap. The validity of the proposed solution is confirmed through comparisons with results from onsite tests and simulations using COMSOL. Pivotal parameters including temperature variations, pile spacing, and the relative stiffness are discussed through examples. Compared with traditional simulation and field test, the proposed solution enables fast and accurate prediction of displacement and load distribution across pile groups, facilitating the safety evaluation of heat exchangers.

The current paper aims to experimentally investigate the thermal performance of geo-energy piles and walls fabricated with Phase Change heat exchangers. Four prototype concrete geo-energy structures (i.e., piles and walls) were tested using two distinct types of heat exchangers, including standard heat exchangers and PCM heat exchangers. The PCM heat exchangers utilized in the current study were filled up with two different types of Phase Change Materials (PCM) with melting points of 26 degrees C and 42 degrees C for geo-energy piles and walls, respectively. The thermal efficiency of the geo-energy piles/walls was experimentally assessed over 100 h of continuous operation under cycles of cooling and heating. The findings illustrated that using PCM heat exchangers led to enhancing the heat transfer efficiency of geo-energy piles by 75 % and 43 % in heating and cooling operations, respectively, compared to those achieved using a standard heat exchanger. Furthermore, the heat transfer performance of geo-energy walls with a PCM heat exchanger was enhanced by 43 % and 32 % in heating and cooling tests, respectively, compared to those achieved using a standard heat exchanger. Moreover, the findings indicated that the inclusion of PCM heat exchangers in geo-energy structures contributed to reducing: the impact on soil temperature and thermal interference radius as well as the potential structural damage due to thermal stress.

The performance of geothermal heat extraction in shallow aquifers depends on both Borehole Heat Exchanger (BHE) and soil or aquifer properties. In this work, an analysis of the thermal yield of a shallow geothermal reservoir was made numerically with the finite element method used to simulate heat and mass transfer in the three-dimensional reservoir. The main parameters for analysis which have been considered are the geometry and physical parameters of the BHE and grout, as well as aquifer matrix and groundwater fluid. Physical parameters are thermal conductivity, flow conductivity, expansion coefficient, porosity, volumetric heat capacity, anisotropy and dispersivity. The numerical tests have been performed in single BHE line source configuration representing numerically modelled thermal response test for the estimation of sustainable heat extraction. The domain size was a 100x100 meter rectangle with a depth of 200 meters. Three main lithological configurations have been modelled: gravel aquifer with low and high convection of groundwater fluid, as well as a shallow geothermal reservoir dominated by clay material without convection. For selected cases, the analysis for temporal and spatial discretization was also made. Three-dimensional transient modelling was made in FEFLOW (R) software with pre- and post-processing done in user-defined Python scripts. The results show the most influential parameters to be considered when setting up the real case simulation of geothermal heating and cooling, as well as optimal temporal and spatial discretization set-up with respect to expected thermal gradients in the reservoir.