Skip to main content

Science – Society – Technology

Thermal resistance of borehole heat exchangers composed of multiple loops and custom shapes

Abstract

Background

The thermal resistance of a borehole can be reduced by employing thermally enhanced grout, increasing the surface area of the loop and locating the legs proximal to the bore wall. Thermal models that are used to predict borehole heat exchange are characterized by either simplified formulations that are restrictive in their application, but utilitarian, or complex multi-dimensional analyses that are cumbersome to implement.

Methods

The borehole resistance methodology presented here offers a straightforward solution that is built on single loop conduction shape factor analysis with thermal shunt accounting and pipe-pipe configuration analysis, to extend to multi-loop borehole configurations and custom kidney extrusions.

Results

The borehole resistance predictions are compared to published data and information listed by manufacturers of multi-loop products in third party thermal tests against standard loops. The results are found to agree within the constraints posed by the model assumptions. The methodology offers a straightforward solution that can be incorporated into popular geothermal loop sizing software such as GLD, GLHEPRO and other system software.

Conclusions

The advantages and challenges of these advanced loop designs are discussed and conclusions drawn. By reducing the bore resistance, one can take advantage of less drilling and proportionally less capital cost for the bore field, while achieving the same loop temperatures.

Background

One of the challenges to the ground source heat pump market is the cost burden of implementing the ground heat exchange. For the case of closed-loop vertical borehole designs, the thermal barrier presented by a loop can be significant. The high-density polyethylene (HDPE) pipe loop contacts the bore wall via a thermal grout. For example, a 7–8 °C temperature difference between the water circulating in the loop and the bore wall reduces the heat pump operating performance. In general, the loop legs are not constrained within the bore cross section, and during installation, the legs can twist and contact one another over a substantial portion of the bore length, thereby degrading performance. Typically, the loop is sealed to the bore wall with a bentonite grout mixture, of which the hydraulic and thermal properties are very important to insure optimal performance of the installation. The use of a thermally enhanced grout, while advantageous to heat transfer between the loop and the bore wall, presents its own challenges. For example, the addition of sand to enhance conduction also increases thermal shunting between the legs of the loop and presents challenges in uniform mixing and deployment through a tremie pipe.

Acknowledging these challenges, a study was undertaken to understand the potential benefit of multi-loop and custom extrusion designs which offer greater surface area for heat transfer, while locating the pipe proximal to the bore wall. This borehole design strategy could reduce the capital cost of the ground heat exchanger by reducing the amount of drilling required to serve a load. The contractor or design engineer has the option of either reducing the amount of drilling, while achieving the same end-of-season temperature as a standard loop, or retaining the same well depth and achieving better heat pump seasonal performance. The first approach reduces capital cost of the loop field but retains the same annual savings, while the second does nothing to lower capital cost but rather improves annual utility savings. For example, the building owner might be more interested in lowering the loop field installation cost by 20 % and accepting less annual operating savings. Another strategic advantage might derive from less land utilization and associated mess.

The purpose of this publication is to provide an analytical basis for calculating borehole thermal resistance (BTR) for the case of multiple loops and custom extrusions. The desire is to incorporate the borehole model into existing well field sizing software packages, such as GLD, GLHEPRO, and system modeling software, such as TRNSYS (Klein et al. 1996; Shonder & Hughes 1997) and HyGCHP, on the market. To accomplish this, the thermal analysis of a single loop (Koenig and Helmke 2014) was chosen as the starting point, acknowledging the contributions of others and developing toward a practical thermal model of a borehole with multiple parallel loops.

Building on the quasi-steady-state modeling (Eskilson 1988) for vertical ground heat exchange, there have been a numerous publications on borehole resistance for a single loop. Previous studies employed finite element numerical techniques (Sharqawy 2009; Sagia 2012; Liao 2012). While these are rigorous, they have also proven to be cumbersome to construct. Others have proposed analytical solutions (Hellström 1991; Gu 1998; Philippe 2010) that are readily applied but are restricted to a limited number of pipe geometries. While these studies support the understanding that increasing the grout thermal conductivity and locating the pipe legs proximal to the bore wall enhances borehole performance, it was left to Redmund (1999) and Beier (2012) to actually demonstrate this fact in field tests. A finite element model for the case of two loops in a borehole was published (Diersch 2010), following the work of Al-Koury et al. (2005). Finally, a recent comparative study of helical pipe and triple U-tube loop designs in foundation pile (Zarrella et al. 2013) is included here, as pertinent to this study.

Borehole thermal resistance for a single loop

A cross-sectional model of the borehole is shown in Fig. 1. The BTR includes contributions from convection inside the pipe (downcomer and upcomer legs), through the pipe walls and through the media that fills the space between the pipes to the surrounding bore wall. The equivalent resistor network model is given in Fig. 2, where the water temperatures inside the pipe loop legs are designated T 1 and T 2, and the surrounding wellbore surface T b. To simplify the analysis, symmetry is assumed in which the pipe spacing is equidistant from the center line of the bore. The resistance for each pipe, R p, is composed of convection at the inside wall and conduction through the pipe wall. The shunt resistance, R s, addresses heat transfer between the pipes through the media, while R f defines the grout filler resistance between pipe surfaces to the surrounding bore wall.

Fig. 1
figure 1

Conduction heat transfer from a pipe loop to the surrounding bore wall

Fig. 2
figure 2

Resistance network model representing well bore heat transfer

The temperatures, in general, can all be considered different, such that T 1 > T 2 > T b for the cooling season or T b > T 2 > T 1 for the heating season.

The three equations that describe the resistive network are as follows:

$$ \left({T}_1 - {T}_{\mathrm{b}}\right) - \left({q}_1 - {q}_3\right)\ {R}_{\mathrm{f}} - {q}_1{R}_{\mathrm{p}} = 0 $$
(1)
$$ \left({T}_2 - {T}_{\mathrm{b}}\right) - \left({q}_2 + {q}_3\right)\ {R}_{\mathrm{f}} - {q}_2{R}_{\mathrm{p}} = 0 $$
(2)
$$ \left({T}_1 - {T}_2\right) - {q}_3{R}_{\mathrm{s}} - \left({q}_1 - {q}_2\right)\ {R}_{\mathrm{p}} = 0, $$
(3)

where q is the heating rate per unit length. The total heat transfer between the pipe loop and the bore wall is:

$$ {q}_1 + {q}_2 = {q}_{\mathrm{t}}=\Delta {T}_{\mathrm{avg}}/\left[\left({R}_{\mathrm{p}} + {R}_{\mathrm{f}}\right)/2\right],\ \mathrm{where}\ \Delta {T}_{\mathrm{avg}}{}^{\circ}\left({T}_{\mathrm{avg}} - {T}_{\mathrm{b}}\right)\ \mathrm{and}\ {T}_{\mathrm{avg}} = \left({T}_1+{T}_2\right)/2. $$
(4)

The total borehole thermal resistance is defined as (R p + R f)/2.

The three equations can be solved simultaneously for the heating rates using matrix methods:

qR = T, defined by the vectors q = {q 1, q 2, q 3} and T = {ΔT 1b, ΔT 2b, ΔT 12}, where ΔT is the temperature difference designated by the subscripts between the pipes and bore wall. The resistance matrix, R, is:

$$ \mathbf{R}=\left|\begin{array}{ccc}\hfill \left({R}_{\mathrm{p}} + {R}_{\mathrm{f}}\right)\hfill & \hfill 0\hfill & \hfill \left( - {R}_{\mathrm{f}}\right)\hfill \\ {}\hfill 0\hfill & \hfill \left({R}_{\mathrm{p}} + {R}_{\mathrm{f}}\right)\hfill & \hfill {R}_{\mathrm{f}}\hfill \\ {}\hfill {R}_{\mathrm{p}}\hfill & \hfill \left( - {R}_{\mathrm{p}}\right)\hfill & \hfill {R}_{\mathrm{s}}\hfill \end{array}\right| $$
(5)

developed from Eqs. (1)–(3). The heating rates, q i, are found by taking the product of the inverse matrix, R 1, with the temperature column vector. The solution requires an iterative approach, where starting values for the ΔTs are refined by operating within the simulation of the borehole and surrounding rock thermal storage model in time. The temporal effects of borehole heat capacity are minor (quasi-steady-state), insuring that little time (approximately a day) is required to elapse until the ΔTs stabilize. Depending on the loop design basis (m/W), the ∆T 12 between the loop legs can be considered to be fixed, e.g., −3.3 or 5.6 °C for the heating and cooling seasons, respectively. For the parameters listed in Table 1, the matrix solution is found to be q i = {−24.7, −8.03, −4.69} and ΔT = {−7.3, −3.96, −3.3}, respectively, for a loop design based on (q 1 + q 2) = −32.73 W/m for the heating season.

Table 1 Parameters chosen for this study

The relative amount of heat shunted between pipe is calculated from the ratio q 3/(q 1 + q 2). This turns out to be an important and useful number, because it defines what fraction of the pipe surface area is effective in heat transport to/from the bore wall. This will be apparent shortly in the shape factor analysis of individual pipe heat transfer.

The values for the pipe resistance, R p, the pipe to bore wall resistance, R f, and the shunt resistance, R s, are developed here assuming steady-state conditions. The pipe resistance is composed of series resistances due to convection at the interior wall from water heat transfer and pipe wall contribution. These can be calculated as follows:

$$ {R}_{\mathrm{p}} = 1/\left(\pi {d}_{\mathrm{i}}{h}_{\mathrm{cv}}\right) + \mathrm{L}\mathrm{N}\left({d}_{\mathrm{o}}/{d}_{\mathrm{i}}\right)/\left(2\pi {k}_{\mathrm{p}}\right), $$
(6)

where the d subscripts (o, i) reflect outside and inside pipe diameters, and h cv is the convection coefficient for water flowing inside the pipe, and k p, the pipe material thermal conductivity. Note that d o/d i = 1/(1 − 2/DR), where DR is the dimension ratio = d o/t, where t is the thickness of the pipe wall.

Similarly, the thermal shunt resistance, R s can be modeled as pipe to pipe conduction shape factor in an infinite medium (Holman 1981), although the borehole is finite.

$$ {R}_{\mathrm{s}} = \cos {\mathrm{h}}^{-1}\left({x}_{\mathrm{pp}}\right)/\left(2\pi {k}_{\mathrm{g}}\right),\ \mathrm{where}\ {x}_{\mathrm{pp}} = \left[2{\left(S\mathrm{p}/d\right)}^2 - 1\right] $$
(7)

Defining dimensionless quantities: d r = d/D and S r = S p/d, the parameter, x pp, can be written as:

$$ {x}_{\mathrm{pp}} = \left[2{S_{\mathrm{r}}}^2 - 1\right]. $$
(7A)

For example, Table 1 lists select parametric values in a typical geo-exchange application:

The grout filler resistance between pipe and bore wall is more complicated. Referring to Fig. 2, the shape factor for heat transfer for a single offset pipe to its surrounding cylinder was compiled by Holman (1981):

$$ {R}_{\mathrm{f}\mathrm{s}} = \cos {\mathrm{h}}^{-1}\left({x}_{\mathrm{f}}\right)/\left(2\pi {k}_{\mathrm{g}}\right),\ \mathrm{where}\ {x}_{\mathrm{f}} = \left({D}^2 + {d}^2 - S{\mathrm{p}}^2\right)/(2Dd). $$
(8)

where D is the bore diameter, d the pipe exterior diameter, S p the center to center (c-c) spacing between pipes, and k g the grout thermal conductivity. Again, this can be defined in dimensionless form as:

$$ {x}_{\mathrm{f}} = \frac{1}{2}\ \left[{d}_{\mathrm{r}}\left(1 - {S_{\mathrm{r}}}^2\right) + 1/{d}_{\mathrm{r}}\right]. $$
(8A)

For the selected parameters in Table 1, for example, the pipe to bore resistance is 0.219 m°C/W. Equation (8), however, is based on a single off-centered pipe conducting to an unobstructed surrounding wall; therefore, one might anticipate that the presence of a second pipe (leg of the loop) would contribute a similar measure, so that the total would simply be half of the single pipe contribution, as suggested by the resistance mesh model given earlier. Unfortunately, this is not quite right in that there is conductive interference due to the presence of the pipe representing the second leg of the loop.

This interference or obstruction is precisely addressed by including the factor, f, defined as:

$$ f = 1-{q}_3/\left({q}_1+{q}_2\right), $$
(9)

so that R f = R fs/f in the denominator of the pipe to bore resistance, R fs, to correct the thermal influence of each pipe. The total borehole resistance, in this case, is 0.172 m°C/W. This is composed of parallel contributions of each pipe individually: 0.089 + 0.219/(1 − 0.143) = 0.345 m°C/W, including the pipe and the grout. For the loop, this value is halved to yield the total borehole resistance stated. Note that the pipe legs are at different temperatures, resulting in the thermal shunt of 14.3 %.

A closed-form solution for the pipe loop to bore resistance, R pb, develops from solving:

$$ \mathrm{R}\mathrm{I}=1/{R}_{\mathrm{f}} = 2\pi {k}_{\mathrm{g}}f/ \cos {\mathrm{h}}^{-1}\left({x}_{\mathrm{f}}\right), $$
(10)

where f = 1 − ΔT/(q t R 3), ΔT is the loop delta-T, q t the total heating rate, and R 3 = (2R p + R s) + (R p R s)/R f.

The latter results from solving the matrix equations for q 3 in terms of the heating rate difference (q 1 − q 2).

The quadratic solution is:

$$ \mathrm{R}\mathrm{I}=-\frac{1}{2}\left(a/b-1/{R}_{\mathrm{f}}\right)\ \left\{1\ \hbox{--}\ \mathrm{SQRT}\left[1-4\left(c/b-a/b\ {R}_{\mathrm{x}}\right)/{\left(a/b-1/{R}_{\mathrm{f}\mathrm{s}}\right)}^2\right]\right\}, $$
(10A)

where a = (2R p + R s), b = R p R s, c = ΔT/(q t R 3), R fs and R 3 as defined in (8) and (10), respectively.

The final borehole resistance (BTR) for a single loop has a closed-form solution that is summarized by:

$$ {R}_{\mathrm{b}} = \left({R}_{\mathrm{p}}+{R}_{\mathrm{f}}\right)/2,\ \mathrm{where}\ {R}_{\mathrm{f}} = 1/\mathrm{R}\mathrm{I}. $$
(11)

The benefit of the closed-form quadratic solution is that it avoids the need for the T information, which requires the full bore simulation, by providing a steady-state value for R f straightaway. The two approaches give the same result, with the matrix solution offering additional information on q.

Comparison to others

In prior work (Koenig and Helmke 2014), which provided a comparison to other reported models, a similar approach was employed to assess single-loop BTR, with the difference being in the use of a geometric obstruction factor, f(β). The comparison, while quite good for larger separations, was shy of reported results for cases in which the two pipes were close to one another. As Fig. 3 demonstrates in a comparison of obstruction factors, the discrepancy is confined to the use of a geometric obstruction factor for pipe separations within 35 % of the bore diameter.

Fig. 3
figure 3

Comparison of obscuration factors for this model (blue) vs. geometric (red)

A comparison of this improved model to the results of Sharqawy, shown in Fig. 4, clearly shows a very good match over the full extent of the pipe separation study. Note, at S p/D ≈ 0.28, the bore resistance experiences a jump as the pipe legs actually contact. This result lends credence to the model proposed here.

Fig. 4
figure 4

Comparison of borehole resistance models

Figure 5 demonstrates the benefit of employing higher conductivity grout (reported in W/m°C) as well as increased pipe separation in lowering borehole resistance. Even though thermal shunting increases with higher grout conductivity, the impact of greater pipe conduction to the bore wall tends to dominate the response. In all cases, however, there is a 25 % net BTR benefit to locating the pipe proximal to the bore wall.

Fig. 5
figure 5

Borehole resistance study as a function of grout thermal conductivity

Methods

Extension to multi-loop borehole geometries

In principle, adding additional loop(s) to a borehole should decrease borehole resistance by increasing the available pipe surface area for heat transfer. It should be noted, however, that the heat transfer effectiveness of an individual loop diminishes with each additional loop. This is because the heat transfer surface area of an individual loop for exchange with the bore wall is occulted by the presence of additional loops.

The method delineated here extends the single-loop BTR model by first defining the boundary in the borehole cross section that separates downcomer and upcomer pipe legs of a common temperature. Pipe loops are assumed to be connected in parallel. The thermal interaction between individual pipes across the boundary is treated independently and, without assigning legs of a common loop, using the single loop methodology presented previously. The total borehole resistance is assembled by summing the individual pipe-pipe resistances of a common spacing weighted by their geometric configuration factor. The configuration factor is developed by noting the number of occurrences of specific pipe-pipe geometries within the borehole. An example is offered to help clarify the procedure.

Fig. 6
figure 6

Borehole resistance (BTR) for 2U design as a function of relative loop leg spacing (Z/D)

Starting with the two-loop borehole, shown in Fig. 8 leftmost diagram, the values of R f are calculated for individual pipe on one side of the boundary line to the two pipes on the other side (diagonal and adjacent configurations). The resistances, R p, R pp, and Rf, are calculated for each configuration. The final (Fig. 6) BTR is composed of the contributions of each configuration. The example in Table 2 for D = 152 mm, d = 33 mm, and k g = 0.78 W/m°C shows the assemblage of the final BTR value from equal contributions of each configuration.

Table 2 Two-loop borehole resistance model results

To conform to the borehole, the ordered pipe planform geometries become square (2U), hexagonal (3U), and octagonal (4U), as indicated in Fig. 8. Pragmatically, in order to accommodate these configurations, the outside diameter of the assembly must remain within a maximum envelope, e.g., 114 mm in a 152-mm bore, or 75 % of the opening. To accomplish this, pipe sizes must shrink from 42 mm (1 1/4 in.) to 33 mm (1 in.) down to 27 mm (3/4 in.) as the number of loops increases. The constraints on allowable geometry impose restrictions on the range of pipe separations that are feasible. Moreover, each of these multi-loop designs requires an assembly to affix the pipe to the desired geometry.

Fig. 7
figure 7

Four-loop borehole configuration

As an example, Fig. 7 shows four loops (eight pipes) arranged evenly around a hollow central form of diameter, d cyl.

The maximum (even) number of pipes N x = 2N L that fit around the form is limited by:

$$ d\ \le\ \left({d}_{\mathrm{cyl}} + d\right)\ \sin \left(\pi /{N}_{\mathrm{x}}\right), $$
(12)

where d cyl is the outside diameter of the form. In addition, there is a practical constraint based on the maximum diameter of the loop assembly that the borehole can comfortably accept,

$$ {d}_{\mathrm{cyl}} + 2d\ \le\ 0.75D $$
(12A)

where 75 % is deemed a reasonable limit. Combining (12) and (12A) leads to the condition:

$$ \left[0.75\left(D/d\right) - 1\right]\ \sin \left(\pi /{N}_{\mathrm{x}}\right)\ \ge \kern0.5em 1 $$
(12B)

So, for example, with a 152-mm (6 in.) bore and a 27-mm (3/4 in.) pipe size, the maximum number of pipes that fit the constraint is ten or, equivalently, five loops. The outside diameter (OD) of this assembly is then (d cyl + 2d), or, in this example, 114 mm (4.5 in.). Other sizes are possibly limited only by the assembly size impact on insertion in the borehole and the ability to insert a grout “tremie” pipe inside the ID of the form, or, ideally, use the form itself as the grout pipe.

To extend the previous model to the case of multi-pipe borehole geometries, under the assumption that the form itself offers negligible thermal isolation, one needs to simply specify the number of pipe loops, the diagonal pipe spacing, and then calculate the pipe-pipe shunt resistance for each configuration. The latter is based on the pipe-pipe spacing for each configuration, as shown in Fig. 8.

Fig. 8
figure 8

Compilation of thermal separations between pipes for two-, three-, and four-loop assemblies

The horizontal line delineates pipe of a common temperature (downcomer vs. upcomer legs). Table 3 summarizes the spacing and number of configurations for each loop geometry, where Z is the fixed diagonal c-c spacing. From this table, the shunt spacing for each separate configuration is identified and the shunt resistance, R s, can be computed for each configuration separately. This value, then, is inserted into the borehole resistance matrix to derive the obstruction factor, f.

Table 3 Multi-loop pipe-pipe spacing for shunt analysis

The final borehole resistance for the multi-loop assembly is calculated from:

$$ \mathrm{B}\mathrm{T}\mathrm{R} = \varSigma\ \left[{R}_{\mathrm{bi}}{\left(\mathrm{C}\mathrm{F}\right)}_i\right]/{N}_{\mathrm{L}} $$
(13)

where R b is the borehole resistance for a single (Z) configuration, multiplied by CF, its configuration fraction, all summed over the number of configurations from i = 1 to N L.

One of the advantages of additional loops is seen in the “N L” factor that appears in the denominator which decreases borehole resistance as loops are added. An additional benefit is that the pipe wall thickness, assuming a fixed DR for pressure rating, decreases as the pipe gets smaller to accommodate the loops.

Results

Using the method outlined in this paper, a compilation of multi-loop borehole resistance is provided in Fig. 9 for specific loop designs and values of the pipe spacing ratio, d/D. The reduction in borehole resistance as additional loops are added is dramatic. This is especially noticeable with the addition of a second loop.

Fig. 9
figure 9

BTR at k g = 0.78 W/m°C as a function of relative pipe spacing and loop designs

By constraining the OD of the multi-pipe loop assembly to 75 % of the bore diameter, the study considered the benefit of additional loops. Fig. 10 shows that the borehole resistance decreases dramatically and was found to be inversely proportional to the number of loops to the power of 0.86.

Fig. 10
figure 10

BTR decrease as loops are added to the borehole

Figure 10 compares BTR values of multiple loop designs reported by specific loop manufacturers to the predictions of this model. In addition to the single-loop match (Fig. 4), there are two additional points of validation:

  • From thermal measurements by Bowman Geothermal of REHAU 1-in. double-loop product in side-by-side installations, they report a BTR of 0.11 (PEX2), a 40 % reduction from single loop, which, after adjustment for HDPE, matches the model with an RMS error of 1.5 %.

  • Agreenability (2014) report a BTR of 0.049 m°C/W in an application of their Twister™ four-loop 3/4-in. pipe technology, which compares favorably to predictions exhibited in Fig. 9.

At some point with addition of multiple loops, the borehole resistance model effectively transitions to a central pipe equivalent diameter, d eq, with a grout gap defined by:

$$ {R}_{\mathrm{b}} = {R}_{\mathrm{p}}/2{N}_{\mathrm{L}} + \mathrm{L}\mathrm{N}\left(D/{d}_{\mathrm{eq}}\right)/\left(2\pi {k}_{\mathrm{g}}\right), $$
(14)

where (Z + d) ≤ d eq ≤ (2N L − 2)d for N L = 2, 3, or (2N L − 1)d for N L = 2. The equivalent diameter is the wetted pipe perimeter limited between the smooth OD of the assembly and the exposed scalloped perimeter of the individual pipe. In this case, the pipe wall between pipes in the assembly becomes adiabatic, except for the end pipe that are at different temperatures.

If the assembly was designed such that the adjacent pipe in the assembly was at different temperatures, then, of course, this would diminish performance because of thermal shunting. Rather, if the assembly could be arranged to segregate pipe of the same temperature (downcomers and upcomers), as suggested in Fig. 11a, then, this might indeed improve performance. An alternate loop design, shown in Fig. 11b, eliminates the intermediate adiabatic walls between pipe in Fig. 11a, in favor of a custom kidney extruded shape. This design by ARB Geowell (2014) eliminates the need for an assembly of smaller pipe in favor of coiled pipe halves that must be joined with a molded transition from kidney pipe to 1 1/2-in. round pipe for connections.

Fig. 11
figure 11

a Segregated pipe assembly. b Extruded kidney pipe deployment

The utility of the proposed borehole thermal model is such that it can be easily extended from simple loops to more complex multi-pipe and custom shape loop configurations. The model addresses only the thermal resistance inside the borehole (BTR), and not the diffusion component into the surrounding formation. The benefit is that it can be incorporated easily into existing software, such as GLHEPRO, to handle multi-loop borehole configurations, avoiding the details of constructing custom numerical models through FEFLOW.

The model produces reasonable results that support the data published by others. The extension of the model to multiple pipe loops in a single borehole is accomplished by (1) defining the mid-plane boundary between downcomer and upcomer pipe, (2) calculating the borehole resistance for each pipe-pipe separation across the boundary, (3) multiplying by the relative number of occurrences of each configuration to yield an average borehole resistance for a representative loop, and finally (4) dividing by the number of loops. Since the loop to bore wall temperature difference is controlled by the product, q(R b), it is noteworthy that it makes no difference whether this is treated as a single loop, in which case (q/N L)R b, or the full multi-loop design, q(R b/N L).

The benefit and practicality of implementing multiple loop borehole designs to improve performance and/or reduce installation cost is arguable. The use of a thermally enhanced grout with a single loop appears to lower borehole resistance by 30 % over neat bentonite grout at standard U-bend separations, despite higher thermal shunting between the legs of the loops. But, unless the pipe legs can be guaranteed to be separated and not twisted and contacting over a substantial length of the bore, this 30 % benefit may never be realized in practice. Physical separation, such as the use of Geo-Clips, is at least as beneficial as thermally enhanced grout. With this in mind, the obvious benefit of multiple loops in reducing borehole resistance is increased surface area for heat transfer and thinner HDPE walls that result from employing smaller diameter pipe loops. Another advantage to multiple loops is in lowering the pressure drop of the ground heat exchanger, which can be as much as 35 %. This benefit develops from dividing the flow into multiple pipe loops of a given cross section as well as shortening the length. In operation, this can save considerable pump energy over the year.

With that said, the question remains whether the assembly of pipe loops can be cost-effective with the benefit of less required drilling. On paper, the savings in drilling more than offsets the extra cost of multiple pipe assembly, which ultimately means that plastic is less expensive than the lifecycle operation of the drill rig. Certainly, there are excellent targets for commercialization of new loop technology, such as commercial applications where space for the loop field is always an issue, and the opportunity to repurpose abandoned water wells that could serve residential customers without extending the depth of the well. Ultimately, driller pricing and competition will prove the winner.

Conclusions

A borehole thermal model for vertical loops was developed from a modified shape factor analysis. An interference factor, f, defined by (1 − q shunt/q), the relative size of the thermal shunt between the loop segments, was employed as the modification required to properly model borehole thermal resistance. The model predictions were compared to reported results and found to match quite nicely over a range of pipe sizes and spacing. Following this success, the model was extended to multi-pipe loop geometries consisting of two-, three-, and four-loop assemblies in a single borehole and a custom kidney extrusion. The borehole resistances were calculated for each of these configurations and presented as a compilation of predicted performance. These predictions were compared to the available data reported by manufacturers of multi-pipe loop assemblies and found to give reasonable results.

Abbreviations

T p :

temperature of pipe loop segment, °C

T b :

temperature of borehole perimeter, °C

ΔT 1n :

temperature difference pipe 1 to (n = b bore; n = 2 pipe 2), °C

q n :

exchange heating rate (pipe to bore wall or pipe to pipe), W/m

R p :

thermal resistance of pipe (convection and pipe wall), m°C/W

R fs :

single pipe thermal resistance (pipe to bore wall), m°C/W

R f :

effective thermal resistance (pipe to bore wall), m°C/W

D :

borehole diameter, m

d :

outside diameter of the loop pipe, m

DR:

dimension ratio: pipe diameter to wall thickness

S p :

center to center spacing between pipe, m

h cv :

convection heat transfer coefficient, W/m2°C

k p :

HDPE pipe thermal conductivity, W/m°C

R s :

pipe-pipe shunt thermal resistance, m°C/W

k g :

grout thermal conductivity, W/m°C

d r :

diameter ratio = d/D

S r :

spacing ratio = S p/d

X pp :

pipe to pipe dimension = [2(S p/d)2 − 1]

X f :

pipe to bore wall dimension = 1/2 [d r(1 − S r 2) + 1/d r]

f :

loop pipe interference factor in thermal exchange

d cyl :

central cylinder diameter, m

N L :

number of loops in borehole

CF:

configuration fraction

R b :

borehole thermal resistance (BTR), m°C/W

d eq :

equivalent diameter of pipe loops actively exchanging heat, m

Z :

diagonal c-c distance between pipe in multi-loop geometry, m

References

  • Al-Koury R, Bonnier PG, Brinkgreve RB (2005) Efficient finite element formulation for geothermal heating systems. Part I: steady-state. Int J Numerical Methods Eng 63(7):933–1013, Pt. II: Transient, IJNME 67 (2006) 5, 725–745

    Google Scholar 

  • Agreenability, J R (2014) TWISTER high efficiency geothermal heat exchanger. IGSHPA Conference, Baltimore, MD

    Google Scholar 

  • ARB Geowell, Koenig A (2014) A new HDPE pipe loop for single well residential opportunity. Paper presented at the PGWA winter conference, Grantsville, PA, 23–24 January 2014

  • Beier RA, Ewbank GN (2012) In-situ thermal response test interpretations: OG and E ground source heat exchange study. IGSHPA, OKSU

    Google Scholar 

  • Crecraft H, Mahlmann J, Rehau in-situ thermal performance testing Bowman Geothermal, LLC, Chantilly, VA: http://www.rehau.com/download/866550/borehole-thermal-conductivity-report.pdf

  • Diersch HJ, Bauer D, Heidemann W, Rühaak W, Schätzl (2010) Finite element formulation for borehole heat exchangers in modeling geothermal heating systems. FEFLOW, DHI-WASY software, White Papers Volume V. 5:5–96.

  • Eskilson P, Claesson J (1988) Simulation model for thermally interacting heat extraction boreholes. Numerical Heat Transfer 13(2):149–165

    Article  Google Scholar 

  • GLD: www.seventhwave.org/hygchp

  • GLHEPRO: https://hvac.okstate.edu/glhepro/overview

  • Gu Y, O’Neal D (1998) Development of an equivalent diameter expression for vertical U-tubes used in ground-coupled heat pumps. ASHRASE Trans 104(2):9

    Google Scholar 

  • Hellström G (1991) Ground heat storage, thermal analyses of duct storage systems, doctoral thesis. Department of Mathematical Physics, University of Lund, Lund, Sweden

    Google Scholar 

  • Holman JP (1981) Heat transfer. McGraw-Hill, New York, pp 68–71

    Google Scholar 

  • HyGCHP: www.ecw.org/hygchp

  • Koenig A, Helmke M (2014) Development of a thermal resistance model to evaluate wellbore heat exchange efficiency. Int J Energy Environ 5(3):297–304

    Google Scholar 

  • Liao Q, Zhou C, Cui W, Jen T (2012) Effective borehole thermal resistance of a single U-tube ground heat exchanger. Numerical Heat Transf 62:197–210

    Article  Google Scholar 

  • Phillippe M, Bernier M, Marchio D (2010) Sizing calculation spreadsheet vertical geothermal borefields. ASHRAE J 7:20–28

    Google Scholar 

  • Redmund CP (1999) Borehole thermal resistance: laboratory and field studies. ASHRAE Trans 105:439–445

    Google Scholar 

  • Sagia Z, Stegou A, Rakopoulos C (2012) Borehole resistance and heat conduction around vertical ground heat exchangers. The Open Chem Eng J 6:32–40

    Article  Google Scholar 

  • Sharqawy MH, Mokheimer EM, Hassan MB (2009) Effective pipe-to-borehole thermal resistance for vertical ground heat exchangers. Geothermics 38:271–277

    Article  Google Scholar 

  • TRNSYS: sel.me.wisc.edu/trnsys/downloads/trnsedapps/gchpdemo/overview.htm; Klein, S.A., et al., TRNSYS Manual, “A Transient Simulation Program”, Solar Energy Laboratory, University of Wisconsin, Version 14.2, for Windows, Sep 1996.

  • Shonder, John and Hughes (1997) Patrick from ORNL (Oak Ridge National Laboratory). Richard Stockton College.

  • Zarrella A, DeCarli M, Galgaro A (2013) Thermal performance of two types of energy foundation pile: helical pipe and triple U-tube. Appl Therm Eng 61(2):301–310

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to AA Koenig.

Additional information

Competing interests

The author declares that he has no competing interests, although ARB Geowell owns the patent for K-Pipe, a kidney pipe custom extrusion. The performance benefit of this pipe was treated objectively within the context of multiple pipe loops. The author has no financial commitment from the company to support the publication of this work nor derives a benefit from its publication.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Koenig, A. Thermal resistance of borehole heat exchangers composed of multiple loops and custom shapes. Geotherm Energy 3, 10 (2015). https://doi.org/10.1186/s40517-015-0029-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40517-015-0029-1

Keywords