Home Effect of buried depth on thermal performance of a vertical U-tube underground heat exchanger
Article Open Access

Effect of buried depth on thermal performance of a vertical U-tube underground heat exchanger

  • Li Yang , Bo Zhang , Jiří Jaromír Klemeš , Jie Liu , Meiyu Song and Jin Wang EMAIL logo
Published/Copyright: June 8, 2021

Abstract

Many researchers numerically investigated U-tube underground heat exchanger using a two-dimensional simplified pipe. However, a simplified model results in large errors compared to the data from construction sites. This research is carried out using a three-dimensional full-size model. A model validation is conducted by comparing with experimental data in summer. This article investigates the effects of fluid velocity and buried depth on the heat exchange rate in a vertical U-tube underground heat exchanger based on fluid–structure coupled simulations. Compared with the results at a flow rate of 0.4 m/s, the results of this research show that the heat transfer per buried depth at 1.0 m/s increases by 123.34%. With the increase of the buried depth from 80 to 140 m, the heat transfer per unit depth decreases by 9.72%.

1 Introduction

Considering global energy emission issues of fossil fuels, renewable energy sources (such as solar and geothermal energies) have been focused on by many researchers. Various phase change materials (PCMs) have been widely used in renewable energy applications [1]. Li et al. [2] proposed a novel solar thermal system coupling with active PCM heat storage wall, which met a time balance between the solar supply and demand. A ground source heat pump (GSHP), like vertical U-tube GSHP, has been widely used for the exploitation of geothermal energy due to its advantages, such as energy-saving and eco-friendly utilizations. A vertical U-tube underground heat exchanger has a high initial investment, which is needed to put forward an accurate heat transfer model for a reduction in initial investment of energy systems [3].

Based on an investigation of fluid–structure coupled simulations, a three-dimensional heat transfer model about the vertical U-tube GSHP will be proposed in this article; and this model will use independent programming development of user-defined functions in the software ANSYS FLUENT 18.0 [4]. To discuss the effects of fluid velocity and buried depth on heat transfer performance of a vertical U-tube ground heat exchanger (GHE), inlet velocities are set to 0.4, 0.6, 0.8, and 1.0 m/s under the same buried depth of 120 m. With the same inlet velocity of 0.6 m/s, the thermal performance of the vertical U-tube GHE is also investigated under different buried depths of 80, 100, 120, and 140 m.

2 Methodologies

The test system consists of a vertical U-tube GHE, ground-coupled heat pump unit, circulation pump, and data acquisition system with measuring devices (like thermocouples, thermometers, and flow meters). These measuring devices mainly collect temperature and flow rates of supply water and return water of the GHE. A diagram is represented in Figure 1.

Figure 1 
               Experimental system. (a) Real picture of the experimental facility. (b) Diagram of the test system.
Figure 1

Experimental system. (a) Real picture of the experimental facility. (b) Diagram of the test system.

It is a complex heat transfer process between the vertical U-tube GHE and the surrounding soil, and some assumptions are made to simplify the model.

For one assumption that the topsoil is not affected by solar radiation, the whole soil has an initial temperature as denoted by equation (1) [5]. Thermal resistance of all contact surfaces is given to be 0.

(1) T = T s + T d Z ,

where T is the temperature at a given depth (K). T s is reference temperature at 0 m, and its value is 282.15 K. T d is temperature gradient of 0.03 K/m, and Z is buried depth (m).

A three-dimensional model has been established with a scale of 1:1. Figure 2 shows the computational domain of the vertical U-tube GHE.

Figure 2 
               Computational domain with geometry dimensions.
Figure 2

Computational domain with geometry dimensions.

The external and internal diameters of the U-pipe are 32 and 26.2 mm, respectively. It is a distance of 100 mm between two branches of the U-pipe. The buried depth of the U-pipe is 120 m, and the soil depth is 121 m. The diameters of the borehole and the surrounding soil are 250 mm and 10 m, respectively.

The computational domain consists of PE pipe, working fluid, backfill, and surrounding soil. Table 1 shows the thermophysical properties of materials, as provided by Tianjin Geothermal Exploration and Development-Design Institute.

Table 1

Thermophysical properties of materials

Term Density (kg/m3) Specific heat (J/kg K) Thermal conductivity (W/m K)
PE pipe 940 2,300 0.65
Water 998.2 4,182 0.6
Backfill material 1,860 840 3
Soil 1,990 1,380 1.7

Initial backfill, PE pipe, and soil are supposed to start at temperature values from equation (1). Based on the boundary conditions during the tests, the inlet velocity is set to 0.6 m/s, and the inlet temperature changes with time, as shown in equation (2), is as follows:

(2) T i = 8.62 × e t 26222.79 4.24 × e t 1008.11 + 303.06 ,

where T i is inlet temperature (K), and t is time (s).

The equation (2) is derived from the experimental data provided by the Tianjin Geothermal Exploration and Development-Design Institute, China. All contact surfaces are set as temperature-coupled walls, and the upper surfaces are set as adiabatic surfaces. Water is employed as a heat carrier in GHEs.

The heat exchange rate was used to analyze the heat transfer performance of vertical U-tube GHEs. An equation of heat exchange rate, as given in Javadi et al. [6], is as follows:

(3) Q = m C p ( T i T o ) ,

where Q is the heat exchange rate (W), m is the mass flow rate (kg/s), C p is the specific heat capacity (J/kg K), T i is the water temperature at the inlet (K), and T o is the water temperature at the outlet (K). The heat exchange rate per meter of borehole depth (Q L) is calculated as follows:

(4) Q L = Q L ,

where L is the depth of a GHE (m).

As shown in Figure 3, model validation and grid independence tests are conducted using cases with four grid meshes from 3 to 24 M. It is found that simulated results are consistent with experimental data tested in summer.

Figure 3 
               Model validation (M for million).
Figure 3

Model validation (M for million).

Figure 3 shows the temperature differences between experimental data and simulated results at different grid numbers. Outlet experimental temperature is 298.42 K. It is found that when the grid number is above 12 M, the outlet temperature is almost constant. Compared to the tested data, temperature variations from the cases with 12 and 24 M grid cells are 3.24 and 3.23%, respectively. The results from 12 M cell case are acceptable and finally used for all the simulations in this research.

3 Results and discussion

The heat transfer performance of a vertical U-tube GHE is investigated under eight operation conditions. Figure 4 shows heat exchange rates with different inlet velocities. It is found that heat exchange rates increase with the increase of the inlet velocity. It is also found that the Q is 3.42 kW at an inlet velocity of 0.4 m/s. An increase of 123.3% is observed when the velocity changes from 0.4 to 1.0 m/s.

Figure 4 
               Heat exchange rates with different inlet velocities.
Figure 4

Heat exchange rates with different inlet velocities.

Heat exchange rates increase with the increase of buried depth, as shown in Figure 5. The heat exchange rate increases by 57% with the increase of the buried depth from 80 to 140 m. Heat exchange rate per meter of borehole depth is directly proportional to the inlet velocity and inversely proportional to the buried depth. Comparing with 28.47 W/m at inlet velocity of 0.4 m/s, Q L increases by 123.34% at 1.0 m/s. Q L decreases by 9.72% with a buried depth increasing from 80 to 140 m. The growth rates of Q are 20.5, 18.4, and 18.1%, when the buried depth changes from 80 to 100 m, 100 to 120 m, and 120 to 140 m, respectively. The growth rate of Q tends to be constant when the buried depth is above 120 m. It is concluded that 120 m is the optimal buried depth, considering thermal performance.

Figure 5 
               Heat exchange rates with different buried depths.
Figure 5

Heat exchange rates with different buried depths.

The construction costs consist of expenses of drilling and grouting. Supplied heat produced by GHEs can be equivalent to electricity power based on ladder electricity price in China. Table 2 shows costs of GHE construction and 70-day electricity saving. Although the construction costs increase with the increase of the borehole depth, the electricity saving using GHEs increases. The recommended borehole depth depends on service time every year and cost recovery time as planned.

Table 2

Costs of GHE construction and electricity saving

Buried depth (m) GHE length (m) Construction cost (USD) 70-day saving bill (USD)
80 160.2 833.9 695.2
100 200.2 1042.1 841.0
120 240.2 1250.3 971.9
140 280.2 1458.5 1100.6

4 Conclusion

In this research, the heat transfer performance of a vertical U-tube GHE was studied numerically. Q increases with the increases of both the inlet velocity and the buried depth. Q L is proportional to inlet velocity and inversely proportional to buried depth. The buried depth of 120 m is optimal in this study. The heat diffusion of turbulent is considered as the main factor affecting the heat transfer performance of the vertical U-tube heat exchanger.


tel: +86-22-60435005

Acknowledgements

This work is supported by the Project of Innovation Ability Training for Postgraduate Students of Education Department of Hebei Province [Grant number CXZZSS2021046] and the Foundation of Key Laboratory of Thermo-Fluid Science and Engineering (Xi’an Jiaotong University), Ministry of Education, Xi’an 710049, P. R. China (Grant No. KLTFSE2018KFJJ01).

  1. Funding information: The authors acknowledge the funding and support by the EU project Sustainable Process Integration Laboratory – SPIL, funded as project No. CZ.02.1.01/0.0/0.0/15_003/0000456, by Czech Republic Operational Programme Research and Development, Education, Priority 1: Strengthening capacity for quality research operated by the Czech Ministry of Education, Youth and Sport under a collaboration agreement with Hebei University of Technology.

  2. Conflict of interest: Authors state no conflict of interest.

References

[1] Yu K, Wang Y, Li Y, Baleta J, Wang J, Sundén B. Effect of phase change materials on heat dissipation of a multiple heat source system. Open Phys. 2019;17:797–807.10.1515/phys-2019-0083Search in Google Scholar

[2] Li H, Li J, Kong X, Long H, Yang H, Yao C. A novel solar thermal system combining with active phase-change material heat storage wall (STS-APHSW): dynamic model, validation and thermal performance. Energy. 2020;201:117610.10.1016/j.energy.2020.117610Search in Google Scholar

[3] Klemeš JJ, Arsenyeva O, Kapustenko P, Tovazhnyanskyy L. Compact heat exchangers for energy transfer intensification. 1st edn. New York: CRC Press; 2015.10.1201/b18862Search in Google Scholar

[4] Khalajzadeh V, Heidarinejad G, Srebric J. Parameters optimization of a vertical ground heat exchanger based on response surface methodology. Energy Build. 2011;43:1288–94.10.1016/j.enbuild.2011.01.007Search in Google Scholar

[5] Samin MY, Faramarzi A, Jefferson I, Harireche O. A hybrid optimisation approach to improve long-term performance of enhanced geothermal system (EGS) reservoirs. Renew Energy. 2019;134:379–89.10.1016/j.renene.2018.11.045Search in Google Scholar

[6] Javadi H, Mousavi Ajarostaghi SS, Pourfallah M, Zaboli M. Performance analysis of helical ground heat exchangers with different configurations. Appl Therm Eng. 2019;154:24–36.10.1016/j.applthermaleng.2019.03.021Search in Google Scholar

Received: 2021-02-01
Revised: 2021-04-15
Accepted: 2021-04-21
Published Online: 2021-06-08

© 2021 Li Yang et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Circular Rydberg states of helium atoms or helium-like ions in a high-frequency laser field
  3. Closed-form solutions and conservation laws of a generalized Hirota–Satsuma coupled KdV system of fluid mechanics
  4. W-Chirped optical solitons and modulation instability analysis of Chen–Lee–Liu equation in optical monomode fibres
  5. The problem of a hydrogen atom in a cavity: Oscillator representation solution versus analytic solution
  6. An analytical model for the Maxwell radiation field in an axially symmetric galaxy
  7. Utilization of updated version of heat flux model for the radiative flow of a non-Newtonian material under Joule heating: OHAM application
  8. Verification of the accommodative responses in viewing an on-axis analog reflection hologram
  9. Irreversibility as thermodynamic time
  10. A self-adaptive prescription dose optimization algorithm for radiotherapy
  11. Algebraic computational methods for solving three nonlinear vital models fractional in mathematical physics
  12. The diffusion mechanism of the application of intelligent manufacturing in SMEs model based on cellular automata
  13. Numerical analysis of free convection from a spinning cone with variable wall temperature and pressure work effect using MD-BSQLM
  14. Numerical simulation of hydrodynamic oscillation of side-by-side double-floating-system with a narrow gap in waves
  15. Closed-form solutions for the Schrödinger wave equation with non-solvable potentials: A perturbation approach
  16. Study of dynamic pressure on the packer for deep-water perforation
  17. Ultrafast dephasing in hydrogen-bonded pyridine–water mixtures
  18. Crystallization law of karst water in tunnel drainage system based on DBL theory
  19. Position-dependent finite symmetric mass harmonic like oscillator: Classical and quantum mechanical study
  20. Application of Fibonacci heap to fast marching method
  21. An analytical investigation of the mixed convective Casson fluid flow past a yawed cylinder with heat transfer analysis
  22. Considering the effect of optical attenuation on photon-enhanced thermionic emission converter of the practical structure
  23. Fractal calculation method of friction parameters: Surface morphology and load of galvanized sheet
  24. Charge identification of fragments with the emulsion spectrometer of the FOOT experiment
  25. Quantization of fractional harmonic oscillator using creation and annihilation operators
  26. Scaling law for velocity of domino toppling motion in curved paths
  27. Frequency synchronization detection method based on adaptive frequency standard tracking
  28. Application of common reflection surface (CRS) to velocity variation with azimuth (VVAz) inversion of the relatively narrow azimuth 3D seismic land data
  29. Study on the adaptability of binary flooding in a certain oil field
  30. CompVision: An open-source five-compartmental software for biokinetic simulations
  31. An electrically switchable wideband metamaterial absorber based on graphene at P band
  32. Effect of annealing temperature on the interface state density of n-ZnO nanorod/p-Si heterojunction diodes
  33. A facile fabrication of superhydrophobic and superoleophilic adsorption material 5A zeolite for oil–water separation with potential use in floating oil
  34. Shannon entropy for Feinberg–Horodecki equation and thermal properties of improved Wei potential model
  35. Hopf bifurcation analysis for liquid-filled Gyrostat chaotic system and design of a novel technique to control slosh in spacecrafts
  36. Optical properties of two-dimensional two-electron quantum dot in parabolic confinement
  37. Optical solitons via the collective variable method for the classical and perturbed Chen–Lee–Liu equations
  38. Stratified heat transfer of magneto-tangent hyperbolic bio-nanofluid flow with gyrotactic microorganisms: Keller-Box solution technique
  39. Analysis of the structure and properties of triangular composite light-screen targets
  40. Magnetic charged particles of optical spherical antiferromagnetic model with fractional system
  41. Study on acoustic radiation response characteristics of sound barriers
  42. The tribological properties of single-layer hybrid PTFE/Nomex fabric/phenolic resin composites underwater
  43. Research on maintenance spare parts requirement prediction based on LSTM recurrent neural network
  44. Quantum computing simulation of the hydrogen molecular ground-state energies with limited resources
  45. A DFT study on the molecular properties of synthetic ester under the electric field
  46. Construction of abundant novel analytical solutions of the space–time fractional nonlinear generalized equal width model via Riemann–Liouville derivative with application of mathematical methods
  47. Some common and dynamic properties of logarithmic Pareto distribution with applications
  48. Soliton structures in optical fiber communications with Kundu–Mukherjee–Naskar model
  49. Fractional modeling of COVID-19 epidemic model with harmonic mean type incidence rate
  50. Liquid metal-based metamaterial with high-temperature sensitivity: Design and computational study
  51. Biosynthesis and characterization of Saudi propolis-mediated silver nanoparticles and their biological properties
  52. New trigonometric B-spline approximation for numerical investigation of the regularized long-wave equation
  53. Modal characteristics of harmonic gear transmission flexspline based on orthogonal design method
  54. Revisiting the Reynolds-averaged Navier–Stokes equations
  55. Time-periodic pulse electroosmotic flow of Jeffreys fluids through a microannulus
  56. Exact wave solutions of the nonlinear Rosenau equation using an analytical method
  57. Computational examination of Jeffrey nanofluid through a stretchable surface employing Tiwari and Das model
  58. Numerical analysis of a single-mode microring resonator on a YAG-on-insulator
  59. Review Articles
  60. Double-layer coating using MHD flow of third-grade fluid with Hall current and heat source/sink
  61. Analysis of aeromagnetic filtering techniques in locating the primary target in sedimentary terrain: A review
  62. Rapid Communications
  63. Nonlinear fitting of multi-compartmental data using Hooke and Jeeves direct search method
  64. Effect of buried depth on thermal performance of a vertical U-tube underground heat exchanger
  65. Knocking characteristics of a high pressure direct injection natural gas engine operating in stratified combustion mode
  66. What dominates heat transfer performance of a double-pipe heat exchanger
  67. Special Issue on Future challenges of advanced computational modeling on nonlinear physical phenomena - Part II
  68. Lump, lump-one stripe, multiwave and breather solutions for the Hunter–Saxton equation
  69. New quantum integral inequalities for some new classes of generalized ψ-convex functions and their scope in physical systems
  70. Computational fluid dynamic simulations and heat transfer characteristic comparisons of various arc-baffled channels
  71. Gaussian radial basis functions method for linear and nonlinear convection–diffusion models in physical phenomena
  72. Investigation of interactional phenomena and multi wave solutions of the quantum hydrodynamic Zakharov–Kuznetsov model
  73. On the optical solutions to nonlinear Schrödinger equation with second-order spatiotemporal dispersion
  74. Analysis of couple stress fluid flow with variable viscosity using two homotopy-based methods
  75. Quantum estimates in two variable forms for Simpson-type inequalities considering generalized Ψ-convex functions with applications
  76. Series solution to fractional contact problem using Caputo’s derivative
  77. Solitary wave solutions of the ionic currents along microtubule dynamical equations via analytical mathematical method
  78. Thermo-viscoelastic orthotropic constraint cylindrical cavity with variable thermal properties heated by laser pulse via the MGT thermoelasticity model
  79. Theoretical and experimental clues to a flux of Doppler transformation energies during processes with energy conservation
  80. On solitons: Propagation of shallow water waves for the fifth-order KdV hierarchy integrable equation
  81. Special Issue on Transport phenomena and thermal analysis in micro/nano-scale structure surfaces - Part II
  82. Numerical study on heat transfer and flow characteristics of nanofluids in a circular tube with trapezoid ribs
  83. Experimental and numerical study of heat transfer and flow characteristics with different placement of the multi-deck display cabinet in supermarket
  84. Thermal-hydraulic performance prediction of two new heat exchangers using RBF based on different DOE
  85. Diesel engine waste heat recovery system comprehensive optimization based on system and heat exchanger simulation
  86. Load forecasting of refrigerated display cabinet based on CEEMD–IPSO–LSTM combined model
  87. Investigation on subcooled flow boiling heat transfer characteristics in ICE-like conditions
  88. Research on materials of solar selective absorption coating based on the first principle
  89. Experimental study on enhancement characteristics of steam/nitrogen condensation inside horizontal multi-start helical channels
  90. Special Issue on Novel Numerical and Analytical Techniques for Fractional Nonlinear Schrodinger Type - Part I
  91. Numerical exploration of thin film flow of MHD pseudo-plastic fluid in fractional space: Utilization of fractional calculus approach
  92. A Haar wavelet-based scheme for finding the control parameter in nonlinear inverse heat conduction equation
  93. Stable novel and accurate solitary wave solutions of an integrable equation: Qiao model
  94. Novel soliton solutions to the Atangana–Baleanu fractional system of equations for the ISALWs
  95. On the oscillation of nonlinear delay differential equations and their applications
  96. Abundant stable novel solutions of fractional-order epidemic model along with saturated treatment and disease transmission
  97. Fully Legendre spectral collocation technique for stochastic heat equations
  98. Special Issue on 5th International Conference on Mechanics, Mathematics and Applied Physics (2021)
  99. Residual service life of erbium-modified AM50 magnesium alloy under corrosion and stress environment
  100. Special Issue on Advanced Topics on the Modelling and Assessment of Complicated Physical Phenomena - Part I
  101. Diverse wave propagation in shallow water waves with the Kadomtsev–Petviashvili–Benjamin–Bona–Mahony and Benney–Luke integrable models
  102. Intensification of thermal stratification on dissipative chemically heating fluid with cross-diffusion and magnetic field over a wedge
Downloaded on 13.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2021-0033/html
Scroll to top button