Circular Image

P.J. Vardon

info

Please Note

182 records found

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. ...
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) - 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) - 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) - 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. ...
Clay-rich geological formations are considered as host rocks for deep geological disposal of radioactive waste. Over the long term, gas will be produced and will migrate through the surrounding geological formation. Gas transport mechanisms have been investigated in laboratory tests. However, the effects of material heterogeneity remain insufficiently explored. This paper presents a stochastic analysis of two-phase flow in clays under gas injection, incorporating spatially correlated porosity. The study evaluates the effects of sample size, gas injection pressure, and the choice between two-dimensional (2D) and two-dimensional (3D) conditions on the statistical outputs, including mean behaviour and variability. The results indicate that larger samples exhibit reduced variability in the degree of saturation and gas permeability due to enhanced averaging effects. Moreover, the variation in results is higher under high gas injection pressure compared to low gas injection pressure. In addition, the variability of results is significantly reduced in 3D simulations compared to 2D, with high-permeability regions more likely to form continuous pathways under 3D conditions, emphasising the necessity of accounting for 3D effects. The findings indicate that sample size is a critical factor in experiments, as it influences the number of tests required to achieve results within a desired level of accuracy. ...
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. ...
Geothermal energy is a key option for decarbonizing heating and cooling in the energy transition. Forecasting geothermal production has inherent uncertainty due to the heterogeneity of geological formations that host the geothermal resource and the limited data available to characterize and quantify these heterogeneities. This uncertainty leads to operational risks such as early thermal breakthrough. Identifying the most valuable monitoring data and data acquisition strategies for operators is key to constraining uncertainties and ultimately de-risking operations in a reliable and cost-effective way. Data-worth analysis quantifies the value of data provided by existing observations or proposed data collection strategies. This study combines Ensemble Smoother with Multiple Data Assimilation and data-worth analysis to constrain uncertainty in production forecasts and reservoir response for a geothermal doublet system located in a clastic, channelized fluvial reservoir. The main monitoring data includes production temperature, injection pressure, and temperature and pressure profiles along the well paths. We show that production temperature and injection pressure alone only can constrain uncertainties in production forecasts. Using observations of well temperature and pressure profiles demonstrates a threefold increase in data worth, which improves both the quantification of production forecasts and reservoir dynamics. At constant injection rates in doublet system, temperature profiles monitored along the wells provide higher data worth compared to pressure profiles. Data from a deviated monitoring borehole is advantageous. Early-time (first year) observations of temperature and pressure profiles along the injector, producer, and monitoring borehole already constrain production uncertainty prior to thermal breakthrough. Data-worth analysis is shown to be most beneficial when conducted across multiple plausible geological scenarios to ensure a more reliable assessment of collection strategies. The findings of this study yield insight into designing informative data-acquisition strategies for direct-use geothermal systems. ...

Mesh dependence and regularisation in the material point method

Conference paper (2026) - José León González Acosta, Miguel A. Mánica, Philip J. Vardon, Michael A. Hicks, Antonio Gens
Much of our effort in numerical analysis within geotechnical engineering is devoted to evaluating the likelihood of failure in a given boundary value problem (BVP). However, certain problems occasionally require studying the post-failure behaviour of the mobilised soil mass and its resulting consequences. These large deformation analyses exceed the capabilities of conventional finite element formulations, requiring the use of specialised numerical techniques, such as the material point method (MPM), which can mitigate mesh distortion issues. However, since MPM is based on the same principles as the finite element method (FEM), it shares many of its limitations, including volumetric locking and the hourglass effect, as well as additional challenges, such as stress oscillations due to material points crossing element boundaries. Furthermore, when combined with a constitutive description exhibiting softening, MPM can lead to non-objective results with a pathological dependence on the adopted mesh and poor convergence properties. Within this context, the present work addresses the importance of regularisation in MPM for the objective simulation of localised deformations in the presence of brittle materials. A nonlocal approach was incorporated within an existing MPM framework and applied to the simulation of a number of simple BVPs with a softening material. As in conventional FEM simulations, results without regularisation showed a more brittle global response and larger strains and displacements as the element size was reduced. Furthermore, and particularly relevant for studying the consequences of a given collapse, run-out distances were shown to depend on the mesh resolution. On the other hand, regularised simulations exhibited consistent behaviour, with a global response and a configuration of localised deformations that were approximately independent of the employed mesh. However, it was demonstrated that stress oscillation issues must also be addressed when softening is considered to prevent numerical instabilities. ...

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. ...

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. ...
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. ...
Journal article (2026) - Matthijs S. van Esch, Martin Bloemendal, Niels Hartog, Philip J. Vardon
The energy transition relies on electrification and due to the increased use of heat pumps this is especially true for the heating and cooling of buildings, causing congestion on the electricity grid. The Aquifer Thermal Energy Storage (ATES) Triplet can decrease the reliance of low-emission heating and cooling on electricity by utilizing on-site energy generation, two heating and cooling supply wells and one well to prevent thermal pollution of the other wells. With this setup the system provides heating and cooling without the use of a heat pump, only by circulating groundwater for the energy supply and asynchronous energy generation. A simulation study is conducted to identify when the ATES Triplet is a more viable option than a standard ATES doublet. Five building energy variables are systematically assessed with a subsurface thermohydraulic dynamics model, i.e. heating demand, cooling demand, injection temperature levels, cutoff temperature levels, and return temperature levels. These scenarios are evaluated based on the calculated electricity requirements, pumping volume, and on-site energy generation requirements. The results show that the Triplet outperforms an average low-temperature ATES doublet in most of the considered cases, especially cases with low sub-surface losses, a high heating ΔT, a low cutoff temperature and a low return temperature after heating. This requires building HVAC (Heating, Ventilation and Air Conditioning) systems to accept as wide a range of supply temperature as possible. Systems with a low heating demand (< 1 TJ) in combination with a high required temperature for heating (> 80 °C) are unfavourable. ...
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. ...
Journal article (2026) - T. M. Grubben, A. Koulidis, M. Bloemendal, N. Hartog, S. Muraro, P. J. Vardon
Boreholes are commonly used in the shallow and deep subsurface for the extraction or injection of fluids, e.g. for drinking water production, geothermal energy extraction or fluidic waste disposal. Many boreholes have been drilled successfully in unconsolidated sediments, however the common approach to evaluate borehole stability predicts failure when boreholes are drilled according to standard practice for groundwater wells, i.e. with little overpressure (< 40 kPa) in unconsolidated sand. To explore this mismatch, this study investigates the role of shear dilation in explaining the existence of stable boreholes in unconsolidated cohesionless sand through experimental and numerical analyses. Firstly, consolidated drained triaxial tests along different stress paths are performed on proxy and cored sand specimens. The results demonstrate the strong and immediate dilative behaviour of sand subjected to lateral unloading. Subsequently, a calibrated SANISAND material model is used in the numerical analysis of a cavity contraction problem. The results indicate that the amount of deformation and extent of a hoop stress reduction zone depend on the initial state of the sand. The early onset of irreversible deformation observed in both the experimental and numerical findings suggests that the standard practice for determining the minimum borehole overpressure is too conservative for unconsolidated sands. It is concluded that the initial state of the sand and associated shear dilation plays a prominent role in the stability problem of boreholes drilled with overpressure pw ∼ 10 kPa. ...
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. ...
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 (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. ...

Results from a full-scale field test

Journal article (2025) - Philip J. Vardon, Marco Gerola, Vincent Leclercq, Korneel de Jong, Jacco Haasnoot, Richard Janssen, Patrick Stoelhorst, Ivaylo Pantev, Jorrit de Vries, More authors...
Energy Quay Walls (EQWs) are innovative energy geostructures which exchange thermal energy with both soil and open water while providing a structural function. A full-scale EQW with thermally activated sheet piles was tested, measuring 8.4 m in length along a 1.75m deep canal, with the sheet piles embedded 13m into the underlying soil. Two different length heat exchangers were used: shallow (3 m length) loops primarily extracting thermal energy from the open water, and deep (15 m length) loops extracting energy also from the soil. The shallow loops demonstrated a high heat extraction rate per activated surface area (∼200 W/m2 at 8 °C water temperature, compared with ∼60 W/m2 for the deep loops), with their performance closely linked to the open water temperature. The shallow loops did not require time to restore surrounding temperatures, indicating stable long-term performance, yet can extract the least energy at the coldest time periods. In contrast, the deeper loops exhibit greater stability across varying open water temperatures and achieve the highest total energy extraction per quay wall length (∼900 W/m at 8 °C water temperature, compared with ∼600 W/m for the shallow loops). Realistic operation of the deep loops lowered the soil temperature by ∼2 °C. ...
Conference paper (2025) - J. Liaudat, P.J. Vardon, M.A. Hicks, A.C. Dieudonné
Gas-induced fracturing in liquid-saturated clay-rich materials presents challenges in understanding and predicting fracture behaviour, due to the complex mechanical and transport properties of clays and the compressibility of gas. This paper introduces a novel experimental device for visualising fluid-driven cracks in clays. The device allows for the induction and observation of two-dimensional cracks in clay-rich, low-permeability materials through the injection of gas or water. The experimental setup comprises precision instrumentation for measuring compression forces, displacement, and fluid pressure, along with high-resolution imaging capabilities. Preliminary tests with Helium gas injection into Boom clay samples demonstrate the device's ability to track fracture evolution. This innovative experimental tool offers insights into the mechanisms governing fluid-driven fractures in clay-rich materials and provides a means to validate numerical models. ...