Cardiff University | Prifysgol Caerdydd ORCA
Online Research @ Cardiff 
WelshClear Cookie - decide language by browser settings

Thermal performance of borehole thermal energy storage for building energy systems in stratified ground with groundwater advection

Gao, Wu, Chu, Pengcheng, Qadrdan, Meysam ORCID: https://orcid.org/0000-0001-6167-2933, Tang, Fujiao, Xu, Wenjie, Li, Jinlong, Zhuang, Duanyang, Zhan, Liangtong and Chen, Yunmin 2026. Thermal performance of borehole thermal energy storage for building energy systems in stratified ground with groundwater advection. Energy and Buildings 369 (Part B) , 118053. 10.1016/j.enbuild.2026.118053

[thumbnail of 3-Manuscript-Clean version-FV.pdf]
Preview
PDF - Accepted Post-Print Version
Download (3MB) | Preview

Abstract

Borehole thermal energy storage (BTES) is widely used in ground source heat pump systems (GSHP) to support building heating and cooling through seasonal thermal storage. However, geological stratification and groundwater flow can strongly affect thermal plume migration, inter-borehole thermal interference, and heat pump operation. This study develops a composite analytical framework to predict BTES thermal performance in stratified ground with groundwater advection. The framework is validated against four independent experimental and numerical benchmarks, yielding RMSE, MAE, and MAPE values spanning 0.163–0.68 °C, 0.003–0.50 °C, and 0.012%–2.02%, respectively. It is then applied to a 90-day heat injection process in a 63-m-deep, three-layered ground system, under L-shaped, linear, and rectangular layouts, each with nine boreholes. Under the baseline condition of 4-m spacing and a groundwater velocity of 5 × 10−7 m/s, the rectangular array shows the least variation in field-averaged outlet temperature with groundwater flow direction (0.06 °C), versus 0.09 °C for the L-shaped and 0.40 °C for the linear arrays. Extended analyses of spacing and velocity show that this robustness is not universal. At 8-m spacing, the L-shaped array becomes comparable to, or slightly more stable than, the rectangular array, while the linear array remains the most sensitive to flow direction. The COP analysis further indicates a trade-off between efficiency and stability, with COP fluctuations of 0.009, 0.014, and 0.070 for the rectangular, L-shaped, and linear arrays, respectively. The vertical thermal discontinuity is quantitatively explained by Peclet numbers, which range from 3.96 to 7.00 in groundwater-bearing layers but fall to zero above the groundwater table. These findings enable more reliable design of BTES for GSHP systems in stratified ground with groundwater advection.

Item Type: Article
Date Type: Publication
Status: Published
Schools: Schools > Engineering
Additional Information: Rights retention policy applied
Publisher: Elsevier
ISSN: 0378-7788
Date of First Compliant Deposit: 14 September 2026
Date of Acceptance: 1 August 2026
Last Modified: 14 Sep 2026 14:30
URI: https://orca.cardiff.ac.uk/id/eprint/188879

Actions (repository staff only)

Edit Item Edit Item

Downloads

Downloads per month over past year

View more statistics