A New Approach for Characterizing Pile Heat Exchangers Using Thermal Response Tests

: Pile heat exchangers offer a cost effective route to implementation of ground-source heat pump systems for many large commercial buildings compared with traditional boreholes. Such pro-jects typically use thermal response tests to determine the key input parameters for system design, namely soil thermal conductivity and heat exchanger thermal resistance. However, this brings chal-lenges for pile heat exchanger based systems, where in situ thermal response tests are known to be less reliable due to the large thermal capacity of the pile. This paper presents a new “black box” resistance capacitive model for applications to pile thermal response tests. The approach is tested against case study data and shown to perform well. Additional test duration savings are shown to be possible if a novel combination of borehole and pile thermal response tests is applied together to determine design parameters.


Introduction
Ground-Source Heat Pump (GSHP) systems can decrease the emission of greenhouse gases resulting from heating, cooling and hot water provision. Energy geostructures, which combine structural and thermal function, are an opportunity to reduce the installation costs of ground heat exchangers (GHE), which are the in-ground component of the GSHP system. Energy geostructures have been installed for several decades [1]. This includes the equipping of foundation piles [2,3], embedded retaining walls [4], or tunnels [5] with plastic heat transfer pipes so they can act as a GHE. Of the types of energy geostructure, pile heat exchangers (PHE), e.g., Figure 1, are the most common. This is partly because most are likely to be associated with an overlying building and hence have ready users of the heat they can supply, but also because their typically cylindrical shape makes them superficially similar in geometry to the more common borehole heat exchanger (BHE). Thermal analysis methods for BHE are well developed and hence are typically drawn on for application to PHE.
However, the radius of a PHE, which can exceed 50 cm, is much greater than the radius of a borehole heat exchanger (BHE) (typically 8 to 10 cm). PHE are also usually much shorter than BHE (typically 10 to 20 m against 100 m to 200 m). Therefore many approaches to the analysis of BHE, which assume the presence of a very long and thin heat source, may not be applicable to PHE. Consequently, more appropriate methods for thermal modeling such as fully discretized models (e.g., finite element analysis) to investigate the pipe arrangements and thermal performance [6,7] or semi-analytical models for thermal dynamic simulation of the whole energy system with hourly time steps [8] are often applied.

Figure 1.
Example of a pile heat exchanger following pile trimming and prior to connection to the building base slab and the ground-source heat pump system. In this particular case the heat transfer pipes are centrally located and still retain their protective sleeves, which were used to avoid their damage during trimming.
During a TRT, a heat-carrying liquid (usually water) circulates in the pipes of the ground heat exchanger with a constant flow volume Qv [m 3 ·s −1 ] while an electric heater supplies a constant power P [W] to the fluid. The entry Tin and exit Tout temperatures of the pile are recorded for the duration of the TRT [9]. The "classical" interpretation of a TRT is based on a purely resistive thermal model developed for BHE [9], i.e., a model overlooking the thermal inertia of the GHE. It assumes that the evolution of the temperature of the heat-carrying liquid Tf is described by the sum of the resistive component and transient infinite line source (ILS) step response, which takes the borehole to be a line emitting power with a constant rate: In Equation (1), T0 is the undisturbed ground temperature, p is the linear power (W·m −1 ) defined as the ratio between the power P and the depth of the exchanger H, is the borehole thermal resistance, is the effective thermal conductivity of the ground, ( ) = ∫ d , t* is the normalized time or Fourier number: The first item is obtained based on an inspection of the initial fluid temperature before heating has commenced, while for 2 and 3 above, an approximation of Equation (1) is generally used, with the Euler constant ( ≈ 0.5773): Equation (3) is valid as soon as a stationary thermal regime in the borehole is reached, i.e., t* > tmin*, which leads to the exclusion of the temperatures measured before tmin of the interpretation. λm and Rb can be identified graphically from the slope and intercept of the curve = (ln( * )). An alternative is to minimize the root mean square error (RMSE) ({ }) (°C) between the measured temperature Tfl and that calculated by the model Tfl,mod: where { } contains the parameters to back-calculate.
Initially developed for BHE, the model described in Equations (1) and (3) do not account for the thermal inertia of the backfilling material. Past work shows that this purely resistive model (Figure 2a) is not suitable for pile heat exchangers of large diameters [11,12], since once λm and Rb have been fitted, the temperature change at small time scales (e.g., 1-10 h) is overestimated by several °C, leading to an underestimation of the transfer capacity of the PHE. Moreover, the duration required to thermally load the backfill material tmin* is barely compatible with the operational constraints of a construction site. Indeed, considering the typical criteria tmin* = 5 [13] and assuming a soil with thermal properties = 1.  Despite these limitations, many authors report application of TRT to PHE to determine the piles and ground properties [11,12]. In these studies, the measured fluid temperature is typically interpreted with the ILS approach, which leads to the exclusion of a significant part of the temperature measurements before tmin* reaches 5. Alternative approaches have been reported in [11,[14][15][16]. Loveridge et al. [16] used pre-defined pile step responses, which lead to values of and Rb similar to that determined by the ILS, while improving the forecast of the temperature changes at the early stage of the TRT. However, the pile step responses had to be chosen from a small number of representative numerically pre-calculated cases for different PHE geometries. When used for TRT back-calculation, selection from these cases is made before the analysis ( [15,17]). It is also not possible to look for the characterization of the heat transfer inside the PHE, since the nature of this is an assumption of the case chosen. This means the approach is not appropriate in all cases and may contain errors related to choice of the closest geometry. Zarrella et al. [14] used the computational Capacity Resistance Model (CaRM) [18] to account for the grouting material inertia. CaRM takes a resistive-capacitive approach and has been developed for specific geometries including for single-U, double-U and coaxial heat exchangers. It has the key limitations that the pipe arrangements in the borehole are predefined (e.g., "pipes close to the borehole wall") and resistances describing the inner thermal transfer must be pre-calculated with a distinct method, such as finite element code. Once the inner parameters have been tuned, CaRM was in very good agreement with a TRT performed on a BHE of 140 mm diameter. Alberdi et al. compared several analytical approaches, as well as a finite element model representing an horizontal PHE cross section [16]. They showed that models that considered both the short length of the PHE and took account of its internal capacity gave the best results when compared with a benchmark 3D numerical simulation.
However, the main approaches, which account for the pile thermal capacity, have limitations as described above, making them applicable only in certain circumstances, or applied with greater error when generalized. Consequently, there remains a need for a flexible model, which can be specific to a given pile geometry, accounting for concrete thermal inertia and is easy to use and can be embedded in a numerical procedure of backcalculation for TRT interpretation.
In this paper we present a simple model of a PHE with thermal inertia, which deals with PHE of large diameters (rb ≫ 10 cm) and is specific to a given pile geometry (Section 2). The model is tested by back analysis of two experimental TRT data sets, where the PHE parameters are back-calculated and compared with the classical ILS interpretation with Equations (1) and (3) (Section 3). We show that the temperature rise at small time scales (i.e., t ≈ 1 h) can be better accounted for with the new simple approach. Additionally, the expected errors are comparable or less compared with standard borehole TRT interpretation. However, the time to convergence of the new model output means that the soil thermal conductivity cannot be determined from a short TRT on a pile. However, we go on to demonstrate that a novel combination of BHE and PHE TRTs interpreted using the new model can lead to a reduction in pile testing time. On this basis we estimate the minimum TRT duration for PHE (Section 4) potentially saving time on practical operations.

Methods
Section 2.1 describes the field test data used in the study and Section 2.2 describes the new resistive capacity model for PHE TRT interpretation. Analysis of the data using the new model and the classical approach introduced in Section 1 is then set out in Section 3.

Experimental Data
Two TRT, summarized in Table 1, and respectively referred to as set B and set C, are analyzed in this paper. Set B was undertaken on a pile of 22.5 cm radius at Richmond, Texas, USA, mostly located in saturated sand [15,19]. Set C was a TRT undertaken on a pile of 30.0 cm radius located in London clay [12]. The inner details of the piles are given in Table 1 and illustrated in Figure 3. Note that for set B, the power was shut down at t = 103.4 h, and restarted 4 h later until t = 140.2 h. The classical interpretation was carried out for t = 103.4 h. For both TRT, the heat-carrier fluid was water, whose properties are assumed to have a density ρfl = 1000 kg·m −3 and a heat capacity Cp,fl = 4180 kJ·kg −1 .K −1 . The temperatures, flow-rate and power were monitored with a time step of 5 min for set B and 1 min for set C. Data of set C were averaged with a time step of 5 min for analysis purposes.
The thermal inertia of the concrete means that the evolution of the fluid temperature goes away from purely linear behavior depending on the logarithm of time foreseen by the ILS model, i.e., Equation (1) (cf. Figure 4). This is even more noticeable for set C, the pile of greater diameter. The temperature measured in set B also exhibits slight periodic variations of 24 h [c.a. 0.5 °C], suggesting the fluid temperature is influenced by the atmosphere. One explanation may be that the pipes connecting the PHE to the TRT module may not have been sufficiently insulated. Note that in what follows, (ρCp)m = 2.4 MJ·K −1 ·m −3 is assumed for both soils based on typical values [10].

Resistive-Capacitive Model for a Pile Heat Exchanger
A model has been developed to take into account the thermal inertia of the backfilling material. The heat transfer inside the concrete relies on a linear capacitance [J·K −1 ·m −1 ] located between two resistances R2 and R3 [K·m·W −1 ] (Figure 2b). The transfer outside the pile considers a hollow cylinder, and relies on the "classical" infinite cylindrical source (ICS) model [20].
The resistive-capacitive model (RC) in Figure 2b is qualified as semi-analytical since the transfer inside the pile is treated numerically whilst the transfer in the ground is treated analytically. The capacitance is simply estimated by: = (5) (ρCp)c is the volumetric calorific capacity of the concrete [J·K −1 ·m −3 ], estimated as 2.11 MJ·K −1 ·m 3 , a value compatible with preceding studies [21,22]. A parameter x represents the position of the capacitive node (Figure 2b) within the borehole total resistance Rb: The definition of the fluid mean temperature Tf remains unchanged: A power balance on the fluid leads to: In what follows, the equations are discretized with a time step Δt. Δt = 15 min has been chosen to remain in line with the experimental measurements.
An analytical solution to the heat equation, the infinite cylindrical source (ICS) ( * ) [20], allows the calculation of the temperature at the borehole wall (i.e., at the ground/pile interface): Under steady-state conditions (constant heat flux), the borehole wall temperature is related to the step response G(t*): During a TRT, the thermal power provided by the test module is constant. However, this power first heats up the concrete before being progressively transferred to the ground. Therefore the steady-state condition cannot be assumed straightaway. The superposition principle is introduced to estimate Tb n at time step n [23]: In Equation (13), the superscripts refer to the time steps, e.g., ≈ ( ∆ ). pb is the linear power received at the borehole wall: The Equations (13) and (14) combine as follows: With: The Equations (7)- (10) and (15) recombine in a matrix system as follows: At each time step n, the vector of the temperatures at the following time step {T} n+1 is determined by an implicit Euler scheme, i.e., with the following approximation: The typical execution time of the numerical model is a fraction of a second.
Note that the model considers that all temperatures are constant along the pile and neglects the axial heat transfer. However, since the yearly variations of temperatures affect the ground up to a few meters, additional research is still required on how these changes, which result in non-uniform initial ground temperature T0 measurements, could affect the TRT result.

Analysis Approach
The initial temperature was determined from the early TRT phases, before heat was injected to the ground. This circulation phase lasted 9 min for set B and 60 min for set C. Averaging the measured temperature yields T0 = 24.97 °C for set B and 14.23 °C for set C. For set B the value was confirmed to be appropriate with reference to three temperature sensors in a nearby borehole. Less than 0.5 °C variation in ground temperature was observed over 18 m depth.
In this paper, the RMSE (Equation (4)) was minimized with the Matlab ® software. { } was determined with the fmincon function, using the SQP (sequential quadratic programming) solver. The active-set solver was tested, but appeared to stop before converging to a local minimum of ({ }) in some cases, providing unreliable results. A second stage has the benefit of using the fitnlm and coefCI functions to fit nonlinear regression models and to determine the 95% confidence intervals (CI) on the fitted parameters. The fitnlm function was run with the results of fmincon as an initial point.
The first step was to interpret the TRT with the classical method described in the introduction. Setting the criterion tmin* = 5, the minimum duration tmin was 60 h for set B and 214 h for set C whilst estimating, respectively, λm at 2.8 W·K·m −1 and 1.4 W·K·m −1 . Indeed, as tmin depends on λm which is sought by the interpretation (see Equation (2)), tmin was determined by a manual trial and error process, increasing or decreasing tmin, then minimizing ({ , }), until the condition on tmin* was reached. The investigation was also undertaken with tmin* = 3. Regarding set C, previous studies [12] reported λm between 1.35 and 1.45 W·K·m −1 depending on the values chosen for tmin* (from 5 to 7) and tmax (from 250 h to 350 h) so this new analysis is consistent with previous work.
Having established tmin, the root mean square error between the field data and the modeled data was minimized over all time steps from tmin until tmax, where tmax was the last data time step used in the analysis and was increased in stages until the end of the test period. In this way, the effect of the data window could be investigated and for an ideal data set convergence of derived parameters (Rb, λm) would be seen with time.
Secondly, to capture the transient phase at the beginning of the TRT, the resistivecapacitive model presented in Section 2.2 was used to back-calculate the ground thermal conductivity λm and the borehole resistances Rb and x parameter by minimizing the misfit Equation (4), i.e., {X} = {λm, Rb, x} with the fitnlm function. Note that the first optimization with fmincon was run with {X} = {λm, R2, R3}. Without the limitations of the stationary assumption for the PHE in the classical method, tmin was fixed to 1 h for both datasets, to capture the transient heat transfer in the pile.

Classical Interpretation
As tmax is increased the root mean square error, thermal conductivity and thermal resistance should converge as the model error is reduced as the PHE thermal capacity is overcome with time. Set B (Figure 5a,c,e) presents a different behavior since the values of λm and Rb do not converge when the duration of the test increases. There are two factors, which may contribute to explaining this result. First, the cyclic daily perturbation of the measured temperature between the pile and the near surface interfere with the interpretation, resulting in a 24 h cyclic pattern in both λm and Rb. It has been observed elsewhere that for TRTs performed on piles, effects like this can lead to increasing errors since any perturbations in temperature due to the environment have a proportionally greater impact later in the tests when the rate of absolute temperature change is lower [24]. Secondly, it has been reported [24] after t* = 10 the pile thermal response will be affected by end effects leading to systematic overestimation of thermal conductivity. For set B, t* was 9.9 at the end of test, so this effect could also impact the late time accuracy.
However, the variation of λm (between 2.7 and 3.2 W·K −1 ·m −1 , that is 15%) remains acceptable in as much as the uncertainty in the interpretation of a TRT is of the order of 10% [23]. The misfit soars for tmax = 105 h since it encompasses temperatures measured after the power shutdown. The daily perturbation of temperature results in narrower intervals of confidence.
For set C, λm and Rb converge if the duration of the TRT tmax is above approximately 220 h for tmin* = 3 (cf. b, d, f). λm and Rb stabilize, respectively, at 1.44 W·K −1 ·m −1 and at 0.136 K·m·W −1 . The interpretation is similar with tmin* = 5 as soon as tmax = 260 h: using tmin* = 3 allows the reduction of the duration of the TRT by only 10%. Beyond this threshold at ≈ 260 h, λm and Rb vary by less than 3% if tmax increases: the thermal transfer at the borehole wall is stationary, which validates the resistive model. tmin was estimated to ≈ 220 h for λm ≈ 1.4 W·K −1 ·m −1 , which means that about 84% (≈(260 − 220)/260) of the TRT duration only served to thermally charge the concrete-and only 40 h was actually used for the interpretation.
The dependence of λm with tmax has been reported in [24] for a TRT exhibiting daily perturbations performed on a 25 m deep BHE, as for set B. Set C does not exhibit such a behavior. One explanation may be that set C PHE is 70% deeper than set B PHE, while λm is half, resulting in a shallower affected thermal zone. A more realistic explanation is that set C measurements do not need to be corrected since specific care has been taken to insulate the pipes. For set B, the perturbation of the fluid temperature, estimated as a deviation from the fitted Jacob approximation of the ILS (Equation (3)), can reach c.a. 0.4 °C, while it remains below 0.1 °C for set A.
The engineer sizing a ground-source heat pump is concerned with the ability of the models representing every system component to forecast the evolution of the fluid temperature at small, medium and large time scales. The capability of the purely resistive model coupled to the ILS (Equation (1)) or the ICS (Equation (11)) has been assessed for both experimental datasets ( Figure 6). The TRNSYS Type 557, which is based on the duct storage model (DST) for ground heat exchangers [25], has been included in the comparison as well. The DST uses superposition of three elements: a heat balance on the fluid, a local resistive process close to an individual pipe (so-called "duct") rather than the borehole, and a global process. All the models use the same values of ground conductivity λm and borehole resistance Rb, obtained by the "classical" interpretation with tmin* = 5 for both sets, tmax = 103 h for set B and tmax = 350 h for set C. Equation (3) is also represented as it is used for the interpretation, though it may lead to a significant underestimation of the temperature at small time scales due to its logarithmic behavior. After 1 h of thermal load the ILS model (Equation (1)) overestimates the temperature of the fluid Tfl by ≈5.0 °C for set B and by ≈3.3 °C for set C. For both datasets, neglecting the thermal inertia of the concrete leads to an underestimation of the pile heat transfer capability, and therefore possibly a cost-ineffective geothermal equipment of many more piles than required. Indeed, the heat pump could exchange a larger amount of thermal energy with the underground than a purely resistive model would forecast.
When parameterizing the duct storage model, the resistance Rb was taken to be equivalent to the "fluid to ground resistance" and for simplicity no account was taken of any additional pipe-to-pipe interactions. The DST is remarkably able to reproduce the temperature evolution. This is, to some extent, surprising since DST is a purely resistive model in terms of internal heat transfer and it does not allow defining the pipe locations in the borehole or pile [25]. However, DST has been primarily designed for borehole thermal energy storage. Consequently, the ground heat exchangers are assumed to be located within a cylinder at the nodes of a hexagonal grid, which will not be realistic for most PHE fields. Meanwhile, thanks to the superimposition principle [21], analytic solutions are far more flexible to handle arbitrary borehole locations. A framework can even be built to select which piles to equip in a field [8]. Figure 7 shows the back-analyzed thermal conductivity for the two datasets using the RC model coupled to the infinite cylindrical source. For set B, the conductivity oscillates in the range 3.1-3.3 W·K −1 ·m −1 as long as t < 103.4 h, so before the power shutdown (Figure 7c). λm exhibits smaller variations than the classical interpretation, suggesting that the new approach is better capable of dealing with daily perturbations in the measured temperature since it uses temporal superposition. However, though the RC model has the ability to deal with time-varying temperature changes, it fails to properly estimate the ground conductivity when encompassing data after the power shutdown, λm then soaring to the range 3.7-3.9 W·K −1 ·m −1 .

Interpretation with the New Resistive-Capacitive (RC) Model
For set C, as soon as tmax > 250 h, the conductivity obtained by the inversion of the RC model converges towards λm = 1.43 W·K −1 ·m −1 , a value almost identical to that obtained by the normal interpretation (Figure7d). The convergence of Rb and x beyond 250 h can also be observed (Figure 7f). Hence, while the RC model appears to allow successful interpretation of set B, which was not possible with the classical approach, the RC interpretation does not allow the reduction of the duration of the TRT compared to the resistive model in this case.

Error Analysis
An in-depth analysis of the error on the back-calculated parameters has been performed by Witte for a conventional TRT on a BHE, given that analytical expressions of λm and Rb are known [23]. For comparison, the estimated error in the present study on the three back-calculated parameters is determined as follows: where account for the input parameters and Δ for the related error on this parameter. The partial derivatives are computed numerically by changing the input value by a small amount (±1%), then minimizing the RMSE to get the updated value of , and computing the relative change of . The measurement error has been computed for a heat injection duration of 100 h for set B and 350 h for set C. The values of Δ are given in Table 2, with the following comments:  According to SIA standards, the volume-specific heat capacities of wet clay and wet sand are, respectively, in the range 2.0-2.8 MJ·K −1 ·m −3 and 2.2-2.8 MJ·K −1 ·m −3 [10]. We considered the error on the ground capacity as the half of these intervals, i.e., 0.3 MJ·K −1 ·m −3 for set B and 0.4 MJ·K −1 ·m −3 for set C. Typical ranges for concrete heat capacity could not be found in the literature, but given the values reported in previous studies [26,27], an error of 0.2 MJ·K −1 ·m −3 was chosen.  The error for the pile diameter and height were determined according to the UK specification for construction tolerances [28]. In this respect, it should be noted that the dimensions of a constructed pile should not be less than the specified dimensions. A tolerance on these dimensions of up to the lesser of 50 mm or 5% is permissible.  For set B, the test was performed with reference to the ASHRAE standard [29]. This states that the accuracy of temperature measurement must be less than 0.3 °C, for power measurements less than 2% and for flow rate measurements less than 5%. These are conservative values, since the test may have been performed with more accurate instruments.  For set C, the client specification had tighter accuracy requirements, which can be reasonably applied. They would be error for temperature of 0.1 °C, flow measurement of 0.01 m·s −1 and power to 5W. For both datasets, the most significant sources of error are related to the heat capacities of the ground and the concrete (refer to Figure 8 which represent the distribution of Δ terms). As expected, temperatures and power play a more significant role for set B, where the standards are less restrictive. Errors on the pile diameter affect the inner parameters Rb and x significantly. The error on the ground thermal conductivity is 0.06 W·K −1 ·m −1 for set C (about 4%) and 0.20 W·K −1 ·m −1 for set B (about 6%) reflecting the more restrictive standards for set C. These values are well within the expected bounds for thermal response testing for boreholes, giving confidence to the new approach. For comparison, for a classical TRT on a BHE, Witte reported a typical error of 5% on the ground thermal conductivity [23].

Pile Thermal Capacity
The interest of the RC model on the resistive model is in the transitory phase. Here we arbitrarily define this phase as t* < 2.5, which results in 27 h for set B and 107 h for set C (Equation (2)). Once their inner parameters have been fitted, the classical and RC methods are compared in Figure 6. Given the similar results between the two methods in terms of conductivity prediction, for comparative purposes, the same input parameters are used in Figure 6. Note that both methods also use the normalized duration tmax* of experimental data, so tmax = 103 h for set B (tmax* ≈ 9.9) and tmax = 350 h for set C (tmax* ≈ 8.5). After 1 h of thermal load, the RC model overestimates the temperature of the fluid Tfl by only ≈ 1.0 °C for set B and ≈1.3 °C for set C. Beyond the transitory phase (t* > 2.5), the two approaches accord well with the measured temperature and both correctly forecast the measured temperature, especially when the t* > 5 criteria is reached. Further, the RC model allows an understanding of the dynamics of heat transfer within the pile (Figure 9). The pile is halfloaded at t* ≈ 0.4 for set B and t* ≈ 0.2 for set C.

Discussion and Recommendations
Whatever the interpretation method, the duration tmax for set C must be of the order of 250 h, i.e., about 10 days. This duration is not necessarily compatible with the management of a construction site. Therefore it is clear that the ground thermal conductivity may be obtained from a much shorter TRT carried out on a special purpose BHE as recommended by [17]. Yet this proposed approach retains a drawback, namely that the pile itself has not been characterized to allow appropriate forward simulation for the design process. No information has been obtained about its thermal resistance. While a closed form analytical expression has now been developed to obtain pile thermal resistance [30] this still requires information about concrete thermal conductivity and assumes accuracy of pipe placements compared with design. Therefore we propose that by using the RC model the combination of a short BHE TRT and an additional short TRT on one of the constructed piles may provide the best solution. One TRT is performed on a BHE having a small radius (e.g., rb = 8 cm) to determine λm within a few days, while simultaneously a TRT is performed on the PHE to obtain the inner parameters for the RC model. This undertaking assumes that the ground is relatively homogeneous and lateral variations of composition and ground properties can be overlooked. This approach is illustrated with respect to datasets B and C. The best-fit values of λm are assumed to have been obtained independently from a BHE TRT. In this case values of 3.24 W/mK and 1.43 W/mK are used, respectively, for sets B and C, based on the results in Section 3. The RC model is then optimized for the pile TRT data with only two rather than three unknowns. The fitted values of x and Rb start to stabilize around tmax ≈ 30-40 h for set B at x ≈ 0.56-0.58 and Rb = 0.080 K·m·W −1 and tmax ≈ 100 h for set C at x ≈ 0.76-0.78 and Rb = 0.121-0.124 K·m·W −1 (Figure 10). For both datasets, the difference, with the value of Rb determined in Section 3, does not exceed 4%.
It is desirable to determine the minimum test duration tmax beyond which the determined values of Rb and x are reliable. Here we use as a criteria the ability of the RC model to forecast the evolution of the fluid temperature at the very end of the data record, from tmin' to from tmax'. A second RMSE indicator ε' is build considering tmin' = 80 h and tmax' = 100 h for set B, and tmin' = 300 h and tmax' = 350 h for set C. For set C, ε' reaches a minimum around tmax* = 100 h (tmax* = 2.4), which suggests this duration is sufficient to estimate the temperature evolution from 300 to 350 h, as confirmed by the comparison to the recorded temperature ( Figure 11).

Conclusions
A resistive-capacitive PHE model has been developed to capture the transient energy transfers within a pile. This model has been shown to perform well and in particular:


Numerical back-calculation of the model parameters on two thermal response tests yield similar values of ground conductivity and thermal resistance as the well-established infinite line source model.


Inclusion of temporal superposition with the model allows reliable results to be obtained even when tests are affected by ambient air interference.  The RC model better represents the transient phase of pile warm-up in the early part of the test (approximately up to a Fourier number t* = 1 to 2).  The errors associated with the calculation of thermal conductivity are all less than 10% and well within expected ranges for boreholes thermal response tests interpreted with the classic infinite line source.  Standard back-calculation using the RC model does not allow to significantly reduce the TRT duration below t* = 5.  However, if the thermal conductivity can be obtained by another means, the time for the RC model to converge is much reduced meaning that pile resistance can be obtained from a pile TRT in a duration corresponding to a Fourier number t* ≈ 2 to 2.5.
Given these characteristics of the RC model the following novel approach for pile characterization is recommended:


Use a borehole at the same site and of the same length as the piles to carry out a BHE TRT to determine the effective soil thermal conductivity using the classical approach.  Carry out a short duration pile TRT according to Fourier number t* ≈ 2 to 2.5, or around 100 h for the cases demonstrated in this paper.  Interpret the pile TRT using the RC model to determine both the pile thermal resistance and the inner resistances of the RC model, which can then be used in forward simulation for design purposes.