Integration of evaporative cooling technique with solar chimney to improve indoor thermal environment in the New Assiut City, Egypt

Cooling buildings in summer is one of the main environmental problems for architects and occupants in many hot dry countries. The summer temperature during these countries reaches peaks of more than 40°C in some. Mechanical air conditioners can solve the problem, but they put a heavy strain on the electricity consumption. Egypt in general has rich sunny and clear skies. Therefore, these conditions encourage to enhance evaporating with natural ventilation and save energy. This paper develops an integration of direct evaporative cooling tower with a solar chimney multi-zone thermal ventilation model. Simulation is done using commercial couple multi-zone airflow under COMIS-TRNSYS software (Madison, WI, USA) to assess natural ventilation and indoor thermal comfort. The results show that the system generates 130.5 m3/h under the effect of solar radiation only and minimum 2 ACH without pressure coefficient which is considered the minimum requirement of ACH. The findings show that the new integrated system interacts with the building envelope and weather conditions to achieve a decrease in indoor temperatures that reach 10°C to 11.5°C compared to outdoor temperatures.


Background
Buildings contribute over 40% of the total global primary energy use corresponding to 24% of the world CO 2 . Heating, ventilation, and air conditioning (HVAC) system is the technology of indoor and vehicular environmental comfort which is responsible for about half of the energy use in buildings [1]. The quest to reduce environmental impacts of conventional energy resources and, more importantly, to meet the growing energy demand of the global population had motivated considerable research attention in a wide range of environmental and engineering application [2]. Effective integration of passive features into the building design can significantly minimize the air conditioning demand while maintaining thermal comfort [3]. Passive evaporative cooling is one of the most efficient and long recognized ways of inducing thermal comfort in predominantly hot dry climates. Historically, evaporative cooling was used extensively in traditional architecture throughout the world's hot arid countries [4][5][6]. In order to extend the use of evaporative cooling, short wind tower cooling was integrated with inclined solar chimney to increase the pressure difference and the system efficiency. This was done after observing the indoor discomfort for one of the Egyptian houses built in the New Assiut City [7].
Many researches have been conducted on using natural ventilation and evaporative cooling strategies for producing cool air and the effect of using a solar chimney on thermal-induced ventilation in buildings. Three important researches were under taken to study these passive systems. Bahadori introduced the idea of capturing the wind within a tower and then pass it through wetted conduit walls [4]. Maerefat and Haghighi also put forward a new solar system employing a solar chimney together with an evaporative cooling cavity. The numerical calculation showed that this integrated system with the proper configurations was capable of providing good indoor conditions during daytime in the living room even at a poor solar intensity of 200 W/m 2 and a high ambient air temperature of 40°C [8]. Alemu et al. developed an integrated model incorporating passive airflow components into a coupled multi-zone ventilation and building thermal model. This model allows an assessment of a combination of passive features such as solar chimney and wind-induced earth-air tunnel for both natural and hybrid ventilation systems at the design stage [9]. A lot of awareness is ongoing worldwide on solar energy utilization [10].
The main objective of this paper is to develop an integrated solar chimney with passive cooling strategies into a coupled multi-zone ventilation and building thermal model. This was done by developing a mathematical model to incorporate passive features into a multi-zone ventilation model incorporating a simple thermal model using COMIS-TRNSYS software (Madison, WI, USA). This development system is used in a room of a single zone to study the performances and the advantages of an integrated system in the numerical model. Also, it focuses on the simulated performances and its impact on the comfort of occupants for the summer period with an integrated proposed system.
A theoretical study of the direct evaporative wet medium An evaporative cooler (EC) is a device that cools air by passing it across a wet medium. As hot air passes over the wet medium, some of its energy goes to evaporate water from the medium. The air comes out of the other side of the device having given up some of its energy and consequently exits at a lower temperature with higher humidity and constant enthalpy. The minimum temperature that can be reached is the thermodynamic wet bulb temperature (TWBT) of the incoming air. It was reported that half of the conventional cooling device in the south west of the USA was replaced by simple direct evaporative air coolers; this would save 18 million barrels of oil per year [5,11].
In this paper, heat and mass transfer from the wet medium into the airflow are taken into account. Based on the theoretical study and initial system conditions, evaporative cooler performance is based on the concept of the adiabatic process. This means no heat is added to or removed from the system; the process of exchanging the sensible heat of the air for latent heat of evaporation from the water is an adiabatic process [12].

A description of the system
The performance of the system is to provide desired comfortable conditions and a suitable rate of ventilation, depending on several parameters such as the ambient conditions (temperature, solar radiation, relative humidity, wind speed, and pressure coefficient). The following dimensions and specifications are applied to the model: The system is located in the New Assiut City, Egypt. New Assiut City is located at 27.30°N latitude and 31.15°E longitude. The solar chimney is oriented towards the south. The calculations were carried out for a single zone and to study the integration of the new proposed system, having a dimension of 4.0 m × 4.0 m × 3.125 m. This system was developed from the reference system of inclined solar chimney attached to a room with an evaporative cooling cavity (ECC) [8]. The old model assumes zero pressure boundaries at the inlet and outlet except for the stack effect. Therefore, it is not taken under consideration the variety of the buildings' static pressure from the atmospheric pressure and the effect of wind pressure coefficient on buildings. Also, there are many reasons that make it very difficult to integrate this system into Egyptian housing. The ECC strategies consumed a lot of water because of water film and a lot of space. Figure 1 shows the diagram of the reference system and the new proposed system.
Nevertheless with the new proposed system, the inclination length of the chimney is 2 m on the inclined angle in order not to extend the vertical height (1.5 m) according to Egyptian Building Regulation Law, and the height of the cooling tower is 1 m. The chimney air gap is 0.2 m, and the inclination of the chimney is 50°based on the angle of the reference system. Table 1 shows the parameters of the solar chimney. Figure 2 shows the location of the room and its dimension in relation to a northern orientation. The performance of the system was monitored at a steady state of condition and then 5 days in the month of June (19, 20, 21, and 23) were selected for studying the performance of the solar chimney between sunrise and sunset in order to study the timewise to focus on the hottest day in the summer season (20 June). Then, parametric studies for solar chimney and wind tower are needed in the future studies in order to choose the best dimensions for a new proposed compact design. Some parameters used in the calculations of solar chimney appear in Table 1.
The evaporative pad provides a large water surface, and the pad is wetted by dripping water from the above source. The EC works by a concentric float valve, which opens when the water level is low in the collecting grill, allowing more water to enter. When the water level returns to the full level, the valve is shut automatically.

Methods
The system is being developed with a solar chimney and a small evaporative cooling wind tower. The experimentation would become time consuming and very expensive. Therefore, the numerical model is made and developed for the system. The developed model is implemented in the COMIS-TRNSYS simulation software. This direct evaporative cooling tower used a component 506d-TESS library which uses a wet medium at the top, and this component was assembled and validated by the Thermal Energy System Specialists, LLC (Madison, WI, USA) and was modified in 2004. The present model includes special airflow component which needs simultaneous prediction of temperature and airflow rates. In the developed multizone ventilation model, a building is idealized as a system of zones, openings, and ducts linked together by discrete airflow paths. Zones are represented by nodes. Chimney and evaporative cooling channels are represented by ducts for resistance and pressure drop calculations. A hydrostatic condition is assumed in zones, and the flow rate in each link is defined as a function of zone pressure, which results in a system of nonlinear equation solver, defined by the mass conservation for each zone [13]. Figure 3 shows the flow chart solving procedure of TRNSYS-COMIS.
The equations are solved numerically to predict zone pressures, and zone pressures are back substituted in the flow equations to predict the airflow rate at each link [13]. The detailed equations used for large openings and ducts are presented with the fundamentals of the multizone airflow model (COMIS, [14,15]. A description of the mathematical models used in each component can be found in the TRNSYS 16 manual handbook. Multi-zone airflow network models deal with the complexity of flows in a building by recognizing the effects of internal flow restrictions. They require extensive information about flow characteristics and pressure distributions [16]. The pressures are calculated in the previous steps as an initial guess for pressures in the next step. Initial zone densities are based on the pressures found in the last step. The solver updates zone densities as it updates zone pressures [13]. Figure 4 shows the flow chart solving procedure for COMIS-TRNSYS.

Pressure calculation and mass flow rate through openings
To estimate the amount of air flow through an opening, it is necessary to know the pressure difference across the opening and its effective flow area. The pressure at an opening can be due to wind, as well as buoyancy, and it is therefore determined by the location of the opening, as well as the internal and external environmental parameters such as wind velocity, pressure coefficient, and outdoor temperatures.
As solar radiation passes through the glazing and absorbed at the wall surface. The air in the chimney is then heated by convection and radiation from the absorber. The decrease in density experienced by the air causes it to rise, whereupon it is then replaced by air from below, i.e., from the attached room. The rate at which air is drawn through the room depends upon the buoyancy force experienced (i.e., depending upon the temperature differential), the resistance of flow through the chimney, and the resistance to the entry of fresh air from the room. While in a wind tower, the opening faces north towards the preferable wind. This compact design with dimensions of 1 m (width) × 0.7 m (depth) and a height of 1 m depends on solar radiation and solar chimney that sucks air from the room and increases the air flow inside the room.  The ventilation rate inside the system is calculated according to six variables (P 2 , P 3 , P 4 , P 5 , P 6 , and P 7 ) which are solved from Equations 1 to 6, taking into account the base pressure on the ground (P = 0) as presented in Figure 3.
where P 1 = (−ρgh 2 + 0.5ρ 1 Cp north ν 2 ) is the pressure at point 1 N/m 2 , α is the discharge coefficient of the opening in the wind tower, and A 1 is the area of the opening (1 m 2 ).
where the local pressure losses in the inlet ζ WetMediumInlet of the evaporative cooling (wet medium), the friction of Figure 2 The location of the room with the integrated system according to north orientation. duct ζ DuctFriction , and local losses of the duct outlet ζ DuctOutlet can be calculated according to the dimensions of dynamic losses due to duct fitting, Re, roughness, duct dimension, and λ. It can be determined in details from COMIS [13]. Also, the indoor velocity equals V A ¼ Q A (m/s).
where the local pressure losses ζ ChimneyInlet in the inlet of a solar chimney, the friction of chimney ζ ChimneyFriction , and local losses of the chimney outlet ζ ChimneyOutlet can be calculated in a similar approach as used in evaporative cooling. The solar chimney is treated in the COMIS program as a rectangular duct connects to the room from one side and a node with volume that equals to zero from the other side to transfer the calculated temperature of the chimney from TRNSYS and then another rectangular duct to outside.
where P7 = − ρg(R r L cos θ) + 0.5ρ 1 Cp south ν 2 , and ρ ¼ 353 θþ273:15 is the air density (kg/m 3 ). Then, the pressure differences across the opening due to buoyancy and wind are required and solved with each others in order to predict the pressure of the room and ventilation rate Q (m 3 /s).

Temperature prediction in the solar chimney
The average air density in the chimney channel ρ 6 is predicted using the thermal resistance network model in TRNbuild-TRNSYS, based on the heat balances of each element of solar chimney. The model used by Ong [17] ang Bassiouny and Koura [18] modified by Alemu et al. [9] is adopted with some assumptions in order to solve the mathematical model. Flow through the chimney was considered in a steady state. The chimney inlet area and exit areas are equal. The thermal network for the physical model considered is shown in Figure 5. The heat balance equations from the thermal network at each element of the solar chimney are described in Equations 7, 8, and 9:

COMIS TRNSYS
The energy balance on the air flowing through the chimney q″ with a length L and width of the air channel: where m • is the mass flow rate of the air component. The mean air temperature T f is given by where γ (mean temperature weighting factor) = 0.75 [19].
Consider the air inlet to the chimney that has a temperature equal to the room average temperature, substitute in the above equation of q″ WL : The heat transfer coefficient from the glass to atmospheric air The radiation heat transfer coefficient between the glass and the sky is where σ = solar absorptance, ε = emissivity, g = gravitational acceleration (9.81), T sky ¼ 0:0552 T 1:5 a : Figure 5 Thermal network for solar chimney with air flow.
The convective heat transfer outside the glass due to the wind is The radiation heat transfer coefficient between wall and glass h rwg is The heat transfer coefficient between the absorber and outside air is Where Δw is the thickness of insulation, k insulation is the thermal conductivity The convective heat transfer coefficient between air and absorber plate or glass (h air.absorber , h air.glass ) where x refers to absorber plate or glass, and d is the where Ra air.x is the Rayleigh number for constant heat flux, Dl is the length. The solar energy absorbed by the glass is The solar energy absorbed by the absorber substituted the above equations into Equations 7, 8, and 9 and the three equations solved simultaneously to give the T g(glass cover) , T f(Air) , T w(absorber) . The mean air temperature T (Air) in the solar chimney channel is used to determine the average air density ρ B which is then substituted in Equation 5; then, the flow rate in the solar chimney is a function of the chimney temperature.

Temperature predictions in an evaporative cooling wind tower
An evaporative cooling wind tower (ECWT) is used to let the hot dry air to enter from the outside and is cooled when it passes over a wet medium. When the air passes through the evaporative cooler, the component calls the TRNSYS psychrometrics to fully determine the state of the inlet air. The psychrometrics returns the inlet air wet bulb temperature and enthalpy backs the inlet temperature, pressure, relative humidity, and humidity ratio. The performance of the evaporative cooler is based on the concept of adiabatic process. The wet bulb temperature of the air entering and exiting the device is assumed to be constant, where the minimum outlet temperature equals the inlet wet bulb temperature. Figure 6 shows evaporative cooler schematic.
The energy balance equation expressing the heat removed from the air to evaporate the water is expressed as The rate of heat transfer from air to water in the wet panel surface, Q sensible (the heat exchange by temperature difference), is given by Δy and the heat gain of latent heat of air (heat exchange by humidity ratio) can be expressed as where w is the humidity ratio. ΔT LM is the log mean temperature difference which is given by [12,21] The performance of the systems with different flow type can be assessed on the saturation efficiency as The air temperature leaving the cooling device can be calculated from (26) The wet bulb temperature of the air entering and exiting the device is assumed to be constant, where the minimum outlet temperature equals the inlet wet bulb temperature.
Since the performance in the evaporative cooler is an adiabatic process, The humidity ratio of the air outlet from the device [22] RH out ¼ X inlet p outside 0:622 þ X inlet ð Þ p saturated:vapour:pressure : ð29Þ The saturate vapor pressure (atm). Bouchahm, 2011 [23] is The corresponding energy removed from the air stream to evaporate water under these circumstances is given by Equation 31: where h enthalpy = C p T + x(r + C P.vapor T), expressed in kJ/ kg (dry) [24]: where C p.air = the specific heat of air (kj/kg.K), x = humidity ratio (kg/kg (dry)), r = evaporative latent heat,

Air flow network and building thermal model
The ventilation model is integrated to the conventional building thermal model in order to have strong influence on the system performance. The internal heat load, solar load, and the heat transfer through walls, windows, roof, ground, and opening are considered and shown in Figure 7a, where Q is the heat transfer. The temperature of the ground at depth of 0.5 m is used in calculation for heat transfer from the ground and derived from Egyptian Typical Meteorological Year (ETMY) which is developed for standards development and energy simulation by US National Climatic Data Center (CA, USA) for a period of 12 years from 1991 until 2003 [25]. Also, the internal load is used in the calculation, where there are four occupants inside the single zone doing light work and seating with light load. Each occupant is assumed to have a sensible gain of 75 W and latent gain of 75 W according to ISO 7730 [26]. The occupants are assumed to be staying in the room for the whole day to apply the maximum load inside. The room light has 13 W/m 2 . It has an energy saving lamp with convective part (visible) of 20%, and one device of 140 W is used. The zones (ECWT, room, and SC) are linked to their ambiance via external nodes and to the north and south facades. Wind pressure data, wind speed, wind velocity profile for the building and for meteo, and orientation are defined. The COMIS program uses wind speed at meteo site in reference height to calculate the speed above the boundary layer using the wind velocity profile of the meteo station. The boundary layer is 60 m high for smooth terrain. Wind pressure coefficients (Cp values) belonging to a certain façade element can be given for several wind directions derived from 'CP-Generator' [13,27] according to the room with the integrated system with maximum height of 4.5 m and around the building. Figure 7b shows the air flow network in COMIS. Table 2 shows the building construction material description based on the ideal thermophysical properties [28] which was used in TRNBuild. The building thermal model is calculated based on the mathematical Equations 33 to 38.
The heat transfer through external walls (and roof ) for room zone and tower zone is predicted as where ne is the number of external walls attached to the zone, T star is star temperature which can be used to calculate a net radiative and convection heat flux from the wall surface, R is the overall thermal resistance value including convective heat transfer resistance on both sides of the surface, and z is the zone [26].
The total gain to chimney zone from all surfaces is the sum of the combined heat transfers: where A s is the inside area of the surface, and S outside surface is the solar radiation for outside surface. The heat transfer across the windows is where U w is the overall heat transfer coefficient across the window (the rate of heat transfer through the window from outside to inside when it is hot per unit area and per unit temperature difference (W/m 2 C), SHGF is the window's solar heat gain factor and its calculation is presented in details [29]. For external surfaces, the long wave radiation is considered using a fictive sky temperature T sky [23]: where f sky is the view factor to the sky, and ε is the long wave emissivity of outside surface = 0.9 for walls.  The sky model used to calculate the radiation on inclined surface for the solar chimney by Perez model, which is considered the best available model.
Heat transfer due to ventilation in a zone is X Then, the total heat load to maintain the zone at specified temperature can be calculated as Then, the zone temperature is iteratively predicted in the calculation.

Steady state condition
The performance of the integrated system has been numerically studied, and the effect of the pressure coefficient on ACH and system performance was studied during the steady state condition. The investigation was at constant wind speed at 4 m/s (at Meteo), wind direction of 360°, and pressure coefficient (C p ) 0.081 for the north façade, and −0.040 for the south façade as shown in Table 3. Also, the effect of the different combinations of the three values of solar radiation, temperature, and relative humidity on indoor environment was studied. The main results show the increase of ACH with an increase of solar radiation. The system is capable of generating 130.5 m 3 /h ventilation rates for a collector area of Figure 9 The variation of the ACH with solar radiation with different combinations. Figure 10 The performance of the system during 3 months of summer using standard weather data. It can be seen from Figure 8 that the chimney cavity temperature is always higher than the ambient temperature with a very high temperature for an aluminum surface (black absorber) that reaches to 64°C at the radiation level of 1,000 W/m 2 , then the surface temperature decreases to 40°C at the radiation level of 400 W/m 2 . It has been observed that the black absorber possesses a high temperature value compared to chimney air temperature; this increase in temperature is due to more captured radiation and storage of high thermal energy. This absorbed energy accelerates the flow of air through the chimney during the day and the night. As a result, this exit air flow will help more fresh air to enter to the space from the opening of the wind tower. Figure 9 shows the increase of ACH with an increase of solar radiation and pressure coefficient with maximum 7 ACH inside the room. Also, the minimum indoor air change rate is 2 ACH with 200 W/m 2 of solar radiation and no wind effect. Thus, the indoor ventilation rate reached zero with no radiation and no wind pressure coefficient. The temperature of the air is higher at the inlet and decreases at the outlet in the evaporative cooler (EC). This is because, initially, the dry air has a large potential to absorb the moisture and gets saturated as the airflow reaches the outlet of the ECWT.
According to ASHRAE 62-2001 standard (Atlanta, GA, USA), ventilation rates dependent on the floor area, whereas the minimum ACH is 0.35, but no less than 15 CFM/person. Also, the Egyptian Energy Code for residential buildings implies a minimum natural ventilation rate of 3 L/s-person (10.8 m 3 /h) for the living room [30]. Therefore, the developed system is able to provide the desired ventilation and thermal needs of the occupants with/without solar radiation as shown in Table 3.

Real weather data
The real weather data of the New Assiut City which is taken from ideal weather data of 20 years [26] was applied to study the system performance during the summer season and the hottest days. In terms of thermal comfort, ASHRAE recognizes that the conditions Figure 11 The performance of the system on 19 to 23 June with focus on the hottest day of June.

Figure 12
The humidity environment of the hottest day (20 June) using psychrometric chart of ASHRAE. Blue: outdoor temperature. Red: indoor temperature. required for thermal comfort in spaces that are naturally conditioned are not necessarily the same as those required for other indoor spaces [31]. So, a new ACS for a naturally ventilated building has been proposed to be integrated within ASHRAE Standard 55, 2004. In the ACS, the mean monthly outdoor air temperature determines the acceptable indoor operative temperature. This relationship is expressed by the following formula: To (comf ) = 0.31Ta (out) +17.8, where To (comf ) is the optimum comfort operative temperature in°C and Ta (out) is the mean monthly outdoor air temperature in°C.
Therefore, most of the outlet temperature after passing through the wetted medium in ECWT is below the acceptable limits of Adaptive Comfort Standard (ACS) as shown in Figure 10. Only 5% of the results data exceed the upper limit of 80% acceptable range of ACS especially for the hottest days. An investigation of the system performance during the 5 days of June was done to focus on the hottest day in the summer season from the weather data file so that a clear pattern is seen in relation to adaptive comfort standard (ACS) as shown in Figure 11.
Arundel et al. proposed that the low or high humidity strongly affects the health and the optimum humidity level which is between 40% and 60% [32]. Therefore, the humidity environment in the living room and the bedroom is very important. Figure 12 shows most of the simulation results for the hottest day, which are located between 40% and 60% except before sunrise between 1 am and 6 am. Thus, this shows the effectiveness of a natural ventilation strategy with a direct evaporative cooling for an indoor environment. Finally, this proposed compact system could be an economical and practical passive alternative to the conventional air conditioning systems in hot and dry climates with no pollution released or energy consumed.

Conclusions
Important conclusions are drawn from the numerical operations of a new integrated system. Pressure coefficient and solar radiation in relation to the area of the collector are the main parameters that affect the air flow rate inside a single zone. When the temperature of the black absorber increases with/without pressure coefficient effects, the ACH increases with a strong relation to solar radiation as shown in Figure 8. The temperature of the indoor environment decreases by 10°C to 11.5°C compared to the outdoor temperature. This helps to achieve comfort during the hottest days of the summer season, with 95% of indoor temperature below the upper range of 80% acceptable comfort range during the period from 21 May to 21 August. The proposed system can be applied during daytime and nighttime, but it can be controlled and its use can be limited during nighttime especially when the inlet outdoor temperature is ≤30.5°C.
Overall, the results show how passive systems interact with the building and weather conditions. This paper studied the effect of a new integrated system for indoor comfort for the New Assiut City as a case study. The authors' future work will be a parametric investigation and optimization of the important operating parameter to achieve compact and high performance design so that guidelines for an efficient optimization process will be done to optimize cases with different locations, dimensions, and number of occupants.

Nomenclature
The list of nomenclature is presented in Table 4.