Study on a Dynamic Numerical Model of an Underground Air Tunnel System for Cooling Applications—Experimental Validation and Multidimensional Parametrical Analysis

: The underground air tunnel system shows promising potentials for reducing energy consumption of buildings and for improving indoor thermal comfort, whereas the existing dynamic models using the computational fluid dynamic (CFD) method show computational complexity and are user ‐ unfriendly for parametrical analysis. In this study, a dynamic numerical model was developed with the on ‐ site experimental calibration. Compared to the traditional CFD method with high computational complexity, the mathematical model on the MATLAB/SIMULINK platform is time ‐ saving in terms of the real ‐ time thermal performance prediction. The experimental validation results indicated that the maximum absolute relative deviation was 3.18% between the model ‐ driven results and the data from the on ‐ site experiments. Parametrical analysis results indicated that, with the increase of the tube length, the outlet temperature decreases with an increase of the cooling capacity whereas the increasing/decreasing magnitude slows down. In addition, the system performance is independent on the tube materials. Furthermore, the outlet air temperature and cooling capacity are dependent on the tube diameter and air velocity, i.e., a larger tube diameter and a higher air velocity are more suitable to improve the system’s cooling capacity, and a smaller tube diameter and a lower air velocity will produce a more stable and lower outlet temperature. Further studies need to be conducted for the trade ‐ off solutions between air velocity and tube diameter for the bi ‐ criteria performance enhancement between outlet temperature and cooling capacity. This study proposed an experimentally validated mathematical model to accurately predict the thermal performance of the underground air tunnel system with high computational efficiency, which can provide technical guidance to multi ‐ combined solutions through geometrical designs and operating parameters for the optimal design and robust operation.


Introduction
As a considerable energy consumer, building energy consumption has become one of the three major energy users all around the world, accounting for more than 30% of the total global energy consumptions [1][2][3]. In China, the energy consumption in buildings increases every year, especially with social development and the improvement of living conditions. The considerable building energy consumptions require high dependence on fossil fuels, which will result in deteriorated environmental concerns [4][5]. The contradiction between the daily increased energy demand and the energy shortage calls for the necessity to deploy renewable energy [6][7][8][9] and to reduce the dependence on fossil fuels. In order to reduce the reliance on traditional fossil fuels without sacrificing economic and social development, the deployment of large-scale renewable energy systems (such as heating, cooling, and electric systems) has become one of the most promising solutions to reduce the building energy consumption and to cover the building demand [10,11]. Depending on the renewable energy sources, the renewable energy systems can be classified into solar-driven [12][13][14][15][16][17], wind-driven, and geothermal-driven energy systems. Compared to the instability in temporal distribution and the unevenness in spatial distribution of solar and wind resources, the geothermal energy shows promising potentials in terms of the abundant resource distribution, the stability of energy supply, the simplicity of system designs and operations, and so on. Geothermal energy has attracted increasing interests with prioritized recommendations for regulating indoor and built environments and for reducing building energy consumptions [18][19][20]. The underlying mechanism for the application the geothermal energy in buildings is that, due to the high thermal inertia of underground soil, independent of the ambient air temperature, solar radiation, and water flows, the soil temperature (3-4 m beneath the ground) normally remains relatively constant during the whole year. The relatively steady soil temperature is promising when being utilized as a "heat sink" or "heating source" for space cooling in summer or as a space heater in winter [21,22].
The underground air tunnel (UAT) system, as one of the most widely used geothermal energy technologies, has attracted more and more interests for reducing energy consumption of buildings [23,24], especially in areas where geothermal energy can be highly recommended. The UAT system is also called the earth-to-air heat exchanger (EAHE), the air-soil heat exchanger (ASHE), or the earth-air tunnel (EAT). A UAT usually consist of a mechanical fan and network buried tubes (commonly made by polyethylene, (PVC), or steel) at a given soil depth. These buried tubes are designed horizontally into underground soil in different geometrical forms, such as ring, serpentine, or parallel tubes. The working principle is that the ambient air is pumped into the buried tubes by the fan through mechanical ventilation. Through the heat exchange with the surrounding soil, the fresh air is cooled by the surrounding soil for space cooling in summer and the outside fresh air is heated for space heating in winter [25]. Depending on the heat-transfer rate and energy quantity, the heat-transfer efficiency is dependent on the input air temperature, air flow rate, tube materials, geometrical dimensions of tubes, buried depth, soil temperature, and soil thermal properties. According to the research results in the latest studies, UAT systems are characterized with the building energy-saving potential and the improvement of indoor thermal comfort [26,27].
To quantitatively investigate the impact of geometrical and operating parameters on the system performance, researchers have mainly focused on the development of mathematical models to accurately predict the thermal and energy performances [28,29]. The commonly used mathematical models can be classified as the steady-state model and dynamic models based on the temperature variation of surrounding soil outside buried tubes with operation time [30]. The multivariables include the tube parameters (i.e., the tube thickness, radius, depth, length, and coefficient of heat conductivity), soil parameters (i.e., thermal conductivity and moisture content), and operating parameters (i.e., air flow rate). The investigated objectives include the outlet temperature and the cooling capacity of UAT systems [31]. With respect to the steady-state model of the UAT system as already developed in academia, the main drawback is ignorance of the variation of soil temperature for complex design and operation parameters [32]. In order to investigate the effects of multivariables on the UAT system performance, Niu et al. [32] developed a one-dimensional steady-state model considering both sensible and latent heat transfer for analyzing the cooling capacity of a UAT system. Results showed that the system's cooling capacity was highly dependent on structural factors, e.g., tube length, diameter, and inlet temperature. In hot-arid and cold climates in Iran, Fazlikhani et al. [33] developed a steady-state model using the MATLAB software to study the influences of various parameters on the cooling/heating capacity of a UAT system. The investigated multivariables included the tube lengths, inlet air, and soil temperatures. Results illustrated that, in both climates, the UAT system can be highly used for preheating in winter and for cooling in summer. As a matter of fact, although the dynamic soil temperature evolution is subtle due to the considerable thermal inertia of the surrounding soil, the ignorance of the temperature evolution will result in the prediction inaccuracy. According to the research results from the cooling application scenario [34], when the soil was located at 0.1 m from the tube wall along the radial direction, the temperature increased from 26.39 °C to 27.27 °C after 5-hour continuous operation. The impact of the soil-temperature rise on the cooling performance of the UAT system is worthy to be investigated, whereas the current literature provides limited progress for the development of mathematical models with the temperature variation of dynamic surrounding soil.
During real operation of the UAT system, temperatures of heat exchange mediums (e.g., flowing air, tube wall, and surrounding soil) are dynamically and constantly changing, depending on the heat-transfer rate and the thermal capacity of the material. The dynamic temperature of each medium can be characterized by the harmonic functions of time, which is one of the effective solutions to complement ignorance of the variation of soil temperature. With the implementation of harmonic function, the outlet temperature can also be considered as a time-harmonic function [35]. Nowadays, the development of dynamic models to accurately predict the dynamic performance of the UAT system has attracted widespread attention. It can be noticed that most dynamic models were based on computational fluid dynamics (CFD) [24,36]. Misra et al. [37] studied the performance of the UAT system using a CFD-based dynamic model. Systematic and parametrical analysis results indicated that the dynamic performance from the thermal transient response was highly dependent on the soil coefficient of heat conductivity and the time-duration for continuous operation. The comparative analysis between the dynamic model and the steady-state model indicated that, compared to the outlet temperature at 18.8 °C in the steady-state model, the predicted outlet temperature was much lower between 16.6 °C and 18.7 °C through the dynamic model. Bhutta et al. [38] presented a comprehensive review on CFD applications in the geometrical design of various heat exchangers. Research results indicated that the CFD was an effective tool for the system design with robust operational performance.
However, with respect to the comprehensive and systematic parametric analysis, the CFD method requires the pretreatment software to structuralize geometric parameters and the grid partition for the development of numerical model. New CFD-based models and mesh generations are required to be conducted to further analyze the impacts of design parameters (e.g., tube diameters, tube thicknesses, and tube lengths) on the system's thermal performance [37]. In addition, the CFD method has been certificated to require higher computing capacities, especially for large-scale systems [39]. The computational complexity during the remodelling process and the simulation process requires the necessity to develop new methodologies.
This study is to cover the above-mentioned scientific gaps. In this paper, a dynamic model of a new UAT system was developed on the MATLAB/SIMULINK platform for the first time. Compared to the most common dynamic simulation method on CFD tools, such as Fluent, the MATLAB/SIMULINK platform is more flexible for multi-domain simulation, together with lower computational complexity for model-based design [34,40]. The MATLAB/SIMULINK platform has been utilized for heating ventilation and air conditioning (HVAC) system operation [41], building envelopes design [42], and hybrid solar heating systems [43]. The implementation of the MATLAB/SIMULINK-based dynamic mathematical model can improve the flexibility for the multidimensional parametrical analysis.
The main contributions of this work include the following:  A dynamic model of the UAT system was developed for the first time through the MATLAB/SIMULINK platform.  Experimental platform was established to calibrate the accuracy of the developed mathematical model.  Based on the developed mathematical model on the MATLAB/SIMULINK platform, multidimensional parametrical analysis has been conducted with high computational efficiency for technical guidance.
The research framework of this paper is as follows: in Section 2, the mathematic model of the UAT system is described. Section 3 shows the validation of the developed model through the experimental results. Multidimensional parametrical analysis was conducted in Section 4 to quantitatively demonstrate the effects of parameters on the outlet temperature and the cooling capacity. Afterwards, conclusions and future studies are drawn in Section 5.

Research Method
In order to investigate the multivariant impacts on the outlet temperature and cooling capacity of the UAT system, parametrical analysis has been conducted on inlet air temperatures, air velocities, and design parameters (e.g., tube depth, diameter, length, material, and wall thickness) using an experimentally calibrated mathematical model. In the mathematical model of the UAT system, the heat conduction was modelled to dynamically calculate the heat transfer within the underground soil. Convective heat-transfer processes along the air flow direction were modelled for the heat exchange between two adjacent air layers and between the flowing air and tube. Multidimensional matrices were developed for different heat exchange mediums, following the principle of energy balance. The multidimensional matrices were programmed in MATLAB language, and then, a block was built for each parameter calculation of thermal performance through the SIMULINK platform using the S-Function, as shown in Figure 1. The S-Function block can expand the function of the SIMULINK platform by computing a subroutine using various programming languages (e.g., MATLAB, C, C++, or Fortran), which can request the input variables for the generation of output results. The developed block can be used to dynamically predict the air, tube, and soil temperature during the continuous operation process.

Model Assumptions
To simplify the dynamic performance prediction of the UAT system without sacrificing the prediction accuracy, some assumptions have been made, as follows: 1. The soil is assumed to be homogeneous, and the thermophysical properties are considered isotropic 2. Air is incompressible with constant density and specific heat 3. The effect of air humidity variation on the heat-transfer process is ignored 4. The heat exchange along the vertical buried tube from the inlet to the buried depth is ignored 5. The thermal interference between two adjacent tubes is ignored 6. The heat conduction between the adjacent soil along the air flowing direction is ignored 7. The effect of convection and radiation on the ground surface on the heat-transfer process of the UAT system is ignored

Mathematical Equations
In this system, the model of air, tube wall, and soil along the flowing air direction can be split into n volume elements and the distance between two-adjacent divided layers is denoted as Δx, as shown in Figure 2. The temperature of each layer is assumed to be the same as the discretized node. Each divided element follows the energy balance law. The energy balance equation of each volume element can be solved at each time step to iteratively calculate the temperature variation. In the mathematical model, the temperatures of each node at time t and the previous time are denoted as T t and T t-dt . The linear mean temperature T = (T t + T t-dt )/2 was adopted for the temperature between time t-dt and t [44]. Based on this, the nodal discretization of the air, tube, and soil located at the j th layer can be illustrated in Figure 3.  (1) Air model: From Equation (1): where ρ a and C a are the air density and the specific heat capacity of air; V a is the divided air volume, V a = πr 2 × Δx, where is the tube inner radius and is the divided length; T a,j , T a,j-1 , and T t,j are the mean air temperatures of the jth and (j-1)th layers and of the tube wall between times t and t-dt, respectively; , is the air temperature at time t-dt; is the air flow rate, m = πr 2 × v × ρ a , where is the air velocity; and R a,j is the heat-transfer resistance between flowing air and tube wall: where ℎ is the convective heat-transfer coefficient, h = 0.023Re 0.8 Pr n λ a 2 r , and n = 0.3 for cooling [32].
is the area of the tube wall, S a = 2 × π × r × . is the coefficient of heat conductivity of the tube wall. is the thickness of the tube wall.

Re is the Reynolds number
is the air coefficient of heat conductivity. Re and Pr are Reynolds number and Prandtl. is the air dynamic viscosity.
(2) Tube wall model From Equation (4), the temperature of node j at time t-dt is shown below: where , , and are the tube density, specific heat capacity, and divided volume, respectively, and where V t = π((r + d) 2 is the tube wall temperature at time t-dt.
, , is the mean temperature of the 1st column soil between time t and t-dt. , , is the heat-transfer resistance between the tube wall and the 1st column soil.
where , is the heat conductivity of the soil. ∆ is the distance of the neighboring divided soil block radius.
(3) Soil model 1) For soil-1 model: From Equation (7), the temperature of node j at time t-dt is shown below: where and are the soil density and specific heat capacity, respectively.
, , is the 1st column soil temperature at time t-dt.
, , and , , are the mean temperatures of the 1st and 2nd column soils between times t and t-dt, respectively.
, , is the heat-transfer resistance between the 1st column and the 2nd column soils: 2) For soil-i model: From Equation (10): From Equation (14):

Overall System of Equations
Based on the abovementioned development of the mathematical model, each divided element or layer has (m+2) equations. The total number of equations is n × (m + 2). A matrix equation was thereafter developed as follows: From Equation (18), the dynamic temperature was calculated by Equation (19). • The air, tube wall, and soil temperature for each divided element at time t can be calculated as follows: where is the vector of the unknown mean temperature between times t and t-dt; is the corresponding matrix of coefficients, which are related to the air parameters, tube parameters, soil parameters, air velocity, etc.; and T t are the air temperature, tube wall temperature, and soil temperature values at times t-dt and t, respectively; and is the column vector of the known terms, e.g., inlet air temperature of T in and outermost soil temperature of T g,j .

Cooling Capacity Calculation
Based on the calculated air temperature, the cooling capacity of the UAT system at time t can be calculated by the following Equation (21) [45].

Q t =m × C a × (T in t -T out t )
where m and C a refer to the air mass flow rate and air specific heat in the tube. T in t and T out t refer to the inlet and outlet temperature at time t, respectively.

Model Validation
The prediction accuracy of the developed dynamic model was compared with the experimental results. The experimental platform of the UAT system was built in Changsha, China. The UAT system consisted of a power equipment and S-shaped buried tube with a depth of 3.5 m and a total length of 34 m. The tube was made of polyethylene (PE) with a diameter of 250 mm and wall thickness of 10 mm. The sensors are PT 100 thermocouple with a permissible error at ± 0.15 °C. All these sensors were calibrated using the ice-water mixture before being connected to the Agilent 34972A for the experimental testing. The airflow velocity in buried tube was measured using the vane probe type anemometer (TSI-9555) with an accuracy of ±3.0%. The schematic diagram of UAT system are presented in Figure 4. The input parameters used for the dynamic model validation are listed in Table  1.  The soil parameters in this paper were estimated from soil types of the experimental site according to Reference [46]. The soil temperature was experimentally measured. The experimental validation was conducted for the cooling mode in summer. The experimental testing in summer was conducted from August 23rd to 24th, 2016, with a continuous operation for 24 hours. The air velocity was 1 m/s during the experimental testing process. Figure 5 shows the comparison results between the experimental testing and simulated results for outlet air temperature and cooling capacity in summer. As shown in Figure 5, the dynamic modeldriven outlet temperature shows a similar trend to the monitored temperature from 15:30 to 7:30, whereas the dynamic model-driven outlet temperature shows a higher increasing trend than the experimental results. The temperature rise of the outlet air temperature is 0.61 °C for the experimental results, and the temperature rise is 1.43 °C for the simulated results. The underlying mechanism is due to the testing errors, resulting from the PT 100 thermocouple and the vane probe type anemometer. Furthermore, the maximum temperature difference between the monitored and simulated air temperature at the outlet is less than 0.81 °C. The corresponding average temperature difference is less than 0.13 °C. Furthermore, the maximum absolute relative deviation (MARE, as shown in Equation (22) [47]) between the monitored and simulated results for the air temperature at the outlet is approximately 3.18%.
where is the number of selected data samples, 1, 2, … , and where , and , are the simulated and monitored air temperatures, respectively. Based on Equation (21), the maximum cooling capacity between the monitored and simulated results can be calculated as 39.59 W and the corresponding average cooling capacity difference is 8.90 W. Furthermore, the maximum absolute relative deviation between the monitored and simulated results for the cooling capacity can be calculated as 6.03%. According to the above comparison results, the developed dynamic model is accurate to further investigate the dynamic thermal performance for multidimensional parameter analysis.

Thermal Performance Analysis
In this section, the multidimensional parameter analysis has been conducted on geometrical parameters and operating parameters of the UAT system. Based on the developed dynamic model, the outlet temperature and cooling capacity of UAT system were predicted with respect to different tube lengths, tube thermal conductivities, tube diameters, and air velocities.

Effect of Tube Lengths
The effect of the tube length on the outlet temperature and cooling capacity was investigated, with the air velocity at 1 m/s and tube inner diameter at 250 mm. Other variables, e.g., inlet air temperatures, tube conductivities, and soil parameters, are kept constant.  Figure 6 shows the impacts of the tube lengths on the cooling capacity and outlet temperature. The average outlet temperatures and cooling capacities are presented in Table 2. As shown in Figure 6 and Table 2, when the tube length is 10 m, the UAT system shows the highest outlet temperature at 28.15 °C and the lowest cooling capacity at 163.29 W. With respect to the increase of the tube length from 10 to 60 m, the outlet temperature decreases from 28.15 °C to 24.40 °C and the cooling capacity increases from 163.29 to 401.95 W. The underlying mechanism is that a longer tube length can enhance the heat exchange between the flowing air and surrounding soil due to the increased timeduration of flowing air in the tube. However, as shown in Figure 6, with the increase of the tube length, the decreasing magnitude of the outlet temperature and the increasing magnitude of the cooling capacity decrease. For instance, when the tube length increases from 10 to 20 m, the decreasing magnitude of the outlet temperature is 1.57 °C and the increasing magnitude of the cooling capacity is 100.15 W. However, when the tube length increases from 50 to 60 m, the decreasing magnitude of the outlet temperature is only 0.23 °C and the increasing magnitude of the cooling capacity is 14.83 W. The underlying mechanism is that, as the decreasing magnitude of the soil temperature decreases with the increase of the tube length, the temperature difference between the flowing air and surrounding soil temperature will be saturated. By contrast, from the perspective of the initial economic investment, the increase of the buried tube length will significantly increase the initial, operational, and maintenance costs due to the difficulties in the system installation and operation processes. Therefore, when systematically considering the outlet temperature, the cooling capacity, and economic performance, the solution for the increase of the buried tube length is not techno-economically competitive, especially with the increase of the buried tube length.

Effect of Tube Thermal Conductivities
The effect of different tube conductivities on outlet temperature and cooling capacity is investigated, when the air velocity, tube length, and tube diameter are set at 1 m/s, 32 m, and 250 mm, respectively. Other variables, e.g., inlet air temperatures and soil parameters, are kept constant.  Figure 7 shows the impacts of tube conductivities on the outlet temperature and cooling capacity. The average outlet temperatures and cooling capacities are listed in Table 3. As shown in Figure 7 and Table 3, the increase of tube conductivities can decrease the outlet temperature and can increase the cooling capacity. For instance, with the increase of the tube coefficient of heat conductivity from 0.5 to 16 W/(m K), the average outlet temperature decreases from 25.30 °C to 25.16 °C and the average cooling capacity increases from 344.43 W to 353.62 W. Furthermore, with the increase of coefficient of heat conductivity of the tube, the decreasing magnitude of the outlet temperature and the increasing magnitude of the cooling capacity decrease. For instance, with the increase of the coefficient of heat conductivity of the tube from 0.5 to 1 W/(m K), the decreasing magnitude of the outlet temperature is 0.07 °C (from 25.30 °C to 25.23 °C), and the increasing magnitude of the average cooling capacity is 4.71 W (from 344.43 to 349.14 W). With the increase of the coefficient of heat conductivity of the tube from 1 to 16 W/(m K), the decreasing magnitude of the outlet temperature is 0.07 °C (from 25.23 °C to 25.16 °C), and the increasing magnitude of the average cooling capacity is 4.48 W (from 349.14 to 353.62 W). In can be concluded that both the outlet temperature and average cooling capacity are less dependent on the coefficient of heat conductivity of tube materials. Therefore, the solution for the increase of the coefficient of heat conductivity of tubes is not quite techno-economically competitive for the performance enhancement.

Effect of Tube Diameters
The effect of different tube diameters on the outlet temperature and the cooling capacity is investigated in the scenario when the air velocity and tube length are 1 m/s and 32 m, respectively. In this study, the selected maximum tube diameter is 400 mm considering the assumption of thermal interference between the adjoining tubes in Section 2.2. As reported in Reference [48], a distance of more than 2D (D is the tube diameter) between the adjoining tubes enables the absolute mean differential of the ratio of the annual exchanged heat rate to be less than 2.1%. Other variables, e.g., inlet air temperatures, tube conductivities, and soil parameters, are kept constant.   Figure 8 shows the impacts of tube diameters on the cooling capacity and the outlet temperature. The average outlet temperatures and cooling capacities are presented in Table 4. As shown in Figure  8 and Table 4, the outlet temperature and cooling capacity of the UAT system increase with the increase of the tube diameter. Specifically, with the tube diameter increasing from 100 to 400 mm, the average outlet temperature increases from 24.05 °C to 26.67 °C and the average cooling capacity increases from 67.84 to 659.29 W. Furthermore, a larger tube diameter is more conducive to improving the cooling capacity. For instance, the increasing magnitude of the average cooling capacity is 78.24 W (from 67.84 to 146.08 W) when the tube diameter increases from 100 to 150 mm and the increasing magnitude of the average cooling capacity is 105.51 W (from 553.78 to 659.29 W) when the tube diameter increases from 350 to 400 mm. However, it is noteworthy that a larger tube diameter will result in a higher outlet temperature, resulting in the increase of the space cooling load for indoor thermal environment. For instance, the maximum outlet temperature can reach 28.96 °C when the tube diameter is 400 mm, and thus, a high temperature has difficulty satisfying the requirements of indoor thermal comfort. In addition, a larger tube diameter will also result in a higher system's initial investment. Thus, an appropriate tube diameter should be carefully selected based on the balance between the outlet temperature and cooling capacity in practical application.

Effect of Air Velocities
The effect of different air velocities on the outlet air temperature and cooling capacity is investigated in the scenario when the tube length is 32 m and the tube diameter is 250 mm. Other variables, e.g., inlet air temperature, tube parameters, and soil parameters, are kept constant.   Figure 9 shows impacts of air velocities on the outlet temperature and cooling capacity. The average outlet temperatures and cooling capacities are presented in Table 5. As shown in Figure 9 and Table 5, both the outlet temperature and cooling capacity increase with the increase of the air velocity. To be specific, with the increase of the air velocity from 0.5 to 4 m/s, the average outlet temperature increases from 24.82 °C to 27.05 °C and the cooling capacity increases from 187.50 to 931.96 W, respectively. The main reason is that a larger air velocity leads to less time-duration for the heat exchange between the surrounding soil and the flowing air. Figure 9 also shows that the outlet air temperature is more stable and lower with the decrease of air velocity. However, the smaller air velocity will lead to a lower cooling capacity. Thus, the trade-off between the lower outlet temperature and larger cooling capacity has to be made for practical applications.

Conclusions and Future Studies
In this study, an underground air tunnel (UAT) system was proposed for the utilisation of geothermal energy for space cooling of buildings. A dynamic mathematical model of the underground air tunnel system was developed on the MATLAB/SIMULINK platform for the dynamic performance prediction, together with an on-site experimental test-rig for the validation. Compared to conventional UAT model developed on the CFD platform, the developed dynamic model on the MATLAB/SIMULINK platform is more flexible for the multi-domain simulation with higher flexibility for the multidimensional parametrical analysis. Based on the developed dynamic model, the multidimensional parametrical analysis has been conducted on geometrical parameters (i.e., tube lengths, tube diameters, and tube thermal conductivities) and operating parameters (i.e., air velocities) to provide technical guidance for the optimal design and robust operation. The main results of this study can be summarized as follows: 1. The on-site experimental validation indicates that the proposed mathematical model is accurate for the dynamic performance prediction. Compared with the on-site experimental data, the developed model is credible and accurate with a maximum absolute relative deviation of 3.18% and 6.03% for outlet air temperature and cooling capacity, respectively 2. The comprehensive and systematic parametrical analysis results indicated that, with the increase of the tube length, the outlet temperature decreases and the cooling capacity increases significantly whereas the decreasing magnitude of outlet temperature and the increasing magnitude of cooling capacity reduce 3. The system performance is slightly dependent on the tube materials, and the technical solution, e.g., the increase of the coefficient of heat conductivity of the tube, is not efficient to improve the system performance 4. Both outlet temperature and cooling capacity of the UAT system increase with the increase of tube diameter and air velocity. The increase of tube diameter and air velocity can improve the system's cooling capacity. However, the decrease of tube diameter and air velocity will result in a more stable and lower outlet temperature, which can reduce the cooling load for indoor thermal comfort The above conclusions are drawn under a 24-hour continuous operation. However, the outlet air temperature and cooling capacity of the UAT system are dependent on the operation time due to the gradual heating of the soil around the buried tubes. Therefore, the assessment of the cooling potential should be carried out over a long-time operation period, e.g., several months, in our following-up study. Further studies need to be conducted for the trade-offs between tube diameter and air velocity for the bi-criteria performance enhancement between cooling capacity and outlet temperature. Furthermore, multivariable optimizations [49] with multi-criteria performances, such as cooling capacity, coefficient of performance (COP), and economic performance and system efficiency, will be conducted in our following-up study. In addition, geometrical configurations of phase change material systems (e.g., displacement position of PCM outside and inside the tube or PCM in the surrounding soil) will be experimentally investigated for technical performance enhancement.