Thermo-Hydraulic Performance of U-Tube Borehole Heat Exchanger with Different Cross-Sections

: For reducing the initial GSHP investment, the heat transfer efﬁciency of the borehole heat exchange (BHE) system can be enhanced to reduce the number or depth of drilling. This paper proposes a novel and simple BHE design by changing the cross-sectional shape of the U-tube to increase the heat transfer efﬁciency of BHEs. Speciﬁcally, in this study, we (1) veriﬁed the reliability of the three-dimensional numerical model based on the thermal response test (TRT) and (2) compared the inlet and outlet temperatures of the different U-tubes at 48 h under the premise of constant leg distance and ﬂuid area. Referent to the circular tube, the increases in the heat exchange efﬁciencies of the curved oval tube, ﬂat oval tube, semicircle tube, and sector tube were 13.0%, 19.1%, 9.4%, and 14.8%, respectively. (3) The heat ﬂux heterogeneity of the tubes on the inlet and outlet sides of the BHE, in decreasing order, is ﬂat oval, semicircle, curved oval, sector, and circle shapes. (4) The temperature heterogeneity of the borehole wall in the BHE in decreasing order is circle, sector, curved oval, ﬂat oval, and semicircle shapes. (5) Under the premise of maximum leg distance, referent to the heat resistance of the tube with a circle shape at 48 h, the heat exchange efﬁciency of the curved oval, ﬂat oval, semicircle, and sector tubes increased 12.6%, 17.7%, 10.3%, and 7.8%, respectively. (6) We found that the adjustments of the leg distance and the tube shape affect the heat resistance by about 25% and 12%, respectively. (7) The ﬂat-oval-shaped tube at the maximum leg distance was found to be the best tube design for BHEs. heat resistance


Introduction
With increasing awareness of environmental protection and the continued development of science and technology, the world energy structure is developing toward multipolarization, from traditional fossil fuels to cleaner wind energy, water energy, tidal energy, geothermal energy, etc. [1][2][3][4][5][6]. As a widely distributed resource, the commercial applications of geothermal energy are rapidly expanding [7]. If the soil temperature is high (greater than 70 • C), geothermal energy can be converted into other forms of easily transportable energy, such as electricity [8]; otherwise, it can be directly used or stored, as achieved using ground source heat pump system (GSHP) [9,10]. The GSHP is a technology that uses the relatively constant temperature of underground soil to convert the low-grade energy in soil into high-grade energy by consuming a small amount of high-grade energy (e.g., electricity), which has environmental advantages and economic benefits. When GSHPs are used for building heating, the radiant heating method provides much higher comfort than the traditional convection heat transfer method of the air conditioning system [11].
The underground heat exchanging aspect of GSHPs is the borehole heat exchanger (BHE) system [12,13], and includes soil, backfill soil, heat exchange tube, and a heat and fast calculation [33]. Hence, we used COMSOL (The COMSOL manufacturers and software companies, COMSOL Inc., Stockholm, Sweden; Suppliers, Kangmodur Software Technology (Shanghai) Co., Ltd., Shanghai, China) to calculate the heat transfer characteristic of BHEs.
Ground initial temperature distribution, soil thermal characteristic and insulation of borehole heat exchanger all played an important role in heating exchange efficiency of BHEs [32,[34][35][36], in which, ground heat conductivity and borehole heat resistance are the two most important parameters [37][38][39]. Ground heat conductivity reflects the heat physical properties of the soil and is related to soil density, water content, groundwater, and other factors [11]  ; the borehole heat resistance reflects the heat transfer effect between the U-pipe and the borehole wall, and is related to factors such as backfill material, U-pipe material, U-pipe shapes, and distribution position. To improve the heat exchange efficiency, scholars have proposed many novel design methods to reduce the heat resistance, such as enhanced heat transfer pipes [14], enhanced heat transfer fluid [40], spiral cloth pipe [41], eccentric U-pipe [42], and others [43,44].
To lower the initial investment of GSHPs, Serageldin et al. [13] proposed a new type of flat oval pipe. In comparison with the traditional circular buried pipe, this new type of oval buried pipe can reduce the internal heat resistance of the borehole by 0.125 m·K/W, and can increase the heat transfer efficiency of the borehole by 18.47%, which is particularly suitable for use in high-density areas such as cities. Jahanbin et al. [12] proposed a curved oval buried pipe and explored the influence of the difference in oval geometry on the heat resistance of the borehole. The results showed that the curved oval cross-section can significantly reduce the heat resistance of the borehole by about 17% compared with the traditional circular cross-section pipes, and this reduction increases with the increase in the oval dimensionless factor. Abuhamdeh et al. [45] studied the influence of semicircular and circular cross-section shapes on the heat transfer efficiency of spiral heat exchangers through numerical simulation. The results showed that spiral heat exchangers with a semicircular cross-section had a higher heat transfer efficiency than those with a circular cross-section. The heat transfer efficiency can be more economically and reliably altered by changing the shape of the cross-section than by adjusting the structure (e.g., spiral shape, heat conducting fin) and material (e.g., enhanced heat transfer backfill material, enhanced heat transfer pipe). However, studies on the shape of the buried pipe cross-section are lacking and the determination of its influence on the heat transfer efficiency of BHEs is necessary.
Based on the above analysis, we employed COMSOL multi-physics to study the heat transfer characteristics of BHEs with different cross-sectional pipes (circle, curved oval, flat oval, semicircle, and sector) through three-dimensional numerical analysis. First, based on the thermal response test (TRT), the rationality of the material parameters, boundary conditions, and geometric dimensions of the analysis model were verified; subsequently, under the premise of the same cross-sectional area and leg distance between pipes, the temperature at inlet and outlet, heat resistance of BHEs, heat flux distribution at inlet and outlet, and temperature distribution at borehole wall were studied. Finally, the heat resistance of BHEs with different cross-sectional pipes were studied with consideration of the maximum leg distance between pipes. The findings in this paper are expected to support the design and improvement of GSHPs.

Physical Model
As shown in Figure 1, the TRT was conducted in Huainan city, Anhui province, China, referring to the ground source heat pump system engineering technical specification of [46]. During the test, the TRT equipment and U-tube form a circulatory system. Then, the rated heating power was set to conduct continuous heat dissipation test for 48 h, and the heating power, water flow, and water temperature were recorded in real-time.
China, referring to the ground source heat pump system en tion of [46]. During the test, the TRT equipment and U-tub Then, the rated heating power was set to conduct continuou h, and the heating power, water flow, and water temperatur

Stratigraphic Descriptions
The landform of the test site is an alluvial landform and t 20 m of the proposed site is the quaternary loose layer, and th and mudstone. The strength of bedrock is generally medium which has higher strength. The specific soil parameters are s

Thermal Response Test (TRT)
As showing in Figure 1, the borehole depth of the TR borehole diameter db was 153 mm; the U-shaped tube was tube with a circular cross-section. The internal diameter dint the tube were 26 and 32 mm, respectively, and the distance b 75 mm (controlled by snap with fixed distance). The backfill m a little sand.

Stratigraphic Descriptions
The landform of the test site is an alluvial landform and the landform is single. Within 20 m of the proposed site is the quaternary loose layer, and the lower bedrock is sandstone and mudstone. The strength of bedrock is generally medium, and the locale is mudstone, which has higher strength. The specific soil parameters are shown in Table 1. As showing in Figure 1, the borehole depth of the TRT was set to 100 m, and the borehole diameter d b was 153 mm; the U-shaped tube was a high-density polyethylene tube with a circular cross-section. The internal diameter d int and external diameter d ext of the tube were 26 and 32 mm, respectively, and the distance between the tube legs d leg was 75 mm (controlled by snap with fixed distance). The backfill material was original soil and a little sand.
To obtain more accurate monitoring data, the flow meter included a high-precision turbine flow sensor of grade 0.5; for temperature measurement, we used class APT100 platinum resistance sensor, and the error was less than 0.15 • C. The measured ground temperature curve is shown in Figure 2, with an average formation temperature of 18.07 • C. The power meter was a WBP114S91 type (Weibo Electronics Co., Ltd., Mianyang, China), and the measurement error was less than 1%. The working voltage of the heater was 220 V and the rated heating power in this test was 3.5 kW (maximum power 8 kW).

Numerical Model
In the numerical analysis model, the calculation model used module, and as the research object, we selected transient heat tra ensure the reliability of the analysis results, the numerical model

•
The thermo-physical properties of the material are constant the change in soil temperature; • The effect of underground water flow on BHEs is neglected; • The soil layer is homogeneous; and • The fluid is incompressible.

Geometry and Meshing
The numerical model was mainly divided into four areas, a cluding soil outside the borehole, backfill soil inside the boreho shape, and heat transfer fluid. A three-dimensional model was es the geometric TRT dimensions. As shown in Figure 3a, the heigh 102 m, and the calculated radius of the model was set to 2 m, wh cording to the heat influence radius in the pre-calculation at 48 borehole db was 153 mm, the inner diameter of the buried tube d diameter dext was 32 mm, and the leg distance of U-tube dleg w parameters of the model are shown in Table 2.

Numerical Model
In the numerical analysis model, the calculation model used a conjugate heat transfer module, and as the research object, we selected transient heat transfer analysis [12,45]. To ensure the reliability of the analysis results, the numerical model is simplified as follows:

•
The thermo-physical properties of the material are constant and do not change with the change in soil temperature; • The effect of underground water flow on BHEs is neglected; • The soil layer is homogeneous; and • The fluid is incompressible.

Geometry and Meshing
The numerical model was mainly divided into four areas, as shown in Figure 3, including soil outside the borehole, backfill soil inside the borehole, HDPE tube with U-shape, and heat transfer fluid. A three-dimensional model was established by referring to the geometric TRT dimensions. As shown in Figure 3a, the height of the model was set to 102 m, and the calculated radius of the model was set to 2 m, which was determined according to the heat influence radius in the pre-calculation at 48 h. The diameter of the borehole d b was 153 mm, the inner diameter of the buried tube d int was 26 mm, the outer diameter d ext was 32 mm, and the leg distance of U-tube d leg was 75 mm. The material parameters of the model are shown in Table 2. Based on the existing research on the cross-sectional shape of U-tubes in BHEs [12,13,45], five analytical models were established. As shown in Figure 4, the cross-sectional shapes were circular, curved oval, flat oval, semicircle, and sector. To better compare the Sustainability 2021, 13, 3255 6 of 20 influence of cross-section shape on the heat transfer efficiency of BHEs, only the crosssection shape of the U-tube was adjusted, and the fluid cross-sectional area in each model was the same, 531 mm 2 . The geometric characteristics of each model are described in Table 3.
cluding soil outside the borehole, backfill soil inside the borehole, HDPE tube with Ushape, and heat transfer fluid. A three-dimensional model was established by referring to the geometric TRT dimensions. As shown in Figure 3a, the height of the model was set to 102 m, and the calculated radius of the model was set to 2 m, which was determined according to the heat influence radius in the pre-calculation at 48 h. The diameter of the borehole db was 153 mm, the inner diameter of the buried tube dint was 26 mm, the outer diameter dext was 32 mm, and the leg distance of U-tube dleg was 75 mm. The material parameters of the model are shown in Table 2.   Based on the existing research on the cross-sectional shape of U-tubes in BHEs [12,13,45], five analytical models were established. As shown in Figure 4, the cross-sectional shapes were circular, curved oval, flat oval, semicircle, and sector. To better compare the influence of cross-section shape on the heat transfer efficiency of BHEs, only the cross-section shape of the U-tube was adjusted, and the fluid cross-sectional area in each model was the same, 531 mm 2 . The geometric characteristics of each model are described in Table 3.

Governing Equations
To determine the transient flow and heat transfer characteristics of the fluid and solid domains, the corresponding governing equations needed to be established, including mass equation, momentum equation, and energy conservation equation.
The mass equation or the continuity equation is: where ρ is density, kg/m 3 ; u is fluid velocity, m/s; and ∇ is the Laplace operator.
The momentum equation is:

Governing Equations
To determine the transient flow and heat transfer characteristics of the fluid and solid domains, the corresponding governing equations needed to be established, including mass equation, momentum equation, and energy conservation equation.
The mass equation or the continuity equation is: where ρ is density, kg/m 3 ; u is fluid velocity, m/s; and ∇ is the Laplace operator. The momentum equation is: where p is static pressure, N/m 2 ; µ is dynamic clay, pa·s; T is the absolute temperature, K; I is the pressure direction vector; and F is the volume force vector, N. The energy conservation equation of the solid domain is: where q = −k∇ 2 T and Q ted is thermoelastic damping, J. The energy conservation equation in the fluid domain is as follows: where Q ρ is the work done by pressure changes and Q vd is the viscous dissipation in the fluid, J.
To better simulate the heat transfer characteristics of fluid flow near the tube wall, k − ε turbulent flow simulation analysis was used: where k is turbulent energy, m 2 /s 2 ; ε is the turbulent energy dissipation ratio, m 2 /s 3 ; P k is turbulence pressure; σ k and σ s are the Prandtl number of k and ε, respectively; c ε1 and c ε2 are experimental values; and µ T is the dynamic eddy viscosity, pa·s.

Initial Conditions and Boundary Conditions
In the COMSOL multi-physics analysis model, the reasonable choice of boundary conditions has a greater influence on the calculation results. Combining results from various references [12,13] and the TRT results, the corresponding heat transfer boundary conditions were established. As shown in Figure 3b, the constant temperature boundary condition was set on the top of the model, and the temperature was 30.85 • C. The left and right sides and bottom of the model were set as adiabatic boundary conditions; the initial temperature of the soil along depth was determined according to the testing temperature, as shown in Figure 2. The initial temperature of the fluid was the inlet temperature T in , which is the sum of the fluid outflow temperature T out and the fluid additional temperature T a heated by the heater in a water tank.

Heating Model by Rated Power
To consider the effect of rated power heating on the inlet temperature in the numerical analysis model, the corresponding heating tank analysis model was established. As shown in Figure 5, the volume of the water tank was 0.045 m 3 , the rated heating power was 3.5 kW, and the volume flow rate was 1.11 m 3 /h. The initial temperature of the inlet was 18 • C. The water enters the water tank from one side of the water tank (temperature T out ), which is heated by a constant power heater (energy Q), and flows into the buried tube (temperature T in ) from the other side. where V is the volume of water in the tank, m 3 ; c f is the heat capacity of water under constant pressure, J/kg·K; m is mass flow rate, kg/h; and Q is the heating power, J.
stainability 2021, 13, x FOR PEER REVIEW Figure 5. Heating water tank model with rated power.
( ) where V is the volume of water in the tank, m 3 ; cf is the heat capacity stant pressure, J/kg·K; m is mass flow rate, kg/h; and Q is the heating At the beginning of the TRT, the temperature of the fluid in the of the soil, which leads to less heat exchange. Hence, the temperatur is negligible, which means Tout is a constant value. In addition, due water tank being fully mixed, we obtain Equation (8): After a few hours, the left side of Equation (7) becomes very sm the right side and can be ignored, and then Equation (7) is adjusted: Thus, Equation (10) can be solved and Tin can be obtained: At the beginning of the TRT, the temperature of the fluid in the BHEs is near to that of the soil, which leads to less heat exchange. Hence, the temperature change rate of T out is negligible, which means T out is a constant value. In addition, due to the water in the water tank being fully mixed, we obtain Equation (8):

Results and Discussion
After a few hours, the left side of Equation (7) becomes very small compared with the right side and can be ignored, and then Equation (7) is adjusted: Considering T out as a constant value, T in (0) = T out . Thus, Equation (10) can be solved and T in can be obtained:

Model Validation
To verify the rationality of the three-dimensional analysis model, the temperature curve of the inlet/outlet obtained by the numerical method was compared with the temperature curve measured in the TRT, as shown in Figure 6. Figure 6 shows that under the combination of rated power heating and cooling of BHEs, the fluid temperature at the tube inlet/outlet continues increasing and gradually reaches an equilibrium state over time, which means the heat of fluid absorbed from the heater and the heat released by the BHEs reaches a dynamic balance. According to the TRT temperature curve, the inlet and outlet temperatures at 48 h are 31.27 and 27.99 • C, respectively, and where the temperature difference is 3.28 • C; the inlet and outlet temperatures are 31.07 and 27.77 • C where temperature difference is 3.30 • C in the numerical analysis model.
Based on the error calculation equation, Where the temperature errors of the inlet, outlet, and temperature differences are 0.6%, 0.8%, and 0.6%, respectively, which meet the requirements of literature [13]. Therefore, we verified that the numerical analysis model is reasonable and can be used to analyze the heat transfer characteristics of BHEs. tainability 2021, 13, x FOR PEER REVIEW Figure 6. Numerical model verification. Figure 6 shows that under the combination of rated power hea BHEs, the fluid temperature at the tube inlet/outlet continues increa reaches an equilibrium state over time, which means the heat of fluid heater and the heat released by the BHEs reaches a dynamic balanc TRT temperature curve, the inlet and outlet temperatures at 48 h are respectively, and where the temperature difference is 3.28 °C; the inle atures are 31.07 and 27.77 °C where temperature difference is 3.30 analysis model.
Based on the error calculation equation, Where the temperature errors of the inlet, outlet, and tempera 0.6%, 0.8%, and 0.6%, respectively, which meet the requirements of li fore, we verified that the numerical analysis model is reasonable and lyze the heat transfer characteristics of BHEs.

Temperature of Inlet and Outlet
To compare the influence of cross-sectional tube shape on the hea of BHEs, the corresponding analysis model was established by referr Section 3. The temperature curves of the inlet and outlet of a tube over in Figure 7.

Temperature of Inlet and Outlet
To compare the influence of cross-sectional tube shape on the heat transfer efficiency of BHEs, the corresponding analysis model was established by referring to parameters in Section 3. The temperature curves of the inlet and outlet of a tube over time are illustrated in Figure 7.  Figure 6 shows that under the combination of rated power heating and co BHEs, the fluid temperature at the tube inlet/outlet continues increasing and gr reaches an equilibrium state over time, which means the heat of fluid absorbed f heater and the heat released by the BHEs reaches a dynamic balance. Accordin TRT temperature curve, the inlet and outlet temperatures at 48 h are 31.27 and 2 respectively, and where the temperature difference is 3.28 °C; the inlet and outlet atures are 31.07 and 27.77 °C where temperature difference is 3.30 °C in the nu analysis model.
Based on the error calculation equation, Where the temperature errors of the inlet, outlet, and temperature differen 0.6%, 0.8%, and 0.6%, respectively, which meet the requirements of literature [13] fore, we verified that the numerical analysis model is reasonable and can be used lyze the heat transfer characteristics of BHEs.

Temperature of Inlet and Outlet
To compare the influence of cross-sectional tube shape on the heat transfer ef of BHEs, the corresponding analysis model was established by referring to param Section 3. The temperature curves of the inlet and outlet of a tube over time are illu in Figure 7.   Figure 7 indicates that the temperature of the inlet and outlet of the buried tube has a non-linear increasing trend, with a larger increase during the initial stage and a smaller increase in the later stage. At the beginning of the experiment, the temperature difference between the fluid and the soil was small that the heat transfer was less, but the heating power was always rated. Therefore, the temperatures of the inlet and outlet of the U-tube increases rapidly at the beginning; with increasing testing time, the fluid temperature increases, the heat transfer efficiency between fluid and the soil increases, and then a dynamic balance is gradually reached between the rated power heater and the cooling power of BHEs. Therefore, the increase rate of the inlet and outlet temperature after 12 h was slow, which tendency is the same as reported by Serageldin et al. [13].
As shown in Figure 8, the temperatures of the inlet and outlet of the circle-shaped tube were 31.07 and 27.77 • C, respectively; 30.64 and 27.34 • C for the curved oval tube, respectively; the outlet temperature decreased by 0.43 • C compared with the circle shaped tube. The inlet and outlet temperatures of the flat oval tube were 30.45 and 27.14 • C, respectively, which is a reduction in the outlet temperature by 0.63 • C compared with the circle-shaped tube. The temperatures of the inlet and outlet of the semicircle tube were 30.76 and 27.46 • C, respectively; the temperature of the outlet was reduced by 0.31 • C compared with the circular tube. The temperatures of the inlet and outlet of the sector tube were 30.58 and 27.28 • C, respectively, which is a decrease in outlet temperature by 0.49 • C compared to the circle-shaped tube. a non-linear increasing trend, with a larger increase during the initial stage and a increase in the later stage. At the beginning of the experiment, the temperature dif between the fluid and the soil was small that the heat transfer was less, but the power was always rated. Therefore, the temperatures of the inlet and outlet of the increases rapidly at the beginning; with increasing testing time, the fluid tempera creases, the heat transfer efficiency between fluid and the soil increases, and the namic balance is gradually reached between the rated power heater and the cooling of BHEs. Therefore, the increase rate of the inlet and outlet temperature after 12 slow, which tendency is the same as reported by Serageldin et al. [13].
As shown in Figure 8, the temperatures of the inlet and outlet of the circletube were 31.07 and 27.77 °C, respectively; 30.64 and 27.34 °C for the curved ov respectively; the outlet temperature decreased by 0.43 °C compared with the circle tube. The inlet and outlet temperatures of the flat oval tube were 30.45 and 27.14 spectively, which is a reduction in the outlet temperature by 0.63 °C compared w circle-shaped tube. The temperatures of the inlet and outlet of the semicircle tub 30.76 and 27.46 °C, respectively; the temperature of the outlet was reduced by compared with the circular tube. The temperatures of the inlet and outlet of the tube were 30.58 and 27.28 °C, respectively, which is a decrease in outlet tempera 0.49 °C compared to the circle-shaped tube. To more accurately evaluate the influence of the cross-sectional shape of a the heat exchange efficiency of BHEs based on outlet fluid temperature, dimens heat exchange efficiency is proposed in combination with the temperature differ the inlet and outlet. The calculation equation is as follows: where Tout is the outlet temperature, °C; Tout,r is the referenced outlet temperature, value is 27.77 °C (temperature of circle-shaped tube); and Tin,r is the referenced in perature, which is 31.07 °C (temperature of circle-shaped tube). Based on Equation (13), the dimensionless heat exchange efficiencies of the shaped, curved oval, flat oval, semicircle, and sector tubes are 100%, 113.0%, 109.4%, and 114.8%, respectively. The descending order is flat oval tube, secto To more accurately evaluate the influence of the cross-sectional shape of a tube on the heat exchange efficiency of BHEs based on outlet fluid temperature, dimensionless heat exchange efficiency is proposed in combination with the temperature difference at the inlet and outlet. The calculation equation is as follows: where T out is the outlet temperature, • C; T out ,r is the referenced outlet temperature, and the value is 27.77 • C (temperature of circle-shaped tube); and T in ,r is the referenced inlet temperature, which is 31.07 • C (temperature of circle-shaped tube). Based on Equation (13), the dimensionless heat exchange efficiencies of the circleshaped, curved oval, flat oval, semicircle, and sector tubes are 100%, 113.0%, 119.1%, 109.4%, and 114.8%, respectively. The descending order is flat oval tube, sector tube, curved oval tube, semicircle tube, and circle shaped tube, which is explained well by the thermal resistance of BHEs with different cross-sectional shapes.

Thermal Resistance of BHEs
The heat resistance is an important parameter for evaluating the heat transfer efficiency of BHEs. Herein, the influence of the cross-sectional shape on the thermal resistance of the BHEs was studied. Considering the complexity of the three-dimensional model, combined with previous research [47,48], a two-dimensional simplified heat resistance analysis model is proposed.
The expression of three-dimensional borehole heat resistance is: where T b is the average temperature of the borehole wall, • C; T p is the average temperature of the U-tube, • C; and q b is the heat flow per unit length of borehole, W/m. The heat transfer between the borehole wall and the tube was only considered in Sharqawy et al.'s method [48] and the corresponding calculating model is shown in Figure  9. Combining the setting law of Sharqawy et al.'s temperature and the ground temperature distribution characteristics in this paper, the boundary temperature of the model was set under the premise of cooling working conditions. The temperature of borehole wall T b , inlet tube T in,int , and outlet tube T out,int were set to 20, 30, and 25 • C, respectively; the material parameters were the same as in Table 2. The temperature in the borehole can reach a stable state within a short time due to the small size of the borehole section, thus, the steady-state analysis was adopted. To more accurately explain the heat exchange efficiency of tubes with different shapes, the model of the heat exchange between tubes (T in,int = 30 • C, T out,int = 25 • C) and the model of heat exchange between the tube and borehole wall (T in,int = 25 • C, T b = 20 • C) were established.

Thermal Resistance of BHEs
The heat resistance is an important parameter for evaluating the heat tran ciency of BHEs. Herein, the influence of the cross-sectional shape on the thermal re of the BHEs was studied. Considering the complexity of the three-dimensiona combined with previous research [47,48], a two-dimensional simplified heat re analysis model is proposed.
The expression of three-dimensional borehole heat resistance is: where Tb is the average temperature of the borehole wall, °C; p T is the average t ture of the U-tube, °C; and qb is the heat flow per unit length of borehole, W/m. The heat transfer between the borehole wall and the tube was only consi Sharqawy et al.'s method [48] and the corresponding calculating model is shown i 9. Combining the setting law of Sharqawy et al.'s temperature and the ground t ture distribution characteristics in this paper, the boundary temperature of the mo set under the premise of cooling working conditions. The temperature of boreh Tb, inlet tube Tin,int, and outlet tube Tout,int were set to 20, 30, and 25 °C, respectiv material parameters were the same as in Table 2. The temperature in the bore reach a stable state within a short time due to the small size of the borehole secti the steady-state analysis was adopted. To more accurately explain the heat excha ciency of tubes with different shapes, the model of the heat exchange between tub = 30 °C, Tout,int = 25 °C) and the model of heat exchange between the tube and boreh (Tin,int = 25 °C, Tb = 20 °C) were established.  Figure 10 plots the heat resistance of BHEs with 3D and 2D models in [48] w ferent section shapes. One can notice that the heat resistance values of the two m are similar. The heat resistances of tubes with different cross-sectional shapes, in ing order, are flat oval, sector, curved oval, semicircle, and circle. The difference resistance with different cross-sectional shapes can be explained in terms of heat e efficiency between two legs and the heat exchange efficiency between the tube a hole wall, which are detailed below.  Figure 10 plots the heat resistance of BHEs with 3D and 2D models in [48] with different section shapes. One can notice that the heat resistance values of the two methods are similar. The heat resistances of tubes with different cross-sectional shapes, in increasing order, are flat oval, sector, curved oval, semicircle, and circle. The differences in heat resistance with different cross-sectional shapes can be explained in terms of heat exchange efficiency between two legs and the heat exchange efficiency between the tube and borehole wall, which are detailed below.
To more accurately explain the heat distribution characteristics inside the borehole, the middle depth of the borehole was selected as a typical borehole section, and the temperature cloud is illustrated in Figure 11. To more accurately explain the heat distribution characteristics inside the boreho the middle depth of the borehole was selected as a typical borehole section, and the te perature cloud is illustrated in Figure 11.  Figure 11 shows the temperature distribution of BHEs with different tube shapes ter 48 h and at a 50 m depth. It can be found that the fluid temperature in the tube is eve distributed and heat exchange is conducted through the tube, backfill soil, borehole w and the soil. Under the same fluid area (531 mm 2 ), a large difference appears at the h exchange perimeter, and the corresponding relationship in decreasing order is flat-o shape (98.9 mm), semicircle shape (94.5 mm), sector shape (92.84 mm), curved oval sha (89.04 mm), and circle shape (81.68 mm). Figure 12 depicts the temperature distribution of BHEs with different tube shapes the given tube temperatures. We found that the temperature distribution of BHEs is d ferent for different tube shapes, which can be used to reflect the heat short-circuit eff between two legs. The description of the heat short-circuit effect in BHEs is that the h exchange occurs between fluid in two legs rather than fluid to soil [15,47], and the ex ence of a heat short-circuit will reduce the heat exchange efficiency of BHEs. Hence, in t study, we recognized the effective evaluation of heat short-circuit between fluid in t legs as the heat flux value between two legs. The heat flux between two legs for cir curved oval, flat oval, semicircle, and sector tubes were 11.133, 11.601, 12.198, 11.836, a 11.327 W/m, respectively, in which the higher heat flux indicates a greater heat short-c cuit effect. The corresponding relationship with the heat short-circuit effect, in decreas order, is flat oval, semicircle, curved oval, sector, and circle tubes. This finding can To more accurately explain the heat distribution characteristics inside the borehole, the middle depth of the borehole was selected as a typical borehole section, and the temperature cloud is illustrated in Figure 11.  Figure 11 shows the temperature distribution of BHEs with different tube shapes after 48 h and at a 50 m depth. It can be found that the fluid temperature in the tube is evenly distributed and heat exchange is conducted through the tube, backfill soil, borehole wall, and the soil. Under the same fluid area (531 mm 2 ), a large difference appears at the heat exchange perimeter, and the corresponding relationship in decreasing order is flat-oval shape (98.9 mm), semicircle shape (94.5 mm), sector shape (92.84 mm), curved oval shape (89.04 mm), and circle shape (81.68 mm). Figure 12 depicts the temperature distribution of BHEs with different tube shapes at the given tube temperatures. We found that the temperature distribution of BHEs is different for different tube shapes, which can be used to reflect the heat short-circuit effect between two legs. The description of the heat short-circuit effect in BHEs is that the heat exchange occurs between fluid in two legs rather than fluid to soil [15,47], and the existence of a heat short-circuit will reduce the heat exchange efficiency of BHEs. Hence, in this study, we recognized the effective evaluation of heat short-circuit between fluid in two legs as the heat flux value between two legs. The heat flux between two legs for circle, curved oval, flat oval, semicircle, and sector tubes were 11.133, 11.601, 12.198, 11.836, and 11.327 W/m, respectively, in which the higher heat flux indicates a greater heat short-circuit effect. The corresponding relationship with the heat short-circuit effect, in decreasing order, is flat oval, semicircle, curved oval, sector, and circle tubes. This finding can be explained by the combined effect of the heat exchange perimeter of the tube and the relative distance between tubes [12,13]. The larger the heat exchange perimeter of a tube and the smaller relative distance between tubes, the more obvious the heat short-circuit effect.  Figure 11 shows the temperature distribution of BHEs with different tube shapes after 48 h and at a 50 m depth. It can be found that the fluid temperature in the tube is evenly distributed and heat exchange is conducted through the tube, backfill soil, borehole wall, and the soil. Under the same fluid area (531 mm 2 ), a large difference appears at the heat exchange perimeter, and the corresponding relationship in decreasing order is flat-oval shape (98.9 mm), semicircle shape (94.5 mm), sector shape (92.84 mm), curved oval shape (89.04 mm), and circle shape (81.68 mm). Figure 12 depicts the temperature distribution of BHEs with different tube shapes at the given tube temperatures. We found that the temperature distribution of BHEs is different for different tube shapes, which can be used to reflect the heat short-circuit effect between two legs. The description of the heat short-circuit effect in BHEs is that the heat exchange occurs between fluid in two legs rather than fluid to soil [15,47], and the existence of a heat short-circuit will reduce the heat exchange efficiency of BHEs. Hence, in this study, we recognized the effective evaluation of heat short-circuit between fluid in two legs as the heat flux value between two legs. The heat flux between two legs for circle, curved oval, flat oval, semicircle, and sector tubes were 11.133, 11.601, 12.198, 11.836, and 11.327 W/m, respectively, in which the higher heat flux indicates a greater heat short-circuit effect. The corresponding relationship with the heat short-circuit effect, in decreasing order, is flat oval, semicircle, curved oval, sector, and circle tubes. This finding can be explained by the combined effect of the heat exchange perimeter of the tube and the relative distance between tubes [12,13]. The larger the heat exchange perimeter of a tube and the smaller relative distance between tubes, the more obvious the heat short-circuit effect. Sustainability 2021, 13, x FOR PEER REVIEW 13 of 21  Figure 13 indicates the temperature distribution of BHEs with different-shaped tube at the given tube and borehole wall temperatures. We found that the temperature distribution and heat transfer direction are all different at different tube shapes. Due to the differences in tube geometry and distribution characteristics, the heat flux values between the borehole wall and circle, curved oval, flat oval, semicircle, and sector tubes were: 28.8, 30.8, 33.2, 30.2, and 30.5 W/m, respectively, in which a higher heat flux indicates a greater heat exchange efficiency of the tube to the borehole wall. The corresponding affecting relationship of the heat exchange efficiency of the tube to borehole wall, in decreasing order, is flat oval, curved oval, sector, semicircle, and circle tubes. These can be explained by the combined effect of the heat exchange perimeter of tube and relative distance of the tube to borehole wall. The larger the heat exchange perimeter of a tube and the smaller relative distance of tube to borehole wall, the higher heat exchange efficiency of tube to borehole wall. Based on Sharqawy et al.'s theory [48] the heat exchange efficiency of BHEs is related to the heat exchange efficiency of the tube to the borehole wall and the heat short-circuit effect between two legs. The greater the heat exchange efficiency of the tube to the borehole wall and the less the heat short-circuit effect between legs, the greater the heat exchange efficiency of BHEs (less heat resistance).
For flat oval tubes, as shown in Figure 14c, the heat flux between the tube and the borehole wall was the greatest, much larger than that of other shapes. Therefore, even if the heat short-circuit effect of the flat oval tube is obvious, it still has a higher heat exchange efficiency, which means a lower heat resistance of 1.251 m·K/W.  Figure 13 indicates the temperature distribution of BHEs with different-shaped tube at the given tube and borehole wall temperatures. We found that the temperature distribution and heat transfer direction are all different at different tube shapes. Due to the differences in tube geometry and distribution characteristics, the heat flux values between the borehole wall and circle, curved oval, flat oval, semicircle, and sector tubes were: 28.8, 30.8, 33.2, 30.2, and 30.5 W/m, respectively, in which a higher heat flux indicates a greater heat exchange efficiency of the tube to the borehole wall. The corresponding affecting relationship of the heat exchange efficiency of the tube to borehole wall, in decreasing order, is flat oval, curved oval, sector, semicircle, and circle tubes. These can be explained by the combined effect of the heat exchange perimeter of tube and relative distance of the tube to borehole wall. The larger the heat exchange perimeter of a tube and the smaller relative distance of tube to borehole wall, the higher heat exchange efficiency of tube to borehole wall.  Figure 13 indicates the temperature distribution of BHEs with different-shaped tube at the given tube and borehole wall temperatures. We found that the temperature distribution and heat transfer direction are all different at different tube shapes. Due to the differences in tube geometry and distribution characteristics, the heat flux values between the borehole wall and circle, curved oval, flat oval, semicircle, and sector tubes were: 28.8, 30.8, 33.2, 30.2, and 30.5 W/m, respectively, in which a higher heat flux indicates a greater heat exchange efficiency of the tube to the borehole wall. The corresponding affecting relationship of the heat exchange efficiency of the tube to borehole wall, in decreasing order, is flat oval, curved oval, sector, semicircle, and circle tubes. These can be explained by the combined effect of the heat exchange perimeter of tube and relative distance of the tube to borehole wall. The larger the heat exchange perimeter of a tube and the smaller relative distance of tube to borehole wall, the higher heat exchange efficiency of tube to borehole wall. Based on Sharqawy et al.'s theory [48] the heat exchange efficiency of BHEs is related to the heat exchange efficiency of the tube to the borehole wall and the heat short-circuit effect between two legs. The greater the heat exchange efficiency of the tube to the borehole wall and the less the heat short-circuit effect between legs, the greater the heat exchange efficiency of BHEs (less heat resistance).
For flat oval tubes, as shown in Figure 14c, the heat flux between the tube and the borehole wall was the greatest, much larger than that of other shapes. Therefore, even if the heat short-circuit effect of the flat oval tube is obvious, it still has a higher heat exchange efficiency, which means a lower heat resistance of 1.251 m·K/W. Based on Sharqawy et al.'s theory [48] the heat exchange efficiency of BHEs is related to the heat exchange efficiency of the tube to the borehole wall and the heat short-circuit effect between two legs. The greater the heat exchange efficiency of the tube to the borehole wall and the less the heat short-circuit effect between legs, the greater the heat exchange efficiency of BHEs (less heat resistance).
For flat oval tubes, as shown in Figure 14c, the heat flux between the tube and the borehole wall was the greatest, much larger than that of other shapes. Therefore, even if the heat short-circuit effect of the flat oval tube is obvious, it still has a higher heat exchange efficiency, which means a lower heat resistance of 1.251 m·K/W.
For circle-shaped tubes, as shown in Figure 14a, the heat exchange efficiency of the tube to the borehole wall was the lowest, being much lower than that of other shapes. Therefore, even if the heat short-circuit effect of the circular tube is the lowest, it still has a smaller heat exchange efficiency, which means a greater heat resistance of 1.433 m·K/W.
For semicircle-shaped tubes, as shown in Figure 14d, the heat exchange efficiency of the tube to the borehole wall and the heat short-circuit effect of the tube were both larger than those of the circular tube. However, the influence of heat short-circuit between two legs was less than that of the heat exchange efficiency of the tube to the borehole wall.
Hence, the tube with semicircle shape still has a higher heat exchange efficiency than the circular tube, resulting in a lower heat resistance of 1.360 m·K/W. For circle-shaped tubes, as shown in Figure 14a, the heat exchange efficiency of the tube to the borehole wall was the lowest, being much lower than that of other shapes. Therefore, even if the heat short-circuit effect of the circular tube is the lowest, it still has a smaller heat exchange efficiency, which means a greater heat resistance of 1.433 m·K/W.
For semicircle-shaped tubes, as shown in Figure 14d, the heat exchange efficiency of the tube to the borehole wall and the heat short-circuit effect of the tube were both larger than those of the circular tube. However, the influence of heat short-circuit between two legs was less than that of the heat exchange efficiency of the tube to the borehole wall. Hence, the tube with semicircle shape still has a higher heat exchange efficiency than the circular tube, resulting in a lower heat resistance of 1.360 m·K/W.
Comparing the curved oval and sector tubes, as shown in Figure 14b,e, we found that the heat exchange efficiency of the sector-shaped tube to the borehole wall was larger than that of the curved oval tube, and the heat exchange efficiency between two legs with a sector-shaped tube was lower than for curved oval tubes. Under the combined effect of the both, the heat exchange efficiency of the sector-shaped tubes in BHEs was greater than that of the curved oval tube. In other words, the heat resistance of the tube with a sector shape (1.318 m·K/W) was lower than that of the curved oval tube (1.344 m·K/W).
As shown in Figure 10, we found that the calculated heat resistance with Sharqawy et al.'s method [48] is similar to that calculated using the 3D model. The calculation errors of the two methods were less than 2%, which means that with constant legs distance, Sharqawy et al.'s method can be used to estimate the heat resistance considering different tube shapes [48].

Temperature Distribution of U-Tube
Serageldin et al. [13,39] stated that during the BHEs cooling process, the drop in fluid temperature is non-uniform, which is related to the soil temperature. Due to the shorter distance between two BHE legs, the heat short-circuit effect affects the temperature distribution between two legs. The temperature distribution of BHEs with different tube shapes along the depth at 48 h is shown in Figure 15.  Comparing the curved oval and sector tubes, as shown in Figure 14b,e, we found that the heat exchange efficiency of the sector-shaped tube to the borehole wall was larger than that of the curved oval tube, and the heat exchange efficiency between two legs with a sector-shaped tube was lower than for curved oval tubes. Under the combined effect of the both, the heat exchange efficiency of the sector-shaped tubes in BHEs was greater than that of the curved oval tube. In other words, the heat resistance of the tube with a sector shape (1.318 m·K/W) was lower than that of the curved oval tube (1.344 m·K/W).
As shown in Figure 10, we found that the calculated heat resistance with Sharqawy et al.'s method [48] is similar to that calculated using the 3D model. The calculation errors of the two methods were less than 2%, which means that with constant legs distance, Sharqawy et al.'s method can be used to estimate the heat resistance considering different tube shapes [48].

Temperature Distribution of U-Tube
Serageldin et al. [13,39] stated that during the BHEs cooling process, the drop in fluid temperature is non-uniform, which is related to the soil temperature. Due to the shorter distance between two BHE legs, the heat short-circuit effect affects the temperature distribution between two legs. The temperature distribution of BHEs with different tube shapes along the depth at 48 h is shown in Figure 15.
tube to the borehole wall was the lowest, being much lower than th Therefore, even if the heat short-circuit effect of the circular tube is th a smaller heat exchange efficiency, which means a greater heat resista For semicircle-shaped tubes, as shown in Figure 14d, the heat exc the tube to the borehole wall and the heat short-circuit effect of the tu than those of the circular tube. However, the influence of heat short-c legs was less than that of the heat exchange efficiency of the tube to Hence, the tube with semicircle shape still has a higher heat exchange circular tube, resulting in a lower heat resistance of 1.360 m·K/W.
Comparing the curved oval and sector tubes, as shown in Figure 1 the heat exchange efficiency of the sector-shaped tube to the borehole w that of the curved oval tube, and the heat exchange efficiency betwe sector-shaped tube was lower than for curved oval tubes. Under the the both, the heat exchange efficiency of the sector-shaped tubes in BH that of the curved oval tube. In other words, the heat resistance of th shape (1.318 m·K/W) was lower than that of the curved oval tube (1.34 As shown in Figure 10, we found that the calculated heat resista et al.'s method [48] is similar to that calculated using the 3D model. Th of the two methods were less than 2%, which means that with constant qawy et al.'s method can be used to estimate the heat resistance consid shapes [48].

Temperature Distribution of U-Tube
Serageldin et al. [13,39] stated that during the BHEs cooling proce temperature is non-uniform, which is related to the soil temperature distance between two BHE legs, the heat short-circuit effect affects the bution between two legs. The temperature distribution of BHEs with d along the depth at 48 h is shown in Figure 15.   Figure 15 shows that the temperature at the inlet side drops more rapidly than on the outlet side with increasing depth, and fewer temperature changes appear when the depth is less than 5 m. In the BHE cooling process, hotter fluid flows into the tube after experiencing the heat exchange with soil, then flows out from the outlet. According to the heat exchange law (heat is always exchanged from high to low temperature) and the temperature distribution results in Figure 15, we concluded that the temperature on the inlet side T in is higher than the temperature on the outlet side T out , and both are higher than the soil temperature T soil . On the inlet side, the fluid heat transfers to the soil and to the fluid on the outlet side; on the outlet side, the fluid heat transfers to the soil and to the fluid on the inlet side, but the higher temperature on the inlet side will reduce the heat transfer efficiency of the fluid on the outlet side due to the heat short-circuit effect. Hence, the heat exchange of the fluid on the inlet side is more efficient than that of the fluid on the outlet side, undergoing a higher temperature drop at the inlet side than on the outlet side. Furthermore, less of a temperature change was observed within a depth of 5 m. This can be explained by the higher soil temperature (shown in Figure 2) that reduces the heat exchange efficiency of fluid to soil. Therefore, when designing the length of BHEs, the depth of soil that is easily affected by seasonal temperature should be considered to maintain the BHE efficiency. Figure 16 depicts the heat flow of two legs with different shapes. We found that the heat flux on the inlet side and outlet side differed due to the heat exchange between two tubes. On the inlet side, the heat fluxes of circle, curved oval, flat oval, semicircle, and sector tubes were 27.93, 28.22, 28.99, 28.54, and 28.13 W/m, respectively; on the outlet side, the heat fluxes were 17.88, 17.59, 16.82, 17.27, and 17.68 W/m, respectively. According to the discussion in Section 4.3 and Serageldin's studies, the larger the heat flux difference between the inlet and outlet sides, the more obvious the heat short-circuit effect. Hence, we concluded that the relationship between the heat short-circuit and tube shapes in decreasing order is flat oval shape, semicircle shape, curved oval shape, sector shape, and circle shape; the tendency is consistent with the discussion in Section 4.3. Figure 15 shows that the temperature at the inlet side drops more rapidly than on outlet side with increasing depth, and fewer temperature changes appear when the d is less than 5 m. In the BHE cooling process, hotter fluid flows into the tube after ex encing the heat exchange with soil, then flows out from the outlet. According to the exchange law (heat is always exchanged from high to low temperature) and the temp ture distribution results in Figure 15, we concluded that the temperature on the inlet Tin is higher than the temperature on the outlet side Tout, and both are higher than the temperature Tsoil. On the inlet side, the fluid heat transfers to the soil and to the flui the outlet side; on the outlet side, the fluid heat transfers to the soil and to the fluid on inlet side, but the higher temperature on the inlet side will reduce the heat transfer ciency of the fluid on the outlet side due to the heat short-circuit effect. Hence, the exchange of the fluid on the inlet side is more efficient than that of the fluid on the o side, undergoing a higher temperature drop at the inlet side than on the outlet side. thermore, less of a temperature change was observed within a depth of 5 m. This ca explained by the higher soil temperature (shown in Figure 2) that reduces the hea change efficiency of fluid to soil. Therefore, when designing the length of BHEs, the d of soil that is easily affected by seasonal temperature should be considered to main the BHE efficiency. Figure 16 depicts the heat flow of two legs with different shapes. We found tha heat flux on the inlet side and outlet side differed due to the heat exchange between tubes. On the inlet side, the heat fluxes of circle, curved oval, flat oval, semicircle, sector tubes were 27.93, 28.22, 28.99, 28.54, and 28.13 W/m, respectively; on the outlet the heat fluxes were 17.88, 17.59, 16.82, 17.27, and 17.68 W/m, respectively. Accordin the discussion in Section 4.3 and Serageldin's studies, the larger the heat flux differ between the inlet and outlet sides, the more obvious the heat short-circuit effect. He we concluded that the relationship between the heat short-circuit and tube shapes in creasing order is flat oval shape, semicircle shape, curved oval shape, sector shape, circle shape; the tendency is consistent with the discussion in Section 4.3.

Borehole Wall Temperature
In the engineering of BHEs, due to the influence of fluid temperature, tube sha and distribution characteristics of the tube, the temperature of the borehole wall is usu non-uniform, which plays an important role in determining the heat affecting rad Hence, the temperature distribution at the borehole wall must be studied considering ferent tube shapes. In the analysis, to more accurately determine the influence of

Borehole Wall Temperature
In the engineering of BHEs, due to the influence of fluid temperature, tube shapes, and distribution characteristics of the tube, the temperature of the borehole wall is usually nonuniform, which plays an important role in determining the heat affecting radius. Hence, the temperature distribution at the borehole wall must be studied considering different tube shapes. In the analysis, to more accurately determine the influence of tube shape on borehole wall temperature, referring to the temperature distribution characteristics in Figure 13, the fluid temperature on the inlet and outlet sides were set to constant values of 30 and 28 • C, respectively, and the corresponding model and results are shown in Figure 17.
shape on borehole wall temperature, referring to the temperature distribution ch istics in Figure 13, the fluid temperature on the inlet and outlet sides were set to values of 30 and 28 °C, respectively, and the corresponding model and results ar in Figure 17.  Figure 17 shows that the temperature distribution along the borehole wal uniform and the tendency is similar for different tube shapes. Comparing the m and minimum temperatures of the borehole wall with different tube shapes, w that temperature difference between the maximum and minimum along the boreh in decreasing order is circle, sector, curved oval, flat oval, and semicircle tubes. T ence of tube shape on borehole wall temperature can be explained in terms of exchange perimeter and the heat resistance of the tube to the borehole wall.

Control of Leg Distance
During the sinking of the U-tube, to ensure the distance between the two leg device is generally used. Figure 16 shows two typical U-tube limit devices. Figure fixed limit device with constant spacing, which was used in the TRTs. In engineeri the drilling is completed, the soil with low strength inside the borehole springs bac the action of lateral soil pressure, which means the borehole diameter will d Hence, if using a fixed limit device, for overcoming the problem of the U-tube st the borehole, the designed distance between two legs of a U-tube usually is less borehole. Figure 18b depicts an elastic limit device that can make U-type legs cont the borehole wall, which maximizes the distance between the two legs. Howev borehole wall is weak or the backfill soil is late, it is possible to lose the desired ef to the elastic force being lower than the lateral soil force. The above two limit de both common in BHEs engineering; hence, to better support engineering of BH maximum leg distance of tubes with different shapes was considered and the corr ing results are shown in Figures 19 and 20.  Figure 17 shows that the temperature distribution along the borehole wall is nonuniform and the tendency is similar for different tube shapes. Comparing the maximum and minimum temperatures of the borehole wall with different tube shapes, we found that temperature difference between the maximum and minimum along the borehole wall in decreasing order is circle, sector, curved oval, flat oval, and semicircle tubes. The influence of tube shape on borehole wall temperature can be explained in terms of the heat exchange perimeter and the heat resistance of the tube to the borehole wall.

Control of Leg Distance
During the sinking of the U-tube, to ensure the distance between the two legs, a limit device is generally used. Figure 16 shows two typical U-tube limit devices. Figure 18a is a fixed limit device with constant spacing, which was used in the TRTs. In engineering, after the drilling is completed, the soil with low strength inside the borehole springs back under the action of lateral soil pressure, which means the borehole diameter will decrease. Hence, if using a fixed limit device, for overcoming the problem of the U-tube sticking in the borehole, the designed distance between two legs of a U-tube usually is less than the borehole. Figure 18b depicts an elastic limit device that can make U-type legs contact with the borehole wall, which maximizes the distance between the two legs. However, if the borehole wall is weak or the backfill soil is late, it is possible to lose the desired effect due to the elastic force being lower than the lateral soil force. The above two limit devices are both common in BHEs engineering; hence, to better support engineering of BHEs, the maximum leg distance of tubes with different shapes was considered and the corresponding results are shown in Figures 19 and 20.     Figure 19 shows that the heat resistance at maximum legs distance is lower than a constant leg distance. Comparing the temperatures in Figures 11 and 20, we concluded that the greater the distance between the two legs, the higher the heat exchange efficiency between the fluid and the smaller the heat short-circuit effect. Hence, heat resistance de creases with the increase in leg distance.

Heat Resistance with Different Leg Distances
With increasing leg distance, the heat resistance reduced from 1. 43 Figure 19 shows that the heat resistance at maximum legs distance is lower than at constant leg distance. Comparing the temperatures in Figures 11 and 20, we concluded that the greater the distance between the two legs, the higher the heat exchange efficiency between the fluid and the smaller the heat short-circuit effect. Hence, heat resistance decreases with the increase in leg distance.

Heat Resistance with Different Leg Distances
With increasing leg distance, the heat resistance reduced from 1.   Figure 19 shows that the heat resistance at maximum legs distance is lower than at constant leg distance. Comparing the temperatures in Figures 11 and 20, we concluded that the greater the distance between the two legs, the higher the heat exchange efficiency between the fluid and the smaller the heat short-circuit effect. Hence, heat resistance decreases with the increase in leg distance.

Heat Resistance with Different Leg Distances
With increasing leg distance, the heat resistance reduced from 1.43 to 1.1 for circle tubes, about 23.2%; reduced from 1.34 to 0.96 for curved oval tubes, about 28.4%; reduced from 1.25 to 0.91 for flat oval tubes, about 27.7%, reduced from 1.36 to 0.99 for semicircle tubes, about 27.4%,; and for sector-shaped tubes, reduced from 1.32 to 1.01, about 23.1%. Hence, the influence of leg distance and tube shapes in reducing order is curved oval, flat oval shape, semicircle, circle, and sector tubes.
Comparing the heat resistance of tubes with different shapes, it can be noted that at maximum leg distance, the heat resistance of flat-oval-shaped tubes was the lowest (0.91 m·K/W), which is a 17.7% reduction compared to circle-shaped tubes; the heat resistance of curved-oval-shaped tube is 0.96 m·K/W, which is a 12.6% reduction compared to the circle-shaped tube; 0.99 m·K/W for the semicircle-shaped tube, which is a 10.3% reduction compared to the circular tube; 1.01 m·K/W for the sector-shaped tube, which is a 7.8% reduction compared to the circle-shaped tube; and the heat resistance of the sector-shaped tube was the biggest, 1.10 m·K/W. Hence, at the maximum leg distance, the heat resistance of the different-shaped tubes in increasing order is flat oval, curved oval, semicircle, sector, and circle.
Based on above analysis, we concluded that the adjustment of the leg distance will affect the heat resistance by about 25% and the adjustment of the tube shape will affect the heat resistance by about 12%. Hence, in engineering design, the leg distance should be considered first, followed by changing tube shape; the flat-oval-shaped tube at maximum leg distance is the best tube design for BHEs of the designs considered.

Conclusions
For improving the heat exchange efficiency of BHEs and reducing initial investment, based on TRTs and COMSOL Multi-physics, the heat transfer characteristics of BHEs with different cross-sectional tubes (circular, curved ellipse, straight ellipse, semicircle, and sector) were studied using 3D numerical analysis. The corresponding conclusions are described as follows: (1) Based on TRTs, we verified the reliability of using COMSOL Multi-physics with three-dimensional models to study the heat exchange characteristic of BHEs. (2) Compared with the inlet and outlet temperatures of the tube with a circle shape at 48 h, the results show that the heat exchange efficiency of the curved oval tube increases by 13.0%, the flat oval tube by 19.1%, the semicircle tube by 9.4%, and sector tube by 14.8%. (3) The heat flux heterogeneity of the tube on the inlet and outlet sides in BHEs with different tube shapes in decreasing order is flat oval shape, semicircle shape, curved oval shape, sector shape, and circle shape. (4) The temperature heterogeneity of the borehole wall in BHEs with different tube shapes in decreasing order is circle shape, sector shape, curved oval shape, flat oval shape, and semicircle shape. (5) Compared with the heat resistance of the circle-shaped tube at 48 h, we found that the heat resistance of the curved oval tube decreased by 12.6%, the flat oval tube by 17.7%, the semicircle tube by 10.3%, and the sector tube by 7.8% at maximum leg distance. (6) The adjustment of the leg distance affects the heat resistance by about 25% and the adjustment of the tube shape affects the heat resistance by about 12%. (7) The adjustment of tube shape can improve the heat exchange efficiency of BHEs, and the flat-oval-shaped tube at maximum leg distance was found to be the best tube design for BHEs.
It is worth noting that the influence of tube cross-section shape on heat transfer characteristics of BHEs at 48 h was studied by combining TRT and a 3D model in this work. To better serve the project about BHEs, long-term BHE heat transfer characteristics (1-10 years) and cyclic loading (cooling and heating) will be investigated according to monitoring data in future research.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.