Next Article in Journal
Geothermal Resources Recognition and Characterization on the Basis of Well Logging and Petrophysical Laboratory Data, Polish Case Studies
Previous Article in Journal
An Interleaved Phase-Shift Full-Bridge Converter with Dynamic Dead Time Control for Server Power Applications
Previous Article in Special Issue
Energy Evaluation of a PV-Based Test Facility for Assessing Future Self-Sufficient Buildings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Influence of Solar Radiation on Heat Delivered to Buildings for Heating

1
Faculty of Environmental Engineering, Lublin University of Technology, Nadbystrzycka 40B, 20-618 Lublin, Poland
2
Faculty of Electrical Engineering and Computer Science, Lublin University of Technology, Nadbystrzycka 38A, 20-618 Lublin, Poland
3
Faculty of Mechanical Engineering, Lublin University of Technology, Nadbystrzycka 36, 20-618 Lublin, Poland
*
Author to whom correspondence should be addressed.
Energies 2021, 14(4), 851; https://doi.org/10.3390/en14040851
Submission received: 30 November 2020 / Revised: 29 January 2021 / Accepted: 2 February 2021 / Published: 6 February 2021
(This article belongs to the Special Issue Buildings Integration of Renewable and Smart Energy Systems)

Abstract

:
Nowadays, the attention of designers and service providers is especially focused on energy efficiency and integration of renewable energy sources (RES). However, the knowledge on smart devices and automated, easily applicable algorithms for optimizing heating consumption by effectively taking advantage of solar heat gains, while avoiding overheating, is limited. This paper presents a simple method for taking into account the influence of solar heat gains in the form of solar radiation for the purposes of forecasting or controlling thermal power for heating of buildings. On the basis of field research carried out for seven buildings (five residential buildings and two public buildings) during one heating season, it was noticed that it was justified to properly narrow down the input data range included in the building energy model calculations in order to obtain a higher accuracy of calculations. In order to minimize the impact of other external factors (in particular wind speed) affecting the heat consumption for heating purposes, it was recommended to consider the data range only at wind speeds below 3 m/s. On the other hand, in order to minimize the impact of internal factors (in particular the impact of users), it was suggested to further narrow down the scope of the input data to an hour (e.g., 10–14 in multi-family residential buildings). During these hours, the impact on users was minimized as most of them were outside the building.

1. Introduction

Regarding climate change and emission of greenhouse gases, the use of renewable energy sources (RES), such as solar radiation, in buildings is very common nowadays [1,2]. Solar radiation may be used in buildings in the form of active systems, such as flat solar collectors [3], but is recommended especially if used passively without special additional investment and exploitation cost.
This may be, for example, passive heating through glazing [4], which may reduce the heat used for heating in the winter season [5].
However, due to solar heat gains, the overheating may occur in buildings not only during heat waves [6] or normal exploitation during the summer period [7,8], but also during the winter season [9].
This may be due to inaccurate control methods of heat supply to a room or building, which does not fully or properly take into account the presence and influence of solar gains on energy consumption, which was shown in the example of co-heating tests [10].
In this light, Danov et al. [11] proposed the method for introducing the solar gain into the calculation of total heat loss coefficient. Knudsen and Petersen [12] underlined that it is not easy to simultaneously minimize overheating and energy costs; therefore, they proposed economic model predictive control for space heating. The use of model predictive control (MPC) may significantly reduce the energy consumption for heating [13], but there are still some barriers in wide application of such control methods [14], connected also with influence of external and internal factors on the building energy consumption [15]. That is why a lot of studies try to address these issues [16,17,18,19].
However, to the best of the authors’ knowledge, the amount of information about smart devices and automated, easily applicable algorithms to optimize heating consumption by effectively taking advantage of solar heat gains, while avoiding overheating, is limited.
This work contributes to this field and presents recommended parameters as well as a simple, widely applicable methodology for proper consideration of solar radiation in building energy model.
Section 2 presents the buildings and methodology of research. Section 3 describes the results and discussion regarding the methodological aspects of considering solar radiation in the creation of a building energy model.

2. Materials and Methods

The field tests were carried out in 7 buildings (Table 1) located approximately 500 m from each other in the Lubelskie Voivodeship in Poland. These included two public buildings (a court building and a store building) and five multi-family residential buildings. The selection of such a set of buildings for this research was aimed at showing the obtained dependencies of the solar radiation influence on heat consumption not only in residential buildings, but also public buildings.
The buildings had thermal insulation in the external walls and the roof. However, the methodology of determining the influence of solar radiation on the amount of thermal power supplied to a building for heating is a universal one, and that is why it may be successfully applied to the buildings without thermal insulation. Besides, because of the simplicity of this method, it is not required to know the emissivity of the building or its parts, solar absorptivity, or surface temperatures of parts of the building (windows, walls, and roof), which necessitates conducting additional measurements and very detailed calculations, and in this way may limit wide usage of this approach.
Each building was equipped with a central hydronic heating system with vertical risers connected to convection radiators. The analyzed buildings were supplied with heat from the district heating network through a thermal node located on the ground floor of each building, which participated in the heat exchange between the heating system and the heating and domestic hot water installation in analyzed buildings.
In each of the heating nodes, a calibrated heat meter was used to measure the heat power supplied to a given building for heating purposes (Qh). The research consisted of the continuous measurement of the heat power supplied to the analyzed buildings for heating (Qh) in the period from 1 October 2015 to 30 April 2016.
In addition to the parameters related to the heat supply to the analyzed buildings, the following weather parameters were also measured using a weather station located in the immediate vicinity of the analyzed buildings: outdoor temperature, wind speed (average and maximum), solar insolation (S_insol in J/cm2), cloudiness (Cloud_1 in %), and cloudiness (Cloud_2 in octants).
Due to the large amount of measurement data, this analysis used only the data related to 1 h.

3. Results

This section presents an analysis of the obtained research results, which was divided into two stages. Section 3.1 shows the analysis of the influence of solar radiation on the heat power supplied to the buildings for heating, in order to properly narrow down the scope of input data for further analysis. Then, Section 3.2 presents the results of the outdoor temperature correction with regard to the parameters characterizing solar radiation.

3.1. Solar Radiation and Heat Power Supplied for Heating

Solar radiation acting on the building envelope has an impact on the temporary demand for heat power during the entire heating season. In turn, the intensity of solar radiation may be influenced by clouding of the atmosphere caused by natural phenomena, such as cloudiness or dust pollution. The influence of this external factor on the amount of heat power supplied to buildings can be taken into account, involving solar insolation or cloudiness. These parameters are available in the databases of weather applications and services, both in the field of historical and forecast data.
However, it should be remembered that, in parallel with solar radiation, many other external factors (in particular, outdoor air temperature and wind speed) and internal factors (user behavior and preferences), which require an appropriate computational approach, also affect the building.
Therefore, in the first place, a detailed analysis related to the impact of external factors related to solar radiation (S_insol, Cloud_1, and Cloud_2) on the heat power supplied to the building for heating purposes was carried out at different hourly ranges of data included in the analysis and wind speed values (please see Table 2), for the purpose of appropriately narrowing down the data scope for further analysis in Section 3.2.
Therefore, in the case of solar insolation (S_insol), the dependencies of the heat power supplied to the analyzed buildings in relation to the outdoor temperature were determined for S_insol from the ranges: at wind speeds (average and maximum) from three different ranges (v < 3 m/s, 3 ≤ v <6 m/s, and for v > 6 m/s) and two hourly intervals, from 6 am to 6 pm and from 10 am to 2 pm.
In the case of cloudiness expressed in % (Cloud_1) and cloudiness expressed in octants (Cloud_2), the dependencies of heat power supplied to the analyzed buildings in relation to the outdoor temperature at wind speeds (average and maximum) from three different ranges (v < 3 m/s, 3 ≤ v < 6 m/s and for v > 6 m/s) and two hourly intervals from 6 am to 6 pm and from 10 am to 2 pm for Cloud_1 from the ranges 0–20%, 20–40%, 40–60%, 60–80%, and 80–100%; and for Cloud_2 in the ranges 0–2 octant (little cloudiness), 3–5 octants (moderate cloudiness), 6–7 (cloudy), and 8 (overcast), were determined.
The results of the analyses for individual cases were developed in the form of graphs and were used to determine the correlation equations together with the determination coefficients (R2), which were presented in the Figure 1, Figure 2 and Figure 3 for building B2 and B5 and in the form of an exemplary summary of the obtained determination coefficients for building B5 (Table 3).
On the basis of Figure 1 and Figure 2, it can be seen that the logical dependencies of the parameters related to the influence of solar radiation on the heat power supplied to buildings are not always obtained, because with increasing cloudiness, the demand for heat power should clearly increase. The reason for this may be the influence of other external factors or users on the heat power needed for heating, as well as the different reactions of individual buildings to a given external factor.
On the other hand, in the case of solar insolation (Figure 3), the logical dependencies of heat power (Qh) on S_insol were obtained for building B5, and the same data range (for the average wind speed below 3 m/s, for 10 am–2 pm) as in the case of Cloud_1. This is because the demand for heat power decreases with increasing solar insolation.
Therefore, apart from the logical analysis of the obtained equations, it is also justified to carry out a detailed analysis of the obtained values of the determination coefficients (example for building B5 in Table 3), which show how accurately the given equation fits the obtained results.
In order to decide which of the wind speed ranges (v < 3 m/s, 3 ≤ v < 6 m/s, v > 6 m/s) is the most reliable and shows the smallest disturbances (has the highest coefficient of determination), the arithmetic averages of coefficients of determination were calculated. The average values were calculated for each of the wind speed ranges, separately for each building. Then, the arithmetic mean of all individual mean values of the determination coefficients of the analyzed buildings were calculated, separately for residential and public buildings, as shown in Figure 4. This analysis does not take into account the division into average and maximum wind speed.
Figure 4 shows that the highest coefficients of determination occur for the wind speeds below 3 m/s. This is the case for cloudiness in octants, solar insolation in J/cm2, and cloudiness in %. The obtained results are logical, because at low wind speed, the influence of this factor (as disturbance) will be the lowest. R2 for the analyzed external factors, in the case of multi-family residential buildings, is 0.8074 (for Cloud_2), 0.9037 (for S_insol), and 0.8482 (for Cloud_1). In turn, for public buildings, the values of the determination coefficients for the wind speed range below 3 m/s were lower (than for residential buildings) and they were 0.8284 (for Cloud_2), 0.8251 (for S_insol), and 0.8003 (for Cloud_1).
On the basis of Figure 4, the mean values of the determination coefficients for cloudiness in octants and cloudiness in percent were also compared. Cloudiness in percent showed visibly higher R2 values for residential buildings, and this was the case for v ˂ 3 m/s, 3 ≤ v ˂ 6 m/s, and v ≥ 6 m/s. For example, for wind speeds below 3 m/s, R2 for cloudiness in percent is greater by 5.05% than R2 for cloudiness in octants.
Then, in order to decide which of the analyzed hourly ranges (6 am–6 pm or 10 am–2 pm) is the most reliable and shows the lowest disturbances (the highest coefficient of determination), the value of arithmetic means of the determination coefficients were also calculated. Figure 5 shows the average coefficients of determination for all the analyzed five multi-family residential buildings and two public buildings, for each of the considered hourly ranges, and for the wind speeds below 3 m/s, without division into their average and maximum values.
From Figure 5, it can be concluded that higher average coefficients of determination occur at 10 am–2 pm. This relationship, although more visible for Cloud_1, is also visible for S_insol for residential buildings. This is due to the minimization of the influence of system users on heat consumption between 10 am and 2 pm, because residents are not at home due to their work or study, and thus the indoor temperature settings are not changed by, for example, adjusting the position of the thermostatic head by thermostatic radiator valves.
This relationship looks different in the case of average determination coefficients for public buildings. In the case of the results for the cloudiness effect, a much higher R2 (by 6.2%) for the 10 am–2 pm interval is visible, but in the case of solar insolation, a higher (by 1.5%) average determination coefficient was obtained for the 6 am–6 pm. This may be caused by a much greater influence of the people staying in public buildings (compared to residential buildings) also at 10 am–2 pm, which may disrupt the operation of the heating system by generating internal heat gains and other behaviors.
Another aspect analyzed is whether the use of average or maximum wind speed in a given hour in the analysis is recommended and has an impact on the obtained determination coefficients of characteristic equations for the analyzed factors determining the influence of solar radiation on heat consumption for heating purposes. Figure 6 summarizes the average determination coefficients for cloudiness (%) and solar insolation (J/cm2), separately for average and maximum wind speeds for the wind speed variant below 3 m/s, and with division into multi-family buildings and public buildings.
On the basis of the results presented in Figure 6, it can be concluded that higher determination coefficients are obtained, taking into account the maximum wind speed values in the analyses, rather than the average values. This is the case for both solar insolation and cloudiness for multi-family residential buildings and public buildings. This is due to the fact that taking into account the maximum wind speed minimizes, to a greater extent, the influence of wind speed as an external factor on disturbances in the process of determining the influence of solar radiation on the heat power supplied to buildings for heating purposes. However, the problem in this case may be that only some weather services have data on the maximum wind speed.
Despite these observations, next stage (in Section 3.2) analyzed both the variants of the average and the maximum wind speed, as they contribute to the high R2 values, especially in the case of multi-family buildings.
On the basis of the above-mentioned detailed results and their analyses, it was found that it was justified to use, in further considerations, only the data from 10 am to 2 pm for the average and maximum wind speed below 3 m/s, which guarantees the best possible consideration of the influence of both solar insolation and cloudiness on the amount of heat power used for heating. This approach to the issue minimizes the influence of both wind speed and system users as disturbing factors, which usually occur in parallel with solar insolation/cloudiness and may result in incorrect consideration of the influence of the analyzed external factors on the demand and prediction of heat power for heating the buildings.
This will allow, using a universal calculation algorithm, development of appropriate equations correcting the outdoor air temperature, separately for each building, which will appropriately take into account the influence of solar insolation or cloudiness, which was presented in Section 3.2.

3.2. Correction of the Outdoor Temperature Due to Solar Radiation

On the basis of the narrowed range of input data (see Section 3.1), the correction of the outdoor temperature due to solar insolation in J/cm2 (S_insol) and cloudiness in % (Cloud_1) was determined separately for each building, because each building may react differently to specific weather conditions. An example of determining the outdoor temperature correction both in terms of wind speed and solar insolation has already been presented by the authors in the article [20].
Nevertheless, for the purposes of this article, an example of determining the outdoor temperature correction due to S_insol on the example of a multi-family building, B5, is presented below.
In this process, the correlation equations were used for four sun exposure ranges (0–35 J/cm2, 35–100 J/cm2, 100–200 J/cm2, and >200 J/cm2) and five cloud cover ranges (0–20%, 20–40%, 40–60%, 60–80%, and 80–100%), which were worked out as part of the calculations in Section 3.1. Using these correlation equations, the values of heat power delivered to buildings were determined for the outdoor temperature range of −20 °C–10 °C (see Table 4).
Then, taking into account the fact that the correlation equations were determined for specific data, the mean insolation was calculated for each of the ranges from which the correlation equations were determined. Lastly, the largest possible logical range was selected (i.e., when an increase in solar insolation leads to a decrease in heat power, and an increase in cloudiness leads to an increase in heat power).
The next step was to determine the increase in the outdoor temperature (N), which was a preliminary correction of the outdoor temperature in the correlation equation for lower solar insolation, enabling us to obtain an identical result of heat power for higher solar insolation. In the case of the example in Table 4, in order to determine the initial correction of the outdoor temperature between the solar insolation ranging from 0–35 J/cm2 to 35–100 J/cm2. This calculation was performed according to the following method (Equations (1) and (2)).
3.0363 · ( t e + N ) + 51.569 = 2.6536   · t e +   49.119
N = 2.6536 · t e + 49.119 51.569 3.0363 t e  
Using Equation (2), the values presented in column three of Table 5 were calculated. Analogous calculations were performed between successive solar insolation intervals and for all values of the outdoor temperature (from −20 °C to 6 °C).
In the next step, the sums of preliminary outdoor temperature corrections between successive variants of the increase in the average solar insolation were determined (see Table 6). The value for solar insolation of 0–35 J/cm2 was now used as a reference value. The process was that for the second average insolation value. The initial temperature corrections from the previous step were copied. For the last two mean solar insolation values, the preliminary temperature corrections presented in Table 5 were summed up. For example, for solar insolation >200 J/cm2 for −3 °C, it was 6.011675827 (see Table 6), because we added preliminary temperature corrections from the solar insolation ranges 35–100 J/cm2, 100–200 J/cm2, and >200 J/cm2, i.e., 1.18502783, 0.972565571, and 3.854082426, respectively, from Table 6. Then, the average values of the temperature correction for the successive variants of the average solar insolation value (in the range 35–100 J/cm2 equal to 62.68 J/cm2, in the range 100–200 J/cm2 equal to 148.31 J/cm2, and in the range >200 J/cm2 equal to 266.52 J/cm2) were calculated and amounted to 1.689194, 2.784762, and 7.836201.
On the basis of the mean values of the outdoor temperature correction determined in this way and the corresponding mean values of solar insolation, Figure 7 was developed, and a correlation equation was determined. This is the correction of the outdoor temperature due to solar insolation for the analyzed multi-family building B5.
Thus, the correction of the outdoor temperature (Teeq) due to S_insol for the example building, B5, can be calculated from Equation (3).
T e e q = t e + ( 0.0309 · n 0.8192 ) ,
where:
Teeq—outdoor temperature correction, taking into account solar insolation (in the case of cloudiness, the formula looks similar, there are only other numerical values), °C;
t e —outdoor air temperature, °C;
n—solar insolation (marked as z in the case of cloudiness), J/cm2 (in the case of cloudiness, %).
The same was done in the case of Cloud_1, and on this basis, corrections of the outdoor temperature due to S_insol and Cloud_1 for the analyzed buildings were determined and are presented in Table 7.
The relationship in the form of the Teeq allows determination of the influence of external factors, in particular, solar insolation/cloudiness on the heat power supplied to the heating system.
Taking into account the obtained results, it can be concluded that the approach that considers the correction equations developed on the basis of the data for solar insolation and average wind speed is characterized by the highest accuracy and logic of the obtained relationships (see Figure 8). This is due to the larger logical ranges of data observed at this stage for the analyzed variants.
Additionally, it should be noted that cloudiness was determined in octants or percentages, which are not precise parameters (compared to solar insolation), but rather subjective and estimated without the use of precise measuring instruments. These parameters do not take into account the position of the Sun in the sky, which affects the amount of solar radiation reaching and generating heat gains in the analyzed buildings. Therefore, as a more reliable parameter (the average value of the determination coefficient for the equations of the outdoor temperature correction due to solar radiation for the analyzed buildings at the level of 0.885), it is recommended to use solar insolation to take into account the influence of solar radiation on the heat power supplied to the building for heating.

4. Conclusions

This work proposed a universal calculation algorithm that will allow determination of the actual impact of solar radiation on the heat power supplied to a room or building for heating purposes, and thus it will be able to counteract the phenomenon of overheating rooms and reduce heat consumption.
Each analyzed building reacted differently to external factors related to solar radiation (in particular, cloudiness or solar insolation), which results from the individual structure of the building, its orientation, ventilation system, tightness, and behavior of occupants.
Therefore, in the first place, a detailed and multi-variant analysis of the impact of external factors related to solar radiation (S_insol, Cloud_1, and Cloud_2) on the heat power supplied to the building for heating purposes was carried out with different ranges of hourly data included in the analysis and wind speed values, to appropriately narrow down the range of input data included in the calculations.
On the basis of the conducted research, it was found that it was justified to narrow down the scope of data included in the calculations (narrowing down to 10 am–2 pm, at wind speeds below 3 m/s) in order to minimize the disturbances caused particularly by wind and building users. This was due to the fact that for such a narrowed range of data, the correlation equations for solar insolation or cloudiness from the outdoor temperature were obtained, which were characterized by high average values of determination coefficients at the level of 0.9043 (for solar insolation) and 0.8827 (for cloudiness in %) for the analyzed residential buildings.
Additionally, when determining the outdoor temperature correction due to the influence of solar radiation, it was found that the obtained results for solar insolation were more accurate (R2 = 0.8850) than for cloudiness (R2 = 0.8472). Therefore, it was recommended to use solar insolation as a parameter as it was more reliable than cloudiness in the calculations of a building energy model.
Besides, it was noticed that because of the simplicity of the method, it was not required to know the emissivity of the building or its parts, solar absorptivity, or surface temperatures of parts of the building (windows, walls, and roof), which necessitates conducting additional measurements and very detailed calculations.
This was connected with the fact that this simple method was based on real data, which were analyzed in a specific way (described above) in order to extract the influence of solar heat gain in the form of solar radiation on thermal power for heating of buildings.
Further research in this area may be focused on the development of methods of controlling heating systems, which, taking into account external factors (in particular, solar insolation) and the factors related to the behavior and preferences of users of heating systems, will allow forecast of the heat power demand and optimize its supply to achieve high energy savings while maintaining the thermal comfort of users.

Author Contributions

Study conception and design: T.C.; Acquisition, analysis, and interpretation of data: T.C., A.M., and A.S.-O.; Drafting of manuscript: T.C., A.M., A.S.-O., A.S., P.M., P.W., Ł.G., M.R.D., and K.Ł.; Critical revision: T.C., A.M., A.S.-O., A.S., P.M., P.W., Ł.G., M.R.D., and K.Ł. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by research project, financed from European Funds (under the Intelligent Development Programme) by the Polish National Centre for Research and Development. Project number: POIR.04.01.02-00-0012/18.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Authors would like to thank Rafał Michałowski (a student of the Lublin University of Technology) for participating in the research. Unfortunately, he passed away before we started to write this article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Frattolillo, A.; Canale, L.; Ficco, G.; Mastino, C.C.; Dell’Isola, M. Potential for Building Façade-Integrated Solar Thermal Collectors in a Highly Urbanized Context. Energies 2020, 13, 5801. [Google Scholar] [CrossRef]
  2. Bianco, L.; Komerska, A.; Cascone, Y.; Serra, V.; Zinzi, M.; Carnielo, E.; Ksionek, D. Thermal and optical characterisation of dynamic shading systems with PCMs through laboratory experimental measurements. Energy Build. 2018, 163, 92–110. [Google Scholar] [CrossRef]
  3. Baccoli, R.; Frattolillo, A.; Mastino, C.; Curreli, S.; Ghiani, E. A comprehensive optimization model for flat solar collector coupled with a flat booster bottom reflector based on an exact finite length simulation model. Energy Convers. Manag. 2018, 164, 482–507. [Google Scholar] [CrossRef]
  4. Ménard, R.; Souviron, J. Passive solar heating through glazing: The limits and potential for climate change mitigation in the European building stock. Energy Build. 2020, 228, 110400. [Google Scholar] [CrossRef]
  5. Kumar, K.; Saboor, S.; Kumar, V.; Kim, K.H.; Babu, T.P.A. Experimental and theoretical studies of various solar control window glasses for the reduction of cooling and heating loads in buildings across different climatic regions. Energy Build. 2018, 173, 326–336. [Google Scholar]
  6. Porritt, S.M.; Cropper, P.C.; Shao, L.; Goodier, C.I. Ranking of interventions to reduce dwelling overheating during heat waves. Energy Build. 2012, 55, 16–27. [Google Scholar] [CrossRef]
  7. Fletcher, M.J.; Johnston, D.K.; Glew, D.W.; Parker, J.M. An empirical evaluation of temporal overheating in an assisted living Passivhaus dwelling in the UK. Build. Environ. 2017, 121, 106–118. [Google Scholar] [CrossRef]
  8. Gamero-Salinas, J.C.; Monge-Barrio, A.; Saanchez-Ostiz, A. Overheating risk assessment of different dwellings during the hottest season of a warm tropical climate. Build. Environ. 2020, 171, 106664. [Google Scholar] [CrossRef]
  9. Brideau, S.; Beausoleil-Morrison, I.; Kummert, M. Collection and storage of solar gains incident on the floor in a house during the heating season. Energy Procedia 2015, 78, 2274–2279. [Google Scholar] [CrossRef] [Green Version]
  10. Stamp, S.; Altamirano-Medina, H.; Lowe, R. Measuring and accounting for solar gains in steady state whole building heat loss measurements. Energy Build. 2017, 153, 168–178. [Google Scholar] [CrossRef]
  11. Danov, S.; Carbonell, J.; Cipriano, J.; Marti-Herrero, J. Approaches to evaluate building energy performance from daily consumption data considering dynamic and solar gain effects. Energy Build. 2013, 57, 110–118. [Google Scholar] [CrossRef]
  12. Knudsen, M.D.; Petersen, S. Economic model predictive control of space heating and dynamic solar shading. Energy Build. 2020, 209, 109661. [Google Scholar] [CrossRef]
  13. Afram, A.; Janabi-Sharifi, F. Theory and applications of HVAC control systems—A review of model predictive control (MPC). Build. Environ. 2014, 72, 343–355. [Google Scholar] [CrossRef]
  14. Wang, Z.; Hong, T. Reinforcement learning for building controls: The opportunities and challenges. Appl. Energy 2020, 269, 115036. [Google Scholar] [CrossRef]
  15. Bilous, I.; Deshko, V.; Sukhodub, I. Parametric analysis of external and internal factors influence on building energy performance using non-linear multivariate regression models. J. Build. Eng. 2018, 20, 327–336. [Google Scholar] [CrossRef]
  16. Hong, T.; Wang, Z.; Luo, X.; Zhang, W. State-of-the-art on research and applications of machine learning in the building life cycle. Energy Build. 2020, 212, 109831. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, R.; Lu, S.; Feng, W. A novel improved model for building energy consumption prediction based on model integration. Appl. Energy 2020, 262, 114561. [Google Scholar] [CrossRef]
  18. D’Amico, A.; Ciulla, G.; Tupenaite, L.; Kaklauskas, A. Multiple criteria assessment of methods for forecasting building thermal energy demand. Energy Build. 2020, 224, 110220. [Google Scholar] [CrossRef]
  19. Turski, M.; Sekret, R. Buildings and a District Heating Network as Thermal Energy Storages in the District Heating System. Energy Build. 2018, 179, 49–56. [Google Scholar] [CrossRef]
  20. Cholewa, T.; Siuta-Olcha, A.; Smolarz, A.; Muryjas, P.; Wolszczak, P.; Anasiewicz, R.; Balaras, C.A. A simple building energy model in form of an equivalent outdoor temperature. Energy Build. 2021, in press. [Google Scholar] [CrossRef]
Figure 1. Dependence of heat power (kW) on the outdoor temperature (°C) and cloudiness (octants), for maximum wind speed below 3 m/s, 6 am–6 pm, and multi-family building B2.
Figure 1. Dependence of heat power (kW) on the outdoor temperature (°C) and cloudiness (octants), for maximum wind speed below 3 m/s, 6 am–6 pm, and multi-family building B2.
Energies 14 00851 g001
Figure 2. Dependence of heat power (kW) on the outdoor temperature (°C) and cloudiness (%), for average wind speed below 3 m/s, 10 am–2 pm, and multi-family building B5.
Figure 2. Dependence of heat power (kW) on the outdoor temperature (°C) and cloudiness (%), for average wind speed below 3 m/s, 10 am–2 pm, and multi-family building B5.
Energies 14 00851 g002
Figure 3. Dependence of heat power (kW) on the outdoor temperature (°C) and solar insolation (J/cm2), for average wind speed below 3 m/s, 10 am–2 pm, and multi-family building B5.
Figure 3. Dependence of heat power (kW) on the outdoor temperature (°C) and solar insolation (J/cm2), for average wind speed below 3 m/s, 10 am–2 pm, and multi-family building B5.
Energies 14 00851 g003
Figure 4. Average determination coefficients for the analyzed external factors related to solar radiation (S_insol, Cloud_1, and Cloud_2), divided into three wind speed ranges for the analyzed multi-family buildings (MB) and public buildings (PB).
Figure 4. Average determination coefficients for the analyzed external factors related to solar radiation (S_insol, Cloud_1, and Cloud_2), divided into three wind speed ranges for the analyzed multi-family buildings (MB) and public buildings (PB).
Energies 14 00851 g004
Figure 5. Average values of the determination coefficients for the analyzed external factors, divided into two hourly ranges (6 am–6 pm and 10 am–2 pm) of the input data range for multi-family buildings (MB) and public buildings (PB).
Figure 5. Average values of the determination coefficients for the analyzed external factors, divided into two hourly ranges (6 am–6 pm and 10 am–2 pm) of the input data range for multi-family buildings (MB) and public buildings (PB).
Energies 14 00851 g005
Figure 6. Average determination coefficients for the analyzed external factors, and for average and maximum wind speed below 3 m/s for multi-family buildings (MB) and public buildings (PB).
Figure 6. Average determination coefficients for the analyzed external factors, and for average and maximum wind speed below 3 m/s for multi-family buildings (MB) and public buildings (PB).
Energies 14 00851 g006
Figure 7. Outdoor temperature correction (°C) depending on the solar insolation (J/cm2).
Figure 7. Outdoor temperature correction (°C) depending on the solar insolation (J/cm2).
Energies 14 00851 g007
Figure 8. Average determination coefficients of the equations of the outdoor temperature correction for solar insolation/cloudiness by average wind speed.
Figure 8. Average determination coefficients of the equations of the outdoor temperature correction for solar insolation/cloudiness by average wind speed.
Energies 14 00851 g008
Table 1. Characteristic of the analyzed buildings.
Table 1. Characteristic of the analyzed buildings.
BuildingBuilding TypeHeated Floor Area (m2)
B1multifamily4195
B2multifamily3970
B3multifamily2910
B4multifamily1845
B5multifamily4264
B6public–courthouse9248
B7public–shop 280
Table 2. Analyzed conditions of external factors related to solar radiation.
Table 2. Analyzed conditions of external factors related to solar radiation.
External Factor Ranges of External Factor Ranges of Wind SpeedHourly Ranges
solar insolation (S_insol) (J/cm2)0–35;
35–100;
100–200;
> 200;
v < 3 m/s;
3 ≤ v < 6 m/s;
V > 6 m/s;
6 am–6 pm;
10 am–2 pm;
cloudiness (Cloud_1) (%) 0–20;
20–40;
40–60;
60–80;
80–100;
v < 3 m/s;
3 ≤ v < 6 m/s;
V > 6 m/s;
6 am–6 pm;
10 am–2 pm;
cloudiness (Cloud_2) (octants)0–2;
3–5;
6–7;
8;
v < 3 m/s;
3 ≤ v < 6 m/s;
V > 6 m/s;
6 am–6 pm;
10 am–2 pm;
Table 3. Summary of the determination coefficients (R2) of the equations developed for multi-family building B5 for 6 am–6 pm and for 10 am–2 pm in parentheses.
Table 3. Summary of the determination coefficients (R2) of the equations developed for multi-family building B5 for 6 am–6 pm and for 10 am–2 pm in parentheses.
Parameter Wind Speedv < 3 m/s3–6 m/sv ≥ 6 m/s
0–2 octantsaverage0.9013 (0.9083)0.8431 (0.7886)- (-)
maximum0.9548 (0.9741)0.9125 (0.9186)0.8593 (0.8782)
3–5 octantsaverage0.8212 (0.8466)0.7675 (0.8303)1 (1)
maximum0.8363 (0.8131)0.7960 (0.8082)0.8026 (0.8579)
6–7 octantsaverage0.8307 (0.8600)0.8011 (0.8201)0.9911 (1)
maximum0.8238 (0.8628)0.8329 (0.8538)0.8062 (0.8304)
8 octantsaverage0.8585 (1.000)0.8342 (-)0.9647 (-)
maximum0.8496 (-)0.855 (-)0.8866 (-)
0–20%average0.9094 (0.9083)0.7674 (0.7886)- (-)
maximum0.951 (0.9741)0.9192 (0.9186)0.8796 (0.9022)
20–40%average0.8055 (0.8526)0.8116 (0.8526)- (-)
maximum0.8612 (0.9118)0.7787 (0.8008)0.8039 (0.8648)
40–60%average0.8673 (0.7428)0.5918 (0.6515)1 (1)
maximum0.9016 (-)0.9306 (0.8934)0.5852 (0.8069)
60–80%average0.7985 (0.8323)0.7565 (0.7698)0.9911 (1)
maximum0.7337 (0.8737)0.8218 (0.8360)0.7438 (0.7687)
80–100%average0.8501 (0.8655)0.8124 (0.8661)0.964 (-)
maximum0.858 (0.8772)- (-)0.8529 (0.8821)
0–35 J/cm2average0.8450 (0.8458)0.8034 (0.8324)0.9385 (1.0000)
maximum0.8562 (0.8847)0.8331 (0.8254)0.8339 (0.8514)
35–100 J/cm2average0.9051 (0.9051)0.7859 (0.8355)0.9858 (1.0000)
maximum0.9329 (0.8662)0.8964 (0.8920)0.8522 (0.8673)
100–200 J/cm2average0.9085 (0.9162)0.7832 (0.7975)- (-)
maximum0.9750 (0.9718)0.8879 (0.9068)0.8831 (0.8806)
>200 J/cm2average0.8172 (0.7841)0.8086 (0.6992)- (-)
maximum- (-)0.9061 (0.8759)0.7424 (0.7472)
Table 4. The results of calculations of the heat power (kW) supplied to the multi-family building B5 for heating, for the data range from 10 am to 2 pm, by an average wind speed of v < 3 m/s.
Table 4. The results of calculations of the heat power (kW) supplied to the multi-family building B5 for heating, for the data range from 10 am to 2 pm, by an average wind speed of v < 3 m/s.
Solar Insolation (J/cm2)
0–3535–100100–200>200
Outdoor Temperature (°C)y = −3.0363x + 51.569y = −2.6536x + 49.119y = −2.572x + 46.783y = −1.8021x + 39.18
1021.20622.58321.06321.159
924.242325.236623.63522.9611
827.278627.890226.20724.7632
730.314930.543828.77926.5653
633.351233.197431.35128.3674
536.387535.85133.92330.1695
439.423838.504636.49531.9716
342.460141.158239.06733.7737
245.496443.811841.63935.5758
148.532746.465444.21137.3779
051.56949.11946.78339.18
−154.605351.772649.35540.9821
−257.641654.426251.92742.7842
−360.677957.079854.49944.5863
−463.714259.733457.07146.3884
−566.750562.38759.64348.1905
−669.786865.040662.21549.9926
−772.823167.694264.78751.7947
−875.859470.347867.35953.5968
−978.895773.001469.93155.3989
−1081.93275.65572.50357.201
−1184.968378.308675.07559.0031
−1288.004680.962277.64760.8052
−1391.040983.615880.21962.6073
−1494.077286.269482.79164.4094
−1597.113588.92385.36366.2115
−16100.149891.576687.93568.0136
−17103.186194.230290.50769.8157
−18106.222496.883893.07971.6178
−19109.258799.537495.65173.4199
−20112.295102.19198.22375.222
Table 5. Preliminary correction of the outdoor temperature, example for a multi-family building B5.
Table 5. Preliminary correction of the outdoor temperature, example for a multi-family building B5.
Solar Insolation (J/cm2)
0–3535–100100–200>200
Outdoor Temperature (°C)y = −3.0363x + 51.569y = −2.6536x + 49.119y = −2.572x + 46.783y = −1.8021x + 39.18
6 0.0506537560.6958094661.160031104
5 0.176695320.7265601451.45937014
4 0.3027368840.7573108231.758709176
3 0.4287784470.7880615012.058048212
2 0.5548200110.818812182.357387247
1 0.6808615750.8495628582.656726283
0 0.8069031390.8803135362.956065319
−1 0.9329447020.9110642153.255404355
−2 1.0589862660.9418148933.55474339
−3 1.185027830.9725655713.854082426
−4 1.3110693941.003316254.153421462
−5 1.4371109571.0340669284.452760498
−6 1.5631525211.0648176064.752099533
−7 1.6891940851.0955682855.051438569
−8 1.8152356491.1263189635.350777605
−9 1.9412772121.1570696415.650116641
−10 2.0673187761.187820325.949455677
−11 2.193360341.2185709986.248794712
−12 2.3194019041.2493216766.548133748
−13 2.4454434671.2800723556.847472784
−14 2.5714850311.3108230337.14681182
−15 2.6975265951.3415737117.446150855
−16 2.8235681591.372324397.745489891
−17 2.9496097221.4030750688.044828927
−18 3.0756512861.4338257468.344167963
−19 3.201692851.4645764248.643506998
−20 3.3277344141.4953271038.942846034
Table 6. Sum of the preliminary outdoor air temperature correction values, example for a multi-family building, B5.
Table 6. Sum of the preliminary outdoor air temperature correction values, example for a multi-family building, B5.
Mean Solar Insolation (J/cm2)
Outdoor Temperature (°C)16.2462.68148.31266.52
6 0.0506540.7464631.906494
5 0.1766950.9032552.362626
4 0.3027371.0600482.818757
3 0.4287781.216843.274888
2 0.5548201.3736323.731019
1 0.6808621.5304244.187151
0 0.8069031.6872174.643282
−1 0.9329451.8440095.099413
−2 1.0589862.0008015.555545
−3 1.1850282.1575936.011676
−4 1.3110692.3143866.467807
−5 1.4371112.4711786.923938
−6 1.5631532.627977.38007
−7 1.6891942.7847627.836201
−8 1.8152362.9415558.292332
−9 1.9412773.0983478.748463
−10 2.0673193.2551399.204595
−11 2.1933603.4119319.660726
−12 2.3194023.56872410.11686
−13 2.4454433.72551610.57299
−14 2.5714853.88230811.02912
−15 2.6975274.039111.48525
−16 2.8235684.19589311.94138
−17 2.9496104.35268512.39751
−18 3.0756514.50947712.85364
−19 3.2016934.66626913.30978
−20 3.3277344.82306213.76591
Average 1.6891942.7847627.836201
Table 7. Equations of the equivalent outdoor temperature corrected due to solar insolation or cloudiness.
Table 7. Equations of the equivalent outdoor temperature corrected due to solar insolation or cloudiness.
Building ParameterWind SpeedEquation
B1solar insolation
(J/cm2)
averagete + (0.0162n + 0.097)
maximumte + (0.0277n − 0.7383)
cloudiness (%)averagete + (−0.0364z + 4.1243)
maximumte + (−0.0533z + 6.1008)
B2solar insolation
(J/cm2)
averagete + (0.0311n − 1.411)
maximumte + (0.0133n − 0.1827)
cloudiness (%)averagete + (−0.0489z + 5.2973)
maximumte + (−0.0671z + 6.7361)
B3solar insolation
(J/cm2)
averagete + (0.0269n − 1.0419)
maximumte + (0.0102n − 0.1478)
cloudiness (%)averagete + (−0.0891z + 8.5884)
maximumte + (−0.0832z + 6.7708)
B4solar insolation
(J/cm2)
averagete + (0.0209n − 0.7183)
maximumte + (0.0377n − 2.2845)
cloudiness (%)averagete + (−0.0169z + 1.8574)
maximumte + (−0.0694z + 5.7688)
B5solar insolation
(J/cm2)
averagete + (0.0309n − 0.8192)
maximumte + (0.0205n − 0.1347)
cloudiness (%)averagete + (−0.0257z + 2.4261)
maximumte + (−0.096z + 7.024)
B6solar insolation
(J/cm2)
averagete + (0.0302n − 4.5632)
maximumte + (0.007n − 0.0983)
cloudiness (%)averagete + (−0.0055z + 0.4776)
maximumte + (−0.0258z + 1.9022)
B7solar insolation
(J/cm2)
averagete + (0.105n − 16.683)
maximumte + (0.0033n − 0.0499)
cloudiness (%)averagete + (−0.1532z + 13.66)
maximumte + (−0.0816z + 8.0567)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cholewa, T.; Malec, A.; Siuta-Olcha, A.; Smolarz, A.; Muryjas, P.; Wolszczak, P.; Guz, Ł.; Dudzińska, M.R.; Łygas, K. On the Influence of Solar Radiation on Heat Delivered to Buildings for Heating. Energies 2021, 14, 851. https://doi.org/10.3390/en14040851

AMA Style

Cholewa T, Malec A, Siuta-Olcha A, Smolarz A, Muryjas P, Wolszczak P, Guz Ł, Dudzińska MR, Łygas K. On the Influence of Solar Radiation on Heat Delivered to Buildings for Heating. Energies. 2021; 14(4):851. https://doi.org/10.3390/en14040851

Chicago/Turabian Style

Cholewa, Tomasz, Agnieszka Malec, Alicja Siuta-Olcha, Andrzej Smolarz, Piotr Muryjas, Piotr Wolszczak, Łukasz Guz, Marzenna R. Dudzińska, and Krystian Łygas. 2021. "On the Influence of Solar Radiation on Heat Delivered to Buildings for Heating" Energies 14, no. 4: 851. https://doi.org/10.3390/en14040851

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