Impact on Thermal Energy Needs Caused by the Use of Different Solar Irradiance Decomposition and Transposition Models: Application of EN ISO 52016-1 and EN ISO 52010-1 Standards for Five European Cities

: To solve the series of heat balances that EN ISO 52016-1 uses to simulate the dynamic hourly energy requirements of a building, detailed climatic data are required as input. Differently from air temperatures, relative humidity and wind speed, which are easily measurable and available in databases, the direct and diffuse solar irradiances incident on the different inclined and oriented surfaces, which are fundamental for the evaluation of solar gains, must be estimated using one of the many regression models available in the literature. Therefore, in this work, the energy needs of buildings were evaluated with the simpliﬁed hourly dynamic method of EN ISO 52016-1 by varying the solar irradiance sets on inclined and oriented surfaces obtained from EN ISO 52010-1 and three other pairs of solar irradiance separation and transposition models. Five European locations and two different window solar transmission coefﬁcients (g gl ) were analysed. The results showed that on average, for the heating period and for both g gl , the use of the different methods causes an average error on the calculation of the annual demand of less or slightly more than 5%; while for the cooling period, the average error on the calculation of the annual demand is 16.4% for the case study with g gl = 0.28 and 25.1% for the case study with g gl = 0.63. On the other hand, analysing the root-mean-square-error of the hourly data, using the model contained in TRNSYS as a benchmark, for most of the cases, when varying window orientations, cities and g gl , the model that diverges furthest from the others is that contained in EN ISO 52010-1.


Introduction
To date, buildings and their construction account for more than a third of the World's final energy consumption and almost 40% of total CO 2 emissions [1]. In fact, around 35% of buildings in the EU are over 50 years old, 75% of the building stock is energy inefficient and the annual renovation rate is just 1% [2]. To promote energy-conscious design, the European Committee for Standardisation (CEN) approved a new set of standards (called 'EN ISO 52000 family' [3]) that provided member states with adequate computational tools to be able to assess the real energy performance of buildings, both at the design and energy requalification stages.
These standards aim at international harmonisation of the methodology for assessing the energy performance of buildings by offering a systemic and holistic approach that considers the dynamic interaction between systems, users and climatic conditions. Of them all, EN ISO 52016-1:2017 [4] is the one that introduces the dynamic hourly method for calculating heating and cooling energy needs, indoor temperatures and sensible and latent heat loads, requiring a limited number of input data and ensuring reproducibility and transparency of calculation [5]. This standard, which is proposed as an alternative to the use of the main simulation software, has recently generated interest from the scientific community, which has begun to evaluate its accuracy. Ballarini et al. [6], comparing the standard's method with Energy Plus, highlighted a discrepancy in terms of thermal loads attributable to the use of different surface heat transfer coefficients and a different modelling of extra thermal radiation to the sky. Zakula et al. [7], through a comparative analysis with TRNSYS, showed a difference in results caused by the use of constant values of the thermal transmittance of glazed surfaces (F w ), the use of constant values of the solar transmission coefficient (g gl ) and to a lesser extent by the heat transfer pattern of opaque elements. With regard to the latter problem, Mazzarella et al. [8] analysed and compared with both TRNSYS and the analytical solution, the model implemented in the Italian National Annex that proposes an alternative method for defining the RC network of opaque elements based on the thermo-physical characteristics of each individual layer that makes up the building element. The use of this model reduced the error on the amplitude of the internal flux by up to 67% and generated an average error on the external flux of only 3% compared to the heat transfer model of opaque surfaces proposed in standard in EN ISO 52016-1. As for the heat transfer model of opaque elements, the problem of using constant values of the solar transmission coefficient (g gl ) was also improved.
In the study by Summa et al. [9], a model (also implemented in the Italian National Annex of EN ISO 52016-1) was tested to calculate an hourly solar transmission coefficient as a function of the angle of incidence of solar radiation, but also as a function of the number of glass panels of the window and the type of glass coating (float glass, low emissivity, solar control, etc.). The results obtained, using TRNSYS as a baseline, showed a 47% improvement in the root-mean-square deviation (RMSE) of summer solar loads (thermal loads caused by solar radiation entering from glazed surfaces) compared to the method proposed by the standard EN ISO 52016-1 with constant g gl . These results were also confirmed in study [10], where it is affirmed that the use of hourly g gl correction factors resulted to be the most sensitive modelling option, among those analysed. As can be observed from the literature studies, EN ISO 52016 [4] shows some discrepancies in terms of energy requirements compared to the main dynamic simulation software (i.e., Trnsys, Energy Plus), but researchers are analysing and proposing improvements to the algorithm, limiting the demand for new input data from designers, which could introduce subjective uncertainty and inaccuracy into the calculation [11].
However, the results of dynamic energy simulations are influenced not only by the calculation models, but also by the accuracy of the hourly variables required as input to EN ISO 52016 [4] and from additional calculation algorithms, such as climatic data. The meteorological parameters, such as dry bulb temperature, relative humidity, solar irradiance, wind speed and direction, are in fact those that have the greatest influence on the energy balance of a building [12,13], and among these, solar irradiance is the most impactful phenomenon [14]. Differently from the other meteorological parameters that are directly measured and present in the Typical Meteorological Years (TMY/TMY2), the measured solar radiation is only the global horizontal radiation (for TMY) and in some cases also the direct normal irradiation (for TMY2). In order to reconstruct diffuse, direct and global irradiance data for all surfaces with arbitrary tilt and orientation, auxiliary calculation methods are needed. To overcome this problem, CEN decided to introduce a single consistent source for all energy performance of buildings standards for the conversion of climate data for energy calculations [15], i.e., EN ISO 52010-1:2017 [16]. This standard allows to calculate the hourly solar irradiance on a surface with arbitrary orientation and inclination from measurements of global horizontal irradiance, temperature, humidity, solar constant, albedo and geographical coordinates of the site. Few studies exist in the literature that attempt to validate EN ISO 52010-1 [16]. In [17], a comparative analysis of four different models for the calculation of direct and diffuse irradiance on arbitrarily inclined and oriented surfaces was performed for five European cities. The results showed that, in comparison to the models implemented in TRNSYS, the model implemented in EN ISO 52010 [16] is the one that mostly overestimates the diffuse solar irradiance on the horizontal plane and underestimates the annual direct irradiance on south, west, east and north facing surfaces by 24.2%, 55.9%, −0.9% and 69.4% respectively. In [18], instead, through the monthly quasi-steady-state method of EN ISO 13790 [19], the energy performance of two typical residential buildings in five Polish cities was calculated, using as input the solar irradiance on inclined and oriented surfaces obtained from 15 different calculation models, including that of EN ISO 52010 [16]. Taking as baseline the TMYs containing irradiances on sloping surfaces with five slope angles and oriented in eight directions, it was seen that the use of EN ISO 52010 [16] introduces a difference on the energy performance indicator for heating, domestic hot water and ventilation (EPH+W) of 1.6 and 1.2 kWh/m 2 for the singlefamily building, respectively, for the use of a g gl = 0.75 and g gl = 0.50, and a difference of 5.1 and 3.9 kWh/m 2 for the multi-family building, varying its energy classification. Therefore, in order to extend the knowledge of these issues, in this research a comparative analysis of the solar loads and heating and cooling energy needs calculated with the simplified hourly dynamic method of EN ISO 52016-1 [4] will be carried out for the first time, varying the solar irradiance sets on inclined and oriented surfaces obtained from EN ISO 52010 [16] and three other pairs of solar irradiance separation and transposition models. The simulations will be conducted for five European locations, on a cubic building, changing each time the orientation of the glazed surface and the solar transmission coefficient (g gl ).
After this Introduction, this paper will follow with Section 2 describing the Methodology, Section 3 illustrating the Case Study, Section 4 discussing the Results and finally Section 5 with the main conclusions and future developments.

Methods
Energy requirements obtained from hourly dynamic simulations performed with EN ISO 52016-1 [4] were compared by varying the input solar irradiances on inclined and oriented surfaces simulated with four different algorithms including the one contained in EN ISO 52010-1:2017 [16].
Both of the above-mentioned standards are part of the new series of standards called EPB 'Energy Performance of Buildings', which aim at international harmonisation of the methodology for assessing the energy performance of buildings, from the identification of input data to the calculation of primary energy (see Figure 1).
EN ISO 52016-1 [4] potentially allows hourly dynamic calculations of heating and cooling energy needs, indoor air temperatures, mean radiant temperatures of the surfaces enclosing the heated thermal zone and operative temperatures of individual thermal zones in the building, in a simplified way, but with accuracy similar to the main simulation software (TRNSYS 17, Energy+ 9.6.0, etc.).
EN ISO 52010-1 [16], starting from the geographical coordinates and the hourly horizontal global irradiance of the site, allows to obtain the direct and diffuse irradiances on any surface with arbitrary orientation and inclination. These irradiances are those required by EN ISO 52016-1 [4] for building energy calculations but can also be used in standards for the design/verification of technical building system components, such as solar thermal collectors and photovoltaic panels. In addition to EN ISO 52010-1 [16], three other algorithms for splitting direct and diffuse radiation, previously compared in [17], were used, i.e., the one implemented: (i) in the work of Summa et al. [17], that makes modifications to the UNI 10349-1:2016 standard [20], (ii) in Meteonorm software [21] and (iii) that implemented in Type 16c of the TRNSYS software [22]. The models used by each algorithm for the split between direct and diffuse radiation are shown in Table 1.  is the clearness index of the atmosphere (relating to diffuse radiation G sol;g is the global horizontal solar irradiance (climatic data set) (W/m 2 ). I ext,h = I ext · sin α sol is the global extra-atmospheric solar irradiance on a horiz plane (W/m 2 ).
I ext is extra-terrestrial radiation (W/m 2 ). α sol is solar altitude angle (°). The Erbs et al. model [25] is used to calculate the diffuse solar irradiance on hor tal plane G sol;d (W/m 2 ) as follows: If k T ≤ 0.22 → G sol;d = 1.0 -0.09·k T ·G sol;g ,   [22] R-P Reindl et al.1990 [28] Perez et al. 1990 [26] The Boland et al. model [23] is used to calculate the diffuse solar irradiance on horizontal plane G sol;d (W/m 2 ) as follows: where k T,d = G sol;g I ext,h is the clearness index of the atmosphere (relating to diffuse radiation). G sol;g is the global horizontal solar irradiance (climatic data set) (W/m 2 ). I ext,h = I ext ·sin α sol is the global extra-atmospheric solar irradiance on a horizontal plane (W/m 2 ).
I ext is extra-terrestrial radiation (W/m 2 ). α sol is solar altitude angle ( • ). The Erbs et al. model [25] is used to calculate the diffuse solar irradiance on horizontal plane G sol;d (W/m 2 ) as follows: where k T = G sol;g I ext is the clearness index of the atmosphere (relating to beam radiation). The Perez et al. 1991 model [27] is used to calculate the normal solar irradiation G sol;n (W/m 2 ) as follows: where I DISC is direct normal irradiation estimate by the DISC model (itself a function of global irradiance and solar zenith angle) [29].
X k t , Z, W, ·k t are obtained from 4-dimensional look-up table based on the bins in Table 1 of paper [30].
where if k t < 0.2 → k b = 0 . otherwise k b = a k t , Z, W, ∆k t ·k t +b k t , Z, W, ∆k t . a k t , Z, W, ∆k t and b k t , Z, W, ∆k t are obtained from 4-dimensional look-up table based on the bins in Table 1 of paper [30].
The Reindl et al. 1990 model [28] is used to calculate diffuse solar irradiance on horizontal plane G sol;d (W/m 2 ) as follows: The Liu et al. 1960 model [24] is used to calculate direct irradiance on tilted/oriented surface I dir (W/m 2 ) as follows: where G sol;b is the direct solar irradiance on horizontal plane (W/m 2 ). θ sol;ic is the solar angle of incidence on the inclined surface. G sol;b sin α sol is calculated using the Formula (18) proposed in [17].
The  and the * model [26] are used to calculate diffuse irradiance on tilted/oriented surface I dif (W/m 2 ) as follows: where a, b are dimensionless parameters. F 1 , F 2 is the circumsolar and horizontal brightness coefficient, function of ε. ε is the dimensionless clearness parameter, at anisotropic sky conditions.
but in EN ISO 52010 [16] the angle θ z was replaced by the angle α sol and the coefficient 1.041 was replaced by 1.014

Case Study
Design choices for the case study are illustrated and justified below, with the aim of showing how the use of different hourly calculation algorithms for the prediction of direct and diffuse solar irradiances on inclined and/or oriented surfaces can vary the outputs of EN ISO 52016-1 [4].

Locations
The five European locations considered in the case study are: De Bilt (Netherlands), Malaga (Spain), Paris (France), Rome (Italy) and Stockholm (Sweden). As already specified in [17], the choice of cities characterised by very different climates (see Table 2) made it possible to test the performance and limits of the four different calculation algorithms as a function of geographical area. For each city, climate data (TMY) were obtained from weather stations officially recognised by the World Meteorological Organization (WMO).

Geometry and Opaque Surfaces
The case study's simplified cube-shaped building has internal dimensions 3.00 m × 3.00 m × 3.00 m. There are no internal partitions and therefore only one airconditioned thermal zone is considered. Table 3 shows the geometric characteristics of the building, identifying the net internal area (A floor_int ), the external area of the walls for each orientation (A wall ), the area of the flat roof (A roof_hor ) and the internal volume (V). The thermal transmittance of opaque elements (external walls and flat roof) was determined by considering the thermal transmittances of the individual locations studied, then choosing the average value of approximately 0.50 W/m 2 K.

Transparent Surfaces
Two window types were simulated, one with standard glass and one with lowemissivity glass, in order to assess the impact of this parameter on the final calculation. Both windows have dimensions of 1.20 m × 2.10 m, a thermal transmittance of 1.10 W/m 2 K and zero frame area. The window with standard double glazing and argon gas cavity has a solar transmission factor (g gl ) of 0.63, while the window with low-emissivity double glazing and argon gas cavity has a g gl of 0.28 (see Table 4). For each window orientation (north, south, east, west, horizontal) a different simulation was performed in order to separate the effects of the direct and diffuse irradiances simulated with the four algorithms on the calculation of the energy demand obtained with EN ISO 52016-1 [4].

Other Assumptions
Thermal bridges, shading caused by internal or external elements and internal gains were not considered. Ground exchange was considered adiabatic. The ventilation rate was set at 0.50 1/h. The heating and cooling system is assumed on 24 h a day with set-point temperatures of 20 • C and 26 • C, respectively.
To summarise, in this work, one building, two types of glazed surfaces, four different calculation procedures for the prediction of direct and diffuse solar radiation on inclined/oriented surfaces, five locations and five window orientations were considered, for a total of 200 simulations (see Figure 2).  The thermal transmittance of opaque elements (external walls and flat roof) was determined by considering the thermal transmittances of the individual locations studied, then choosing the average value of approximately 0.50 W/m 2 K.

Transparent Surfaces
Two window types were simulated, one with standard glass and one with low-emissivity glass, in order to assess the impact of this parameter on the final calculation. Both windows have dimensions of 1.20 m × 2.10 m, a thermal transmittance of 1.10 W/m 2 K and zero frame area. The window with standard double glazing and argon gas cavity has a solar transmission factor (ggl) of 0.63, while the window with low-emissivity double glazing and argon gas cavity has a ggl of 0.28 (see Table 4). For each window orientation (north, south, east, west, horizontal) a different simulation was performed in order to separate the effects of the direct and diffuse irradiances simulated with the four algorithms on the calculation of the energy demand obtained with EN ISO 52016-1 [4].

Other Assumptions
Thermal bridges, shading caused by internal or external elements and internal gains were not considered. Ground exchange was considered adiabatic. The ventilation rate was set at 0.50 1/h. The heating and cooling system is assumed on 24 h a day with set-point temperatures of 20 °C and 26 °C, respectively.
To summarise, in this work, one building, two types of glazed surfaces, four different calculation procedures for the prediction of direct and diffuse solar radiation on inclined/oriented surfaces, five locations and five window orientations were considered, for a total of 200 simulations (see Figure 2).

Statistical Descriptors
The simulated energy needs for the different cities and with the different solar irradiance data sets were compared using the average, the standard deviation, the mean error and the RMSE (root-Mean-Square Error).
The average of the energy demand value (heating or cooling) for each city is: where x i : i-th annual energy demand value (heating or cooling) determined using as input the set of solar irradiances calculated with the models described in Section 2.
N: number of models used (N = 4). The standard deviation is used to express the dispersion of the data around the position index, which in our case is the arithmetic mean of the annual energy needs for heating and cooling for each city where x i : i-th annual energy demand value (heating or cooling) determined using as input the set of solar irradiances calculated with the models described in Section 2.
x: average value of the annual energy demand (heating or cooling) determined using as input the set of solar irradiances calculated with the models described in Section 2.
N: number of models used (N = 4). The mean error or uncertainty for each city is: The RMSE, on the other hand, being proportional to the size of the squared error, is also sensitive to outliers [33,34] and in this case indicates the distance in terms of hourly energy needs of the B-L, E-P* and P-P models compared to the R-P model of TRNSYS.
where P i : i-th energy need value determined using as input the set of solar irradiances calculated with TRNSYS type 16 (R-P).
O i : i-th energy need value determined using as input the solar irradiance sets calculated with the B-L, E-P* and P-P models. n: number of hours when there is an energy need for heating/cooling. The solar irradiances calculated with TRNSYS Type 16c are taken as benchmark, because this type contains the two split models considered most valid in the literature, namely Reindl [28] and Perez [26]. In addition to Perez's well-known anisotropic sky model, Reindl's model is the only one of those analysed that includes multiple predictor variables in the calculation of direct solar irradiance in the horizontal plane, i.e., the clarity index, solar altitude angle, air temperature and relative humidity. Furthermore, Reindl states that although the clarity index is the most important variable for low or medium clearness index (k T ) values (cloudy sky), in the case of high k T values (clear sky) the solar altitude angle becomes predominant, improving the irradiance prediction compared to Erbs' model by 26% [28] (model used by EN ISO 52010-1 [16]).

Impact of Solar Decomposition and Transposition Models on Calculated Energy Demand
The annual energy needs for heating and cooling were compared with the respective solar loads by varying the orientation of the glass surface for the five European locations studied. The solar loads are calculated from the solar irradiances (input data) simulated with the four different calculation models described in Section 2. Each point in Figures 3 and 4  Looking at the case of the g gl = 0.28, it is evident from Figure 3 that the winter solar loads are less dispersed than the summer ones, despite varying location, window orientation and calculation method.

Impact of Solar Decomposition and Transposition Models on Calculated Energy Demand
The annual energy needs for heating and cooling were compared with the respective solar loads by varying the orientation of the glass surface for the five European locations studied. The solar loads are calculated from the solar irradiances (input data) simulated with the four different calculation models described in Section 2. Each point in Looking at the case of the ggl = 0.28, it is evident from Figure 3 that the winter solar loads are less dispersed than the summer ones, despite varying location, window orientation and calculation method. Focusing on the winter period (see Table A1), it can be observed that on buildings with a rooflight (horizontal window placed on the ceiling), the four irradiance calculation models lead to similar standard deviations (henceforth dev.s) for all five locations, both in terms of incoming solar loads (1.2 kWh/year) and in terms of energy needs (16.0 kWh/year). In fact, taking as baseline the hourly heating needs obtained using the solar irradiances calculated with the TRNSYS R-P split method, the RMSEs of the other three models are very low, with values of 6.1 Wh, 6.3 Wh and 9.5 Wh for the B-L, P-P and E-P* methods, respectively (see Table 5). Focusing on the winter period (see Table A1), it can be observed that on buildings with a rooflight (horizontal window placed on the ceiling), the four irradiance calculation models lead to similar standard deviations (henceforth dev.s) for all five locations, both in terms of incoming solar loads (1.2 kWh/year) and in terms of energy needs (16.0 kWh/year). In fact, taking as baseline the hourly heating needs obtained using the solar irradiances calculated with the TRNSYS R-P split method, the RMSEs of the other three models are very low, with values of 6.1 Wh, 6.3 Wh and 9.5 Wh for the B-L, P-P and E-P* methods, respectively (see Table 5). As in the case of ggl = 0.28, it is evident from Table 8 that the south and west orientations show the highest average dev.s with errors of 5.5% (78.9 kWh/year) and 5.2% (80.2 kWh/year) on the calculation of the heating energy needs, respectively.
Focusing on the cities analysed, comparing the dev.s with the average annual winter demand estimated with EN ISO 52016-1 [4] when varying the solar radiation split methods, the average error among all orientations on the winter demand estimate is less than 5% for all locations except for Malaga, where the average error is 13.9% and the maximum error is 19.6% (south). Different from the case with ggl = 0.28, where all cities showed an error on all orientations of less than 5%, in the case with ggl = 0.63, De Bilt and Paris show an error on the annual demand estimate of 6.1% and 6.4% in the west, while Rome shows an error of 7.5% in the south.
Analysing, instead, the hourly heating energy needs, it can be seen that for all five orientations, the highest RMSE values are those obtained using as input the irradiances calculated with the E-P* method of EN ISO 52010-1 [16] with values of 11.6 Wh, 76.1 Wh, 41.0 Wh, 76.0 Wh and 37.3 Wh, respectively, for the horizontal, south, east, west and north window orientation (see Table 9). The methods that show better RMSE, on the other hand, are the B-L for the horizontal, east, west and north orientation with values lower than the maximum ones of the E-P* method of 24.0%, 34.0%, 35.7% and 64.2%, and the P-P method for the south orientation with lower RMSE than the maximum of 51.3%.
In the summer season (see Figure 4 and Table A2), the dev.s of the total solar loads for buildings with south, east, west and north windows is 15.6 kWh/year, 73.4 kWh/year, 142.9 kWh/year and 26.3 kWh/year, respectively, while the average dev.s of the annual cooling energy demand for buildings with south, east, west and north windows (see Table  9) results to be 19.7 kWh/year, 75.1 kWh/year, 117.1 kWh/year and 26.9 kWh/year with an average increase over the case study of ggl = 0.28 of 434.9%, 283.6%, 289.0% and 281.1%.  Table 5. Case study with g gl = 0.28. Average RMSE of heating and cooling demand when varying the solar irradiance transposition method. The blue boxes indicate the worst results for heating and the red boxes the worst results for cooling. The asterisk * used beside "E-P*" indicates the change in EN ISO 52010 [16] from the classical method of  on the calculation of the "dimensionless clarity parameter" (see Formula (11)). Remaining still on the winter period of the graph in Figure 3, considering the buildings with south, east, west and north glazed surfaces, the dev.s. of the total solar loads between the different models is similar and equal, respectively, to 29.4 kWh/year, 29.0 kWh/year, 26.6 kWh/year and 26.3 kWh/year, while the dev.s. of the annual heating energy demand is more variable (see Table 6), showing higher values for the south and west orientations with a corresponding average dev.s. among all the cities of the case study of 54.4 kWh/year and 56.8 kWh/year, equal in both cases to about 3.2% on the calculation of the heating energy demand. Focusing on the cities analysed, comparing the dev.s with the average annual winter energy needs estimated with the EN ISO 52016-1 [4] standard when varying the radiation splitting methods, an average error on the energy needs estimate of less than 5% is obtained for all locations, except for Malaga, where the average error is 6.7% and the maximum error 10.6% (south).

B-L E-P*
Analysing the hourly heating energy needs, it is noted that for all four orientations, the highest RMSE values are those obtained using as input the irradiances calculated with the E-P* method of EN ISO 52010 [16] with values of 41.0 Wh, 23.3 Wh, 45.4 Wh and 19.3 Wh, respectively, for the south, east, west and north orientation (see Table 5). On the other hand, the RMSE values of the B-L and P-P methods were similar and lower than the E-P* method of 47.9%, 16.7%, 33.6% and 54.8% for the south, east, west and north orientation, respectively.
Analysing the summer period (see Figure 3 and Table A1), unlike the winter period, the points appear to be more scattered. The dev.s of the total solar loads for the buildings with south, east, west and north glazed surfaces is 5.6 kWh/year, 18.6 kWh/year, 33.6 kWh/year and 8.7 kWh/year, respectively, while the average dev.s of the annual cooling energy demand (see Table 6) shows higher values mainly for the west orientation with average dev.s values of 30.1 kWh/year, which is 16.4% error in the calculation of the cooling energy demand.
Focusing on the analysed cities, comparing the dev.s with the average annual cooling energy requirements estimated with EN ISO 52016-1 [4], varying the radiation split methods produces an error on the energy requirement estimate significantly higher than 5% for all locations in the case study. In particular, the worst estimate occurs for the De Bilt location where the average error is 23.8% and the maximum error is 47.8% for the west orientation. These high percentages are undoubtedly the result of very low cooling energy demands due to both the location and the low emissivity of the glass; but also for Malaga, where the energy demand is higher, the average error is 8.2% and the maximum error is 16.7% for the western orientation.
Looking from annual to hourly cooling demand, once again the E-P* method appears to be the one with the highest RMSE for the south, west and north orientations with values of 28.9 Wh, 106.3 Wh and 34.0 Wh, respectively; while for the east orientation, the method with the highest RMSE appears to be the P-P with a value of 55.1 Wh (see Table 5).
The methods showing better RMSE, on the other hand, are the B-L for east, west and north orientation with values 30.4%, 46.0% and 33.0% lower than the maximum values, and the P-P method for south orientation with an RMSE 39.3% lower than the maximum values.
To conclude, calculating the maximum delta of heating and cooling requirements for each orientation of the glazed surface and for each location studied (Table 7), varying the irradiance split method leads to annual energy variations of up to 188. 8 kWh/year (21.0 kWh/year/m 2 ) for heating and up to 184.3 kWh/year (20.5 kWh/year/m 2 ) for cooling, i.e., a variation of +8.5% and −37.2% compared to the average of the energy needs obtained with the four irradiance split methods. By decreasing the window performance and switching from a low-e double glazing unit with g gl of 0.28 to a float double glazing unit with g gl of 0.63 (see Figure 4), the dev.s. between the models analysed increases for both heating and cooling requirements, except for the heating requirements for buildings with horizontal glazing surfaces, where the dev.s. between the models is lower by 2.4 kWh/year (see Table A2).
In the heating season (Table 8), for buildings with south, east, west and north glazed surfaces, the dev.s of solar loads increases compared to the case study with g gl = 0.28, with average values of 30.3 kWh/year (+3.3%), 35.3 kWh/year (+21.7%), 26.7 kWh/year (+0.3%) and 41.2 kWh/year (+54.4%). As a consequence, the dev.s of energy needs also increases compared to the case study with g gl = 0.28, with values of 78.9 kWh/year (+45.0%), 49.8 kWh/year (+56.8%), 80.2 kWh/year (41.2%) and 33.7 kWh/year (+126.4%) for the south, east, west and north orientations, respectively. Table 8. Case study with g gl = 0.63. Yearly heating and cooling energy demand obtained by averaging the results of the simulated solar irradiance datasets with the four calculation models described in Section 2 and their respective standard deviations. The blue boxes indicate the worst results for heating and the red boxes the worst results for cooling. As in the case of g gl = 0.28, it is evident from Table 8 that the south and west orientations show the highest average dev.s with errors of 5.5% (78.9 kWh/year) and 5.2% (80.2 kWh/year) on the calculation of the heating energy needs, respectively.
Focusing on the cities analysed, comparing the dev.s with the average annual winter demand estimated with EN ISO 52016-1 [4] when varying the solar radiation split methods, the average error among all orientations on the winter demand estimate is less than 5% for all locations except for Malaga, where the average error is 13.9% and the maximum error is 19.6% (south). Different from the case with g gl = 0.28, where all cities showed an error on all orientations of less than 5%, in the case with g gl = 0.63, De Bilt and Paris show an error on the annual demand estimate of 6.1% and 6.4% in the west, while Rome shows an error of 7.5% in the south.
Analysing, instead, the hourly heating energy needs, it can be seen that for all five orientations, the highest RMSE values are those obtained using as input the irradiances calculated with the E-P* method of EN ISO 52010-1 [16] with values of 11.6 Wh, 76.1 Wh, 41.0 Wh, 76.0 Wh and 37.3 Wh, respectively, for the horizontal, south, east, west and north window orientation (see Table 9). The methods that show better RMSE, on the other hand, are the B-L for the horizontal, east, west and north orientation with values lower than the maximum ones of the E-P* method of 24.0%, 34.0%, 35.7% and 64.2%, and the P-P method for the south orientation with lower RMSE than the maximum of 51.3%. In the summer season (see Figure 4 and Table A2), the dev.s of the total solar loads for buildings with south, east, west and north windows is 15.6 kWh/year, 73.4 kWh/year, 142.9 kWh/year and 26.3 kWh/year, respectively, while the average dev.s of the annual cooling energy demand for buildings with south, east, west and north windows (see Table 9) results to be 19.7 kWh/year, 75.1 kWh/year, 117.1 kWh/year and 26.9 kWh/year with an average increase over the case study of g gl = 0.28 of 434.9%, 283.6%, 289.0% and 281.1%.
Furthermore, focusing on the summer period, the variation in the decomposition and transposition models imply a variation on the cooling energy needs substantially greater than 5% for each city. In particular, the worst estimate is for the De Bilt location where the average error is 27.3% and the maximum error is 50.3% for the west orientation. These high percentages are undoubtedly due to very low cooling energy demands, but even for Malaga where the energy demand is higher, the average error is 12.7% and the maximum error is 23.8% for the western orientation.
Regarding the cooling hourly energy needs, the highest RMSEs are given again by the E-P* method with values of 19.2 Wh, 86.4 Wh, 241.9 Wh and 79.5 Wh, respectively, for the horizontal, south, west and north window orientation, while for the east orientation, the worst RMSE value is given by the P-P method equivalent to 132.9 Wh (see Table 8). The best RMSEs are obtained from the P-P models for the horizontal and southern orientation of the glazed surface, with values 27.1% and 51.4% lower than the maximum values of the E-P* method; while the B-L model gives lower RMSEs of 33.1%, 42.4% and 39.6% for the east, west and north orientations.
Calculating the maximum delta of heating and cooling demand for each orientation of the glazing surface and for each location (Table 10), varying the irradiance split method leads to annual energy variations up to 270.2 kWh/year for heating and up to 565.2 kWh/year for cooling, representing a variation of 13.9% and −52.7% compared to the average demand obtained with the four irradiance split methods. These results appear to be congruent with Michalak's study [18] where the use of EN ISO 52010-1 for the calculation of solar irradiance leads to a significant change in the energy performance indicator for heating and domestic hot water (EP H+W ) compared to that provided in the Polish TMYs, impacting the energy classification of a building. This result is shown to be very important because the improvement of the building's energy class is derived only from the choice of the calculation method of decomposition and transposition of solar irradiance and not from a real efficiency intervention.

Conclusions
Defining the boundary conditions of the conditioned environment, such as climate data, is crucial when assessing the energy performance of buildings. Most of these data (such as temperature, humidity and horizontal global irradiance) are available in special databases, while others, such as direct and diffuse solar irradiances on horizontal and inclined/oriented planes, must be estimated with special models obtained from regressions of experimental data gathered in different parts of the world, therefore not perfectly adequate to provide irradiances on every point of the earth's surface. In this study we evaluated four pairs of models for calculating the direct and diffuse irradiance of horizontal and inclined/oriented surfaces (see Section 2) and calculated the uncertainty that these models introduce in the estimation of the energy needs of a building, simulated with the simplified hourly dynamic method of EN ISO 52016-1 [4] on five European locations, also varying the solar transmission coefficient of the glazed surface (g gl ). The main results follow: The dev. of solar loads for the winter period on the various orientations and as the solar transmittance coefficient of glass changes does not cause an excessive increase in incoming solar loads, since winter solar radiation is generally low. For the summer period, on the other hand, the increase in the solar transmittance coefficient of the glass causes a significant increase in incoming solar loads on all orientations with values exceeding 178%.

•
The dev.s of the heating energy demands on the various orientations and for both case studies (ggl), show average values of uncertainty on the calculation of the total energy demand for heating below or slightly above 5%, thus acceptable.

•
On the other hand, when calculating the cooling demand, in both case studies, the largest errors are obtained for the west orientation of the glazed area, with average uncertainty values on the calculation of the total energy demand for cooling of 16.4 percent. In contrast to the heating needs, considering also the other orientations, the average errors on the estimation of the cooling demand are considerably higher than 5% and therefore not acceptable.

•
Focusing on the individual localities, while the estimation of the annual winter demand by varying the solar radiation split methods shows an average error of less than 5%, the individual city of Malaga shows errors on the winter demand estimation up to 10.6% with g gl = 0.28 and up to 19.6% with g gl = 0.63. This indicates that some of the split methods used are less appropriate to predict irradiance for the city of Malaga.

•
On the other hand, comparing the dev.s with the respective average annual summer demand when varying the radiation split methods, an error on the estimation of the energy demand that is 5% higher for all locations is obtained for both the case study with g gl = 0.28 and with g gl = 0.63. This indicates that the solar radiation split methods for estimating cooling requirements severely affect the results of EN ISO 52016-1 [4], especially for De Bilt where the average error considering all orientations is 27.3%.

•
Calculating the RMSE of both hourly heating and cooling needs and taking as baseline the results obtained using the R-P solar radiation split method contained in TRNSYS, for both g gl = 0.28 and g gl = 0.63, the model that deviates the most from these results is the one contained in EN ISO 52010-1 [16], i.e., the E-P* model.
Possible future developments will be to compare models of solar irradiance splitting the components on inclined and oriented surfaces with real data collected experimentally with pyranometers placed vertically on each main orientation. Furthermore, the experi-mentally measured solar loads will be compared with the solar loads identified with the main simulation software. Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to data from licensed software.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Table A1. Case study ggl = 0.28. Heating/cooling energy demand and solar loads obtained using the four different solar irradiance decomposition and transposition models. The asterisk * used beside "E-P*" indicates the change in EN ISO 52010 [16] from the classical method of  on the calculation of the "dimensionless clarity parameter" (see Formula (11) Table A2. Case study ggl = 0.63. Heating/cooling energy demand and solar loads obtained using the four different solar irradiance decomposition and transposition models. The asterisk * used beside "E-P*" indicates the change in EN ISO 52010 [16] from the classical method of  on the calculation of the "dimensionless clarity parameter" (see Formula (11)).