Abstract
As part of the energy design synthesis method, complex dynamic building simulation database was created with IDA ICE code for all family house building configurations for a considered problem. In this paper, the annual heat energy demand output parameter is considered to serve as basis of a building energy design investigation. The sensitivity analysis performed by Morris' elementary effect method was used. As the result of the sensitivity analysis of the output parameter, the most important input parameters can be identified, that influence the buildings' energy efficiency, that can support further building designs.
1 Family house
Design of buildings has a key role in the world's energy consumption and due to the negative environmental impacts caused by the building industry sustainable buildings frequently gain focus recently. In this relation, the process of building energy and climate responsive design represents a key factor. For energy issues see for example [1], and for nature-based solutions consider for example [2]. Unfortunately, conventional building design method in industry practice generates only one or a very limited number of concepts, based on previous experience. This is often supported by the fact, that architects consider the artistic side of design more important than complex building physics simulations and complicated mathematical models to be evaluated from the beginning of the design. Thus, when the architectural plan is ready, i.e., the building body shape together with the space organization or layout are ready as the most fundamental design features, and the technical apparatus is designed subsequently.
Energy design method [3] integrates some of these high-level engineering calculations to implement a sustainable architectural design. The method includes some heuristic building simulations, quantifying the chosen design concepts. As its extension, Energy Design Synthesis (EDS) method is a unique technique ensuring optimal buildings performing highest energy efficiency while offering best comfort. A family house geometry case was considered, where the house was constructed from 6 building blocks. First geometry generation rules were determined to generate all feasible and potentially optimal building configurations for a family house class. For these 167 configurations as accepted geometries, 5010 complex dynamic thermal simulations were performed, i.e., all the configurations were modeled and the family houses were fully generated within IDA ICE code. Please note that during the simulations two different structures including isolation as main building property were considered: one meeting the minimum standards of a family house, while the other is close to a passive house. Further, the simulations included 3 wall-window ratios and 5 orientations. From the simulation data, the annual heating energy demand as output parameter and some input parameters were selected for further investigations. Based on the results of the simulations, the most influential parameters and their dependencies were identified.
2 Sensitivity analysis
Generally, sensitivity analysis is used to identify the important parameters together with their relations, i.e., the individual importance of the selected design variables as well as their joint-effect is measured and evaluated. In other words, it is the study of how uncertainty in the output of a model can be apportioned to different sources of uncertainty in the model input. A sensitivity analysis determines the contribution of the individual design variable to the total performance of the design solution. As a result the dynamics of the variables can be investigated. Note that the effects may vary widely, therefore the order of importance is also significant. Further, individual, as well as aggregate influence occurs. For building design it is of importance, since architects will understand the underlying coherence between the input parameters of a building and its output annual energy demand. Therefore, already in the first building shape and layout design step expert decisions can be made that has direct significant influence on energy consumption sustainability.
Sensitivity Analysis (SA) methods can be classified for example in the following way:
- (i)global sensitivity methods consider multiple design parameters;
- (ii)local sensitivity methods evaluate the output variability based on the variation of a single design parameter; and
- (iii)screening methods consider usually two extreme values on both sides of the standard value for selected design parameters to evaluate, which design parameters the building performance is significantly sensitive to.
Some local SA applications consider model parameters as varying inputs, and aim at assessing how their uncertainty impacts model performance [4]. Applications often apply One-At-a-Time (OAT) approach where the sample is constructed in a way that between two sampling points the influence of only one input variable is considered while all other variables remain constant; individual influences of the output variable are examined. In other words OAT design is called when only one parameter changes values between consecutive simulations, i.e., how the output value changes in case only one input variable changed and the other variables remained or considered to be constant. On the other hand, SA applications often consider interactions between the input variables also when examining the influence of the input variables on the output. In general, global SA sampling is more often performed by All-At-a-Time (AAT) approach, where all the input factors are varied simultaneously. As a consequence, the sensitivity to each input variable considers the direct influence of that factor as well as their joint influence.
It is one of the most important elements of the models to consider the distribution of the input variables. Frequently, well-known distributions describe the input variables, for example uniform, normal, lognormal, or Weibull distribution. However, in practice many variables describing the models are not continuous and their distribution cannot be described in this manner.
Here, Morris Elementary Effect (MEE) method is used [5], when the effect of the change along the adequately scale of OAT randomly selected variable is measured. This can be well applied with a larger number of variables also as a sampling based method.
High result of the absolute expected value shows great influence of the input variable to the output. A low standard deviation means that effects are constant, namely the output is linear on this variable. High result of its standard deviation shows that there is a non-linear effect, or the interaction between the input variable and other input variables is significant. Therefore, significant parameters are depicted in the section of the
3 Applied method
The selection of the right type of the input variables is crucial from the sampling point of view. Morris method requires uniform distribution on the [0,1] interval, nevertheless in practice this cannot be fulfilled. In many cases, the variables are discrete and their distribution cannot be known. To overcome this problem, some special sampling methods are applied, and the selected method serves as the heart of the sensitivity analysis afterwards.
For example, some processes are based on the Latin hypercube sampling techniques; other are based on distance measure and uniformity. A potential solution can be the following: the levels of the parameters are given by quantiles, then in the samples the even distribution of the variables is reached; for example [8]. There the sampling is not directly done, but based on the quantile values and thus the levels of the trajectories of the Morris method appear as equal distances.
There are other processes to reach uniformity also. For example [9], which is called sampling for uniformity method, or [10], which is known as the enhanced Sampling for Uniformity method, eSU in short. This latter is used hereinafter. Note that this is a key point, since here some variables cannot be considered as continuous and in case of eSU, the sampling method defines grids that correspond to the values of the variables.
4 Model input parameters
The following total of 15 buildings' design variables was considered. As building design parameters structure, wall window ratio, and orientation were selected. Further, a set of input parameters related to the various building configurations was selected: roof surface (r), surface connected to the ground (g), balcony (b), external walls (w); together with specific edges and vertices. Further, a complex descriptor selected by the architects, namely envelope surface/Area/net floor Surface (A/S) ratio, was also considered as input variable.
5 Results
The purpose of the sensitivity analysis here was to determine the contribution of the selected individual building design variables to the total energy performance. In other words, the order of the variables with the most influence on the annual energy demand is to be determined to support efficient or optimal building design.
Here, Morris method was considered with a level of 6, i.e., for the variable models the trajectories were moved on six levels. Note that in the literature the minimum level is 4. The calculations were based on 30 trajectories.
Table 1 demonstrates the sensitivity analysis results of the full input data; i.e., the order of the considered input variables is given according to Eq. (3). Note that the table shows clearly the difference in the interpretation of the sensitivity measures discussed earlier. The corresponding unit of the investigated annual heating energy demand is kWh/a, and therefore the unit of the sensitivity measures
Order of the variables
# | Parameters | |||
1 | structure | 3,982 | −3,982 | 424 |
2 | orientation | 725 | 312 | 841 |
3 | b | 668 | 452 | 694 |
4 | a_negative_edge | 643 | 599 | 588 |
5 | g | 541 | 532 | 389 |
6 | r | 473 | 460 | 410 |
7 | g_edge | 460 | 457 | 313 |
8 | window_rate | 378 | −69 | 475 |
9 | r_positive_edge | 374 | 366 | 312 |
10 | a_positive_air_positive_vertex | 334 | 318 | 261 |
11 | r_negative_edge | 254 | 33 | 344 |
12 | g_positive_positive_vertex | 220 | 203 | 251 |
13 | a_positive_edge | 200 | 4 | 269 |
14 | a_per_s | 107 | 103 | 131 |
15 | w | 98 | 12 | 146 |
Table 2 demonstrates the sensitivity analysis results of the input data related to structure 1, where the minimum standards were considered as the necessary building properties, while in Table 3 is demonstrated the sensitivity analysis results of the input data related to structure 2, close to a passive house’ properties. The order of the considered input variables is given according to Eq. (3) in both cases. It is important to mention that there is only a slight difference in the order as well as in the data of the results, which can be considered as not significant and may be a result of the random generation steps. Obviously, the structure and orientation are the most important parameters. Then, variables describing the various building configurations follow for every separate study.
Order of the variables for structure 1
# | Parameters | |||
1 | orientation | 834 | 322 | 1,005 |
2 | b | 715 | 515 | 715 |
3 | a_negative_edge | 663 | 609 | 622 |
4 | g | 558 | 556 | 391 |
5 | r | 462 | 454 | 395 |
6 | g_edge | 458 | 446 | 333 |
7 | window_rate | 412 | −371 | 390 |
8 | r_positive_edge | 403 | 388 | 342 |
9 | a_positive_air_positive_vertex | 343 | 324 | 267 |
10 | r_negative_edge | 286 | 98 | 359 |
11 | g_positive_positive_vertex | 221 | 201 | 260 |
12 | a_positive_edge | 203 | 80 | 267 |
13 | a_per_s | 113 | 108 | 144 |
14 | w | 98 | 16 | 144 |
Order of the variables for structure 2
# | Parameters | |||
1 | orientation | 708 | 290 | 844 |
2 | a_negative_edge | 639 | 603 | 557 |
3 | b | 638 | 410 | 655 |
4 | g | 583 | 582 | 435 |
5 | r | 491 | 483 | 422 |
6 | g_edge | 445 | 445 | 326 |
7 | r_positive_edge | 352 | 345 | 286 |
8 | window_rate | 319 | 258 | 293 |
9 | a_positive_air_positive_vertex | 298 | 274 | 230 |
10 | r_negative_edge | 272 | 26 | 361 |
11 | g_positive_positive_vertex | 207 | 197 | 212 |
12 | a_positive_edge | 190 | −39 | 257 |
13 | a_per_s | 93 | 90 | 113 |
14 | w | 84 | −26 | 120 |
The
Figure 3 corresponds to monotonicity, the variables on the
6 Conclusions
The sensitivity measures
Acknowledgements
This work was partially supported by the 2019-2.1.11-TÉT Bilateral Scientific and Technological Cooperation.
Many thanks for Drs. Yogesh Khare and Rafael Muñoz-Carpena for development of the Matlab/Octave packages for Morris' Elementary Effect methods using enhanced Sampling for Uniformity.
References
- [1]↑
M. S. Albdour, B. Baranyai, and M. M. Shalby, “Overview of whole-building energy engines for investigating energy-related systems,” Pollack Period., vol. 18, no. 1, pp. 36–41, 2022.
- [2]↑
Q. He and A. Reith, “Using nature-based solutions to support urban regeneration: A conceptual study,” Pollack Period., vol. 17, no. 2, pp. 139–144, 2022.
- [3]↑
B. Baranyai, B. Póth, and I. Kistelegdi, “Planning and research of smart buildings and constructions with the ‘Energy design Roadmap’ method,” Pollack Period., vol. 8, no. 3, pp. 15–26, 2013.
- [4]↑
F. Pianosi, K. Beven, J. Freer, J. W. Hall, J. Rougier, D. B. Stephenson, and T. Wagener, “Sensitivity analysis of environmental models: A systematic review with practical workflow,” Environ. Model. Softw., vol. 79, pp. 214–232, 2016.
- [5]↑
M. D. Morris, “Factorial sampling plans for preliminary computational experiments,” Technometrics, vol. 33, no. 2, pp. 161–174, 1991.
- [6]↑
F. Campolongo, J. Cariboni, and A. Saltelli, “An effective screening design for sensitivity analysis of large models,” Environ. Model. Softw., vol. 22, no. 10, pp. 1509–1518, 2007.
- [7]↑
D. G. Sanchez, B. Lacarrière, M. Musy, and B. Bourges, “Application of sensitivity analysis in building energy simulations: Combining first- and second-order elementary effects methods,” Energy Build, vol. 68, pp. 741–750, 2014.
- [8]↑
A. Saltelli, M. Ratto, T. Andres, F. Campolongo, J. Cariboni, D. Gatelli, M. Saisana, and S. Tarantola, Global Sensitivity Analysis. The Primer. John Wiley & Sons, Ltd, 2007.
- [9]↑
Y. P. Khare, R. Muñoz-Carpena, R. W. Rooney, and C. J. Martinez, “A multi-criteria trajectory-based parameter sampling strategy for the screening method of elementary effects,” Environ. Model. Softw., vol. 64, pp. 230–239, 2015.
- [10]↑
J. Chitale, Y. Khare, R. Muñoz-Carpena, G. S. Dulikravich, and C. Martinez, “An effective parameter screening strategy for high dimensional models,” in Proceedings on ASME International Mechanical Engineering Congress and Exposition, vol. 7, Tampa, Florida, USA, Nov. 3–9, 2017, Art no. V007T09A017.