Refurbishment of a Social Interest Building in Mexico Using Earth-to-Air Heat Exchangers

: The refurbishment of a social interest building using Earth-to-Air Heat Exchanger (EAHE) was studied in representative dry climatic conditions of Mexico (dry, very dry, temperate, and sub-temperate). A simulation method that uses both computational ﬂuid dynamics (CFD) and building energy simulation (BES) was used to analyze the inﬂuence of the EAHE on the indoor conditions of a room. First, CFD simulations of the EAHE were performed using climatic data and soil properties of the four representative cities, and then the results were loaded into the TRNSYS software to estimate the indoor air temperature and the building room’s thermal loads. When connected to a building room on a warm day, the EAHE reduced the indoor air temperature by a factor ranging between 1.7 and 3.2 °C, while on a cold day, the EAHE increased the indoor air temperature of the room by between 1.0 and 1.9 °C. On the other hand, the EAHE reduced the daily cooling load of the room by a factor between 2% and 6%. The EAHE also reduced the daily heating load by between 0.3% and 11%. Thus, EAHE as a refurbishment technology can beneﬁt social interest buildings in Mexico.


Introduction
Currently, daily human life is very attached to technology. For instance, we use air conditioning systems to achieve thermal comfort. This technology allows activities to be carried out more comfortably, especially in zones with extreme weather conditions. However, air conditioners are large electricity consumers; thus, energy refurbishment technologies currently seek to condition indoor spaces with the lowest possible energy consumption. One of these refurbishment alternatives is the Earth-to-Air Heat Exchanger (EAHE), which supplies airflow to the indoor environment of buildings [1][2][3]. This airflow improves indoor air conditions without consuming a high amount of electricity because of the use of the thermal inertia of the soil. In an EAHE, the ambient air is supplied by a fan to one or more pipes buried in the soil at a depth where the outdoor environment has a small influence. Because the soil temperature is different from the ambient temperature, the airflow in the pipes exchanges heat with the soil, which acts as a heat source or sink depending on the season or climatic zone [4]. Thus, the air from the buried pipes is cooled or heated and sent to the buildings to benefit their indoor thermal comfort. The amount of heat exchanged between the air in the pipe and the soil depends on various factors such as weather conditions, pipe length, diameter, and soil properties. Several researchers have studied the EAHE by performing experimental tests. Those studies have provided information about the effectiveness of EAHE in reducing or increasing airflow temperature. The capacity of EAHE to cool or heat the air depends on its design and the thermophysical properties of the soil. The following points can be highlighted regarding the design of EAHE: • An EAHE installed in soil with a high humidity level presents better thermal performance and better air temperature drops than EAHE installed in dry soil [5,6]. • An airspeed between 1 and 5 m/s provides good airflow at the EAHE outlet; however, this reduces the time of interaction between the ground and the air, reducing the temperature difference between the inlet and the outlet [7]. • There are no significant changes in results concerning the material of the pipes [8]. • The temperature of the soil surrounding the pipe increases as time passes, which can be between 6 and 12 h of continuous operation of the system, and this increase in temperature decreases the capacity of the EAHE to reduce the temperature of the air at the outlet [9][10][11]. • A diameter between 25 and 30 cm reduces pressure drops in the pipe [12]. • The drop or rise in temperature to the air caused by the EAHE can reach values between 2 and 20°C and the cooling or heating capacity ranges between 2000 and 8000 W. Those values vary depending on weather conditions, construction parameters, and soil type [13][14][15][16][17][18][19][20].
The previous findings have allowed a better understanding of how an EAHE behaves; however, studying an EAHE by performing experimental tests requires considerable time and might be expensive. Thus, EAHE has also been studied theoretically to subject it to different conditions and understand the influence of various parameters on its thermal behavior. Many of these studies have certain similarities in developing their mathematical models, such as the following: constant soil and air thermophysical properties, constant air velocity, and laminar, transitional, or turbulent flow. These models are usually solved using computational fluid dynamics (CFD) [21]. On the other hand, most analytical studies consider one-dimensional models, generally whose formulations are solved using the thermal resistance method. Another notable difference is that the analytical studies consider steady-state conditions, while most numerical studies consider EAHE in an unsteady state [22]. These differences make analytical studies easier to develop but limited in representing reality, while numerical studies have higher accuracy but are relatively complex. Other computational tools, such as TRNSYS, implement components (type 952 and 997) where, by means of finite differences, they model the heat conduction between the pipe and the surrounding soil.
Some works have analyzed the influence of the soil temperature on the behavior of EAHE [12,[23][24][25][26][27]. Other studies have considered different passive systems coupled to the EAHE [28][29][30][31][32]. Some studied different EAHE configurations [2,[33][34][35] or different weather and soil conditions in which the EAHE interacts [36]. The results of these studies have allowed us to understand which soil conditions benefit the EAHE. For instance, a high humidity level in the soil or soil moisture benefits the thermal potential of the EAHE [37,38]. The typical configuration of an EAHE is a horizontal pipe. However, for its optimal operation, long lengths are sometimes needed, so various designs can allow long lines in a reduced space, such as the spiral configuration, serpentine, and others known as Z because of the way the fluid travels [39][40][41][42][43]. In the case of climatic conditions, the system works better for extreme weather conditions and deficient ambient humidity levels compared to mild and humid climates [34,44,45].
Other research works have mentioned that the EAHE can increase the ventilation potential of buildings. Thus, some researchers have evaluated the EAHE system connected to a room. Maytorena and Hinojosa [46] numerically analyzed the thermal behavior of an EAHE connected to a building room in the dry zone of Sonora, Mexico. The research was divided into four cases, where the authors determined the temperature inside the room while thermally insulating some walls and allowing heat gain in others. They showed that the EAHE maintained the indoor air between 27 and 29°C regardless of where the heat gain was activated. Other authors, such as Ahmed et al. [15], experimentally studied the conditioning capacity of the EAHE connected to a room in a warm humid zone in Rockhampton, Australia. The authors found that the EAHE reduced the indoor air temperature by 2.1°C, while the humidity level increased by 1.2%. Other works such as Skotnicka-Siepsiak [47] have studied the energy capacity of the EAHE connected to a room in Olsztyn, Poland. The results showed that the system reduced 45% of the energy consumption of the room by providing 75% of the cooling energy required by the room.
Other studies were conducted to analyze different passive systems coupled to a room. For instance, Long et al. [48] studied the benefits of connecting a solar chimney with phase change material (PCM) with an EAHE. The results showed that the airflow of the EAHE increases by 50% when coupling this system to a solar chimney with PCM and with 0.8°C lower indoor room temperature than the case without PCM.
An EAHE may have excellent cooling or heating potential, which is defined as the air temperature difference between the outlet and the inlet. Still, when the supply airflow to the room is insufficient, it could not benefit the indoor environment regardless of the supply air temperature. Previous works of the authors of the current research have numerically analyzed the performance of an EAHE for different cities in Mexico [49][50][51]. Although those studies showed that the EAHE had an acceptable thermal performance for most climatic zones, they do not reveal how the airflow from the EAHE could benefit the indoor environment of a building. Therefore, this work aims to estimate the influence of an EAHE on the indoor air temperatures and heating and cooling loads of a room that belongs to a social interest building situated in regions representing the different dry climate conditions of Mexico. For this purpose, a coupling was carried out between computational fluid dynamics (CFD) and TRNSYS (BES). First, a CFD model was used to perform simulations of the EAHE in four cities, and then the results were loaded into the TRNSYS V16 software through a data reader to estimate a building room's cooling and heating loads.

Physical Model
The EAHE was modeled with a typical horizontal configuration. The system consists of a 5 m long PVC D30 pipe buried at 2 m. The dimensions of the EAHE are due to the space available in low-income housing, whose dimensions are between 30 and 50 m 2 . A 0.05 m insulation layer was considered to minimize energy losses or gains in the vertical outlet pipe of the EAHE. It is supposed that the air that is introduced to the building room has the same temperature and velocity as the air at the outlet of the EAHE. Thus, it is assumed that the elbow that connects the EAHE with the room is adiabatic and does not interrupt the airflow or cause additional head losses. Figure 1 shows the geometry of the EAHE, while its dimensions are shown in Table 1.
The assumptions for numerical modeling are shown in Figure 2. For the soil section (Figure 2a), we assume the following: • Two-dimensional heat transfer. • Constant and temperature-independent thermophysical properties of soil. • The temperature of the soil is determined by solving the heat conduction equation affected by the boundary conditions of the surface exposed to the outdoor environment.
For the pipe (Figure 2b), we assume the following: • The pipe of the circular cross-section is modeled as a square cross-section of equivalent area [52].
• The heat conduction of the pipe wall is not considered because the thickness is very thin. • There is no evaporation or condensation in the pipes of the EAHE. • Convection heat transfer is the dominant mechanism in the pipes. • There is laminar and transitional airflow inside the pipes because it is considered that a low-power fan will supply the air.

Room
In Mexico, homes are classified according to their surface area and price, with the residential plus being the largest and most expensive, and those with the lowest price and size are called social interest. The latter is between 30 and 50 m 2 in size [53], and consists of two bedrooms, a living-dining room, a kitchen, and a bathroom, as shown in Figure 3.  The room selected for the coupling of the system is the largest bedroom of the building from Figure 3 (room B), with dimensions of 3 × 2.85 m 2 . The room is highlighted in the Figure and the orientation of the building is also shown. The room's height was established according to the housing building code, which specifies that it must be 2.7 m [53]. Table 2 shows the thermophysical properties of the construction materials of the analyzed building used in the simulation together with the thicknesses of each one, whose values were obtained from the TRNSYS materials library. These materials were used to form the walls, roof, and floor configuration, which were defined as a composition of one or several layers of material, from the inside to the outside; the building envelope components are composed as follows: • Exterior walls (exposed to outdoor environment): a layer of gypsum plaster, a layer of mortar, brick, and a layer of mortar. • Interior walls (exposed to indoor environment): a layer of mortar, brick, and a layer of mortar. • Floor: a layer of marble, lime cement mortar, and a layer of light concrete. • Roof: a layer of gypsum plaster, a layer of mortar, hollow block, lime cement mortar, and a waterproofing layer.

Mathematical Model of the EAHE
The system analyzed is formed by the soil and EAHE. The mathematical model considers the heat conduction in the soil and heat transfer by convection in laminar and transitional regime inside the pipes. The governing equations of the heat transfer phenomenon are continuity, momentum, and energy, which are discretized using the finite volume method. These equations can be solved in a transient state; however, it requires a lot of resources and computational time. Thus, the quasi-transient state was used, which consists of modeling steady states every ten minutes. The governing equations are as follows: The system is subject to four boundary conditions, three adiabatic and one in which an energy balance occurs, proposed by Mihalakakou et al. [22]; these conditions are detailed below.
(a) Left and right boundaries (x = 0, x = H x ) are considered adiabatic: For the bottom boundary, the current work uses the concept of penetration depth used by some research works available in the literature [54,55]. This concept assumes that the daily heat wave (the daily oscillation in temperature of the outdoor air entering the EAHE) does not penetrate the soil beyond a certain radius, and this radius is lower than 1 m in the y downward direction. Therefore, this boundary can be considered as an adiabatic boundary after a certain depth. The bottom boundary can be expressed as follows: (b) Bottom boundary (y = 0): (c) Upper boundary (y = H y ). To represent the heat transfer in the northern boundary on the soil surface, the following energy balance proposed by Mihalakauko et al. [22] was used: (i) CE is the convective energy exchanged between the air and the soil surface: T amb is the temperature of the environment according to the day, and T sur is the temperature of the soil surface. h sur is the convection heat transfer coefficient on the soil surface and it is determined from the following equations [56]: (ii) SR is the short wavelength solar radiation absorbed from the soil surface, which is calculated by Mihalakauko et al. [22]: κ is the absorptance of the soil and G is the solar radiation that reaches the surface of the soil.
(iii) LR is the long wavelength radiation, which can be calculated as follows: ε is the emittance of the soil surface and ∆R is a term that depends on the relative humidity of the earth and air at the surface. A value of 63 W/m 2 has proven to be a good approximate value for this variable, shown in Badescu [56].
(iv) LE is the latent heat flux from the soil surface due to evaporation, which can be calculated as follows [22]: RH is the relative humidity of the ambient air and f is a number that is a function of the type of land cover. This fraction can take the following values: . It is considered that the temperature of air entering the pipes is equal to the temperature of the outdoor environment at a constant velocity, which is governed by the Reynolds number value.
For the numerical solution of the EAHE model, the governing equations of the system (soil + pipe thermal insulation + air) were solved as if it all was fluid, and the thermophysical properties were assigned according to the position in the system. In the solid region (soil and pipe thermal insulation), a blocking-off method was used, which consists of setting up the velocity components equal to zero. In this way, the hydrodynamic effect in the soil and insulation is canceled, and heat transfer is then governed only by heat conduction. In the fluid region (inside the pipes), such a blocking-off method was not used; therefore, the velocity and temperature fields were obtained for this region. In the current research, the soil temperature was not considered to be equal to the annual average outdoor air temperature, as the bottom boundary condition of the model. Instead, the soil temperature is determined by solving the mathematical model and its boundary conditions. The solution provides the temperature distribution of the soil, which depends on the type of soil, the upper boundary condition, which is affected by the outdoor environment, and the interaction with the pipes of the EAHE.

Numerical Modeling of the EAHE
The finite volume technique divides the physical domain into small control volumes. In this work, the finite volume technique uses a Cartesian grid and does not specify a z direction. Heat transfer in the z direction is considered negligible. A computational node is placed where the desired variable is calculated in each volume. Interpolation schemes were implemented to discretize the governing equations, using the central scheme for the conductive terms and a hybrid scheme for the convective ones. The SIMPLE algorithm was implemented for coupling the momentum and continuity equations [57,58]. A convergence criterion of 10 −10 was established for all variables to ensure good results and precision.
The optimal number of control volumes into which the domain can be divided varies in each case study. Thus, it is necessary to carry out a mesh independence study to find the number of volumes that allow the best precision in the results without consuming extra computational resources. Figure 4 shows the division of the physical domain in a computational mesh. The spaces separating the pipe volumes are more refined than those corresponding to the soil. The meteorological conditions for the mesh independence study were those of Juárez City because it is the place with the most extreme climate. The temperature was 46.6°C with solar radiation of 786 W/m 2 , and the relative humidity was 5% with a wind speed of 1.5 m/s.  The mesh independence study consisted of varying the number of nodes representing the pipe length (N x3 ) from 71 to 111 nodes, while the rest of the sections remained constant.
The results showed that after 101 nodes, there are no significant changes in the velocity and temperature of the air in the outlet. Once the number of nodes was established for the length of the pipe, the number of nodes varied in the sections of N x2 , N x4 , and N y2 from 50 to 90 nodes. From 71 nodes, it was observed that the outlet air temperature and speed no longer varied significantly. The number of nodes in the physical domain is 286 in the x-direction and 215 in the y-direction; the distribution is presented in Table 3. The computational code was verified by solving the exercise available in House et al. [59]. The exercise determines the velocities and temperatures of the air laminar flow regime in steady state by natural convection in a cavity with an embedded solid in the center. The cavity has a vertical hot wall (T H = 25°C) and a vertical cold wall (T C = 15°C). The dimension of cavity H is determined based on the Rayleigh number for conductivity ratios of 0.2 and 0.5 in which the Nusselt number was compared. The verification results showed that the absolute percentage difference between the results of this study and results obtained by House et al. [59] is 0.6% and 0.7% for the conductivity ratios mentioned above. Thus, the developed computational code provides reliable results. Other verification tests performed for the computational code used here can be found in previous publications of the authors of the current work [50,60].

BES-CFD Coupling
The BES-CFD coupling is a technique that allows us to take advantage of both parties in heat transfer analysis. The BES-CFD coupling can be of a dynamic type when both programs work together, exchanging information at each moment, or static when the programs work using the information received from the other, but they work independently. These benefits include reduced time and computational resources of convergence in CFD, more exact heat transfer coefficients used in BES, and better estimation of building energy requirements compared to standalone BES simulations [61]. In this study, a static coupling was implemented where the TRNSYS building modeling software implements the data obtained from the CFD numerical modeling through the following procedure: • From the numerical modeling of EAHE, a .txt file is obtained, which contains the air temperature at the outlet of the EAHE, the air changes per hour, and the air humidity of each modeled case. • The VENTILATION component is activated in the TRNbuild interface, and the variables are named. • The .txt file is read through a data reader in the simulation studio, which supplies the TYPE 56 component with the variables required by the VENTILATION component.
The BES-CFD coupling is carried out using a data reader and TYPE 56, as shown in Figure 5.
The room simulation in TRNSYS was carried out under free-floating conditions to analyze the influence of the EAHE on thermal comfort. The thermal loads for each chosen city's warmest and coldest days were also calculated to explore the effect of the EAHE on the thermal loads.

Climatic Conditions and Soil Properties
For the numerical modeling or CFD simulations, the ambient temperature, wind speed, and relative humidity of warmest and coldest days of four representative cities with dry climate present in Mexico were used as input data in 10 min intervals. These cities were selected based on their prevailing climate. The climatic data were provided by the National Meteorologic Service-National Commission of Water (SMN-CONAGUA). However, it is not possible to show all of the input data used in simulations, so here are presented only the average, maximum, and minimum values of the climatic data used in simulations (Table 4). These data were obtained from the warmest and coldest days of the year 2020. The same climatic data used for the CFD simulations were used to perform the simulations of the social interest building using TRNSYS. To perform the CFD simulations it is also necessary to have the properties of the soil in each of the representative cities in Mexico. Thus, according to previous research [62,63], the types of soil in each city are Mexico City (silt), Juarez City (sand), Zacualtipan (limestone rock), and Monterrey (limestone rock). Table 5 presents the thermophysical properties of the different types of soil considered in the current work.

Results
This section is divided into two parts. The first part presents the running schedule of the EAHE, which is determined by comparing the temperature of the air at the outlet of the EAHE and the temperature of the indoor air in the room. The second part presents the influence of connecting an EAHE on the indoor air temperatures and cooling and heating loads of the building room.

Determination of the EAHE Running Schedule
The numerical modeling of the heat exchanger was carried out for different air velocities, according to the Reynolds number in laminar and transitional flow. This section compares the temperature of the indoor air in the room, determined with TRNSYS, and the temperature of air supplied by the EAHE in each city, which was obtained with CFD simulation. This comparison allowed us to discriminate which schedule the EAHE works well with and which it does not. Therefore, if the EAHE does not comply with its duty of cooling or heating as needed, it is supposed that the air from EAHE does not enter the room. Figure 6a shows a graphical comparison between the indoor air temperature of the room without EAHE and the air temperature at the outlet of the EAHE on the warmest day in Monterrey. Different air velocities for the air into the pipes are presented in the figure as a function of the Reynolds number (Re). On this day, the indoor air temperature of the building room ranges between 29.5 and 45°C. The air from the EAHE is colder than the indoor air of the room the whole day for any value of Re. During the first seven hours, between 00:00 h and 07:00 h, the air from the EAHE was up to 15°C colder than the indoor air. This temperature difference tends to decrease as solar radiation increases; however, it remains below room temperature, with a difference of 7°C between 7:00 and 13:00 h for small Reynolds values. When the radiation decreases, around 16:00 h, the air temperature from the outlet of the EAHE decreases.
On the coldest day in Monterrey, the air temperature supplied by the EAHE and the indoor air temperature in the room are shown in Figure 6b. The indoor air temperature of the building room ranges between 6 and 18.6°C. Thus, the room in this city requires heating to reach thermal comfort, especially the period before 14:00 h. Figure 6b shows that the airflow provided by the EAHE is warmer than the air in the room during most of the day at high Re. As solar radiation increases, this heating effect tends to be greater from 08:00 h, with a temperature difference of about 11°C. However, as solar radiation decreases, the EAHE loses the heating effect. The supplied air temperature becomes lower than that of the room from 16:00 h, with a temperature difference of about 11°C. This behavior makes it unfeasible to keep the system running when the room reaches a temperature close to 15°C; thus, the EAHE could only provide 17 h of heating.  Figure 7a shows the comparison between the temperature of the air supplied by the EAHE and the indoor air temperature of the room in Juárez City. Throughout the warmest day, the air temperature supplied by the EAHE is lower than that reached by the room for all Re numbers. But, as Re is smaller, the temperature of the air supplied by the EAHE becomes colder than the indoor air. During the first hours, this temperature difference is close to 15°C when Re = 100, which decreases as solar radiation increases after 9:00 h. However, this temperature difference reaches between 7 and 10°C between 12:00 and 16:00 h. As solar radiation decreases, the supplied air recovers its cooling capacity for the rest of the day, reaching about 27°C below the indoor room temperature, which provides up to 24 h of cooling effect.
On the coldest day in Juarez City, the EAHE supplies air at a temperature higher than the temperature of the room for most of the day for high Reynolds numbers (Figure 7b). During the first hours, a cooling effect is obtained by the system towards the room of up to 10°C; however, as solar radiation increases, the heating effect of the system increases for all Re values, reaching a temperature difference of 15°C between 9:00 h and 14:00 h. This temperature difference tends to decrease as the radiation decreases after 14:00 h, until a cooling effect is obtained again after 18:00 h, so a heating effect of between 9 and 18 h is obtained depending on the speed at which the air is driven by the EAHE.  4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 Figure 8a presents the temperature of the air in the room and the air supplied by the EAHE in Zacualtipan. During the first hours of the day, a significant difference between both temperatures occurs (15°C). But, during the hours with solar radiation, this difference decreases between 8:00 and 13:00 h, reaching a value between 3 and 5°C, increasing again after 13:00 h. As solar radiation decreases, the effect of the air provided by the EAHE is cooling for the rest of the day, reaching a temperature difference of 20°C. Thus, the supplied air temperature is colder than the air in the room the whole day. Figure 8b shows that the air supplied by the EAHE during the first 8 h has a temperature lower than that of the room on the coldest day in Zacualtipán. As solar radiation increases, a heating effect begins after 8:00 h, reaching a temperature up to 11°C greater than the indoor air between 8:00 and 13:00 h. Subsequently, once the radiation starts to decrease, a cooling effect occurs again from 16:00 h, with a temperature difference of between 6 and 10°C, which can be harmful considering that the air in the room reaches values between 10 and 15°C. The maximum temperature increase is 10.5°C at 12:50 h, with a heating effect between 8:00 and 16:00 h.  4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23  The temperature of the air supplied by the EAHE and the indoor air temperature on the warmest day in Mexico City are presented in Figure 9a. The temperature of the air provided by the EAHE during the first hours of the day is between 7 and 12°C lower than the room air temperature. However, in the hours of maximum solar radiation, between 8:00 and 15:00 h, the cooling effect decreases until the supplied air temperature is higher, between 5 and 10°C, than the indoor air. As the day progresses and the radiation decreases, the cooling effect of the air supplied by the EAHE is recovered; the air supplied by the EAHE had a temperature between 10 and 16°C lower than the indoor air.
On the coldest day, Figure 9b shows that when solar radiation is available, the temperature of the air supplied by the EAHE is warmer, between 12 and 18°C (from 7:00 to 17:00 h), than the indoor air. On the other hand, for hours without solar radiation, the temperature of the supplied air is lower, between 2 and 5°C, than that reached by the indoor air, which causes a cooling effect. The EAHE provided a heating operation time of 10 h.  Table 6 summarizes the running schedules of the EAHE in the four cities analyzed in this work for the different values of Re. The table presents the hours in which the EAHE works well on the warmest day and the hours in which the EAHE works well on the coldest day. These data are presented in the column labeled working time (W.T.). The table also presents the schedules in which the EAHE complies with its cooling and heating duty for the warmest and coldest days, respectively.

Results of the Simulations of the EAHE Connected to the Building Room
This section presents the analysis of the EAHE connected to a room. The effect of EAHE on the indoor air temperature of the room and its cooling and heating thermal loads is presented here. Such results were obtained by coupling each city's CFD model of the EAHE and the BES room model.
As mentioned in previous sections, three input data are needed to feed the VENTI-LATION component in TRNSYS: rate of air changes, the temperature of the air from the EAHE that is injected into the house, and relative humidity. The air change ratio (ACR) was calculated according to [64]: The volumetric flow rate is the product of the average velocity of the fluid at the outlet of the EAHE (v ave ) and the area of the pipe (A). This volumetric flow is divided by the total volume of the room, in this case, a space of 20.52 m 3 . Table 7 reports the ACR values for each Reynolds number with different numbers of pipes implemented for a room. It is worth mentioning that pipe layout is such that no thermal interference between pipes occurs. The value obtained with one pipe was multiplied by the number of lines to obtain the ACR for more pipes. As expected, the ACR increases as the number of pipes and the Reynolds number increase. The peak ACH = 3.1 1/h is achieved for 20 pipes and Re = 4000. The effect of EAHE with 5, 10, 15, and 20 pipes on the thermal loads of the room was simulated. The maximum number of pipes was selected according to the area of the room (8.55 m 2 ). As the number of pipes increases, the airflow into the room increases. This is done to increase the ACR and improve the indoor conditions of the room. The airflow from the EAHE is supplied to Room B, as shown in Figure 3. The results from simulations of the EAHE connected to the room are presented below. It is worth mentioning that on a warm day, the EAHE only provides air to the room when the supplied air temperature is lower than that of indoor air, while for the coldest day, the EAHE only provides air to the room during the hours when the supplied air is warmer than the indoor air. It is supposed that the outlet of the EAHE is closed when this system does not comply with its duty of cooling or heating when required; in this way, it is avoided by introducing warm air on the warmest day and cold air on the coldest day. Due to restrictions of the numerical code, it is suitable for laminar and transitional flow only, but if it is required to use high air velocities or turbulent regime to achieve a high air change ratio (ACR), it could not provide reliable results. Thus, it was decided to increase the number of pipes to determine the amount of ACR necessary to generate a significant air change ratio within the room. Implementing 20 pipes in an area of 8.55 m 2 to obtain an air change value that generates a significant effect within the room could be somewhat difficult to carry out; however, within a real system, it is possible to obtain the same maximum value of ACR reported in the study, modifying specific EAHE parameters such as air velocity and pipe diameter to increase volumetric flow, which results in fewer pipes for the same ACR value. On the other hand, according to the research developed by Minaei and Rabani [65], it is considered that a separation distance greater than 15 cm does not significantly affect the temperature at the outlet of the EAHE, so if someone wants to implement 20 pipes in the room, 10 pipes can be grouped on the south wall and 10 pipes on the east wall with a separation distance of 15 cm, covering a length of 2.4 m on each side. Another factor that is worth mentioning is the soil thermal saturation. When the EAHE is working, the temperature of the soil surrounding the pipes could increase as time goes on, then after a period between 6 and 12 h of continuous operation, the soil becomes thermally saturated and it could affect the performance of the EAHE [9]. The soil thermal saturation was not addressed in this work, but it should be considered when intermittent EAHE use is recommended. Table 8 presents the maximum and minimum indoor air temperatures of the building room without EAHE obtained with TRNSYS and both temperatures when the EAHE is connected to the building room obtained with a CFD-BES coupling. The table presents results for airflow with Re = 4000, and it is supposed that twenty pipes are connected to the building room for all cities. Further, the temperature difference (∆T) between the indoor air temperature of the room without EAHE and with EAHE is also presented. A high Reynolds value was chosen because a high velocity of the airflow supplied by the EAHE favors its mixing with the room's indoor air.
On the warmest day, the values of ∆T for the different cities range between 1.7 and 3.2°C. The indoor air temperature reduction provided by the EAHE appears to be small, but it is worth mentioning that the room is closed and does not have other infiltration. Although infiltration is a relevant parameter, especially for having a healthy indoor environment, here, infiltration, other than the air introduced by the EAHE into the room, was not considered, because this research aims to focus on the influence of the EAHE on the indoor air temperature of the room. Thus, the indoor air temperature of the room without EAHE and the corresponding temperature of a room with EAHE were compared. Thus, leaving out other parameters could give a better understanding of the behavior of the air in the room with the two conditions mentioned above. In Monterrey, the EAHE was able to reduce the indoor air temperature by 3.2°C, which is the maximum temperature reduction of all cities. The second more significant temperature reduction occurred in Juárez City with a value of 2.8°C. Then, Zacaltipán had a temperature reduction of 2.6°C. Finally, in Mexico City, the EAHE reduced the indoor air temperature by up to 1.7°C. Thus, on warm days, buildings situated in cities in Mexico benefit from installing EAHE because the indoor air temperature is decreased.
On the other hand, Table 8 also shows the values of ∆T for the coldest day of the different cities. The values of ∆T ranged between 1.2 and 3.5°C. In Juárez City, where indoor air reached the lowest value, the EAHE increased indoor air temperature by 3.5°C. The second more significant temperature increase occurred in Mexico City with a value of 2.3°C. Finally, the indoor air temperature increase provided by the EAHE in Zacaltipán (1.2°C) was very similar to that in Monterrey (1.3°C). The EAHE connected to a building room can influence the thermal load required to condition the room space. Because the most significant effect on temperature is obtained with the airflow supplied by 20 pipes for the different Reynolds numbers, the cooling and heating loads of the room with EAHE with 20 pipes are presented here. Table 9 shows the cooling load of the warmest day for each city under study. The values of the cooling load indicate that the thermal load decreases when airflow increases with velocity. The most significant cooling load reduction occurred in Zacualtipán, 2.39 kWh, which is equal to a 6% reduction. The second more significant cooling load reduction was obtained for the room in Monterrey and Juárez City, where the cooling load reductions were similar, with values of 1.64 kWh and 1.61 kWh, respectively. Both values also represent a 2% cooling load reduction. Finally, a minor decrease occurred in Mexico City, with 0.86 kWh (3%). The cooling load reductions are relatively small, but it should be remembered that these values belong to just one day. Depending on the number of days of the year that need cooling in each type of climate, the EAHE could reduce significantly the annual cooling load. Although the previous reductions are small, the year-round savings could be important depending on the type of climate in which the EAHE is installed. Future research should consider the year-round performance of EAHE connected to the building room.  Table 10 presents the heating loads of the building room on the coldest day of the year 2020. The table demonstrates an EAHE with 20 pipes, and Re = 4000 provides the most significant room heating load reduction. When the EAHE is installed in the building room in Mexico City, it provided the most significant heating load reduction; it was reduced by 3.51 kWh (11%). In Zacualtipán, the EAHE reduced the heating load by 2.11 kWh, equivalent to 5%. In Juárez City, the heating load was reduced by 0.49 kWh (2%). In Monterrey, the EAHE reduced the heating load by just 0.16 kWh (0.3%). As in the cooling loads, the heating load reductions might be small, but these values belong to just one day. Depending on the number of days of the year that need heating, the EAHE could reduce significantly the annual heating load.

Conclusions
A simulation method that couples BES and CFD was used to simulate the thermal behavior of an EAHE connected to a building room in four cities with different climates in Mexico. Simulations were performed for the warmest and the coldest days of 2020, and the following is concluded: • The comparison of the indoor air temperature and the temperature of the air at the outlet of the EAHE allowed us to understand the hours in which the EAHE works well and in which it does not work for each of the four cities. On the warmest day, the EAHE could supply air with a lower temperature than indoor air of the room for 24 h in Monterrey, Juárez City, and Zacualtipán. The EAHE worked well for 18 h in Mexico City. • On the coldest day, Juárez City, Monterrey, Mexico City, and Zacualtipán require heating, and the EAHE could supply air with a higher temperature than indoor air for 24 h, 17 h, 10 h, and 8 h, respectively. • It was found that the indoor air temperature conditions improved due to installing the EAHE with 20 pipes. On the warmest day, the EAHE decreased the indoor air temperature by a factor between 1.7 and 3.2°C. On the coldest day, the EAHE increased the indoor air temperature between 1.5 and 3.5°C. • The effect of connecting an EAHE on the daily heating and cooling loads of the building room was also analyzed. On the warmest day, the EAHE reduced the cooling load of the room between 2% and 6%. On the coldest day, the EAHE reduced the heating load by a factor between 0.3% and 11%. Although the previous reductions are small, the year-round savings could be important depending on the type of climate in which the EAHE is installed.
Future research should consider the year-round performance of EAHE connected to the building room in the different types of climates in Mexico. Further, the payback period of the EAHE should be found in order to determine in which locations it is feasible to install EAHE for improving indoor environment conditions. Finally, the soil thermal saturation was not addressed in this work, but it should be considered when intermittent EAHE use is recommended. memory of our mentor and friend Jesús Xamán, who guided us to develop this work and encouraged everybody when we needed it.

Conflicts of Interest:
The authors declare no conflicts of interest.