Next Article in Journal
Planning Tourism in Protected Natural Areas: Safety, Soft Law and Conflict Management between Beach Users. The Case of Surf in Aljezur, Portugal
Next Article in Special Issue
Housing Evaluation Methodology in a Situation of Social Poverty to Guarantee Sustainable Cities: The Satisfaction Dimension for the Case of Mexico
Previous Article in Journal
Ferulated Pectins from Sugar Beet Bioethanol Solids: Extraction, Macromolecular Characteristics, and Enzymatic Gelling Properties
Previous Article in Special Issue
A Research Methodology for Mitigating Climate Change in the Restoration of Buildings: Rehabilitation Strategies and Low-Impact Prefabrication in the “El Rodezno” Water Mill
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of a Combination of Thermal Insulation and Cool Roof for the Refurbishment of Social Housing in Southern Spain

by
Carlos-Antonio Domínguez-Torres
1,
Helena Domínguez-Torres
2 and
Antonio Domínguez-Delgado
3,*
1
Escuela Técnica Superior de Arquitectura, Universidad de Sevilla, Avda. Reina Mercedes 2, 41012 Sevilla, Spain
2
Facultad de Ciencias Económicas y Empresariales, Universidad de Sevilla, Avda. Ramón y Cajal 1, 41018 Sevilla, Spain
3
Department of Applied Mathematics 1, Escuela Técnica Superior de Arquitectura, Universidad de Sevilla, Avda. Reina Mercedes 2, 41012 Sevilla, Spain
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(19), 10738; https://doi.org/10.3390/su131910738
Submission received: 24 August 2021 / Revised: 20 September 2021 / Accepted: 22 September 2021 / Published: 27 September 2021

Abstract

:
Social housing built in the middle of the last century in Spain suffers from poor thermal insulation conditions that cause situations of discomfort and energy poverty. For this reason, the energetic refurbishment of the envelope of this social building stock is necessary to overcome these situations and reduce energy consumption aimed at achieving interior comfort for its occupants. The goal of this work is to optimize a constructive solution that combines cool roof techniques with the use of thermal insulation applied to the refurbishment of the roof of buildings belonging to a quarter of social housing in Seville, Spain. The optimization analysis is based on the computation of the energy performance of the roofs when the energy retrofitting measure is applied, considering a variety of combinations of solar reflective coatings and insulation layer thickness, performing a dynamic analysis that accounts for the aging effect of the cool coats on the monthly roof energy performance and on the economic balance for the whole life cycle (LC) span. Economic and energy optimization analysis show that a suitable combination of cool roof emissivity and insulation layer thickness produces significant savings in the operational energy and in the economic profitability of the proposed retrofitting measure: the optimum combination obtained provides for the entire life cycle timespan an energy savings of 5.71 GJ/m 2 and a cost savings equivalent to the 63.1 % of the total costs when compared to the non-refurbished roof. The application of a time-dependent pattern for the changes on time produced by the aging effect on the cool roof emissivity, and its effects on the optimization of the combination of cool roof and insulation layer, can be considered novel in literature, both from an energy and an economic point of view.

1. Introduction

The building sector is responsible for around 40 % of final energy consumption in Europe and is a major source of CO 2 emissions. In order to reduce these emissions, in the last few decades, efforts have been made, consisting largely of the promulgation of regulations for the implementation of energy savings and efficiency measures, such as the Spanish norm DB-HE-2019 [1], which establishes new standards for the reduction of energy consumption in buildings, in accordance with the directives set by the EU for the H2030 [2], which establishes as goals for 2030 an improvement in energy efficiency equal to or greater than 32.5 % , cuts in greenhouse gas emissions equal to or greater than 40%, and a share for renewable energy equal to or greater than 32 % , all from 1990 levels [2].
In climates characterized by mild winters and very long hot seasons, such as the case of southern Spain and much of the Mediterranean area, a large amount of energy consumption is necessary to obtain the internal comfort of buildings during the cooling season. This poses a serious problem in the case of the social park housing built in southern Spain in the mid-twentieth century, in which envelopes often lack thermal insulation.
The energy efficiency of roofs has been object of a wide research focused on the enhancement and the optimization of its energy performance. This is largely justified by the fact that the roof is the envelope component that receives the most solar radiation in summer and in winter is responsible of an important percentage of the heat losses. Two widely used types of techniques aim to improve the thermal performance of roofs: those that modify their thermal properties, how to add or change thermal insulation, and those who carry out roof surface treatments, how to install radiant barriers or cool roofs. Taking into account that, in the geographical context considered here, in summer, most of the time, the sky is clear, and, in winter, the temperatures are not too cold, the cool roof technique seems advisable to lower the temperature of the roof, while the combination with a layer of insulation can help to achieve indoor comfort conditions and reduce energy consumption for this purpose.
Yu et al. [3] analyzed the optimum thickness of the thermal insulation installed on the roof of a residential building, and they studied different colors for the outer roof surface; to carry out the analysis, a life-cycle cost analysis (LCCA) was performed by considering solar-air degree-hours and four cities representing different China geographical regions; authors concluded that the effect of the colors on the optimal insulation thickness is different for each city. Gentle et al. [4] studied the energy behavior of a roof, analyzing different combinations of solar albedo, thermal emittance, and thermal resistance of a roof; a cost benefit assessment of the different combinations was performed that concluded the importance of high values for the roof albedo. Ramamurthy et al. [5] used the minimization of yearly costs for heating and cooling in order to determine the optimal combination of roof solar albedo and insulation thickness; they concluded that an albedo greater than 0.7, combined with a polyiso foam insulation of 18 cm, significantly reduces thermal flux through the roof when applied to new constructions; likewise, they found that, if installations costs are included in the analysis, the optimal combination is reached for a 4-inch insulation thickness with a payback period of 13 years. Arumugam et al. [6] studied the combination of solar reflective coating, radiant barriers, and thermal insulation for a roof by using EnergyPlus; a total of five climatic zones in India and 88 different roof combinations were analyzed; through an economic analysis, they found that the use of solar reflective coating and radiant barriers on the roof reduced the need for insulation layer for all the climatic regions analyzed.
Other studies have focused on optimizing the properties of the roof. Piselli et al. [7], by using dynamic energy computation, developed a method to perform an optimization analysis; the objective of this analysis was to optimize the solar reflectivity of the roof surface to minimize yearly energy loads for conditioning the building; the analysis was carried out by considering five cities of Italy located in different climate zones as the case study; from the performed optimization analysis, it was concluded that the the climate context is the most influential factor on the optimum roof solar reflectance.
Shi and Zhang [8] studied the effect on the reduction of annual energy loads of different solar reflectance and long-wave emissivity combinations for the outer surface of the building envelope; the study was carried out in 35 different climatic zones of the world and concluded that the greatest potential for energy savings is obtained for tropical climates and for high solar reflectance and high long-wave emissivity values.
Furthermore, a large number of studies deal with usefulness of cool roofs to improve indoor comfort conditions and to reduce energy consumption [9,10,11,12,13,14,15,16,17,18,19,20,21,22], while other studies analyze the benefits derived from the the attenuation of the heat island effect in cities, as well as achieving greater comfort in them [13,14,17].
It is worth noting that the performances of the building in summer and winter condition depend on the global behavior of the building, contending with global the performance of the roof and the infills. This means that the evaluation of thermal lack must also be extended to the performances of the infills. On this last topic, the literature is wide and complex, particularly for the fact that the nonstructural elements, such as the infills, must maintain the integrity under earthquake effect. Some works related with this topic are Kibert [23], Vailati and Monti [24], Manfredi and Masi [25], and Vailati et al. [26].
Likewise, some other recent works related to the research presented here and that may be of interest are the developed by Zingre et al. [27], Hernández-Pérez [28], and De Masi et al. [29].
Regarding the economic perspective, Levinson and Akbari [30] demonstrated the ability of cool roofs to achieve significant monetary savings when used in commercial building retrofits in the USA. Boixo et al. [19] stated that a large-scale implementation of cool roofs in the whole region of Andalusia, southern Spain, would have the ability to save yearly around 59 million euros in electricity costs, considering for this calculation only dwellings with flat roofs and the use of electricity as energy source, which would imply a reduction of about 136,000 metric tons on CO 2 emissions.
Nevertheless, as stated in Reference [31], in the literature, there are few studies that analyze the cost-effectiveness of cool roofs and, therefore, the combination of cool roof and insulating layer, based on an LCCA. Jo et al. [32] performed a research to investigate the impact on energy consumption and savings of replacing, in a commercial building, the original dark material of the roof with a solar reflective coating; after conducting a cost benefit analysis over a 20-year period, including installation and maintenance costs, the authors found monthly reductions of 1.3–1.9% and 2.6–3.8% in electricity consumption after a roof retrofit of 50 % and 100 % , respectively. Sproul et al. [33] performed an economic comparative analysis on the effect of some colors, specifically white, green, and black flat roofs in the United States; after a 50-year life-cycle analysis (LCA) that included installation, replacement, and maintenance costs, they concluded that the best economics results are obtained for the white color. Zhang et al. [34], after performing an LCA, demonstrated that a combination of roof ventilation and cool paint produces an annual energy savings for cooling of 109.94 USD/m 2 , with a payback period of 2 and 6 months for unventilated and ventilated roofs, respectively, under the weather conditions in Singapore. Hernández-Pérez [28] carried out a comparison of the cost-effectiveness of reflective white and colored roofs versus conventional gray roofs in six cities in Mexico. A 10-year LCCA concluded that, for the case of no thermal insulation, gray roofs are less cost effective than white and colored reflective roofs for all the considered locations. Yuan et al. [35] performed an optimization analysis of the combination of thermal insulation thickness and surface reflectivity through a 10-year LCCA; as result of the analysis, optimum combinations were proposed for each of the six considered regions in Japan. Saafi and Daouas [31] performed a 20-year LCCA in order to find the optimum combination of solar reflective coatings and insulation thickness for roofs of residential buildings in Tunisia; authors demonstrated that the use of these combinations are cost-effective and lead to economic savings of up to 44.53 Tunisien dirhams per m 2 with a payback period equal to 3.4 years; moderate values of the insulation thickness, combined with values of the solar reflectivity as high as possible, are recommended for use in the Tunisian climate; in this study, the effect of loss of reflectivity of the cool roof is taken into account when carrying out the life-cycle cost analysis; nevertheless, this effect is assumed to be constant throughout the whole life-cycle timespan.
The optimization analysis carried out in the present work tries to establish the cost-effectiveness of an Energy Saving Refurbishment Measure (ESRM) consisting of the combination of an insulation layer and a roof cool coat when used to refurbish roofs belonging to the social housing built in Spain in the middle of the past century and to optimize this combination both from an energy and economic point of view under the climate conditions of Southern Spain.
Taking into account the obsolete energy conditions of the building park considered, the present analysis is performed with the objective of providing guidelines to decision-makers for the implementation of refurbishment actions on the social housing aimed at reducing situations of internal thermal distress, energy poverty, and associated unhealthy situations. Thus, the improvements on energy efficiency that result from the optimization of the combination of thermal insulation and cool roofs for the refurbishment of social housing can be framed in Buildings Regeneration Strategy to Climatic Change Mitigation, Energy, and Social Poverty.
Likewise, the present optimization analysis is proposed in order to appraise in a deeper way the combined effect of insulation layer thickness variability and different values of solar reflectivity for the roof coat and the effect on this combined system, together with the aging effect of the cool paint on the energy performance of the refurbished roof, as well as their implications on the economic and energy LCCA of these combinations when they are implemented in the retrofitting of the social housing park under consideration.
To carry out the optimization analysis, the integration of the aging pattern for the cool coat with the proposed variability in the combinations of insulation layer thickness and roof reflectivity requires a specific simulation code allowing the accurate calculation on thermal loads for the refurbished roofs. Usually, commercial packages for building energy simulation (BES) do not have this ability. This software limitation may be the cause of the lack of studies that address the time-dependent loss of solar reflectivity due to the aging effect of the cool coat. To the best of our knowledge, there exists a gap in the literature with regard to studies that perform optimization analysis by considering annual variable costs due to the aging effect, when performing the LCCA of the proposed combination. Such a gap is especially important in the geographic area analyzed in this research where cooling demand is usually very high in the hot season.
This way, the thermal performance of a combination of insulation layer and cool roof is analyzed through a comprehensive numerical study when such a combination is used to refurbish roofs from residential buildings belonging to the social housing built in Seville, Spain, in the middle decades of the 21st century before the enactment of the first Spanish legislation intended at regulating energy demand in buildings, the NBE-CT-79 [36], in 1979.
The main objective of this work is to perform an energy and economic optimization of the proposed refurbishment measure within a broad framework considering the effect on the LCCA of the real weather of Seville, a wide range of combinations for the insulation layer thickness, and the cool coat reflectivity, as well as the aging effect of the cool coating, the energy source, and equipment type used to condition the buildings and current standard economic parameters in Spain and Europe.
Finally, it is important to note that the findings of the present research can be extended to other places with climatic conditions similar to those existing in the south of Spain, as are as are almost all southern European regions.

2. Methodology

Compared to previous research aimed to estimate the optimum combination of cool roof reflectivity and insulation thickness, the methodology followed in this work takes into account the time-dependent loss of reflectivity of cool paint produced by the aging effect and the cyclical restoration of the solar reflectivity by power washing of the cool roof. In this way, in order to estimate the impact on thermal loads of the evolutionary change in roof reflectivity in combination with different insulation thickness, a dynamic finite difference model is developed that calculates the monthly thermal loads for the entire LCA timespan, taking into account the monthly updated reflectivity estimated following the patterns shown in Reference [37].
The LCCA is performed by taking into account the costs of surface reflective coating material, its installation, and the power washing when done, as well as the costs of the insulation material and its installation. To compute the costs for the whole LCA span, the energy cost is computed monthly for every studied scenario, considering the changes in the energy consumption produced by the cool coating reflectivity changes along the LCA span, and then the costs are estimated under the assumption of the use of three different sources of energy to obtain thermal comfort: heat pumps for heating and cooling, air conditioners to cool and gas boilers to heat, and, finally, air conditioners to cool and electrical cheating.
The thermal load calculations are used to find out the most advantageous combination of roof surface reflectivity and insulation thickness in order to reduce the energy load needed to obtain indoor comfort in the studied buildings and the considered climate framework for the whole LCA time period. This way, the economic analysis based in the energy results is used to conclude the optimum combination that provides the minimum total cost for the LCA timespan.

2.1. Energy Roof Model

In the energy analysis, the hypothesis of one-dimensional heat conduction through the roof is considered. It is assumed that the roof is made up of M parallel layers meeting the following hypothesis: the material of every layer is homogeneous and isotropic, thermal properties of materials remain constant under temperature changes, and there is no internal heat generation. Then, for the multilayer heat conduction through the roof, the linear heat conduction equation,
T j ( x , t ) t = α j 2 T j ( x , t ) x 2 , j = 1 , , M ,
can be applied in each layer j, with T j being the temperature in layer j, and α j the thermal diffusivity of the layer material j [38].
Equation (1) requires the definition of boundary conditions. The boundary conditions for the external and internal surfaces of roof are obtained from the energy balance equation corresponding to each surface.
The outer surface of the roof is affected by the solar radiation I S W , the long-wave exchange with the sky q L W , the heat exchange by convection between the outer surface and the ambient air q c , e x t , and, finally, by the heat conduction through the outer layer. Thus, the energy balance for this surface can be written as:
k T n + q c , e x t + α S W · I S W + q L W = 0 ,
where n is the outward normal vector to the roof outside surface, k is the thermal conductivity of the outer layer material, and α S W is the solar absorptivity of the external roof surface.
The external convection heat transfer q c , e x t is computed as
q c , e x t = h c , e x t ( T a T ) ,
where T a is the outdoor air dry-bulb temperature, T the temperature of the outer surface of the roof, and the correlation
h c , e x t = 8.18 + 2.28 V R [ W / m 2 K ]
proposed by Hagishima and Tanimoto [39], based on experiments carried out on a roof, has been adopted. Here, V R is the wind velocity above the roof in m/s.
The long-wave radiation exchange between the roof external surface and the sky is calculated as:
q L W = Q s k y L W Q w L W ,
with Q w L W being the thermal radiation intensity emitted by the external surface of the roof, and Q s k y L W being the sky down-welling long-wave radiation.
For computing Q w L W , the law of Stefan–Boltzmann establishes that
Q w L W = ε σ T 4 ,
where ε is the surface emissivity, T the surface temperature in Kelvin degrees, and σ = 5.67 × 10 8 [W/m 2 K 4 ] is the Stefan–Boltzmann constant.
Q s k y L W is estimated through the expression:
Q s k y L W = ε s k y σ T a 4 ,
where ε s k y is the sky emissivity that, following Walton [40] and Clark et al. [41], has been computed as:
ε s k y = ( 0.787 + 0.764 ln ( T d p 273 ) ) ( 1 + 224 10 4 n 35 10 4 n 2 + 28 10 5 n 3 ) ,
where T d p is the dew point temperature in Kelvin, and n is the fraction of sky covered by clouds in tenths.
Other sources of thermal radiation, such as plant masses or nearby buildings higher than the building where the studied roof is placed, have not been considered in the present study, for the sake of brevity.
For the roof inside surface, the energy balance is given by
κ T n + q c r , i n t = 0 ,
where T is the roof interior surface temperature, n is the outward normal vector to the roof interior surface, κ is the thermal conductivity of the roof internal layer, and q c r , i n t is the combined convective-radiant heat transfer balance between the roof interior surface and the indoor temperature computed as:
q c r , i n t = h c r , i n t ( T T r o o m ) , ( W / m 2 ) .
Here, h c r , i n t is the combined convective-radiant heat transfer coefficient defined as the sum of the convective h c o n v , i n t and the radiant h r a d , i n t heat transfer coefficients. Following the correlations established in Reference [42] for horizontal surfaces, the convective heat transfer coefficients are taken as h c o n v , i n t = 0.92 (W/m 2 ) for downward heat flow and as h c o n v , i n t = 4.04 (W/m 2 ) for upward heat flow. The radiant heat transfer coefficient is taken as h r a d , i n t = 5.72 ε (W/m 2 ), where ε is the emissivity of the roof inside surface. These values are recommended by ASHRAE [43], and they are often used in the calculation of heat transfer in buildings [44,45].
Finally, the M roof layers are assumed to be in good contact; therefore, the interlayer thermal resistance is neglected.

2.2. Changes in Reflectivity Due to Aging Effect and Maintenance

To consider realistic conditions, the aging effect of the cool coating must be considered. Some authors have pointed out that the factors that result in a loss of reflectivity in cool paints are mainly the weatherization of the cool roofs, the characteristics of the paint, and the material on which the paint is applied [22,46]. Other studies quantified the lost reflectivity and the rate at which this loss occurs. This way, in Reference [47], a loss of the solar reflectivity for white roofs of about 10–30% was found; this reduction was found mostly during the first year, whilst, in Reference [48], it was established that the solar reflectance of a cool paint decreased about 11% after being subjected to 400 h of artificial accelerated weathering. In Bretz and Akbary [49], it is stated that the decrease in the solar reflectivity of the the cool coat reflectivity occurred mostly in the first year, and, very likely, in the first months, and, in Reference [50], from an in-field study, it was shown that, for all reflective roofing types analyzed, in the first two years, about 95% of the loss of solar reflectivity occurred, and about 98% of the total aging effect took place in the interval of three years after the installation.
Finally, from in-field measurements [49,50], it has been found that a power washing of the cool roof could restore up to 90–100% of initial values for the solar reflectivity of the cool coat. In Reference [49], the cold coat was washed one year after its installation, and, in Reference [50], no data is provided on the period of time in which the wash was performed.
According to the aging behavior mostly supported by previous literature, to compute the heat flux through the roof, the following pattern for the reflectivity loss of the cool coating has been considered: during the first year, the solar reflectivity suffers a loss of the 20 % of its original value, with most of this reduction occurring in the first few months; in the second year, a reflectivity loss equal to 8.5 % is taken, and, during the third year, a loss equal to 0.9 % is considered. After these first three years, the solar reflectivity decrease is assumed to stabilize [49], until a loss of 30 % is reached for the whole LC timespan. With these assumptions, the aging pattern stated in Reference [50] is fulfilled. On the other hand, according to the aforementioned literature, the solar reflectivity is considered to be restored to 90% of its initial value after a power washing and then following the same rate of aging loss that occurs after the first application of the cool layer.
Finally, the pattern established in Reference [37] for the time evolution of the solar reflectivity, that collects the foregoing, is used in the computations. This pattern is illustrated in Figure 1 and Figure 2. In Figure 1, the evolution of the solar reflectivity of a cool roof with an initial value equal to 0.9 is shown for a triennial washing maintenance.
In Figure 2a, the monthly and accumulated loss of the solar reflectivity is shown for a period of 42 months, and, in Figure 2b,c, the assumed evolution of the solar reflectivity without maintenance and for a decennial washing, respectively, are shown for a period of 20 years.
An immediate consequence of the evolution over time of the exterior reflectivity is the need to calculate the heat transfer through the roof on a monthly basis for the whole LCA timespan. For this, the reflectivity value in Equation (2) must be updated for each moment of time, following the aging pattern described above. That is, the heat transfer through the roof cannot be considered the same for all the years of the LCA timespan, which is a frequent assumption in the literature.

2.3. Thermal Load Calculations

The key point of the developed analysis is the accurate computation of the heat transfer through the roof to the building indoor and the associated thermal loads to obtain indoor comfort. To do that, the load at time t is calculated using the expression entered above:
q ( t ) = h c r , i n t ( T ( t ) T r o o m ) ( W / m 2 ) ,
which takes into consideration the convective and radiant thermal flux from the internal roof surface to the building indoor. q ( t ) yields positive values when the heat flux enters to the building, and negative values when heat flows out of the building.
Then, if weather hourly values are considered, as it is usual in building energy software, the hourly load is derived from (7) for every hourly interval. This way, daily, monthly, and yearly loads are straightforwardly obtained by integrating the hourly loads on the timespan of interest.
The thermal loads due to the heat transfer through the roof, including the cooling and heating loads, are computed varying the absorptivity of the cool paint from 0.1 to 0.5 , with this being the usual value range to consider a roof as a cool roof [51]; for computations, steps equal to 0.1 are considered for the absorptivity. The insulation thickness ranges from 0 to 10 cm, with the latter being the maximum value usually used for the roof retrofitting in southern Spain. For the computations, steps of 1 cm for the thickness of the insulation layer are considered.
We assume a fixed value of the indoor temperature in each season of the year. The conditioning system is considered to be in continuous mode with a fixed temperature set point at 24 C in the cooling season, at 20 C in the heating season, and at 22.5 C for the intermediate seasons, according with the comfort temperatures established in the Spanish regulations for thermal installations in buildings [52].
In order to estimate the energy loads due exclusively to the heat transfer through the roof, we follow the approach from References [4,31,53], and, then, the zone beneath the roof is assumed to have the envelope adiabatic, except for the roof.
The input data needed to compute the thermal loads are:
  • Meteorological data: Air temperature, solar radiation, cloud cover, wind speed, and relative humidity. These data can be hourly or any other frequency smaller than one hour.
  • Roof data: solar and thermal reflectivity of the outer roof surface and geometrical and thermophysical properties of the roof layers.
  • Indoor data: indoor air temperature.
On the other hand, with regard to the aim of assessing the energy and economic performance of the refurbishment measure through the lifetime period of 20 years according to the service time of the roof coating, taken into consideration is the effect stemming from the action of environmental agents as rain, dust, air particles, moisture, or sun, as well as their own aging, on the initial reflectivity properties of the cool paint, as described in Section 2.2.
To estimate the thermal loads, the following cool coating maintenance scenarios were analyzed: no washing, and decennial, quinquennial, and triennial washing. Taking into account the social character of the buildings in which the roof reform is applied, it seems appropriate to consider an annual roof wash unlikely; realistically, only no maintenance, or washing every five or ten years, scenarios should be considered.
The thermal loads were estimated using the open source FreeFem++ software [54]. Hourly loads, including heating and cooling loads, were computed for the whole LC time interval.

2.4. Case Study

The case study considered is a building in the city of Seville. The city is located in the Guadalquivir valley in southern Spain; see Figure 3.
Its location in a valley with a low altitude above the sea, its meridional latitude, and its dry climate, together with the high levels of solar radiation, cause a long hot season, spanning months beyond summer which, together with its soft winters, make this city an appropriate place for the application of cool roofing techniques.
The climatic chart of Seville is shown in Figure 4 from data from the Spanish State Meteorological Agency [55]. The average annual temperature is 19.2 C, with maximum average temperatures up to 40 C in the months of July and August and a minimum average of 5.7 C in January. The normal incident solar radiation reaches a maximum daily mean equal to 8.3 kWh/m 2 in July and a minimum daily mean equal to 2.3 kWh/m 2 in December. Winters are mild while summers are warm, dry, and sunny. On the other hand, autumns and springs are characterized both by moderate temperatures and by being the seasons in which most of the rainfall occurs. According to the described characteristics, the studied area climate is classified as Mediterranean, following the Köppen-Geige climate classification.
It is noteworthy to observe the temperature differential between the temperature of the sky and the ambiance temperatures. As can be noted in Figure 4, this difference is up to 15 C throughout the year with respect the mean ambiance temperature, which is an indicator of the potential of the longwave heat exchange with the sky to reduce the roof surface temperature.
The studied building belongs to a quarter of social housing called La Juncal Figure 5 built in the sixties of the last century that exhibits the constructive characteristics common in the social housing built in Seville before 1979, when the first Spanish legislation designed to regulate energy demand in buildings was enacted.
In Figure 6a, the layout of a typical roof belonging to the social housing studied is shown, and, in Table 1, the dimensions and thermophysical values of the components of the roofs are reported. This constructive configuration for roofs is the usual in the social housing built in the sixties of the past century at southern Spain [56]. For the outer surface of the roof, the absorptivity solar radiation value is taken equal to 0.8, with the latter being a typical value of the material making up this layer. This roof will henceforth be called the reference roof.

2.5. Energy Saving Refurbishment Measures

In order to improve the energy performance of the reference roof, a combined refurbishment measure was applied to the outside layer of the roof. This measure consisted of the installation of an EPS insulation and the application of an external layer of white elastomer cool paint to the outer surface. In Figure 6b, the layout of the retrofitted roof is shown, and, in Table 2, the dimensions and thermophysical values of the components of the retrofit are reported. The combination of the insulation thickness and the value of the cool paint absorptivity is then analyzed to optimize the energy and economical performance of the roof, taking into consideration the time variation of the cool coating reflectivity discussed above.

2.6. Model Validation

In order to verify the reliability of the numerical model used in the present study, a twofold process of validation is performed: first, the model outputs are compared to the experimental values obtained from the monitoring of a full-scale outdoor test cell; and, second, the calculated values from the numerical model are compared to the values provided by a well known and validated building energy simulation (BES) tool, EnergyPlus, when applied to the case-study used here and described in Section 2.4.
The indices used to check the reliability of the numerical model were the root mean square error (RMSE), the mean of the residuals ε m e a n , the standard deviation of the residuals S D ( ε m e a n ) , and the r 2 index, which are often used in the process of validation of building energy simulation models [57].
For both validation procedures, a winter week, from 7 to 14 February 2017, and a summer week, from 17 to 24 August 2017, were considered. The climatic data for both weeks were local data recorded by a meteorological station placed on the roof of a test cell in the city of Seville, as explained in Section 2.6.1. The indoor conditions assumed for the thermal load calculations, as in Section 2.3, are considered: An ambient interior temperature kept constant to 20 C in the winter week, and to 24 C in the summer week. In Figure 7, the ambiance temperature and solar radiation for both weeks are shown.

2.6.1. Experimental Validation

For the experimental validation, an outdoor full-scale test cell placed in the center of the city of Seville was used. The interior dimensions of the test cell are 2.40 m wide, 3.20 m deep, and 2.70 m high. The thermophysical values of the materials and the roof stratigraphy are presented in Table 3. The values shown in the table are the nominal ones provided by the Spanish Technical Building Code [58] and the material specifications sheets from the manufacturers. The outer surface of the roof is painted red, and a solar absorptivity of 0.5 and a thermal emissivity of 0.9 are considered for it.
The test cell is fully instrumented to measure the temperatures of the internal surfaces of the envelope, and data are recorded every 10 min. Thermocouples with an accuracy of ±0.75 C and an operating range of 250 to 350 C are used to monitor the temperature of the inner surfaces. Finally, a local meteorological station located on the roof of the cell is used to record weather data. More details on the test cell and its monitoring equipment can be found in Reference [59].
In Table 4, the values of the statistical indices computed for the temperature of the inside roof surface, are shown.
The results obtained for the indicators shown in Table 4 allow for asserting that the presented model is able to compute the temperatures of the internal roof surface with good precision; the values obtained for the different indicators are similar to those obtained in other works for the validation of energy numerical models for a building’s envelope [57].

2.6.2. Comparison with EnergyPlus Software

The EnergyPlus [60] software is a well known BES tool that has been widely used for the calculation of the energy performance of buildings, proven to be very efficient in the computation of the thermal behavior of buildings. An example of the application of EnergyPlus to the analysis of the thermal behavior of roofs can be found in Reference [31].
To carry out the comparison between the numerical model presented here and the EnergyPlus tool, the reference roof of the case study in Section 2.4 and described in Table 1 is used, and the values for the temperature of interior roof surface provided by EnergyPlus are considered as benchmark.
In Table 5, the values of the statistical indices computed for the temperature of the indoor roof surface, are shown.
The values of the indices shown in Table 5 allow to stablish the good fit between the values from the numerical method and from the EnergyPlus tool. As in the case of the experimental validation, these values are similar to those obtained in other validation processes of numerical thermal models for building envelopes [57]. On the other hand, it can be observed that the indicators for the comparison between the models have a slightly better performance than for the experimental validation, a fact also reported in the literature.

3. Energy Optimization Analysis

In this section, an analysis of the energy performance of a variety of roof solar reflectivity and insulation thickness combinations is carried out. The objective of this section is to determine the combination that provides the minimum energy consumption for the whole LCA timespan and to calculate the energy savings of each combination when compared to the reference case.
The analysis is carried out considering a time dependent change of the cool coating solar reflectivity due to the aging effect and different maintenance patterns, as described in Section 2.2.
The solar absorptivity of the cool coat varies from 0.1 to 0.5, with steps equal to 0.1. The insulation thickness ranges from 0 to 0.1 m, with steps equal to 0.1. For the outer surface of the retrofitted and the reference roofs, the thermal emissivity taken is equal to 0.9. The considered maintenance for the cool roof was a power washing every three, five, or ten years. The case where no power washing is performed was also included in the analysis.
The transmission loads through the roofs are calculated for the entire LC timespan taking into account the annual variability of the thermal loads during that period of time due to the changes in time of the cool layer reflectivity produced by the effects of aging and maintenance patterns.

3.1. Thermal Loads Results

In Figure 8 and Figure 9, the total thermal loads (cooling plus heating loads) for the whole LC timespan are shown for different values of solar absorptivity, thickness insulation, and maintenance protocols, together with the total load for the reference roof. As can be seen in these figures, the effect of low absorptivity values on the total thermal loads is very noticeable as it is observed when comparing with the reference case, even in the case of the retrofitted roof without insulation. For all the combinations and all maintenance protocols considered, the total loads are lower than for the reference case, although there is a significant difference between the loads for the uninsulated roof and the insulated ones.
On the other hand, it is observed that the lower the solar absorptivity, the lower is the total load. It is also observed that increasing the thickness of the insulation reduces the loads. So, for every value of the solar absorptivity, the total load decreases as the insulation thickness increases. However, regarding the maintenance protocol, only little differences are found among the total loads, with a slight tendency to decrease as the frequency of the cool coat washes increases. Thus, it can be drawn up that the role played by periodic washing in the total thermal loads is relatively modest.
Analyzing the heating and the cooling loads separately, it is clear that, for all the cases of absorptivity, thick insulation, and maintenance protocols, the smaller the cool roof absorptivity, the smaller the cooling load, as can be seen in Figure 10 and Figure 11. In addition, it can be noted that the smallest cooling loads are obtained for the case of the triennial wash, which is in accordance to the fact that the original roof solar absorptivity values are restored with a higher frequency than in the other maintenance protocols. Therefore, the cooling loads gradually increase as the frequency of cool roof maintenance decreases and reaches the highest values for the case of no power washing.
On the other hand, it can be seen that increasing the thickness of the insulation has a strong effect in reducing the cooling loads, so that the thicker the insulation layer, the lower the cooling load. It can be concluded that, for all the maintenance protocols, the lowest load is obtained for the thickest insulation and the lowest solar absorption. Regarding the maintenance protocol, the most favorable case is the triennial power washing.
Regarding the heating loads, the behavior is the opposite to that described for the cooling loads, as can be observed in Figure 12 and Figure 13. Now, due to solar gains decrease, the retrofitted cool roof without insulation has higher heating loads than the reference case. Clearly, the winter penalty effect on solar gains due to the cool roof is the cause of this. For the roofs with insulation, the only case where the heating loads are greater or close to that of the reference roof is when solar absorptivity equals 0.1. For values greater than 0.1 of the roof absorptivity, the heating loads are lower when compared to the reference case. This result holds for all the insulation thickness and maintenance protocols. Therefore, it can be stated that, although the heating loads increase as the solar absorptivity of the roof decreases, the insulating layer is able to compensate the penalty effect due to the low values of absorptivities, except for the value of the absorptivity equal to 0.1, as discussed before. Again, the role of the insulation layer in reducing the heat flux is obvious, i.e., the thicker the insulation, the lower the heating load. Finally, it is again observed that the effect of the cool coat washing is quite small on the heating loads over the LC timespan.
This opposite behavior of the heating and cooling loads described for the retrofitted roofs is the cause of the small differences observed among the total thermal loads for the whole LC timespan when considering different maintenance scenarios, while the changes in solar reflectivity produced by the aging of the cool coat act in the opposite direction. That is, on the one hand, when power washing is carried out, it produces a decrease in the cooling loads and an increase in the heating loads, while the aging of the cool coat produces the opposite effect in such a way that all these processes combined they result in a certain equilibrium in the resulting total loads. To take into account the aging of the cool coat, the aging pattern discussed in Section 2.2 is used for the thermodynamic computations by using the monthly reflectivity given by the pattern.
In Figure 8 and Figure 9, the most remarkable fact, from a physical point of view, is the difference between roofs equipped with thermal insulation and roofs without it; this is due to the higher heat flux values through roofs without thermal insulation, which implies, in the cold season, a greater loss of heat from the interior and, in the hot season, a greater gain, so that these effects produce a significant thermal load increasing for the roofs without insulation when compared to the equipped with insulation ones. On the other hand, it can be observed that, for the roof without insulation, the reduction of heat gain due to the solar reflective coat compensates the penalty effect on winter; this means that the reduction of thermal conduction through the roof in summer, due to the lower absorption of solar radiation, compensates, in terms of energy consumption to obtain indoor comfort, the increase in energy consumption necessary to achieve indoor comfort conditions in winter, due to the aforementioned lower absorption of solar radiation. This effect is clear in Figure 10, Figure 11, Figure 12 and Figure 13, where it can be observed the strong reduction in cooling loads for the roof with cool coat without insulation when compared to the reference case, while, regarding the heating loads, although higher for the reference roof, it has a lower difference on thermal loads when compared to the cool roof; this difference is collected by the total load, which, as said before, is lower for the roof with cool coat.
On the other hand, noteworthy is the energy performance of the roofs with insulation. These roofs are able of reducing thermal conduction in both directions, from inside to outside and vice versa, which produces a significant reduction in thermal loads compared to roofs without thermal insulation, as can been observed in Figure 8 and Figure 9. It is observed, in Figure 10, Figure 11, Figure 12 and Figure 13, that the higher the thickness of the insulation layer, the less the thermal loads, for both cooling and heating, which is in accordance with the reduction of heat transfer through the roof caused by the insulation layer. On the other hand, as expected, it can be observed the progressive reduction of the cooling loads as the solar absorptivity increases, obviously due to the decreasing solar radiation absorbed by the roof, whilst, for the same reason, the opposite is observed for heat loads. However, what is concluded, as it is shown in Figure 8 and Figure 9, is that, for the climatic conditions and the geographical framework in which the study has been developed, the total load decreases as the reflectivity increases due to the fact that the reduction in thermal flux to the interior achieved in summer, and which impacts on cooling loads, is greater than the increase in thermal flux from the inside to the outside, which occurs in winter, and that is reflected in an increase in heating loads.

3.2. Energy Optimization Results

This section studies the optimum combination of solar reflectivity and insulation thickness in terms of thermal loads for the entire LC time interval. Energy savings of the different combinations when compared to the reference case are also reported.
From the results shown in Figure 8 and Figure 9, it can be concluded that the optimum combination to reduce energy consumption is absorptivity equals 0.1, and insulation thickness equals 0.1. This result holds for all the maintenance scenarios, although the results for the different maintenance scenarios are very close to each other. To clarify this, in Table 6 and Table 7, the values of the total thermal loads for the LC timespan are shown.
The similarity of the loads among the different maintenance protocols for the cool roof can be observed in all those cases for which the absorptivity values and the thickness of the insulation layer are the same. Note that, for the optimum combination established, the total thermal load difference between the optimum maintenance scenario, i.e., triennial washing, and the worst scenario, i.e., no washing, is only of 0.01676 GJ/m 2 , equivalent to a 1.55 % difference.
In Table 8 and Table 9, the savings for the total thermal loads of every retrofitted roof when compared to the reference roof are shown. In accordance with the energy results presented above, all the combinations are found to provide energy savings compared to the reference case.
From the above, it is straightforward that the greatest energy savings are obtained by the retrofitting combination that uses a cool roof with an absorptivity equal to 0.1 and an insulation layer thickness equal to 0.1 m, with the triennial washing of the cool coat being the case that gives the highest energy saving; on the other hand, the energy savings decrease as the absorptivity increases, and the insulation thickness decreases, according to the previous discussion on thermal loads.

4. Cost Optimization Analysis

The objective of this analysis is to assess the cost-effectiveness of the energy retrofitting measures proposed when compared to the reference case, as well as to find out the combination of cool coating reflectivity and insulation layer thickness that yields the minimum LC cost under the climate conditions considered.
The optimum combination of cool coating reflectivity and insulation thickness is a function of the costs of the cool paint and the insulation material, the application of the cool paint, the insulation layer installation, the cool coat maintenance, the thermal loads due to the energy transfer through the roof, the coefficient of performance (COP) of the heating and cooling systems, the retrofitting materials lifetime, the energy prices, the energy inflation rate, and the monetary discount rate.
Despite the aforementioned tendency of heating and cooling loads to decrease when the thickness of the insulation layer increases, a rise in the reflectivity of the cool coat tends to reduce the cooling load, as well as the solar gain, in winter. Furthermore, the aging and the washing of the cool coat give rise on their own to contrary effects on the loads. So, this time-varying evolution of the thermal loads, combined with the effect of the costs of the retrofitting material and its maintenance, makes it necessary to carry out a detailed analysis of the combined effect on the LCA cost of the different retrofitting combinations in order to obtain the optimum results. The main difference between the present analysis and the usual cost analysis is that, here, the costs are computed monthly for the whole LC time interval, with the goal of assessing the effect on the evolution of costs, along time, produced by changes in the thermal loads associated to the cool coat reflectivity variation, throughout the LC timespan.

4.1. Methodology

First, the assessment of the cost-effectiveness of every ESRM is done. To that aim, the methodology from Reference [61] is used. This methodology performs an LCCA, where all anticipated costs are estimated and discounted to their present worth.
The procedure consists in computing the costs for each future time period using the future value of the variables involved and then discounting each future cost to its present worth. Thereby, the sum of all the present values of the costs yields the life-cycle cost.
Once the LC cost has been determined for each of the considered ESRMs, including different scenarios for maintenance, the cost-effectiveness of each case is assessed by its LC cost: the lower this cost, the more cost effective the ESRM, and the one having the lowest LC cost is selected as the optimum.
The lifetime for the life-cycle analysis is taken as equal to n = 20 years, which is the mean service life of the used cool coat and the insulation material [62]. Then, for every considered ESRM, the life-cycle cost, the associated savings, and the payback period are computed.
In the present LCCA, the method P1-P2 introduced in Reference [61] is used with some little changes to account for the time changing character of the energy operating costs associated to the variation of the cool coat solar reflectivity due to the aging and washing of the coat. We use the expression
P W k = 1 / ( 1 + d ) k
to compute the present worth ( P W ) of one monetary unit belonging to the future time interval k (usually years), where d is the market discount rate (fraction per time period).
In the present analysis, k, for k = 1 , , n , is taken as the year in which the costs are being computed. Then, if we note as Q c ( k ) and Q h ( k ) the cooling and heating loads at year k for each ESRM, the energy payment for the whole year k is given by:
C e ( k ) = Q c ( k ) × P e , c C O P c × ( 3.6 × 10 6 ) + Q h ( k ) × P e , h C O P h × ( 3.6 × 10 6 ) ( 1 + i ) k 1 ,
where P e , c and P e , h are the current prices of energy used to cool and heat for cooling and heating, respectively, C O P c and C O P h are the coefficients of performance for cooling and heating, and, finally, i is the considered inflation rate for energy costs.
Then, in accordance with expression (8), the present worth of the future energy payment C e ( k ) at the period k is given by
P W k ( C e ( k ) ) = C e ( k ) / ( 1 + d ) k .
It is noteworthy that, although C e ( k ) represents the energy costs for year k, it is computed by taking into account the changing values of energy consumption for every specific year k according to the monthly changes in the solar reflectivity of the cool coat.
To compute C e ( k ) , the energy source considered for cooling is electricity, while, for heating, the use of electricity or gas have been considered, given that they are the usual sources of energy used in the geographical framework under analysis.
To include the maintenance in the cost computation, we name as i M the inflation rate for this kind of costs. Then, again applying (8), the present worth of the future maintenance payment C M ( k ) at the period k is given by
P W k ( C M ( k ) ) = C M ( k ) / ( 1 + d ) k = δ ( k ) C M ( 1 + i M ) k 1 / ( 1 + d ) k ,
where C M is the current price of maintenance per unit area of roof surface, and δ ( k ) is defined as:
δ ( k ) = 1 , i f   t h e   m a i n t e n a n c e   i s   d o n e   a t   p e r i o d   k ; 0 , i f   t h e   m a i n t e n a n c e   i s   n o t   d o n e   a t   p e r i o d   k .
Then, the present worth of the energy total cost for the life-cycle timespan is computed as C e = k = 1 n P W k ( C e ( k ) ) , whilst the present worth value of maintenance life-cycle total cost is computed as C M = k = 1 n P W k ( C M ( k ) ) . Now, if we note as C I , c c and C I , i n s the initial costs of installation of the cool coat and the insulation layer, the total initial cost for installation is given by C I = C I , c c + C I , i n s . Finally, the present worth of the life-cycle total cost per unit area for each considered ESRM and maintenance scenario is given by:
C t = C e + C M + C I .
Once C t is computed for every case, the optimum case is considered to be the one that yields the lower value for C t .
On the other hand, to compute the possible savings of the proposed ESRMs, the P W of the LC total cost must be estimated for the reference case, too. Then, the P W of the energy costs for the reference case at period k, P W k ( C e ( r e f ) ( k ) ) , is computed by:
P W k ( C e ( r e f ) ( k ) ) = C e ( r e f ) ( k ) / ( 1 + d ) k ,
where C e ( r e f ) ( k ) is the energy future cost at period k for the reference roof calculated by:
C e ( r e f ) ( k ) = Q c ( r e f ) ( k ) × P e , c C O P c × ( 3.6 × 10 6 ) + Q h ( r e f ) ( k ) × P e , h C O P h × ( 3.6 × 10 6 ) ( 1 + i ) k 1 ,
with Q c ( r e f ) ( k ) and Q h ( r e f ) ( k ) being the cooling and heating loads, respectively, for the reference roofs at the time period k.
Now, we note by C I ( r e f ) and C M ( r e f ) the installations and maintenance costs, respectively, of the minimum retrofitting action to be done on the reference roof and which is deemed to be the bituminous paint retrofitting.
Then, the difference between the costs for the cool roofs and the reference case for the whole LC period is estimated through the net savings (NS):
N S = k = 1 n P W k ( C e ( r e f ) ( k ) ) P W k ( C e ( k ) + ( C M ( r e f ) C M ) + ( C I ( r e f ) C I ) .
A positive value for N S means that the use of the ESRM under consideration yields economic savings when compared to the reference case and that the retrofitting system is cost-efficient. Due to the outlined thermal loads change over time, the net savings (11) cannot be reduced to a simple analytical expression, and its computation must be performed through the yearly computation, according to the yearly changes in thermal loads; this represents a distinguishing feature of our analysis with respect to the usual analysis of cost-effectiveness for retrofitting measures found in the literature.
To compute the period of investment return or payback period (PB), we follow Reference [63] and take PB as the time horizon t, for which the value of NS becomes zero. To compute this value of t, since the terms of the series in Equation (11) are not constant on k, the net savings accumulated for every time horizon k 0 is computed first:
N S ( k 0 ) = k = 1 k 0 P W k ( C e ( r e f ) ( k ) ) P W k ( C e ( k ) C M C I .
Then, if N S ( k 0 ) < 0 and N S ( k 0 + 1 ) > 0 , the value of the payback period t is estimated as:
t = N S ( k 0 ) N S ( k 0 + 1 ) N S ( k 0 ) + k 0 .
The value of t given by (13) is consistent with the value provided by the usual formulas for the PB found in the literature used under the consideration of constant energy consumption for all the years of the LC period [31,61]. In Appendix C, a comparison between the formula given by (13) and the formula for the computation of the PB from References [31,61,63] is performed.

4.2. Variables for the LCA Economic Analysis

The input variables needed to carry out the economic calculations involved in the LCCA are shown in Table 10.
The costs of the cool and bituminous paint application, their maintenance, and the insulation material and its installation are taken from Reference [62]. The electricity and natural gas prices are the average final prices for household consumers in Spain during the year 2020 as reported by Eurostat Office [64]. The values of the energy inflation and the discount rates have been taken considering the average values of these parameters under Spanish and European economic conditions, with such values being previously used for the LCA cost analysis in former research [37,65]. Regarding the conditioning equipment, it has been considered that cooling is always done by air conditioning machinery using electricity as the only energy source, while, for heating, we have considered the following cases:
  • Heating by air (heat pumps).
  • Heating by natural gas.
  • Heating by electricity radiators.
The machinery is considered common in the local market and it has a mid-range efficiency. The efficiency indices of the conditioning devices are listed in Table 11.

4.3. Cost Optimization Results

In this section, the economic results from the LCCA for the considered ESRMs are shown. The objective of this analysis is twofold: first, the cost-effectiveness of the ESRMs will be assessed; second, the optimum case for the combination of insulation layer, roof solar absorptivity, and maintenance protocol will be determined.
The NS and PB were calculated by using the economic variables introduced in Section 4.1 and considering the aging and maintenance scenarios presented in Section 2.2. For each ESRM, a positive NS means real savings and, conversely, a negative NS implies that the ESRM has a higher LCA cost than the reference case and, therefore, no real savings. Finally, once C t is computed for every considered ESRM and for every cool roof maintenance protocol, the lowest value of C t determines the optimum combination.
On the other hand, the cost of the roof retrofitting must be paid by all the owners of the building, while the benefit of the retrofitting would only be noticeable in the housing under the roof; however, in the event of damage, according Spanish law, the cost of restoring the roof is paid by all the owners of the building; so, since the owners and tenants of the dwellings taken as the case study are characterized by having low incomes, a retrofitting action on the roof is unlikely, except for the existence of some deterioration of the roof that makes such retrofitting necessary; thus, to perform the analysis, it is assumed that the reference roof needs some kind of renovation or maintenance. This way, the application of a bituminous paint is chosen since this is considered to be the most usual retrofitting measure in the socioeconomic context within which the study is performed.
In Table 12, the optimum values with regard to the use of pump for heating and cooling are shown for every protocol maintenance and for the whole LC timespan. These values have been extracted from Table A4 in Appendix B, where all the LC total costs are shown for the different values of both the insulation thickness layer and the roof reflectivity, as well as for all the colorredconsidered roof maintenance protocols. For this case of considering electricity as the only energy source to conditioning, the estimated LC total cost for the reference case is equal to 174.21/m 2 . Then, this value is used to obtain the economic savings and the payback periods shown in the table.
As can be observed in Table 12, when a pump is used for heating and cooling, the cool roof absorptivity that provides the minimum total cost is equal to 0.1, and the insulation thickness is equal to 0.08 for all the maintenance scenarios. The optimum result is obtained for the non-maintenance scenario with a total cost equal to 64.27/m 2 , a total savings equal to 109.94/m 2 , and a payback period of 4.54 years.
In Table 13, the optimum values regarding the use of gas for heating and of a pump for cooling are shown for every protocol maintenance and for the whole LC timespan. These values have been sourced from Table A5 in Appendix B, where all the LC total costs are shown for the different values of both the insulation thickness layer and the roof reflectivity, as well as for all the roof maintenance protocols. Within the framework of this energy source for conditioning, the LC total cost for the reference case is equal to 180.07/m 2 . For the cases of no-maintenance and decennial maintenance, the combination that provides the maximum savings is that which has a solar absorptivity equal to 0.1 and an insulation thickness equal to 0.08, while, for the quinquennial maintenance, the combination yielding the maximum savings is that which has a solar absorptivity equal to 0.1 and an insulation thickness equal to 0.09 and, finally, for the triennial maintenance, the combination that provides the maximum savings is that which has a solar absorptivity equal to 0.2 and an insulation thickness equal to 0.09. For this scenario of mixed energy used to condition the building, the optimum result is also obtained for the non-maintenance scenario with a total cost equal to 68.27/m 2 , a total savings equal to 111.80/m 2 , and a payback period of 4.05 years.
In Figure 14, Figure 15, Figure 16 and Figure 17, the insulation cost, the cool coating cost (including maintenance), the energy consumption cost, and the total cost are displayed with respect to the insulation thickness for the optimum absorptivities shown in Table 12 and Table 13, for all the maintenance scenarios and for the two considered sources of energy.
As it can be seen from these figures, there is a rapid decrease in energy costs when the insulation thickness begins to increase from zero. This also results in a rapid decrease in total costs. However, for values of insulation thickness greater than 0.07 m, energy costs decrease at a much weaker rate according to the decrease in the rate of reduction of thermal loads. Thereby, the decrease in energy costs is not balanced by the increase in insulation costs, with this fact being much more evident for the higher frequency washing of the cool coating. This way, when using electricity as the only energy source for conditioning, the curve of total costs has a minimum for the insulation thickness equal to 0.08 m that corresponds to the optimum insulation thickness for all the cool roof maintenance cases; when using gas for heating and electricity for cooling, in the case of quinquennial cool roof maintenance, the optimum is reached for an insulation thickness equal to 0.09 m, whilst, for the cases of no-maintenance, decennial, and triennial cool roof maintenance, the curve of total costs exhibits a minimum for an insulation thickness equal to 0.08 m, even though, for the triennial maintenance, the minimum total cost is reached for a solar absorptivity equal to 0.2, as is shown in Table A5 in Appendix B.
As it can be observed in Table A4 and Table A5 in Appendix B, even if, as it is reported in the energy analysis, the minimum energy consumption was obtained for triennial washing when absorptivity equals 0.1 and insulation thickness equals 0.1 m, energy cost savings are not enough to compensate the rise of cost derived of the increase in insulation thickness and of the cool roof maintenance. On the other hand, as it is deduced from those tables, all the studied ESRMs are cost-effective and produce positive savings. This result holds for all the different cases of maintenance of cool roofs and of machinery for conditioning.
Thus, from the performed economic analysis, it can be concluded that the optimum combination is that obtained for a cool coat solar absorptivity equal to 0.1, an insulation thickness equal to 0.08 m, no washing of the cool coat, and when using a pump for heating and cooling.
Finally, it is noteworthy to observe the limited scope of the cool roof maintenance through power washing as it can be deduced from the energy and economic analysis. This fact has been discussed formerly in the literature; this way, in Reference [50], it is pointed out that the differences induced by the washing of the cool roofs are very small for the LCA total period, which calls into question its suitability, especially when considering the possible damage to the cool coating during maintenance.

5. Conclusions

The thermal performance of social building roofs under Southern Spain climate and its refurbishment by means of combining a cool roof coat and a thermal insulation layer have been analyzed in the present paper. The thermal dynamic of the roofs was computed through a numerical model using a finite difference method in order to assess the impact of different combinations of cool coating reflectivities and insulation thickness when used to retrofit the roofs of social buildings built in the sixties of the last century at the city of Seville, Spain. The time-dependent impact of the aging effect on the roof coat reflectivity was also incorporated into the model.
Combinations of cool coating, for which values of solar absorptivity were equal to 0.1, 0.2, 0.3, 0.4, and 0.5, with insulation layer thickness ranging from 0 to 0.1 m, were analyzed, and the results obtained were compared to those reported for the initial reference roof belonging to the social housing. In order to perform the analysis under realistic conditions, the aging effect of the cool coatings was taken into account by considering a pattern for the aging effect and its impact on the solar reflectivity of the cool coat. This aging effect was estimated a monthly basis, and it was introduced in energy computations for the whole service lifetime of the refurbished roof to perform an LCCA in order to appraise the cost-effectiveness of the proposed combinations.
Then, a comprehensive study was conducted by taking into consideration different values for initial roof reflectivity and insulation thickness, different sources of energy for the conditioning of the building, and different maintenance patterns to restore the cool coating reflectivity.
The results obtained for the different combinations here studied point to significant savings in the operational energy. Such savings are larger as the absorptivity of the external coating of the retrofitted roof decreases and the insulation layer thickness increases. This way, a decrease in the annual total loads by 84.56 % is found for the best case that is given by insulation thickness equal to 0.1, initial roof absorptivity equal to 0.1, and triennial washing.
Likewise, it has been concluded the limited impact of the cool roofs maintenance on the total loads during the LCCA period of time. This way, the values of total loads are similar for the considered maintenance scenarios, with the maximum difference of around 1.55 % for each combination being considered. However, regarding the heating and cooling loads separately, some variations in the evolution of the loads can be observed for every maintenance scenario.
In order to evaluate the economical value of the different combinations studied for the refurbishment of the roofs belonging to the social housing considered, a 20-year LCCA was conducted. The LCCA pointed out that all the combinations of insulation layer thickness, roof solar reflectivity, and maintenance scenarios are cost-efficient. Again, the low impact of the maintenance on the total costs for the whole LCCA timespan is clear. However, the curves of costs result are affected by installation costs, so that the combination that produces the lowest operational cost is not the one with the lowest total cost for the entire period of time analyzed. Thus, it can be concluded from the LCCA performed that the optimum combination is that obtained for a cool coat solar absorptivity equal to 0.1, an insulation thickness equal to 0.08 m, no washing of the cool coat, and when using a pump for heating and cooling. Such a combination allows a total savings equal to 109.94/m 2 , equivalent to the 63.10 % of the reference costs, as well as a payback period of 4.54 years.
According to the results obtained in the present research, the use of an adequate combination of insulation and cool roof for the refurbishment of energy obsolete roofs is profitable in terms of energy consumption, economic costs, and environmental benefits. Because of all this, its use is recommended for the refurbishment of roofs belonging to the social housing in the climatic zone considered.

Author Contributions

Conceptualization, C.-A.D.-T., H.D.-T. and A.D.-D.; methodology, all the authors; software, all the authors; formal analysis, all the authors; investigation, all the authors; writing—original draft preparation, H.D.-T. and C.-A.D.-T.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Economy and Competitiveness of the Spanish Government and the European Regional Development Fund through the research and development project “Parametric Optimization of Double Skin Facades in the Mediterranean Climate to Improve Energy Efficiency Under Climate Change Scenarios” (ref BIA2017-86383-R).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Thermal Loads Results

Table A1. Total loads (GJ/m 2 ) for the LC timespan.
Table A1. Total loads (GJ/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.05.0352.88322.62432.38792.17333.19182.91212.65052.41172.1951
0.031.52071.38031.25111.13191.02251.53591.39411.26361.14321.0327
0.051.14691.04130.94430.83570.74801.15841.05170.95370.84410.7555
0.070.91900.83530.75830.68620.62210.92820.84360.76590.69310.6283
0.10.71830.66100.59680.54370.48900.72550.66760.60280.54910.4939
Quinquennial WashingTriennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.03.26922.97032.69522.44792.16193.31823.02742.75552.50732.2820
0.031.58411.43781.30321.17911.06511.59681.44931.31361.18851.0736
0.051.19471.08470.98360.87050.77921.20421.09340.99150.87750.7854
0.070.95730.87010.78990.71480.64800.96490.87710.79620.72050.6532
0.10.74820.68850.62170.56640.50940.75420.69410.62660.57090.5135
Table A2. Heating loads (GJ/m 2 ) for the LC timespan.
Table A2. Heating loads (GJ/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.03.16022.88322.62432.38792.17333.19182.91212.65052.41172.1951
0.031.52071.38031.25111.13191.02251.53591.39411.26361.14321.0327
0.051.14691.04130.94430.83570.74801.15841.05170.95370.84410.7555
0.070.91900.83530.75830.68620.62210.92820.84360.76590.69310.6283
0.10.71830.66100.59680.54370.48900.72550.66760.60280.54910.4939
Quinquennial WashingTriennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.03.26922.97032.69522.44792.16193.31823.02742.75552.50732.2820
0.031.58411.43781.30321.17911.06511.59681.44931.31361.18851.0736
0.051.19471.08470.98360.87050.77921.20421.09340.99150.87750.7854
0.070.95730.87010.78990.71480.64800.96490.87710.79620.72050.6532
0.10.74820.68850.62170.56640.50940.75420.69410.62660.57090.5135
Table A3. Cooling loads (GJ/m 2 ) for the LC timespan.
Table A3. Cooling loads (GJ/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.01.87482.30572.76393.24573.74821.76262.16772.59863.05153.5240
0.030.75040.93531.12911.33051.53940.80980.99221.18411.38951.5783
0.050.56620.70450.84971.00691.19340.58860.69630.83921.01111.1608
0.070.46450.57680.69440.81500.94500.44670.56300.68010.80540.9304
0.10.36180.46970.53610.63080.72970.34450.43340.52540.62380.7196
Quinquennial WashingTriennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.01.69852.08892.50412.94053.39581.60241.97072.36232.77413.2036
0.030.80430.99101.18181.38541.57640.75150.94741.13821.33691.5448
0.050.58550.68080.81920.98331.13700.56370.65220.78980.94691.1105
0.070.42620.54340.66080.78680.91380.40130.51720.63570.76000.8935
0.10.32800.42170.51330.60960.70740.30920.40300.49530.59190.6921

Appendix B. Life Cycle Total Costs

Table A4. Cooling and heating by pump. Total costs (/m 2 ) for the LC timespan.
Table A4. Cooling and heating by pump. Total costs (/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.0129.32131.66135.21139.77145.31132.30134.28137.93142.46146.94
0.01104.76106.01109.27112.66125.55104.83106.54109.11113.82126.63
0.0291.3592.5595.0298.21105.2393.9395.9198.07103.17106.63
0.0383.5184.5686.3088.3290.9085.4086.4988.2190.5792.53
0.0477.2078.0479.6980.9783.6978.8379.9881.0583.9284.63
0.0571.7472.5873.8074.8777.0173.3973.8375.0376.4978.64
0.0667.1468.9470.1171.5574.8269.7770.0871.7972.2675.63
0.0764.6166.0167.4168.9670.4167.0767.9268.6369.9372.03
0.0864.2765.4066.4967.7168.8266.5167.1067.7168.7469.63
0.0964.4764.7865.1665.7366.0166.6767.2167.5367.6067.63
0.164.8164.8664.9565.0365.1167.0867.4167.6467.7167.74
Quinquennial WashingTriennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.0135.76138.16141.76146.01150.50140.60143.02146.13149.32153.54
0.01106.21109.31112.66116.03120.76109.76112.83116.21119.96124.15
0.0295.3896.8598.22100.92103.7698.82100.23101.97104.08106.49
0.0388.2089.7790.1892.3794.8892.1693.5595.1997.2898.68
0.0480.8682.1083.5684.9886.5285.1585.2086.6288.2189.91
0.0576.1576.6977.8179.0280.3179.9680.6481.7582.7384.07
0.0671.7472.3974.1075.2276.8175.7976.6077.5978.3979.11
0.0769.3770.1071.2272.4673.7473.7474.4375.3776.2477.04
0.0868.3368.3968.4668.7868.9972.8273.1973.5873.9074.38
0.0968.4068.5568.6168.8369.1172.8973.2673.8174.2174.61
0.168.7268.7168.8969.3869.5773.1173.3974.0774.4874.71
Table A5. Cooling by pump and heating by gas. Total costs (/m 2 ) for the LC timespan.
Table A5. Cooling by pump and heating by gas. Total costs (/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.0141.20142.47145.06148.73153.473143.75144.84147.56151.26155.93
0.01108.10109.51112.27112.66125.55111.63112.84114.27117.02122.63
0.0295.7096.6599.8798.21105.2398.6399.31100.87102.63104.11
0.0388.2288.9190.9992.5894.7490.9391.5192.8194.7295.68
0.0480.6181.3482.0980.9783.6983.6383.9485.3686.0286.73
0.0576.0476.4877.34678.2279.177.5877.7778.6079.7880.48
0.0671.1471.5472.6171.5574.8273.6973.9874.6975.8676.33
0.0769.1969.6170.3571.3672.7470.5770.9271.4972.5573.78
0.0868.2768.7069.1970.9468.8269.8170.0770.6171.5471.63
0.0968.5468.8869.4670.4366.0170.0870.1370.8371.2071.55
0.169.2769.6370.1570.5965.1170.3470.4970.9471.3371.49
0.0147.44147.90150.42153.99158.57152.81151.93154.03157.35161.78
0.01114.09115.78118.66121.53124.34118.95118.23120.57122.17126.03
0.02101.76102.99104.22106.01108.16106.12107.13108.55110.28111.01
0.0393.7994.3195.5396.7397.8197.9798.2699.32100.76101.8
0.0485.6586.1087.1688.6189.9289.9590.0391.1192.1693.88
0.0580.3180.3681.1282.2983.3684.2284.2584.9085.8887.15
0.0675.5575.7876.3077.1277.9879.7879.7080.4881.3982.39
0.0773.0573.6174.1374.8975.5377.6277.2577.6278.579.73
0.0871.6372.1172.6673.1874.8075.5274.9975.3875.9676.98
0.0971.371.8572.7173.0374.0675.5375.0675.5175.7176.41
0.171.4572.0772.9373.1273.9775.6875.2775.8676.0176.28

Appendix C. Comparison of Formulas to Compute the Payback Period

In this appendix, the formula given by Equation (13), used here to compute the payback period, is compared to the one used in former works as in References [31,63], where a constant value of the yearly energy consumption for the whole LCCA timespan is considered. The formula used in these reference to calculate the payback period ( P B ), when applied to the framework of the present work, is given by:
P B = L n 1 ( d i ) C I , i n s + P 2 C I , c c A s L n 1 + i 1 + d i d P B = ( 1 + i ) C I , i n s + P 2 C I , c c A s i = d ,
where A s is the annual energy saving, taken constant, and P 2 = 1 + P 1 M s , with M s being a term that accounts for the present worth value of the cool coat maintenance, and P 1 is defined as
P 1 = u = 1 n ( 1 + i ) ( u 1 ) ( 1 + d ) u .
More details about these formulas can been found in the cited works.
To perform the comparison between both formulas, the optimum combination of insulation layer and cool coat reflectivity found in the developed analysis, is considered. Because the formula given by Equation (A1) is only valid for constant values of yearly energy consumption, to compute A s , the value of the energy consumption for the whole LCCA timespan of the optimal combination is divided by n, the number of years considered for the LCCA. Then, both formulas are applied for these values of energy consumption and for three different cases of i and d, the one used here, and two more combinations, in order to have some extra comparison values. The results of the comparison are shown in Table A6. In this table, the PB formula that uses constant values of yearly energy consumption is named Constant, and the presented here is named Variable.
Table A6. Comparison of formulas to compute the payback period.
Table A6. Comparison of formulas to compute the payback period.
PB Formulabbb i = 2 % , d = 1 . % bb i = 3 % , d = 1.5 % bb i = 1.5 % , d = 3 %
Constantbbb4.1834   [years]bb4.1721   [years]bb4.4442   [years]
Variablebbb4.1827   [years]bb4.1710   [years]bb4.4460   [years]
As it is observed in Table A6, the values provided by both formulas are very close. It is noteworthy to observe that the PB obtained for the optimal combination, when yearly energy consumption is considered constant, differs of the value of PB obtained when considering the energy consumption variable due to the aging effect, as in Table 12, which supports the stated need to take this effect into account when conducting the LCCA and the optimization analysis.

References

  1. DB-HE2019. Basic Energy Saving Document. Spanish Royal Decrete 732/2019. Available online: https://www.codigotecnico.org/DocumentosCTE/AhorroEnergia.html (accessed on 15 June 2021).
  2. European Commission. DG Energy. Available online: https://ec.europa.eu/clima/policies/strategies/2030-en (accessed on 2 July 2021).
  3. Yu, J.; Tian, L.; Yang, C.; Xu, X.; Wang, J. Optimum insulation thickness of residential roof with respect to solareair degree-hours in hot summer and cold winter zone of China. Energy Build. 2011, 43, 2304–2313. [Google Scholar] [CrossRef]
  4. Gentle, A.R.; Aguilar, J.L.C.; Smith, G.B. Optimized cool roofs: Integrating albedo and thermal emittance with R-value. Sol. Energy Mater. Sol. Cells 2011, 95, 3207–3215. [Google Scholar] [CrossRef] [Green Version]
  5. Ramamurthy, P.; Sun, T.; Rule, K.; Bou-Zeid, E. The joint influence of albedo and insulation on roof performance: A modeling study. Energy Build. 2015, 102, 317–327. [Google Scholar] [CrossRef] [Green Version]
  6. Arumugam, R.S.; Garg, V.; Ram, V.V.; Bhatia, A. Optimizing roof insulation for roofs with high albedo coating and radiant barriers in India. J. Build. Eng. 2015, 2, 52–58. [Google Scholar] [CrossRef] [Green Version]
  7. Piselli, C.; Saffari, M.; Gracia, A.; Pisello, A.L.; Cotana, F.; Cabeza, L.F. Optimization of roof solar reflectance under different climateconditions, occupancy, building configuration and energy systems. Energy Build. 2017, 151, 81–97. [Google Scholar] [CrossRef]
  8. Shi, Z.; Zhang, X. Analyzing the effect of the longwave emissivity and solar reflectance of building envelopes on energy-saving in buildings in various climates. Sol. Energy 2011, 85, 28–37. [Google Scholar] [CrossRef]
  9. Akbari, H.; Levinson, R.; Rainer, L. Monitoring the energy-use effects of cool roofs on California commercial buildings. Energy Build. 2005, 37, 1007–1016. [Google Scholar] [CrossRef] [Green Version]
  10. Romeo, C.; Zinzi, M. Impact of a cool roof application on the energy and comfort performance in an existing non-residential building. A Sicilian case study. Energy Build. 2013, 67, 647–657. [Google Scholar] [CrossRef]
  11. Akbari, H.; Bretz, S.; Kurn, D.M.; Hanford, J. Peak power and cooling energy savings of high-albedo roofs. Energy Build. 1997, 25, 117–126. [Google Scholar] [CrossRef]
  12. Akbari, H.; Levinson, R. Evolution of cool-Roof standards in the US. Adv. Build. Energy Res. 2008, 2, 1–32. [Google Scholar] [CrossRef]
  13. Akbari, H.; Konopacki, S. Calculating energy-saving potentials of heat-island reduction strategies. Energy Policy 2005, 33, 721–756. [Google Scholar] [CrossRef]
  14. Akbari, H.; Matthews, H.D. Global cooling updates: Reflective roofs and pavements. Energy Build. 2012, 55, 2–6. [Google Scholar] [CrossRef]
  15. Xu, T.; Sathaye, J.; Akbari, H.; Garg, V.; Tetali, S. Quantifying the direct benefits of cool roofs in an urban setting: Reduced cooling energy use and lowered greenhouse gas emissions. Build. Environ. 2012, 48, 1–6. [Google Scholar] [CrossRef]
  16. Bozonnet, E.; Doya, M.; Allard, F. Cool roofs impact on building thermal response: A French case study. Energy Build. 2011, 43, 3006–3012. [Google Scholar] [CrossRef]
  17. Santamouris, M. Cooling the cities—A review of reflective and green roof mitigation technologies to fight heat island and improve comfort in urban environments. Sol. Energy 2014, 103, 682–703. [Google Scholar] [CrossRef]
  18. Seifhashem, M.; Capra, B.R.; Milller, W.; Bell, J. The potential for cool roofs to improve the energy efficiency of single storey warehouse-type retail buildings in Australia: A simulation case study. Energy Build. 2018, 158, 1393–1403. [Google Scholar] [CrossRef]
  19. Boixo, S.; Diaz-Vicente, M.; Colmenar, A.; Castro, M.A. Potential energy savings from cool roofs in Spain and Andalusia. Energy 2012, 38, 425–438. [Google Scholar] [CrossRef]
  20. Pisello, A.L.; Santamouris, M.; Cotana, F. Active cool roof effect: Impact of cool roofs on cooling system efficiency. Adv. Build. Energy Res. 2013, 7, 209–221. [Google Scholar] [CrossRef]
  21. Kolokotroni, M.; Gowreesunker, B.L.; Giridharan, R. Cool roof technology in London: An experimental and modelling study. Energy Build. 2013, 67, 658–667. [Google Scholar] [CrossRef] [Green Version]
  22. De Masi, R.F.; Ruggiero, S.; Vanoli, G.P. Acrylic white paint of industrial sector for cool roofing application: Experimental investigation of summer behavior and aging problem under Mediterranean climate. Sol. Energy 2018, 169, 468–487. [Google Scholar] [CrossRef]
  23. Kibert, C.J. Sustainable Construction: Green Building Design and Delivery; John Wiley and Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  24. Vailati, M.; Monti, G. Earthquake-Resistant and Thermo-Insulating Infill Panel with Recycled-Plastic Joints. In Earthquakes and Their Impact on Society. Springer Natural Hazards; DÁmico, S., Ed.; Springer: Cham, Switzerland, 2016. [Google Scholar] [CrossRef]
  25. Manfredi, V.; Masi, A. Seismic Strengthening and Energy Efficiency: Towards an Integrated Approach for the Rehabilitation of Existing RC Buildings. Buildings 2018, 8, 36. [Google Scholar] [CrossRef] [Green Version]
  26. Vailati, M.; Monti, G.; Di Gangi, G. Earthquake-Safe and Energy-Efficient Infill Panels for Modern Buildings. In International Conference on Earthquake Engineering and Structural Dynamics; Rupakhety, R., Olafsson, S., Eds.; Springer: Cham, Switzerland, 2018; Volume 44. [Google Scholar] [CrossRef]
  27. Zingre, K.T.; Kiran Kumar, D.E.V.S.; Wan, M.P.; Chao, C.Y.H. Effective R-value approach to comprehend the essence of integrated opaque passive substrate properties. J. Build. Eng. 2021, 44, 102865. [Google Scholar] [CrossRef]
  28. Hernández-Pérez, I. Influence of Traditional and Solar Reflective Coatings on the Heat Transfer of Building Roofs in Mexico. Appl. Sci. 2021, 11, 3263. [Google Scholar] [CrossRef]
  29. De Masi, R.F.; Gigante, A.; Ruggiero, S.; Vanoli, G.P.; Ruggiero, S. Impact of weather data and climate change projections in the refurbishment design of residential buildings in cooling dominated climate. Appl. Energy 2021, 303, 117584. [Google Scholar] [CrossRef]
  30. Levinson, R.; Akbari, H. Potential benefits of cool roofs on commercialbuildings: Conserving energy, saving money, and reducing emission ofgreenhouse gases and air pollutants. Energy Effic. 2010, 3, 53–109. [Google Scholar] [CrossRef] [Green Version]
  31. Saafi, K.; Daouas, N. A life-cycle cost analysis for an optimum combination of cool coating and thermal insulation of residential building roofs in Tunisia. Energy 2018, 152, 925–938. [Google Scholar] [CrossRef]
  32. Jo, J.H.; Carlson, J.D.; Golden, J.S.; Bryan, H. An integrated empirical and modeling methodology for analyzing solar reflective roof technologies on commercial buildings. Build. Environ. 2010, 45, 453–460. [Google Scholar] [CrossRef]
  33. Sproul, J.; Wan, M.P.; Mandel, B.H.; Rosenfeld, A.H. Economic comparison of white, green, and black flat roofs in the United States. Energy Build. 2013, 71, 20–27. [Google Scholar] [CrossRef]
  34. Zhang, Z.; Tong, S.; Yu, H. Life Cycle Analysis of Cool Roof in Tropical Areas. Procedia Eng. 2016, 169, 392–399. [Google Scholar] [CrossRef]
  35. Yuan, J.; Farnham, C.; Emura, K.; Alam, M.A. Proposal for optimum combination of reflectivity and insulation thickness of building exterior walls for annual thermal load in Japan. Build. Environ. 2016, 103, 228–237. [Google Scholar] [CrossRef]
  36. Spanish Royal Decrete 2429/1979. Approving the Basic Building on Thermal Conditions in Buildings NBE-CT-79; BOE Nro. 253; Norma Básoca de la Edificación: Barcelona, Spain, 1979; pp. 24524–24550. [Google Scholar]
  37. Domínguez-Delgado, A.; Domínguez-Torres, C.A.; Domínguez-Torres, H. Energy and Economic Life Cycle Assessment of Cool Roofs Applied to the Refurbishment of Social Housing in Southern Spain. Sustainability 2020, 12, 5602. [Google Scholar] [CrossRef]
  38. Hickson, R.I.; Barry, S.I.; Mercer, G.N.; Sidhu, H.S. Finite difference schemes for multilayer diffusion. Math. Comput. Model. 2011, 54, 210–220. [Google Scholar] [CrossRef]
  39. Hagishima, A.; Tanimoto, J. Field measurements for estimating the convective heat transfer coefficient at building surfaces. Build. Environ. 2003, 38, 873–881. [Google Scholar] [CrossRef]
  40. Walton, G.N. Thermal Analysis Research Program Reference Manual; NBSSIR 83-2655; National Bureau of Standards: Washington, DC, USA, 1983; p. 21. [Google Scholar]
  41. Clark, G.; Allen, C. The Estimation of Atmospheric Radiation for Clear and Cloudy Skies. In Proceedings of the 2nd National Passive Solar Conference (AS/ISES), Philadelphia, PA, USA, 16–18 March 1978; pp. 675–678. [Google Scholar]
  42. Kusuda, T. Heat transfer in buildings. Handbook of Heat Transfer Applications, Rohsenow, W.M., Hartnett, J.P., Ganic, E.N., Eds.; McGraw-Hill: New York, NY, USA, 1985. [Google Scholar]
  43. ASHRAE Fundamentals Handbook (S.I.); American Society of Heating, Refrigerating and Air Conditioning Engineers Inc.: Atlanta, GA, USA, 1997.
  44. Sanjuan, C.; Suárez, M.J.; Blanco, E.; Heras, M.D.R. Development and experimental validation of a simulation model for open joint ventilated facades. Energy Build. 2011, 43, 3446–3456. [Google Scholar] [CrossRef]
  45. Daouas, N. Impact of external longwave radiation on optimum insulation thickness in Tunisian building roofs based on a dynamic analytical model. Appl. Energy 2016, 177, 136–148. [Google Scholar] [CrossRef]
  46. Mastrapostoli, E.; Santamouris, M.; Kolokotsa, D.; Vassilis, P.; Venieri, D.; Gompakis, K. On the ageing of cool roofs: Measure of the optical degradation, chemical and biological analysis and assessment of the energy impact. Energy Build. 2016, 114, 191–199. [Google Scholar] [CrossRef]
  47. Eilert, P. High Albedo (Cool) Roofs: Codes and Standards Enhancement (CASE) Study; Pacific Gas and Electric Company: San Francisco, CA, USA, 2000. [Google Scholar]
  48. Xue, X.; Yang, J.; Zhang, W.; Jiang, L.; Qua, J.; Xu, L.; Zhang, H.; Song, J.; Zhang, R.; Li, Y.; et al. The study of an energy efficient cool white roof coating based on styrene acrylate copolymer and cement for waterproofing purpose. Part I: Optical properties, estimated cooling effect and relevant properties after dirt and accelerated exposures. Constr. Build. Mat. 2015, 98, 176–184. [Google Scholar] [CrossRef]
  49. Bretz, S.E.; Akbari, H. Long-term performance of high-albedo roof coatings. Energy Build. 1997, 25, 159–167. [Google Scholar] [CrossRef] [Green Version]
  50. Durability, Pausing for Reflection on Aging of the Cool Roof. J. Archit. Coat. 2008, 63, e5.
  51. Synnefa, A.; Cool Roofs Standars and the ECRC Product Rating Program. European Cool Roofs Council. 2016. Available online: https://coolroofcouncil.eu/wp-content/uploads/2019/05/ECRC-Product-rating-manual_2017.pdf (accessed on 23 March 2021).
  52. Reglamento de Instalaciones Térmicas en los Edificios. Ministerio para la Transición Ecológica y el Reto Demográfico. Available online: https://energia.gob.es/desarrollo/EficienciaEnergetica/RITE/Paginas/InstalacionesTermicas.aspx (accessed on 2 March 2021).
  53. Rossi, M.; Rocco, V.M. External walls design: The role of periodic thermal transmittance and internal areal heat capacity. Energy Build. 2014, 68, 732–740. [Google Scholar] [CrossRef]
  54. Auliac, S.; Le Hyaric, A.; Morice, J.; Hecht, F.; Ohtsuka, K.; Pironneau, O.; FreeFem++. , 3rd ed. Version 3.31-2, 2014. 2017. Available online: http:www.freefem.org (accessed on 5 June 2021).
  55. Agencia Estatal de Meteorología de España. Guía Resumida del Clima en España (1981–2010). 2010. Available online: http://www.aemet.es/es/conocermas/recursos_en_linea/publicaciones_y_estudios/publicaciones/detalles/guia_resumida_2010 (accessed on 2 June 2021).
  56. Domínguez Amarillo, S.; Sendra Salas, J.J.; Oteiza San José, I. La Envolvente Térmica de la Vivienda Social. El Caso de Sevilla, 1939 a 1979; 2016 Csic.; Editorial CSIC: Madrid, Spain, 2016; p. 174. ISBN 978-84-00-10124-4. [Google Scholar]
  57. Cattarin, G.; Causone, F.; Kindinis, A.; Pagliano, L. Empirical and comparative validation of an original model to simulate the thermal behaviour of outdoor test cells. Energy Build. 2018, 158, 1711–1723. [Google Scholar] [CrossRef]
  58. Código Técnico de la Edificación; Ministerio de Fomento del Gobierno de España: Madrid, Spain, 2009.
  59. León, A.L.; Suárez, R.; Bustamante, P.; Campano, M.A.; Moreno, D. Design and performance of test cells as an energy evaluation model of facades in a mediterranean building area. Energies 2017, 10, 1816. [Google Scholar] [CrossRef] [Green Version]
  60. US Department of Energy; Energy Efficiency and Renewable Energy Office; Building Technology Program. EnergyPlus 8.0.0. Available online: http://apps1.eere.energy.gov/buildings/energyplus/ (accessed on 10 June 2021).
  61. Duffie, J.A.; Beckman, W. Solar Engineering of Thermal Processes; John Wiley and Sons Inc.: Hoboken, NJ, USA, 2006. [Google Scholar]
  62. CYPE Engineers. CYPE Ingenieros Construction Price Generator. 2017. Available online: http://www.cype.es/ (accessed on 22 March 2021).
  63. Vincelas, F.F.C.; Ghislain, T.; Robert, T. Influence of the types of fuel and building material on energy savings into building in tropical region of Cameroon. Appl. Therm. Eng. 2017, 122, 806–819. [Google Scholar] [CrossRef]
  64. Available online: https://ec.europa.eu/eurostat/statistics-explained/index.php/Natural-gas-price-statistics (accessed on 2 March 2021).
  65. Nydahl, H.; Andersson, S.; Astrand, A.P.; Olofsson, T. Energy and Buildings Including future climate induced cost when assessing building refurbishment performance. Energy Build. 2019, 203, 109428. [Google Scholar] [CrossRef]
Figure 1. Evolution of solar reflectivity for an initial value equal to 0.9 and a triennial washing maintenance.
Figure 1. Evolution of solar reflectivity for an initial value equal to 0.9 and a triennial washing maintenance.
Sustainability 13 10738 g001
Figure 2. (a) Percent loss of solar reflectivity, monthly and accumulated. (b) Solar reflectivity evolution for the aged case without maintenance. (c) Solar reflectivity evolution for the aged case with decennial washing.
Figure 2. (a) Percent loss of solar reflectivity, monthly and accumulated. (b) Solar reflectivity evolution for the aged case without maintenance. (c) Solar reflectivity evolution for the aged case with decennial washing.
Sustainability 13 10738 g002
Figure 3. Geographical framework and location of the city of Seville (Spain).
Figure 3. Geographical framework and location of the city of Seville (Spain).
Sustainability 13 10738 g003
Figure 4. Climatic chart of Seville.
Figure 4. Climatic chart of Seville.
Sustainability 13 10738 g004
Figure 5. La Juncal (Seville): Aerial views.
Figure 5. La Juncal (Seville): Aerial views.
Sustainability 13 10738 g005
Figure 6. (a) Layout of the reference roof. (b) Layout of the retrofitted roof.
Figure 6. (a) Layout of the reference roof. (b) Layout of the retrofitted roof.
Sustainability 13 10738 g006
Figure 7. Weather conditions for: (a) a winter week (7 to 14 February 2017); (b) a summer week (17 to 24 August 2017).
Figure 7. Weather conditions for: (a) a winter week (7 to 14 February 2017); (b) a summer week (17 to 24 August 2017).
Sustainability 13 10738 g007
Figure 8. Total loads (GJ/m 2 ) for the whole LCA timespan: (a) no washing; (b) decennial washing.
Figure 8. Total loads (GJ/m 2 ) for the whole LCA timespan: (a) no washing; (b) decennial washing.
Sustainability 13 10738 g008
Figure 9. Total loads (GJ/m 2 ) for the whole LCA timespan: (a) quinquennial washing; (b) triennial washing.
Figure 9. Total loads (GJ/m 2 ) for the whole LCA timespan: (a) quinquennial washing; (b) triennial washing.
Sustainability 13 10738 g009
Figure 10. Cooling loads (GJ/m 2 ) for the whole LCA timespan: (a) no washing; (b) decennial washing.
Figure 10. Cooling loads (GJ/m 2 ) for the whole LCA timespan: (a) no washing; (b) decennial washing.
Sustainability 13 10738 g010
Figure 11. Cooling loads (GJ/m 2 ) for the whole LCA timespan: (a) quinquennial washing; (b) triennial washing.
Figure 11. Cooling loads (GJ/m 2 ) for the whole LCA timespan: (a) quinquennial washing; (b) triennial washing.
Sustainability 13 10738 g011
Figure 12. Heating loads (GJ/m 2 ) for the whole LCA timespan: (a) no washing; (b) decennial washing.
Figure 12. Heating loads (GJ/m 2 ) for the whole LCA timespan: (a) no washing; (b) decennial washing.
Sustainability 13 10738 g012
Figure 13. Heating loads (GJ/m 2 ) for the whole LCA timespan: (a) quinquennial washing; (b) triennial washing.
Figure 13. Heating loads (GJ/m 2 ) for the whole LCA timespan: (a) quinquennial washing; (b) triennial washing.
Sustainability 13 10738 g013
Figure 14. No washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Figure 14. No washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Sustainability 13 10738 g014
Figure 15. Decennial washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Figure 15. Decennial washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Sustainability 13 10738 g015
Figure 16. Quinquennial washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Figure 16. Quinquennial washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Sustainability 13 10738 g016
Figure 17. Triennial washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Figure 17. Triennial washing life-cycle costs (GJ/m 2 ): (a) cooling and heating by pump; (b) heating by gas and cooling by pump triennial washing.
Sustainability 13 10738 g017
Table 1. Thermophysical characteristics of the reference roof.
Table 1. Thermophysical characteristics of the reference roof.
LayerDescriptionThickness
(m)
Density
(kg/m3)
Specific Heat
(J/kg K)
Conductivity
(W/m K)
1 (Out.)Ceramic tiles0.00520008001.00
2Mortar0.01200010001.40
3Protective Layer0.015115010000.23
4Mortar0.01200010001.40
5Carbon cinders0.16406571.40
6Concrete vault0.3133010001.32
7 (Ins.)Gypsum plaster0.01100010000.32
Table 2. Thermophysical characteristics of the refurbished roof.
Table 2. Thermophysical characteristics of the refurbished roof.
LayerDescriptionThickness
(m)
Density
(kg/m3)
Specific Heat
(J/kg K)
Conductivity
(W/m K)
1(Out.)White elastomer0.0015115010000.23
2Regularization layer0.01200010001.40
3EPS0 to 0.13012100.04
4–10 (Ins.)Same as layers 1 to 7 of the reference roof (Table 1)
Table 3. Thermophysical characteristics of the roof materials.
Table 3. Thermophysical characteristics of the roof materials.
LayerDescriptionThickness
(m)
Density
(kg/m3)
Specific Heat
(J/kg K)
Conductivity
(W/m K)
1(Out.)Galvanized corrugated
sheet panels0.002780046050.2
2Ventilated air chamber0.2001.18410070.0255
3Sandwich panel0.200401884.150.017
2Wool rock panel0.1601008400.046
3(In.)Sandwich panel0.100401884.150.017
Table 4. Statistical indicators for the experimental validation.
Table 4. Statistical indicators for the experimental validation.
RMSE ( C) ε mean ( C) SD ( ε ) ( C) r 2
Winter0.2930.1790.2130.955
Summer0.3010.2510.3090.961
Table 5. Statistical indicators for the comparison with EnergyPlus.
Table 5. Statistical indicators for the comparison with EnergyPlus.
RMSE ( C) ε mean ( C) SD ( ε ) ( C) r 2
Winter0.2630.0840.1070.985
Summer0.2870.1210.2260.981
Table 6. Total loads (GJ/m 2 ) for the LC timespan.
Table 6. Total loads (GJ/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.05.03505.18895.38825.63355.92164.95445.07985.24915.46325.7190
0.032.27122.31562.38022.46252.56192.34582.38632.44772.53282.6110
0.051.71311.74591.79391.84261.94141.74701.74801.79291.85511.9163
0.071.38351.41211.45271.50121.56711.37491.40661.44591.49841.5588
0.11.08021.13071.13291.17451.21881.07971.10101.12811.17291.2136
Table 7. Total loads (GJ/m 2 ) for the LC timespan.
Table 7. Total loads (GJ/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.04.96785.05925.19935.38845.55774.92064.99815.11785.28135.4856
0.032.38842.42872.48502.56452.64152.34832.39672.45192.52552.6184
0.051.78021.76551.80281.85381.91621.76791.74561.78131.82441.8959
0.071.38351.41351.45071.50161.56181.36621.39431.43191.48051.5467
0.11.07631.11031.13491.17601.21681.06341.09701.12201.16281.2056
Table 8. Total load savings (GJ/m 2 ) for the LC timespan.
Table 8. Total load savings (GJ/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.02.11812.05531.91571.73591.51712.09922.00291.86231.68351.4685
0.034.61684.57244.50784.42554.32614.59604.54744.48064.39074.3049
0.055.17495.14215.09415.04544.94665.18165.14635.09555.02614.9663
0.075.50455.47595.43535.38685.32095.50845.47815.44025.39105.3296
0.15.80785.75735.75515.71355.66925.80835.7875.75995.72295.6753
Table 9. Total load savings (GJ/m 2 ) for the LC timespan.
Table 9. Total load savings (GJ/m 2 ) for the LC timespan.
No WashingDecennial Washing
Insulation
Thickness
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
Abs. =
0.1
Abs. =
0.2
Abs. =
0.3
Abs. =
0.4
Abs. =
0.5
0.02.10202.01381.89001.72771.58652.05611.98631.86421.69701.4892
0.034.58704.53694.47314.38684.30504.58514.53554.47864.40354.3206
0.055.17395.14335.09795.03434.96845.15485.14045.10135.04624.9750
0.075.50645.47905.44445.39635.33285.49955.47895.44825.40485.3446
0.15.81175.77775.75315.72705.68065.82465.79105.76635.73005.6862
Table 10. Input variables used in the LCA.
Table 10. Input variables used in the LCA.
VariableValue
EPS cost116€/m 3
Insulation installation cost13.15€/m 2
Cool paint cost9.45€/m 2
Washing cost1.63€/m 2
Bituminous paint cost8.01€/m 2
Bituminous paint maintenance0.4€/m 2 yearly
Electricity cost0.2403kWh
Gas cost0.0736kWh
Energy inflation rate 3 %
Discount rate 1.5 %
Lifetime20years
Table 11. Conditioning equipment efficiency.
Table 11. Conditioning equipment efficiency.
DeviceEfficiency
Cooling pumpCOP = 3.4
Heating pumpCOP = 3
Gas natural heating η s = 0.8
Electrical Radiators heating η s = 1
Table 12. Optimum results for cooling and heating by pump.
Table 12. Optimum results for cooling and heating by pump.
Cool Roof
Washing
Solar
Absorptivity
Optimum
Thickness (m)
Minimum Total
Cost (€/m 2 )
Maximum total
Savings (€/m 2 )
Payback Period
(Years)
No0.10.0864.27109.944.54
Decennial0.10.0866.51107.704.78
Quinquennial0.10.0868.33105.885.01
Triennial0.10.0872.82101.395.48
Table 13. Optimum results for cooling by pump and heating by gas.
Table 13. Optimum results for cooling by pump and heating by gas.
Cool Roof
Washing
Solar
Absorptivity
Optimum
Thickness (m)
Minimum Total
Cost (€/m 2 )
Maximum total
Savings (€/m 2 )
Payback Period
(Years)
No0.10.0868.27111.804.05
Decennial0.10.0869.81110.264.48
Quinquennial0.10.0971.30108.774.71
Triennial0.20.0874.99105.084.97
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Domínguez-Torres, C.-A.; Domínguez-Torres, H.; Domínguez-Delgado, A. Optimization of a Combination of Thermal Insulation and Cool Roof for the Refurbishment of Social Housing in Southern Spain. Sustainability 2021, 13, 10738. https://doi.org/10.3390/su131910738

AMA Style

Domínguez-Torres C-A, Domínguez-Torres H, Domínguez-Delgado A. Optimization of a Combination of Thermal Insulation and Cool Roof for the Refurbishment of Social Housing in Southern Spain. Sustainability. 2021; 13(19):10738. https://doi.org/10.3390/su131910738

Chicago/Turabian Style

Domínguez-Torres, Carlos-Antonio, Helena Domínguez-Torres, and Antonio Domínguez-Delgado. 2021. "Optimization of a Combination of Thermal Insulation and Cool Roof for the Refurbishment of Social Housing in Southern Spain" Sustainability 13, no. 19: 10738. https://doi.org/10.3390/su131910738

APA Style

Domínguez-Torres, C. -A., Domínguez-Torres, H., & Domínguez-Delgado, A. (2021). Optimization of a Combination of Thermal Insulation and Cool Roof for the Refurbishment of Social Housing in Southern Spain. Sustainability, 13(19), 10738. https://doi.org/10.3390/su131910738

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop