Numerical Analysis of the Factors Influencing a Vertical U-Tube Ground Heat Exchanger

The development of a three-dimensional, unsteady state model, which couples heat transfer with groundwater seepage for a vertical U-tube ground heat exchanger (GHE) is presented. The influence of underground soil thermal properties, grout materials, inlet water temperature and velocity, and groundwater seepage on heat transfer in the GHE is examined. The results indicate that before the heat in the borehole is saturated, the heat flux in the GHE is directly proportional to the thermal conductivity coefficient of the grout materials. The radius of the thermal effect of the GHE and the recovery rate of the temperature in the soil are also proportional to the thermal diffusion coefficient of the soil. In cooling mode, the increase of the inlet water temperature of the GHE results in enhanced heat transfer. However, this may cause issues with heat buildup. The increase of the inlet water velocity in the GHE enhances heat convection in the tube. The effect of thermal-seepage coupling in groundwater can reduce the accumulated heat, thus resulting in the effective enhancement of the heat transfer in the GHE.


Introduction
Ground-source heat-pump (GSHP) systems have been widely applied due to their advantages in efficiency, energy conservation, and emission reduction.The GSHP systems use geothermal energy from shallow earth layers as a heat source or a heat sink, and absorb or emit heat to the underground through ground heat exchangers (GHE) to provide refrigeration or heat supply to buildings [1][2][3][4].Their first application was in the USA in 1945.The GSHP, as an innovative device, succeeded in revealing its great advantages [5].The ability of a GHE mainly depends on the structure and thermal properties of the surroundings as well as on the operating conditions.To ensure the performance of a GSHP, the careful design of GHEs is essential [6,7].
Various types of GHEs are available on the market.The vertical U-tubes, which are inserted into a borehole and often backfilled with cement, bentonite and sand, are very efficient, but the drilling cost is rather high.Thus, a lot of researchers have aimed to improve the performance of GHEs and reduce the borehole depth [8].Jalaluddin and Miyara [9] built three GHE models, including U-tube, double-tube, and multi-tube, to investigate the effect of tube number.Based on simulation results, during cooling mode, the heat flux of GHE at 22 h increased by 23.7% for the double-tube and 54.2% for the multi-tube.This indicated that the increase in tube number provided a possibility to improve the thermal performance and reduce the drilling cost of the GHEs.
The horizontal and helically coiled GHEs are also good choices for high-efficiency and low-cost [10,11].Congedo et al. [12] investigated three different tubes (linear, helical and slinky) that influence the performance of GHE by CFD simulations, and the analysis revealed that the helical heat exchanger arrangement performed the best.
The shank spacing is also one of the structural factors on which this study is focused.In order to investigate the factors influencing thermal conductivity and thermal borehole resistance, a sensitivity study based on numerical simulations, considering the effects of the shank spacing, initial thermal distribution and thermal dispersivity, was performed by Wagner et al. [13].It was demonstrated that the effects of the first two parameters were just below 10%.On the other hand, thermal dispersivity plays a key role because its influence overestimates the effective thermal conductivity by a factor of 1.2-2.9 as compared to the saturated condition.
Bouhacina et al. [14] added longitudinal fins to the inner surface of the U-tube to obtain a novel GHE, and then the effects and benefits of this novel GHE were studied numerically.Two models were built to simulate the thermal and dynamic behaviors of a traditional GHE and the novel GHE.The result showed that under the same inlet condition, a greater velocity was noticed in the case of the novel GHE, and the heat exchange was up to 7% more efficient.
Thermal properties have also been considered.Three parameters of the GHE, including subsurface thermophysical properties, pipe materials and operating modes, were investigated by Lous et al. [15].The results showed that the porosity, thermal conductivities, and geothermal flux play an important role in the heat transfer.Furthermore, the parameters that affect thermal comportment are primarily the thermal conductivity of the U-tube, the discharge rate of the circulating water and the heat flux.
Several studies have taken the influence of convection into consideration.A CFD simulation of the U-tube was considered by Gustafsson et al. [16].A temperature gradient was obtained by setting a lower temperature on the far-field wall and a higher temperature on the tube wall, which induced a velocity gradient in the underground water.This resulted in an increased heat flux.Yang et al. [17] developed a GHE model that accounts for groundwater advection in a porous media to investigate the influence of groundwater on the thermal distribution.It was discovered that better diffusivity conditions can increase the rate of soil heat diffusion.Thus, the heat buildup can be alleviated.In addition, groundwater advection can efficiently remove the underground heat buildup.Dehkordi and Schincariol [18] used a fully discretized finite-element model to evaluate the thermal effect and hydrogeological properties of a GHE.The sensitivity of major parameters and thermo-hydrogeological factors that influence the system's performance, as well as their impact, were analyzed over 6-month and 25-year operation periods.The results indicated that the groundwater flow (>10 −7 m/s), thermal conductivity and thermal gradient between the background and inlet result in a more efficient heat exchange.However, the effect of thermal interference between the branch tubes was overlooked while analyzing the grout thermal conductivity.It can be proved that the thermal interference would be enhanced, and then the GHE efficiency reduced by using thermally enhanced grout.
However, these studies did not thoroughly consider the influencing factors; therefore a comparative analysis is required.In this paper, the vertical GHE applied under different scenarios is evaluated.A three-dimensional unsteady state model for heat transfer in the vertical U-tube GHE is developed, and the influence of underground soil thermal properties, grout materials, inlet water temperature and velocity, and groundwater seepage on the heat transfer in the GHE is analyzed.In the discussion on heat flux, the factors are investigated from different aspects in detail such as the thermal-seepage coupling, thermal interference between branch tubes, turbulent flow in tubes, and recovery of ground temperature.

Model Description
As shown in Figure 1, a U-tube with an outer diameter of 37 mm and thickness of 2.5 mm is inserted into a 150 mm borehole at a depth of 100 m, while the distance between the two tubes is 50 mm.The assumptions that have been applied in the model are: (i) the physical properties of both soil and grout are homogenous and isotropic; (ii) the contact thermal resistances of joints are negligible; (iii) the soil is a saturated porous medium and (iv) the underground water seepage flow is horizontal.All the simulations were performed using the commercial Fluent 6.3 CFD software (Fluent, Inc., New York, NY, USA) and the Gambit 2.0 preprocessor (Fluent, Inc.).and grout are homogenous and isotropic; (ii) the contact thermal resistances of joints are negligible; (iii) the soil is a saturated porous medium and (iv) the underground water seepage flow is horizontal.All the simulations were performed using the commercial Fluent 6.3 CFD software (Fluent, Inc., New York, NY, USA) and the Gambit 2.0 preprocessor (Fluent, Inc.).It is reasonable to refine the mesh along the horizontal plane and coarsen that along the vertical plane.As shown in Figure 2, the ground and grout domain is meshed using triangular elements, the fluid domain is meshed using quadrangular elements, the U-tube is meshed using boundary layer elements, the curving section is meshed with 15 interval counts and the overall volume grid is generated by sweeping the model with a space of 0.5 m.A mesh convergence study was performed by comparing the simulations and the Chen et al. [19] experiment.The results in Table 1 show the difference in heat flux based on three types of mesh density.The finer mesh produced the minimum error (0.6%) and was adopted to simulate the GHE.The model is able to predict the heat flux very accurately.It is reasonable to refine the mesh along the horizontal plane and coarsen that along the vertical plane.As shown in Figure 2, the ground and grout domain is meshed using triangular elements, the fluid domain is meshed using quadrangular elements, the U-tube is meshed using boundary layer elements, the curving section is meshed with 15 interval counts and the overall volume grid is generated by sweeping the model with a space of 0.5 m. (iii) the soil is a saturated porous medium and (iv) the underground water seepage flow is horizontal.All the simulations were performed using the commercial Fluent 6.3 CFD software (Fluent, Inc., New York, NY, USA) and the Gambit 2.0 preprocessor (Fluent, Inc.).It is reasonable to refine the mesh along the horizontal plane and coarsen that along the vertical plane.As shown in Figure 2, the ground and grout domain is meshed using triangular elements, the fluid domain is meshed using quadrangular elements, the U-tube is meshed using boundary layer elements, the curving section is meshed with 15 interval counts and the overall volume grid is generated by sweeping the model with a space of 0.5 m.A mesh convergence study was performed by comparing the simulations and the Chen et al. [19] experiment.The results in Table 1 show the difference in heat flux based on three types of mesh density.The finer mesh produced the minimum error (0.6%) and was adopted to simulate the GHE.The model is able to predict the heat flux very accurately.A mesh convergence study was performed by comparing the simulations and the Chen et al. [19] experiment.The results in Table 1 show the difference in heat flux based on three types of mesh density.The finer mesh produced the minimum error (0.6%) and was adopted to simulate the GHE.The model is able to predict the heat flux very accurately.

The Governing Equations
The generalized equation can be expressed as [20]: The k and ε realizable equations, which characterize the turbulent flow, can be expressed as follows [21]: where µ, ν, and µ t are the dynamic viscosity, kinematic viscosity, and turbulent viscosity, respectively; G k and E are the production terms of k and ε; σ k and σ ε are the Prandtl number of k and ε with values of 1.0 and 1.2, respectively; C 1 and C 2 are empirical constants.
The momentum equation of the soil has an additional source item, which is the viscosity loss of underground water [22]; the equation can be expressed as: where S u is the source term; α and u are the soil permeability and groundwater flow velocity, respectively.
The energy equations of the soil can be divided into two parts: One part is the solid phase of soil: The other part is the fluid phase of underground water: The total equation can be expressed as [23]: where ϕ is the soil porosity.

Data Processing
The initial and boundary conditions are set as follows: (1) The initial temperature is a constant of 290.6 K.
(2) The far-field and bottom surfaces are set as undisturbed boundaries.
(3) The top ground surface is set as a convective boundary with an air temperature of 289.0 K and a heat transfer coefficient of 10 W/(m 2 •K).(4) The inlet is set as a velocity-inlet boundary, and the outlet is set as a pressure-outlet boundary.(5) The reference parameters for simulation are listed in Table 2.The overall heat transfer capacity of GHE can be expressed as Equation ( 8): where q l is the heat flux per unit borehole depth; ρc p f is the volumetric heat capacity; v in is the inlet flow; T in and T out are the temperatures of inlet and outlet, respectively; H is the borehole depth.

Influence of Grout Materials
Cement, bentonite, or sand are generally used as grout materials in engineering practice; the thermophysical properties of these materials vary significantly.Three different grout materials are selected for the analyses and their parameters are shown in Table 3.The heat flux histories for the different grout materials are shown in Figure 3.During the initial period, heat transfer in the GHE is concentrated in the inside part of the borehole.The heat flux is determined by the thermophysical properties of the grout materials, and is found to be the highest when cement is used for grout.After some time, the heat is accumulated to a certain quantity and the heat flux is determined by the thermophysical properties of the soil outside the borehole.Therefore, the difference in the heat flux is reduced after 10 h.
As shown in Figure 4, the fluid temperature in the U-tube decreases along the fluid path.The reduction in the temperature reaches its highest value when the cement grout is used.The temperature of sand and cement grout increases near the outlet tube, which has an unfavorable effect on the heat transfer because the heat in the inlet tube is transferred to the outlet tube by the grout materials within a narrow distance (50 mm).Therefore, the temperature in the outlet tube rises again and its recovery depends on how high the thermal conductivity coefficient is of the grout materials.

Influence of Soil
Three different soils are selected and the parameters are listed in Table 4.The heat flux histories are shown in Figure 5.The greater the thermal conductivity coefficient of the soil, the higher the heat flux of the GHE.The heat flux reaches its highest value when the ground materials are composed of granites.At the beginning, the difference in the heat flux is small and it increases with time as the heat diffuses towards the soil regions.

Influence of Soil
Three different soils are selected and the parameters are listed in Table 4.The heat flux histories are shown in Figure 5.The greater the thermal conductivity coefficient of the soil, the higher the heat flux of the GHE.The heat flux reaches its highest value when the ground materials are composed of granites.At the beginning, the difference in the heat flux is small and it increases with time as the heat diffuses towards the soil regions.

Influence of Soil
Three different soils are selected and the parameters are listed in Table 4.The heat flux histories are shown in Figure 5.The greater the thermal conductivity coefficient of the soil, the higher the heat flux of the GHE.The heat flux reaches its highest value when the ground materials are composed of granites.At the beginning, the difference in the heat flux is small and it increases with time as the heat diffuses towards the soil regions.The field temperature at a depth of 20 m under the ground after 30 days of continuous operation in the borehole and another 10 days of a natural recovery period is shown in Figure 6.The heat is transferred from the center of the borehole to its far ends, with the highest temperature at the center of the borehole.The higher the value of the thermal diffusion coefficient of soil, the larger its radius of thermal effect.After 30 days, the radius of thermal effect values for silt, sand, and granite are 1.68 m, 1.78 m, and 1.95 m, respectively.
The temperature at the center decreases during the natural recovery period.The temperature in the soil gradually recovers and its rate can be estimated using Equation ( 9): where is initial soil temperature, is the soil temperature under continuous operation, is the soil temperature after recovery.The recovery rate of the field temperature in silt, sand, and granite are 0.773, 0.778, and 0.789, respectively.Therefore, the higher the thermal diffusion coefficient of the soil, the faster the natural recovery of the ground temperature.The field temperature at a depth of 20 m under the ground after 30 days of continuous operation in the borehole and another 10 days of a natural recovery period is shown in Figure 6.The heat is transferred from the center of the borehole to its far ends, with the highest temperature at the center of the borehole.The higher the value of the thermal diffusion coefficient of soil, the larger its radius of thermal effect.After 30 days, the radius of thermal effect values for silt, sand, and granite are 1.68 m, 1.78 m, and 1.95 m, respectively.
The temperature at the center decreases during the natural recovery period.The temperature in the soil gradually recovers and its rate can be estimated using Equation ( 9): where T o is initial soil temperature, T is the soil temperature under continuous operation, T is the soil temperature after recovery.The recovery rate of the field temperature in silt, sand, and granite are 0.773, 0.778, and 0.789, respectively.Therefore, the higher the thermal diffusion coefficient of the soil, the faster the natural recovery of the ground temperature.The field temperature at a depth of 20 m under the ground after 30 days of continuous operation in the borehole and another 10 days of a natural recovery period is shown in Figure 6.The heat is transferred from the center of the borehole to its far ends, with the highest temperature at the center of the borehole.The higher the value of the thermal diffusion coefficient of soil, the larger its radius of thermal effect.After 30 days, the radius of thermal effect values for silt, sand, and granite are 1.68 m, 1.78 m, and 1.95 m, respectively.
The temperature at the center decreases during the natural recovery period.The temperature in the soil gradually recovers and its rate can be estimated using Equation ( 9): where is initial soil temperature, is the soil temperature under continuous operation, is the soil temperature after recovery.The recovery rate of the field temperature in silt, sand, and granite are 0.773, 0.778, and 0.789, respectively.Therefore, the higher the thermal diffusion coefficient of the soil, the faster the natural recovery of the ground temperature.

Influence of Inlet Water Temperature
Three different inlet water temperatures of 298 K, 303 K, and 308 K are selected, and the heat flux histories are shown in Figure 7; with the increase of the inlet temperature, the temperature potential to initial ground temperature increases and the heat transfer capacity of the GHE is higher.The heat flux in the GHE is the highest when inlet water temperature is 308 K, with an average value of 87.8 W/m.

Influence of Inlet Water Temperature
Three different inlet water temperatures of 298 K, 303 K, and 308 K are selected, and the heat flux histories are shown in Figure 7; with the increase of the inlet temperature, the temperature potential to initial ground temperature increases and the heat transfer capacity of the GHE is higher.The heat flux in the GHE is the highest when inlet water temperature is 308 K, with an average value of 87.8 W/m.The field temperature at a depth of 20 m under the ground is shown in Figure 8; the higher the value of the inlet water temperature, the higher the rise in the underground field temperature.When the inlet water temperature values are 298 K, 303 K, and 308 K, the borehole wall temperature values are 294 K, 297.3 K, and 301.2 K, respectively.The radius of thermal effect values are 0.26 m, 0.35 m, and 0.42 m, respectively.The heat in the underground will be accumulated if the inlet water temperature is too high, which has a negative effect on long-term operation of GHEs.

Influence of Inlet Water Velocity
Three different inlet water velocities of 0.2 m/s, 0.3 m/s and 0.6 m/s are selected and the heat flux histories are shown in Figure 9; the heat flux increases with the increase in the inlet water velocity.When the inlet water velocity is 0.6 m/s, the heat flux in the GHE is the highest, with an average value of 76 W/m.This is mainly because the heat-carrying fluid flow increases with the increase of the flow velocity.The field temperature at a depth of 20 m under the ground is shown in Figure 8; the higher the value of the inlet water temperature, the higher the rise in the underground field temperature.When the inlet water temperature values are 298 K, 303 K, and 308 K, the borehole wall temperature values are 294 K, 297.3 K, and 301.2 K, respectively.The radius of thermal effect values are 0.26 m, 0.35 m, and 0.42 m, respectively.The heat in the underground will be accumulated if the inlet water temperature is too high, which has a negative effect on long-term operation of GHEs.

Influence of Inlet Water Temperature
Three different inlet water temperatures of 298 K, 303 K, and 308 K are selected, and the heat flux histories are shown in Figure 7; with the increase of the inlet temperature, the temperature potential to initial ground temperature increases and the heat transfer capacity of the GHE is higher.The heat flux in the GHE is the highest when inlet water temperature is 308 K, with an average value of 87.8 W/m.The field temperature at a depth of 20 m under the ground is shown in Figure 8; the higher the value of the inlet water temperature, the higher the rise in the underground field temperature.When the inlet water temperature values are 298 K, 303 K, and 308 K, the borehole wall temperature values are 294 K, 297.3 K, and 301.2 K, respectively.The radius of thermal effect values are 0.26 m, 0.35 m, and 0.42 m, respectively.The heat in the underground will be accumulated if the inlet water temperature is too high, which has a negative effect on long-term operation of GHEs.

Influence of Inlet Water Velocity
Three different inlet water velocities of 0.2 m/s, 0.3 m/s and 0.6 m/s are selected and the heat flux histories are shown in Figure 9; the heat flux increases with the increase in the inlet water velocity.When the inlet water velocity is 0.6 m/s, the heat flux in the GHE is the highest, with an average value of 76 W/m.This is mainly because the heat-carrying fluid flow increases with the increase of the flow velocity.

Influence of Inlet Water Velocity
Three different inlet water velocities of 0.2 m/s, 0.3 m/s and 0.6 m/s are selected and the heat flux histories are shown in Figure 9; the heat flux increases with the increase in the inlet water velocity.When the inlet water velocity is 0.6 m/s, the heat flux in the GHE is the highest, with an average value of 76 W/m.This is mainly because the heat-carrying fluid flow increases with the increase of the flow velocity.Meanwhile, with the increase in the flow velocity, the diffusion rate of kinetic energy increases, the heat transfer in turbulent flow is strengthened, and the surface heat transfer coefficient is higher.The heat transfer coefficient, Reynolds and Nusselt numbers are calculated by the method of surface integral weighted average as: where ∅ is the variable through surface A. The results of interior flow are shown in Table 5; when the inlet water velocity values are 0.2 m/s, 0.3 m/s, and 0.6 m/s, the computed average surface heat transfer coefficients are 17.97 W/(m 2 •K), 21.06 W/(m 2 •K), and 25.95 W/(m 2 •K), respectively.

Influence of Groundwater Flow
The change in the heat flux when the Darcy velocities are set at 2 × 10 −7 m/s, 2 × 10 −6 m/s, and 2 × 10 −5 m/s is shown in Figure 10.During the initial period, the seepage in the ground can remove some of the heat and reduce the borehole wall temperature.Therefore, the heat flux in the case of seepage is higher than that in the case of no seepage.However, after the heat in the borehole is accumulated for a certain period, e.g., after 5 h, groundwater with the Darcy velocity of 2 × 10 −7 m/s can no longer reduce the borehole wall temperature effectively, and the heat transfer intensity in the GHE is close to that in the case of no seepage.When the Darcy velocity is 2 × 10 −6 m/s, the same case described above occurs after 20 h.When the Darcy velocity is increased to 2 × 10 −5 m/s, it has medium water permeability.The positive effect induced by the thermal-seepage coupling exceeds the negative effect caused by heat buildup, and the heat transfer process is improved.
As shown in Figure 11, when a Darcy velocity is 2 × 10 −7 m/s, it basically belongs to an impermeable stratum, and the temperature distribution is generally the same as that in the case of no seepage.When the Darcy velocity is 2 × 10 −6 m/s, the ability of heat diffusion is improved.When the Darcy velocity is 2 × 10 −5 m/s, the thermal-seepage coupling effect causes the contour line to shift towards the seepage direction, and the radius of the thermal effect around the borehole is significantly reduced.Meanwhile, with the increase in the flow velocity, the diffusion rate of kinetic energy increases, the heat transfer in turbulent flow is strengthened, and the surface heat transfer coefficient is higher.The heat transfer coefficient, Reynolds and Nusselt numbers are calculated by the method of surface integral weighted average as: where ∅ is the variable through surface A. The results of interior flow are shown in

Influence of Groundwater Flow
The change in the heat flux when the Darcy velocities are set at 2 × 10 −7 m/s, 2 × 10 −6 m/s, and 2 × 10 −5 m/s is shown in Figure 10.During the initial period, the seepage in the ground can remove some of the heat and reduce the borehole wall temperature.Therefore, the heat flux in the case of seepage is higher than that in the case of no seepage.However, after the heat in the borehole is accumulated for a certain period, e.g., after 5 h, groundwater with the Darcy velocity of 2 × 10 −7 m/s can no longer reduce the borehole wall temperature effectively, and the heat transfer intensity in the GHE is close to that in the case of no seepage.When the Darcy velocity is 2 × 10 −6 m/s, the same case described above occurs after 20 h.When the Darcy velocity is increased to 2 × 10 −5 m/s, it has medium water permeability.The positive effect induced by the thermal-seepage coupling exceeds the negative effect caused by heat buildup, and the heat transfer process is improved.
As shown in Figure 11, when a Darcy velocity is 2 × 10 −7 m/s, it basically belongs to an impermeable stratum, and the temperature distribution is generally the same as that in the case of no seepage.When the Darcy velocity is 2 × 10 −6 m/s, the ability of heat diffusion is improved.When the Darcy velocity is 2 × 10 −5 m/s, the thermal-seepage coupling effect causes the contour line to shift towards the seepage direction, and the radius of the thermal effect around the borehole is significantly reduced.

Conclusions
The vertical U-tube GHE is thoroughly evaluated in this study.A three-dimensional unsteady state model coupling heat transfer with groundwater seepage is developed, and the model has been validated with GHE experimental data.A large number of simulations are employed to investigate factors influencing the vertical GHE.The analysis is compared based on the heat flux in the borehole, fluid temperature along the U-tube, thermal interference between branches, turbulent flow in tubes, and temperature recovery in the ground.
The results clearly show the importance of underground soil thermal properties, grout materials, inlet water temperature and velocity, and groundwater seepage on thermal performance.The conclusions can be drawn as: (1) Before the heat in the borehole is saturated, the heat flux in the GHE is directly proportional to the thermal conductivity coefficient of the grout materials.(2) The greater the thermal diffusion coefficient of the soil materials, the larger the radius of the thermal effect of the GHE and the faster the recovery rate of the temperature in the soil.

Conclusions
The vertical U-tube GHE is thoroughly evaluated in this study.A three-dimensional unsteady state model coupling heat transfer with groundwater seepage is developed, and the model has been validated with GHE experimental data.A large number of simulations are employed to investigate factors influencing the vertical GHE.The analysis is compared based on the heat flux in the borehole, fluid temperature along the U-tube, thermal interference between branches, turbulent flow in tubes, and temperature recovery in the ground.
The results clearly show the importance of underground soil thermal properties, grout materials, inlet water temperature and velocity, and groundwater seepage on thermal performance.The conclusions can be drawn as: (1) Before the heat in the borehole is saturated, the heat flux in the GHE is directly proportional to the thermal conductivity coefficient of the grout materials.(2) The greater the thermal diffusion coefficient of the soil materials, the larger the radius of the thermal effect of the GHE and the faster the recovery rate of the temperature in the soil.

Conclusions
The vertical U-tube GHE is thoroughly evaluated in this study.A three-dimensional unsteady state model coupling heat transfer with groundwater seepage is developed, and the model has been validated with GHE experimental data.A large number of simulations are employed to investigate factors influencing the vertical GHE.The analysis is compared based on the heat flux in the borehole, fluid temperature along the U-tube, thermal interference between branches, turbulent flow in tubes, and temperature recovery in the ground.
The results clearly show the importance of underground soil thermal properties, grout materials, inlet water temperature and velocity, and groundwater seepage on thermal performance.The conclusions can be drawn as: (1) Before the heat in the borehole is saturated, the heat flux in the GHE is directly proportional to the thermal conductivity coefficient of the grout materials.(2) The greater the thermal diffusion coefficient of the soil materials, the larger the radius of the thermal effect of the GHE and the faster the recovery rate of the temperature in the soil.(3) Increasing the inlet water temperature strengthens the heat transfer in the GHE, but it can cause heat buildup problems.(4) Increasing the inlet water velocity can enlarge the heat convection in the tube; when the inlet water velocity values are 0.2 m/s, 0.3 m/s, and 0.6 m/s, the computed average surface heat transfer coefficients are 17.97 W/(m 2 •K), 21.06 W/(m 2 •K), and 25.95 W/(m 2 •K), respectively.(5) The thermal-seepage coupling effect in groundwater can remove the accumulated heat, and therefore, can effectively enhance the heat transfer in the GHE.

Figure 2 .
Figure 2. The finite elements of GHE.
homogenous and isotropic; (ii) the contact thermal resistances of joints are negligible;

Figure 2 .
Figure 2. The finite elements of GHE.

Figure 2 .
Figure 2. The finite elements of GHE.

Figure 3 .Figure 4 .
Figure 3. Heat flux variation for the different grout materials.

Figure 4 .
Figure 4. Temperature in both inlet and outlet tubes for the different grout materials.

Figure 5 .
Figure 5. Heat flux variation for the different soils.

Figure 6 .
Figure 6.Temperature distributions under the ground for the different soils.(a) After a continuous operation of 30 days; (b) after a natural recovery of 10 days.

Figure 5 .
Figure 5. Heat flux variation for the different soils.

Figure 5 .
Figure 5. Heat flux variation for the different soils.

Figure 6 .
Figure 6.Temperature distributions under the ground for the different soils.(a) After a continuous operation of 30 days; (b) after a natural recovery of 10 days.

Figure 6 .
Figure 6.Temperature distributions under the ground for the different soils.(a) After a continuous operation of 30 days; (b) after a natural recovery of 10 days.

Figure 7 .
Figure 7. Heat flux variation for the different inlet temperatures.

Figure 8 .
Figure 8. Temperature distributions under the ground for different inlet temperatures.

Figure 7 .
Figure 7. Heat flux variation for the different inlet temperatures.

Figure 7 .
Figure 7. Heat flux variation for the different inlet temperatures.

Figure 8 .
Figure 8. Temperature distributions under the ground for different inlet temperatures.

Figure 8 .
Figure 8. Temperature distributions under the ground for different inlet temperatures.

Figure 9 .
Figure 9. Heat flux variation for the different inlet velocities.

Figure 9 .
Figure 9. Heat flux variation for the different inlet velocities.

Figure 10 .
Figure 10.Heat flux variation for the different groundwater seepage values.

Figure 11 .
Figure 11.Temperature distributions under the ground for different groundwater seepage values.

( 3 )
Increasing the inlet water temperature strengthens the heat transfer in the GHE, but it can cause heat buildup problems.(4) Increasing the inlet water velocity can enlarge the heat convection in the tube; when the inlet water velocity values are 0.2 m/s, 0.3 m/s, and 0.6 m/s, the computed average surface heat transfer coefficients are 17.97 W/(m 2 •K), 21.06 W/(m 2 •K), and 25.95 W/(m 2 •K), respectively.(5) The thermal-seepage coupling effect in groundwater can remove the accumulated heat, and therefore, can effectively enhance the heat transfer in the GHE.

Figure 10 .Figure 10 .
Figure 10.Heat flux variation for the different groundwater seepage values.

Figure 11 .
Figure 11.Temperature distributions under the ground for different groundwater seepage values.

( 3 )
Increasing the inlet water temperature strengthens the heat transfer in the GHE, but it can cause heat buildup problems.(4) Increasing the inlet water velocity can enlarge the heat convection in the tube; when the inlet water velocity values are 0.2 m/s, 0.3 m/s, and 0.6 m/s, the computed average surface heat transfer coefficients are 17.97 W/(m 2 •K), 21.06 W/(m 2 •K), and 25.95 W/(m 2 •K), respectively.(5) The thermal-seepage coupling effect in groundwater can remove the accumulated heat, and therefore, can effectively enhance the heat transfer in the GHE.

Figure 11 .
Figure 11.Temperature distributions under the ground for different groundwater seepage values.

Table 2 .
The reference parameters of GHE for simulation.

Table 3 .
The thermophysical properties of the grout materials.

Table 4 .
The thermophysical properties of the soils.

Table 4 .
The thermophysical properties of the soils.

Table 4 .
The thermophysical properties of the soils.

Table 5 .
The results of interior flow.

Table 5 .
The results of interior flow.