Abstract

Thermal storage capacity and airflow rate of a solar chimney combined with different PCMs are numerically studied during nighttime. PCMs with phase change temperatures of 38°C, 44°C, 50°C, and 63°C are selected in this numerical study. Results show that the maximum average ventilation rate of 610 kg/m2 and maximum thermal storage of 4750 kJ/m2 are achieved at the phase change temperature of 38°C. However, for phase change temperature of 63°C, night ventilation does not occur under the identical conditions. The findings reveal that a lower phase change temperature can increase the chargeability (and therefore the dischargeability) of a solar chimney, since a higher phase change temperature demands higher solar radiation intensity and longer charging time for a solar chimney. For PCM with a phase change temperature of 44°C, most of the heat stored in PCM is lost to ambient through glass cover by radiation and only a small portion is used for heating the air within air channel.

1. Introduction

A solar chimney is a natural draft system that has already been applied in the building ventilation widely and has attracted the interest of many scholars around the world. Extensive numerical or/and experimental studies in terms of evaluating ventilation performance [15] and optimizing design of solar chimney [610] have been performed. A solar chimney consisted of roof solar collector, and Trombe wall was experimentally investigated by Khedari et al. [11]. Results show that in June~July of a hot and humid climate, the temperature, air velocity, and flow rate per area of a solar chimney inside the room varies from 35–37°C, 0.02–0.08 m/s, and 0.01–0.02 m3/s per 1 m2, respectively. Hirunlabh et al. [12] performed a numerical study on four new configurations of roof solar collector under Thailand weather conditions. By using the new configuration of roof solar collector, the highest volume flow rate of air is 0.072 m3/s or 0.0206 m3/s per 1 m2 of a solar chimney. Chungloo and Limmeechokchai found that with a solar chimney, the indoor temperature can be reduced by 1.0–3.5°C at high ambient temperature and high solar intensity (32.0–40.01°C) in the daytime [13]. Miyazaki et al. [14] numerically investigated the effects of solar chimneys on thermal load mitigation of office buildings in Japan. The results showed that in the natural ventilation mode, a larger chimney area was required to reduce the passive cooling load of the building. Between 10 a.m. and 12 p.m. of an average day in May, when the chimney width is 4 m, the cooling load was less than that of the no solar chimney design.

Studies reveal that the height [6], chimney position [7], type of absorber [8], width and depth of cavity [15], inclination angle [16], and the insulation or thermal mass in the solar chimney [17] are the primary factors that should be considered when designing a solar chimney. The above analysis indicated that solar chimney possesses the potential of cooling effect and cooling load reduction, and the potential can be enhanced by optimizing the designing factors. The available studies mainly concentrated on the performance characteristics of solar chimney during the daytime. If the night ventilation or all-day ventilation can be achieved, the solar energy utilization as well as the indoor thermal comfort would be greatly improved. Amori and Mohammed [18] experimentally investigated the effect of integrating phase change material (PCM) in a solar chimney and found that integrating a solar chimney with PCM yield longer ventilation period after the sunset. In this paper, a PCM unit is integrated into the solar absorber of the chimney to take advantage of PCM’s isothermal phase change characteristic as well as high latent heat density. The PCM absorbs the heat radiated from the sun and stores heat during the daytime and then releases the absorbed heat during the nighttime. By this means, night ventilation can be achieved without consuming any additional energy sources. This has been indicated by Liu and Li’s study [19], in which the performance of a conventional vertical solar chimney and a PCM-based one with the same geometry was compared. In their study, the paraffin wax with phase change temperature of 38–42°C was selected as PCM. It was found that the effective ventilation time of the PCM-based solar chimney can be as long as 13 h 50 min, while the conventional one’ ventilation time is only 1 h 20 min. Similarly, Zhou and Pang [20] experimentally studied the thermal behavior of a system with a collector-storage wall using PCM with phase change temperature 26 ± 1°C. The experimental results indicated that by integrating PCM in the wall, the whole ventilation period of the system is as long as 17.5 h.

It is should be noticed that the PCM can significantly influence the ventilation performance of the PCM-based solar chimney, especially the phase change temperature of the PCM, which is a vital factor among the physical properties of the PCM. It directly affects the absorber temperature and consequently affects the air flow rate and thermal storage capacity.

Based on the findings from above literatures, a mathematical model is developed to examine the ventilation performance and thermal storage capacity of a solar chimney using different PCMs under hot summer climatic conditions. The results of this paper could help to enhance the performance of such a PCM-based solar chimney by optimizing the phase change temperature.

2. Description of Mathematical Model

The schematic diagram of the PCM-based solar chimney is shown in Figures 1 and 2. The system mainly consists of glazing cover, air cavity, heat absorber plate, PCM container, and extruded insulation panels.

2.1. Calculation of the Energy Stored during the Day

The inlet and outlet of this system are closed when PCM is storing the solar energy during the day. Solar radiation passes through the glazing cover and is absorbed by the absorber. Then, the absorbed energy is used to raise the temperature of PCM. When the PCM temperature reaches its phase change temperature, the PCM starts to melt and stores latent heat simultaneously. During this period, the temperature of the PCM remains relatively constant.

Applying the energy balance concept on the glazing cover and the absorber plate, the following equations are given:

Theoretically, the glass receives heat from solar radiation and radiative heat from the absorber plate while it loses energy to ambient. The absorber plate, made by a thin stainless steel with high thermal conductivity, is used to absorb solar energy. Hence, it is reasonably assumed that there is no thermal gradient across the thickness of the absorber plate. Consequently, the absorber plate absorbs solar radiation through the glass cover while releases heat to the glass cover and the PCM.

In this paper, the PCM is sealed in a metallic container, which is assumed to be in good thermal contact with the PCM. And the absorber plate serves as a heat exchanger between the air and the container. Hence, the control volume method for discrete one-dimensional phase change problems is adopted. The PCM container is divided into () cells along the thickness direction. The one-dimensional energy equation for each internal node “” can be written as follows:

The heat transfer between different components of this system is shown in Figure 1, where qs is the solar radiation intensity. (1)The thermal storage capacity of the PCM(2)The radiative heat transfer coefficient between the absorber and glass cover has been obtained from [21](3)The overall heat loss coefficient from the glass cover to ambient includes the convection caused by wind and radiative heat transfer from glass cover to sky. It can be written as

The radiative heat transfer coefficient from the glass cover to the sky may be obtained from [21]

The sky temperature and the convective heat transfer coefficient between the glazing cover and the surrounding air, which is affected by wind speed, can be calculated by [21, 22] (4)The overall heat transfer coefficient between the PCM and the absorber plate is given by(5)The heat that the PCM transfers to indoor is given by

2.2. Calculation of the Air Flow Rate during the Night

The discharging period begins when the chimney starts ventilating during the night. The air exchanges heat with the absorber plate as it passes through the air channel and then flows to ambient through the outlet. For various components of the system, the energy balance equations for the discharging period are presented below.

Applying the energy balance concept on the glass cover, absorber plate, and air flow in the air channel, the following equations are yielded:

Assuming that the natural convection of the air occurs between the glazing cover and the absorber plate, the convective heat transfer coefficient hcv,a-g between the glazing cover and the air in the channel can be calculated as follows [23]: where is the inclination angle of the absorber plate ( in this study) and the Grashof number is

It is assumed that the mean air temperature inside the chimney is equal to the average value of the absorber temperature and the glazing cover temperature at the absorber side.

The relationship between inlet air temperature and outlet air temperature can be described as follows [24, 25]: where , the mean temperature approximation coefficient, was found in the literature to be 0.74 [21, 22] and the inlet air temperature can be regarded as the outdoor air temperature.

The airflow rate as a result of the buoyancy effect can be calculated as follows [26]: where is found in the literature to be 0.57 [26].

During the daytime, the initial state of all the PCM is in its solid phase and the whole system’s temperature is assumed to be the same as the ambient temperature. For the night ventilation period, the temperature of each system component and the PCM temperature depend on the final charging conditions. The inlet air temperature is equal to the ambient temperature. The temperatures of , , , and as well as the airflow rate are obtained with a program solved in “MATLAB.” The air outlet temperature can be obtained with the mean air temperature and (14). With the developed code, the performance of the PCM-based solar chimney as well as the heat transfer between various components of the solar chimney has been investigated.

3. Analysis and Discussion of the Numerical Results

The input parameters of the proposed system are given in Tables 1 and 2. The ventilation time is scheduled for 19:00–7:00. The monthly average solar radiation at the horizontal in Chongqing, China is close to the total solar radiation of June, and the solar radiation on June 26 approaches to the daily average solar radiation of June. Therefore, the climatic parameters of June 26 in Chongqing obtained from [27] are applied in this study, and the solar radiation intensity is given in Figure 2.

This paper analyzes the effect of phase change temperature on the ventilation rate of a solar chimney under Chongqing climatic conditions. The PCMs used in this numerical study are myristoyl, dodecylic acid, myristelaidic acid, and palmitic acid, and their corresponding phase change temperatures are 38°C, 44°C, 50°C, and 63°C, respectively.

Figure 3 shows the stored energy and ventilation rate for different PCMs. It is seen that the higher the phase change temperature, the lower the energy storage and the ventilation rate per unit area of absorber. That is to say, the PCM with phase change temperature of 38°C has the best performance among the studied PCMs, while the PCM with phase change temperature of 63°C has the poorest performance. Specifically, for the PCM with phase change temperature of 38°C, the solar chimney achieves the maximum energy storage of 4750 kJ/m2 and the maximum ventilation rate of 610 kg/m2. However, for the PCM with phase change temperature of 63°C, there is little energy stored in PCM, and consequently the airflow is not observed during 19:00–07:00.

Figure 4 shows the absorber temperature and the PCM temperature for four different PCMs when storing heat during daytime. As shown, the thermal behavior characteristics of PCM under the identical conditions are different. The temperature variations of PCM display a typical fully charging process for phase change temperatures of 38°C and 44°C, respectively. The melting times for these two PCMs are 5 and 7.5 h, respectively. However, for the phase change temperatures of 50 and 63°C, the maximum PCM temperatures is 49.6°C and 53°C at the end of charging process accordingly, indicating that PCMs are not fully melting. Especially for phase change temperature of 63°C, the PCM temperature even decreases after the peak solar radiation intensity. This is because the absorber surface temperature descends quickly with solar radiation intensity weakens due to the high surface temperature. Consequently, the sensible heat stored in PCM transfers to the absorber surface.

This indicates that a lower phase change temperature can increase the chargeability and dischargeability of the solar chimney. The larger amount of energy stored and discharged by the solar chimney and consequently a higher ventilation rate is achieved. A higher phase change temperature demands higher solar radiation intensity and longer charging time and results in poor chargeability of the chimney.

Figures 57 show the variations of air temperature and the transient ventilation rate with time.

It can be seen from Figures 57 that the changing trend of ventilation rate is similar to that of the air temperature difference between inlet and outlet. Higher phase change temperature leads to poor chargeability of the chimney and has the shorter ventilation time. The ventilation times corresponding to different PCMs are shown in Table 3.

As shown in Table 3, when the phase change temperature of the PCM is 38°C, the longest ventilation time of 12 hours was achieved while the PCM with phase change temperature of 50°C has the shortest ventilation time which is only 4.33 hours. This is because more solar energy was stored in the PCM with lower phase change temperature.

Figures 8 and 9 show the energy storage and release of different components of chimney combined with the dodecylic acid (phase change temperature of 44°C) during the day and night, respectively.

It can be seen from Figure 8 that the available solar energy is as high as 11,000 kJ/m2, while the energy stored in PCM is only 3150 kJ/m2, which accounts for 28.6% of the . The heat transferred by radiation from glass cover to ambient is 5090 kJ/m2 that accounts for 46.3% of . It means that nearly 50% of the solar energy is lost through glass cover to ambient during the daytime period. The solar energy utility efficiency is still low for PCM of dodecylic acid.

It can be seen from Figure 9 that most of energy stored in PCM is transferred to absorber surface, accounts for 89% of . However, 80% of loses to ambient through glass cover during night ventilation. The heat absorbed by air is only about 13% of . This means that most of heat stored in PCM is lost to ambient through glass cover by radiation during the night ventilation; only a small portion is used for heating the air within air channel.

4. Validation of the Numerical Model

The numerical model is validated by comparing numerical predictions of PCM temperature and air mass flow rate with the experimental results obtained from Liu and Li’s study [28]. This experiment investigated the thermal performance of a solar chimney integrated with PCM for a given heat flux. The PCM temperature during charging and discharging processes and air flow rate were investigated. As shown in Figure 10, a rectangular PCM container with dimensions of 1000 mm wide × 1600 mm high × 40 mm deep was constructed of stainless steel. A total of 50 kg RT42 was used to fill the container, and the physical properties of RT42 was given in [28]. In their experimental testing, the openings of inlet and outlet were closed during the charging period. The discharging period happened right after the charging period. A number of K-type thermocouples with the measurement range of −50 to 250°C and accuracy of ±0.3°C were distributed on the glass cover, absorber plate, inside the PCM, and the air channel. Air velocities inside the air channel are measured by a TSI8455 air flow probe with the measurement range of 0–50 m/s and accuracy of ±3% of reading. The setup of the developed model in this study follows the experimental conditions. Figure 11(a) shows the PCM temperatures obtained by numerical simulation and experiment testing. It can be seen that there was a good agreement between the numerical data and experimental data. The coherence indicates the validity of the developed numerical model. Figure 11(b) shows that the predicted air flow rate varies with the measured one within a small divergence. Therefore, it has the capability to predict the mass flow rate accurately.

5. Conclusions

In this present paper, a numerical model is developed for optimizing the phase change temperature for a PCM-based solar chimney by considering the amount of energy stored in PCM and ventilation rate. Based on the numerical results, the following conclusions can be made: (1)The higher the phase change temperature, the lower the energy storage and ventilation rate per unit area of absorber. In this study, the PCM with phase change temperature of 38°C has the best performance among the studied PCMs while the PCM with phase change temperature of 63°C has the poorest performance. Specifically, for the PCM with phase change temperature of 38°C, the solar chimney achieves the maximum energy storage of 4750 kJ/m2 and the maximum ventilation rate of 610 kg/m2. However, for PCM with phase change temperature of 63°C, there is little energy stored in PCM and the ventilation does not occur during the scheduled ventilation time.(2)A lower phase change temperature can increase the chargeability and dischargeability of the solar chimney. The larger amount of energy stored and discharged by the solar chimney and consequently a higher ventilation rate is achieved. A higher phase change temperature demands higher solar radiation intensity and longer charging time, and results in poor chargeability of the chimney.(3)During daytime, the energy stored in PCM with phase change temperature of 44°C, only accounts for 28.6% of the available solar energy, while nearly 50% of the solar energy is lost through glass cover to ambient by radiation during daytime period. The solar energy utility efficiency is low.(4)During night ventilation period, the heat absorbed by air is only about 13% of . This means that most of the heat stored in PCM is lost to ambient through glass cover by radiation; only a small portion is used for heating the air within air channel.(5)It can be concluded that most of the stored energy is lost through glass cover to ambient by radiant heat transfer. Therefore, reducing the radiation heat transfer coefficient of glass cover can increase the night ventilation effectively.

Nomenclature

A:Surface area (m2)
Ain, Aout:Cross sectional area of inlet and outlet (m2)
Ar:Aspect ratio of Ain/Aout
Cd:Coefficient of discharge of air channel inlet
C:Specific heat (J·kg−1 K−1)
h:Convective heat transfer coefficient
Kw1:Thermal conductivity of brick wall (W·m−1 K−1)
Kw2:Thermal conductivity of extruded insulation panels (W·m−1 K−1)
L:Length (m)
m:Mass flow rate (kg·s−1)
Q:Total heat (kJ)
q:Heat transfer (W)
qs:Solar radiation intensity (W·m−2)
T:Temperature (K)
:Overall heat transfer coefficient from PCM to indoor (W·m−2 K−1)
:Overall heat transfer coefficient from glass cover to ambient (W·m−2 K−1)
V:Volume (m3)
v:Velocity (m·s−1).
Greek Symbols
Δw1:Thickness of brick wall
Δw2:Thickness of extruded insulation panels
Δt:Temperature difference (K)
α:Absorptivity
ρ:Density (kg·m−3)
:Gravitational acceleration (m·s−2)
σ:Stefan-Boltzmann constant
ε:Emissivity
β:Inclination angle measured from the horizontal plane
δa:Gap between absorber wall and glass
γ:Constant for mean temperature approximation
λ:Thermal conductivity (W·m−1 K−1)
τ:Glass transmittance
ν:Kinematic viscosity (m2·s−1)
ζ:Emissivity.
Subscripts
amb:Ambient conditions
a:Air in the air flow channel
a-g:From air to glass
b:From PCM to indoor
c:Chimney
cv,w-a:Convective from black wall to air
cv,a-g:Convective from air to glass
in:Inlet
out:Outlet
g:Glass
g-sky:From glass to sky
i:From PCM to indoor
p:PCM
p-w:From PCM to black wall
r:Room
r,g-sky:Radiative from glass to sky
r,w-g:Radiative from black wall to glass
s:Sky
t:Air flow channel of the solar roof
w:Absorber wall
w-a:From absorber wall to air
w-g:From absorber wall to the glass
w-p:From absorber wall to PCM
x:Storage capacity of the PCM.
Dimensionless Terms
Nu:Nusselt number
Pr:Prandtl number
Ra:Rayleigh number
Gr:Grashof number.

Conflicts of Interest

The authors declare no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the support from the National Science Foundation of China (no. 51478058), the 111 Project (no. B13041), and the Fundamental Research Funds for the Central Universities (no. 106112015CDJXY210006).