Multi-Objective Optimisation of the Energy Performance of Lightweight Constructions Combining Evolutionary Algorithms and Life Cycle Cost

This paper discusses the thermal and energy performance of a detached lightweight building. The building was monitored with hygrothermal sensors to collect data for building energy model calibration. The calibration was performed using a dynamic simulation through EnergyPlus® (EP) (Version 8.5, United States Department of Energy (DOE), Washington, DC, USA) with a hybrid evolutionary algorithm to minimise the root mean square error of the differences between the predicted and real recorded data. The results attained reveal a good agreement between predicted and real data with a goodness of fit below the limits imposed by the guidelines. Then, the evolutionary algorithm was used to meet the compliance criteria defined by the Passive House standard for different regions in Portugal’s mainland using different approaches in the overheating evaluation. The multi-objective optimisation was developed to study the interaction between annual heating demand and overheating rate objectives to assess their trade-offs, tracing the Pareto front solution for different climate regions throughout the whole of Portugal. However, the overheating issue is present, and numerous best solutions from multi-objective optimisation were determined, hindering the selection of a single best option. Hence, the life cycle cost of the Pareto solutions was determined, using the life cycle cost as the final criterion to single out the optimal solution or a combination of parameters.


Introduction
The building sector represents 40% of the total European Energy balance [1]. In Portugal, the residential sector is responsible for approximately 30% of the total primary energy consumption [2]. To comply with the first Energy Performance Building Directive (EPBD) [1] requirements, new national codes were published in 2006. Following the EPBD recast publication in 2010 [1], these documents were reviewed and republished [3,4] with the goal of renovating buildings incorporating the Near Zero-Energy Building (NZEB) requirements.
The EPBD recast [1], had the ambitious goal that all new buildings should have nominal energy demands near to zero by 2020. Among the major strategies pointed out to reduce the energy consumption and greenhouse gas (GHG) emissions in the building sector [5], two key strategies Figure 1. Scheme of covariance matrix adaptation evolution strategies (CMA-ES) and hybrid differential evolution (HDE) algorithms (adapted from Reference [19]).
In the decision-making process, the economic dimension plays a key role. One of the most well accepted methodologies to evaluate and compare rehabilitation alternatives is life cycle cost (LCC) analysis. LCC is the sum of the present value of investment and operating costs for the building and service systems, including those related to maintenance and replacement, over a specified life span. For instance, Hasan et al. [23] have used LCC, combined with simulation, on the optimisation of the U-values of typical Finnish constructions. In the context of PH design, several authors applied LCC as an optimisation criteria [24][25][26]. The difficulties concerning the selection of a single optimum solution after a multi-objective optimisation were discussed in Reference [27]. The authors used LCC as a final criterion in the decision-making process. This paper is based on a lightweight building case study erected in the scope of an industrial research project "MODIKO Passive House"-in which a steel frame system was developed and optimised in order to comply with the requirements defined by the PH standard. The project included the construction of a real model in the Aveiro region (central coast of Portugal). The objective of this paper was to present and discuss the results of the optimisation procedure applied for the selection of the constructive solutions, including the adaptability of the PH standard for other climate regions of the Portuguese mainland. The combination between lightweight construction and PH design was tested, using thermal comfort, energy efficiency and LCC as performance criteria.

Methodology
The aim of this research was to present a methodology for the optimisation of cost effective envelope constructive solutions and heating ventilating and air conditioning (HVAC) systems to be applied to lightweight buildings in order to comply with PH comfort and energy requirements. To achieve this goal, a building energy model (BEM) was created using EnergyPlus (EP) and calibrated with monitoring data collected in situ.  Figure 1. Scheme of covariance matrix adaptation evolution strategies (CMA-ES) and hybrid differential evolution (HDE) algorithms (adapted from Reference [19]).
In the decision-making process, the economic dimension plays a key role. One of the most well accepted methodologies to evaluate and compare rehabilitation alternatives is life cycle cost (LCC) analysis. LCC is the sum of the present value of investment and operating costs for the building and service systems, including those related to maintenance and replacement, over a specified life span. For instance, Hasan et al. [23] have used LCC, combined with simulation, on the optimisation of the U-values of typical Finnish constructions. In the context of PH design, several authors applied LCC as an optimisation criteria [24][25][26]. The difficulties concerning the selection of a single optimum solution after a multi-objective optimisation were discussed in Reference [27]. The authors used LCC as a final criterion in the decision-making process. This paper is based on a lightweight building case study erected in the scope of an industrial research project "MODIKO Passive House"-in which a steel frame system was developed and optimised in order to comply with the requirements defined by the PH standard. The project included the construction of a real model in the Aveiro region (central coast of Portugal). The objective of this paper was to present and discuss the results of the optimisation procedure applied for the selection of the constructive solutions, including the adaptability of the PH standard for other climate regions of the Portuguese mainland. The combination between lightweight construction and PH design was tested, using thermal comfort, energy efficiency and LCC as performance criteria.

Methodology
The aim of this research was to present a methodology for the optimisation of cost effective envelope constructive solutions and heating ventilating and air conditioning (HVAC) systems to be applied to lightweight buildings in order to comply with PH comfort and energy requirements. To achieve this goal, a building energy model (BEM) was created using EnergyPlus (EP) and calibrated with monitoring data collected in situ.
The first task of the methodology was collecting temperature data from the monitoring campaign and the parameters of the real construction building to define the BEM. Secondly, the collected data were used to calibrate and validate the BEM using a hybrid evolutionary algorithm to minimise the deviation of the Root Mean Squared Error (RMSE) between the measured and simulated indoor air temperature, suiting uncertain input parameters (design variables). The accuracy of the BEM was assessed by the goodness of fit (GOF) and Coefficient of Variation Root Mean Squared Error (CV RMSE) hourly criteria, according to the methodology proposed in [21]. The hybrid evolutionary algorithm was used to instruct the engine calculation software (EP). The calibration procedure was only carried out for Aveiro as no real data were available for the other regions. The results found in the calibration may have been different if data from a different city was used. Despite this limitation, the authors believe that using real data increases the accuracy of the model and adds value to the research.
Thirdly, the optimisation procedure was applied in six regions of Portugal: Bragança, Oporto, Aveiro, Lisbon, Évora and Faro. These locations were selected because they represent the different climatic regions of the country. The performance functions included the assessment of the indoor thermal comfort while minimising heating demand and overheating rate. The optimisation was carried out using the hybrid evolutionary algorithm for a multi-objective optimisation. The opaque constructive solutions, the windows, the air infiltration rate and the HVAC system were chosen as variables for the optimisation. The multi-objective optimisation procedure provided the Pareto front of optimum solutions. In Section 5, the multi-objective optimisation using evolutionary algorithm is presented in detail. For each optimisation, a total of 10,000 simulations were carried out. The optimisation was performed using a server with Intel Core i7 5820K with six dedicated cores working at 3.3 GHz.
Finally, an economic analysis based on the LCC was implemented for the selection of the optimum solutions based on the multi-objective optimisation. This methodology is depicted in Figure 2. The first task of the methodology was collecting temperature data from the monitoring campaign and the parameters of the real construction building to define the BEM.
Secondly, the collected data were used to calibrate and validate the BEM using a hybrid evolutionary algorithm to minimise the deviation of the Root Mean Squared Error (RMSE) between the measured and simulated indoor air temperature, suiting uncertain input parameters (design variables). The accuracy of the BEM was assessed by the goodness of fit (GOF) and Coefficient of Variation Root Mean Squared Error (CV RMSE) hourly criteria, according to the methodology proposed in [21]. The hybrid evolutionary algorithm was used to instruct the engine calculation software (EP). The calibration procedure was only carried out for Aveiro as no real data were available for the other regions. The results found in the calibration may have been different if data from a different city was used. Despite this limitation, the authors believe that using real data increases the accuracy of the model and adds value to the research.
Thirdly, the optimisation procedure was applied in six regions of Portugal: Bragança, Oporto, Aveiro, Lisbon, Évora and Faro. These locations were selected because they represent the different climatic regions of the country. The performance functions included the assessment of the indoor thermal comfort while minimising heating demand and overheating rate. The optimisation was carried out using the hybrid evolutionary algorithm for a multi-objective optimisation. The opaque constructive solutions, the windows, the air infiltration rate and the HVAC system were chosen as variables for the optimisation. The multi-objective optimisation procedure provided the Pareto front of optimum solutions. In Section 5, the multi-objective optimisation using evolutionary algorithm is presented in detail. For each optimisation, a total of 10,000 simulations were carried out. The optimisation was performed using a server with Intel Core i7 5820K with six dedicated cores working at 3.3 GHz.
Finally, an economic analysis based on the LCC was implemented for the selection of the optimum solutions based on the multi-objective optimisation. This methodology is depicted in Figure 2.

Building Location and General Characterisation
The building used as case study was built in 2015 in a suburban area of the city of Aveiro on the north coast of Portugal mainland (latitude 40 • 60 , longitude 8 • 60 ) about 15 km away from the Atlantic coast and 50 m above sea level. The building had a contemporary architecture design inserted in a residential area with one family dwelling buildings of similar geometry, as displayed in Figure 3.

Building Location and General Characterisation
The building used as case study was built in 2015 in a suburban area of the city of Aveiro on the north coast of Portugal mainland (latitude 40°60′, longitude 8°60′) about 15 km away from the Atlantic coast and 50 m above sea level. The building had a contemporary architecture design inserted in a residential area with one family dwelling buildings of similar geometry, as displayed in Figure 3. The one-family building had two floors and a steel structure, shown in Figure 4. The building had 148 m 2 of treated floor area and 318 m 2 of exterior surface area, of which 198.4 m 2 were opaque external wall. The glazing to wall ratio was 15.4%, approximately. This ratio was quite different for the four external walls, with a ratio of 26% in the north and 60% in the south external wall. To prevent overheating, a shading overhang and external venetian shutters were installed in the south-oriented windows as shading devices.  The building structure was founded on a massive ground floor concrete slab with thermal insulation underneath which supported the cold formed steel structure. The exterior wall panels consisted of a thick thermal insulation layer assembled with a secondary steel structure. The flat roof was composed of thermal insulation panels fixed to zinc sheeting over the main steel structure. The constructive solutions are displayed in Table 1, including the corresponding U-values. The one-family building had two floors and a steel structure, shown in Figure 4. The building had 148 m 2 of treated floor area and 318 m 2 of exterior surface area, of which 198.4 m 2 were opaque external wall. The glazing to wall ratio was 15.4%, approximately. This ratio was quite different for the four external walls, with a ratio of 26% in the north and 60% in the south external wall. To prevent overheating, a shading overhang and external venetian shutters were installed in the south-oriented windows as shading devices.

Building Location and General Characterisation
The building used as case study was built in 2015 in a suburban area of the city of Aveiro on the north coast of Portugal mainland (latitude 40°60′, longitude 8°60′) about 15 km away from the Atlantic coast and 50 m above sea level. The building had a contemporary architecture design inserted in a residential area with one family dwelling buildings of similar geometry, as displayed in Figure 3.  The one-family building had two floors and a steel structure, shown in Figure 4. The building had 148 m 2 of treated floor area and 318 m 2 of exterior surface area, of which 198.4 m 2 were opaque external wall. The glazing to wall ratio was 15.4%, approximately. This ratio was quite different for the four external walls, with a ratio of 26% in the north and 60% in the south external wall. To prevent overheating, a shading overhang and external venetian shutters were installed in the south-oriented windows as shading devices.  The building structure was founded on a massive ground floor concrete slab with thermal insulation underneath which supported the cold formed steel structure. The exterior wall panels consisted of a thick thermal insulation layer assembled with a secondary steel structure. The flat roof was composed of thermal insulation panels fixed to zinc sheeting over the main steel structure. The constructive solutions are displayed in Table 1, including the corresponding U-values. The building structure was founded on a massive ground floor concrete slab with thermal insulation underneath which supported the cold formed steel structure. The exterior wall panels consisted of a thick thermal insulation layer assembled with a secondary steel structure. The flat roof was composed of thermal insulation panels fixed to zinc sheeting over the main steel structure. The constructive solutions are displayed in Table 1, including the corresponding U-values.  The linear thermal bridges were analysed in accordance with EN ISO 10211:2007 [28] and all comply with the Ψvalue ≤ 0.01 W/(m 2 ·°C) required by PH.
The thermal properties of the windows and doors were: Thermal transmission coefficient Uw,installed = 1.72 W/(m 2 ·°C) and Uw,installed = 1.69 W/(m 2 ·°C) for the north and south windows, respectively; solar heat gain coefficient (SHGC) = 0.57; and external doors Uw,installed = 1.90 W/(m 2 ·°C). These values take into to account the frame (Uf) and glass edge thermal bridge (Ψg) in accordance with ISO 10077 [29] and the thermal bridge due to installation (ΨInstall) in accordance with EN ISO 10211 [28]. For the HVAC system, a compact heat pump unit was chosen, incorporating: Ventilation and passive recovery unit, as well as a supplementary energy efficient heating system, to offset the air temperature when the heat exchanger cross flow is not sufficient.

Weather Files Considered for BEM
The World Map of Köppen-Geiger Climate Classification [30] establishes an organisation of the world map classified in five main climates. Portugal is assembled in the Csa (south of the country) and Csb (centre and north of the country) regions. In this study, six locations were selected, representative of different winter and summer conditions and severity, which together characterise the climatic regions of the country, shown in Figure 5. Bragança and Évora are representative of the north and south interior, respectively. Oporto, Lisbon and Faro represent the north, centre and south coast of Portugal's mainland, respectively. Aveiro is the location of the case study.  The linear thermal bridges were analysed in accordance with EN ISO 10211:2007 [28] and all comply with the Ψvalue ≤ 0.01 W/(m 2 ·°C) required by PH.
The thermal properties of the windows and doors were: Thermal transmission coefficient Uw,installed = 1.72 W/(m 2 ·°C) and Uw,installed = 1.69 W/(m 2 ·°C) for the north and south windows, respectively; solar heat gain coefficient (SHGC) = 0.57; and external doors Uw,installed = 1.90 W/(m 2 ·°C). These values take into to account the frame (Uf) and glass edge thermal bridge (Ψg) in accordance with ISO 10077 [29] and the thermal bridge due to installation (ΨInstall) in accordance with EN ISO 10211 [28]. For the HVAC system, a compact heat pump unit was chosen, incorporating: Ventilation and passive recovery unit, as well as a supplementary energy efficient heating system, to offset the air temperature when the heat exchanger cross flow is not sufficient.

Weather Files Considered for BEM
The World Map of Köppen-Geiger Climate Classification [30] establishes an organisation of the world map classified in five main climates. Portugal is assembled in the Csa (south of the country) and Csb (centre and north of the country) regions. In this study, six locations were selected, representative of different winter and summer conditions and severity, which together characterise the climatic regions of the country, shown in Figure 5. Bragança and Évora are representative of the north and south interior, respectively. Oporto, Lisbon and Faro represent the north, centre and south coast of Portugal's mainland, respectively. Aveiro is the location of the case study.  The linear thermal bridges were analysed in accordance with EN ISO 10211:2007 [28] and all comply with the Ψvalue ≤ 0.01 W/(m 2 ·°C) required by PH.
The thermal properties of the windows and doors were: Thermal transmission coefficient Uw,installed = 1.72 W/(m 2 ·°C) and Uw,installed = 1.69 W/(m 2 ·°C) for the north and south windows, respectively; solar heat gain coefficient (SHGC) = 0.57; and external doors Uw,installed = 1.90 W/(m 2 ·°C). These values take into to account the frame (Uf) and glass edge thermal bridge (Ψg) in accordance with ISO 10077 [29] and the thermal bridge due to installation (ΨInstall) in accordance with EN ISO 10211 [28]. For the HVAC system, a compact heat pump unit was chosen, incorporating: Ventilation and passive recovery unit, as well as a supplementary energy efficient heating system, to offset the air temperature when the heat exchanger cross flow is not sufficient.

Weather Files Considered for BEM
The World Map of Köppen-Geiger Climate Classification [30] establishes an organisation of the world map classified in five main climates. Portugal is assembled in the Csa (south of the country) and Csb (centre and north of the country) regions. In this study, six locations were selected, representative of different winter and summer conditions and severity, which together characterise the climatic regions of the country, shown in Figure 5. Bragança and Évora are representative of the north and south interior, respectively. Oporto, Lisbon and Faro represent the north, centre and south coast of Portugal's mainland, respectively. Aveiro is the location of the case study.
The linear thermal bridges were analysed in accordance with EN ISO 10211:2007 [28] and all comply with the Ψ value ≤ 0.01 W/(m 2 · • C) required by PH.
The thermal properties of the windows and doors were: Thermal transmission coefficient U w,installed = 1.72 W/(m 2 · • C) and U w,installed = 1.69 W/(m 2 · • C) for the north and south windows, respectively; solar heat gain coefficient (SHGC) = 0.57; and external doors U w,installed = 1.90 W/(m 2 · • C). These values take into to account the frame (U f ) and glass edge thermal bridge (Ψ g ) in accordance with ISO 10077 [29] and the thermal bridge due to installation (Ψ Install ) in accordance with EN ISO 10211 [28]. For the HVAC system, a compact heat pump unit was chosen, incorporating: Ventilation and passive recovery unit, as well as a supplementary energy efficient heating system, to offset the air temperature when the heat exchanger cross flow is not sufficient.

Weather Files Considered for BEM
The World Map of Köppen-Geiger Climate Classification [30] establishes an organisation of the world map classified in five main climates. Portugal is assembled in the Csa (south of the country) and Csb (centre and north of the country) regions. In this study, six locations were selected, representative of different winter and summer conditions and severity, which together characterise the climatic regions of the country, shown in Figure 5. Bragança and Évora are representative of the north and south interior, respectively. Oporto, Lisbon and Faro represent the north, centre and south coast of Portugal's mainland, respectively. Aveiro is the location of the case study.

Ground Floor Slab External Walls Flat Roof
The linear thermal bridges were analysed in accordance with EN ISO 10211:2007 [28] and all comply with the Ψvalue ≤ 0.01 W/(m 2 ·°C) required by PH.
The thermal properties of the windows and doors were: Thermal transmission coefficient Uw,installed = 1.72 W/(m 2 ·°C) and Uw,installed = 1.69 W/(m 2 ·°C) for the north and south windows, respectively; solar heat gain coefficient (SHGC) = 0.57; and external doors Uw,installed = 1.90 W/(m 2 ·°C). These values take into to account the frame (Uf) and glass edge thermal bridge (Ψg) in accordance with ISO 10077 [29] and the thermal bridge due to installation (ΨInstall) in accordance with EN ISO 10211 [28]. For the HVAC system, a compact heat pump unit was chosen, incorporating: Ventilation and passive recovery unit, as well as a supplementary energy efficient heating system, to offset the air temperature when the heat exchanger cross flow is not sufficient.

Weather Files Considered for BEM
The World Map of Köppen-Geiger Climate Classification [30] establishes an organisation of the world map classified in five main climates. Portugal is assembled in the Csa (south of the country) and Csb (centre and north of the country) regions. In this study, six locations were selected, representative of different winter and summer conditions and severity, which together characterise the climatic regions of the country, shown in Figure 5. Bragança and Évora are representative of the north and south interior, respectively. Oporto, Lisbon and Faro represent the north, centre and south coast of Portugal's mainland, respectively. Aveiro is the location of the case study.  The weather files for these regions were obtained from the National Laboratory for Energy and Geology (LNEG) [32]. The files were built from weather data collected by the Portuguese Sea and Atmosphere Institute (IPMA) using measurements carried out between 1971 and 2000. The weather data of the regions used in the simulations are summarised in Figures 6 and 7. The weather files for these regions were obtained from the National Laboratory for Energy and Geology (LNEG) [32]. The files were built from weather data collected by the Portuguese Sea and Atmosphere Institute (IPMA) using measurements carried out between 1971 and 2000. The weather data of the regions used in the simulations are summarised in Figure 6 and Figure 7. Figure 6 and Figure 7 show that for all regions the months with both the highest temperatures and highest irradiances occur from June to August with the lowest relative humidity. Évora and Faro represent the hottest regions of the country. The lowest temperature and irradiance occur from December to February, with the highest relative humidity. Bragança, Oporto and Aveiro, located in north of the country, have more severe winter conditions. The locations with higher solar irradiance are Lisbon and Faro, both on the coastline.

Monitoring of the Building
The real data collected during the hygrothermal monitoring campaign were fundamental to achieving an accurate BEM. To that the end, thermo-hygrometer sensors were installed in six compartments to record the air temperature and relative humidity. The temperature probe had an accuracy of 0.5 °C and a resolution of 0.1 °C while the humidity probe had an accuracy of 3% and a  The weather files for these regions were obtained from the National Laboratory for Energy and Geology (LNEG) [32]. The files were built from weather data collected by the Portuguese Sea and Atmosphere Institute (IPMA) using measurements carried out between 1971 and 2000. The weather data of the regions used in the simulations are summarised in Figure 6 and Figure 7. Figure 6 and Figure 7 show that for all regions the months with both the highest temperatures and highest irradiances occur from June to August with the lowest relative humidity. Évora and Faro represent the hottest regions of the country. The lowest temperature and irradiance occur from December to February, with the highest relative humidity. Bragança, Oporto and Aveiro, located in north of the country, have more severe winter conditions. The locations with higher solar irradiance are Lisbon and Faro, both on the coastline.

Monitoring of the Building
The real data collected during the hygrothermal monitoring campaign were fundamental to achieving an accurate BEM. To that the end, thermo-hygrometer sensors were installed in six compartments to record the air temperature and relative humidity. The temperature probe had an accuracy of 0.5 °C and a resolution of 0.1 °C while the humidity probe had an accuracy of 3% and a  Figures 6 and 7 show that for all regions the months with both the highest temperatures and highest irradiances occur from June to August with the lowest relative humidity. Évora and Faro represent the hottest regions of the country. The lowest temperature and irradiance occur from December to February, with the highest relative humidity. Bragança, Oporto and Aveiro, located in north of the country, have more severe winter conditions. The locations with higher solar irradiance are Lisbon and Faro, both on the coastline.

Monitoring of the Building
The real data collected during the hygrothermal monitoring campaign were fundamental to achieving an accurate BEM. To that the end, thermo-hygrometer sensors were installed in six compartments to record the air temperature and relative humidity. The temperature probe had an accuracy of 0.5 • C and a resolution of 0.1 • C while the humidity probe had an accuracy of 3% and a resolution of 0.1%. The building was monitored from the 22 October 2017 to 21 December 2017. The sensor's positions inside the compartments were chosen in order to avoid direct sun exposure from the glazed areas. The monitoring acquisition system was logged at 10-minute intervals and averaged hourly. The sensors were distributed in accordance with ISO 7726 recommendations [33], positioned approximately at a 1-m height. Figure 8 shows the position of the thermo-hygrometer (TH) sensors in the compartments. resolution of 0.1%. The building was monitored from the 22 October 2017 to 21 December 2017. The sensor's positions inside the compartments were chosen in order to avoid direct sun exposure from the glazed areas. The monitoring acquisition system was logged at 10-minute intervals and averaged hourly. The sensors were distributed in accordance with ISO 7726 recommendations [33], positioned approximately at a 1-m height. Figure 8 shows the position of the thermo-hygrometer (TH) sensors in the compartments.
(a) (b) The exterior dry bulb temperature, relative humidity, solar irradiation, and wind speed and direction were collected from a nearby weather station and used in the BEM for calibration purposes. Direct and diffuse irradiance were converted from the monitored global horizontal solar radiance using the operational model developed by Cipriano et al. [34]. Figure 9 presents the variation in the air temperature in the six compartments and in the exterior during the monitoring period. These results show that the indoor temperature practically complies with the thermal comfort range (20-25 °C) established by PH requirements during almost the entire period. The air infiltration was evaluated by a blower door test, shown in Figure 10a. The air change rate at a pressure difference of 50 Pa (n50) was 0.9 h −1 which indicates a very airtight building envelope and represents 10% below the limit value of 1.0 h −1 required by the PH standard for European warm climates, as presented in the Passive-On [35] study.
The air change rate (ACR) due to the Mechanical Ventilation with Heat Recovery (MVHR) system was evaluated with discrete measurements of the air velocity at the ventilation grids of each compartment, shown in Figure 10b  The exterior dry bulb temperature, relative humidity, solar irradiation, and wind speed and direction were collected from a nearby weather station and used in the BEM for calibration purposes. Direct and diffuse irradiance were converted from the monitored global horizontal solar radiance using the operational model developed by Cipriano et al. [34]. Figure 9 presents the variation in the air temperature in the six compartments and in the exterior during the monitoring period. These results show that the indoor temperature practically complies with the thermal comfort range (20-25 • C) established by PH requirements during almost the entire period.

201
The exterior dry bulb temperature, relative humidity, solar irradiation, and wind speed and 202 direction were collected from a nearby weather station and used in the BEM for calibration purposes.

203
Direct and diffuse irradiance were converted from the monitored global horizontal solar radiance 204 using the operational model developed by Cipriano et al. [34].   The air infiltration was evaluated by a blower door test, shown in Figure 10a. The air change rate at a pressure difference of 50 Pa (n 50 ) was 0.9 h −1 which indicates a very airtight building envelope and represents 10% below the limit value of 1.0 h −1 required by the PH standard for European warm climates, as presented in the Passive-On [35] study.
The air change rate (ACR) due to the Mechanical Ventilation with Heat Recovery (MVHR) system was evaluated with discrete measurements of the air velocity at the ventilation grids of each compartment, shown in Figure 10b. The air velocity was logged with an anemometer with an accuracy of 0.2 m/s and a measuring range of 0.6 to 40 m/s. In each compartment, the final velocity was the average of nine points around the ventilation grids in accordance with Reference [36], as illustrated in Figure 10c. was the average of nine points around the ventilation grids in accordance with Reference [36], as illustrated in Figure 10c.

Numerical Model Definition
Based on the parameters defined in the previous sections (3.1-3.3), a dynamic thermal simulation model was created in OpenStudio ® (Version 2.6.0, National Renewable Energy Laboratory, Golden, CO, USA) to be simulated in EnergyPlus ® software. Figure 11 presents the Northeast and Southwest views of the BEM designed through the SketchUp ® and OpenStudio ® plugin where the adjacent constraints were included in the BEM geometry as shading surfaces. A complete description of the BEM can be found in [20]. The model includes nine thermal zones (TZ), corresponding to the internal compartments of the building, shown in Figure 11. Thermal zone TZ01-Living room, which includes the hall and the staircase, establishes the connection with the elevated floor and thus is considered in both. The remaining TZ are as follows: TZ02-kitchen; TZ03-ground floor suite bedroom; TZ04-ground floor bathroom; TZ05-first floor suite bedroom; TZ06-first floor bathroom and TZ07/TZ08-first floor single room. Real schedules were defined for the occupancy, equipment, electric devices and artificial lighting hours. The building was occupied by a family composed by two adults and two children (100% of occupation level/density) and the schedule routine was based on the behaviour of the residents, as presented in Table 2. During the remaining hours the zones were unoccupied. The remaining internal gains (lighting and electric equipment) were included in the model, also in accordance with real schedules provided by users.

Numerical Model Definition
Based on the parameters defined in the previous sections (3.1-3.3), a dynamic thermal simulation model was created in OpenStudio ® (Version 2.6.0, National Renewable Energy Laboratory, Golden, CO, USA) to be simulated in EnergyPlus ® software. Figure 11 presents the Northeast and Southwest views of the BEM designed through the SketchUp ® and OpenStudio ® plug-in where the adjacent constraints were included in the BEM geometry as shading surfaces. A complete description of the BEM can be found in [20]. The model includes nine thermal zones (TZ), corresponding to the internal compartments of the building, shown in Figure 11. Thermal zone TZ01-Living room, which includes the hall and the staircase, establishes the connection with the elevated floor and thus is considered in both. The remaining TZ are as follows: TZ02-kitchen; TZ03-ground floor suite bedroom; TZ04-ground floor bathroom; TZ05-first floor suite bedroom; TZ06-first floor bathroom and TZ07/TZ08-first floor single room.
Energies 2018, 11, x FOR PEER REVIEW 9 of 23 was the average of nine points around the ventilation grids in accordance with Reference [36], as illustrated in Figure 10c.

Numerical Model Definition
Based on the parameters defined in the previous sections (3.1-3.3), a dynamic thermal simulation model was created in OpenStudio ® (Version 2.6.0, National Renewable Energy Laboratory, Golden, CO, USA) to be simulated in EnergyPlus ® software. Figure 11 presents the Northeast and Southwest views of the BEM designed through the SketchUp ® and OpenStudio ® plugin where the adjacent constraints were included in the BEM geometry as shading surfaces. A complete description of the BEM can be found in [20]. The model includes nine thermal zones (TZ), corresponding to the internal compartments of the building, shown in Figure 11. Thermal zone TZ01-Living room, which includes the hall and the staircase, establishes the connection with the elevated floor and thus is considered in both. The remaining TZ are as follows: TZ02-kitchen; TZ03-ground floor suite bedroom; TZ04-ground floor bathroom; TZ05-first floor suite bedroom; TZ06-first floor bathroom and TZ07/TZ08-first floor single room. Real schedules were defined for the occupancy, equipment, electric devices and artificial lighting hours. The building was occupied by a family composed by two adults and two children (100% of occupation level/density) and the schedule routine was based on the behaviour of the residents, as presented in Table 2. During the remaining hours the zones were unoccupied. The remaining internal gains (lighting and electric equipment) were included in the model, also in accordance with real schedules provided by users. Real schedules were defined for the occupancy, equipment, electric devices and artificial lighting hours. The building was occupied by a family composed by two adults and two children (100% of occupation level/density) and the schedule routine was based on the behaviour of the residents, as presented in Table 2. During the remaining hours the zones were unoccupied. The remaining internal gains (lighting and electric equipment) were included in the model, also in accordance with real schedules provided by users. The MVHR was always in operation to guarantee adequate indoor air quality and thermal comfort. During the monitored period, the external venetian shading devices were manually controlled by the occupants. Once again, the shading device's schedule was defined according to the occupants' behaviour: (i) closed-from 18:00 to 8:00 and (ii) open-from 8:00 to 18:00.

Definition of the Unknown Parameters
In the first step of the model calibration, 15 input parameters (x0 to x14) were identified as the ones with the highest uncertainty. Consequently, a range of variation was defined for each one, presented in Table 3. Parameters x0 to x3 (insulation thickness) were selected as uncertainty variables to take into to account the thermal bridge's impact and, therefore, their range of variation is rather narrow. Parameters x4 and x5 are related to the exterior windows U value coefficient, with a range limit (5%) of the value given by window manufacturer. Parameter x6 (n 50 ) has a significant impact on the indoor air temperatures, as well as on the space heating or cooling energy needs in buildings. The range defined (0.3 to 1.5 h −1 ) was based on a blower door test carried out with a 50 Pa pressure difference. Parameters from x7 to x14 represent the air change rate (m 3 /h) provided by the mechanical ventilation system. The ACR was directly evaluated from the ventilation grids using different parameters divided by thermal zone. All uncertainties were considered as continuous variables with the box constraint limits defined in Table 3.

Accuracy Criteria for the Building Energy Model
The model was validated according to the criteria (or the limit values) defined by the following standards: (i) American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) guideline 14 [37], (ii) the International Performance Measurement and Verification Protocol (IPMVP) [38] and (iii) the Federal Energy Management Program (FEMP) [39]. The strategy used to calibrate the model consists of achieving the best match between measured and simulated indoor air temperatures. For the BEM calibration, two indexes were evaluated to assess the model's accuracy: (1) normalised mean bias error (NMBE) and (2) coefficient-of-variation of the root mean square error (CV RMSE). The assumption of standardised statistical indexes is proposed by the authors in Reference [40] to represent the performance of a model defined in References [40][41][42].
The NMBE index establishes the fit of the simulated and monitored data. NMBE is defined by Equation (1) and calculated as the sum of errors between the measured and simulated data for each timestep (hourly data) of the considered period.
The calibration process, according to the methodology proposed, defines the root mean square (RMSE) coefficient as the objective function in the evolutionary algorithm, which establishes the basis to CV RMSE evaluation. The RMSE is a measure of the variability of the data defined by Equation (2). The difference in paired data points is calculated and squared hourly.
The CV RMSE index enables the determination of how well a model fits the data by capturing offsetting errors between measured and simulated data and acting as a measure of the GOF index. The CV RMSE index is calculated using Equation (3): In this equation, M i is the measured data and S i is the simulated data for the time i, which refers to the hour of the day. The M i is the average value of the measured data points for the considered period n (hourly).
The GOF and CV RMSE indexes were the selected criteria to validate the model's accuracy. A consolidated way to define an index based on NMBE and CVRMSE is to calculate GOF. The GOF index is determined by Equation (4):

Evaluation of the Calibration Procedure
The EnergyPlus outputs were programmed using the Energy Management System (EMS) application to calculate the RMSE index. Then, the evolutionary algorithm was used to minimise the objective function (RMSE index) aimed at finding a trade-off between the input design parameters (defined as variables in Table 3) according to monitored indoor air temperatures. Figure 12a shows the coefficient of determination (r 2 ) between the real data and simulated results derived from the best id. attained. The optimal solution is shown in Figure 12b. For temperatures between 20 • C and 22 • C, the deviation between simulated results and monitored data reveals a tendency of higher temperatures for the monitored data. From 22 • C to 24 • C, the observed deviation can be considered well fit, with small discrepancies in this range. In an overall analysis, the coefficient of determination of 0.99 reveals the accuracy of model-fitting.  The CV RMSE values attained, shown in Table 4, are below the limit values imposed in the ASHRAE Guidelines [37], the IPMVP [38] and the FEMP [39] standards. The validation of the BEM follows the criteria of these standards, as shown in Table 4 According to the GOF index, ASHRAE Guidelines recommend a GOF below 11% and Cipriano et al. [34] suggests a GOF below to 3% for trial's acceptance. In the presented work, the attained GOF is lower than the values suggested. In this study the three guidelines identified were verified, meaning that there is an acceptable agreement between measured and simulated data.

Multi-Objective Optimisation Using Evolutionary Algorithm
Finding the best design solution for a PH building is a complex problem in the sense that there are several good solutions resulting from the combination of parameters which satisfy the energy demand limits and comfort requirements. Thus, the use of evolutionary algorithms with multiple conflicting objectives could be a possible tool to use for assessing design actions.
The optimisation process starts from the calibrated model described in the previous section, followed by the multi-objective optimisation process using the hybrid evolutionary algorithm to develop the optimised BEMs.
The objective functions are the goals of the optimisation process and are discussed in Section 5.2. The calculation of the objective functions was performed using the Energy Management System (EMS) application of EnergyPlus ® . The hybrid evolutionary algorithm was used to instruct the energy simulation program and to select a set of optimised alternatives. The CV RMSE values attained, shown in Table 4, are below the limit values imposed in the ASHRAE Guidelines [37], the IPMVP [38] and the FEMP [39] standards. The validation of the BEM follows the criteria of these standards, as shown in Table 4 According to the GOF index, ASHRAE Guidelines recommend a GOF below 11% and Cipriano et al. [34] suggests a GOF below to 3% for trial's acceptance. In the presented work, the attained GOF is lower than the values suggested. In this study the three guidelines identified were verified, meaning that there is an acceptable agreement between measured and simulated data.

Multi-Objective Optimisation Using Evolutionary Algorithm
Finding the best design solution for a PH building is a complex problem in the sense that there are several good solutions resulting from the combination of parameters which satisfy the energy demand limits and comfort requirements. Thus, the use of evolutionary algorithms with multiple conflicting objectives could be a possible tool to use for assessing design actions.
The optimisation process starts from the calibrated model described in the previous section, followed by the multi-objective optimisation process using the hybrid evolutionary algorithm to develop the optimised BEMs.
The objective functions are the goals of the optimisation process and are discussed in Section 5.2. The calculation of the objective functions was performed using the Energy Management System (EMS) application of EnergyPlus ® . The hybrid evolutionary algorithm was used to instruct the energy simulation program and to select a set of optimised alternatives. A multi-objective optimisation process, to optimise buildings parameters, involves identifying the optimisation parameters: envelope solutions and control systems features presented and discussed in Section 5.3.
After optimisation parameter analysis and the objectives function calculation process mentioned above, the detailed implementation of the multi-objective optimisation process using an evolutionary algorithm for solving BEM problem was performed. Finally, results generated after the multi-objective optimisation process can be reported using the Pareto-optimal fronts indicating all non-dominated solutions. These results are presented in Section 5.4.
For a thorough understanding of the whole process of the multi-objective optimisation, a workflow is summarised in Figure 13. A multi-objective optimisation process, to optimise buildings parameters, involves identifying the optimisation parameters: envelope solutions and control systems features presented and discussed in Section 5.3.
After optimisation parameter analysis and the objectives function calculation process mentioned above, the detailed implementation of the multi-objective optimisation process using an evolutionary algorithm for solving BEM problem was performed. Finally, results generated after the multiobjective optimisation process can be reported using the Pareto-optimal fronts indicating all nondominated solutions. These results are presented in Section 5.4.
For a thorough understanding of the whole process of the multi-objective optimisation, a workflow is summarised in Figure 13.

Objective Functions
The overheating rate (with a maximum limit of 10%) and heating demand (with a maximum limit of 15 kW·h/(m 2 ·a)) were settled as the objective functions (PH compliance criteria).
The overheating rate is the sum of the discomfort time for the entire building according to the PH requirements and to EN 15251 [43] for category II relative to a normal comfort level adjusted for new buildings. The 25 °C limit is defined by the PH standard, however a more permissible limit of 26 °C was also evaluated taking into account the warm climate context. The overheating rate upper limit established by the adaptive comfort analysis of the EN 15251 was the third criteria chosen. In a second approach for the climatic regions where the overheating rate criteria is not fulfilled, an alternative cooling criterion was tested, adding a cooling unit into the model. In these cases, the cooling demand (with a maximum limit of 15 kW·h/(m 2 ·a)) was used as an objective function.

Optimisation Parameters
In the optimisation process, the identification of the parameters in use puts forward a set of alternative measures for building design. Table 5 presents the list of parameters for the input database.
The minimum of the reference parameters (Uvalue for opaque and translucent surfaces) was defined from Portuguese standards (REH) [3]. In the BEM, a schedule was created to control the blinds activation (automatic shading device with horizontal slats regulated between 5° to 90°) in the summer season triggered by indoor air temperature values. The shading system used allows for the possibility to achieve the optimal activation temperature for each climate region studied, also including optimisation of the slat angle of the blinds. The external blinds worked in automatic mode

Objective functions Optimisation parameters Run optimisations and iterative simulations
• Overheating rate: • Heating demand:

Objective Functions
The overheating rate (with a maximum limit of 10%) and heating demand (with a maximum limit of 15 kW·h/(m 2 ·a)) were settled as the objective functions (PH compliance criteria).
The overheating rate is the sum of the discomfort time for the entire building according to the PH requirements and to EN 15251 [43] for category II relative to a normal comfort level adjusted for new buildings. The 25 • C limit is defined by the PH standard, however a more permissible limit of 26 • C was also evaluated taking into account the warm climate context. The overheating rate upper limit established by the adaptive comfort analysis of the EN 15251 was the third criteria chosen. In a second approach for the climatic regions where the overheating rate criteria is not fulfilled, an alternative cooling criterion was tested, adding a cooling unit into the model. In these cases, the cooling demand (with a maximum limit of 15 kW·h/(m 2 ·a)) was used as an objective function.

Optimisation Parameters
In the optimisation process, the identification of the parameters in use puts forward a set of alternative measures for building design. Table 5 presents the list of parameters for the input database.
The minimum of the reference parameters (U value for opaque and translucent surfaces) was defined from Portuguese standards (REH) [3]. In the BEM, a schedule was created to control the blinds activation (automatic shading device with horizontal slats regulated between 5 • to 90 • ) in the summer season triggered by indoor air temperature values. The shading system used allows for the possibility to achieve the optimal activation temperature for each climate region studied, also including optimisation of the slat angle of the blinds. The external blinds worked in automatic mode from the 1 June to 31 October. The parameters used in the optimiser are presented in Table 5, as well as the box constraints. Two different pieces of mechanical equipment for heating and ventilation were combined as a discrete variable to simulate two scenarios: (i) traditional HVAC system and (ii) heat pump unit with heat recovery ventilation (MVHR). The MVHR is an essential component to achieve the PH requirements, however, in some regions their use could be unnecessary. This conclusion was presented in [20], which defined that the heat recovery feature of the mechanical ventilation system may not be necessary in regions with low heating demand, such as the southern regions of the Portuguese mainland.
The parameters defined for each region were combined with six different windows solutions, shown in Table 6, integrated as a discrete variable in the optimisation process. These solutions resulted from a survey on the typical solutions used in southern Europe. The data for this survey were collected from manufacturers with windows expertise.

Numerical Results
The results are presented as a Pareto front, which represents the optimal set of solutions after 10,000 runs for each region. Pareto fronts are frequently used for analysing BEM optimisation results with multi-objective solutions, defining a set of design solutions as an optimal solution set. For these possible solutions, a unique situation occurs where a single objective adversely affects other objectives [44].
The PH features optimisation, shown in Tables 5 and 6, for different climate regions in the Portuguese mainland were performed and are displayed as data plots shown in Figure 14. The weather files used for each region were presented in detail in Section 3.2. Building stakeholders can assess the trade-offs between optimised design solutions of the Pareto front for three comfort criteria defined and referred to in Section 5.3.
The multi-optimisation was computed, defining the optimum constructive solutions to fulfil the PH requirements, shown in Figure 14. The results show an exceedance of the overheating risk rate when considering the EN15251 standard approach. For the others approaches, and considering the maximum admissible value of 10% for the overheating risk, an active cooling system is needed for air conditioning, excluding the Oporto and Aveiro regions, shown in Figure 14a,d, respectively, where the indoor air temperatures are lowest.
Oporto, Aveiro, Lisbon, Évora and Faro, characterised by moderate temperatures felt during the winter season, provide an easy way to comply with the heating demand limit for the criteria under study. For Bragança, the heating demand is higher, which is characterised by a more severe winter, however, the heating demand limit was accomplished. Figure 14 presents the Pareto front results for the three comfort criteria defined in Section 5.3 associated with the overheating rate and cooling demand limit. These results, which are severely influenced by the weather location conditions, are in line with others research studies presented in References [10,45]. when considering the EN15251 standard approach. For the others approaches, and considering the maximum admissible value of 10% for the overheating risk, an active cooling system is needed for air conditioning, excluding the Oporto and Aveiro regions, shown in Figure 14a,d, respectively, where the indoor air temperatures are lowest. Oporto, Aveiro, Lisbon, Évora and Faro, characterised by moderate temperatures felt during the winter season, provide an easy way to comply with the heating demand limit for the criteria under study. For Bragança, the heating demand is higher, which is characterised by a more severe winter, however, the heating demand limit was accomplished. Figure 14 presents the Pareto front results for the three comfort criteria defined in Section 5.3 associated with the overheating rate and cooling demand limit. These results, which are severely influenced by the weather location conditions, are in line with others research studies presented in References [10,45].

Discussion
In summary, these results suggest that using the EN 15251 criteria, the PH requirement for both limits was easily achieved for all climate regions. Moreover, the majority of the scenarios in the Pareto front accomplish the heating demand.
For the overheating rate limits, considering the 25 °C and 26 °C maximum temperature criteria, the PH requirements were achieved only for the Oporto and Aveiro regions. Thus, an active cooling system was needed to reduce the overheating rate for the other climate regions under study. The use of an active system allows the compliance of the cooling demand criteria, for Lisbon, Évora and Faro, preventing the overheating risk with a lower cooling demand.
In order to adapt the PH concept for buildings in southern European warm climates, it is plausible to consider the maximum indoor air temperature limit 26 °C and therefore, this was chosen as the main criteria, using the descriptive statistics shown next. The results were statistically analysed and organised by region climate, shown in Table 7. The coefficient of variation (CV), mean and standard deviation of each set of parameters are included in the table.

Discussion
In summary, these results suggest that using the EN 15251 criteria, the PH requirement for both limits was easily achieved for all climate regions. Moreover, the majority of the scenarios in the Pareto front accomplish the heating demand.
For the overheating rate limits, considering the 25 • C and 26 • C maximum temperature criteria, the PH requirements were achieved only for the Oporto and Aveiro regions. Thus, an active cooling system was needed to reduce the overheating rate for the other climate regions under study. The use of an active system allows the compliance of the cooling demand criteria, for Lisbon, Évora and Faro, preventing the overheating risk with a lower cooling demand.
In order to adapt the PH concept for buildings in southern European warm climates, it is plausible to consider the maximum indoor air temperature limit 26 • C and therefore, this was chosen as the main criteria, using the descriptive statistics shown next. The results were statistically analysed and organised by region climate, shown in Table 7. The coefficient of variation (CV), mean and standard deviation of each set of parameters are included in the table. One of the main features of the PH concept is the thermal building envelope defined with the purposes of reducing heat losses. Regarding the thermal insulation of the opaque envelope, high thickness of thermal insulation is fundamental as an average value close to the maximum allowed was attained for parameters #01 and #02. Concerning the insulation of the ground floor slab (#00), no clear trend was identified in the results, nor its dependence on the exterior climate. In fact, the variability of the insulation of the ground floor slab is related with the higher annual thermal amplitude, since higher values were found in cities with severe winter and summer. In colder climates, the insulation of walls (#01) is close to the maximum, while in warmer climates (Lisbon, Évora and Faro), more variability was observed, with the results pointing to higher CV. The maximum thickness was obtained for the roof insulation (#02), regardless of the location. There was a clear trend for window solutions (#03 and #04), with triple glazing windows (S6) for North and South orientation without significant variability between climates. Hence, the strongly insulated opaque envelope and double to triple glazing envelope is desirable, with the benefit of minimising the indoor air temperature asymmetry [46].
Regarding the shadings, the slat angle parameter (#05) tends to be low, which can be explained by an irrelevant importance of solar gains. Although, Bragança, Lisboa, Évora and Faro had higher values, with a slat angle at 90 • positioned (closed) in some cases as they have a considerable CV. The activation temperature of blinds (#06) is relatively equal for all climate cases, observing a tendency to increase in hotter climates, with an ideal activation temperature about 22.5 • C for all climate regions.
The bypass activation temperature (#07) reveals no obvious relationship with climate severity. In addition, bypass air flow rate (#08) showed no clear relationship with climate regions, excluding the fact that there is a greater variability in climate regions with greater annual thermal amplitudes. The statistical results show no clear trend of the air change rate (#09 to #16) and a preference for the use of a ventilation system with heat recovery in relation to the common mechanical system (#17).
The statistical results reveal the need to have a strongly insulated thermal envelope and very good glazing solutions are needed to achieve the PH requirements for the Portuguese mainland. Other input parameters showed lower variability, demonstrating that these are stabilised values, leading to consensual choices for the designers. However, the input parameters allow a set of combinations of solutions (input parameters) for the input database meeting the PH compliance criteria. Although, an additional criterion, as an economic referee, is needed to simplify the options for the decision maker.

Life Cycle Costs Analysis
This section is an essential complement to the previously discussed descriptive statistical analysis of the best proposals. The main goal was to provide the building stakeholders and designers with the crucial information on which to base decision making. The multi-objective optimisation with differing objectives (two or more) results always range of optimum solutions and therefore, theoretically, it is not possible to define what solution is better than others. Introducing an additional criterion, besides the objective functions defined in the optimisation process, can be a successful strategy if one intends to achieve a single solution. Thus, by including an economic analysis, the authors attempt to identify the most cost-effective scenario from the set of Pareto optimum solutions. Traditionally, the result of an economic analysis is included in the multi-objective optimisation as an additional objective. However, the intention of this paper was to explore an alternative manner of taking into account this goal. With this approach, the relative importance given to this aspect is higher and it is possible to achieve a single optimum solution.
In this study, the investment evaluation was carried out using the LCC methodology developed for the analysis of energy conservation projects [47]. The LCC analysis was fulfilled evaluating three types of cost: investment; operation and maintenance (O&M) and energy costs. Thus, to calculate the LCC of an energy conservation solution, the following Equation (3) was used [47]: In this equation, I is the present value (PV) of investment costs, O&M is the present value of operating and maintenance, and E is the energy costs present value.
The methodology requires the calculation of the present value of the entire cash flow throughout the period of analysis. A discount rate of 3% was chosen and the LCC analysis was carried out for a lifetime of 20 years. The energy related costs were calculated a price of 0.18 €/kW·h, with a corresponding annual escalation rate (e) of 0.20%.
The main objective of LCC analysis is to define an optimum solution from the solutions/results that define the Pareto front. This method was also implemented by the authors of Reference [27], using the economic optimisation as a decision factor, to select the best constructive solution amongst the dataset of the Pareto front. For the LCC calculations, the input data of variable investment includes the thermal insulation and the window solution prices. The remaining inputs used in the optimisation process have no influence over the investment cost. The installation cost of the variables investment is presented in Table 8. The LCC analysis was carried out for all regions in the study. The LCC results are presented in Table 9. The lowest LCC values of the Pareto front solutions are shown in Table 9, highlighting the best solution (id#) and their variable's values separately for each climate regions. Faro, Évora and Lisbon are the regions with the lowest LCC values. The Aveiro region presents the highest LCC value, representing an increase of approximately 44% of the LCC value, when compared with Faro, the region with lowest LCC. The higher LCC values occur in the cities located in the north region, while the lowest LCC values occur in the south region cities with warmer climates.

Conclusions
This study focuses on a dynamic thermal simulation of a light steel frame building using the EP software as a calculation engine. An evolutionary algorithm was used to calibrate the building energy model and to find the best constructive solutions and features of the active systems to fulfil the PH requirements for different regions of the Portugal mainland. Regarding the calibration process, the differences between real data and simulated results were minimised and the attained results allow us to conclude that the building energy model is calibrated according to the standards ASHRAE guideline 14 [37], IPMVP [38] and FEMP [39].
The constructive technology of lightweight steel frame construction offers a significant cost reduction compared to the traditional massive construction solutions and allows a higher quality control of construction procedures. However, the building technology studied herein presents a low thermal inertia and increased consequential risk of overheating [48]. The Passive House concept applied to lightweight constructive systems is viable for Portugal, however some adaptations of the constructive solutions are needed. The overheating risk obliges the use of simple upgrades in the constructive components and features, such as the use of automated blinds, to attain even better results with respect to the overheating reduction.
For southwest European climates, the overheating criteria imposed by the Passive House standard was not achieved for all climate regions of Portugal mainland and the adaptive comfort model is recommended. For some regions, an active cooling system is required, however, with low cooling demand to comply with the 25 • C criterion limit.
The features optimisation allowed the definition of several models that represent good solutions that comply with the Passive House requirements. However, stakeholders and designers must select the best solution that suits economic feasibility and high thermal comfort. Several best options are attained with two conflicting objectives, because there is no single dominated solution. Thus, a life cycle cost analysis was carried out to calculate the minimum life cycle cost value as a criterion of the decision to find and single out a best solution assuring the viability of the best combination designs. The methodology proposed in this study increases the energy efficiency of building construction, revealing a two-step approach to support decision making of the best building features to assure Passive House compliance criteria.