Next Article in Journal
Enhancement of Hosting Capacity with Soft Open Points and Distribution System Reconfiguration: Multi-Objective Bilevel Stochastic Optimization
Next Article in Special Issue
Residential Buildings’ Foundations as a Ground Heat Exchanger and Comparison among Different Types in a Moderate Climate Country
Previous Article in Journal
Partial Discharge Detection Based on Anomaly Pattern Detection
Previous Article in Special Issue
Reviewing the Modeling Aspects and Practices of Shallow Geothermal Energy Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Explicit Multipole Formula for the Local Thermal Resistance in an Energy Pile—The Line-Source Approximation

1
Building Physics, Lund University, 221 00 Lund, Sweden
2
Building Services, Lund University, 221 00 Lund, Sweden
*
Author to whom correspondence should be addressed.
Energies 2020, 13(20), 5445; https://doi.org/10.3390/en13205445
Submission received: 22 September 2020 / Revised: 13 October 2020 / Accepted: 15 October 2020 / Published: 19 October 2020
(This article belongs to the Special Issue Computational Geothermal Energy Applications)

Abstract

:
This paper presents a closed-form quite handy formula for the local thermal resistance Rb between the temperature of the bulk heat-carrier fluid in the pipes, equally spaced on a concentric circle inside a circular energy pile, and the mean temperature at the periphery of the pile. The so-called multipole method is used to calculate the temperature field. An important improvement of the multipole method is presented, where Cauchy’s mean value theorem of analytical functions is used. The formula for thermal resistance Rb0 for the zero-order approximation (J = 0), where only line heat sources at the pipes are used, is presented. The errors using zeroth-order approximation (J = 0) are shown to be quite small by comparisons with eight-order approximation (J = 8) with its accuracy of more than eight digits. The relative error for the local thermal resistance Rb0 for the zero-order approximation (J = 0) lies below 5% for a wide range of input parameter values. These ranges are judged to cover most practical cases of application. The smallest local thermal resistance Rbmin is, with some exceptions, obtained when the pipes lie directly in contact with the pile periphery. A neat formula for this minimum is presented.

1. Introduction

Ground-source heat pump (GSHP) systems are among the most energy-efficient and environmentally clean heating and cooling technologies available today [1]. Due to their higher efficiency, they offer substantial potential for energy savings and carbon emissions reduction in buildings [2]. The most common application of GSHP systems is with vertical borehole heat exchangers [3]. Boreholes are routinely drilled to 100- to 400-m depths and require a specific drilling area as well as specialized drilling equipment. The high capital cost of installing borehole heat exchangers is one of the underlying reasons that hamper the widespread deployment of GSHP systems in many places around the world. For new constructions, an alternative to using borehole heat exchangers is to use piled foundations with embedded heat exchange elements as ground heat exchangers. The thermally active foundation piles, commonly referred to as energy or thermal piles, are dual-function structural elements that transmit mechanical loads through unstable soil layers to lower ground levels with high bearing capacity as well as exchange heat with the surrounding ground to provide heating and cooling to the building. They do not require the capital cost and land requirements of borehole heat exchangers and have the potential to lower the environmental footprint of GSHP systems by using fewer materials with low embodied energy.
Energy piles are constructed in many different configurations, as shown in Figure 1, and have typical diameters varying between 0.1 m and up to 3.0 m and lengths ranging from 10 to 60 m [4]. The number and arrangement of heat transfer pipes in the pile depend on the size and construction method of the energy piles. Driven piles of steel or concrete usually have one to two U-pipes arranged in a series or parallel configuration. The heat transfer pipes are placed towards the center of the pile in the hollow driven piles and are positioned around the reinforcement bars in the solid driven piles. Continuous flight auger (CFA) piles have one to three U-pipes installed in the center of the pile. Bored piles have one to several U-pipes mounted on the reinforcement frame close to the boundary of the pile. The differences in construction methods and geometrical arrangements lead to different thermal characteristics of energy piles.
The current design and analysis methods of energy piles are primarily based on the methods developed for borehole heat exchangers [5]. This is done under the assumption that the thermal characteristics of energy piles and borehole heat exchangers are similar, which is not entirely accurate but simplifies the analysis. Energy piles usually have a greater number of heat transfer pipes and a far smaller aspect (length to diameter) ratio, which differentiates their thermal response from borehole heat exchangers. A key characteristic in this regard is the thermal resistance, which characterizes the heat transfer inside the ground heat exchangers. The thermal resistance Rb of a ground heat exchanger is the ratio of the temperature difference between the circulating fluid in the heat transfer pipes and the surrounding ground to the heat flux per unit length of the ground heat exchanger. The larger number of pipes and their many possible arrangements in energy piles complicate the heat transfer and make the calculation of thermal resistance calculation more complex than for borehole heat exchangers.
There are two main approaches for calculating the thermal resistance of energy piles. In the first approach, thermal resistance is estimated using in situ thermal response testing of energy piles [6,7]. This approach is not generally useful for design purposes as it can only estimate thermal resistance of an existing energy pile. In the second approach, thermal resistance is computed using mathematical models based on conductive heat transfer. Several mathematical models have been proposed and used in this regard. Readers are referred to studies [8,9,10] for a detailed survey and comparison of these models. Some of the commonly used models for calculating the thermal resistance of energy piles include those of Loveridge and Powrie [11], Paul [12], Sharqawy et al. [13], Bauer et al. [14], and Diao et al. [15]. A great majority of these models only work well under certain conditions and perform poorly under others [10]. Moreover, among the existing methods, only the Loveridge and Powrie model can calculate the thermal resistance of energy piles with more than two U-pipes. In the absence of suitable analytical methods, the common practice has been to use a conservative estimate of thermal resistance when sizing energy piles with three or more U-pipes.
One of the most robust and accurate methods for calculating the thermal resistance of ground heat exchangers is the multipole method. It is a mathematically rigorous and thoroughly validated analytical method used for calculating the thermal resistances of any number of arbitrarily placed pipes in a composite region. The model, originally developed in the 1980s by Bennet, Claesson, and Hellström [16], uses a combination of line heat sources and so-called multipoles to yield an exact solution of the steady-state temperature field and thermal resistances for arbitrarily placed pipes in a circular region surrounded by an infinite region. The accuracy of the model increases with the number of multipoles used for the calculation. The 10th-order multipole provides a greater than 8-decimal digit accuracy. The original model with 10th-order multipole was provided as a FORTRAN code of approximately 600 lines. It has since been incorporated into several building energy simulation software and ground heat exchanger design programs for calculating the thermal resistance of borehole heat exchangers with single and double U-pipes [17,18]. A principal modification to the model was made in 2011 by Claesson and Hellström [19], who simplified the algorithm by putting a condition on the mean temperature around the periphery of the ground heat exchanger instead of prescribing a constant temperature condition at a suitable radial distance outside the ground heat exchanger in the original model. The modified multipole model is available as a Mathcad program.
The multipole method has a highly sophisticated mathematical formulation and a rather complex algorithm. Its implementation in design and simulation programs is intricate and requires considerable programming expertise and effort. This has led to the development of closed-form multipole formulas that, under certain conditions, can be as accurate as the original multipole method in calculating the thermal resistance of ground heat exchangers. So far, exact multipole formulas have only been developed for single and double U-pipe configurations. Claesson and Javed [20] presented zeroth-, first-, second-, and third-order multipole formulas for single U-pipe ground heat exchangers with symmetrically positioned pipes. Later, the authors also presented zeroth- and first-order multipole formulas for double U-pipe ground heat exchangers with symmetrically positioned pipes [21]. The accuracy of the presented formulas depends upon the order of their multipoles. The zeroth-order multipole formulas only account for line sources and do not consider any multipoles. Javed and Spitler [10] noted that for a single U-pipe borehole heat exchanger the accuracy of the zeroth-order multipole formula varies from case to case. In some cases, the accuracy of the zeroth-order formula is quite good; in other cases, it is not as good. However, the first- and second-order multipole formulas provide an accuracy of better than 3% and 1%, respectively, compared to the original multipole method under all practical conditions [20]. Similar trends have also been observed for double U-pipe ground heat exchangers by Claesson and Javed [21]. The multipole formulas for single and double U-pipe ground heat exchangers are also applicable to energy piles with a corresponding number of U-pipes.
In this paper, an important improvement of the multipole method is presented. Cauchy’s so-called mean-value theorem for analytical functions is used to determine the value of certain complex-valued integrals around the pipe periphery (see Section 3). The objective of this paper is to derive a closed-form zeroth-order multipole formula for any number of evenly distributed pipes in a ground heat exchanger. The derived formula is applicable for energy piles with multiple U-pipes equally spaced on a circle. The accuracy of the presented formula is established by a detailed parametric study that brackets all or almost all real-world energy piles with a circular cross-section.

2. Energy Pile with Pipes Equally Spaced on a Circle

Figure 2 shows the considered type of energy pile. The N (=8) vertical heat exchanger pipes are equally spaced on the circle r = rc. Typical values for N may be in the range from 2 to 12. The left-hand figure shows a horizontal cross-section of the energy pile and the surrounding ground, and the right-hand figure shows in greater detail the pipe with the center at (rc, 0) on the positive x-axis.
The radius of the energy pile is rb. The index b for borehole is kept from other multipole papers. The thermal conductivity in the ground outside the energy pile, r > rb, is λ (W/(m·K)), and it is λb inside the circular energy pile region outside the pipes. The outer radius of the pipes is rp, and Rp denotes the thermal resistance from bulk fluid in the pipes to the outer pipe wall (per meter pipe, (K·m)/W). The pipe resistance Rp is the sum of the thermal resistances of the pipe wall (an annulus) and the fluid boundary layer.
The heat injection from each pipe to the ground is q0 (W/m). The value of q0 is negative for heat extraction from the ground. The required fluid temperature Tf for the heat carrier fluid to achieve the prescribed heat flux q0 is to be calculated.
The (steady-state) thermal problem requires that an outer ground temperature is prescribed. As will be discussed in Section 6, it is sufficient to prescribe the average temperature Tbav around the periphery r = rb of the energy pile. This way is close to the thermal problem that was first presented in [19].
The thermal problem has a high degree of symmetry. The temperature field in the wedge indicated in the left-hand Figure 2 around the x-axis with the opening angle φ = 2π/N is repeated around each of the N pipes. The temperature field is symmetric in y, so the upper and lower halves of the wedge mirror each other.
The primary input data are:
N , q 0 , λ b , λ , r p , r c , r b , R p , T bav .
The thermal resistance Rp ((K·m)/W) may be represented by a nondimensional parameter β, and the ratio between the two conductivities by the parameter σ:
β = 2     π   λ b R p , β 0 , σ = λ b λ λ b + λ , 1 < σ < 1 .
The positions of the N pipes on the circle r = rc become using the complex exponential ze:
z n = r c z e n , z e = e 2 π i N = cos ( 2   π   /   N ) + i sin ( 2   π   /   N ) , n = 1 ,     2 ,       ,     N ; z 1 = r c z e 1 , z 2 = r c z e 2 , , z N = r c z e N = r c ; | z e | = 1 , arg ( z e ) = 2   π N .
The positions of consecutive pipes are obtained by the complex rotation ze or increase of the polar angle by 2π/N. The last pipe N lies on the positive x-axis at z = rc.
The circles around the pipes and the pile must not intersect. This gives the following restrictions on rp, rc, and rb for any N:
N 2 : r p sin ( π   /   N ) r c r b r p ; N = 1 : 0 r c r b r p , N π arcsin ( r p / ( r b r p ) ) .
The lower limit is the case when the pipes are squeezed in the center so that consecutive pipes touch each other. In the upper limit, the pipes touch the wall of the energy pile. The left-hand inequalities give an upper limit for N depending on rp/(rbrp) as indicated to the right. The pipes may be the two shanks of U-pipes and then N is an even number. However, all formulas in the paper are valid for any N including N = 1.
The thermal resistance Rb ((K·m)/W) between the heat carrier fluid in the pipes and the periphery of the energy pile with the average temperature Tbav is defined by the relation:
N     q p R b = T f T bav , R b = 1 K b N     q p = K b ( T f T bav ) .
The inverse of the thermal resistance is the corresponding thermal conductance Kb (W/(K·m)) of the energy pile.

3. The Multipole Method Applied to the Energy Pile

This section presents the application of the multipole method to an energy pile with multiple U-pipes equally spaced on a circle.

3.1. Temperature Field from the N Line Heat Sources

Complex notations and Cartesian and polar coordinates are used:
z = x + i y = r e i φ = r cos ( φ ) + i r sin ( φ ) ,   | z | = r = x 2 + y 2 .
The complex-valued logarithm of any complex function f (z) is defined by:
ln ( f ( z ) ) = ln ( | f ( z ) | ) + i arg ( f ( z ) ) , π < arg ( f ( z ) ) π .
Explicit expressions for the temperature field are derived in Appendix A1, Appendix A2 and Appendix A3. The temperature field has the following form:
T ( x , y ) = T c ( z ) = T bav + q 0 2   π   λ b Re [ W c ( z , r c ) ] , z = x + i y
where Wc (z, rc) is obtained from the sum of the N complex-valued line sources in Figure 2, left. The following compact expressions are derived in detail in Appendix A1, Appendix A2 and Appendix A3, Equation (A12):
W c ( z , r c ) = { ln ( r b N z N r c N ) + σ ln ( r b 2 N r b 2 N z N r c N ) | z | r b   λ b   λ [ N ln ( r b z ) + 2   λ λ b + λ ln ( z N z N r c N ) ] | z | r b .
The temperature field involves the absolute value of two polynomials in z:
P c ( z ) = | z N r c N | , P σ ( z ) = | r b 2 N z N r c N | .
The formulas for the temperature in the pile and ground regions become:
T ( x , y ) = T bav + q 0 2   π   λ b [ ln ( r b N P c ( z ) ) + σ ln ( r b 2 N P σ ( z ) ) ] , | z | r b
T ( x , y ) = T bav + q 0 2   π   λ [ N ln ( r b | z | ) + 2   λ λ b + λ ln ( | z | N P c ( z ) ) ] , | z | r b .
This is the line-source solution without using multipoles. The above temperature field satisfies the following conditions:
  • Laplace equation in the pile and ground regions,
  • The heat flux of q0 (W/m) from each pipe,
  • The internal boundary conditions of continuous temperature at the pile periphery,
  • The internal boundary conditions of continuous radial heat flux at the pile periphery (the thermal conductivity is λb in the pile region and λ in the outside ground),
  • The average temperature of Tbav at the pile periphery.
The temperature field in polar coordinates becomes from (11) and (12):
T c ( r     e i φ ) = T pol ( r , φ ) = T bav + { q 0 2   π   λ b [ ln ( r b N P c ( r   e i φ ) ) + σ ln ( r b 2 N P σ ( r   e i φ ) ) ] r r b N q 0 2   π   λ ln ( r b r ) + q 0 π   ( λ b + λ ) ln ( r N P c ( r   e i φ ) ) r r b .
The two denominators in the logarithms become:
P c ( r   e i φ ) = | r N   e i N φ r c N | = r 2 N 2     r N r c N   cos ( N   φ ) + r c 2 N , P σ ( r   e i φ ) = | r b 2 N r N   r c N   e i N φ   | = r b 4 N 2     r b 2 N   r N r c N   cos ( N φ ) + r 2 N   r c 2 N .
The temperature around the wall of the energy pile becomes:
T pol ( r b , φ ) = T bav + q 0   π   ( λ b + λ ) ln ( r b N P cb ( φ ) ) , P cb ( φ ) = r b 2 N 2     r b N r c N   cos ( N   φ ) + r c 2 N .
P cb ( 0 ) = r b N r c N , P cb ( 0.5   π   /   N ) = r b 2 N + r c 2 N , P cb ( π   /   N ) = r b N + r c N .
The wall temperature varies around the average Tbav roughly following cos(N·φ).

3.2. Boundary Condition at the Pipes

The boundary condition at the periphery of any pipe remains to be fulfilled. It is sufficient to consider pipe N at zN = rc due to the rotational symmetry, where the temperature field is repeated for every increase of the polar angle φ by 2π/N. See Figure 2 and the discussion at Equation (A13).
Let ρ be the local radial distance from the center of pipe N and ψ the local polar angle from this center z = rc as indicated in Figure 2, right:
z = r c + ρ e i ψ , 0 < ρ r b r c , π < ψ π .
The circular region around pipe N inside the energy pile area extends at least to the pipe radius rp. Let hp (W/(m2·K)) denote the heat transfer coefficient (per unit pipe area) between the heat carrier fluid in the pipes and the outer side of the pipe wall. The pipe thermal conductance Kp (W/(m·K)) is obtained by multiplication by the length of the circumference:
K p = 2     π   r p h p , R p = 1 2     π   r p h p , λ b   h p = λ b 2     π   r p R p = β r p .
In the right-hand formula, (2), left, is used. The boundary condition involves a heat balance between the radial heat flux at the outside pipe wall and, the difference between the bulk fluid temperature Tf and the temperature outside the pipe multiplied by the heat transfer coefficient hp. The temperature outside the pipe and the heat flux vary around the pipe with the local polar angle ψ:
h p [ T f T c ( r c + r p e i ψ ) ] = q p N ( ψ ) = ( λ b ) ρ [ T c ( r c + ρ e i ψ ) ] | ρ = r p ,   π < ψ π .
It should be noted that the above two heat fluxes have the dimension of W/m2.
The above exact boundary condition at each point of the periphery may be written in the following way:
T BC ( ψ ) = T c ( r c + r p e i ψ ) + q p N ( ψ ) h p = T f , π < ψ π .
The sum of the temperature outside the pipe and the radial heat flux over the pipe divided by hp shall be equal to (a constant) Tf for all ψ. The subscript BC is used to denote this left-hand boundary-condition function of ψ. In a more precise solution, the boundary condition may be fulfilled with increasing accuracy by adding multipole components to the temperature field with suitably adjusted strength factors. This is, however, left to a sequel paper. Here, only line heat sources are used, which means that the boundary condition is fulfilled approximately as a mean heat balance for the considered pipe.
The integral of qpN (ψ) around the pipe is equal to the heat flux q0. Integration of (20) gives:
π π q p N ( ψ ) r p   d ψ = q 0 , π π T c ( r c + r p e i ψ ) r p   d ψ + 1 h p q 0 = 2     π   r p T f .
Using (18), center, the following relation between the fluid temperature Tf0 and the mean temperature outside the pipe is obtained in the present approximation without multipoles (J = 0):
1 2     π π π T c ( r c + r p e i ψ )     d ψ + R p q 0 = T f 0 .
The subscript 0 is added for the fluid temperature to underline that it is the value for J = 0. The remaining task is now to calculate this mean temperature around the pipe from the temperature field (8) and (9), upper line.

3.3. Mean Temperature at the Outer Periphery of the Pipes

The temperature Tc (z) consists of the line source at z = rc and a remaining part that is an analytic function of z near z = rc. The direct contributions from the N pipes in (9), upper left, is rewritten in the following way:
ln ( r b N z N r c N ) = ln ( r b z r c ) + ln ( r b N z N r c N z r c r b ) = ln ( r b z r c ) + ln ( r b N 1 S ( z , r c ) ) , z N r c N z r c = S ( z , r c ) = n = 1 N z N n r c n 1 , S ( r c , r c ) = N r c N 1 .
The full solution Wc (z, rc) consists of the line source with its singularity at z = rc and an analytic part Wa (z, rc):
W c ( z , r c ) = ln ( r b z r c ) + W a ( z , r c ) , W a ( z , r c ) = ln ( r b N 1 S ( z , r c ) ) + σ ln ( r b 2 N r b 2 N z N r c N ) .
The integral of the line source becomes:
1 2     π π π Re [ ln ( r b r c + r p e i ψ r c ) ]     d ψ = 1 2     π π π Re [ ln ( r b r p e i ψ ) ]     d ψ = ln ( r b r p ) .
There is a nice mean-value theorem for analytic functions by Cauchy, which is directly applicable to the integral of Wa (z, rc) around the pipe at z = rc. Let f (z) be a function that is analytic on and inside a circle with the radius rp around a point z0 in the complex z-plane. Then, the mean value of the function on the circle is equal to the value of f (z) at the center of the circle:
1 2     π π π f ( z 0 + r p e i ψ )     d ψ = f ( z 0 ) , f ( z )   analytic   in     | z z 0 | r p .
This formula may be found in any standard book on analytical functions, for example [22]. An important idea in the improved new approach to the multipole technique presented here is to use this mean-value theorem for the analytic function Wa (z, rc):
1 2     π π π W a ( r c + r p e i ψ , r c )     d ψ = W a ( r c , r c ) .
The integral of Tc around the pipe in (22) gets a logarithmic contribution (25) from the line source and another one (27) from the analytical part:
1 2     π π π T c ( r c + r p e i ψ )     d ψ = T bav + q 0 2   π   λ b [ ln ( r b r p ) + Re [ W a ( r c , r c ) ] ] .
The value of Wa (z, rc) for z = rc is from (23) and (24):
W a ( r c , r c ) = ln ( r b N 1 N r c N 1 ) + σ ln ( r b 2 N r b 2 N r c 2 N ) .
The fluid temperature now becomes from (22), (27), and (29):
T f 0 = T bav + R p q 0 + q 0 2   π   λ b [ ln ( r b r p ) + ln ( r b N 1 N r c N 1 ) + σ ln ( r b 2 N r b 2 N r c 2 N ) ] .

4. Thermal Resistance of the Energy Pile

The main goal of this paper is to derive a formula for the (local, steady-state) thermal resistance between the heat carrier fluid in the pipes and the periphery of the energy pile with its prescribed average temperature Tbav. This resistance becomes from (5) and (30):
R b 0 = T f 0 T bav N   q 0 = R p N + 1 N 2   π   λ b [ ln ( r b N N   r p   r c N 1 ) + σ ln ( r b 2 N r b 2 N r c 2 N ) ] .
The formula is based on the line-source approximation for the temperature field. An important task is then to determine the accuracy of the formula for various input parameter values; see Section 8.

4.1. The Thermal Resistance Rb0 as a Function of rc

The thermal resistance Rb0 of the energy pile depends on the input parameters (1) and in particular on the radial distance rc from the center z = 0 to the pipes. Equation (31) may be written as:
R b = R p N + 1 N 2   π   λ b   [ ln ( r b N   r p ) + R ( r c ) ] .
The dependence on rc is expressed by the function R’(rc):
R ( r c ) = ( N 1 ) ln ( r b r c ) + σ ln ( r b 2 N r b 2 N r c 2 N ) , r c r b r p .
The range of allowed values for rc is limited by the restrictions that the circles of pipes and energy pile must not intersect each other, (4).
An interesting question is now what is the best choice for rc? The thermal resistance shall be as small as possible. The resistance should normally decrease as the distance to the pile periphery decreases. The derivative of R’(rc) is readily obtained from (33):
d R d r c = N 1 r c + σ 2   N r c 2 N 1 r b 2 N r c 2 N = 2   N r c 2 N 1 r b 2 N r c 2 N [ N 1 2   N r b 2 N r c 2 N r c 2 N σ ] .
The derivative is negative for all rc if σ is negative, i.e., λbλ. Then, the minimum occurs at the right end rc = rbrp. The derivative becomes positive if the right-hand expression within the square brackets is negative. This may happen for larger ratios of λb/λ and rc close to rbrp, but for N greater than, say 5, the minimum lies very close to the right-hand limit, and this limit can be used quite safely in the present applications to determine the minimum of Rb0: R bmin = R b 0 | r c = r b r p .
For λb larger than λ, it may be better to distance the pipes slightly from the energy pile periphery to avoid too much heat flow in the ground region with its lower conductivity. The effect on the minimum is, however, negligible in the present applications. For the cases N = 1, 2, 3, 4, the effect is somewhat greater and worth considering.

4.2. Minimum Thermal Resistance for rc = rb − rp

The minimum value for the thermal resistance of the energy pile occurs exactly or with good accuracy when the pipes touch the periphery of the energy pile. The formula for the minimum is obtained from (31):
R bmin = R p N + 1 N 2   π   λ b [ ln ( r b N N   r p   ( r b r p ) N 1 ) + σ ln ( r b 2 N r b 2 N ( r b r p ) 2 N ) ] .
The minimum resistance depends on the ratio rp/rb in the following way:
R bmin = R p N + 1 2   π   λ b   N R ( r p ) , r p = r p r b , r c = r b r p , R ( r p ) = ln ( 1 N   r p ) + ( N 1 ) ln ( 1 1 r p ) + σ ln ( 1 1 ( 1 r p ) 2 N ) .
The ratio rp/rb is rather small. The following approximation (derived from suitable Taylor expansions of the logarithms) shows the behavior for small ratios:
R ( r p ) ln ( 1 N   r p ) + ( N 1 ) ( r p + 0.5 ( r p ) 2 ) + σ [ ln ( 1 2     N   r p ) + ( N 0.5 ) ( r p 2 N 5 12 ( r p ) 2 ) ] .
The two quadratic terms in rp′ may be neglected for say rp′ ≤ 0.1. The minimum resistance of the energy pile when the pipes touch the pile wall, rc = rbrp, becomes, with good accuracy:
R bmain R p N + 1 N 2   π   λ b { ( 1 + σ ) [ ln ( r b N   r p ) +   ( N 1 ) r p r b ] σ [ ln ( 2 ) 0.5 ]     r p r b } .

5. Formulas and Graphs for the Temperature Field

The temperature field in Cartesian and polar coordinates is given in Section 3.1, Equations (11)–(14). Here, the temperature field is plotted and discussed for a reference case.

5.1. Reference Case

The following data are chosen as the reference case:
N = 8 , q 0 = 10   W / m , λ b = 1.5   W / ( m K ) , λ = 3   W / ( m K ) , 2 r b = 0.6   m , r c = r b r p , 2 r p = 0.032   m , β = 0.75 , T bav = 0   ° C .
The pipe resistance Rp and the conductivity parameter σ become:
R p = β 2   π   λ b = 0.0796   m K / W ,   σ = λ b λ λ b + λ = 1 3 .
The average temperature Tbav at the pile wall is put to zero in all formulas and graphs below for simplicity. The resistance of the energy pile and the fluid temperature become (31):
R b 0 = 0.024   m K / W ,   K b 0 = 1 / R b 0 = 41.7 , T f 0 = R b 0 N   q p = 1.916 .

5.2. Plots of the Temperature Field

The overall temperature field T (x, y) = Tpol (r, φ) for the reference case (39) is presented in Figure 3, left. Six isotherms, T = 1, 0, −1, −2, −3, −4, are shown. The black circle r = rb shows the periphery of the energy pipe. The isotherm T = 0 winds around the pile periphery along which the average temperature is zero. The pipes with the fluid temperature Tf0 = 1.916 in heat carrier fluid lie inside the N = 8 isothermal circles for T = 1. The temperatures at a few particular points become:
T ( 0 , 0 ) = 0.465 , T pol ( r b , 0 ) = 0.733 , T pol ( r b , π N ) = 0.352 , T ( r c + r p , 0 ) = 0.733 , T ( r c , r p ) = 1.16 , T ( r c r p , 0 ) = 1.41 .
The isotherm T = −1 has a wavy form influenced by the position of the pipes. The deviation from a circle decreases outwardly. The isotherm T = −3 is very close to a perfect circle.
Figure 3, right, shows a surface plot of the temperature field for the reference case. The pipes are clearly seen as temperature spikes. The very flat region around the center z = 0 with values close to T(0,0) = 0.465 is to be noted.
Figure 4 shows radial profiles Tpol (r, φ) from (13) and (14) in the interval 0.2 < r < 0.5 for φ = 0, π/(5N), (3π)/(5N), and π/N. The red top curve for φ = 0 tends to infinity at the position r = rc of the line source. The bottom curve shows the radial profile for the radial line φ = π/N, which lies in the middle between pipe n = N and pipe n = 1 in Figure 2. The black dotted line shows the radial average temperature Tpolav (r), which is defined below in (46). The four curves converge to the radial average temperature for r greater than 0.4. The difference between the top and bottom curve is 0.091, 0.015, 0.004, and 0.00006 at r = 0.4, 0.5, 0.6, and 1, respectively.
The four curves lie very close to the value T (0, 0) = 0.465 at the center for 0 < r < 0.2. The deviation from temperature at r = 0 is less than 0.0002 for r < 0.1 and less than 0.06 for r < 0.2. The temperature field is very flat in the center region of Figure 3 (left) and 3 (right). This means that changes in the thermal conductivity or introduction of a return pipe with the same fluid temperature in the central region r < 0.2 will have a negligible, or perhaps marginal, effect on the value of the thermal resistance.

6. Total Radial Heat Flux: Thermal Resistance to Any Radius r > rb

The Laplace equation in polar coordinates for a thermal conductivity that depends on r reads:
1 r r ( r λ ( r ) T pol r     ) + 1 r 2 φ ( λ ( r ) T pol φ     ) = 0 , λ ( r ) = { λ b r < r b λ r > r b .
The equation is valid everywhere except at the N heat sources on the circle r = rc. Let Tpolav (r) denote the average temperature obtained from integration over −π < φ < π. Integration of (43) in φ gives:
T polav ( r ) = 1 2     π π π T pol ( r , φ )     d φ 1 r d d r [ r λ ( r ) d T polav d r       ] + 1 r 2 [ λ ( r ) 2     π T pol φ     ] φ = π φ = π = 0 .
The heat flux in the φ-direction is the same at φ = π and φ = −π. The right-hand term above vanishes and the average temperature Tpolav (r) satisfies the one-dimensional radial heat conduction equation except at r = rc, where the total radial heat flux qtotrad (r) increases from zero to N·q0:
d d r [ r λ ( r ) d T polav d r       ] = 0 , r r c ; q totrad ( r ) = 2     π   r λ ( r ) d T polav d r       = { 0 r < r c N   q 0 r > r c .
The temperature Tpolav (r) becomes constant for r < rc. It is zero at r = rb (for Tbav = 0) and it varies logarithmically as ln (rb/r) for r > rc, inversely proportional to the thermal conductivities λb and λ:
T polav ( r b ) = 0         T polav ( r ) =   N   q 0 2     π   λ ln ( r b r ) , r r b ; T polav ( r ) =   N   q 0 2     π   λ b ln ( r b r ) , r c r r b , T polav ( r ) = N   q 0 2     π   λ b ln ( r b r c ) , 0 r r c .
This average temperature is shown by the dotted line in Figure 4. It should be noted that the derivative of the average temperature is discontinuous at r = rb following λ (r), (43) right, while the heat flux is continuous. It should also be noted that Tpolav (0) from (46), right, gives the same value T (0) = 0.465 as in (42).
The thermal resistance Rb refers to the resistance between the fluid in the pipes and the average temperature at the periphery of the energy pile, (5). The corresponding thermal resistance between the fluid in the pipes and the average temperature at any radius r0 > rb is obtained in the following way:
r 0 > r b : T f T polav ( r 0 ) = T f T polav ( r b ) + T polav ( r b ) T polav ( r 0 ) = N   q 0 R b N   q 0 2     π   λ ln ( r b r 0 ) .
The resistance from fluid to r = r0 becomes equal to the sum of Rb and the thermal resistance of the annular region rb < r < r0 with its thermal conductivity λ:
R f   to   r 0 = R b + 1 2     π   λ ln ( r 0 r b ) , r 0 > r b .
This simple thermal network is illustrated in Figure 5.

7. Deviation from the Exact Boundary Condition at the Pipes

The boundary condition at the pipes is discussed in Section 3.2. Equation (20) gives the exact boundary condition at the periphery of the pipes:
T BC ( ψ ) = T c ( r c + r p e i ψ ) + q p N ( ψ ) h p = T f , π < ψ π .
The temperature at the outside of the pipe as a function of z is given by (8) and (9), top line. The radial heat flux qρN (ψ) is defined in (19), right. It is fairly straight-forward to determine this radial derivative with respect to ρ from the formula for the temperature:
q ρ N ( ψ ) = N   q 0 2   π   Re [ z ( ψ ) N 1 e i ψ z ( ψ ) N r c N σ z ( ψ ) N 1 e i ψ r c N r b 2 N z ( ψ ) N   r c N ] , z ( ψ ) = r c + r p e i ψ .
This gives an explicit formula for TBC (ψ).
Figure 6 shows the function (49) for the exact boundary condition around the periphery of the pipes for the reference case (39). The chosen fluid temperature Tf0 is the average or mean value of the boundary function over ψ as derived above in Section 3.2, Equations (28) and (30). The difference between the two curves represents the deviation from the exact boundary condition. The largest difference 0.092 between TBC (ψ) and the chosen Tf0 occurs at ψ = ±π. The relative error lies below 5%. The temperature at the outside of the pipe and the radial heat flux qpN, (50), divided by the heat transfer coefficient hp are also shown. These curves vary more than TBC (ψ), but they counteract each other so that the deviation in the boundary condition is much smaller.

8. Comparison with Exact Results from Multipoles Solutions of Order J = 8

Higher-order multipole solutions (J = 1, 2, …) give the thermal resistance with increasing accuracy. The general solution for any order of multipoles J will be reported in a sequel paper. Multipoles up to J = 8 are used in the comparisons of accuracy for the thermal resistance of the energy pile.

8.1. Error for Rb (J) as a Function of J

The error for Rb (J) compared to Rb (8) is:
Δ R b ( J ) = R b ( J ) R b ( 8 ) R b ( 8 ) = T f ( J ) T f ( 8 ) T f ( 8 ) , J = 0 ,   1 ,     7 .
Figure 7 shows the absolute relative error |∆Rb (J)| for J = 0, 1, … 7 for the reference case.
The relative error for J = 0 is 0.0069 or 0.69%. There is a dramatic increase in accuracy with increasing J. This is typical for the multipole method. The multipole method gives the solution and in particular the thermal resistance Rb (8) with 8 digits of accuracy. The error for thermal resistance is always much smaller than the relative deviations for the boundary conditions at the pipes; see Figure 6.

8.2. Relative Error ∆Rb0 (N, rb, rc, rp, σ, β) for a Suitable Set of Parameter Values

The input parameters (1) determine the temperature field and the thermal resistance of the energy pile. The temperature field (11) and (12) is proportional to q0/λb. The non-dimensional relative error ∆Rb0 = ∆Rb (0) with J = 0 in (51) will depend on the number of pipes N, two dimensionless parameters, and the three radii:
Δ R b 0 = Δ R b 0 ( N , σ , β , r p , r c , r b ) ; r p sin ( π   /   N ) r c r b r p , β 0 , 1 < σ < 1 .
It may be noted that the relative error is unchanged when the three radii are scaled with the same factor k > 0:
Δ R b 0 ( N , σ , β , r p , r c , r b ) = Δ R b 0 ( N , σ , β , k r p , k r c , k r b ) .
To investigate the performance of the proposed formula, a parametric study is made, varying the number of pipes N, the pile and ground thermal conductivities, the dimensionless pipe thermal resistance, and the pile geometry. The parameters and levels for the full parametric study, representing 6 × 5 × 3 × 3 × 4 = 1080 cases, are summarized in (54):
N = 2 ,   4 ,   6 ,   8 ,   10 ,   12 ; 2   r b = 0.16 ,   0.3 ,   0.6 ,   1.2 ,   2.4 ; r c = r c , min ,   2   r b 3 ,   r b r p ; 2   r p = 0.032 ; λ b   /   λ = 0.5 ,   1 ,   2 ,     or     σ = 1 3 ,   0 ,   1 3 ; β = 0.25 , 0.5 ,   1 , 2 .
These first 3 parameters represent 90 different geometries; in all cases, the outer diameter of the U-pipe is 0.032 m. The pile diameters are between 0.16 and 2.4 m. The smallest diameter is smaller than a typical precast driven thermal pile and the largest diameter is closer to what might be found in a large CFA or bored thermal pile. Taken together, the considered diameters bracket almost all energy pile installations. The three rc radii correspond to close, moderate, and wide pipe spacings. The pipes touch each other and the wall of the energy pile for rc equal to rc,min (close spacing), and rc equal to rbrp (wide spacing), respectively. These cases are typical of CFA and hollow driven piles, respectively. The case rc equal to 2rb/3 (moderate spacing) is representative of a bored energy pile where heat exchanger pipes are mounted on the reinforcement ring. Again, the three spacings virtually bracket all possible spacings. Similarly, the three levels of sigma bracket the most likely ratios of the borehole and ground thermal conductivities found in practice.
Table 1 shows the relative errors ∆Rb0 (N, rb, rc, rp, σ, β) in percent. The grey cells in the table indicate relative errors of 10% or more. More than 50% of the errors are marked as zero, which means that the error is below 0.5%. For rb greater than or equal to 0.3 m, the relative error ∆Rb0 (N, rb, rc, rp, σ, β) is less than or equal to 5% for 643 of the 648 cases and less than 10% for all cases. For rb = 0.15 m, the relative error ∆Rb0 (N, rb, rc, rp, σ, β) is less than or equal to 5% for 156 of the 216 cases and less than 10% for all except a few cases. The relative error exceeds 10% only for 8 cases when N is greater than or equal to 8 and β is equal to 0.25 or 0.50 (i.e., for very low Rp and Rb values). These cases are very unlikely in practice and are possible only under the implausible conditions of extremely low thermal conductivity of energy piles and laminar flow in the heat exchanger pipes. Moreover, although the relative error in these cases is large, the absolute error is small and of the order of 0.002 (K·m)/W) or less. For rb equal to 0.08 m, the relative error ∆Rb0 (N, rb, rc, rp, σ, β) is greater than 10% for over one-third of the considered cases. With N equal to 12 and rp equal to 0.016 m, the rc,min becomes larger than 2rb/3, and hence the moderate spacing condition does not hold. The relative error for these cases with rb equal to 0.08 m is classified as not applicable (n/a) in Table 1. An interesting observation from Table 1 is that the relative error ∆Rb0 (N, rb, rc, rp, σ, β) is positive for β less than 1 and is negative for β larger than 1. This suggests that the zero-order formula overestimates the thermal resistance of energy piles for β less than 1 and underestimates the thermal resistance for β greater than 1. For all cases with β equal to 1, the relative error is fairly close to zero. This information is valuable when setting margins of safety for design.
One parameter that can affect the accuracy of the multipole Formula (39) is the pipe radius. To demonstrate the impact of different pipe radii on the relative errors of the parametric study (54), Table 2 and Table 3 present relative errors ∆Rb0 (N, rb, rc, rp, σ, β) for rb equal to 150 mm and rp equal to 12.5 mm, and rb equal to 300 mm and rp equal to 20 mm, respectively. Comparing the results of Table 1 and Table 2 indicates that for energy piles with rb equal to 0.150 m, the relative errors ∆Rb0 (N, rb, rc, rp, σ, β) become considerably lower for rp equal to 0.0125 m than for rp equal to 0.016 mm. For rp equal to 0.0125 m, the relative error exceeds 10% only when N is greater than or equal to 10, and β equal to 0.25. Comparison of the results of Table 1 and Table 3 shows that for energy piles with rb equal to 0.3 m, the relative errors ∆Rb0 (N, rb, rc, rp, σ, β) are marginally higher for rp equal to 0.02 m than for rp equal to 0.016 m. Nevertheless, the relative error for rp equal to 0.02 m exceeds 10% only for the single case of N equal to 12, σ equal to −0.33, and β equal to 0.25.

9. Further Discussion and Concluding Summary

The newly developed zeroth-order multipole formula provides a fast, flexible, and accurate method for calculating the thermal resistance of energy piles. The newly presented zero-order multipole Formula (31) for Rb0 may be used in virtually all engineering applications for energy piles with rb greater than or equal to 0.3 m and rp less than or equal to 0.02 m with pipes equally spaced on a circle. The relative error ∆Rb0 (N, rb, rc, rp, σ, β) may exceed 10% in a few odd cases, but the absolute error in such cases is very small. Furthermore, the zero-order multipole Formula (31) may also be used for all energy piles with rb between 0.15 and 0.3 m and rp less than or equal to 0.0125 m if N is less than or equal to 10, and β is greater than or equal to 0.25. For applications with rp equal to 16 mm and rb equal to 0.15 to 0.3 m, the zero-order multipole Formula (31) may be used if N is less than or equal to 8. The zero-order multipole formula may also be used for all energy pile applications if β is equal to 1 or N is equal to 2.
To the best of the authors’ knowledge, zeroth-order multipole Formula (31) is the first, and the only closed-form formula published to date, that can be used to calculate the thermal resistance of energy piles with any number of heat exchanger pipes. The proposed formula marks a substantial improvement on the current practice of sizing energy piles with three or more U-pipes based on a guessed value of thermal resistance, which can lead to improperly sized systems. The explicit formula will facilitate designers to step beyond the limits of using a conservative thermal resistance estimate, and will enable them to optimize the energy pile design by making informed decisions about design choices, such as the position and the number of heat exchanger pipes in the pile. The proper sizing of energy piles will, in turn, allow for the overall system size to be reduced and the system performance to be improved.
A possible direction for future research is the development of higher-order multipole formulas for calculating the local thermal resistance of energy piles. The higher-order formulas will further improve the accuracy of the thermal resistance calculation, especially for small-diameter energy piles (less than 300 mm). The derivation of such formulas is expected in the near future and they will be presented soon in a sequel paper.
To conclude, in this paper, a closed-form Formula (31) for the thermal resistance of energy piles with any number of heat exchanger pipes evenly spaced on a concentric circle inside the pile periphery was derived. The proposed formula is based on the zero-order multipole approximation (J = 0) and considers only line heat sources at the pipes. It was demonstrated that the newly derived formula can calculate the thermal resistance of large-diameter energy piles (rb ≥ 0.3 m) with an accuracy better than 10% essentially under all conditions. The relative errors are less than or equal to 5% for over 99% of the 648 cases analyzed in this study. It was also demonstrated that the zero-order formula can also be used for smaller-diameter energy piles (0.15 ≤ rb < 0.3 m) for most conditions, except when the dimensionless pipe thermal resistance β is very small and the number N of heat exchanger pipes in the energy pile tops 8. For these specific cases and for energy piles with particularly smaller diameters (rb < 0.15 m), the use of higher-order multipoles is recommended—closed-form formulas for which will be presented in our forthcoming papers. This paper also presents a handy Formula (35) for the minimum thermal resistance of energy piles when the pipes touch the periphery of the pile.

Author Contributions

Both authors contributed equally to the design, concept, and preparation of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

hpHeat transfer coefficient between fluid and outer side of pipe wall, W/(m2·K), see Equation (19)
JNumber of multipoles, J = 0, 1, 2 …
NNumber of pipes in the energy pile; N = 2 for single U-tube
q0Heat rejection rate from each pipe to the ground, W/m
RbLocal thermal resistance between fluid in the U-tube to the energy pile wall, (m·K)/W
Rb0Local thermal resistance for the zero-order, (m·K)/W, see Equation (31)
RbminMinimum local thermal resistance, (m·K)/W, see Equation (38)
RpTotal fluid-to-pipe resistance for a single pipe i.e., one leg of the U-tube, (m·K)/W
rbRadius of the energy pile, m
rcRadius of the concentric circle, m, see Figure 2
rpOuter radius of heat exchanger pipe, m
TTemperature, °C
TbEnergy pile wall temperature, °C
TbavAverage temperature at the energy pile wall, °C
TfMean fluid temperature inside the U-tube, °C
βDimensionless thermal resistance of one U-tube leg, see Equation (2)
Rb0Relative error in local thermal resistance for zero-order, %
λThermal conductivity of the ground, W/(m·K)
λbThermal conductivity of the energy pile, W/(m·K)
σThermal conductivity ratio, dimensionless, see Equation (2)

Appendix A. Detailed Mathematical Derivations for the Temperature Field

Appendix A.1. Complex-Valued Line Heat Source

The complex-valued line heat source function W0 (z, z0) is the basic solution and a starting point for the multipole method:
W 0 ( z , z 0 ) = { ln ( r b z z 0 ) + σ ln ( r b 2 r b 2 z ¯ z 0 )   ,     | z | r b , ( 1 + σ ) ln ( r b z z 0 ) + σ 1 + σ 1 σ ln ( r b z ) ,   | z | r b .
The real part of the function W0 (z, z0) gives the temperature field for a line heat source at z = z0 with the strength q = 2πλb (W/m). It will satisfy the Laplace equation in pile and ground regions since it is the real part of an analytical function except for the right-hand top logarithm involving z ¯ . However, this latter part is the complex-conjugate of the analytic logarithm function defined by (A3), and its real part will satisfy the Laplace equation. This is a bit tricky and must be kept in mind in the analyses and derivations below.
The following complex notations are used:
z = x + i y ,   | z | = x 2 + y 2 ,   z ¯ = x i y ,   z 0 = x 0 + i y 0 ,   | z 0 | < r b r p .
The complex-valued logarithm of any complex function f(z) is defined by (7):
ln ( f ( z ) ) = ln ( | f ( z ) | ) + i arg ( f ( z ) ) , π < arg ( f ( z ) ) π
The upper expression of (A1) is valid in the circular energy pile region and the lower one in the ground outside the pile. The first top logarithm represents an ordinary line heat source, while the second term represents a line source at the “mirror” point z = r b 2 / z ¯ 0 . The expressions in the two regions are constructed so that the real part, and the ensuing radial heat flux at the pile and ground sides, become equal at the pile radius |z| = rb This means that the real solution satisfies the two internal boundary conditions at the pile wall. The average temperature from the real part of W0 (z, z0) (which will determine Tbav) taken around the pile periphery is zero.
The real part of W0 (z, z0) will be used in the temperature fields below. The complex conjugate of z times z0 may be replaced by z times the complex conjugate of z0:
Re [ ln ( r b 2 r b 2 z ¯ z 0 ) ] = ln ( r b 2 | r b 2 z ¯ z 0 | ) = ln ( r b 2 | r b 2 z z ¯ 0 | ) = Re [ ln ( r b 2 r b 2 z z ¯ 0 ) ] .
The expression (A1) outside the pile (lower line) may be written in the following alternative way:
| z | r b :     W 0 ( z , z 0 ) = 1 + σ 1 σ ln ( r b z ) + ( 1 + σ ) ln ( z z z 0 ) , 1 + σ 1 σ = λ b λ , 1 + σ = 2   λ b λ b + λ .
The first part presents a simple line source at the center z = 0 with the radial logarithmic temperature decrease ln(rb/r). The second part tends to zero as r = |z| increases and it determines the φ-dependence for r > rb.

Appendix A.2. Temperature Field from N Line Heat Sources on a Circle

The (steady-state, real-valued) temperature field Tc (z) for the pipes on the circle is to be determined for the prescribed heat flux q0 (W/m) from all pipes, and any prescribed average temperature Tbav around the pile periphery. From this, the required fluid temperature Tf is to be determined from the boundary condition at the pipes. The thermal resistance Rb of the energy pile is then determined in accordance with (5).
The temperature field consists of a linear combination of the real part of the line heat source solution taken at the center of each pipe:
T c ( z ) = T bav + q 0 2   π   λ b Re [ W c ( z , r c ) ] , z n = r c z e n .
W c ( z , r c ) = { n =   1 N [ ln ( r b z z n ) + σ ln ( r b 2 r b 2 z z ¯ n ) ]   ,     | z | r b , λ b λ [ N ln ( r b z ) + 2   λ λ b + λ n =   1 N ln ( z z z n ) ] ,   | z | r b ,   z ¯ n = r c z e n .
where Wc (z, rc) is the sum of the N complex-valued line sources (A1) for the energy pile and ground regions. The “conjugate” substitution (A4) is used. The prescribed average temperature is to be added since the average temperature at the pile periphery is zero for each term in the sum. The factor before the sum ensures that the heat flux from each pipe is q0. The total heat flux from all pipes to the surrounding ground is N·q0 (W/m).
Equation (A6) is the line-source solution without using multipoles. The above temperature field satisfies the five conditions listed in Section 3.1 (after Equation (12)).

Appendix A.3. New Formulas for the Temperature Field

The sums in (A7) for Wc (z, rc) of the N complex-valued line sources involve two sums of logarithms:
Σ 1 = n = 1 N ln ( r b z z n ) , Σ 2 = n = 1 N ln ( r b 2 r b 2 z z ¯ n ) , z ¯ n = r c z e n .
The first sum of logarithms may be rewritten as a logarithm of products:
Σ 1 = ln [ r b N ( z r c z e 1 ) ( z r c z e 2 ) ( z r e z φ N ) ] = ln ( r b N z N r c N ) .
The roots of the polynomial zNrcN = 0 are z = zn, n = 1, … N. The two polynomials in the denominator of order N have the same roots. Hence, they are identical. The other sum of logarithms may, in the same way, be rewritten as:
Σ 2 = ln [ r b 2 N ( r b 2 z r c z e 1 ) ( r b 2 z r c z e N ) ] = ln ( r b 2 N r b 2 N z N r c N ) .
where, again, the two polynomials in the right-hand denominator of the logarithms are identical since they have the same roots:
r b 2 N ( z r c ) N = 0         z r c z e n = r b 2 , n = 1 ,   2 ,     N .
It may be noted that the two left-hand expressions involving logarithms in (A9) and (A10) may differ by 2πi times an integer depending on the choice of arguments in the logarithms. However, this does not influence the real part used for the temperature field. This comment is applicable to the addition of logarithms in all other cases in this paper.
The complex-valued temperature field (A7) may now be written in the following concise form inserting the sums (A8)–(A10):
W c ( z , r c ) = { ln ( r b N z N r c N ) + σ ln ( r b 2 N r b 2 N z N r c N )   ,   | z | r b λ b λ [ N ln ( r b z ) + 2   λ λ b + λ ln ( z N z N r c N ) ] ,   | z | r b .
These formulas are the starting point in Section 3.1, Equation (9). The above temperature fields depend on zN. Let the point z be “rotated” m times by ze.
( z e m z ) N = ( e 2 π i N m ) N z N = ( e 2 π i N N ) m z N = z N Re [ W c ( z e m z , r c ) ] = Re [ W c ( z , r c ) ] .
The real part of the temperature field repeats itself for any rotation (ze)m. This means that the temperature field is the same in every wedge with the opening angle Δφ = 2π/N as indicated in Figure 2, left. The temperature field is also symmetric in y. It is sufficient to consider the temperature field in the wedge 0 ≤ φ ≤ π/N, r ≥ 0.
Figure 3 (left) and (right), in Section 5.2 show the temperature field from (A6) and (A12) for the reference case (39). The black circle r = rb = 0.3 shows the periphery of the energy pile. The positions of the N pipes and the rotational symmetry (N = 8) are clearly seen.

References

  1. Rees, S. (Ed.) Advances in Ground-Source Heat Pump Systems; Woodhead Publishing: Cambridge, UK, 2016. [Google Scholar]
  2. Garber, D.; Choudhary, R.; Soga, K. Risk Based Lifetime Costs Assessment of a Ground Source Heat Pump (GSHP) System Design: Methodology and Case Study. Build. Environ. 2013, 60, 66–80. [Google Scholar] [CrossRef]
  3. Rivoire, M.; Casasso, A.; Piga, B.; Sethi, R. Assessment of Energetic, Economic and Environmental Performance of Ground-Coupled Heat Pumps. Energies 2018, 11, 1941. [Google Scholar] [CrossRef] [Green Version]
  4. Vieira, A.; Alberdi-Pagola, M.; Christodoulides, P.; Javed, S.; Loveridge, F.; Nguyen, F.; Cecinato, F.; Maranha, J.; Florides, G.; Prodan, I.; et al. Characterisation of Ground Thermal and Thermo-Mechanical Behaviour for Shallow Geothermal Energy Applications. Energies 2017, 10, 2044. [Google Scholar] [CrossRef] [Green Version]
  5. Loveridge, F.; Powrie, W. Temperature response functions (G-functions) for single pile heat exchangers. Energy 2013, 57, 554–564. [Google Scholar] [CrossRef]
  6. Loveridge, F.; Powrie, W.; Nicholson, D. Comparison of two different models for pile thermal response test interpretation. Acta Geotech. 2014, 9, 367–384. [Google Scholar] [CrossRef] [Green Version]
  7. Park, H.; Lee, S.R.; Yoon, S.; Choi, J.C. Evaluation of thermal response and performance of PHC energy pile: Field experiments and numerical simulation. Appl. Energy 2013, 103, 12–24. [Google Scholar] [CrossRef]
  8. Lamarche, L.; Kajl, S.; Beauchamp, B. A review of methods to evaluate borehole thermal resistances in geothermal heat-pump systems. Geothermics 2010, 39, 187–200. [Google Scholar] [CrossRef]
  9. Javed, S.; Spitler, J.D. Calculation of borehole thermal resistance. In Advances in Ground-Source Heat Pump Systems, 1st ed.; Rees, S.J., Ed.; Woodhead Publishing: Cambridge, UK, 2016; pp. 63–95. [Google Scholar]
  10. Javed, S.; Spitler, J.D. Accuracy of borehole thermal resistance calculation methods for grouted single U-tube ground heat exchangers. Appl. Energy 2017, 187, 790–806. [Google Scholar] [CrossRef]
  11. Loveridge, F.; Powrie, W. 2D thermal resistance of pile heat exchangers. Geothermics 2014, 50, 122–135. [Google Scholar] [CrossRef]
  12. Paul, N.D. The Effect of Grout Thermal Conductivity on Vertical Geothermal Heat Exchanger Design and Performance. Master’s Thesis, South Dakota University, Brookings, SD, USA, July 1996. [Google Scholar]
  13. Sharqawy, M.H.; Mokheimer, E.M.; Badr, H.M. Effective pipe-to-borehole thermal resistance for vertical ground heat exchangers. Geothermics 2009, 38, 271–277. [Google Scholar] [CrossRef]
  14. Bauer, D.; Heidemann, W.; Müller-Steinhagen, H.; Diersch, H.J. Thermal resistance and capacity models for borehole heat exchangers. Int. J. Energy Res. 2011, 35, 312–320. [Google Scholar] [CrossRef]
  15. Diao, N.R.; Zeng, H.Y.; Fang, Z.H. Improvement in modeling of heat transfer in vertical ground heat exchangers. HVAC R Res. 2004, 10, 459–470. [Google Scholar] [CrossRef]
  16. Bennet, J.; Claesson, J.; Hellström, G. Multipole Method to Compute the Conductive Heat Flows to and between Pipes in a Composite Cylinder; Notes on Heat Transfer 3; University of Lund: Lund, Sweden, 1987. [Google Scholar]
  17. Earth Energy Designer (EED), Version 4.2; BLOCON: Lund, Sweden, 2019.
  18. Ground Loop Heat Exchanger Professional (GLHEPro), Version 5.0; International Ground Source Heat Pump Association (IGSHPA): Stillwater, OK, USA, 2016.
  19. Claesson, J.; Hellström, G. Multipole method to calculate borehole thermal resistances in a borehole heat exchanger. HVAC R Res. 2011, 17, 895–911. [Google Scholar]
  20. Claesson, J.; Javed, S. Explicit multipole formulas for calculating thermal resistance of single U-tube ground heat exchangers. Energies 2018, 11, 214. [Google Scholar] [CrossRef] [Green Version]
  21. Claesson, J.; Javed, S. Explicit multipole formulas and thermal network models for calculating thermal resistances of double U-pipe borehole heat exchangers. Sci. Technol. Built Environ. 2019, 25, 980–992. [Google Scholar] [CrossRef] [Green Version]
  22. Nehari, Z. Introduction to Complex Analysis; Allen and Bacon Inc.: Boston, MA, USA, 1961. [Google Scholar]
Figure 1. Common types of energy piles.
Figure 1. Common types of energy piles.
Energies 13 05445 g001
Figure 2. Energy pile with N (=8) heat exchanger pipes equally spaced on the circle r = rc.
Figure 2. Energy pile with N (=8) heat exchanger pipes equally spaced on the circle r = rc.
Energies 13 05445 g002
Figure 3. Temperature field (left) and surface plot of the temperature field (right) for the reference case (39).
Figure 3. Temperature field (left) and surface plot of the temperature field (right) for the reference case (39).
Energies 13 05445 g003
Figure 4. Radial temperature profiles for φ = 0, π/(5N), 3π/(5N), and π/N. The dotted curve shows the average temperature Tpolav (r), (46).
Figure 4. Radial temperature profiles for φ = 0, π/(5N), 3π/(5N), and π/N. The dotted curve shows the average temperature Tpolav (r), (46).
Energies 13 05445 g004
Figure 5. Thermal resistance network between the fluid in the pipes and the average temperature at any radius r0 > rb.
Figure 5. Thermal resistance network between the fluid in the pipes and the average temperature at any radius r0 > rb.
Energies 13 05445 g005
Figure 6. Boundary-condition function TBC (ψ) (49), and its average Tf0. The dotted curves show the boundary temperature and the radial heat flux divided by the heat transfer coefficient (50).
Figure 6. Boundary-condition function TBC (ψ) (49), and its average Tf0. The dotted curves show the boundary temperature and the radial heat flux divided by the heat transfer coefficient (50).
Energies 13 05445 g006
Figure 7. Absolute relative error ∆Rb for J = 0, 1, 2, …7 for the reference case (39).
Figure 7. Absolute relative error ∆Rb for J = 0, 1, 2, …7 for the reference case (39).
Energies 13 05445 g007
Table 1. Relative error ∆Rb0 (N, rb, rc, rp, σ, β) in percent for N = 2 to 12, and rc = rc,min (close spacing), 2rb/3 (moderate spacing), and rbrp (wide spacing) a. The pipe radius rp is 0.016 mm.
Table 1. Relative error ∆Rb0 (N, rb, rc, rp, σ, β) in percent for N = 2 to 12, and rc = rc,min (close spacing), 2rb/3 (moderate spacing), and rbrp (wide spacing) a. The pipe radius rp is 0.016 mm.
rb (mm)σβCloseModerateWide
246810122468101224681012
80−0.330.255131721316528172533n/a31227466585
0.503791216281481216n/a1511182632
1.00000−1−1−20000−1n/a0000−1−2
2.00−3−8−10−13−16−210−2−5−9−15n/a−1−2−5−9−14−21
0.000.255131721316417162533n/a1822416182
0.503791216280381216n/a039172532
1.00000−1−1−20000−1n/a0000−1−2
2.00−3−8−10−13−16−210−2−5−9−15n/a0−2−5−9−14−21
0.330.255131721316307162533n/a0518375879
0.503791216270381216n/a028162431
1.00000−1−1−20000−1n/a0000−1−2
2.00−3−8−10−13−16−210−2−5−9−15n/a0−1−4−8−14−21
150−0.330.253899910024689249142129
0.502555560123451246912
1.00000000000000000000
2.00−2−6−7−7−8−80−1−2−3−5−70−1−2−3−4−5
0.000.25389991002468901491623
0.502555560123450124710
1.00000000000000000000
2.00−2−6−7−7−8−80−1−2−3−5−700−1−2−3−5
0.330.25389991002468900261219
0.50255556012345001359
1.00000000000000000000
2.00−2−6−7−7−8−80−1−2−3−5−700−1−1−3−4
300−0.330.25356555001122123579
0.50133333000111111234
1.00000000000000000000
2.00−2−5−5−5−5−50000−1−10−1−1−1−1−2
0.000.25356555001122001235
0.50133333000111000112
1.00000000000000000000
2.00−2−5−5−5−5−50000−1−10000−1−1
0.330.25356555001122000012
0.50133333000111000001
1.00000000000000000000
2.00−2−5−5−5−5−50000−1−100000−1
600−0.330.25244443000000112233
0.50133222000000011111
1.00000000000000000000
2.00−2−4−4−4−4−4000000000−1−1−1
0.000.25244443000000000011
0.50133222000000000000
1.00000000000000000000
2.00−2−4−4−4−4−4000000000000
0.330.25244443000000000000
0.50133222000000000000
1.00000000000000000000
2.00−2−4−4−4−4−4000000000000
1200−0.330.25233333000000111122
0.50122222000000000111
1.00000000000000000000
2.00−1−3−3−3−3−3000000000000
0.000.25233333000000000000
0.50122222000000000000
1.00000000000000000000
2.00−1−3−3−3−3−3000000000000
0.330.25233333000000000000
0.50122222000000000000
1.00000000000000000000
2.00−1−3−3−3−3−3000000000000
a Relative error larger than 10% shown with a grey background.
Table 2. Relative error ∆Rb0 (N, rb, rc, rp, σ, β) in percent for N = 2 to 12, and rc = rc,min (close spacing), 2rb/3 (moderate spacing), and rbrp (wide spacing) a. The pipe radius rp is 0.0125 m.
Table 2. Relative error ∆Rb0 (N, rb, rc, rp, σ, β) in percent for N = 2 to 12, and rc = rc,min (close spacing), 2rb/3 (moderate spacing), and rbrp (wide spacing) a. The pipe radius rp is 0.0125 m.
rb (mm)σβCloseModerateWide
246810122468101224681012
150−0.330.2537777701235613691419
0.50244444001233112468
1.00000000000000000000
2.00−2−5−6−6−7−700−1−1−2−30−1−1−2−3−3
0.000.253777770123560125913
0.50244444001233001246
1.00000000000000000000
2.00−2−5−6−6−7−700−1−1−2−300−1−1−2−3
0.330.253777770123560012510
0.50244444001223000134
1.00000000000000000000
2.00−2−5−6−6−7−700−1−1−2−3000−1−1−2
a Relative error larger than 10% shown with a grey background.
Table 3. Relative error ∆Rb0 (N, rb, rc, rp, σ, β) in percent for N = 2 to 12, and rc = rc,min (close spacing), 2rb/3 (moderate spacing), and rbrp (wide spacing) a. The pipe radius rp is 0.02 m.
Table 3. Relative error ∆Rb0 (N, rb, rc, rp, σ, β) in percent for N = 2 to 12, and rc = rc,min (close spacing), 2rb/3 (moderate spacing), and rbrp (wide spacing) a. The pipe radius rp is 0.02 m.
rb (mm)σβCloseModerateWide
246810122468101224681012
300−0.330.253666660112341347913
0.50144444001122112345
1.00000000000000000000
2.00−2−5−5−6−6−6000−1−1−20−1−1−1−2−2
0.000.25366666011234001358
0.50144444001122001124
1.00000000000000000000
2.00−2−5−5−6−6−6000−1−1−2000−1−1−2
0.330.25366666011234000125
0.50144444001122000012
1.00000000000000000000
2.00−2−5−5−6−6−6000−1−1−20000−1−1
a Relative error larger than 10% shown with a grey background.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Claesson, J.; Javed, S. Explicit Multipole Formula for the Local Thermal Resistance in an Energy Pile—The Line-Source Approximation. Energies 2020, 13, 5445. https://doi.org/10.3390/en13205445

AMA Style

Claesson J, Javed S. Explicit Multipole Formula for the Local Thermal Resistance in an Energy Pile—The Line-Source Approximation. Energies. 2020; 13(20):5445. https://doi.org/10.3390/en13205445

Chicago/Turabian Style

Claesson, Johan, and Saqib Javed. 2020. "Explicit Multipole Formula for the Local Thermal Resistance in an Energy Pile—The Line-Source Approximation" Energies 13, no. 20: 5445. https://doi.org/10.3390/en13205445

APA Style

Claesson, J., & Javed, S. (2020). Explicit Multipole Formula for the Local Thermal Resistance in an Energy Pile—The Line-Source Approximation. Energies, 13(20), 5445. https://doi.org/10.3390/en13205445

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop