Abstract
The present study reports air-based novel energy saving compact bank type earth air heat exchanger system (EAHE) (having only 25m2 area) fabricated from Polyvinylchloride pipes of nominal diameter 0.203 mm (8 inches) for inlet and outlet headers and main pipes and 0.150 mm (6 inches nominal diameter) for branched pipes, which is installed in Shaheed Bhagat Singh State University, Ferozepur in Punjab state, India (North West border area of India) and studied for the duration from April to August 2021 using the full factorial design with four inlet variables at multilevel to find a model equation between inlet and outlet variables, i.e., heat exchanger effectiveness and achieved temperature difference. An induced airflow mode has been used to provide near-uniform flow through all corresponding pipes of symmetrical-shaped earth air heat exchangers. The values of Dry bulb temperature (DBT) and wet bulb temperature (WBT) have been measured for air inlet, for air outlet, as well as for ambient air using resistance temperature detectors (RTD) installed at the requisite locations. Earth's undisturbed temperature (EUT) was also noted by installing RTD at a depth of 2 m, and the average measured value of EUT is found to be 28.5°C. The current system has imparted the temperature difference variation from 0.4 oC to 9.4 oC, and the effectiveness varied from 0.16 to 0.82 during the whole season. This system could give a cooling potential of 0.0523 kJ/s to 1.587 kJ/s using a mass flow rate of 0.163 kg/s.Experiments have been designed methodically to apply the full factorial technique. Most favorable parameters have been found for hot and dry and hot and humid weather. The current study is novel in terms of significantly improving the effectiveness of EAHE and addressing the central issue of space limitation in urban areas.
Keywords
Introduction
An earth air heat exchanger helps regulate the temperature of ambient air flowing through it before being fed to the conditioned space. These are placed 1–3 m below the ground surface, and the air may pick up or lose heat to the surrounding soil depending on weather conditions. The concept behind this system is that the temperature below the earth's surface beyond a certain depth becomes independent of ambient air temperature, which is called Earth's Undisturbed Temperature (EUT). The EUT is higher and lower than the ambient temperature in the winter and summer seasons, respectively. When ambient air is passed through the tubes buried underground at a considerable depth, the air becomes warmer in winter and cooler in summer than ambient air. So it is considered to be a promising technique for energy savings. Although it is a very energy-efficient option, it has not been used on a large scale as it requires considerable space. EAHE is a very inert cooling technique with no environmental problems and consumes minimal electrical energy. India's Energy Conservation Building Code (ECBC 2017) has considered the EAHEs, the best air conditioning technique for the buildings sector. Almost 40% of the world's high-grade energy is consumed in the heating and cooling of buildings. Combating global warming through reduced energy consumption in buildings is one of the European Union strategies. The use of renewable energy sources is being encouraged across the globe; 1 thus, EAHEs ideally provide comfortable temperatures for preheating and pre-cooling in the winter and summer seasons, respectively.
Climatic conditions affect the performance of EAHE. Studies have been carried out on EAHEs across the globe under different climatic conditions.2–5 These studies have reported that EAHEs work best in semi-arid and temperate climates but worst in humid climates.5,6 The places with alternate heating and cooling demands have been shown to utilize this technology optimally. 3 Soil composition also affects the performance; hence the effect of soil composition has been studied,7–11 and it was concluded that sandy soil gives the best performance as far as heat transfer is concerned. 8 The soil with a higher humidity level performs better than dry soil. 12 It has been reported that the addition of Bentonite to soil can improve the performance of EAHE; thus, the length of the tube can be reduced, and the increase in energy gain achieved is 13–14%. 11 The thermal conductivity of dense soil is high and thus gives better performance. 13 The material of pipes does not have much impact on the performance of EAHE. 14
Extensive testing has been done on EAHEs in both physical and simulative modes. A few studies have reported actual fabrication and installation of EAHE in full scale. Many studies have been conducted which have used simulation techniques15–19 using software, i.e., CFD, ANSYS, TRANSYS. Numerical modeling20–23 and analytical modeling at different parametric values using different configurations have been done to a large extent. Most of the experiments have been performed with pipes laid in a straight line24–35 and serpentine style36–39 in both open mode 40 and closed mode 41 to study their usefulness. Large space requirement is a significant constraint for that. So pipes were laid in meandering style by other researchers 27 to overcome this problem. U, S, and Z configurations have also been used.22,39,42–44 Temperature difference ranging from 5°C to 20°C has been reported. In one of the studies, a higher temperature difference of 21.3 has been achieved with a water jacket on the outlet pipe. 45 The longer the EAHE, the more the effectiveness as the outlet temperature gets closer to EUT. Vertical borehole heat exchangers have also been studied with groundwater flow. A new convective ground resistance was added for a peclet number greater than 0.01. 46
Some studies are based on combined systems where EAHE is coupled with evaporative coolers, desiccant dehumidification systems (experimentally47,48 as well as numerically,49–52) and vapor compression systems42,53,54 to improve its competence in the air conditioning industry. Energy savings of 71%, 49 74%, 55 and 75% 56 and CO2 emissions savings of 76 to 91% 52 have been reported in coupled systems. Indoor air quality can be improved by coupling EAHE with a solar-assisted desiccant cooling system also. 57
The ground heat source pump was examined for cooling and heating purpose and the co-efficient of performance was found to be 4.29 for cooling and 3.97 for heating purpose in Korea. 58
Bank type EAHE can be easily fitted in while laying the foundation of buildings as compact and does not meddle with the civil foundation restraints. EAHEs with long straight pipes can be practicable in countryside areas, but metropolitan cities are unsuitable due to space constraints. The cost of land, as well as quarrying, is high in urban areas. So for the present study, a compact bank type EAHE(requiring only 25m2 area) was conducted in Shaheed Bhagat Singh State University, Ferozepur in Punjab state, India, from April to August 2021.
Materials and methods
Bank shaped earth air heat exchanger consisting of two headers (inlet and outlet) with four parallel distribution lines has been decided. EAHE was fabricated with Polyvinylchloride pipes of nominal diameter 0.203 mm (8 inches) for inlet and outlet headers and main pipes and 0.150 mm (6 inches nominal diameter)for branched pipes. Leak-proof joints were made between the lines. The symmetrically shapedEAHE was installed in a trench, as shown in Supplementary Figure 1. The diameter of the main pipes is 1.33 times the diameter of connecting lines, as the literature confirms that the primary pipe diameter should be 1.2 to 1.8 times the parallel branched pipes. 43 The labeled top and front views of the header connected to the main pipe, with all dimensions, are shown in Supplementary Figure 2.
Two similar setups for experimentation were created and installed at the depths of 1 m and 2 m below the earth's surface. Induced draught of airflow is supplied using a circular duct fan. The values of Dry bulb temperature (DBT) and wet bulb temperature (WBT) have been measured at (a) air inlet (b) air outlet (c) ambient air, with the help of resistance temperature detectors (RTD) installed at the requisite locations. Earth's undisturbed temperature (EUT) was also noted by installing RTD at a depth of 2 m, away from the effect of EAHE. Instruments used of the required specifications are given in Table 1, and their pictures are given in supplementary Figure 3.
Instrumentation and component used.
Many observations have been recorded daily from April to August 2021 at two different depths and three airflow velocity values. The values of DBT and WBT have been measured with the help of RTDs with an uncertainty level of ± 0.1 °C, and velocity of air has been measured with Vane type anemometer with the measurement uncertainty of ± 0.1 m s − 1. The airflow rate has been calculated using airflow velocity, and the uncertainty falls within ± 2% of the measured value. The specific humidity of air was calculated from the DBT and WBT of inlet and outlet air, respectively. The uncertainty in the calculated values from the observed values is found to be ± 2.67 × 10–4 kg kg − 1 of dry air. Output response variables are the effectiveness of the heat exchanger and temperature difference. They have been calculated from the expressions given below.
Cp is the specific heat at constant pressure.
The detail of inlet parameters and response variables, along with their minimum and maximum values, is given in Table 2.
Comparison of results with previous studies.
Results and discussion
Bank type EAHE has not been explored much. Previous studies have also not explained the effect of specific humidity as well as the effects of 2-factor or 3-factor interactions. The current study is very much novel in terms of the shape of EAHE and its compact design as it occupied an area of only 25m2 and the values of effectiveness (
Effect of factors interactions on dT
The
It has been noted from the graph (Figure 1) that the value of

Variation of dT with Tdi for three different velocities at a specific humidity of 0.015 kg/kg of dry air (a) at h = 1m (b) at h = 2 m.

Variation of dT with Tdi for two different depths at specific humidity of 0.015 kg/kg of dry air and velocity of 3.2 m/s.

Variation of dT with Tdi for three different specific humidities at a velocity of 3.2 m/s (a) at h = 1m (b) at h = 2 m.
The similar trends as shown in Figures 1 and 2 were obtained in earlier literature studies. In one such study, temeprature difference of 6.2 oC was obtained at a depth of 1.5 m with a length of 60 m 40 and, in another study, temeprature difference of 8.4 oC was obtained at a depth of 1.5 m with a length of 60 m 30 both conducted in Algeria. One more study has reported a temperature difference of 8 oC at a depth of 1 m with a length of 39 m 38 in India. However, the current study has reported temeprature difference of 9.4 oC at a depth of 2 m at Ferozepur, India. Moreover, the detailed comparison of present study with previous studies is shown in Table 2. However the results reported in current study are better in terms of temperature difference as well as size and shape.
The
To the best of knowledge, no study has reported the values of specific humidity in literature for EAHE system, which makes the current study novel.
It is noted that the lowest value of 0.32°C of
Effect of factors interactions on effectiveness (
)
Figure 4 shows the variation of

Variation of ɛ with Tdi at two different depths (a) at w = 0.0105 kg/kg of dry air and v = 2.4 m/s (b) w = 0.015 kg/kg of dry air and v = 3.2 m/s ( c) w = 0.0206 kg/kg of dry air and v = 4.0 m/s.
Similar trends were obtained in literature and one such study has reported effectiveness of 0.79 at a velocity of 3.5 m/s and inlet temeprature of 43.5 oC 37 in Algeria and in another study, effectiness of 0.8 at a velocity of 3.5 m/s and inlet temeprature of 49 oC is reported 27 in Algeria. The current study has reported efectivess of 0.82 at a velocity of 2.4 m/s with inlet temeprature of 43 oC and specific humidity of 0.105 kg/kg of dry air.
Statistical modeling
Full Factorial design with 04 inlet parameters (Table 3) at multilevel has been used to analyze 02 response variables (Table 2). All input, as well as outlet parameters, are numeric and continuous. Depth of buried pipes as inlet parameter has been taken at two levels, whereas other three parameters, i.e., DBT, specific humidity, and airflow velocity, have been taken at three levels. So a set of 54 experiments has been taken from the whole set of readings. Seven extra design points have been taken to improve the model's accuracy.
Detail of parameters along with their minimum and maximum values.
The data analysis has been done with the help of software - Design Expert 13.0.5.0. Analysis of data and prediction of response variables was tried with the help of linear, 2FI, quadratic, cubic, quartic, fifth, and sixth-order models, and it was seen that the reduced 2FI model fitted well for the analysis of effectiveness as well as temperature difference as response variables. Though the higher-order models gave better fitting polynomial, the highly complex equation became cumbersome to use. These models have an acceptable range of residual errors. Significant terms were identified after running AICc forward regression on reduced 2FI models. Some nominal terms were added to maintain the hierarchy of variables to make the analysis of variance viable. Factor coding was coded, and the Sum of squares was Type III – Partial. The Model F-value proves the significance of the model. Only those model terms are significant that have P-values less than 0.0500. The difference between Predicted R2and Adjusted R2 is less than 0.2. Adeq Precision measured the signal-to-noise ratio. A ratio greater than 4 is desirable. All these values for both models have been tabulated in Table 4. The values conform to the requirements for valid models. So these models can be used to navigate the design space.
Models for response variables .
The normal plots of residuals versus various terms were plotted, and it was noted that maximum points lie on the straight line, and some points were scattered. No specific S pattern could be observed, proving that a normal distribution pattern was followed and no further transformation was required. The final equations for predicting temperature difference and effectiveness are given below.
Final equation in terms of actual factors:
The graphs indicate that the model equations were appropriate for predicting the values of

Predicted vs. actual values (a) for response 1: DT (b) for response 2:
Point prediction and model validation:
The comparison has been made between predicted and observed output parameters (Table 5) at given inlet parameters that were not used for modeling.
Predicted Vs. Observed Values.
The predicted values are almost consonant with the observed values, confirming the models' validity. Hence the present study has explored an energy-saving natural compact air conditioning system that can be commercialized for urban buildings.
Response Surface Model has been used in literature for Parametric optimization of vertical ground heat exchanger 59 and parametric optimization of horizontal earth to air heat exchanger.27,60 No study has been found to the best of knowledge for the optimization of EAHE based on Full Factorial design which makes the current study novel.
Conclusion
The compact-shaped bank type EAHE is very effective in pre-cooling the air that can be directly fed to the conditioned space or the air handling unit of the conventional air conditioning system, depending on weather conditions. The effectiveness was found to be 0.82 in the present study, which is comparatively better than many studies giving a significant improvement in effectiveness. As the value of temperature difference decreased with the increase in moisture content of the air, it was found that the system was more effective in hot and dry weather. The temperature difference varied from 0.4°C to 9.4°C, and the effectiveness varied from 0.16 to 0.82 during the whole season. The system can provide a cooling potential of 0.0523 kJ/s to 1.587 kJ/s at a mass flow rate of 0.163 kg/s. The airflow velocity had a minimal effect, and lower velocity increased the temperature difference. Increased depth of EAHE increased the temperature difference by 2–3 oC. A polynomial has been fitted between the inlet parameters and response variables using full factorial design, taking 63 experimental readings fitted into the design table. The model equations have given −5% to 5% error between the experimental and predicted values. So the model is considered to be appropriate. Hence, the bank type compact-sized EAHE is seen as a promising and suitable energy-saving alternative or to support conventional air conditioning systems, especially in urban areas where space availability is a constraint.
Supplemental Material
sj-docx-1-pie-10.1177_09544089221096052 - Supplemental material for Experimental studies on natural bank type heating/cooling energy saving Air conditioning system
Supplemental material, sj-docx-1-pie-10.1177_09544089221096052 for Experimental studies on natural bank type heating/cooling energy saving Air conditioning system by Vaishali Goyal, Arun Kumar Asati, Rajeev Kumar Garg and Amit Arora in Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering
Footnotes
Acknowledgements
The authors are highly grateful to the university's civil and electrical maintenance staff for helping with the experimental setup's fabrication and installation.
Disclosure statement
There is no conflict of interest with anyone.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Supplemental material
Supplemental material for this article is available online.
Appendix
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
