Next Article in Journal
Representing Small Commercial Building Faults in EnergyPlus, Part II: Model Validation
Previous Article in Journal
Probabilistic Seismic Assessment of Existing Masonry Buildings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Algorithm Applied to Multi-Criteria Selection of Thermal Insulation on Industrial Shed Roof

by
Michel Nikolaos Stamoulis
,
Gerson Henrique dos Santos
*,
Wagner Barth Lenz
and
Angelo Marcelo Tusset
Federal University of Technology – Paraná – UTFPR, Mechanical Engineering Graduate Program (PPGEM) Doutor Washington Subtil Chueire, 330 - CEP 84017-220 Ponta Grossa, Brasil
*
Author to whom correspondence should be addressed.
Buildings 2019, 9(12), 238; https://doi.org/10.3390/buildings9120238
Submission received: 25 October 2019 / Revised: 17 November 2019 / Accepted: 19 November 2019 / Published: 21 November 2019

Abstract

:
The rational use of energy has motivated research on improving the energy efficiency of buildings, which are responsible for a large share of world consumption. A strategy to achieve this goal is the application of optimized thermal insulation on a building envelope to avoid thermal exchanges with the external environment, reducing the use of heating, ventilation and air-conditioning (HVAC) systems. In order to contribute to the best choice of insulation applied to an industrial shed roof, this study aims to provide an optimization tool to assist this process. Beyond the thermal comfort and cost of the insulation, some hygrothermic properties also have been analysed to obtain the best insulation option. To implement this optimization technique, several thermo-energetic simulations of an industrial shed were performed using the Domus software, applying 4 types of insulation material (polyurethane, expanded polystyrene, rockwool and glass wool) on the roof. Ten thicknesses ranging from 0.5 cm to 5 cm were considered, with the purpose of obtaining different thermal comfort indexes (PPD, predicted percentage dissatisfied). Posteriorly, the best insulation ranking has been obtained from the weights assigned to the parameters in the objective function, using the technique of the genetic algorithm (GA) applied to multi-criteria selection. The optimization results showed that polyurethane (PU) insulation, applied with a thickness of 1 cm was the best option for the roof, considering the building functional parameters, occupant metabolic activity, clothing insulation and climate conditions. On the other hand, when the Brazilian standard was utilized, rock wool (2 cm) was considered the best choice.

1. Introduction

Heating, ventilation and air-conditioning (HVAC) systems are the main causes of high-electric energy consumption in buildings in the world [1,2]. In order to reduce this energy consumption, the application of thermal insulation on the envelope is an efficient strategy to avoid thermal exchanges with the external environment, thereby increasing the energy efficiency of buildings. The concept of energy efficiency in buildings refers to providing thermal comfort for the occupants with low energy consumption.
Therefore, the design of an energy-efficient building should take into account the external environmental conditions such as geographic location, position, and the climate of the region. In addition, other internal factors are of great importance, such as the purpose of the building, periods of use, and activities of the occupants. Besides these conditions, the materials used on the envelope play a crucial role in the thermal performance of the buildings. In this context, the choice of the best insulation is difficult due to the various parameters that should be analyzed such as performance, cost, uncertainties related to hygrothermal properties, sustainability issues, etc. In this line of research, Adamczyk and Dylewski [3] proposed a methodology for assessing the environmental and economic benefits of thermal insulation. The results indicated that this investment is beneficial for ecological and economic reasons. Soares et al. [4] provided a survey of the main methods to measure the thermal transmittance and thermal behavior of construction elements such as the insulation materials.
In order to promote the most adequate choice, other scientific works found in the literature can also be cited. Firstly, a literature review on determining the optimal thickness of the thermal insulation material in a building envelope and its effect on energy consumption has been carried out by Kaynakli [5]. Kaynakli [6] also investigated some parameters that affect the optimal insulation thickness, such as inflation and discount rates, lifetime, energy costs, the heating/cooling loads, and the properties of the insulation material. The optimal thickness and payback period was determined through an economic model based on life-cycle cost analysis. Akyüz et al. [7] analyzed the application of thermal insulation on the walls and roof of an airport terminal from energy, environment, and cost aspects. The optimal insulation thickness was also calculated and by applying thermal insulation on the walls and roof of the terminal building, a reduction of heat loss of 48% and 56% was observed, respectively.
Lin et al. [8] implemented a multi-linear regression model and an artificial neural network model to predict the building thermal load and the number of discomfort degree hours considering the variable thermal mass, insulation, absorptance of solar radiation, and glazing ratio. The models were coupled with a multi-objective genetic optimization algorithm to minimize the building thermal load and to improve the thermal comfort for a green building in China. An extensive review of the optimization methods and their application in energy-efficient architectural building design can also be found in Kheiri [9]. In addition, decision-making approaches were discussed for multi-objective optimizations.
Dubrow and Krarti [10] presented a simulation–optimization tool that couples a genetic algorithm to a whole-building energy simulation engine to select optimal values of some parameters associated with the envelope to minimize energy use for residential buildings. Different building shapes were investigated and some features were considered in the optimization analysis including wall and roof constructions, foundation types, insulation levels, and window types and areas. Alsayed and Tayeh [11] determined an optimal insulation thickness, which considers climatic conditions, wall structure, type and cost of insulation, cost of energy, and other economic parameters using an economic model, based on life-cycle cost analysis over a building lifetime of 20 years. Energy savings and payback periods analysis were performed for expanded polystyrene (EPS) and polyurethane foam (PU) and the results showed better performance of polyurethane foam for the cases studied.
Grygierek and Grygierek [12] presented a multi-objective optimization of the selected design parameters in a single-family building in temperate climate conditions with natural ventilation. The influence of four types of windows, their size, building orientation, insulation of external wall, roof and ground floor, and infiltration on the life-cycle costs and thermal comfort were analyzed. A controller was proposed, with the task of imitating occupants in controlling the natural ventilation airflow. Baniassadi et al. [13] carried out studies of economic optimization of thermal insulation thickness and phase-change material (PCM) layer using ENERGYPLUS software coupled to a genetic algorithm. Results showed that the optimal PCM thickness was zero for all cases and the optimum insulation thickness varied significantly by climatic regions of the country studied, Iran.
Gelegenis and Axaopoulos [14] developed a complete work in order to indicate the optimal insulation thicknesses of walls’ insulation. A multi-parametric mathematical analytical approach is presented incorporating the existence of glazing at the wall, the thermal absorptance of the wall, the base temperature of the heated space, the thermal characteristics, and the heat and solar gains of the heated space. Other relevant work was undertaken by Rosasco and Perini [15] who investigated which factors are considered by a panel of experts in the fields of design and architecture for the selection of a building roof solution. Results of the study identified their priority and through a sustainability-based multicriteria analysis, the role played by each one in the decision process. Different approaches in the analysis of roof insulation can also be verified in literature. Taylor et al. [16] examined the thermal losses from various roof systems in terms of the economic cost of lost insulation due to thermal bridging. Gürlich et al. [17] analyzed the daylight usage in buildings and quantified the daylight and energy performance of a complex multi-layer textile membrane roof. In this work a genetic algorithm (GA) was been chosen as the optimization tool. The ability of the GA to search simultaneously different regions of a solution space makes it possible to find a diverse set of solutions for difficult problems with non-convex, discontinuous, and multi-modal solution spaces [18]. In building area, Caldas and Norford [19] can be cited as one of the first authors who used the GA applied as a generative and search procedure to look for optimized design solutions in terms of thermal and lighting performance in a building. Although the literature contains works [20] that evaluated the building properties by combining life-cycle assessment methodology and genetic algorithm, this work applied a multi-criteria selection to obtain the best insulation for an industrial shed roof.
In some emerging countries like Brazil, the application of a regulation for energy efficiency in buildings is not yet mandatory and also most industrial buildings do not have air-conditioning systems. Moreover, thermal comfort is not taken into account at the design stage, and the use of a simple insulated roof could improve the well-being and productivity of the occupants. Thus, this work focused on a simple and objective method that can directly contribute to the best choice of thermal insulation in industrial buildings. In this study, four different insulation materials (polyurethane foam—PU, expanded polystyrene—EPS, rockwool and glass wool) with a thickness range from 0.5 cm to 5 cm were considered in order to verify the thermal confort indexes [21]. Beyond the thermal comfort and cost of the insulation, some hygrothermic properties such as water vapor diffusion resistance factor, and thermal diffusivity have also been analysed to obtain the best insulation option. To implement this optimization technique, thermo-energetic simulations of an industrial shed were performed using Domus software [22]. The climate of Curitiba, Brazil (Lat: −25.43 and Long: −49.27) was investigated due to the high daily thermal amplitude presented. Recife, Brazil (Lat: −8.05, Lon: −34.88) was also verified due to its tropical climate. Climatic archive test reference year (TRY) was used in the simulations [23].
This study is in line with the work presented by Civic and Vucijak [24] that used a case study and also proposed the VIKOR method (VlseKriterijumska Optimizacija I Kompromisno Resenje in Serbian) for multi-criteria optimization of insulation. However, in this work the focus is directed to the thermal comfort (higher weight in the objective finction) in which it is one of the most important aspects of the indoor environmental quality due to its effects on well-being of the people and building energy saving, while [24] focused on the emissions aspect of avoiding carbon dioxide emissions during fabrication of insulant materials. The VIKOR method also relies on a fixed dimensionless matrix to evaluates the objective function. On the other hand, GA uses the coefficients and it is more flexible, additionally, GA prevails when non-linearities are involved, because of the genetic operations that can overcome rough transitions.

2. Problem Description

In order to promote the optimal selection of insulation material in the range of the Brazilian market applied to the roof of an industrial shed, the decision criterion has been based in four parameters that should be maximized or minimized:
  • Cost of insulation ($/cm), criteria to be minimized,
  • Predicted percentage dissatisfied (PPD) thermal comfort index (%), criteria to be minimized,
  • Water vapour diffusion resistance factor, criteria to be maximized,
  • Thermal diffusivity (m2/s), criteria to be minimized.
Insulation costs vary according to the density and an average cost of the materials has been adopted to simplify the optimization methodology, as described in Table 1.
The PPD index, created for predicting the general thermal sensation and degree of discomfort, was elaborated by Fanger [21] using heat-balance equations applied to the body human. Environmental parameters as air temperature, mean radiant temperature, air velocity and relative humidity are taken into account for the thermal sensation calculation as well as the clothing insulation and human metabolic rate.
The water vapor diffusion resistance factor is the ratio of the water vapor permeability of the air to the water vapor permeability of the material. The water vapor permeability is the quantity of water vapor transmitted per unit time through a unit area of the material per unit of vapor pressure. Therefore, this parameter is important for the prediction of mold growth and durability of the insulation. On the other hand, the thermal diffusivity measures the ability of a material to conduct thermal energy relative to its ability to store. The internal temperature of buildings with the envelope formed by materials with small thermal diffusivity varies more slowly in response to external climatic conditions.
To simulate the thermal comfort index PPD, an industrial shed with area of 300 m2 (Figure 1) has been modeled in the Domus software [22] developed for the simulation of multizone buildings and analysis of both thermal comfort and energy use. Two cities were considered in the simulations: Curitiba, Brazil (subtropical climate) and Recife, Brazil (tropical climate).
Table 2 presents the data that was considered for the thermal-load simulation of the building model.
Simulations were performed using four types of resistive insulation material. Ten thickness ranging from 0.5 cm to 5 cm were considered, for the purpose of obtaining different thermal comfort indexes (PPD). Table 3 shows the hygrothermal properties of insulations.
Internal thermal gains (equipment and lighting) have not been considered in order to ensure that the thermal load can be attributed only to the external climatic conditions. In the evaluation of thermal comfort indexes, a person performing a moderate activity was considered, with metabolic rate of 1.6 met. A clothing thermal resistance of 0.691 clo (1 clo = 0.155 m2K/W) and air velocity of 0.01 m/s have been adopted.
In the optimization procedure, MATLAB software [26] was chosen to implement the GA. This evolutionary computing technique operates on a population of candidates in parallel. Different areas of the solution space can be searched, allocating an appropriate number of candidate solutions for the search in various regions. Thus, after the objective function is implemented in the GA, the generation, crossing and mutation operations provide the result of the best insulator for the criteria and coefficients applied. For this analysis, the multicriteria decision analysis method (MCDA) was used to perform the classification. MCDA is a common human dilemma that involves decision making strategy different outcomes and different solutions for the same problems. In many cases, the best output it is not explicit and some variables are difficult to quantify. In these complex cases, such as the selection of material for insulation a common option is to select methods of optimization to elucidate the answer. In this case, before implementing the GA, the input data was set, after being normalized, and a weighting factor was arbitrarily attributed to each parameter.
In order to define the objective function, the equation must be associated with the unsatisfied people index (PPD) minimization, insulation cost minimization, vapor diffusion resistance factor maximization and thermal diffusivity minimization, as described by Equation (1):
Fobj = 0.60 PPD + 0.25 Cost + 0.10 μ L + 0.05 α
where µ is the vapor diffusion resistance factor, L, the insulation thickness and α, the thermal diffusivity. A mathematical optimization problem is one in which some function is either maximized or minimized relative to a given set of alternatives. The function to be minimized or maximized is called the objective function and the set of alternatives is called the feasible region or constraint region. In this normalized function (Equation (1)), where the sum of the weighting factor is equal to 1, a factor of 0.6 has been considered for the PPD index, 0.25 for the insulation cost, 0.1 for the vapor-diffusion resistance factor and 0.05, for the thermal diffusivity. These weighting factors can be chosen arbitrarily according to user needs. The cost factor, for example, in underdeveloped countries is more important than in developed countries. On the other hand, the factor associated to vapor diffusion resistance is more important in humid regions than in dry regions.
In Figure 2, the flowchart for obtaining the optimal parameters is presented. Initial conditions are random in the algorithm, generating the initial gens. Consequently, the Domus database returns to the objective function that will check the stop criteria and repeat the operation.
For the variables input in the objective function, normalization has been performed, in which the variables were divided by the maximum value of each criteria. After, the normalizad PPD values (daily average) were grouped in 4 matrices corresponding to the different insulation materials (12 × 10), defined by 12 months and 10 different thicknesses. As the maximum number of population generated is 40 individuals, which refers to the combination of 4 materials and 10 thicknesses, individuals were randomly selected to repeat and compose the initial population. In order to ensure the convergence of the method and adopt parameters suitable for small samples for numerical simulations, the following parameters were defined: a mutation factor of 1% of the individuals; 5% for the survival factor; and 94% for the species’ crossover. According to [27,28,29,30], these parameters are appropriate for a problem similar to that considered in this paper.

3. Results and Discussions

Figure 3 presents the results of the daily mean values of PPD for different thicknesses of the PU insulation, using the Curitiba, Brazil, climatic archive (TRY [23]).
In the months of June to August, the results showed the worst performance, with values between 26% and 39% of PPD. April and November showed the best thermal comfort index, with PPD values varying between 12% and 13%. In relation to the thickness effect, a behavior inversely proportional to the PPD values to the months of winter and summer was observed. In the case of PU insulation, the increase in thickness from 0.5 cm to 5 cm resulted in a maximum reduction of 14% in PPD (from 39% to 26%) in the winter months (South hemisphere) and an increase of 5% in the summer (from 18% to 23%). This effect can be attributed to the mild cold climate of Curitiba, Brazil, and the daily average used, where the nocturnal PPD values have a significant influence on the results. Other insulations had similar behavior to the PU (not shown) due to similar thermal diffusivity.
In order to verify the effect of the insulation for a city with a hot climate (tropical), Figure 4 shows the PPD (%), using the PU insulation (0.5 cm and 5 cm) to the city of Recife, Brazil. Unlike Curitiba, Recife presents less oscillatory results, with a significant improvement of thermal comfort over the year. To verify the clothing effect in the PPD index, PU was used (0.5 cm and 5 cm) on the roof and the Curitiba climatic archive was considered in the simulations. Clothing insulation of 0.69 clo and 1.18 clo also were utilized in the analysis.
Figure 5 shows that when a high clothing insulation (1.18 clo) is adopted, the effect of the insulation on the roof is not significant in the winter months for the PPD values. However, for a low clothing insulation (0.69 clo) in this period, PPD showed a reduction of 10% with increasing insulation thickness on the roof. For the summer, the increase in thickness of the insulation decreased the PPD by 8% (1.18 clo). In the case of low clothing insulation, the opposite behavior was observed, suggesting the need for an increase in the thermal load for a PPD reduction.
In the optimization process, a minimal population was analyzed, shown in Figure 6.
As shown in Figure 6, after 4 generations there was not much improvement from the initial conditions and the local best prevail. The minimum requirement for convergence was 5 population that allowed the population to have genetic operations along the generations to improve the results.
The objective function had the best result when its value tended to a minimum. In the PU case (Figure 7), the best values were found between March and April, for thicknesses of 0.5 cm to 1.0 cm, and in December, for the same thicknesses. From June to August, the objective function values presented the worst results, with behavior similar to the PPD values, shown in Figure 3. This characteristic was expected due to the PPD weight of 60% attributed in the objective function.
Table 4 presents the maximum and minimum values obtained in the objective function for all insulations. According to the results, PU and the rock wool had the lowest values for the objective function. However, taking into account the other proposed criteria and following the convergence criteria, it is possible to establish a ranking of the insulators with the respective thickness.
Table 5 shows that the PU insulation with a thickness of 1 cm presented the lowest value of the objective function. In sequence, the rock wool presented the best results (some thicknesses), followed by fiber glass and EPS, respectively.
NBR 15220 (Brazilian Standard [25]) establishes an overall heat transfer coefficient lower than 2 W/m2K for the roof. Table 6 presents the approximate thickness of each insulation obtained through the equivalent thermal resistance.
Through the ranking presented in Table 5, the position and the classification for each insulation is shown in Table 7. In the choice, the nearest and lowest thickness has been considered. As observed in Table 7, the rock wool (2 cm) is the best choice for the application in an industrial shed roof, considering the criteria of selection, building configuration and climate conditions.

4. Summary and Conclusions

This work aimed to provide an optimization tool able to assist in the choice of the best thermal insulation to be applied on the roof of an industrial shed. Beyond thermal comfort and costs, the vapor diffusion resistance factor and thermal diffusivity were also considered in the insulation selection.
Thermal comfort indexes (PPD) have been obtained through thermo-energetic simulations using the Domus software. In the analysis, 4 thermal insulations (polyurethane (PU), expanded polystyrene (EPS), rock wool and glass wool) were considered and 10 different thicknesses (0.5 to 5 cm) were verified. The best choice was obtained using the genetic algorithm (GA) technique applied to the multicriteria selection. The optimization results showed that the polyurethane (PU) insulation applied at a thickness of 1 cm was the best option for installation on the roof, considering the functional parameters of the building, occupant metabolic activity, clothing insulation and climate conditions.
Taking into account the Brazilian standard (NBR 15220) which establishes an overall heat-transfer coefficient lower than 2 W/m2K for insulated roofs, a rock wool insulation with 2 cm of thickness presented the best classification.
The ranking of the insulations was obtained through the weighting factors chosen arbitrarily for the objective function. Other weighting factors could lead to a different ranking. For example, in underdeveloped countries, the term cost would have a greater weight in objective function. On the other hand, in humid regions, the vapor diffusion resistance factor could be more significant. In addition, other clothing insulation or human metabolic rates used to calculate the PPD could modify the results presented. However, the use of the genetic algorithm proved to be capable of generating a diverse set of solutions and is able to capture extreme solutions to any objective function. For future work, analysis of an adaptive thermal comfort model in the choice of the best insulation is suggested. Other parameters of a building also can be analyzed in order to increase the energy efficiency of buildings, such as: building orientation, window area and localization, paint color and finishes, among others.

Author Contributions

This research has been developed within the framework of M.N.S. dissertation under the guidance of G.H.d.S. (Supervisor) and jointly supervised by A.M.T. In this work, the contributions of the authors were: writing and simulation, M.N.S.; supervision and review, G.H.d.S.; software and investigation, W.B.L. and A.M.T.

Funding

The authors thank the Brazilian funding agencies CAPES and CNPq for research support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berardi, U. Building Energy Consumption in US, EU, and BRIC Countries. Procedia Eng. 2015, 118, 128–136. [Google Scholar] [CrossRef]
  2. International Energy Agency (IEA). Transition to Sustainable Buildings: Strategies and Opportunities to 2050; International Energy Agency: Paris, France, 2013. [Google Scholar]
  3. Adamczyk, J.; Dylewski, R. The impact of thermal insulation investments on sustainability in the construction sector. Renew. Sustain. Energy Rev. 2017, 80, 421–429. [Google Scholar] [CrossRef]
  4. Soares, N.; Martins, C.; Gonçalves, M.; Santos, P.; Da Silva, L.S.; Costa, J.J. Laboratory and in-situ non-destructive methods to evaluate the thermal transmittance and behavior of walls, windows, and construction elements with innovative materials: A review. Energy Build. 2019, 182, 88–110. [Google Scholar] [CrossRef]
  5. Kaynakli, O. A review of the economical and optimum thermal insulation thickness for building applications. Renew. Sustain. Energy Rev. 2012, 16, 415–425. [Google Scholar] [CrossRef]
  6. Kaynakli, O. Parametric Investigation of Optimum Thermal Insulation Thickness for External Walls. Energies 2011, 4, 913–927. [Google Scholar] [CrossRef]
  7. Akyüz, M.K.; Altunta, O.; Sögüt, M.Z. Economic and Environmental Optimization of an Airport Terminal Building’s Wall and Roof Insulation. Sustainability 2017, 9, 1849. [Google Scholar] [CrossRef]
  8. Lin, Y.; Zhou, S.; Yang, W.; Li, C.-Q. Design Optimization Considering Variable Thermal Mass, Insulation, Absorptance of Solar Radiation, and Glazing Ratio Using a Prediction Model and Genetic Algorithm. Sustainability 2018, 10, 336. [Google Scholar] [CrossRef]
  9. Kheiri, F. A review on optimization methods applied in energy-efficient building geometry and envelope design. Renew. Sustain. Energy Rev. 2018, 92, 897–920. [Google Scholar] [CrossRef]
  10. Dubrow, D.T.; Krarti, M. Genetic-algorithm based approach to optimize building envelope design for residential buildings. Build. Environ. 2010, 45, 1574–1581. [Google Scholar] [CrossRef]
  11. Alsayed, M.F.; Tayeh, R.A. Life cycle cost analysis for determining optimal insulation thickness in Palestinian buildings. J. Build. Eng. 2019, 22, 101–112. [Google Scholar] [CrossRef]
  12. Grygierek, K.; Ferdyn-Grygierek, J. Multi-Objective Optimization of the Envelope of Building with Natural Ventilation. Energies 2018, 11, 1383. [Google Scholar] [CrossRef]
  13. Baniassadi, A.; Sajadi, B.; Amidpour, M.; Noori, N. Economic Optimization of PCM and insulation layer thickneess in residential buildings. Sustain. Energy Technol. Assess. 2016, 14, 92–99. [Google Scholar]
  14. Gelegenis, J.; Axaopoulos, P. A Multi-Parametric Mathematical Approach on the Selection of Optimum Insulation Thicknesses in Buildings. Buildings 2017, 7, 15. [Google Scholar] [CrossRef]
  15. Rosasco, P.; Perini, K. Selection of (Green) Roof Systems: A Sustainability-Based Multi-Criteria Analysis. Buildings 2019, 9, 134. [Google Scholar] [CrossRef]
  16. Taylor, T.J.; Willits, J.; Hartwig, C.A.; Kirby, J.R. Optimizing Single-Ply Low-Slope Roofing Assemblies for Insulation Value. Buildings 2018, 8, 64. [Google Scholar] [CrossRef]
  17. Gürlich, D.; Reber, A.; Biesinger, A.; Eicker, U. Daylight Performance of a Translucent Textile Membrane Roof with Thermal Insulation. Buildings 2018, 8, 118. [Google Scholar] [CrossRef]
  18. Konak, A.; Coit, D.W.; Smith, A.E. Multi-objective optimization using genetic algorithms: A tutorial. Reliab. Eng. Syst. Saf. 2006, 91, 992–1007. [Google Scholar] [CrossRef]
  19. Caldas, L.G.; Norford, L.K. A design optimization tool based on a genetic algorithm. Autom. Constr. 2002, 11, 173–184. [Google Scholar] [CrossRef]
  20. Pushkar, S.; Becker, R.; Katz, A. A methodology for design of environmentally optimal buildings by variable grouping. Build. Environ. 2005, 40, 1126–1139. [Google Scholar] [CrossRef]
  21. Fanger, P.O. Thermal Comfort: Analysis and Applications in Environment Engineering; McGraw Hill Inc.: New York, NY, USA, 1972. [Google Scholar]
  22. Mendes, N.; Oliveira, R.C.L.F.; Santos, G.H. Domus 2.0: A Whole-Building Hygrothermal Simulation Program. In Proceedings of the Eighth International Conference on Building Performance Simulation (IBPSA 03), Eindhoven, The Netherlands, 11–14 August 2003. [Google Scholar]
  23. CIBSE. Chartered Institution of Building Services Engineers. UK. Available online: https://www.cibse.org/weatherdata (accessed on 15 February 2019).
  24. Civic, A.; Vucijak, B. Multi-criteria Optimization of Insulation Options for Warmth of Buildings to Increase Energy Efficiency. Procedia Eng. 2014, 69, 911–920. [Google Scholar] [CrossRef]
  25. Associação Brasileira de Normas Técnicas (ABNT). NBR 15220 (Desempenho Térmico de Edificações); Associação Brasileira de Normas Técnicas: Rio de Janeiro, Brazil, 2003. (In Portuguese) [Google Scholar]
  26. MATLAB. The Language of Technical Computing; The MathWorks, Inc.: Natick, MA, USA, 2018. [Google Scholar]
  27. Lynch, M.; Conery, J.; Bürger, R. Mutation Accumulation and the Extinction of Small Populations. Am. Nat. 1995, 146, 489–518. [Google Scholar] [CrossRef]
  28. Goldberg, D.E. Genetic Algorithms in Search, Optimization & Machine Learning; Addison-Wesley: Boston, MA, USA, 1989. [Google Scholar]
  29. Conn, A.R.; Gould, N.I.M.; Toint, P. A Globally Convergent Augmented Lagrangian Algorithm for Optimization with General Constraints and Simple Bounds. SIAM J. Numer. Anal. 1991, 28, 545–572. [Google Scholar] [CrossRef]
  30. Conn, A.R.; Gould, N.I.M.; Toint, P.L. A Globally Convergent Augmented Lagrangian Barrier Algorithm for Optimization with General Inequality Constraints and Simple Bounds. Math. Comput. 1997, 66, 261–288. [Google Scholar] [CrossRef]
Figure 1. Industrial shed model.
Figure 1. Industrial shed model.
Buildings 09 00238 g001
Figure 2. Flowchart applied in optimization by genetic algorithm (GA).
Figure 2. Flowchart applied in optimization by genetic algorithm (GA).
Buildings 09 00238 g002
Figure 3. PPD (%) for 10 thicknesses of the PU insulation used on the roof.
Figure 3. PPD (%) for 10 thicknesses of the PU insulation used on the roof.
Buildings 09 00238 g003
Figure 4. PPD (%) for 2 thicknesses of the PU insulation used on the roof in Recife, Brazil.
Figure 4. PPD (%) for 2 thicknesses of the PU insulation used on the roof in Recife, Brazil.
Buildings 09 00238 g004
Figure 5. Clothing insulation effects on PPD.
Figure 5. Clothing insulation effects on PPD.
Buildings 09 00238 g005
Figure 6. Population sensibility.
Figure 6. Population sensibility.
Buildings 09 00238 g006
Figure 7. Objective function for the PU insulation.
Figure 7. Objective function for the PU insulation.
Buildings 09 00238 g007
Table 1. Insulation costs to cover an area of 1 m2 of roof.
Table 1. Insulation costs to cover an area of 1 m2 of roof.
InsulationTypeDensity (kg/m3)Insulation Cost ($/cm)
EPSSlab401.25
PUSlab324.25
Fiber GlassRoll500.58
Rock WoolRoll500.19
Table 2. Summary of building input data.
Table 2. Summary of building input data.
Geographic LocationCuritiba (Lat: −25.43, Lon: −49.27) and Recife (Lat: −8.05, Lon: −34.88)
Building OrientationNorth
Area300 m2
Geometric Dimensions12 m × 25 m × 4 m
Roof Area300 m2
Configuration WallsMortar (20 mm) + Brick (240 mm) + Mortar (20 mm)
Dimension Doors (2)18 m2 (6 m × 3 m)
Window to Wall Ratio20%
Windows Area (8)5 m2 (5 m × 1 m)
Table 3. Hygrothermal property of roof insulation layers [25].
Table 3. Hygrothermal property of roof insulation layers [25].
MaterialsDensity (kg/m3)Specific Heat (J/kgK)Water Vapor Diffusion Resistance Factor (-)Thermal Conductivity (W/mK)Thermal Diffusivity (m2/s)
EPS401420300.0407.04 × 10−7
PU321670440.0356.55 × 10−7
Fiber Glass507001.50.0451.29 × 10−6
Rock Wool507501.30.0451.20 × 10−6
Table 4. Maximum and minimum values in the objective function for the 4 isolations.
Table 4. Maximum and minimum values in the objective function for the 4 isolations.
Objective Function ValuePUEPSRock WoolFiber Glass
Maximum0.4450.6390.4350.586
Minimum0.3750.5910.3850.532
Table 5. Ranking of the optimal configuration of the insulations.
Table 5. Ranking of the optimal configuration of the insulations.
RankInsulationThickness (cm)GA Value
1PU1.00.375
2PU0.50.379
3Rock Wool0.50.385
4Rock Wool1.00.389
5Rock Wool1.50.394
6Rock Wool2.00.403
7Rock Wool3.00.406
8PU2.00.408
9Rock Wool3.50.409
10PU1.50.413
11PU3.00.415
12Rock Wool2.50.42
13Rock Wool4.00.421
14PU3.50.424
15PU2.50.425
16Rock Wool4.50.429
17PU4.00.431
18Rock Wool50.435
19PU4.50.438
20PU50.445
21Fiber Glass3.00.532
22Fiber Glass3.50.544
23Fiber Glass0.50.546
24Fiber Glass1.00.550
25Fiber Glass4.00.554
26Fiber Glass1.50.555
27Fiber Glass4.50.559
28Fiber Glass2.50.584
29Fiber Glass2.00.584
30Fiber Glass50.586
31EPS0.50.591
32EPS1.00.598
33EPS1.50.603
34EPS2.00.608
35EPS2.50.614
36EPS3.00.619
37EPS3.50.624
38EPS4.50.629
39EPS4.00.629
40EPS50.639
Table 6. Insulation thickness calculated according to NBR 15220 [25].
Table 6. Insulation thickness calculated according to NBR 15220 [25].
InsulationThickness (cm)Thermal Conductivity (W/mK)
EPS2.000.040
PU1.750.035
Fiber Glass2.250.045
Rock Wool2.250.045
Table 7. Ranking and classification of the insulations following the NBR 15220 standard.
Table 7. Ranking and classification of the insulations following the NBR 15220 standard.
Rock WoolPUFiber GlassEPS
Thickness (cm)2.01.52.02.0
GA value 0.40300.41300.58400.6080
Ranking Position6102934
Classification1234

Share and Cite

MDPI and ACS Style

Stamoulis, M.N.; Santos, G.H.d.; Lenz, W.B.; Tusset, A.M. Genetic Algorithm Applied to Multi-Criteria Selection of Thermal Insulation on Industrial Shed Roof. Buildings 2019, 9, 238. https://doi.org/10.3390/buildings9120238

AMA Style

Stamoulis MN, Santos GHd, Lenz WB, Tusset AM. Genetic Algorithm Applied to Multi-Criteria Selection of Thermal Insulation on Industrial Shed Roof. Buildings. 2019; 9(12):238. https://doi.org/10.3390/buildings9120238

Chicago/Turabian Style

Stamoulis, Michel Nikolaos, Gerson Henrique dos Santos, Wagner Barth Lenz, and Angelo Marcelo Tusset. 2019. "Genetic Algorithm Applied to Multi-Criteria Selection of Thermal Insulation on Industrial Shed Roof" Buildings 9, no. 12: 238. https://doi.org/10.3390/buildings9120238

APA Style

Stamoulis, M. N., Santos, G. H. d., Lenz, W. B., & Tusset, A. M. (2019). Genetic Algorithm Applied to Multi-Criteria Selection of Thermal Insulation on Industrial Shed Roof. Buildings, 9(12), 238. https://doi.org/10.3390/buildings9120238

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