Next Article in Journal
Bioremoval of Methylene Blue from Aqueous Solutions by Green Algae (Bracteacoccus sp.) Isolated from North Jordan: Optimization, Kinetic, and Isotherm Studies
Previous Article in Journal
Restorative Effects of Park Visiting on Physiology, Psychology, and Society and the Factors Influencing Park Visiting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplified Calculation of Tsol Based on Dynamic Numerical Simulation of Tsky in Diverse Climates in China

1
School of Architecture & Urban Planning, Shenzhen University, Shenzhen 518060, China
2
Shenzhen Center for Human-Oriented Environment and Sustainable Design (CHESD), Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(1), 839; https://doi.org/10.3390/su15010839
Submission received: 16 November 2022 / Revised: 27 December 2022 / Accepted: 28 December 2022 / Published: 3 January 2023

Abstract

:
An accurate calculation of sol-air temperature (Tsol) is very important for urban environments and building energy consumption. There are various methods that can be used to calculate Tsol by considering sky radiation effects. Climate conditions are vital factors affecting sky temperature (Tsky). In this paper, in order to select an appropriate calculation method to determine long-wave radiation, a theoretical analysis was carried out based on the effect of Tsky on the thermal gain of building envelopes due to long-wave radiation. Typical annual meteorological data were selected to calculate Tsol for 10 meteorological stations covering five building thermal zones in China. The application of the Tsol model was studied using MBE as the measurement standard, and a linear regression equation for the calorific value of the envelope obtained via the Tsky estimation method and the Tsky dynamic calculation method was established. The results show that relative humidity is the key meteorological factor that affects the application of the Tsol model and that the Tsky dynamic calculation should be used to calculate long-wave radiation in regions with low relative humidity. A thermal correction equation for buildings was obtained for use in areas lacking meteorological data and to provide a basis for sustainable building design.

1. Introduction

As buildings serve as living and working spaces in various climates, studying the physical laws of the interaction between buildings and natural climates provides a basis for the sustainable development of urban environments [1]. First of all, an evaluation of the cooling energy of buildings needs to consider the heat transport through the building envelope in the transient state, based on the periodic cooling load temperature difference of the walls and roofs. However, the outdoor temperature near a building surface varies significantly throughout the day, depending on various environmental factors and the thermal properties of the surface materials.
To take these parameters into account in the calculation of the environmental temperature and cooling load, the sol-air temperature (Tsol) equation or Stewart’s tabulated total equivalent temperature difference/time averaging (TETD/TA) value was proposed in 1944, and it was later adapted into the ASHRAE Guide and Data Book in 1961 (ASHRAE 1961) [2]. Many existing standards, such as the CIBSE guidelines, define Tsol as a hypothetical external surface temperature index that combines four related variables, namely, short-wave radiation, long-wave radiation exchange, air temperature, and wind speed on the surface of the envelope [3].
Two objects with different temperatures and in close proximity to each other tend to exchange energy via LW radiation (infrared) in order to reach equilibrium. These exchanges happen inside and outside of buildings. It is difficult to study these exchanges, as all surrounding factors must be taken into account, including the sky. These exchanges must be calculated between each and every surface of a building, the temperature and emissivity of which must be known. Moreover, these calculations are complicated by the surface temperature having an exponent of 4. Infrared calculations thus require the thermal calculation of a larger system and the introduction of non-linear terms; therefore, a simplified solution has been proposed [4].
Based on the building location, it is possible to define the climatic area and, consequently, identify a suitable correlation. ASHRAE proposed a simple correlation, assuming a difference between sky and ambient temperatures equal to 6 °C [2]. Considering the heat transfer between a building and the environment, ISO 13790 allows for the calculation of the temperature of the sky. If the sky temperature calculation is not available, the equations listed in Table 1 can be used. ISO 13790 provides a simplified computation for determining the annual energy requirements for heating in residential and non-residential buildings.
To provide guidance for the thermal design of civil buildings, the thermal design code for civil buildings (GB 50176-2016) proposed calculations and specified values [6] Due to the lack of sky radiation data in most parts of China, it is difficult for most architects to master the complex calculation of long-wave radiation. It takes the minimum limit index of thermal design as the calculation target and proposes a T s o l model, as shown in Equation (1), which considers the total solar radiation and the outdoor dry ball temperature, ignoring long-wave radiation and ground short-wave radiation.
Many scholars have studied the influence of meteorological factors on building load. Lauzet’s [7] research shows that, without considering the urban context and climate, 15–89% of heating energy consumption is ignored, and 131–200% of cooling energy consumption is ignored. Jian Hang [8] used scale experiments to compare surface temperature, air temperature, albedo, and SEB components in deep and shallow valleys under clear, partially cloudy, and cloudy sky conditions in a humid subtropical climate. Sky conditions had a significant effect on diurnal variations in the surface temperatures of street canyons. The temperatures of different canyon surfaces peaked at similar times under partially cloudy sky conditions, while they peaked at different times under clear sky conditions. Luca [9] assessed the impact of sky temperature under different climatic conditions on energy demand for building years. It was found that the annual energy demand for heating and cooling was affected significantly, ranging from −10% to + 19% in a tropical climate, from −10% to + 13% in a dry climate, from −19% to + 28% in a mild climate, and finally from −43% to + 83% in snow conditions.
The sensitivity of the sky temperature to building energy requirements and climatic factors suggests that there are large differences between clear sky and cloudy conditions. The clearness index can be regarded as an indicator of cloud cover. Due to the influence of geographical latitude, altitude, the terrain barrier, and other factors in western Inner Mongolia, eastern Xinjiang, the western Qinghai–Tibet Plateau, and other inland areas, the annual amount of water vapor available is very small, and precipitation is scarce, so the total monthly clear sky index remains high throughout the year [10]. Liu studied the sky temperature related to long-wave radiation heat transfer [11]. According to the monthly average data of 82 stations in China from 1960 to 1970, the sky temperature at each station and its differences from the air temperature were calculated. The results of Tsol model calculations that ignore the effect of sky radiation will be biased in areas with strong sky radiation [12].
Parameters such as the reflected radiation intensity and the ground blackness coefficient in the above model need to be accurately measured after the building is completed; therefore, the model is not conducive to thermal design. Many researchers have determined the external surface temperature of the envelope under the influence of long-wave radiation by means of actual measurements and thermal analyses [13,14]. Simulation tools adapted to long-wave radiation calculations have been proposed. Conductive heat gains and losses through building envelopes are transient in nature because of the dynamism associated with hygro-thermal boundary conditions. Simulation analyses have become the main research method used in this field, as well as multi-software comparative analyses of the sensitivity of climate parameters and outdoor environmental differences to long-wave radiation impacting building load [15,16,17].
It can be seen that this link between sky radiation and Tsol is only examined in a few research projects, and there are even fewer dynamic comparisons of Tsol in different climate regions in China.
In this paper, a dynamic calculation method for the sky temperature is introduced by examining the heat gain of a building in a simulation in order to evaluate the impact of long-wave radiation on building heat gain in different climatic regions, and a regression model for long-wave radiation estimation is established to provide a reference for the thermal design of the envelope and a calculation for the energy saving of a building under different climatic conditions, as shown in Figure 1.
To this end, this paper is organized as follows: Section 2 presents the methodology used in this study and illustrates the model used for the simulations. Section 3 compares the calculation results of long-wave radiation in different climatic regions. Section 4 adjusts the numerical solution model. Finally, Tsol is calculated and evaluated in different scenarios.

2. Long-Wave Radiation Calculation and Data Analysis Methods

The Tsol models are divided into two categories according to the level   of influence of the model accuracy on the actual load and the calculation parameters: one is the simplified Tsol model adopted by the current code for the purpose of convenient calculations in actual engineering, and   the other is a long-wave radiation Tsol model with a higher accuracy considering long-wave radiation heat transfer.

2.1. The Simplified Tsol Model

The   thermal design code for civil buildings (GB 50176-2016) [6] proposes an equivalent temperature value that increases the solar radiation relative   to the outdoor air temperature for the convenience of calculation. The calculation formula ignores the ground-reflected radiation and the atmospheric long-wave radiation heat transfer, and   it is shown in Equation (1):
T sol = t e + ρ s I α e
where Tsol is the   sol -air temperature, °C; t e is the outside air temperature, °C; I is the solar irradiance projected on the external surface of the building envelope, W / m 2 ; ρ s is the solar radiation absorption coefficient of the external surface; and   α e is the external surface heat transfer coefficient, W / m 2 · K .

2.2. The Simplified Sky Temperature Calculation

A simplified calculation method for estimating the sky temperature is derived from the radiation and heat balance relationship between the air near the ground and the atmosphere, and   it is defined as Equations (2) and (3):
σ T sky 4 = Q s k y = Q a i r = ε a i r σ T a 4
T sky = ε air 4 T a
where T a is the air temperature at 1.5~2.0 m from the ground, K; ε a i r is the emissivity of the air near the ground, which can be calculated with ε a i r = 0.741 + 0.0062 t d p ; and   t d p is the air dew point temperature near ground, °C.

2.3. The Dynamic Sky Temperature Calculation

For short- term or instantaneously metered long-wave radiant heat, the surface temperature should be used as the calculation parameter. According to the observation data of the ground radiation balance of a meteorological station, the value of Tsky can be calculated. Various formulas have been proposed to estimate the long-wave radiation emitted downward by the atmosphere and upward by the ground.
Some researchers have carried out statistical analyses on the long-term observation data of 82 meteorological stations in China and proposed Equation (4) for the sky temperature, which fits well with the measured data [14]:
T sky = [ 0.9 T g 4 ( 0.32 0.026 P q ) ( 0.30 + 0.70 S ) T a 4 ] 1 / 4
where T g is the surface temperature, K; T a is the air temperature at 1.5~2.0 m from the ground, K; P q is the partial pressure of the water vapor in the air near the ground, mbar; and S is the sunshine rate, that is, the ratio of actual sunshine hours to possible sunshine hours for the whole day.
Considering the diurnal difference caused by the influence of cloud cover on the insolation rate, the long-wave irradiance of the atmosphere on the horizontal plane during the day [18] is written as Equation (5):
T sky = σ T a 4 [ 0.904 ( 0.304 0.061 P w ) S h 0.005 P w ]
where T s k y is the daytime sky radiation on a horizontal surface, W/m2; σ is the Stefan–Boltzmann constant, 5.67 × 10 8 W/(m2K4); T a is the air temperature, K; P w is the water vapor partial pressure, hPa; and S h is the hourly sunshine rate. S h is replaced by 1 ( N h / 8 ) at night, where N h is the hourly cloud cover. Since the S h value is not stable during the first and last hours of the day, these two hours are considered nighttime periods.

2.4. Methods for Statistical Data

The   differences between the calculation results of the simplified T s o l model and the long-wave radiation T s o l model were   analyzed , and the Root Mean Square Error (RMSE) and Mean Bias Error (MBE) were   used to determine applicability. The level   of influence of the meteorological factors on the models ’ adaptability was analyzed using the correlation coefficient (Pearsons r). A linear regression equation for long-wave radiation estimation and the calculation results of the dynamic method were established based on the results of the dynamic calculation of the sky temperature. The coefficient of determination ( R 2 ) and the standard error of the regression line (STE) were used to evaluate the calibration equation.

2.5. Overview of the Study Area

In order to adapt the thermal engineering design of a building to the regional climate and to ensure the basic indoor thermal environment requirements, the building energy-saving design standards put forward the limit requirements of the prescribed indicators according to different climatic regions. The T sol model mainly involves the insulation characteristics of winter and summer. This paper selects 10 typical cities for the analyses, namely, Harbin, Xining, Turpan, Lhasa, Shanghai, Chongqing, Guangzhou, Yuanjiang, Guiyang, and Kunming, as shown in Table 2. The selection principle is mainly to select typical cities that cover all five thermal zones in China and that are located in different climate zones with different degrees of dryness and humidity.

3. Long-Wave Radiation Impact Analysis

3.1. Typical Daily T s o l Calculation

According to the simplified model (Equation (1)) and the long-wave radiation model (Equation (2)), the air temperature estimation method is used to determine the long-wave radiation in order to calculate the hourly values of T s o l in winter and summer, and the meteorological parameter data are typical meteorological year data in Chinese Standard Weather Date (CSWD) format [20]. The external surface heat transfer coefficient is 23 in winter and 19 in summer, and the solar radiation absorption coefficient of the external surface is 0.74.
The calculated value of the long-wave radiation T s o l model is set to the base value, and the Mean Bias Error (MBE) of the simplified T s o l model for each city is calculated, as shown in Table 3. A typical day in summer (July 22) and a typical day in winter (January 22) are selected to analyze differences hour by hour, as shown in Figure 2. The Root Mean Square Error (RMSE) is used as the evaluation index to analyze the applicability of the simplified T sol model under various climatic conditions, as shown in Figure 3.
It can be seen in Table 3 that, compared with the long-wave radiation T sol model, the simplified T sol model overestimates the winter T sol value overall (3.01 °C ≤ MBE ≤ 4.26 °C), because the formula does not take into account long-wave radiation cooling, so the calculated value is generally high. However, in summer, the simplified T sol model generally underestimates T sol (−1.58 °C ≤ MBE ≤ 2.73 °C) at all stations, except for the Xining, Turpan, and Lhasa stations. This is because the long-wave radiation T sol model uses the air temperature estimation method to calculate the long-wave radiation. Affected by the relative humidity, the long-wave radiation value of the building envelope is lower than the calculated value of the sky long-wave radiation. In summer, the long-wave radiation cooling effect of the building envelope coupled with the low relative humidity is significant.
In Figure 2a, it can be seen that the simplified T sol model overestimates all of the T sol values of a typical day in winter. The temperature difference ranges from 2.96 to 4.94 °C, the daily fluctuation of the temperature difference is small (0.14 to 1.12 °C), and the mean value of the difference at each station is relatively similar (3.02 to 4.34 °C). The temperature differences of the models in Lhasa, Xining, Turpan, and Harbin are relatively high from 12:00 to 22:00. Among them, the temperature difference of the model is the largest in Lhasa at night, while in the other areas, the temperature difference shows a decreasing trend from 12:00 to 22:00. This is because a clear sky at night is more similar to a theoretical black body than a cloudy sky, so a clear sky at night can act as a good absorber of radiant heat and cause the temperature of the structure to drop significantly. It is worth noting that the temperature difference of the model in Harbin is smaller than that in areas with high outdoor temperatures; that is, the temperature difference of the model is not consistent with the thermal zone in winter.
In Figure 2b, it can be seen that, in summer, the simplified T sol model underestimates the T sol (−3.86~3.96 °C) of all regions, except for Turpan and Lhasa. The daily fluctuation of the temperature difference is larger in summer than in winter (0.48~2.93 °C). The temperature differences of the models show the same overall downward trend from 8:00 to 24:00. Among them, the nocturnal temperature differences in Shanghai and Guangzhou are the highest at −3.86 and −3.18 °C, respectively, and the average temperature difference in the other regions is around −1 °C. The average temperature differences in Turpan and Lhasa are 3.02 and 0.83 °C, respectively, among which, the temperature difference in Turpan reaches a maximum of 3.96 °C at night; this means that long-wave radiation has an obvious cooling effect at night in Turpan.
Overall, in areas with high humidity, the T sol calculated via the simplified T sol model is overestimated in winter and underestimated in summer. In winter, the calculated value of the simplified T sol model is 3.01~4.26 °C higher. Therefore, when calculating the heating load of buildings in winter, there will be adverse effects, especially in Lhasa, Xining, and Turpan, and the nighttime T sol value will be quite different from the actual value. Considering the small differences in the MBE of the winter models at different stations, it is recommended to correct the T sol value by 3∓4 °C when using the simplified model to calculate the heating load. The MBE of the summer models at most stations is about −1 °C, which has little impact on the calculation of the cooling energy consumption of air conditioners. The calculated values of the simplified T sol model at night in Shanghai and Guangzhou are relatively low, but the thermal insulation design of the envelope is limited by the maximum outdoor temperature. Therefore, the temperature differences of the nighttime model have relatively little effect on thermal insulation design, and it is safe for thermal insulation design that the calculated value of the simplified T sol model in summer is relatively large. The MBE in the Turpan area reaches 2.73 °C in summer, so the influence of long-wave radiation should be fully considered in thermal insulation design.

3.2. Model Applicability Analysis in Winter and Summer

Taking the Root Mean Square Error (RMSE) as the evaluation index, the models of ten stations are selected for an applicability analysis, as shown in Figure 3. Based on a comparison of the total solar radiation, relative humidity, and the average outdoor dry-bulb temperature in winter and summer at each station, Pearson coefficients of climatic factors are shown in Table 4, Table 5 and Table 6, the key climatic factors affecting the model selection are analyzed, as shown in Figure 4.
The calculation results of the simplified T sol model and the long-wave radiation T sol model are relatively similar in summer. Among them, Turpan, Shanghai, and Guangzhou have higher RMSE values than the other regions (1.44 °C ≤ RMSE ≤ 2.91 °C). The simplified T sol model in Harbin, Shanghai, Chongqing, and Guiyang has a high calculation accuracy in winter, and the RMSE is less than 3.5 °C. Overall, all of the regions, except for Harbin and Guiyang, should use the long-wave radiation T sol model for calculation, especially Turpan and Lhasa, as the simplified T sol model will produce values that largely differ from those of the actual situation.
The reason for the difference in the RMSE values at each station is related to the calculation of long-wave radiation for each meteorological parameter. Figure 4 shows the outdoor dry-bulb temperature, total solar radiation, and outdoor relative humidity of the 10 typical cities in winter and summer. The correlation coefficients for the winter RMSE are 0.30, 0.80, and −0.82, respectively, and the correlation coefficients for the summer RMSE are 0.74, 0.07, and −0.52, respectively.
The F value of the dry-bulb temperature is much greater than 1 in winter, so the null hypothesis can be rejected. The ability of the independent variable to explain the dependent variable is 73.3%. The summer relative humidity values have the highest F values, so the null hypothesis can be rejected. Small residuals between the heat exchange, good fitting, and high interpretation can allow for the prediction of the building radiation heat exchange using statistical methods.

4. Analysis of the Influence of Climate on Building Heat Gain

The long-wave radiation calculation values discussed above are obtained using the air temperature estimation method, which does not consider the influence of the insolation rate or cloud cover on the heat gain of a building. In order to propose a simplified calculation method of long-wave radiation suitable for different climatic regions, the F test is carried out on the annual heat gain of the buildings, which is calculated using the estimation method at each station based on the results of the dynamic calculation. A regression equation for sky temperature estimation and the calculation results of the dynamic method are established, and the applicability of the simplified long-wave radiation calculation method is analyzed using STE as the judgment index.

4.1. Model Settings

The Apache Sim anisotropic solar radiation model is selected, and the temperature of the building is determined before the simulation. There are 10 preprocessing days, direct shading and internal sun tracking are calculated using Sun Cast, the number of time steps is set to 10 min, and the average output time interval of the results is 60 min [2]. A typical multi-story residential building is used as a geometric model to analyze and calculate the heat gain of the envelope. A standard floor plan of such a building is shown in Figure 5. The total height of the building is 18.8 m, and the shape coefficient is 0.28. The long-wave radiation rate of the roof is 0.9, and the solar radiation absorption coefficient is 0.7. The long-wave radiation rate of the exterior wall is 0.9, and the solar radiation absorption coefficient is 0.5. The geometric model used for calculation is shown in Figure 6.

4.2. Heat Gain Analysis of the Envelope

The heat gain of the envelope determined using the dynamic method is set to calculate the long-wave radiation as the base value, and the hour-by-hour difference is analyzed for the whole year for each station using the estimation method, as shown in Figure 7.
The F test is carried out on the heat gain determined using the estimation method and the dynamic calculation. The results show that the p values are less than 0.05, except for at the Harbin station; that is, there is a significant difference between the estimation method and the dynamic calculation. When the building heat gains need to be calculated accurately, the estimation method cannot be used to calculate the sky temperature. It can be seen in Figure 7 that the medians of the differences in Harbin, Xining, Lhasa, and Chongqing are relatively low, 0.23, 0.47, 0.30, and 0.31 kW/h, respectively. The medians of the differences in Turpan, Guangzhou, and Kunming are all higher than 1 KW/h. The annual difference in the Chongqing area has the smallest change (Q3 − Q1 = 0.36 kW/h), while the annual difference is relatively high in Kunming (Q3 − Q1 = 1.48 kW/h).

4.3. Correction Equation for Building Heat Gain Calculation

The calculation results of the long-wave radiation estimation method and the dynamic method at each station are fitted, as shown in Table 7. The fitted equations are evaluated using R2 and STE, as shown in Figure 8.
As shown in Figure 8, the R2 of the linear regression equations of the long-wave radiation estimation method and the dynamic method are above 0.99 for all areas, except for Guangzhou and Kunming. The statistical test results of the fitting equation are all p ≤ 0.01, and the fitting results are good (0.2 ≤ STE ≤ 0.88). The lower the total solar radiation, the smaller the STE. When the long-wave radiative heat gain needs to be accurately calculated in areas where outdoor meteorological parameter data are lacking, the appropriate equation of the air temperature estimation method should be used.

5. Conclusions

The long-wave (LW) radiation exchange with the sky and surrounding surfaces depends on the form factors, temperature, and emissivity of a building’s surfaces, as well as those of the surrounding buildings. Most of the current studies are based on thermal equilibrium studies with ideal conditions, ignoring the urban environment and building surfaces with the same temperatures. When considering the local microclimate and heat exchange flux between a building and its surrounding environment, it should be noted that the geometric features are more complex in urban environments. First, surrounding buildings create a shadow effect that reduces the surface temperature from solar radiation during the day. Second, the reduction in the sky perspective factor caused by building density limits the radiation cooling of the sky. There are also secondary reflections of long-wave radiation that interact with each other depending on the layout of adjacent buildings. Third, the local air flow transformation of buildings caused by urban forms has an impact on long-wave radiation heat transfer. That is, the building layout is one of the reasons for the gap between the building performance and the measured performance.
Aysan [13] considers the long-wave radiation heat transfer between the external surface of the envelope and the surrounding environment and the long-wave radiation heat transfer between the external surface of the envelope and the sky atmosphere separately, and they also consider the reflected radiation in the total solar radiation intensity. The calculation results show that the exchange process of urban surfaces under different microclimates is strongly influenced by the characteristics of the surface and the urban configuration. In this pattern, experimental and computational results are presented demonstrating the impact of the geometry of the inner courtyard and the covering materials on the thermal condition near the building envelopes and, consequently, on T sol .
Therefore, when simulating the heat and mass flow in and around a building caused by long-wave radiation in different climate regions, the influences of the local microclimate and adjacent buildings should be considered to accurately predict the energy performance of the building using the T sol model. The conclusions obtained should be limited to the author’s example and should not be extended.
Through a comparative analysis of the T sol models and their long-wave radiation calculation methods under different climatic conditions, the following conclusions are drawn:
  • In areas with high relative humidity, the simplified T sol model has higher T sol values in winter and lower T sol values in summer. When calculating the heating load in winter, each station should be corrected by 3~4 °C. In the summer in Turpan, the calculated value is conservative, and the long-wave radiation effect at night should be fully utilized for thermal insulation design.
  • Except for in Harbin and Guiyang, the long-wave radiation T sol model should be used in all areas, especially in Lhasa and Turpan; there will be large errors in calculating T sol when using the simplified T sol model.
  • The relative humidity is a key meteorological factor affecting the applicability of the T sol model in summer; in areas with higher relative humidity, the simplified T sol model is more suitable. When the air temperature is used to estimate the long-wave radiation in winter, the emissivity of the air near the ground ε air should be corrected for application in different climatic regions and seasons.
  • When calculating the long-wave radiation, the air temperature estimation method can be used in Harbin and Chongqing, but the dynamic calculation of the sky temperature should be adopted in the other areas. When the long-wave radiation heat transfer needs to be accurately calculated in areas where the outdoor meteorological parameter data are lacking, the appropriate equation of the air temperature estimation method should be used.
In this study, the long-wave radiation heat transfer is determined by using the calculation parameters and the sky temperature with higher precision to conduct calculations and simulations, and the correct equation for the air temperature estimation is obtained. A correction method and a theoretical basis are provided for the regions where the accuracy of the long-wave radiation heat transfer calculation needs to be high. In future research, we will focus on areas where the standard error of the regression equation is large, and we will determine appropriate correction coefficients using actual measurements combined with dynamic simulation methods so as to provide a simpler and more accurate calculation method for practical engineering applications.

Author Contributions

Conceptualization, J.C.; methodology, J.C.; software, J.C.; validation, J.C.; formal analysis, M.W.; investigation, J.C.; resources, J.C.; data curation, J.C.; writing—original draft preparation, J.C.; writing—review and editing, J.C. and Y.F.; visualization, J.C.; supervision, J.C.; project administration, J.C.; funding acquisition, Y.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number no. 52178020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available on request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. O’Callaghan, P.W.; Probert, S.D. Sol-air temperature. Appl. Energy 1977, 3, 307–311. [Google Scholar] [CrossRef]
  2. ASHRAE. ASHRAE Guide and Data Book: Fundamentals and Equipment; American Society of Heating and Air-Conditioning Engineers: Atlanta, GA, USA, 1961. [Google Scholar]
  3. Elangovan, R.; Kumar, A.; Alur, R. Thermal Performance of Building Envelops. Energy Build. 2017, 8, 169–188. [Google Scholar]
  4. CIBSE. Chartered Institution of Building Services. In Environmental Design: CIBSE Guide A; CIBSE: London, UK, 2006. [Google Scholar]
  5. Thomas, O.; Ohlsson, K.E.A.; Ronny, Ö. Measurement of the environmental temperature using the sol-air thermometer. Energy Procedia 2017, 132, 357–362. [Google Scholar]
  6. GB 50176-2016; Code for Thermal Design of Civil Building. Ministry of Housing and Urban-Rural Development of the People’s Republic of China: Beijing, China, 2016.
  7. Nicolas, L.; Auline, R.; Marjorie, M. How Building Energy Models Take the Local Climate into Account in an Urban Context—A Review; Elsevier: Amsterdam, The Netherlands, 2019; Volume 116. [Google Scholar]
  8. Zhang, K.; Chang, Q. Modeling of downward longwave radiation and radiative cooling potential in China. J. Renew. Sustain. Energy 2019, 11, 66501. [Google Scholar]
  9. Evangelisti, L.; Guattari, C.; Asdrubali, F. On the sky temperature models and their influence on buildings energy performance: A critical review. Energy Build. 2019, 183, 607–625. [Google Scholar] [CrossRef]
  10. Mehdi, Z.; Yogi, G.D.; Elias, S. A review of clear sky radiative cooling developments and applications in renewable power systems and passive building cooling. Sol. Energy Mater. Sol. Cells Int. J. Devoted Photovolt. Photothermal Photochem. Sol. Energy Convers. 2018, 178, 115–128. [Google Scholar]
  11. Lu, X.; Xu, P.; Wang, H.; Yang, T.; Hou, J. Cooling potential and applications prospects of passive radiative cooling in buildings: The current state-of-the-art. Renew. Sustain. Energy Rev. 2016, 65, 1079–1097. [Google Scholar] [CrossRef]
  12. Hang, J.; Wang, D.; Zeng, L.; Ren, L.; Shi, Y.; Zhang, X. Scaled outdoor experimental investigation of thermal environment and surface energy balance in deep and shallow street canyons under various sky conditions. Build. Environ. 2022, 225, 109618. [Google Scholar] [CrossRef]
  13. Aysan, F. Comparative analysis of sol-air temperature in typical open and semi-closed courtyard spaces. Build. Simul. 2022, 15, 957–973. [Google Scholar]
  14. Liu, S.Y.; Huang, Y.F. Discussion on Effective Sky Temperature. Acta Energ. Sol. Sin. 1983, 1, 6. [Google Scholar]
  15. Naveros, I.; Bacher, P.; Ruiz, D.P. Setting up and validating a complex model for a simple homogeneous wall. Energy Build. 2014, 70, 303–317. [Google Scholar] [CrossRef]
  16. Li, D.H.; Yang, L.; Lam, J.C. Impact of climate change on energy use in the built environment in different climate zones—A review. Energy 2012, 42, 103–112. [Google Scholar] [CrossRef]
  17. Liu, D.L.; Jia, X.W.; Yang, J.L.; Dong, G.M. Analysis of Urban Radiative Field Simulations. J. Tsinghua Univ. Sci. Technol. 2019, 3, 6. [Google Scholar]
  18. Al-janabi, A.; Kavgic, M.; Mohammadzadeh, A.; Azzouz, A. Comparison of EnergyPlus and IES to model a complex university building using three scenarios: Free-floating, ideal air load system, and detailed. J. Build. Eng. 2019, 22, 262–280. [Google Scholar] [CrossRef]
  19. GBT 17297-1998; Names and Codes for Climate Regionalization in China-Climate Zones and Climatic Regions. State Administration for Market Regulation: Beijing, China, 1998.
  20. EnergyPlus Weather Data. Available online: https://energyplus.net/weather (accessed on 15 March 2022).
Figure 1. Technical roadmap of this study.
Figure 1. Technical roadmap of this study.
Sustainability 15 00839 g001
Figure 2. (a) Temperature variation values of sol-air temperature model on a typical day in winter. (b) Temperature variation values of sol-air temperature model on a typical day in summer.
Figure 2. (a) Temperature variation values of sol-air temperature model on a typical day in winter. (b) Temperature variation values of sol-air temperature model on a typical day in summer.
Sustainability 15 00839 g002
Figure 3. RMSE of sol-air temperature model.
Figure 3. RMSE of sol-air temperature model.
Sustainability 15 00839 g003
Figure 4. The mean values of meteorological parameters in winter and summer.
Figure 4. The mean values of meteorological parameters in winter and summer.
Sustainability 15 00839 g004
Figure 5. Typical floor plan of simulated building.
Figure 5. Typical floor plan of simulated building.
Sustainability 15 00839 g005
Figure 6. Simplified building mode.
Figure 6. Simplified building mode.
Sustainability 15 00839 g006
Figure 7. The difference in the heat gain of the envelope.
Figure 7. The difference in the heat gain of the envelope.
Sustainability 15 00839 g007
Figure 8. Standard error of regression line.
Figure 8. Standard error of regression line.
Sustainability 15 00839 g008
Table 1. Direct models according to ISO 13790 [5].
Table 1. Direct models according to ISO 13790 [5].
CorrelationSite
Tsky = Tamb − 11Temperate areas
Tsky = Tamb − 9Sub-polar areas
Tsky = Tamb − 13Tropical areas
Table 2. Outdoor climatic parameters of each region [19].
Table 2. Outdoor climatic parameters of each region [19].
Dividing RegionClimate Sub-RegionCityAltitude/mClimate RegionCalculation PeriodAverage Relative Humidity
(%)
Average Outside Temperature
(°C)
Total Solar Radiation
(W/m2)
Severe cold zoneI(B)Harbin143Sub-humid area (B)summer (Jun.–Aug.)72.7821.29190.57
winter
(Dec.–Feb.)
74.90−16.0867.72
I(C)Xining2296Sub-arid area (C)summer (Jun.–Aug.)66.6516.44208.66
winter
(Dec.–Feb.)
46.25−6.05107.17
Cold zoneII(B)Turpan37Extremely arid area (E)summer (Jun.–Aug.)33.1130.97223.11
winter
(Dec.–Feb.)
51.45−3.5481.40
II(A)Lhasa3650Sub-arid area (C)summer (Jun.–Aug.)59.8215.76239.73
winter
(Dec.–Feb.)
29.13−0.12142.21
Hot summer and cold winter zoneIII(A)Shanghai3Humid area (A)summer (Jun.–Aug.)82.6026.28165.74
winter
(Dec.–Feb.)
71.676.0898.69
III(B)Chongqing259Humid area (A)summer (Jun.–Aug.)78.0126.96118.04
winter
(Dec.–Feb.)
84.329.1738.65
Hot summer and warm winter zoneIV(B)Guangzhou41Humid area (A)summer (Jun.–Aug.)84.2028.01148.52
Winter
(Dec.–Feb.)
70.4314.53101.42
IV(B)Yuanjiang401Sub-humid area (B)summer (Jun.–Aug.)74.4328.34188.62
winter
(Dec.–Feb.)
64.1518.03146.64
Warm zoneV(A)Guiyang1224Humid area (A)summer (Jun.–Aug.)76.4923.25161.99
winter
(Dec.–Feb.)
80.616.5458.27
V(A)Kunming1887Humid area (A)summer (Jun.–Aug.)79.0519.93168.53
winter
(Dec.–Feb.)
66.209.21139.23
Table 3. The MBE of T sol model in each station (°C).
Table 3. The MBE of T sol model in each station (°C).
TimeHarbinXiningTurpanLhasaShanghaiChongqing
Summer−0.100.612.731.02−1.30−1.09
Winter3.013.703.654.263.473.28
Table 4. Annual Pearson calculation.
Table 4. Annual Pearson calculation.
PearsonHarbinXiningTurpanLhasaShanghaiChongqingGuangzhouYuanjiangGuiyangKunming
Global radiation (W/m2)0.093 −0.046 −0.253 −0.053 −0.204 −0.247 −0.136 −0.326 −0.082 0.110
Dry-bulb temperature (°C) 0.645 0.141 0.083 −0.088 −0.094 −0.560 0.007 −0.148 −0.075 0.076
External relative humidity (%) −0.117 0.087 0.020 −0.060 0.117 0.228 0.046 0.208 0.057 −0.017
Table 5. Results of the Pearson calculation in summer.
Table 5. Results of the Pearson calculation in summer.
PearsonHarbinXiningTurpanLhasaShanghaiChongqingGuangzhouYuanjiangGuiyangKunming
Global radiation (W/m2)0.069 0.207 −0.741 0.395 −0.611 −0.325 −0.763 −0.755 −0.212 −0.004
Dry-bulb temperature (°C) 0.214 0.294 −0.683 0.578 −0.431 −0.262 −0.707 −0.734 −0.074 0.054
External relative humidity (%) −0.180 −0.274 0.436 −0.504 0.405 0.202 0.718 0.695 0.083 −0.037
Table 6. Results of the Pearson calculation in winter.
Table 6. Results of the Pearson calculation in winter.
PearsonHarbinXiningTurpanLhasaShanghaiChongqingGuangzhouYuanjiangGuiyangKunming
Global radiation (W/m2)−0.244 −0.504 −0.336 −0.667 −0.532 −0.557 0.270 0.225 −0.404 0.162
Dry-bulb temperature (°C) 0.625 0.145 0.242 −0.380 −0.203 −0.738 0.178 0.495 −0.208 −0.015
External relative humidity (%) −0.037 −0.106 0.016 0.250 0.511 0.580 −0.193 −0.366 0.360 0.016
Table 7. The regression equation of long-wave radiation heat transfer at each station.
Table 7. The regression equation of long-wave radiation heat transfer at each station.
StationRegression EquationR2STE/(kW/h)
Harbiny = 1.03x + 0.7387R2 = 0.99820.52
Xiningy = 1.0118x + 0.8398R2 = 0.99480.74
Turpany = 1.0103x + 1.2301R2 = 0.99640.72
Lhasay = 1.001x + 0.6759R2 = 0.99540.70
Shanghaiy = 0.9984x + 0.6881R2 = 0.99580.56
Chongqingy = 0.9835x + 0.3093R2 = 0.99880.20
Guangzhouy = 0.9852x + 1.0088R2 = 0.98920.70
Yuanjiangy = 0.9956x + 0.9566R2 = 0.99190.63
Guiyangy = 0.9887x + 0.8529R2 = 0.99080.66
Kunmingy = 0.9793x + 0.9423R2 = 0.98910.88
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, J.; Fan, Y.; Wang, M. Simplified Calculation of Tsol Based on Dynamic Numerical Simulation of Tsky in Diverse Climates in China. Sustainability 2023, 15, 839. https://doi.org/10.3390/su15010839

AMA Style

Chen J, Fan Y, Wang M. Simplified Calculation of Tsol Based on Dynamic Numerical Simulation of Tsky in Diverse Climates in China. Sustainability. 2023; 15(1):839. https://doi.org/10.3390/su15010839

Chicago/Turabian Style

Chen, Jie, Yue Fan, and Menghan Wang. 2023. "Simplified Calculation of Tsol Based on Dynamic Numerical Simulation of Tsky in Diverse Climates in China" Sustainability 15, no. 1: 839. https://doi.org/10.3390/su15010839

APA Style

Chen, J., Fan, Y., & Wang, M. (2023). Simplified Calculation of Tsol Based on Dynamic Numerical Simulation of Tsky in Diverse Climates in China. Sustainability, 15(1), 839. https://doi.org/10.3390/su15010839

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop