1. Introduction
Since the early 1970s, building energy performance simulation (BEPS) tools have played a crucial role in the assessment of the environmental and energy-related impacts of buildings, which represent the largest energy consuming sector in the developed world [
1]. Specifically, due to the growing concern about the energy efficiency of buildings, a significant number of simulation tools has been developed and adopted to support efficient design or redesign, construction or refurbishment, and operation of energy efficient buildings. Such tools also allow dealing with all the most important phenomena occurring in the building itself [
2,
3,
4] and among it and the applied energy efficiency techniques or renewable energy sources [
5]. As a result, dynamic simulation is a well-recognized method for the assessment of the building energy performance. The use of dynamic simulations is particularly crucial during the very early stages of the building design, as well as during the refurbishment process of existing buildings (where full scale testing campaigns are costly and/or not practicable). The adoption of BEPS tools is more and more crucial and highly recommended for the design of the next generation of buildings (NZEBs), for the promotion of the integration of renewable energy technologies and for the building energy diagnosis [
6,
7].
Although BEPS tools have been in use throughout the building energy community for decades, their adoption has been boosted by the recent advances in computational methods and computer calculation power, that have provided key opportunities also for the enhancement of such tools and for the development of new ones [
8]. In fact, despite of the availability of commercial BEPS tools and codes, more and more up to date tools (mainly conceived for research applications) have been developed (e.g., [
9,
10,
11]) with the intention of providing greater design flexibility while accomplishing the challenging multiple aim of [
12]: (i) investigating new unreleased energy efficiency measures and strategies, (ii) meeting the highest expectations for the comfort of occupants, (iii) increasing models’ robustness and fidelity also toward the implementation of innovative control strategies [
13]. Moreover, novel BEPS tools, also coupled to commercial ones, are constantly being developed to help the definition and development of the next generation of buildings [
5], to be designed through suitable computer-based analyses [
14,
15].
Despite the number of new developed tools and codes and the wider and wider adoption of dynamic BEPS tools, building energy simulation is still nowadays a complicated process that requires modeling and analytical skills for both the development and use phases. As an example, the use of BEPS tools and the analysis of the related results can be considered as a challenge for building designers and practitioners who are sometimes bewildered by the selection of a suitable BEPS tool and by the uncertainty about the fidelity of the calculation results [
12]. Therefore, novel in-house developed BEPS tools must be validated in order to ensure reliable and accurate energy analyses and to prevent a negative impact on the acceptance of the obtained results. Nevertheless, as well known, despite of the advancements in validation methodologies for building energy simulation tools, the validation process is still time consuming and rather difficult to accomplish [
16]. It is also worth noting that while the validation of a single mathematical model for the energy performance analysis of a new building technology, component or material, can be often easily assessed, the experimental validation of a whole BEPS code is not always feasible and currently practicable. This is particularly true for the new generation of dynamic building simulation codes that must be capable of simulating a variety of innovations which would require extensive extremely costly and time consuming testing procedures [
5,
17]. Note that, accurate energy analyses often could not be ensured by simply adding to existing BEPS codes new subroutines, suitably developed for simulating innovative technologies and investigating their effectiveness. In any case, at the end of the entire modeling process a validation of the whole BEPS code is always necessary and recommended [
12].
In this regard, in order to validate novel BEPS tools, several general criteria and standard procedures have been developed and are available in [
18]. They consist of comprehensive and integrated suites of building energy analysis tool tests, involving analytical, comparative, and empirical methods. Their adoption has been recently emphasized by the energy performance building directive (EPBD) issued by the European Union that underlines the need of new certified tools to be developed for decision-makers and practitioners in order to provide integrated building design applications and to obtain the compliance with higher energy efficiency standards.
One of the most important international standard tools conceived in order to increase the confidence of BEPS code outcomes is the Building Energy Simulation TEST (BESTEST) procedure [
19]. This validation method was purposely developed for testing, diagnosing and validating the capabilities of new BEPS tools. The BESTEST procedure consists of suitable analytical techniques and tests to be carried out for the comparison of the results obtained thorough a novel BEPS program or design tool with those obtained by the current state-of-the-art codes. Recently, a number of the BESTEST cases have been incorporated into ANSI/ASHRAE Standard 140-2011 [
20]. Note that, the International Energy Agency (IEA) assumed the BESTEST as the official standard procedure for the new BEPS codes validation. The BESTEST procedure was also adopted by the European Committee for Standardization (CEN) as the test for checking the reference cooling load and energy calculation methods, based on the requirements of several standards addressing different aspects of the EPBD [
21]. Through the BESTEST and similar procedures, many new BEPS codes were validated [
22,
23,
24].
Although different mathematical models and assumptions adopted in the available BEPS codes (e.g., for assessing the occurring building physic phenomena), by following the above mentioned validation approach all the occurring bugs can be detected and fixed. At the end of the validation procedure, reliable predictions of the system behavior can be obtained in terms of: (i) energy use for heating, cooling, lighting, etc.; (ii) operating temperatures of the building elements and indoor air; (iii) occupants’ indoor comfort.
This paper is focused on this framework. In particular, it describes the most important steps and outcomes obtained by following the BESTEST validation procedure, purposely adopted for the diagnosis and validation of a dynamic BEPS tool, named DETECt and developed for research aims by Buonomano and Palombo [
25]. This code is able to dynamically predict the thermal behavior of mono and multi-zone (and multistory) buildings and to assess the benefits of different and advanced building envelope techniques in case of different weather locations, envelope materials (also including integrated phase change materials and solar based technologies), building shapes, orientations and geometries. DETECt was implemented in a computer code written in MatLab [
26]. The tool and the code are not yet available for use and inspection by any users. Nevertheless, since DETECt was expressly developed for research aims, it has been used and inspected by several researchers (with different backgrounds) for the implementation of novel features. From this point of view, by means of DETECt, novel building envelope technologies or plant/systems and control strategies were implemented and studied by the authors and by several researchers [
27,
28,
29]. The authors are also working on a future freely available release of the tool with a suitable user friendly interface.
The BESTEST procedure adopted for the validation of DETECt (e.g., in case of the basic mono-zone building) was selected on the basis of several criteria, such as completeness, accuracy, reproducibility and cost effectiveness of the test suite. In fact, this procedure consists of widely available high quality empirical validation data sets and it includes detailed and unambiguous documentations of the input data for a selected number of representative design conditions.
The outcomes of this paper refer to several main topics, such as: (i) comparison of DETECt results vs. the mandatory BESTEST qualification cases (in terms of heating and cooling annual energy demand and integrated peak-load, annual transmitted and incident solar radiation, annual hourly zone temperature); (ii) results and discussion regarding the assessment of some hourly variables and their comparison vs. the reference trends; (iii) description of the iterative refinement process carried out for minimizing the sources of the divergences, detected during the validation tests and occurred among the results obtained by the presented code and the reference ones.
DETECt surpassed the BESTEST reliability test. In particular, low or very low deviations of the results obtained by the novel presented code
vs. those provided by the validation procedure are observed. Note that an additional code to code validation test, purposely conceived for several commercial buildings located in different weather locations, was also carried out [
25]. The development of the presented model and the correspondent code validation were obtained after several years of research work and the activity of modeling innovative technologies for building envelopes and HVAC systems, to be implemented in DETECt, is still today extant.
In the paper, some original examples concerning the capability of the presented code to assess the building energy performance are also reported. Some graphics about the temperature and heat fluxes dynamic trends and spatial temperature gradients into walls are also reported.
3. Innovative Capabilities of DETECt
DETECt has been developed for research purposes with the aim of assessing the energy performance of innovative building integrated technologies and configurations, and particularly to analyze their impact on the building energy demands. The development of new in-house building dynamic simulation models, capable to simulate innovative building energy efficient measures, is well encouraged in literature and supported by the recent advances in analysis and computational methods and calculation power [
8]. In fact, despite of the availability of commercial BEPS codes (with high level of flexibility and complete user interfaces/data libraries), new in-house simulation models are developed for research aims and, more in general, for implementing new unreleased energy efficiency measures, which are particularly necessary and recommended for the design of the next generation of buildings, such as NZEBs [
5,
14,
15].
An innovative building design requires the suitable integration of different innovative elements, whose energy effectiveness must be preventively assessed in order to avoid unexpected/ undesired performances and the design failure [
4]. In this regard, computer-based analyses are well recognized to support the understanding of complex systems (such as NZEBs) by providing feedbacks on the performance of all the investigated innovative design scenarios [
12]. Based on high order RC networks, several new dynamic simulation models have been developed for analyzing, as an example, the impact on the building energy demand of glazed/opaque surfaces and internal solar radiation [
50], the thermal behavior of multi glazing façades [
51], the thermo-physical behavior of a Phase Change Material (PCM) layer in combination with transparent and opaque materials [
52], the performance of Building Integrated Solar Thermal and Photovoltaic systems (BIST-BIPVT) [
53]. New BEPS computer-based models are also required to be flexible and suitable for the calculation of household actions and to effectively account model occupants’ behavior [
54]. Nevertheless, the available literature remarks a lack of knowledge regarding the availability of studies focused on the energy performance obtained through the building integration and combination of such innovative technologies [
55,
56]. Thus, an additional research effort is still required in order to obtain detailed dynamic simulation models for the energy analysis of the whole building-plant system [
57]. In this regard, very few studies available in literature deal with: (i) building integrated PV systems coupled to PCM panels; (ii) passive and active effects on the overall building energy performance of BIST and BIPV systems, also coupled to PCMs wallboards [
58].
Therefore, DETECt was developed and has been enhanced with the aim of obtaining a reliable tool capable of carrying out comprehensive whole building dynamic simulations while assessing the energy interaction of the abovementioned building integrated technologies and novel strategies. It is also worth noting that, although commercial and open source analysis tools (e.g., TRNSYS, Energy-Plus,
etc. [
59]) or generic programming environments [
9,
60] are regularly adopted among building designers and researchers for swift numerical analyses, such programs are less flexible with respect to the controller development and the swift implementation of specific mathematical models (e.g., for describing the thermal dynamic of innovative materials or technologies by means of whole building simulation environment). DETECt has been conceived in this framework, with the aim to fill the above mentioned lack of knowledge regarding the need of comprehensive whole building energy performance analyses for innovative building integrated measures for energy saving. Examples of the use of DETECt to analyze innovative building integrated technologies and implemented control strategies (also toward the NZEB goal) are reported in [
27,
28,
49].
In the following, a short description of the main assumptions about the simulation models related to the innovative capabilities of DETECt is reported. Such models are developed in order to perform comprehensive and whole building energy performance analyses, by also properly taking into account the passive and active effects due to the building integration of innovative technologies. Specifically, with the help of DETECt it has been possible to:
Investigate the application of both single solid phase and phase change materials (PCM), for opaque and glazing elements [
28]. Thermo-physical properties of the building envelope (
i.e., of the RC network) can be considered as constant (in case of single solid phase materials) and time-dependent (in case of phase change materials, as a function of the layers temperature (effective specific heat capacity model technique [
61]).
Take into account the mutual energy interaction between building integrated PV and PVT systems (BIPV and BIPVT) and indoor building spaces. The tool also enables the simulation of BIPV and BIPVT coupled to PCM sub layers. Such systems can be integrated into the building envelope, such as the roof or façade, and/or in parts of them (differently from some commercial codes, where BIPV and BIPVT often cannot be simulated as only partly integrated in the simulated surfaces) [
28].
Simulate the thermal behavior of attached sunspaces and complex glazing systems, by accurately taking into account the beam solar radiation that, passing through them, enter the interior (and adjacent) thermal zones [
28].
Meet the building daylighting and shading requirements by controlling the use of artificial lights and windows solar shadings with the aim to optimize the visual comfort, the heating and cooling demands, and the need for artificial lighting. To this aim, the adopted simplified daylighting model (based on the lumen method [
62]) allows one to simulate the switching on/off of artificial lights and the modulation of external solar shadings (obtained by varying their shading coefficients) which can be optionally scheduled during any simulation time interval [
28].
Implement advanced occupant schedules models in order to accurately take into account the occupants dynamic behavior and their interaction with the building (opening of windows solar shadings, variation of natural ventilation rates, switching on/off HVAC plant, etc.).
Implement advanced controllers (e.g., based on the online adaptation of the control gains) for the optimal calculation of the sensible and latent energy demands [
27].
Calculate the energy requirements and thermal-hygrometric comfort of multi-enclosed thermal zones, featuring special environments (e.g., indoor hospitals spaces (including multiple infant-incubators), laboratory temperature control rooms,
etc.) [
49].
Perform sensitivity and/or parametric analyses with a single (so called “one-click”) simulation run [
28], remarkably useful to: (i) assess the energy effectiveness and the economic and financial implications of innovative building features and technologies, also toward the NZEB goal; (ii) provide additional insight on specific parameters influencing the building performance, especially at early design steps and in case of non-standard solutions.
4. Model Validation
The aim of the paper is to present the steps followed to complete a code-to-code validation of DETECt. As it is well known, there are only a few ways to check the accuracy of a whole-BEPS model [
63]: empirical validation, analytical verification and comparative testing. The difference among such approaches depends on the way in which the calculated outputs or results obtained by a program, subroutine, algorithm or software object are compared to some data considered as reference ones. By the empirical validation, the calculated results under exam are compared to the measured data obtained from a real building, test cell, or laboratory experiment. Through the analytical verification procedure, the calculated output data are compared to the results achieved through a standard analytical solution or a generally accepted numerical method for simulating the heat transfer mechanisms under very simple and highly constrained boundary conditions. By means of the comparative test procedure, the calculated results are compared to those provided by several reference programs that may be considered validated or more detailed and, presumably, more reliable than the code under test [
64]. Details about advantages and disadvantages of these procedures are reported in [
63,
65].
In this paper in order to check the DETECt reliability, the comparative test procedure (code-to-code validation) is used. The reason for this choice is based on the unavailability of an experimental set-up. Several detailed code-to-code validation procedures, developed in the last years, are available in literature. In general, they incorporate several different test cases for analyzing building models accuracy, under different operating conditions. From this point of view, the IEA commissioned a number of projects to develop proper validation methodologies for building energy models [
30,
63,
66,
67]. Among them, the BESTEST procedure is the one adopted and utilized for assessing the accuracy of the presented simulation code.
4.1. BESTEST Validation Procedure
The BESTEST validation procedure includes several test cases, organized in diagnostic and qualification (mandatory) series, which allow the user to analyze the influence of different physical phenomena on the numerical results of building energy models, covering a high fraction of the physical phenomena occurring within the built environment. The built tests for the iterative diagnostic procedure allow codes to be examined over a broad range of parametric interactions based on a variety of output types, minimizing the possibility of compensation errors (which cannot be easily assessed by experimental campaigns, often time consuming or unfeasible) [
19]. From this point of view, by following most of the tests each upon other, it is possible the check the effect in the calculation results of many model features, such as thermal mass, solar and internal heat gains, window-shading devices, infiltration and setback thermostat control.
The features of the BESTEST reference building respond to a compromise between those of U.S. and European constructions. A simple rectangular plant geometry and a single thermal zone with no interior partitions are taken into account,
Figure 4. Although the test cases explore different combinations of parameters and settings, the BESTEST results are related to a single annual hourly Typical Meteorological Year (TMY1) weather data file, related to the Denver, CO (USA). Its climate features a cold clear winter and a hot dry summer, with large diurnal temperature variations throughout the year. In order to check the accuracy of a simulation code under exam, for each case test the BESTEST procedure provides several results obtained by means of several codes considered as reference ones. Such programs represent, for the building energy performance analysis, the best state of the art [
19]. The results provided by these reference codes delimit the reference range of the case test and the confidence interval in which the results of the simulation code under test have to be included. Note that the existence of such reference range means that no truth standard results exists in building energy performance simulations [
19].
DETECt has surpassed the BESTEST reliability test. About three years were necessary for the model development and the related optimization and diagnostic processes. During such period a high number of simulations was carried out in order to satisfy the BESTEST validation procedure. In particular, by following the BESTEST diagnostic procedure and the test cases logic method [
19], the algorithms or existing bugs causing discrepancies have been, step by step, fixed and eliminated. The complete code calibration was also achieved by iteratively improving the calculation procedure of some model subroutines. From this point of view, in the early version of the presented code, the divergences found between the DETECt results and the BESTEST ones were mainly due to the oversimplified handling of the solar radiation effect on the interior and exterior building surfaces. Note that through this diagnostic procedure it was possible to ungroup the particular physical effects simultaneously simulated by means of different code algorithms, detecting the related errors that in a standard simulation process could balance themselves. Moreover, such tests also allowed checking the code against itself, with the aim of analyzing the model sensitivity by varying the operating parameters taken into account in the mathematical algorithms.
4.2. BESTEST Test Cases and DETECt Assumptions
In order to definitely validate DETECt, comparisons with the mandatory BESTEST qualification cases were necessary [
19]. Through such qualification cases, the DETECt reliability was checked by changing the building features and several parameters, such as direct solar gain through windows, window orientation and shadings, heating and cooling running hours, set point temperatures, building air change per hour (ACH), the set-back thermostat control.
In such tests the energy performance of a basic building envelope (
Figure 4), modeled by taking into account two different inertial masses, is dealt with. Thus, two series of qualification cases are taken into consideration for the light (cases from 600 to 650) and heavyweight (cases from 900 to 950) building [
19]. The BESTEST building geometry for both qualification cases 600 and 900 is a simple rectangular plant insulated room with two south facing windows. All the related geometric features are shown in
Figure 4a. Note that, the two south windows are considered as an equivalent single one. With respect to such cases, in cases 610 and 910 a unique horizontal overhang of 1 m width above the windows is added to the building,
Figure 4b. Shifting to cases 620 and 920 the two windows are modeled on the west and east building façades,
Figure 4c. With respect to such cases, in cases 630 and 930 for each window a horizontal overhang and two lateral vertical fins of 1 m width are taken into account,
Figure 4d.
The thermophysical properties of the BESTEST building materials are reported in
Table 1. In
Table 2 the envelope layers and the material thicknesses are reported. All the features of the mandatory qualification test cases for both light and heavy building are summarized in
Table 3. For the calculation of the window solar transmission (depending on the angle of direct solar radiation), a solar heat gain coefficient at normal incidence equal to 0.787 is taken into account. For all the mandatory BESTEST cases, the exterior and interior emissivity of the opaque surfaces are set to 0.9, while their short wave absorptance is set to 0.6. Glass infrared emissivity and absorption coefficient are set to 0.84 and 0.10, respectively. More details about such BESTEST cases are reported in [
19].
The main assumptions adopted in DETECT are described. Specifically, for opaque and transparent surfaces, the internal surface unitary convection heat transfer coefficient is assumed to be equal to 3.16 W/m
2·K for horizontal heat transfer on vertical surfaces and to 4.13 and 1.00 W/m
2·K for upward and downward heat transfer on horizontal surfaces, respectively. The external surface unitary convection heat transfer coefficient is calculated by assuming a wind velocity of 5 m/s in case of opaque surfaces and it is set to 16.37 W/m
2·K in case of windows. Note that, eventual corrections of glass surface coefficients according to window overhangs and fins are disregarded. A combined radiative and convective coefficient of the window air gap of 6.297 W/m
2·K is assumed. For the investigated test cases, the DETECt simulation sample time is set to 300 s in order to obtain satisfactory dynamic building behavior and suitable system peak loads [
68]. The number of nodes in which each building element is split range from 10 to 30 depending on the related weight, according to [
69]. The heavier the mass element, the higher the number of nodes. Finally, in order to assure the result stability and the solution convergence, the tolerance and maximum time step solver parameters are set to 0.015 and 60 s, respectively. Here, the integration interval is set to 1 s. All the numerical DETECt results are obtained by the MatLab ode23t ODE solver [
39].
4.3. Comparison with the BESTEST and Code Calibration Discussion
In the following, the main results obtained through the BESTEST validation procedure are summarized. Note that about thousand dynamic simulations were necessary for completing the calculation of the heating and cooling demands and temperature profiles of all the considered buildings. For all the BESTEST cases, with the exception of those in free floating regime, the comparison results related to the yearly heating and cooling energy requirements are shown in
Figure 5 and
Figure 6, respectively.
For the same cases, the obtained heating and cooling peak loads and their comparison
vs. the BESTEST ones are reported in
Figure 7 and
Figure 8, respectively. As it is possible to observe in such figures, the differences between the results achieved by means of DETECt are always included in the confidence intervals allowed by the considered validation procedure. The BESTEST lower and upper limits reported in such figures delimit, case by case, the provided reference results ranges. The DETECt heating and cooling energy demands and loads, for all the given cases, are always included in the allowed confidence intervals.
In the following, a discussion about the refinement iterative process carried out for minimizing the sources of the detected differences, occurred during the validation tests, is reported. At the beginning of such procedures these differences were mostly due to the model algorithms related to the assessment of the solar flux absorbed by the building outside surfaces and distributed inside the indoor space. From this point of view, for the sake of simplicity, in the first versions of the presented code, no distinction between convective and radiative heat transfer on the inside and outside surfaces was assumed. Such phenomena were successively modeled through different and distinguished algorithms. In a second version of the code the solar flux transmitted by the windows toward the indoor zone was still assumed to be absorbed only by the floor. As a result of all such assumptions, an initial remarkable disagreement of the results vs. those of the reference programs was obtained. Note that, such discrepancy was detected and fixed by simulating different additional BESTEST diagnostic tests, such as case 270 (whose description and results are not reported here for the sake of brevity). Thus, in a subsequent code version, such radiation was modeled by directing it only toward the floor and partly reflecting it toward the other interior surfaces, which absorb the incoming solar radiation as a function of their absorption coefficients. Although such improvements, the long wave internal radiation was still disregarded in the model. Such simplification resulted in remarkable discrepancies particularly for the heavyweight buildings cases. In fact, the validation process showed how important is the effect of the solar gain when the building zone is characterized by thermally massive walls storing the heat gained by the absorbed solar radiation.
For these reasons, in the ultimate DETECt version, internal long wave radiation heat exchanges and multiple indoor surfaces reflections, considered as diffuse, were taken into account [
70]. However, after this enhancement, both the low and high mass qualification tests for shaded cases showed light residual results differences caused by the solar radiation transmission on the east and west window surfaces. This problem was linked to the modeling of the horizontal overhang and the vertical fins located around the windows (
Figure 4d), disregarding the shadow on the adjoining opaque wall surface. For the cases 630 and 930,
Figure 6, the results obtained for the annual cooling requirements are close to the related minimum range bandwidth. Therefore, such results could not be linked to the disregarded estimation of the shades opaque surfaces. After several tests, such results were finally linked to the calculation of the transmitted annual solar radiation on the east and west building surfaces. For the shaded cases, such incoming radiation is resulted lower than the minimum one yielded by the reference programs,
Table 4. Here, for the unshaded test cases, the annual transmitted solar radiations (diffuse and direct) through west- and south-facing windows (concerning cases 620 (or 920) and 600 (or 900), respectively) are reported.
In the same table, for the same oriented surfaces, the diffuse and direct shaded annual transmitted solar radiations (concerning cases 630 (or 930) and 610 (or 910)) are also shown. Note that, concerning the transmitted solar radiation, such good agreements have been obtained by means of a variable direct shortwave transmission coefficient of windows. This, in the first version of the code, was kept constant without accounting for the hourly sunlight incident angle dependence. In
Table 4 the annual incident solar radiations on all the vertical and horizontal surfaces are also reported. In this table a good agreement among the DETECt results and those provided by the BESTEST procedure is always achieved.
As written above, the BESTEST validation procedure requires the assessment of some hourly variables and their comparison
vs. the reference trends [
19]. In particular, the calculated hourly free floating temperatures trends and their annual minimum and maximum values must to be checked. From this point of view, in
Table 5 the annual minimum and maximum temperatures for the free floating cases (cases 600FF/900FF and 650FF/950FF) are reported. In this table, it is possible to observe that almost all the DETECt results are included in the BESTEST allowed ranges. For the same cases, the comparison
vs. the BESTEST data of the simulated free floating temperatures trends are reported for the winter day, selected for the BESTEST procedure, in
Figure 9a and for the summer one in
Figure 9b. Here, the calculated profiles are resulted always included between the minimum and maximum reference ones. In
Figure 10a for case 600 and in
Figure 10b for case 900 the heating and cooling requirements are reported for the winter day selected by the BESTEST procedure. Also in this case the calculated profiles are always included between the minimum and maximum reference ones.
Note that, according to the BESTEST results no cooling energy demand is required by the heavyweight building,
Figure 10b. This is also due to the considered cooling set point (27 °C,
Table 3). Comparing the calculated free floating temperature trend (
Figure 9a) to the heating and cooling requirements one of the same winter day (
Figure 10a), it is possible to observe that only for the lightweight case no correspondence between the related hourly values is obtained. In fact, to the calculated temperature close to the reference upper limit should correspond a calculated heating energy requirement close to the reference lower limit.
Such an outcome, probably due to the adopted thermostat controller, was also noted in the results achieved by some BESTEST reference programs. Since the calculated hourly energy requirements and temperature trends for the considered days, reported in
Figure 9 and
Figure 10, agree well with the reference results, concerning the remaining results no additional problems regarding the thermal dynamics (as well as the thermostat controller) can be observed.
5. Outputs
In the following some examples about the potentiality of the presented code are reported. In particular, for some of the above described test cases the time history and the spatial trends of some thermophysical properties useful in a building energy analysis are reported.
In
Figure 11,
Figure 12,
Figure 13,
Figure 14 and
Figure 15 some dynamic outputs obtained through DETECt are reported. By means of such a graphical tool the main thermal phenomena occurring in an investigated building can be detected. In addition, design and operating guidelines can be easily obtained through each energy performance simulation. In particular, in
Figure 11 for the test cases 600 (
L: lightweight) and 900 (
H: heavyweight), by switching off the HVAC system (free floating, FF, regime), the dynamic trends of the indoor air temperature (
Tin,FF), the internal surfaces temperatures of the window glass (
Tgl,FF) and the roof (
Trf,FF) are depicted for several summer days (August 21st–25th).
In
Figure 12 for the same summer days and the same internal surfaces and indoor air temperatures of
Figure 11, the dynamic profiles of the free floating regime are depicted for the test cases 620 (
L: lightweight) and 920 (
H: heavyweight). In the same graphs the profiles of the outdoor air temperature (
Tout) and of the beam solar radiations incident on the south vertical wall (
),
Figure 11, and on the east and west vertical walls (
and
),
Figure 12 are also shown. In the high left side of these figures the yearly dynamic trends
Tin,FF, for both light and heavyweight envelopes, are also reported. Note that, in the case of free floating regime, the thermal fluxes through walls are directed from inside to outside, most of the time. In
Figure 11 and
Figure 12 the thermal inertia effect of the building envelopes can be observed by means of the peaks of
Tin and
Trf, which are always delayed with respect to the
Tgl ones. In
Figure 11, the
Tgl are obviously the closest to the peaks of
. Note that for the lightweight envelopes, higher interior and indoor temperature are always achieved by the building with facing south windows
vs. that one with windows located on the east and west façades. Obviously, the indoor air temperature fluctuation during the whole day is more remarkable for the lightweight envelope building. Note also that, without switching on the HVAC system, during the daytime the indoor air temperatures are far (
H: heavyweight) or very far (
L: lightweight) from the comfort ones (
Tin,FF always higher than 27 °C). During the night time, such occurrence is detected only for the lightweight envelope building (
Tin,FF lower than 20 °C). In the investigated summer days, for both the envelope cases the indoor air temperature is always much higher than the outdoor one. At last, by means of such graphical tools the average radiant temperature of the investigated thermal zone can be easily assessed.
In
Figure 13 and
Figure 14, in case of switched on HVAC system, for both the light (
L) and heavy (
H) envelopes of test cases 600/900 and 620/920, the temperature dynamic trends of indoor air (
Tin) are shown. Such profiles are referred to the HVAC system activation beyond the set point temperatures of 20 and 27 °C for the same previously accounted summer days of the FF regime figures.
In the same figures, the profiles of
Tout and of the total solar radiations on vertical south (
),
Figure 13, and east and west (
and
),
Figure 14, external surfaces of the buildings are also reported. The dynamic trends of the related sensible cooling heat flux (
) are shown too. In the high left side of such figures the dynamic profiles of all these parameters are reported for the whole year.
is obviously accounted during the
Tin control, which continually occurs for the whole day. Note that the cooling energy demand is higher for the lightweight envelope building and the heating load is obtained only for the lightweight one. The short time lag between
(
L) and
(
H) is due to the global low envelope weight difference. Note that the light and heavy wall (floor) weights are 17 and 145 kg/m
2 (26 and 122 kg/m
2), respectively, whereas the same roof weight of about 21 kg/m
2 is simulated for both the light and heavy envelopes.
In
Figure 15 on the left side for cases 600FF and 900FF, respectively, the dynamic profiles of the simulated internal and external roof surfaces,
Trf,int and
Trf,ext, respectively, and the dynamic profiles of the indoor air,
Tin, temperatures obtained through DETECt and a reference simulation code, REF, are reported for a summer day (June 12th). In the same figures the
Tout profile is also depicted. Here, as reference software, for comparison purposes, was assumed TRNSYS 17, considered as a BEPS tool well representative of the recent state of the art. Note that, the results of the BESTEST cases achieved by using such newer tool version can show significant differences
vs. the original previous one included in the BESTEST procedure. This is due to the code improvement (as consequence of the remarkable developments in building sciences and technologies achieved from 1980s) and to the user interpretations of standards and simulation assumptions during the BESTEST cases reproduction [
71]. Concerning the comparison between DETECt and REF dynamic profiles, a good agreement is observed. The reference software seems to slightly underestimate the
Tin assessment. In particular, it was verified that the calculated heating and cooling demands and peaks obtained through TRNSYS 17 for the BESTEST cases 600 and 900, differ from those calculated by means of previous TRNSYS versions [
72]. In addition, whereas the heating requirements obtained through TRNSYS 17 still comfortably fall within the reference present BESTEST acceptability ranges, the cooling ones slightly fall below the correspondent confidence lower limits. At the right side of
Figure 15 the spatial temperature gradients across the roof stratification are reported for some time intervals of the day and for both light and heavyweight envelopes. On the
x-axis of such graphs the considered thickness between the interior and exterior surfaces is reported. It must be noted that for both light and heavyweight buildings, the material and thickness of the roof layers are the same. Here, an incoming heat flux through such building component is detected for most of the time. In particular, although the very short time lag between indoor air temperatures of light and heavyweight envelopes, an outgoing heat flux from the building through the roof occurs only from 7:00 to 13:00 for case 600FF and only from 12:00 to 17:00 for case 900FF. Note that higher temperature gradients are obtained across the central insulation layer, as expected.
6. Conclusions
In this paper the main features of a dynamic building energy performance simulation model, purposely developed for research purposes [
25], are described. The building simulation model, called DETECt, is capable of predicting the dynamic thermos-hygrometric behavior of mono- and multi-zone, as well multi-story buildings and assessing the benefits of different and advanced building envelope integrated techniques (e.g., building integrated solar thermal systems, phase change materials,
etc.), solar gain controls and advanced daylighting strategies in case of different weather locations, envelope materials, building shapes, orientations and geometries. The code is also able to carry out suitable parametric analyses, which can be easily performed through a single iterative simulation run, without the need of re-entering, step by step, the building models and the features to be varied.
The steps and results of the validation procedure, based on the BESTEST, carefully followed to validate the presented dynamic simulation model, are here described. A brief description of the mathematical models and algorithms are also reported in order to provide insights about the engine logic and to show the model assumptions that primarily influence the fidelity and accuracy of the simulation results. Specifically, the physical building model implemented in DETECt is based on a high order resistive-capacitive (RC) thermal network obtained through distributed parameters. Capabilities and outputs of the code, including few application examples, are finally presented.
DETECt exceeded the carried out validation procedure. The simulation results obtained through DETECt, compared to those provided by the BESTEST, are always included within the provided confidence intervals. In order to show the potentiality of the simulation code, the most important outputs referred to different BESTEST cases are reported and suitably discussed. They refer to heating and cooling demands, calculated over the whole year or a selected time horizon, heating and cooling loads, temperatures of indoor air, building elements surfaces and nodes.
DETECt allows one to carry out comprehensive and whole building energy performance and comfort analyses, by properly taking into account the passive and active effects due to the building integration of innovative building integrated technologies, such as BIPVT-BISTS, PCM,
etc. [
28]. The code also includes a graphical tool for plotting the dynamic profiles of building temperatures and heat fluxes and the time-variant spatial trend of temperature in the building elements. DETECt can be also linked to other additional computer subroutines for the performance simulation of building plant systems also supported by innovative and renewable energy technologies. Thanks to the code, any innovative technology can be implemented within the simulation code, depending on the occurring research needs. Finally, with the help of all the results obtained through the presented code, useful guidelines for the energy design of new buildings (such as NZEBs) or for the energy refurbishment of the existing ones can be obtained.