Optimization of Multilayered Walls for Building Envelopes Including PCM-Based Composites

This work proposes a numerical procedure to simulate and optimize the thermal response of a multilayered wallboard system for building envelopes, where each layer can be possibly made of Phase Change Materials (PCM)-based composites to take advantage of their Thermal-Energy Storage (TES) capacity. The simulation step consists in solving the transient heat conduction equation across the whole wallboard using the enthalpy-based finite element method. The weather is described in detail by the Typical Meteorological Year (TMY) of the building location. Taking the TMY as well as the wall azimuth as inputs, EnergyPlusTM is used to define the convective boundary conditions at the external surface of the wall. For each layer, the material is chosen from a predefined vade mecum, including several PCM-based composites developed at the Institut für Werkstoffe im Bauwesen of TU Darmstadt together with standard insulating materials (i.e., EPS or Rockwool). Finally, the optimization step consists in using genetic algorithms to determine the stacking sequence of materials across the wallboard to minimize the undesired heat loads. The current simulation-based optimization procedure is applied to the design of envelopes for minimal undesired heat losses and gains in two locations with considerably different weather conditions, viz. Sauce Viejo in Argentina and Frankfurt in Germany. In general, for each location and all the considered orientations (north, east, south and west), optimal results consist of EPS walls containing a thin layer made of the PCM-based composite with highest TES capacity, placed near the middle of the wall and closer to the internal surface.


Introduction
According to the International Energy Agency (IEA), the buildings and construction sector across the world was responsible for 36% of final energy use and 39% of carbon dioxide (CO 2 ) emissions in 2018 [1]. In Germany-and in Europe as well-buildings consume almost 40% of the total energy [2]. Far from Europe, in Argentina, a similar situation is faced: buildings are the largest energy consumers (33% of the total energy) [3].
Engineers, scientists, and decision makers all around the world are being appealed to promote best construction practices for using the least, cleanest, and/or most economical resources. The optimization The paper is structured as follows. Section 2 introduces the optimization-based methodology for the design of multilayered walls for building envelopes to be optimal in the light of a given performance index (as example, a purely thermal criterion is defined here). Section 3 describes the thermal modeling of a multilayered wall in a building envelope, pointing out the influence of local weather conditions on boundary conditions and determining the temperature as a function of the design. Section 4 is devoted to describing the vade mecum of materials for building envelopes, with emphasis on the PCM-based composites it contains. Section 5 shows and discuss the results of applying this methodology to walls with different orientations located in two different climates. Finally, Section 6 addresses the concluding remarks and future work.

Optimization-Based Design of Wallboards for Building Envelopes
Let us consider the N-layer wallboard with total thickness L. As shown in Figure 1, ∈ [0, ] denotes the distance to the external surface and the layer = 1,2, … , has thickness Δ and lies between and = + Δ . To optimize the performance of the wall, each layer is allowed to be made of a different material. Let the integer identify the material at layer = 1,2, … , . In the context of optimization, plays the role of a categorical design variable that can take a finite number of levels. Each level is an integer that identifies a candidate material. In general, each can take its own levels; in other words, the layer can have its particular set of material choices; for instance, layer 1 (the external one) may only allow hydrophobic and non-flammable materials.
Each possible design of the -layered wall is identified by the set = [ , , … , ]. Then, we propose to determine the optimal design, say opt , by solving the following optimization problem (Equation (1)): where is the cost function representing the performance of the wall. Several choices for can be found in the literature: the energy consumption for comfort in airconditioned rooms and the degree-hours of discomfort in naturally-ventilated rooms, which can be considered either as multiple objectives [28,29] or combined in a weighted sum as a unique objective [27], the life cycle cost [30], the environmental impact [30,31], and a sum of the initial investment and the energy bill minus the resale value [32]. In these works, the cost function involves the whole building, but they can perfectly serve to characterize the performance of the wallboards enveloping these buildings.
In general, the problem in Equation (1) is subject to inequality constraints that serve to prescribe specific thresholds for the thermal transmittance (U-value) [32], the structural compliance, the weight of the whole wall, etc. The objective and these may be interchangeable. For instance, one could either minimize the energy demand without exceeding a given budget or minimize the budget without exceeding a given energy demand. To optimize the performance of the wall, each layer is allowed to be made of a different material. Let the integer d i identify the material at layer i = 1, 2, . . . , N. In the context of optimization, d i plays the role of a categorical design variable that can take a finite number of levels. Each level is an integer that identifies a candidate material. In general, each d i can take its own levels; in other words, the layer i can have its particular set of material choices; for instance, layer 1 (the external one) may only allow hydrophobic and non-flammable materials.
Each possible design of the N-layered wall is identified by the set d = [d 1 , d 2 , . . . , d N ]. Then, we propose to determine the optimal design, say d opt , by solving the following optimization problem (Equation (1)): where C is the cost function representing the performance of the wall. Several choices for C can be found in the literature: the energy consumption for comfort in air-conditioned rooms and the degree-hours of discomfort in naturally-ventilated rooms, which can be considered either as multiple objectives [28,29] or combined in a weighted sum as a unique objective [27], the life cycle cost [30], the environmental impact [30,31], and a sum of the initial investment and the energy bill minus the resale value [32]. In these works, the cost function involves the whole building, but they can perfectly serve to characterize the performance of the wallboards enveloping these buildings.
In general, the problem in Equation (1) is subject to inequality constraints that serve to prescribe specific thresholds for the thermal transmittance (U-value) [32], the structural compliance, the weight of the whole wall, etc. The objective C and these may be interchangeable. For instance, one could Materials 2020, 13, 2787 4 of 13 either minimize the energy demand without exceeding a given budget or minimize the budget without exceeding a given energy demand.
For example, let us consider a purely thermal performance criterion adopting as cost function the total undesired heat loads all along the 8760 h of a typical meteorological year (TMY) (Equation (2)): where h int is the heat convection coefficient between the wall and the indoor environment, T (h) surfint is the temperature at the internal wall surface at hour h, T max tgt and T max tgt are the maximum and minimum target (ideal) indoor temperatures (maybe time-dependent), and u = (u + |u|)/2 is the ramp function such that there is no contribution to the undesired heat gains C gain if T In the above equation, h int , T max tgt , and T min tgt are assumed to be given, while T (h) surfint is determined by the design d in a way to be defined in the next section.

Temperature Evolution Across the Wall
The temperature T at a distance x ∈ [0, L] from the external surface of the wall at the time instant t > 0 is governed by the heat conduction Equation (3) subject to the initial conditions (Equation (4)) T(x, 0) = T 0 (4) and the boundary conditions (Equations (5) and (6)) In the above equations . T = ∂T/∂t is the temperature rate, ρc eff is the effective heat capacity, k is the thermal conductivity, q ext is the heat flux from the outdoor environment through the external wall surface (which evolves in time following the local weather conditions as described in Section 3.1), h int is the heat convection coefficient at the internal wall surface (the same as in Equation (2)) and T room is the indoor room temperature.
Given the multilayered nature of the wall, the physical properties k and ρc eff are layer-wise defined: at a distance x ∈ (x i−1 , x i ) from the external surface, they are those of the material in layer i, that is the material d i in the vade mecum (see Figure 1).
In general, material properties depend on temperature (Equations (7)). Further, the effective heat capacity in PCMs is rate-dependent during phase changes [33]. Thus, The temperature dependence of material properties makes the heat Equation (3) nonlinear. This nonlinearity becomes severe in presence of PCMs due to the peaks in their heat capacity during Materials 2020, 13, 2787 5 of 13 phase changes. Then, it is crucial to develop a robust solver of the heat conduction Equation (3). Here, recourse is made to the enthalpy-based finite element formulation proposed by Morgan et al. [24]. For the sake of conciseness, the reader interested in the finite element implementation is referred to the just cited work.
Finally, since the thermal properties k and ρc eff at each layer i are dependent on d i , it becomes apparent that the temperature T at any point x across the wall at any time instant t depends on the , making explicit the influence of the design d on the performance of the wall as represented by the cost function C(d) given by Equation (2).

External Boundary Conditions
The heat flux from the outdoor environment is defined as follows (Equations (8)): where q solar is the absorbed short-wave (direct and diffuse) solar radiation, α is the solar absorbance (assumed equal to 0.6 h ext is the external heat transfer coefficient T out is the outdoor temperature, and T surfext (t) = T(0, t, d) is the temperature at the external surface. The heat flux q ext depends on t not only via T surfext (t) but also through T out , h ext and q solar , which change following the instantaneous local weather conditions as defined by the Typical Meteorological Year (TMY).
Regarding T out , this coincides with either the dry bulb temperature T db for h ext < 1000 W/(m 2 K) or the wet-bulb temperature T wb for h ext = 1000 W/(m 2 K), with T db and T wb directly taken from the TMY, where they vary on an hourly basis. On the other hand, q solar and the convection coefficient h ext are computed using EnergyPlus™ [34]. To take into account not only the local weather conditions but also the influence of the surface facing angle on these variables, EnergyPlus is applied to a square 4-m-wide, 4-m-deep, and 3-m-high thermal zone with walls facing north, east, south, and west. Regarding h ext , it is computed using the so-called AdaptiveConvectionAlgorithm in EnergyPlus, which takes into account the wind direction and magnitude and sets h ext to an arbitrarily high value (1000 W/(m 2 K)) at the wall exposed to the wind when it is raining, forcing to assume T out = T db at that instant.

Vade Mecum of Materials for Building Envelopes
To improve the performance of a multilayer wall, let each layer be built of a material thoughtfully chosen from the henceforth called vade mecum of materials for building envelopes. This is a database that should contain a wide choice of materials in terms of thermal properties (conductivity, specific heat, and thermal energy storage capacity, etc.) as well as non-flammability, water and air tightness, weight, cost, and embedded energy, among others, such that it offers a large enough design space for optimization purposes. Mathematically speaking, each material in the vade mecum is a level of the categorical design variable d i . Note that a vade mecum containing M materials gives rise to N M different designs for an N-layered wall.
As a first step, the current vade mecum is built on the basis of purely thermal criteria. Particularly, we are interested in using two cement-based PCM-composites, henceforth referred to as MPCM-p and RBA-p, with p related to the PCM content.
The MPCM-p is a concrete with w/c = 0.45, 70 vol% of normal aggregates (granitic crushed stones) and 30 vol% of a PCM paste containing p =0, 10, or 20 vol% of microencapsulated PCM. The PCM is Micronal ® DS 5038 X type, which is a powder of microencapsulated paraffin wax developed by BASF, with a melting point of around 26 • C and a heat storage capacity of 145 kJ/kg [35]. The effective heat capacity of MPCM-p mixtures as distinct temperature-dependent functions for either heating or cooling are shown on the left of Figure 2 (see [25] for more details on these mixtures).
The RBA-p is a concrete with w/c = 0.50, 30 vol% of surrounding cement paste, and 70 vol% of recycled brick aggregates (RBA), which are filled with p = 0, 65, and 80 vol% of PCM. In this case, the PCM is the non-encapsulated paraffin wax RT 25 HC ® developed by RUBITHERM ® , with a melting point of around 25 • C and a heat storage capacity of 210 kJ/kg [36]. The effective heat capacity of MPCM-p mixtures are also temperature-dependent functions that differ for heating and cooling, and are depicted on the right of Figure 2. More details on RBA-p composites are given in [17].
To enlarge the choice of thermal properties, the vade mecum also contains two widely used insulating materials: rockwool and expanded polystyrene (EPS). The current vade mecum is summarized in Table 1. Let us note that this is a first version, to be continuously enriched with additional properties (cost, embedded energy, etc.) and enlarged all along the recently launched project NRG-STORAGE, mainly to account for other performance indexes (cost, environmental aspects, etc.). Further, the wall must exhibit fire safety and water tightness, among other essential requirements for real life applications. Considering fire safety for instance, PCM-based concretes offer a non-flammable choice; on the contrary, EPS and wools exhibit a high flammability, which is a well-known critical issue yet to be solved. Thus, although indispensable, these additional requirements may be detrimental to the performance of the wall; therefore, in the optimization problem in Equation (1), they are not represented by the cost function C but by inequality constraints. requirements for real life applications. Considering fire safety for instance, PCM-based concretes offer a non-flammable choice; on the contrary, EPS and wools exhibit a high flammability, which is a well-known critical issue yet to be solved. Thus, although indispensable, these additional requirements may be detrimental to the performance of the wall; therefore, in the optimization problem in Equation (1), they are not represented by the cost function but by inequality constraints.

Numerical Results
This section reports the results of the optimization of multilayered systems for external walls, considering four different orientations (N, S, W, and E) and two different locations (Sauce Viejo in

Numerical Results
This section reports the results of the optimization of multilayered systems for external walls, considering four different orientations (N, S, W, and E) and two different locations (Sauce Viejo in Argentina and Frankfurt in Germany). Sauce Viejo has a humid subtropical climate, which is considered as Cfa in the Köppen-Geiger classification; the average temperatures along typical summer and winter weeks are 24.4 and 13.8 • C, respectively. Frankfurt has a warm temperate climate, Cfb in the Köppen-Geiger classification, with average temperatures of 18.0 and 1.7 • C, respectively, for typical summer and winter weeks. For later discussions concerning orientation, let us keep in mind that Sauce Viejo and Frankfurt are in the Southern and Northern hemispheres, respectively.
For an accurate evaluation of the performance of a building envelope all along a year, the weather at each location is described by its typical meteorological year (TMY). This TMY is a database with relevant weather variables (including dry-bulb and dew-point temperatures, wind speed and direction, relative humidity, total sky cover, ceiling height, atmospheric pressure, global horizontal solar radiation, diffuse and direct normal solar radiation, precipitation, etc.) given on an hourly basis along twelve typical meteorological months, which were chosen from different years following the Sandia method [37]. The TMYs for Sauce Viejo and Frankfurt can be downloaded for free from Climate.OneBuilding [38]; that at Sauce Viejo was recently generated by Bre and Fachinotti [39].
As pointed out in Section 2, the cost function C(d) (to be minimized) represents the undesired heat loads along a TMY defined by Equation (2). We further assume that: (i) the room temperature T room appearing in the boundary condition at the internal surface (Equation (6)) is ideally maintained at 24 • C; (ii) the maximal and minimal temperatures involved in the definition of C(d) (Equation (2)) are set to the same value, i.e., T min tgt ≡ T min tgt ≡ T room = 24 • C; and (iii) the convection coefficient of the internal surface in Equations (2) and (6) is set to the typical value h int = 8.24 W/ m 2 K , as adopted by Biswas et al. [21,22] for non-reflective vertical interior walls.

Reference Solutions on Homogenous Walls
For comparison purposes, let us start by considering a homogeneous 20-cm-thick external wall made of one of the candidate materials found in the vade mecum (Table 2). The transient heat conduction Equation (3) was solved using the enthalpy-based finite element method [24] with Euler-backward (implicit) time stepping. A previous analysis was carried out to determine the best deal between and computational cost, from which we decided to use four linear finite elements per layer and constant time steps of half-an-hour (i.e., 17,520 time steps along the TMY). To take due account of the steep variation of the effective heat capacity in PCM-based composites, six Gauss points were used in the corresponding finite elements, while two Gauss points (as usual) were used in the remaining finite elements.
Then, 64 problems were solved, one per each location, material, and orientation. From now on, for the sake of simplicity, let us refer simply as N to the case of a wall facing N, and so on for E-, S, and W-facing cases.
At Sauce Viejo, as shown in Figure 3, the undesired heat gains C gain and losses C loss are balanced. Using the PCM-concrete composites, C gain is maximal for the N and minimal for S; for W and E, it is almost as prejudicial as for N. Regarding C loss , its maximum and minimum occur at S and N, respectively. Once again, E and W are closer to the worst case (S). In general, considering the total undesired loads C = C gain + C loss , the worst case is W, followed by E, N and S, in that order. Using insulating materials (either rockwool or EPS), the total C considerably reduced. Best performances are achieved using EPS (the least conductive material of the vade mecum). For EPS, C attains its minimum value for S. At Frankfurt, as shown in Figure 4, loss is considerably higher than gain . This made the conclusions for loss also valid for the total heat loads . By using the PCM-concrete composites, the performance of the wall is not as sensitive to the orientation as it is for Sauce Viejo. In general, N is the worst case, but E, W and S are not much better. Once again, the performance is greatly improved by using insulating materials.

Optimal Multilayered Walls
Given the 20-cm-thick external wall of the preceding section, let it be made of 20 equally thick layers, and each layer is allowed to be made of one the eight materials in the vade mecum (Table 1). At Frankfurt, as shown in Figure 4, C loss is considerably higher than C gain . This made the conclusions for C loss also valid for the total heat loads C. By using the PCM-concrete composites, the performance of the wall is not as sensitive to the orientation as it is for Sauce Viejo. In general, N is the worst case, but E, W and S are not much better. Once again, the performance is greatly improved by using insulating materials. At Frankfurt, as shown in Figure 4, loss is considerably higher than gain . This made the conclusions for loss also valid for the total heat loads . By using the PCM-concrete composites, the performance of the wall is not as sensitive to the orientation as it is for Sauce Viejo. In general, N is the worst case, but E, W and S are not much better. Once again, the performance is greatly improved by using insulating materials.

Optimal Multilayered Walls
Given the 20-cm-thick external wall of the preceding section, let it be made of 20 equally thick layers, and each layer is allowed to be made of one the eight materials in the vade mecum (Table 1).

Optimal Multilayered Walls
Given the 20-cm-thick external wall of the preceding section, let it be made of 20 equally thick layers, and each layer is allowed to be made of one the eight materials in the vade mecum (Table 1). In this case, a design is actually a stacking sequence defined by d = [d 1 , d 2 , . . . , d 20 ], where d i identifies the material in layer i.
Then, the optimal stacking sequence is determined by solving the nonlinear, integer programming problem given by Equation (1) by using genetic algorithms (GA). Note that a homogeneous wall made of the material m in the vade mecum represents a possible design d of the current multilayered wall, where all the layers are made of the same material m, i.e., d i = m = constant for i = 1, 2, . . . , N. The best homogeneous wall from the previous section (that made of EPS for all the cases) has been included in the initial population for GA.
The optimal solutions for both locations and the four orientations are shown in Table 2. Not surprisingly, only two materials, among the eight possible ones in the vade mecum, are present in all the optimal solutions: that with the lowest conductivity (level 8 = EPS) and that with the highest effective thermal energy storage capacity (level 6 = RBA-80). The great majority of the layers are made of EPS (16 or 18 for Sauce Viejo and 19 for Frankfurt).
At Sauce Viejo (see Figure 5), despite the little use of RBA-80, the energy performance of the optimal multilayered wall is considerable better than that of the EPS (insulating only) wall: C is reduced between 23.4% and 45.6%, bounds corresponding to S and N, respectively. current multilayered wall, where all the layers are made of the same material , i.e. = = for = 1,2, … , . The best homogeneous wall from the previous section (that made of EPS for all the cases) has been included in the initial population for GA.
The optimal solutions for both locations and the four orientations are shown in Table 2. Not surprisingly, only two materials, among the eight possible ones in the vade mecum, are present in all the optimal solutions: that with the lowest conductivity (level 8 = EPS) and that with the highest effective thermal energy storage capacity (level 6 = RBA-80). The great majority of the layers are made of EPS (16 or 18 for Sauce Viejo and 19 for Frankfurt). At Sauce Viejo (see Figure 5), despite the little use of RBA-80, the energy performance of the optimal multilayered wall is considerable better than that of the EPS (insulating only) wall: is reduced between 23.4% and 45.6%, bounds corresponding to S and N, respectively.  At Frankfurt, as shown in Figure 6, only Layer 11 is made of RBA-80 for all the considered orientations. For this optimal wall, C gain is greatly reduced (from 70.6% for W to 78.9% for N). However, the weight of C gain into C is considerably lower than that of C loss , for which there is a modest improvement: from 0.4% for N to 8.3% for S. At the end, C is reduced between 4.8% for N to 18.8% for S.
In general, for Sauce Viejo (in the Southern hemisphere) as well as for Frankfurt (in the Northern one), the most beneficial effect of RBA-80 is observed at the most sun-exposed wall. Furthermore, the use of RBA-80 leads to considerably higher improvements in Sauce Viejo, where the temperatures remain a longer time in the phase-change range of RBA-80 (from 20 • C to 26 • C, approximately).
Finally, let us remark that Table 2 is the initial version of a new vade mecum of multilayered wall systems for building envelopes. For the time being, it offers a quick answer to the question of minimizing the undesired thermal loads at locations climatically close to those analyzed here. Once again, it is a goal of the recently launched NRG-STORAGE project to continuously enlarge this vade mecum to offer quick solutions for more climates considering different performance indexes.
At Frankfurt, as shown in Figure 6, only Layer 11 is made of RBA-80 for all the considered orientations. For this optimal wall, gain is greatly reduced (from 70.6% for W to 78.9% for N). However, the weight of gain into is considerably lower than that of loss , for which there is a modest improvement: from 0.4% for N to 8.3% for S. At the end, is reduced between 4.8% for N to 18.8% for S. Figure 6. Undesired heat loads in a 20-cm-thick external wall at Frankfurt (Germany) for different orientations: comparison between the EPS insulating wall and the optimal multilayered wall.
In general, for Sauce Viejo (in the Southern hemisphere) as well as for Frankfurt (in the Northern one), the most beneficial effect of RBA-80 is observed at the most sun-exposed wall. Furthermore, the use of RBA-80 leads to considerably higher improvements in Sauce Viejo, where the temperatures remain a longer time in the phase-change range of RBA-80 (from 20 °C to 26 °C, approximately).
Finally, let us remark that Table 2 is the initial version of a new vade mecum of multilayered wall systems for building envelopes. For the time being, it offers a quick answer to the question of minimizing the undesired thermal loads at locations climatically close to those analyzed here. Once again, it is a goal of the recently launched NRG-STORAGE project to continuously enlarge this vade mecum to offer quick solutions for more climates considering different performance indexes.

Conclusions
This article introduces an optimization-based methodology to improve the performance of multilayer walls to be used in building envelopes. This is done in the following steps: i. A vade mecum of materials is built to be used in building envelopes, including particularly materials with thermal energy storage capacity and insulating properties, which should offer a wide choice of material responses. ii.
Given the location and the wall azimuth, EnergyPlus is used to translate the local typical meteorological year (TMY) into time-dependent boundary conditions for the heat conduction equation. iii.
Given a multilayered wall "design", i.e. a specific stacking sequence of materials chosen from the vade mecum of Step i, the heat conduction equation, subject to the boundary conditions from Step ii, is solved along a whole Typical Meteorological Year (TMY) using a robust finite element method to determine the temperature across a given multilayered wall. iv.
The temperature evolution at the internal wall surface, resulting from Step iii, serves to determine the thermal performance index of the current design, here defined by a cost function representing the undesired thermal loads along a TMY. v.
Genetic algorithms are used to make the designs evolve until achieving optimal performance. vi.
Steps ii-v must be repeated first for the remaining orientations at the same location and then for different locations.

Conclusions
This article introduces an optimization-based methodology to improve the performance of multilayer walls to be used in building envelopes. This is done in the following steps: i. A vade mecum of materials is built to be used in building envelopes, including particularly materials with thermal energy storage capacity and insulating properties, which should offer a wide choice of material responses. ii. Given the location and the wall azimuth, EnergyPlus is used to translate the local typical meteorological year (TMY) into time-dependent boundary conditions for the heat conduction equation. iii. Given a multilayered wall "design", i.e., a specific stacking sequence of materials chosen from the vade mecum of Step i, the heat conduction equation, subject to the boundary conditions from Step ii, is solved along a whole Typical Meteorological Year (TMY) using a robust finite element method to determine the temperature across a given multilayered wall. iv. The temperature evolution at the internal wall surface, resulting from Step iii, serves to determine the thermal performance index of the current design, here defined by a cost function representing the undesired thermal loads along a TMY. v. Genetic algorithms are used to make the designs evolve until achieving optimal performance. vi. Steps ii-v must be repeated first for the remaining orientations at the same location and then for different locations.
Here, the methodology was applied to optimize thermal performance of envelopes in Sauce Viejo (Argentina) and Frankfurt (Germany), having humid subtropical and warm temperate climates, respectively. Further, the walls were assumed to face north, east, south, and west. Considering a 20-cm-thick external wall, the optimal solution in any case is mostly made of EPS (the best insulating material in the current vade mecum), including a 1-4-cm-thick layer of a PCM-based composite (the material with the highest thermal energy storage capacity in the current vade mecum). In general, This PCM-based layer is placed next to the middle of the wall, closer to the internal surface. In this way, this methodology defines not only the proper placement but also the proper thickness of the PCM-based layer considering weather and orientation.
Furthermore, despite the little use of PCMs, the undesired heat loads were reduced in comparison to a 20-cm-thick EPS wall: up to 18.8% and 45.6% for Frankfurt and Sauce Viejo, respectively. The better performance at Sauce Viejo is explained by the fact that local temperatures remain for longer periods closer to the phase-change temperatures of the PCM-based composites available in the current vade mecum.
Further steps which follow from this research will address the enrichment of the vade mecum, including particularly PCM-based composites with various phase-change temperature ranges and materials.
Once the vade mecum of materials proves to be wide enough, optimal solutions can be obtained by applying the current methodology within a wide range of climates and different performance indexes, which will be taken as inputs for a new vade mecum of building envelopes.