Methodological Approach for the Development of a Simpliﬁed Residential Building Energy Estimation in Temperate Climate

: Energy ratings and minimum requirements for thermal envelopes and heating and air conditioning systems emerged as tools to minimize energy consumption and greenhouse gas emissions, improve energy e ﬃ ciency and promote greater transparency with regard to energy use in buildings. In Latin America, not all countries have building energy e ﬃ ciency regulations, many of them are voluntary and more than 80% of the existing initiatives are simpliﬁed methods and are centered in energy demand analysis and the compliance of admissible values for di ﬀ erent indicators. However, the application of these tools, even when simpliﬁed, is reduced. The main objective is the development of a simpliﬁed calculation method for the estimation of the energy consumption of multifamily housing buildings. To do this, an energy model was created based on the real use and occupation of a reference building, its thermal envelope and its thermal system’s performance. This model was simulated for 42 locations, characterized by their climatic conditions, whilst also considering the thermal transmittance fulﬁlment. The correlation between energy consumption and the climatic conditions is the base of the proposed method. The input data are seven climatic characteristics. Due to the sociocultural context of Latin America, the proposed method is estimated to have more possible acceptance and applications than other more complex methods, increasing the rate of buildings with an energy assessment. The results have demonstrated a high reliability in the prediction of the statistical models created, as the determination coe ﬃ cient (R2) is nearly 1 for cooling and heating consumption.


Introduction
Energy labelling of buildings and minimum requirements for insulation, solar control and performance of Heating, Ventilation and Air Conditioning (HVAC) systems emerged as tools to minimize the energy consumption and greenhouse gas emissions, improve the energy efficiency and promote greater transparency with regard to energy use in buildings [1].
The energy consumption of buildings is established by means of the relation between the energy demand and the performance of the mechanical space conditioning systems. The energy demand varies according to climatic conditions, the characteristics of the thermal envelope and the operative conditions of the building (occupancy, household appliance usage habits, equipment and lighting).
The first energy efficiency of building regulations emerged in the 1950s, in France, Switzerland and Germany, due to the concerns of the government entities to decrease the high heating energy Figure 1. Development of building rating systems in temperate countries. Own elaboration adapted from the Köppen-Geiger climatic classification [7].
In order to know the energy rating context in countries with temperate climates, 47 initiatives for regulation and energy rating, distributed among 27 countries, were revised (Table 1). Almost half (48%) of the 47 initiatives revised are compulsory and cover the national territory. It has to be highlighted that, while all developed countries studies have energy efficiency regulations and energy rating schemes, not all emerging countries have the tools to assess the energy efficiency of their buildings. Figure 1. Development of building rating systems in temperate countries. Own elaboration adapted from the Köppen-Geiger climatic classification [7].
In order to know the energy rating context in countries with temperate climates, 47 initiatives for regulation and energy rating, distributed among 27 countries, were revised (Table 1). Almost half (48%) of the 47 initiatives revised are compulsory and cover the national territory. It has to be highlighted that, while all developed countries studies have energy efficiency regulations and energy rating schemes, not all emerging countries have the tools to assess the energy efficiency of their buildings. Due to the analysis of the application of these regulations and schemes, two contexts are highlighted: the European Union and Latin America.
In the European Union, policies and strategies of each country member are conducted by the European Directives, the directives 2002/91/CE [8] and 2010/31/UE [9] being the most influential in the implementation of energy efficiency regulations and of the building energy rating scheme. They require European member states to achieve a particular result without dictating the means of achieving that result. Directives normally leave member states with a certain amount of leeway as to the exact rules to be adopted. Directives can be adopted by means of a variety of legislative procedures depending on their subject matter.
In Latin America, each country develops and implements the policies and strategies of its territory, as well as the regulatory framework for buildings. Not all Latin American countries have regulations focused on the energy assessment and energy efficiency of buildings. Furthermore, many regulations and energy rating schemes are voluntary [10].
The main compulsory initiatives in Latin America are the Brazilian Programme for Energy Labelling of Buildings (PBE Edifica) and the case of Rosario (Argentina). The first one requires an A level for new public buildings; the second one requires the fulfilment of certain values for thermal transmittance, risk of condensation, air permeability and cooling and heating thermal loads, established in the Ordinance 8757/11, in order to obtain the necessary licenses [10].
Regarding the assessed values of the building energy efficiency initiatives revised Table 1, in Latin America, 62% assess the achievement of certain requisites (thermal transmittance, solar factor, thermal capacity . . . ), 19% evaluate the energy demand and 19% focus on energy consumption. In Europe and the rest of the world, these proportions are different, with the main proportion of these initiatives centered on energy consumption, followed by the assessment of the energy demand, shown in Figure 2.
Due to the analysis of the application of these regulations and schemes, two contexts are highlighted: the European Union and Latin America.
In the European Union, policies and strategies of each country member are conducted by the European Directives, the directives 2002/91/CE [8] and 2010/31/UE [9] being the most influential in the implementation of energy efficiency regulations and of the building energy rating scheme. They require European member states to achieve a particular result without dictating the means of achieving that result. Directives normally leave member states with a certain amount of leeway as to the exact rules to be adopted. Directives can be adopted by means of a variety of legislative procedures depending on their subject matter.
In Latin America, each country develops and implements the policies and strategies of its territory, as well as the regulatory framework for buildings. Not all Latin American countries have regulations focused on the energy assessment and energy efficiency of buildings. Furthermore, many regulations and energy rating schemes are voluntary [10].
The main compulsory initiatives in Latin America are the Brazilian Programme for Energy Labelling of Buildings (PBE Edifica) and the case of Rosario (Argentina). The first one requires an A level for new public buildings; the second one requires the fulfilment of certain values for thermal transmittance, risk of condensation, air permeability and cooling and heating thermal loads, established in the Ordinance 8757/11, in order to obtain the necessary licenses [10].
Regarding the assessed values of the building energy efficiency initiatives revised Table 1, in Latin America, 62% assess the achievement of certain requisites (thermal transmittance, solar factor, thermal capacity…), 19% evaluate the energy demand and 19% focus on energy consumption. In Europe and the rest of the world, these proportions are different, with the main proportion of these initiatives centered on energy consumption, followed by the assessment of the energy demand, shown in Figure 2. According to Krstić and Teni [11], building performance evaluation models can be classified as black, grey or white box. The black box is the model with the least information required to the user and the white box is the model with the highest information required to perform the energy evaluation of the building. A black box model uses a simple mathematical or statistical model which relates a set of influential input parameters. The energy prediction under this model is based on a statistical approach according to a relevant database. According to Krstić and Teni [11], building performance evaluation models can be classified as black, grey or white box. The black box is the model with the least information required to the user and the white box is the model with the highest information required to perform the energy evaluation of the building. A black box model uses a simple mathematical or statistical model which relates a set of influential input parameters. The energy prediction under this model is based on a statistical approach according to a relevant database. A white box model is based on building physics, and requires a detailed description of the building and its thermal systems. What is known as the general method, based on energy simulations using DOE 2 or EnergyPlus, which calculates the dynamic energy performance of the building [ 12], can be included in this classification. This model uses complex tools, requiring a high level of expertise, and consumes lots of time and resources [13]. A grey box model mixes both above-mentioned methods, as it requires certain key parameters identified from a physical model and the energy prediction is based on statistical methodologies. In these models, a rough description of the building geometry is enough.
The input data for simplified methods (black or grey box models) are diverse. The analysis of 43 simplified methods focused on energy consumption (Table 2) shows that the most common input data is weather data (16), followed by characteristics of the building (13), energy bills (13), occupancy loads (11), characteristics of the thermal envelope (8), characteristics of the HVAC systems (5), building typology (5) and other inputs. In total, 23 of those methods are destined to residential buildings, 14 for office buildings, 5 for commercial buildings and only one is for any use. The simplified method of Petersdoff, Boernans and Hamish [14] is based on five building typologies with eight different levels of thermal insulation simulated on three climatic regions of Europe: cold, moderate and warm. Wong, Wan and Lam [15] proposed a model for office buildings with daylighting controls in subtropical climates, considering four climatic variables, four envelope variables and the type of day. Korolija, et al. [16] developed a model to predict the energy consumption of HVAC in office buildings in the United Kingdom, from four building typologies, seven envelope variables, five HVAC systems and five operational variables. Nath and Lamberts [17] established a simplified model to estimate the cooling energy consumption of an office building by simulating it on 18 Brazilian cities. They modified 15 building parameters, 12 HVAC system variables and three occupancy variants. The cities were selected based on a classification of the Brazilian local climates-cities-on 18 groups based on an index built from degree-hours and enthalpy.
The presence of general methods in Latin America is very low and their application is very scarce. The lower amount of required information and the lower time consumption of the simplified methods have allowed many Latin American national initiatives to be based on simplified methods. In fact, more than 80% of the Latin American initiatives are simplified methods and are centered on energy demand limitation or in the establishment of minimum requirements of the thermal envelope. However, the application of these tools is reduced.
It has to be highlighted that the energy assessment based on the energy demand or the characteristics of the thermal envelope does not integrate all the present heat interchange processes and does not account for energy consumption and greenhouse emissions, as the performance of HVAC systems are not taken into account.
The main objective is to develop a simplified calculation model to estimate energy consumption in order to assess the environmental impact of the building stock and to guarantee a higher acceptance and application in Latin America, with the objective of increasing the number of buildings with energy evaluations. This work is focused on multifamily housing buildings due to their high representativeness, especially in cities with a population greater than 200,000 habitants [18].

Materials and Methods
In order to create a simplified energy consumption estimation model for residential buildings in temperate climates for Latin America, the next steps were followed, as shown in Figure 3: 1.
Analysis of current energy building performance initiatives.

2.
Data collection: case study, user energy performance and present environmental conditions. 3.
Development and calibration of the energy model.

4.
Simulation scenarios: locations and U-value thresholds and proposals for walls and roofs.
Equations for energy consumption estimation.
The aim of the first step is the finding of the common parameters of the thermal envelope of the residential building energy performance initiatives for those countries of Latin America with cities located in temperate climates, in order to get the threshold values, which will conform the simulation scenarios.
The second step aims at gathering all necessary information needed to define and characterize the case study in order to build the energy model: the geometrical and constructive definition of the case study, the elaboration of the users' profile and the present environmental performance of the case study [19]. This case study, a residential building, is representative of this type of residence, as shown in the national office of statistics of each selected country.
1. Analysis of current energy building performance initiatives. 2. Data collection: case study, user energy performance and present environmental conditions. 3. Development and calibration of the energy model. 4. Simulation scenarios: locations and U-value thresholds and proposals for walls and roofs. 5. Energy consumption and linear regression. 6. Equations for energy consumption estimation. The aim of the first step is the finding of the common parameters of the thermal envelope of the residential building energy performance initiatives for those countries of Latin America with cities located in temperate climates, in order to get the threshold values, which will conform the simulation scenarios.
The second step aims at gathering all necessary information needed to define and characterize the case study in order to build the energy model: the geometrical and constructive definition of the case study, the elaboration of the users' profile and the present environmental performance of the case study [19]. This case study, a residential building, is representative of this type of residence, as shown in the national office of statistics of each selected country. The third phase compares the monitored environmental performance of the case study with the simulated results in order to calibrate the energy model. Once the energy model was calibrated, the simulation scenarios were defined which are the function of the locations of the selected countries with temperate climates and the legislative threshold values for common thermal envelope parameters found in the first step. Finally, five constructive scenarios representing the variety of these threshold values are defined.
In order to get the energy consumption and its variation from the simulation scenarios, the thermal model of the original case study is simulated for each location of each country and for each location changing the constructive scenario according to its national requisites.
The correlation between the energy consumption data and the climatic variables is studied due to the development of linear regression statistical models where the most relevant variables are highlighted.
Finally, once these variables are identified, the correlational equations for heating and cooling estimation are defined, forming a simplified method for energy estimation based on few climatic data, easily available for building agents.
These steps can be grouped in two: development of the experiment and statistical analysis. The first group, covering the steps 1 to 4, is exposed in this chapter and the statistical analysis of the energy consumption and development of the correlational equations, steps 5 to 7, are shown in the results chapter.

Analysis of Current Energy Performance of Buildings Initiatives
Focused on Latin American context, the actual initiatives for energy performance of buildings centered on residential buildings were prioritized, which include objective requisites about comfort and energy efficiency. The initiatives for residential buildings from countries with at least five locations with temperate climate and the availability of weather files for EnergyPlus (epw) were selected.
The analysis also includes the Spanish scheme for being a referent in Latin America due to historical linkages and the use of a common language, among other factors. Furthermore, climatic similarities with the temperate regions of Latin America and their weather files for energy simulations are available, enabling a comparison between the Spanish minimum requisites and the Latin American threshold values [20][21][22]. Table 3 shows the National Building Energy Rating schemes, the mandatory standard that limits the characteristics of the thermal envelope and the parameters limited. In bold are the common parameters of the thermal envelope of the selected schemes.  [26].
In the Argentinian context there is an energy rating scheme for residential buildings, the IRAM 11900:2017 standard, but this is of voluntary fulfilment and is based on the requirements of the Law 13059:2003 which forces different IRAM standards related to the characteristics of the thermal envelope to be applied.
The Brazilian Programme includes the Technical Regulation of Quality for the Energy Efficiency Level of Residential Buildings (RTQ-R [27]), which limits the values of the different parameters of the thermal envelope according to the bioclimatic zone of the location. The Chilean context is similar; the Building Energy Efficiency Rating is based on the limited values determined in a specific legislative document called Thermal Regulation [28]. The Mexican Building Energy Rating is mandatory for new housing and is based on the calculation program of PassivHaus, which requires the minimum levels specified in the Official Mexican Standards (NOM), also compulsory. The Spanish context is similar to the Brazilian, Chilean and Mexican contexts: there is a Building Energy Rating scheme which also allows for the certification of the fulfilment of the thermal requirements of the thermal envelope.
The common parameter of the thermal envelope that is limited by the selected schemes is the thermal transmittance (U-value) of the external walls (vertical and horizontal, i.e., walls and roofs). So, this parameter is selected to be modified and then to build different simulation scenarios.

Case Study
Nowadays, more than 80% of the population of Latin America and the Caribbean live in cities [29]. Data from the Economic Commission for Latin America and the Caribbean (ECLAC) [30] show that the percentage of the population that lives in cities with more than 20,000 inhabitants in Argentina, Brazil, Chile and Mexico is 79.8%, 70.6%, 80.3% and 70.2% respectively. Furthermore, according to the Economic Commission for Latin America and the Caribbean (CEPALSTAT) [31], the percentage of urban housing in these countries is 90%, 83%, 86% and 76% respectively, with approximately 75% on property.
Housing buildings are the most representative building typology in countries with temperate climates and are the greatest energy consumers during their useful life [18]. In Spain, 68% of residential units are located in multi-family buildings, being the majority in cities with more than 5000 inhabitants. It was found that 56.3% of residential units have no thermal insulation [32]. This percentage is lower in Latin America, but still of relevance; 57% of the Chilean residential buildings are located in metropolitan cities; of this buildings, 81% are multi-family residential buildings [32]. In Argentina, 21.38% of residential units are located in multi-family buildings in cities with more than 20,000 inhabitants [33]. In Brazil, the average amount of multi-family buildings is 35% [ 34]. A lower presence of multi-family buildings can be found in Mexico [35].
The case study is a multi-family building based on traditional building fabric, designated to middle class families according to the Brazilian Programme "Minha Casa Minha Vida" Figure 4.

Case Study
Nowadays, more than 80% of the population of Latin America and the Caribbean live in cities [29]. Data from the Economic Commission for Latin America and the Caribbean (ECLAC) [30] show that the percentage of the population that lives in cities with more than 20,000 inhabitants in Argentina, Brazil, Chile and Mexico is 79.8%, 70.6%, 80.3% and 70.2% respectively. Furthermore, according to the Economic Commission for Latin America and the Caribbean (CEPALSTAT) [31], the percentage of urban housing in these countries is 90%, 83%, 86% and 76% respectively, with approximately 75% on property.
Housing buildings are the most representative building typology in countries with temperate climates and are the greatest energy consumers during their useful life [18]. In Spain, 68% of residential units are located in multi-family buildings, being the majority in cities with more than 5000 inhabitants. It was found that 56.3% of residential units have no thermal insulation [32]. This percentage is lower in Latin America, but still of relevance; 57% of the Chilean residential buildings are located in metropolitan cities; of this buildings, 81% are multi-family residential buildings [32]. In Argentina, 21.38% of residential units are located in multi-family buildings in cities with more than 20,000 inhabitants [33]. In Brazil, the average amount of multi-family buildings is 35% [34]. A lower presence of multi-family buildings can be found in Mexico [35].
The case study is a multi-family building based on traditional building fabric, designated to middle class families according to the Brazilian Programme "Minha Casa Minha Vida" Figure 4. The structure is of reinforced concrete; the external walls are based on brick masonry with a plaster layer on both sides; the windows are simple glazed with aluminum frame without solar protection; the roof is based on a concrete slab and aluminum-zinc tiles. The U-values are 2.51 W/m 2 K for the external walls, 1.96 W/m 2 K for roofs, 5.80 W/m 2 K for windows and 2.88 W/m 2 K for the internal floors.
The building is located in a plot of land at the center of Criciúma, Santa Catarina, Brazil, being 28°41′ South its latitude. According to Köppen classification, its climate is Cfa (humid subtropical climate), characterized by hot and humid summers, and mild winters [7]. According to the Brazilian Standard NBR 15.220: Thermal Performance of Buildings [36], Criciúma is located in climatic zone II, characterized by an annual average daily temperature higher than 18 °C, as shown in Figure 5. The urban land is mainly plane, the surrounding is based on isolated multi-family buildings, so there is no vegetation or other buildings to obstruct the solar incident radiation or exposure to wind. The structure is of reinforced concrete; the external walls are based on brick masonry with a plaster layer on both sides; the windows are simple glazed with aluminum frame without solar protection; the roof is based on a concrete slab and aluminum-zinc tiles. The U-values are 2.51 W/m 2 K for the external walls, 1.96 W/m 2 K for roofs, 5.80 W/m 2 K for windows and 2.88 W/m 2 K for the internal floors.
The building is located in a plot of land at the center of Criciúma, Santa Catarina, Brazil, being 28 • 41 South its latitude. According to Köppen classification, its climate is Cfa (humid subtropical climate), characterized by hot and humid summers, and mild winters [7]. According to the Brazilian Standard NBR 15.220: Thermal Performance of Buildings [36], Criciúma is located in climatic zone II, characterized by an annual average daily temperature higher than 18 • C, as shown in Figure 5. The urban land is mainly plane, the surrounding is based on isolated multi-family buildings, so there is no vegetation or other buildings to obstruct the solar incident radiation or exposure to wind.

User Energy Performance
In order to define the occupancy and usage profile to be incorporated to the energy model, a survey of the inhabitants of the case study was carried out. A questionnaire was developed as a google formulary and sent by email to the responsible of each residential unit. From a population of 480 inhabitants, a response rate of 11% was recorded.
The focus of the survey was the understanding of those phenomena that could influence perceived comfort levels by recollecting information related to occupancy habits, HVAC preferences, possible modification of the internal covering, as well as perceptions about thermal comfort inside their dwellings.
The structure of the survey was composed of two introductory sections and four sections asking about occupancy, activities, timing, thermal comfort perception, air conditioning equipment and actions followed to reach thermal comfort Figure 6. The introductory sections were informative and exposed the objective of the survey, its extent and the expected time to fulfil it. The first group of questions was focused on the identification of the users' opinion about their perceived thermal comfort under two situations: using or not using HVAC equipment. The following section asked about their habits: how they occupy and use their dwelling and how they are dressed while they stay at home. The next section focused on the preferences for using HVAC equipment: type of HVAC equipment, frequency of use, operation mode, natural

User Energy Performance
In order to define the occupancy and usage profile to be incorporated to the energy model, a survey of the inhabitants of the case study was carried out. A questionnaire was developed as a google formulary and sent by email to the responsible of each residential unit. From a population of 480 inhabitants, a response rate of 11% was recorded.
The focus of the survey was the understanding of those phenomena that could influence perceived comfort levels by recollecting information related to occupancy habits, HVAC preferences, possible modification of the internal covering, as well as perceptions about thermal comfort inside their dwellings.
The structure of the survey was composed of two introductory sections and four sections asking about occupancy, activities, timing, thermal comfort perception, air conditioning equipment and actions followed to reach thermal comfort Figure 6.

User Energy Performance
In order to define the occupancy and usage profile to be incorporated to the energy model, a survey of the inhabitants of the case study was carried out. A questionnaire was developed as a google formulary and sent by email to the responsible of each residential unit. From a population of 480 inhabitants, a response rate of 11% was recorded.
The focus of the survey was the understanding of those phenomena that could influence perceived comfort levels by recollecting information related to occupancy habits, HVAC preferences, possible modification of the internal covering, as well as perceptions about thermal comfort inside their dwellings.
The structure of the survey was composed of two introductory sections and four sections asking about occupancy, activities, timing, thermal comfort perception, air conditioning equipment and actions followed to reach thermal comfort Figure 6. The introductory sections were informative and exposed the objective of the survey, its extent and the expected time to fulfil it. The first group of questions was focused on the identification of the users' opinion about their perceived thermal comfort under two situations: using or not using HVAC equipment. The following section asked about their habits: how they occupy and use their dwelling and how they are dressed while they stay at home. The next section focused on the preferences for using HVAC equipment: type of HVAC equipment, frequency of use, operation mode, natural ventilation and other passive strategies followed by the users. The last section asked about the constructive state of the dwelling in order to know if the thermal envelope had been modified.
The survey was based on multiple choice questions. An example can be found in Figure 7. Once the answers were received, they were analyzed in order to get statistical information about user behavior (an example is shown in Figure 8). Based on this statistical information, the user profile could be created.  The introductory sections were informative and exposed the objective of the survey, its extent and the expected time to fulfil it. The first group of questions was focused on the identification of the users' opinion about their perceived thermal comfort under two situations: using or not using HVAC equipment. The following section asked about their habits: how they occupy and use their dwelling and how they are dressed while they stay at home. The next section focused on the preferences for using HVAC equipment: type of HVAC equipment, frequency of use, operation mode, natural ventilation and other passive strategies followed by the users. The last section asked about the constructive state of the dwelling in order to know if the thermal envelope had been modified.
The survey was based on multiple choice questions. An example can be found in Figure 7. Once the answers were received, they were analyzed in order to get statistical information about user behavior (an example is shown in Figure 8). Based on this statistical information, the user profile could be created.  The surveys showed that two-bedroom dwellings were mostly occupied by two inhabitants; meanwhile four inhabitants were the majority occupancy for three-bedroom dwellings. The daily time duration of each user activity was 8 h for sleeping, 8 h for working, 2 h for cooking and eating, 2 h for cleaning and personal care and 4 h for entertainment.
Bathrooms, kitchen and common circulation areas were spaces where ventilation was the main strategy for thermal conditioning, meanwhile bedrooms and the living room used air conditioning equipment (if any). The HVAC systems on dwellings are air-to-air heat pump splits, fans and radiators. During summer, 100% of inhabitants used fans to reach thermal comfort at any time during the day, meanwhile 63% of the population use the split heat pump. However, 63% opened the window every time they wanted to ventilate the dwelling. Only 16% of the users maintained closed windows. During spring and autumn, 46% used fans, 9% used the split and 4% needed to use electric radiators. In these seasons there was a great portion (37%) which did not use any mechanical space conditioning system, and 24% kept their windows closed. During winter, those occupants with split systems used them (63%), 22% had and used electric radiators and 20% of the population did not use  The surveys showed that two-bedroom dwellings were mostly occupied by two inhabitants; meanwhile four inhabitants were the majority occupancy for three-bedroom dwellings. The daily time duration of each user activity was 8 h for sleeping, 8 h for working, 2 h for cooking and eating, 2 h for cleaning and personal care and 4 h for entertainment.
The users adapted their clothing according to seasonal variations, being tropical for summer (clo = 0,3), intermediate for spring and autumn (clo = 0,5) and light coated for winter (clo = 1,0). Bathrooms, kitchen and common circulation areas were spaces where ventilation was the main strategy for thermal conditioning, meanwhile bedrooms and the living room used air conditioning equipment (if any). The HVAC systems on dwellings are air-to-air heat pump splits, fans and radiators. During summer, 100% of inhabitants used fans to reach thermal comfort at any time during the day, meanwhile 63% of the population use the split heat pump. However, 63% opened the window every time they wanted to ventilate the dwelling. Only 16% of the users maintained closed windows. During spring and autumn, 46% used fans, 9% used the split and 4% needed to use electric radiators. In these seasons there was a great portion (37%) which did not use any mechanical space conditioning system, and 24% kept their windows closed. During winter, those occupants with split systems used them (63%), 22% had and used electric radiators and 20% of the population did not use any heating equipment or did not have it. In this season, 54% of the inhabitants kept their windows closed.
The dwellings are equipped with a fridge (24 h functioning), two PCs and two TVs that are utilized for 2 h each, and the shower is used for 1 h. Daylight is the main lighting source from 8:00 h to 18:00 h and there is a fluorescent lamp in each space which functions from 18:00 h to 22:00 h.
Surveys verified that inhabitants used adaptive comfort practices. The operation hypothesis, thus, has been assembled bearing in mind the strategies and actions that are allowed per adaptive The surveys showed that two-bedroom dwellings were mostly occupied by two inhabitants; meanwhile four inhabitants were the majority occupancy for three-bedroom dwellings. The daily time duration of each user activity was 8 h for sleeping, 8 h for working, 2 h for cooking and eating, 2 h for cleaning and personal care and 4 h for entertainment.
The users adapted their clothing according to seasonal variations, being tropical for summer (clo = 0.3), intermediate for spring and autumn (clo = 0.5) and light coated for winter (clo = 1.0).
Bathrooms, kitchen and common circulation areas were spaces where ventilation was the main strategy for thermal conditioning, meanwhile bedrooms and the living room used air conditioning equipment (if any). The HVAC systems on dwellings are air-to-air heat pump splits, fans and radiators. During summer, 100% of inhabitants used fans to reach thermal comfort at any time during the day, meanwhile 63% of the population use the split heat pump. However, 63% opened the window every time they wanted to ventilate the dwelling. Only 16% of the users maintained closed windows. During spring and autumn, 46% used fans, 9% used the split and 4% needed to use electric radiators. In these seasons there was a great portion (37%) which did not use any mechanical space conditioning system, and 24% kept their windows closed. During winter, those occupants with split systems used them (63%), 22% had and used electric radiators and 20% of the population did not use any heating equipment or did not have it. In this season, 54% of the inhabitants kept their windows closed.
The dwellings are equipped with a fridge (24 h functioning), two PCs and two TVs that are utilized for 2 h each, and the shower is used for 1 h. Daylight is the main lighting source from 8:00 h to 18:00 h and there is a fluorescent lamp in each space which functions from 18:00 h to 22:00 h.
Surveys verified that inhabitants used adaptive comfort practices. The operation hypothesis, thus, has been assembled bearing in mind the strategies and actions that are allowed per adaptive comfort model, that is, the changes of clothing for inhabitants and the operation of windows in order to get the dwelling ventilated, shown in Table 4 and Figure 9.   Figure 9. Occupancy hypothesis.
The environmental temperature range for comfort and HVAC was taken from the adapted Givoni's chart, with a minimum value of 18 °C and a maximum of 27 °C.
The air change rate was set to 2 ac/h due to the ventilation strategy of the users and the bad quality of permeability of the window frames. This value is in accordance with the suggested values of the standard IRAM 11.604/01 [37].

Present Environmental Conditions
In order to get the present environmental conditions of the building for the purpose of calibrating the energy model, seven dwellings were monitored. Dry bulb air temperature (insideoutside), relative humidity (inside-outside) and global exterior solar radiation were recorded for a time-step of 30 min from the 14th to 22th in January and from 10th to 18th in July, representing winter and summer conditions, both in 2016.
Inside thermal conditions were recorded by two HOBO UX100-003 data loggers per dwelling (in the main bedroom and living room), located 150 cm from the floor and protected from direct solar radiation, shown in Figure 10. The external temperature, relative humidity and global radiation were recorded with another HOBO data logger. The sensors were protected from direct solar radiation and precipitation. The environmental temperature range for comfort and HVAC was taken from the adapted Givoni's chart, with a minimum value of 18 • C and a maximum of 27 • C.
The air change rate was set to 2 ac/h due to the ventilation strategy of the users and the bad quality of permeability of the window frames. This value is in accordance with the suggested values of the standard IRAM 11.604/01 [37].

Present Environmental Conditions
In order to get the present environmental conditions of the building for the purpose of calibrating the energy model, seven dwellings were monitored. Dry bulb air temperature (inside-outside), relative humidity (inside-outside) and global exterior solar radiation were recorded for a time-step of 30 min from the 14th to 22th in January and from 10th to 18th in July, representing winter and summer conditions, both in 2016.
Inside thermal conditions were recorded by two HOBO UX100-003 data loggers per dwelling (in the main bedroom and living room), located 150 cm from the floor and protected from direct solar radiation, shown in Figure 10. The external temperature, relative humidity and global radiation were recorded with another HOBO data logger. The sensors were protected from direct solar radiation and precipitation. The measured values were drawn in the adapted Givoni's chart, as it takes into account the adaptability of the user to different temperatures and relative humidity ranges [38]. Although it is designed to include only outdoor air conditions, indoor measured values were drawn in order to confirm the strategies followed by users, as they had reflected on the surveys. The first comfort range (zone E) considers relative humidity values varying from 30% to 50%; the second range (zone F) considers an adaptation to relative humidity values of up to 80%. The temperature for a comfort zone with low humidity varies from 18.5 °C to 27.5 °C but varies from 17 °C to 26 °C if the relative humidity is high, shown in Figure 11. It can be observed that during the summer week, 3.24% of the measured hours in the living room were in thermal comfort zone E without any air conditioning strategy. Considering the natural ventilation and adaptation to high humidity, the percentage of hours in comfort increased to 28.70%. In the main bedroom, those values were 0.46% and 10.19%, respectively. It can be seen that the main The measured values were drawn in the adapted Givoni's chart, as it takes into account the adaptability of the user to different temperatures and relative humidity ranges [38]. Although it is designed to include only outdoor air conditions, indoor measured values were drawn in order to confirm the strategies followed by users, as they had reflected on the surveys. The first comfort range (zone E) considers relative humidity values varying from 30% to 50%; the second range (zone F) considers an adaptation to relative humidity values of up to 80%. The temperature for a comfort zone with low humidity varies from 18.5 • C to 27.5 • C but varies from 17 • C to 26 • C if the relative humidity is high, shown in Figure 11. The measured values were drawn in the adapted Givoni's chart, as it takes into account the adaptability of the user to different temperatures and relative humidity ranges [38]. Although it is designed to include only outdoor air conditions, indoor measured values were drawn in order to confirm the strategies followed by users, as they had reflected on the surveys. The first comfort range (zone E) considers relative humidity values varying from 30% to 50%; the second range (zone F) considers an adaptation to relative humidity values of up to 80%. The temperature for a comfort zone with low humidity varies from 18.5 °C to 27.5 °C but varies from 17 °C to 26 °C if the relative humidity is high, shown in Figure 11. It can be observed that during the summer week, 3.24% of the measured hours in the living room were in thermal comfort zone E without any air conditioning strategy. Considering the natural ventilation and adaptation to high humidity, the percentage of hours in comfort increased to 28.70%. In the main bedroom, those values were 0.46% and 10.19%, respectively. It can be seen that the main It can be observed that during the summer week, 3.24% of the measured hours in the living room were in thermal comfort zone E without any air conditioning strategy. Considering the natural ventilation and adaptation to high humidity, the percentage of hours in comfort increased to 28.70%.
In the main bedroom, those values were 0.46% and 10.19%, respectively. It can be seen that the main bioclimatic strategy followed by users to reach the comfort zone during summer is ventilation, as air movement reduces the perceived temperature.
During the winter week, 6.49% of the measured hours in the living room were in comfort zone E without passive strategies. As users in this period use artificial heating systems, and let solar radiation enter the dwelling but keep the windows closed, the humidity increases but within a certain range, and comfort hours increased this percentage to 49.80%. In the case of the main bedroom, those values were 4.24% and 38.89%, respectively. It can be observed how the users' strategies, obtained from the surveys, during winter have a higher impact on reaching the comfort zone, even with high humidity, than in the summer time.

Development and Calibration of the Energy Model
An energy model of the case study was created in DesignBuilder (version 5.0.0.137), based on the widely respected EnergyPlus simulation program [39]. The input climate data were the energy plus weather data (epw) [40], modified with the external measured data.
Geometric, constructive, occupancy, equipment, HVAC systems and usage characteristics were created according to the obtained data from the documentation of the building and the user profile created from surveys.
In order to validate the thermal model, a comparison between simulated internal temperature on free running and the measured values from monitoring Figure 12. bioclimatic strategy followed by users to reach the comfort zone during summer is ventilation, as air movement reduces the perceived temperature. During the winter week, 6.49% of the measured hours in the living room were in comfort zone E without passive strategies. As users in this period use artificial heating systems, and let solar radiation enter the dwelling but keep the windows closed, the humidity increases but within a certain range, and comfort hours increased this percentage to 49.80%. In the case of the main bedroom, those values were 4.24% and 38.89%, respectively. It can be observed how the users' strategies, obtained from the surveys, during winter have a higher impact on reaching the comfort zone, even with high humidity, than in the summer time.

Development and Calibration of the Energy Model
An energy model of the case study was created in DesignBuilder (version 5.0.0.137), based on the widely respected EnergyPlus simulation program [39]. The input climate data were the energy plus weather data (epw) [40], modified with the external measured data.
Geometric, constructive, occupancy, equipment, HVAC systems and usage characteristics were created according to the obtained data from the documentation of the building and the user profile created from surveys.
In order to validate the thermal model, a comparison between simulated internal temperature on free running and the measured values from monitoring Figure 12. It can be observed that, in summer, the maximum measured temperature was 31.1 °C, meanwhile the maximum estimated temperature was 35.3 °C. In winter, the maximum inside measured temperature was 20.3 °C and the estimated was 26.2 °C with a thermal time lag of 11 h. These values indicate that the energy model overestimates the temperature values, so it was necessary to identify the sources of uncertainty in order to calibrate the model.
Regarding the measured values, it was observed that the heterogeneity of the values was the main uncertain source. The total typical uncertainty of the measured values was 0.15 °C (95% confidence interval, coverage factor k = 2) indicating that the uncertainty is low and that the air stratification during the most critical summer and winter week is negligible.
The geometry of the model was the actual geometry, as it was measured from the case study. The properties of the materials were manually introduced based on the real construction. The weather file is a statistical compilation, so there is a difference between measured data and the data It can be observed that, in summer, the maximum measured temperature was 31.1 • C, meanwhile the maximum estimated temperature was 35.3 • C. In winter, the maximum inside measured temperature was 20.3 • C and the estimated was 26.2 • C with a thermal time lag of 11 h. These values indicate that the energy model overestimates the temperature values, so it was necessary to identify the sources of uncertainty in order to calibrate the model.
Regarding the measured values, it was observed that the heterogeneity of the values was the main uncertain source. The total typical uncertainty of the measured values was 0.15 • C (95% confidence interval, coverage factor k = 2) indicating that the uncertainty is low and that the air stratification during the most critical summer and winter week is negligible.
The geometry of the model was the actual geometry, as it was measured from the case study. The properties of the materials were manually introduced based on the real construction. The weather file is a statistical compilation, so there is a difference between measured data and the data contained in the weather files. Occupancy and internal loads were also a source of uncertainty, as the user profile was created based on the survey, which implies certain subjectivity.
A sensitivity analysis was carried out between three outputs and the parameters of thermo-physical properties of the materials, the general conditions of the building and its surrounding in order to assess the uncertainty of the energy model. The operative temperature was found to be the output which better reflects the modification of the parameters of the model. To adjust the model, the most influential parameters were the air change rate and the thermal conductivity of concrete in roofs and of the brick layer in walls.
The result of the energy balance showed that internal loads and solar heat gains were a significant heat source that provoked the estimated temperatures to be higher than the measured ones. Internal equipment loads and lighting were eliminated from the model and, due to the results from the sensitivity analysis, the air change rate was modified. Thermal properties of the roof and walls were not modified in order to keep the representativeness of the simulation scenarios. Calibration of the model was performed according to guidelines from ASHRAE Guideline 14 [41]. This guideline states that a model could be considered validated if its mean bias error (MBE) is no larger than 10%, and if the coefficient of variation of the root-mean-square error (CV(RMSE)) is not larger than 30% when the hourly data is used for the validation. The validation was based on the monitored and simulated dry bulb indoor air temperature, which was measured in 30-min intervals.
In each monitoring period, the MBE and CV(RSME) satisfied the 10% and 30% limits respectively for indoor air temperature, as recommended in ASHRAE Guideline 14.

Simulation Scenarios
In order to build a simplified calculation method and assess the energy consumption di fference due to the fulfilment of the minimum requisites of thermal transmittance defined by the Argentinian, Brazilian, Chilean, Mexican and Spanish initiatives, two simulation series were carried out. The first one consists of a simulation of the case study built in different locations with temperate climates. The second one consists of a simulation of the case study in these locations but modifying the thermal transmittance of the opaque elements of the thermal envelope according with the legal requisites of the country.

Locations and U-value Thresholds
In order to establish the locations with temperate climates for the first set of simulations, cities from Argentina, Brazil, Chile, Mexico and Spain with temperate climates and a population higher than 200,000 inhabitants were selected, totaling 307 locations. The information about population was obtained from national censuses.
As energy consumption is highly related to climate rigor, and this is related to degree days [ 42], cities with heating and cooling degree days higher than 1500 for each regime were selected, resulting in 42 locations, giving priority to those locations with a bigger population, more distance between them and with available weather files.
The selected locations account for 19.05% of the population of Argentina, 2.52% of the population of Brazil, 46.18% of the population of Chile, 13.53% of the population of Mexico and 8.69% of the population of Spain, so the results would cover a significant population from Argentina, Chile and Mexico. Table 5 shows the selected locations, heating and cooling degree days (HDD and CDD), national climatic zone classification and the minimum thermal transmittance for façades and roofs according to their national requirements for thermal envelopes.

Proposals for Walls and Roofs
Five couples of thermal transmittances (façade and roof) were created in order to deal with the requisites of each country. The first couple was established to be the U-values of the case study as the façade and roof materials and configuration are traditionally used in Latin America.
Normative thermal transmittance for the façade and roof for each location were identified Table 5. In Argentina, Brazil, Chile and Spain U-values vary depending on the climatic zone [22,28,36,43]. In Argentina, there is a difference between winter and summer requirements, so the most restricted values were chosen to fulfil both situations. In Mexico, U-values are only defined for certain locations as the normative does not specify a climatic classification [44]; furthermore, there is a change in U-values for buildings of more than three stories high, so the most restricted values were selected.
Finally, the admissible U-values were grouped based on similar ranges, excluding those values corresponding to the case study. Four couples of U-values for the façade and roof were defined in Table 6. These additional U-values present different materials in the main structure of each element. Excepting Wall II, all additional elements incorporate a thermal insulation layer. Walls are covered on both sides with a plaster layer, the lower side of the roofs are also covered with a plaster layer and the upper side of the slab is a layer of expanded clay.
In order to control the appearance of pathologies, the possibility of condensation risk was analyzed for each wall-roof couple [45]. It was verified that, under normal climatic conditions, there is no superficial or interstitial condensation risk in external walls and roofs.

Results
In this chapter, the statistical analysis of the simulation results and the development of the correlational equations are shown. Finally, in order to contrast the results of the equations, they have been compared with the energy consumption resulting from the energy rating schemes of the countries selected.

Energy Consumption and Linear Regression
The improved energy model was simulated in the 42 locations (scenario R) and also in the 42 locations modifying the U-values of the façade and roof according to national requisites, grouped in five couples (scenario N). The following energy consumption results were obtained, shown in Table 7. It can be observed that, in general, the modification of the original U-values for the walls and roofs in the requisite values of the national legislation, reduces the energy consumption for cooling and heating.
It also can be observed that, although locations within the same climatic classification generally trend to perform in a similar way (i.e., more heating than cooling consumption), there are cases where different trends can be found in the same climatic zone.
A multiple linear regression study was carried out in order to identify the climatic variables more related to heating and cooling energy consumption. Geographic and climatic data for each location were extracted by means of an analysis of the weather file with the Climate Consultant software including the altitude, latitude, monthly average temperature (Tavg), relative humidity (RH), global radiation (RAD), wind speed (Wsp) and sky cover percentage (SKcv).
Furthermore, degree days (DD) were calculated for each location by using the average temperature (Tavg) of each location considering a base temperature of 18 • C [46]. The minimum and maximum design temperature (Td) were also calculated. The maximum design temperature (Td max) was calculated by adding 3.5 • C to the average maximum temperature (Tmax avg) corresponding to the warmest month, and the minimum design temperature (Td min) was calculated by subtracting 4.5 • C from the average minimum temperature (Tmin avg) corresponding to the coldest month according to indications from IRAM 11.603 [45].
Monthly average temperature, relative humidity, global radiation, wind speed and sky cover percentage values were delimited for two different situations in order to determine the set of data more correlated with the variation of energy consumption. The first one considers the annual average; the second considers the average of the three warmest months (average summer-as) to be correlated to cooling consumption and the three coldest months (average winter-aw) to be correlated to heating consumption, shown in Table 8. The linear regression analysis was performed by SPSS software (version 15.0). The interpretation of the results was based on the consideration of the values of the determination coefficient. R 2 presents null correlation between their variables if the value is 0; very low correlation if the value is between 0.01 and 0.19; low correlation if the value is between 0.2 and 0.39; moderated correlation if the value is between 0.4 and 0.69; high correlation if the value is between 0.70 and 0.89; very high correlation if the value is between 0.9 and 0.99; and a perfect correlation if the value is 1 [47]. Positive values indicate that dependent the variable increases as the independent variable increases; negative values indicate that the dependent variable decreases as the independent variable increases.
The adjusted determination coefficient (R 2 ) was used to assess the reliability of the mathematical model. The relevance of the statistical model and the significance of the variables were verified for each case by means of the null p-value hypothesis test in the ANOVA and the probability values. A value ≤0.05 was adopted to test the hypothesis that the analyzed variable is of significance.
The method of multiple linear regression analysis was to add the variables forward, in which the software identifies the variables most correlated with the dependent variable (energy consumption), becoming the main variable. Lately, more variables are progressively added in order to increase the R 2 , provided that they are influencers of the dependent variable and improve the statistical model, becoming secondary variables. The latest model corresponds to the highest R 2 between the dependent variable and the main and secondary variables, becoming the optimal statistical model. The independent variables not included are neither significant nor influential on the dependent variable.
Results from the multiple linear regression analysis are shown in Tables 9 and 10. It can be observed that degree days are the main variable for heating and cooling energy consumption. The determination coefficients for both situations are 0.86 for heating energy consumption (winter) and 0.87 for cooling energy consumption (summer), indicating a strong variation of the energy consumption with degree days. The inclusion of secondary variables increased the correlation between energy consumption and climatic variables.
The final statistical model presents R 2 = 0.92 for summer and R 2 = 0.93 for winter (heating energy consumption). These results show that cooling energy consumption depends almost completely (92%) on the variation of cooling degree days, latitude and cooling degree days for the three warmest months. Heating energy consumption depends almost completely (93%) on the variation of heating degree days, global radiation of the three coldest months, average minimum temperature and average wind speed for the three coldest months.
Although heating and cooling degree days are the most influential variables on heating and cooling energy consumption, it is observed that not all of the five selected countries use degree days to define climate rigor, shown in Table 11, this is the case for Argentina and Mexico, which use the maximum design temperature and the average maximum temperature, respectively. In Mexico, furthermore, there is no standard for winter conditions, and in Chile there is no regulation on summer conditions. It is observed by analyzing the cooling energy consumption from the Chilean locations, shown in Table 7, that the climate in these cities has been demonstrated to be slightly rigorous during summer, so there is no need for a summer condition standard. However, heating energy consumption in Mexico was shown to be significant so it would be convenient to include any kind of regulation for winter conditions.

Equations for Energy Consumption Estimation
Once the statistical models were improved, the equations for the energy consumption estimation of multi-family residential buildings in temperate climates were defined (Equations (1) and (2) (1) where HEC is the heating energy consumption, CDD is the cooling degree days (base 18 • C), L is latitude (m) and CDD_s is the average cooling degree days for the three warmest months (base 18 • C). CEC = 394.912 + (0.164 * HDD) − (871.390 * RAD aw ) − (12.793 * T amin ) + 3.233 * W spaw (2) where CEC is the cooling energy consumption, HDD is the heating degree days (base 18 • C), RAD_aw is the average global radiation for the three coldest months (W/m 2 ), T_amin is the average minimum temperature ( • C) and W_spaw is the average wind speed for the three coldest months (km/h).
Once the equations are defined, the correlation between the results given by the equations and those obtained by the simulation, shown in Table 7, were compared in order to visualization of their adjustment, shown in Figure 13. It can be observed that the R 2 values are 0.92 for cooling and 0.94 for heating, indicating very high reliability in the estimation of energy consumption by the equations.

Contrasting Predictions from Equations and Established Methods
In order to identify the main similarities and differences between the equations and the existing methods to validate and optimize the model, a comparison between the equations and the Brazilian, Chilean and Spanish existing methods was carried out. The Argentinian and the Mexican method were excluded, as their procedures do not assess thermal loads.
Energy consumption results from equations were compared to the values obtained from the simplified calculation template provided by the PBE Edifica from Brazil and by the CEV from Chile. In the Spanish case, there are some tools, simplified and general methods. For this analysis, the simplified tool CE3X was used.
In the Spanish and Chilean cases, the methods enabled the energy demand to be obtained, so Equation 3 was applied in order to find the energy consumption. The seasonal average performance of the systems has been defined as 1.02 for summer and 1.45 for winter, according to the Institute for the Diversification and Saving of Energy (IDAE), from the Government of Spain [48]: where C is the final energy consumption (kWh/m 2 year), D is the energy demand (kWh/m 2 year) and η is the seasonal average performance (%) It was observed that the estimated heating energy consumption by the equation has a difference of 11.2% from CEV results, 21% from CE3x and 210% from PBE Edifica. The estimated cooling energy consumption by the equation has a difference of 43% and 596% from CE3X and PBE Edifica, respectively, shown in Table 12. It can be observed that the R 2 values are 0.92 for cooling and 0.94 for heating, indicating very high reliability in the estimation of energy consumption by the equations.

Contrasting Predictions from Equations and Established Methods
In order to identify the main similarities and differences between the equations and the existing methods to validate and optimize the model, a comparison between the equations and the Brazilian, Chilean and Spanish existing methods was carried out. The Argentinian and the Mexican method were excluded, as their procedures do not assess thermal loads.
Energy consumption results from equations were compared to the values obtained from the simplified calculation template provided by the PBE Edifica from Brazil and by the CEV from Chile. In the Spanish case, there are some tools, simplified and general methods. For this analysis, the simplified tool CE3X was used.
In the Spanish and Chilean cases, the methods enabled the energy demand to be obtained, so Equation (3) was applied in order to find the energy consumption. The seasonal average performance of the systems has been defined as 1.02 for summer and 1.45 for winter, according to the Institute for the Diversification and Saving of Energy (IDAE), from the Government of Spain [48]: where C is the final energy consumption (kWh/m 2 year), D is the energy demand (kWh/m 2 year) and η is the seasonal average performance (%) It was observed that the estimated heating energy consumption by the equation has a difference of 11.2% from CEV results, 21% from CE3x and 210% from PBE Edifica. The estimated cooling energy consumption by the equation has a difference of 43% and 596% from CE3X and PBE Edifica, respectively, shown in Table 12. The results indicate that the proposed method is closer to heating energy consumption calculated by the Chilean and Spanish tools and is farther from the cooling energy consumption calculated by the Spanish and Brazilian tools. When comparing the estimated energy consumption results from the equations (Table 7) with the degree days required for the 42 locations (Table 5), it is observed that although there is a very high correlation, for the Brazilian and Spanish cases, the energy consumption calculated by their national tools presents low or very low correlation with degree days, indicating that the proposed method is more consistent with the variation of the local climate rigor ( Figure 14).  The results indicate that the proposed method is closer to heating energy consumption calculated by the Chilean and Spanish tools and is farther from the cooling energy consumption calculated by the Spanish and Brazilian tools. When comparing the estimated energy consumption results from the equations (Table 7) with the degree days required for the 42 locations (Table 5), it is observed that although there is a very high correlation, for the Brazilian and Spanish cases, the energy consumption calculated by their national tools presents low or very low correlation with degree days, indicating that the proposed method is more consistent with the variation of the local climate rigor ( Figure 14).

Discussion
It is observed that in the Brazilian and Spanish systems, climatic conditions are introduced by means of a selection of the climatic zone (statistically) providing the same energy consumption for different locations within the same climatic zone (i.e., Blumenau, Chapecó, Criciúma, Florianópolis and Porto Alegre in Brazil); the input data in the proposed method are specific climatic data of the location, so it provides further approximation to real conditions. In this regard, the minimum difference to the calculated energy consumption by the Chilean system demonstrates a very high correlation between energy consumption and degree days for each selected Chilean city.
The geometry, building fabric, HVAC systems and user profile characteristics required by the energy rating schemes of the studied countries are similar to the input data of the energy model. However, in the national tools these are input data, while in the proposed method this information is implicit to the equations, so is not possible to modify the building characteristics.
It has to also be highlighted that while this method was built from a user profile developed from a real profile, in the Spanish tool the occupancy and usage profile is predefined for residential buildings with an occupancy thermal load much higher than that of the real profile obtained by surveys. The difference between these profiles is, without doubt, one of the main sources of difference between energy consumption results from the proposed method and the national tools.
The developed methodology is very close to the methods used by Petersdorff, Boermans and Harnisch [14], and Nath and Lamberts [17]. In both cases, an energy model was simulated under different climatic conditions in order to obtain energy consumption data for a statistical analysis. However, the proposed model works with 42 locations when the revised studies worked with three and 18 locations, respectively. Furthermore, the thermal envelope characteristics in the proposed model is centered in the modification of insulation thickness but Petersdor ff, Boermans and Harnisch [ 14] also vary building typologies, and Nath and Lamberts [17] modified 15 building characteristics without modifying insulation levels. HVAC characteristics are of higher relevance in the model of Nath and Lamberts [17], as it is centered on office buildings.

Conclusions
This research describes the methodology developed to create a simplified method to estimate the energy consumption of multi-family housing buildings located in temperate climates, whose input data are just a few climatic variables.
The estimated energy consumptions from the equations present very high values of the determination coefficient (R 2 = 0.92 and R 2 = 0.94), demonstrating the viability of the application of this method as a tool in quantifying energy consumption.
The proposed method is mainly manual, being easier to use than energy simulation software. It can be implemented in other tools as a spread sheet, in order for even easier use. It can be used by professionals during the building design or reconstruction stages helping to make decisions as it predicts the energy consumption for both heating and cooling, by research centered on the energy performance of buildings, allowing an assessment and qualification of the buildings, and by governmental actors to transfer information to populations in a simple manner, generating awareness about energy efficiency in buildings.
In contrast to other schemes that included climatic conditions based on climatic zones, the proposed method requires some specific local climatic variables as input data, so the energy consumption results are more consistent with the local climatic variations. There is a certain difference in the results from those calculated by national schemes based on climatic zones as they output the same energy consumption for different locations within the same climatic zone. In contrast, results from the proposed method and the Chilean one, as they vary with degree days, are very close.
The multiple linear regression analysis demonstrates that, in general, degree days is the most influential variable on energy consumption in residential buildings. In summer conditions, latitude and average degree days for the three warmest months are also influential. In winter conditions, the average global solar radiation for the three coldest months, average minimum temperature and average wind speed for the coldest months have to be taken into account.