Circular Image

P.J. Vardon

info

Please Note

175 records found

Geothermal energy has the potential to decarbonize heating, cooling, and power production. However, managing the efficient and sustainable exploitation of geothermal resources is challenging due to the limited data availability, which restricts our ability to characterize and quantify the multi-scale, hierarchical geological structures of the hosting reservoirs. In this study, we propose a scenario-based data assimilation framework that enables the efficient modelling of multiple complex geological scenarios and is linked to flow and heat transfer simulations for subsequent uncertainty analysis. This framework is based on an ensemble smoother with multiple data assimilation (ESMDA) and demonstrated on a channelized fluvial geothermal reservoir. By improving the open-source Rapid Reservoir Modelling (RRM) tool, we efficiently create multiple deterministic fluvial geothermal reservoir scenarios that honors facies along well paths in a probabilistic manner by randomly selecting, cropping, and stacking channelized layers from the layer template library. Petrophysical properties for each scenario are then modelled using geostatistics to generate a geologically plausible and sufficiently diverse ensemble of reservoir realizations. The multiple scenarios and corresponding ensemble realizations are then subjected to heat and fluid flow simulations using the open-source Delft Advanced Research Terra Simulator (open-DARTS) to quantify the uncertainty of production temperatures and reservoir pressures. Finally, ESMDA is employed to assimilate temperature and pressure profiles at the injection well, monitoring borehole, and production well across all members of the ensemble realizations for the different geological scenarios. We demonstrate the applicability of our framework using a synthetic, yet geologically consistent, case study of a low-enthalpy geothermal system where heat is produced from a geothermal doublet located in a channelized fluvial sandstone reservoir. The framework enables the falsification of geological scenarios with poor data assimilation performance that is unlikely to reflect the actual reservoir architecture, and supports the identification of plausible geological scenarios that are more likely to represent the subsurface geology based on the deviation of modelled and observed well temperature and pressure profiles. The workflow offers an efficient way to constrain geological uncertainties inherent to geologically complex geothermal reservoirs and improve the forecasting of production temperatures and pressure differences. ...
High-temperature aquifer thermal energy storage (HT-ATES) can play a key role in the energy transition. For well completion of conventional low-temperature ATES and groundwater wells, grout and/or clay pellets are typically utilised as annular materials to ensure the long-term well integrity. It is not yet known if such materials can also be used in HT-ATES working conditions. In this work, a novel approach to evaluate the sealing performance for such completion materials is proposed and tested over multiple thermal heating and cooling cycles representative of the conditions of HT-ATES operation. The experimental framework utilises a novel experimental design to test the apparent transmissivity of the annular material, followed by micro-CT scanning. During each test, up to 11 thermal cycles are applied, with temperature variations between 22oC and 90oC. For grouts after 7 days of curing, micro-CT scans reveal debonding and the occurrence of micro-annuli with an equivalent diameter of approximately 26% of the original cross-section. After 28 days of curing, the thermal cycles had a much reduced impact on micro-annulus formation. The corresponding apparent transmissivity decreased up to 80% for samples containing a high percentage of cementitious minerals and a low water-to-grout ratio. The clay pellets, saturated with fresh water, demonstrated effective sealing capacity and an impermeable behaviour. However, clay pellets saturated with 0.25 mol/L NaCl, showed up to an 85% decrease in swelling capacity yet still exhibited impermeable behaviour. The results indicated that thermal cycles affect the integrity of grouts, while clay pellets show resilience to them. Furthermore, longer curing periods and specific chemical compositions improve sealing performance and provide resilience to thermal cycles. ...
Journal article (2026) - David Geerts, Wen Liu, Alexandros Daniilidis, Philip J. Vardon, Gert Jan Kramer
Integrating renewable energy into district heating creates a heat supply–demand mismatch that High-Temperature Aquifer Thermal Energy Storage (HT-ATES) can help address. However, the potential greenhouse gas emission reduction and financial benefits of HT-ATES have received limited attention. Additionally, the interplay between the demand, supply components, and HT-ATES has been overlooked, while the assessment of integrating HT-ATES into a district heating system is crucial to understanding the benefits of the HT-ATES implementation. This study evaluates the integration of HT-ATES into a district heating system, focusing on both economic and environmental performance indicators. It novelly accounts for the dynamic operational interactions between HT-ATES and other system components, enabling a more realistic assessment of operational choices. The model is applied to a case study of a simplified district heating system. The results show that the relative size of the heat supplier compared to heat demand is a key determinant of the cost-effectiveness of HT-ATES. In the case study, a geothermal doublet reduced the levelized cost of heat by 25–37 €/MWh compared to a gas boiler, while also reducing reliance on fossil fuels. In contrast, HT-ATES had a limited impact on total system costs, regardless of whether it operated when stored heat was available or was used for peak shaving. Nevertheless, HT-ATES increased the renewable energy share by 9%–18% across all scenarios. Furthermore, the optimal geothermal capacity differed depending on whether HT-ATES was included. Finally, while a high renewable energy share can be cost-effective, achieving 100% renewable heat was found to be highly cost-ineffective in this case. These results support informed decision-making on HT-ATES implementation under appropriate system design conditions. ...

A mechanistic approach to understand its water resistance behaviour

Cow-dung is a widely used stabiliser applied in traditional earthen buildings with one objective to improve water resistance. However, most research has focused on explaining its mechanical strength, with only one study suggesting water resistance mechanism via formation of insoluble compounds at high pH, a phenomenon uncommon in natural cow dung and soil mixtures. This article investigates the water-resistance behaviour of cow-dung stabilised compressed earthen blocks (CD-CEBs) through an extensive experimental programme to understand the influence of cow-dung and soil related factors and to characterise the components of cow-dung responsible for its water resistance. It was found that the small-sized microbial aggregates (SSMA) present in cow-dung, which are negatively charged hydrophobic aggregates of low specific surface area, are responsible for enhanced water resistance of CD-CEBs. The insights gained from experiments are compiled to recommend the following strategies for improved performance of CD-CEBs: (i) The use of wet cow-dung is advised over dry cow-dung as it provided over 80 times better water resistance; (ii) Adopting a higher compaction liquid content (by 3%) improved the water resistance by over 40 times; (iii) The water resistance of CD-CEBs was improved over 30 times by using soils rich in low-swelling clay minerals such as kaolinite. A case study applying these findings demonstrates the successful scaleup from the lab to field showcasing potential of cow-dung and soil in low-carbon construction. ...
Despite the advantages of using Bayesian networks for probabilistic risk assessment, adoption in practice has been limited due to the lack of realistic, facility-scale studies. Scaling up from systems to facility-level safety assessments poses challenges in (i) integrating external hazards and their cascading effects, and (ii) resolving non-homogeneity of various technical and human reliability models. The novelty of the study is in formalising risk integration using Bayesian networks, at facility scale, and demonstrating its effectiveness in addressing associated challenges. A Bayesian network-based multi-hazard risk framework is introduced and demonstrated for a nuclear power plant subject to flooding and earthquake hazards, capturing dependencies among hazards and consequences. Individual reliability models – conventionally extraneous to facility-wide risk models – are included as subnetworks by using Bayesian network-based surrogate models for technical systems and a Bayesian networks approach for human reliability modelling. Two approaches are used for subnetwork integration – object-oriented and unified Bayesian networks. The unified approach allows for prediction, diagnostics and inter-causal reasoning since Bayesian inference is bi-directional. Conversely, in the object-oriented approach, diagnostics are limited to within individual subnetworks and as a consequence the model can potentially neglect dependencies between objects. However, the object-oriented model requires only 50 % of the computational memory and consumes less than 25% of the runtime as the unified network, while improving visual clarity of the risk model. The model reveals key insights – for example, variations in operator stress or available response time during a hazard event can result in up to a 77 % change in top event probability – demonstrating its effectiveness in capturing critical relationships in complex, facility-scale risk scenarios. These findings can be used to suitably allocate resources towards risk mitigation and plant safety management. ...
Soft stimulation technologies have been proposed as a means to reduce the breakdown pressure and mitigate the risk of induced seismicity during geothermal reservoir stimulation. Yet, the underlying mechanisms remain poorly understood due to the complexity of the coupled thermo-hydro-mechanical (THM) processes. In this work, a fully coupled THM model is developed to evaluate and compare the performance of different stimulation scenarios (monotonic, stepwise injection rate, cyclic injection rate or temperature, and stepwise combined with cyclic injection rate stimulation) on a synthetic, highly permeable reservoir with near-borehole clogging. Simulation results show that stepwise injection rate stimulation yields the most favourable outcomes, followed by the stepwise injection rate combined with cyclic injection rate stimulation. On the other hand, fatigue effects are seen to play a negligible role in the improved performance since the tensile stress at the fracture tip is relaxed with the continuous fracture growth. In addition, cyclic injection temperature stimulation is generally neither better nor worse than monotonic stimulation, but has slightly different characteristics, creating more local damage controlled by the period of the injection cycle. Cyclic injection rate stimulation can slightly reduce the peak pressure, compared with monotonic stimulation, but only when the injection rate is low. The reduction in peak pressure occurs due to the combination of thermally-induced stresses associated with cooling and incremental damage rather than any influence of fatigue. Stepwise or low-frequency cyclic injection rate stimulation are suggested rather than a high-frequency cyclic injection rate stimulation, while injection with cyclic temperatures is suggested when more local damage is wanted. ...
Journal article (2026) - Wen Luo, Anne Catherine Dieudonné, Josselin Ouf, Florian Amann, Philip J. Vardon
Rocks can undergo fatigue failure when subjected to cyclic mechanical, hydraulic, or thermal loadings, or a combination of these. Therefore, accounting for possible fatigue damage is important for subsurface engineering projects, such as the cyclic stimulation of geothermal reservoirs. However, existing models do not simultaneously account for degradation of both tensile strength and stiffness under varying-amplitude loading and coupled thermo-hydro-mechanical (THM) conditions. To address this, a new cohesive zone model is developed to account for the effect of fatigue on tensile strength and stiffness. The model is then used within the framework of zero-thickness interface elements to simulate the response of pre-existing or new fractures. Hydraulic and thermal processes are included in both the cohesive interface elements and the continuum elements, allowing the consideration of coupled thermo-hydro-mechanical processes. The fatigue damage variable is set to evolve with the number and magnitude of cycles according to Palmgren-Miner's rule. The proposed method is validated against three laboratory tests from the literature, including cyclic Brazilian test, cyclic hydraulic fracturing test and cyclic thermal stimulation test. All three validation results show that the fatigue damage or reduced breakdown pressure can be well reproduced. Mesh sensitivity based on the simulation of the Brazilian test, in which interface elements are inserted in-between all the continuum elements, highlights the influence of the mesh orientation and mesh density on the simulation results. In addition, stabilisation of the method is demonstrated by increasing the mechanical viscosity, which must be used with care to avoid predicting a longer fatigue life. The ability of the method to handle varying-amplitude cyclic loading is demonstrated by the simulation of a synthetic cyclic loading scheme based on the Brazilian test. The proposed method can be used to support the design of cyclic thermal stimulation campaigns for geothermal (or other) reservoirs, by being able to simulate the reduction in strength due to fatigue, and thus reducing stimulation pressures needed. ...
Journal article (2026) - Mouadh Rafai, Merita Tafili, Yuepeng Dong, Philip J. Vardon
Realistic numerical modeling of energy piles in soft soil requires advanced constitutive relationships capable of capturing the inherent thermo-plastic behavior of the surrounding ground. In this study, a newly developed rate-dependent, thermo-plastic constitutive model, called AVISA-T, is employed within the Plaxis Finite Element (FE) code to simulate the response of a well-instrumented energy pile embedded in multilayered soft soils subjected to thermo-mechanical loading. Following material parameter calibration and model prediction validation using non-isothermal laboratory tests on the soils surrounding the pile, the model was employed to simulate full-scale in-situ tests. In these simulations, the pile was initially subjected to either 0 % or 60 % of its bearing capacity and then exposed to continuous cooling over a period of up to three months. The AVISA-T model effectively reproduces the development of contractive and expansive strains, as well as compressive and tensile stresses that coexist along the pile shaft, including the accumulation of residual strains and stresses. In the absence of axial mechanical load, both residual contractive and expansive strains were observed, accompanied by irreversible uplift of the pile head, primarily attributed to non-uniform, unrecovered temperature changes. Moreover, under higher mechanical loading, the model captures dragdown effects resulting from thermal shrinkage of the surrounding soil, which contributes to the accumulation of permanent strains, stresses, and settlements. A comparison between simulations using the common Modified Cam Clay (MCC) model, the AVISA model without thermal effects and the AVISA-T model highlights the importance of using models including thermal plasticity for engineering practice. ...
Journal article (2026) - David Bruhn, Hemmo A. Abels, Patrick Fulton, Virginie Harcouët-Menou Harcouët-Menou, Ernst Huenges, Stefan Jansen, Alexis Koulidis, Susanne Laumann, Haiyan Lei, Joseph Moore, Paula Rulff, Thorben Schöfisch, Auke Barnhoorn, Evert Slob, Philip J. Vardon, Liliana Vargas Meleza, Denis Voskov, Claire Bossennec, Aoife K. Braiden, Maren Brehme, Romain Chassagne, Alexandros Daniilidis, Mathieu Darnet, Guy Drijkoningen
Low-enthalpy geothermal heat production is becoming increasingly common, which leads to the potentially competitive use of the available subsurface space, especially in densely populated urban areas. A specific challenge presented by the high density of different geothermal systems is understanding the details of convective and conductive heat flow processes and detailed monitoring of properties and processes in the subsurface.

On the TU Delft campus, we aim to drill a borehole of around 4.5 km depth to be used for the exploration, observation, and monitoring of subsurface processes that will be part of a larger research infrastructure under development. This so-called urban energy laboratory includes – in addition to the deep multi-use borehole – a well-instrumented geothermal doublet drilled in 2023, reaching to a depth of 2.2 km; a local seismic monitoring system (installed in 2022); an ultra-sensitive portable seismic monitoring array; and a high-temperature aquifer heat storage system (HT-ATES), for which a pilot well was drilled in 2024. With this urban energy laboratory, we want to tackle problems and better understand processes related to multiple and/or competing subsurface uses in urban environments. The deep exploration and monitoring borehole is designed specifically to monitor fluid and/or flux movement in 3D with unprecedented precision, aiming to understand the propagation of the geothermal cold front and reservoir pressures.

During the 3 d International Continental Scientific Drilling Program (ICDP)-sponsored UrbEnLab workshop, 75 scientists from 17 countries met in Delft, the Netherlands, in June 2024 to prioritize the scientific ambitions of the deep exploration and monitoring borehole and to discuss potential techniques that could be applied to tackle them. Assessing the life cycle of a geothermal system situated in a complex heterogeneous sedimentary system was defined as the broad aim, with revealing the detailed flow field established being a key priority. ...

Monitoring Borehole DEL-HTO-P01

Report (2026) - Alexis Koulidis, Martin van der Schans, Philip J. Vardon, Martin Bloemendal
This document describes the drilling, completion and testing of the pilot borehole (DEL-HTO-P01), which will later be utilized as a monitoring well for High Temperature Aquifer Thermal Energy Storage (HT-ATES). The work includes all materials, components, tools, and services that are related to the preparation, delivery, installation, measurement, and completion of the components belonging to the test drilling and monitoring well. The following analysis combines various data sources to characterize the subsurface for the Delft Demonstration site (WP1). The location of the DEL-HTO-P01 is demonstrated in Figure 1, which is located at RD-coordinates X = 85,257 m, Y = 445,757 m.

All data presented in this report have been published via TU Delft institutional data repository 4TU.ResearchData under the CC BY 4.0 license. The initial borehole dataset and the CT scan data can be found via https://doi.org/10.4121/1694ba82-db41-4017-8d1c-9de3ce1a785e and https://doi.org/10.4121/9902ebc7-dbd6-43b1-a3c3-85ff9ab645a9. ...
Journal article (2026) - M. Rafai, M. Tafili, Y. Dong, P. J. Vardon
In this study, a newly developed rate-dependent thermo-plastic constitutive model was enhanced to incorporate thermally accelerated creep and implemented into the Plaxis finite-element code, enabling the simulation of the behaviour of a well-instrumented energy pile in multilayered soft soils under thermomechanical loads. First, the model was validated against non-isothermal laboratory tests on soils surrounding the pile, and then against simulations of field tests. The results revealed that the inclusion of thermally accelerated creep improves the prediction of irreversible pile settlement, which is primarily attributed to the accumulation of volumetric contraction in the surrounding soil after each thermal cycle. The analysis also distinguishes between drag down effects resulting from thermo-elastic and thermo-plastic behaviour, as well as those induced by long-term creep. ...
Journal article (2026) - Merita Tafili, Mouadh Rafai, Mohammad Javad Salimi, Mohammadsadegh Ashrafi, Phil Vardon, Torsten Wichtmann
Energiepfähle werden zunehmend als multifunktionale Gründungselemente eingesetzt, die tragende Funktionen mit der Nutzung geothermischer Energie kombinieren. Während das mechanische Verhalten in der Bemessung üblicherweise berücksichtigt wird, werden thermische Einwirkungen häufig in vereinfachter elastischer Form behandelt. Langzeitiges Heizen und Kühlen kann jedoch irreversible Bodenverformungen und Spannungsumlagerungen hervorrufen, insbesondere in weichen und geschichteten Böden. Dieser Beitrag untersucht die Relevanz thermoplastischen Bodenverhaltens für die Bemessung von Energiepfählen unter Verwendung eines thermo-visko-hypoplastischen Stoffmodells für feinkörnige Böden (AVISA-T), das zyklische Belastung berücksichtigt. Das Modell wird in einem Finite-Elemente-Programm implementiert und anhand großmaßstäblicher Energiepfahlversuche in Delft (NL) unter unterschiedlichen axialen Lastniveaus und langzeitiger thermischer Beanspruchung evaluiert. Die Ergebnisse zeigen, dass thermoplastische Effekte maßgeblich die Pfahlkopfverformungen sowie die Umlagerung der Normalkräfte entlang des Pfahls bestimmen können, was mit konventionellen thermo-elastischen oder thermo-elasto-plastischen Ansätzen nicht erfasst wird. ...
Conference paper (2025) - G. Song, S. Geiger, D. Voskov, H. Abels, P. Vardon
Long-term geothermal production is subject to considerable uncertainty due to limited data availability and inherent geological heterogeneity. While observation and data acquisition improve our understanding of the reservoir, they also contribute significantly to project costs. It is essential to identify the most informative observation strategy. In this study, we apply a previously developed scenario-based data assimilation framework that integrates rapid geological modelling, efficient numerical simulation, and Ensemble Smoother with Multiple Data Assimilation (ESMDA) to constrain uncertainties in reservoir properties and production forecasts to a synthetic but geologically realistic fluvial geothermal system and conduct a data worth analysis to evaluate the impact of different observations (production temperature and injection pressure, well temperature and pressure profiles, etc.) on uncertainty reduction. Results show that production temperature and injection pressure alone, though cost-effective, are insufficient to significantly reduce uncertainties in reservoir performance forecasts. In contrast, well temperature and pressure profiles exhibit substantially higher data worth, leading to much better-constrained predictions. Moreover, incorporating a monitoring borehole further constrains uncertainty by capturing subsurface dynamics between the injector and producer. These findings underscore the importance of monitoring pressure and temperature profiles in the wells of a geothermal doublet. ...

A conditional Generative Adversarial Network for geotechnical subsurface schematisation

Journal article (2025) - F. A. Campos Montero, B. Zuada Coelho, E. Smyrniou, R. Taormina, P. J. Vardon
Subsurface schematisations are a crucial geotechnical problem which generally consists of filling substantial gaps in subsurface information from the limited site investigation data available and relying heavily on the engineer’s experience and occasionally geostatistical tools. To address this, schemaGAN, a conditional Generative Adversarial Network (GAN) to generate geotechnical subsurface schematisations from site investigation data is introduced. This novel method can learn complex underlying rules that govern the subsurface geometries and anisotropy from a big database of training cross-sections, and can produce subsurface schematisations from Cone Penetration Tests (CPT) in an insignificant timeframe. To test and demonstrate the performance of schemaGAN, a database of 24,000 synthetic geotechnical cross-sections with their corresponding CPT data was created, including spatial variability and gradually spatially varying layers. After training, the effectiveness of schemaGAN was compared against several interpolation methods, and it is seen that schemaGAN outperforms all other methods, with results characterised by clear layer boundaries and an accurate representation of anisotropy within the layers. SchemaGAN’s superior performance was confirmed through a blind survey, and in two real case studies in the Netherlands, where the model demonstrates better predictive accuracy for known CPT data. ...
Efficient geothermal resource development remains challenging due to inherent geological uncertainty and limited subsurface data. A proof-of-concept for a digital twin for a fluvial geothermal reservoir, similar to the Delft campus geothermal project, is presented. This digital twin has the aim to integrate geological scenario modeling, production simulation, uncertainty analysis, and data assimilation to mitigate operational risks, reduce maintenance costs, extend reservoir longevity, and enhance the overall sustainability of this project. In this contribution, we assess the efficiency of the ensemble smoother with multiple data assimilation (ESMDA) for subsurface property inversion of a fluvial geothermal system. First, we developed an efficient method that allows for the swift creation of multiple geological scenarios of channelized reservoir geometries, fully constrained to well information, using Rapid Reservoir Modeling (RRM). Next, we generated an ensemble containing multiple geological realizations for a given scenario representing the geothermal system using stochastic reservoir modelling. For a single scenario and its ensemble of stochastically generated property distributions, heat flow and production rates were simulated using the Delft Advanced Research Terra Simulator (DARTS). One of the ensemble members and its simulated production data were taken as the “truth” (or reference) case. ESMDA was then employed to invert the property distribution within the fluvial channels of all other ensemble members, using the “observed” temperature and pressure data along the injection and production well from the “truth” case. We also consider the presence of a monitoring borehole to analyze how additional monitoring data impacts the convergence of ESMDA. The simulation results of the posterior models demonstrated a significant reduction in root mean square error for temperature and pressure data which align more closely with the “observations” compared to the prior models. This outcome confirms the feasibility of applying ESMDA for data assimilation in fluvial geothermal systems, such as the Delft campus geothermal project. ...
Journal article (2025) - Mouadh Rafai, Diana Salciarini, Philip J. Vardon
The effect of the load level on long-term thermally induced pile displacements and the impact of cyclic thermal loads on the bearing capacity of energy piles are investigated via a full-scale in situ test in Delft, The Netherlands. The pile was loaded to a specific target of 0, 30, 40, or 60% of its calculated ultimate bearing capacity. At the end of each loading step, up to ten cooling–natural heating cycles were applied. The pile behavior during monotonic cooling and cyclic cooling–natural heating in terms of the displacement along the pile is reported, with a focus on permanent displacements. During monotonic (pile/ground) cooling, a settlement of the pile head and an uplift of the pile segment near the pile tip were observed in all four tests. In addition, under higher mechanical load, the pile head displacement was larger while the uplift was lower due to the imposed mechanical load. During cyclic thermal load, under zero mechanical load, pile head displacement was fully reversible while permanent uplift of the lowest pile segment was observed and attributed mainly to the permanent dragdown of the surrounding soil. Under moderate mechanical loads (30 and 40%), thermal cycles induced an irreversible pile head settlement, which stabilized with an increasing number of cycles. In addition, a permanent pile settlement along the pile was observed at the end of these tests. Under high mechanical load (60%), the irreversible settlement along the pile continued to increase with only a slight reduction in rate, being higher compared to moderate mechanical loads. In this test, a normalized pile head settlement of 0.124% was observed after ten thermal cycles. The permanent settlement of the pile under thermo-mechanical loads was mainly attributed to the contraction of sand beneath the pile tip and thermal creep at the soil–structure interface. The pile bearing capacity was observed to increase after thermo-mechanical tests, mainly due to the residual/plastic pile head displacement, which in turn densified sand leading to an increase in tip resistance. ...

Performance analysis and optimisation

Journal article (2025) - Marco Gerola, Francesco Cecinato, Vincent Leclercq, Philip J. Vardon
Energy Quay Walls (EQWs) are innovative energy geostructures with the unique capability to exchange heat with both soil and open water. Although previous laboratory testing demonstrated a promising energy efficiency for this type of system, its novelty necessitated thorough research to advance comprehension of its thermal behaviour and optimise energy efficiency. This paper conducts an in-depth examination of EQWs, employing numerical models validated against real data from a full scale test in Delft, The Netherlands.

Two Finite Element numerical models were developed to (i) reconstruct the undisturbed (i.e. pre-geothermal activation) temperature profile within the soil and (ii) conduct a comprehensive (3D) analysis of heat exchange processes in an EQW application (i.e. during geothermal activation), calibrating relevant parameters with field test data, providing valuable insights into its energy efficiency. Following validation, the geothermal activation model was employed to assess the impact of the flow regime within the heat exchanger pipes and the velocity of the open water on the energy efficiency of the EQW system. Additionally, the contributions of soil, water, and air to the energy gain are investigated. The results indicate that the primary source of energy gain is from open water, and the dominance of this contribution is further increased by the presence of turbulent flow within the heat exchanger pipes. However, the soil can play a key role in short term energy delivery. Furthermore, this study emphasises the importance of the open water movement, revealing a 48% reduction in energy extraction for fully stationary water scenarios. ...
Journal article (2025) - Philip J. Vardon, Anne Catherine Dieudonné, John S. McCartney, Jean Michel Pereira, David Smeulders, Guillermo Narsilio
It is clear that to address climate change, an energy transition which makes a large-scale use of the subsurface is needed. The subsurface will play a critical role in this transition, serving as a resource for new sources of energy production and storage, a foundation for energy infrastructure, and a repository for waste by-products from energy production (e.g., radioactive waste disposal, CO2 geo-sequestration). Furthermore, there are challenges in understanding material behaviour due to complex coupled phenomena, measuring material properties and upscaling the physical phenomena to engineering scale structures. Uncertainties, heterogeneities and long timescales offer additional challenges, as does bringing technology ever closer to dense populations. This is the topic of Energy Geotechnics. In the next decade and decades, society needs to complete the energy transition, and to do so the already substantial changes need to be vastly accelerated. This brings many challenges, which academics, consultants, contractors and authorities need to address together. [...] ...
Journal article (2025) - Mouadh Rafai, Diana Salciarini, Philip J. Vardon
This paper presents quantitative data from a field test on a new type of energy pile, called a displacement cast in situ energy pile. The test pile was installed in a multilayered soft soils and subjected to a continuous cooling for 3 months, with no mechanical load. Afterwards, the pile was loaded to a specific target of 20 or 60 % of its calculated ultimate bearing capacity and then subjected to up to five thermal cycles. Under zero mechanical load, the results revealed that the compressive/tensile stresses coexist along the pile. Under low mechanical load (20 %), thermal cycles induced irreversible residual contractive strains and stresses as well as a limited pile head settlement. Under high mechanical load (60 %) and extreme operating conditions, i.e., negative temperatures which could have indicated a frozen interface, further irreversible settlements observed at the end of this test. Mechanical pile tests however indicated no impact of stress history (including the freezing test) on the shaft resistance and the overall pile-bearing capacity. ...
Abstract (2025) - Amirhossein Hashemi, Martin Bloemendal, Philip Vardon, Pieter Goverse, Guus De Rechter
Underground Thermal Energy Storage (UTES) technologies are essential for advancing low-carbon heating and cooling systems, particularly in urban areas where space constraints and retrofitting challenges pose significant barriers. In this study the performance of a system of novel coaxial diagonal borehole heat exchangers (BHE) is analyzed during September–December 2024.

The Home Smart Energy (HSE) system, implemented in Medemblik, Netherlands, features a nine-borehole diagonal array arranged in a circular configuration. The boreholes are drilled at a 60° or 45° angle to depths of up to 40 meters, operating in a closed-loop coaxial setup. A brine mixture of water, operates with a flow rate of 3100 l/h, and 14% glycol lowers the freezing point below 0°C, allowing the system to supply higher capacities. The heat pump extracts the heat from the BHE’s, supported by solar thermal collectors to charge the BHE’s in summer, ensuring efficient year-round heating. An extensive monitoring framework, including Distributed Temperature Sensing (DTS), provides detailed insights into system performance during operation.

The HSE system demonstrated consistent performance under varying configurations and conditions. With all nine boreholes active, the system achieved a seasonal Coefficient of Performance (COP) ranging from 3.8 to 5.2, with daily energy outputs averaging 125 to 220 kWh/day. During December 2024, tests were conducted using three boreholes in different configurations at a reduced flow rate of 2800 l/h. These tests showed that borehole arrangement moderately influenced system performance, with the adjacent configuration achieving slightly higher energy outputs and COP, compared to the dispersed configuration.

The system also demonstrated significant energy cost savings of €954 during November and December 2024, attributed to a reduction in gas consumption by over 700 m³ compared to the previous year. These findings confirm that diagonal shallow co-axial borehole arrays are a scalable and sustainable UTES solution, offering substantial energy savings and CO₂ reductions in dense urban settings. ...