Modelling the Thermal Energy Storage of Cementitious Mortars Made with PCM-Recycled Brick Aggregates

This paper reports a numerical approach for modelling the thermal behavior and heat accumulation/liberation of sustainable cementitious composites made with Recycled Brick Aggregates (RBAs) employed as carriers for Phase-Change Materials (PCMs). In the framework of the further development of the fixed grid modelling method, classically employed for solving the well-known Stefan problem, an enthalpy-based approach and an apparent calorific capacity method have been proposed and validated. More specifically, the results of an experimental program, following an advanced incorporation and immobilization technique, developed at the Institut für Werkstoffe im Bauwesen for investigating the thermal responses of various combinations of PCM-RBAs, have been considered as the benchmark to calibrate/validate the numerical results. Promising numerical results have been obtained, and temperature simulations showed good agreement with the experimental data of the analyzed mixtures.


Introduction
The energy demand for heating and cooling the global building stock represents a massive part of the total energy consumption around the world (≈ 40%) [1]. In the EU, it accounts for about half of all energy consumption [2]. To attenuate this number, thermal efficiency of construction and building elements, like walls, roofs, and floors, has become the most important measure to enhance energy savings of the new and existing building stock. Recently, with the introduction of the 2019/2021 EU Buildings Directive, it has been additionally emphasized that new buildings of the EU-Member States have to be designed as "Nearly-Zero-Energy Buildings" from the beginning of 2021 onward [3]. For public non-residential buildings, this obligation should become active in 2019 and will be permanently anchored in future building legislation. Moreover, the German Energy Saving Regulation (EnEV) 2014 has proclaimed 25% stricter requirements for energy savings in their 2016 issue, which affects the set of most important regulations of new and existing German buildings. These regulations define the maximum annual primary energy demand and the maximum permissible loss of transmission heat through the building envelope, on the basis of a reference building. In order to

Basic Principles: Thermodynamics and Thermal Energy Storages
The enthalpy H of a thermodynamic system can be defined as with U the internal energy, p the pressure and V the volume of the system. Then, by introducing the first law of thermodynamics, for closed systems and an infinitesimal process, the following equation can be stated as i.e., the variation of the internal system energy (dU) is equal to an infinitesimal amount of heat added (δQ) minus an infinitesimal amount of work performed (δW). Under the hypothesis that only p and V spend work (quasi-static process), δW = p dV. It follows that dU = δQ − p dV.
For a constant pressure, which represents most cases of PCM-concrete systems of construction and building applications, it can be assumed that a small variation of enthalpy is equal to a small amount of heat added: dH = δQ.

Enthalpy Description and Apparent Calorific Capacity Method (ACCM)
In accordance with the fixed grid modelling method and with its discretization through the ACCM approach (see reference [17]), the classical equation for describing a heat conduction problem can written as follows: where Q is the heat of the system, t the time, λ = λ(T, x) the thermal conductivity of the material, which depends on the temperature T and position vector x (of the considered body Ω), .
q v is the possible source term, while ∇. and ∇ are the divergence and gradient tensorial operators. Substituting Equation (5) into Equation (6) leads to the following equation: which is the mostly adopted equation for solving phase changes in simulations for construction and building applications, and it is known in literature as the enthalpy-based method.
Applying the chain rule to ∂H ∂t of Equation (7), ∂T ∂t (8) and by introducing the concept of the Apparent Calorific Capacity Method (ACCM), the following temperature-dependent (apparent of effective) heat capacity expression can be written [17]: = ρC e f f (T) (9) Thus, Equation (7) modifies into the following heat transfer equation .
Finally, the following initial condition (IC) and boundary conditions (BCs) complete the description of the phase change problem: BCs: where Γ T , Γ q , Γ c and Γ r are the essential and natural (flux, convection and radiation) boundary conditions of the domain Ω. T D is the specified temperature imposed at the essential boundary, q the heat flux, h the heat convection and κ the radiation coefficients, respectively, T env the environmental temperature, and, finally, T rad the external radiation source temperature.

Overview of the Experimental Test Data
This section reports the results of an experimental program performed for characterizing the thermal-energy response of PCM-RBA mortar mixtures and their components. The experimental data are used for validation of the numerical simulations discussed in Section 5.
Six mixtures were considered, having a w/c ratio of 0.5 and various amounts of PCM-RBA volume fractions (Table 1). All mixtures were prepared according to EN 196-1 [27]. Recycled bricks (labelled "SB") and high porosity Poroton®fired-clay blocks (labelled "PB") have been considered [28]. These building materials, processed in the form of medium/coarse aggregates and provided by a local company (SHW GmbH, Messel, Germany), were used as carriers (containers) for storing a predefined volume of PCM. RT 25 HC [29] paraffin waxes were used as PCM. They are characterized by a high crystallinity and possess an excellent heat store capacity during phase changes, from solid to liquid and vice versa. The thermo-physical properties of the paraffin wax, considered in this research, are a melting temperature of 25 • C, a storage capacity of approximately 230 kJ/kg, a latent heat capacity of almost 200 kJ/kg, thermal conductivity (in both phases) of 0.20 W/m × K, and densities of 880 (liquid)/770 (solid) kg/m 3 . The PCM-RBAs were thus produced following an advanced encapsulation technique, proposed and patented by the Institut für Werkstoffe im Bauwesen of TU-Darmstadt (Darmstadt, Germany) [30]. The mixture names (labels) in Table 1 aim at providing the key information on the amount of PCM, filling the RBA's open porosity (expressed in volume fraction of the open capillary pore space), and the type of RBAs considered in the mixture. For example, the label "REF-SB" refers to the reference mixture (without PCM) using SB type of bricks; or "SB-65" indicates a mixture using SB-bricks and a filling degree of PCM of 65 V.-% of the total SB-RBA open capillary porosity. The complete description (materials, methods, results and discussion) of the experimental campaign is available in Mankel et al. [28].

DSC Measurements: Aggregates, Paste, and PCMs
Differential Scanning Calorimetry (DSC) tests were performed for each component used in the investigated PCM-RBA mortar mixtures. Their heat storage capacity has been expressed in terms of bulk density times the specific heat capacity, i.e., ρ × C p . Three samples per each component were analyzed and the mean value for each of them was presented in this section. Particularly, the DSC thermograms under either heating or cooling, for the cement paste with a w/c of 0.50 and for both SBand PB-RBAs, have been shown in Figure 1a. They were done within the temperature range of 10 • C to 40 • C and using a heating/cooling rate of 10 K × min −1 . From these results, it can be observed that the sensible behavior of both RBAs is almost similar, whereby a slightly higher sensible heat storage capacity can be detected for the cement paste.  In Figure 1b, the DSC results of the used paraffin wax (Rubitherm® RT25HC, Rubitherm Technologies GmbH, Berlin, Germany) is shown. The test procedure was conducted in accordance with the IEA DSC 4229 PCM Standard [31] to determine the final heating/cooling rate for the considered dynamic DSC measurements. The adopted heating/cooling rate was 0.25 K × min −1 on the final measured results. This value represented a compromise between accuracy of the results and acceptable mitigation of heating rate measurements, fulfilling the IEA DSC 4229 requirements [31].
DSC curves of the considered PCM, in Figure 1b, shows a sensible heat storage character of the material in those temperature ranges, far from the phase change responses (e.g., in the solid and liquid stages), and shows a pronounced latent peak in the region close to the temperature where the phase change occurs (Tm = 24.5 °C for heating and Tm = 22.95 °C for cooling).  In Figure 1b, the DSC results of the used paraffin wax (Rubitherm®RT25HC, Rubitherm Technologies GmbH, Berlin, Germany) is shown. The test procedure was conducted in accordance with the IEA DSC 4229 PCM Standard [31] to determine the final heating/cooling rate for the considered dynamic DSC measurements. The adopted heating/cooling rate was 0.25 K × min −1 on the final measured results. This value represented a compromise between accuracy of the results and acceptable mitigation of heating rate measurements, fulfilling the IEA DSC 4229 requirements [31].  Figure 1b, shows a sensible heat storage character of the material in those temperature ranges, far from the phase change responses (e.g., in the solid and liquid stages), and shows a pronounced latent peak in the region close to the temperature where the phase change occurs (T m = 24.5 • C for heating and T m = 22.95 • C for cooling).

Thermal Conductivity of the PCM-RBA Mortars
The thermal conductivity of PCM-RBA mortar mixtures was determined using the Hot-Disk transient plane source method [28]. For this aim, three samples of 150 mm × 150 mm × 80 mm cuboids were tested. The measurements were done with a 9.9 mm diameter sensor, in three different specimen sides. Steady-state conditions with a temperature of 20 • C were considered. In Figure 2, it can be observed that all mixtures have quite comparable thermal conductivities, which range between 0.696 (min. value) and 0.846 (max. value) W/mK. As a general trend, it can be observed that the mixtures with PB RBAs deal with slightly lower conductivities than the SB RBA ones. This can be attributed to the higher porosity of the SB RBAs, which affect the overall conductivity of the composites. It could also be observed that for higher PCM contents, the thermal conductivity is lower.
In Figure 1b, the DSC results of the used paraffin wax (Rubitherm® RT25HC, Rubitherm Technologies GmbH, Berlin, Germany) is shown. The test procedure was conducted in accordance with the IEA DSC 4229 PCM Standard [31] to determine the final heating/cooling rate for the considered dynamic DSC measurements. The adopted heating/cooling rate was 0.25 K × min −1 on the final measured results. This value represented a compromise between accuracy of the results and acceptable mitigation of heating rate measurements, fulfilling the IEA DSC 4229 requirements [31].
DSC curves of the considered PCM, in Figure 1b, shows a sensible heat storage character of the material in those temperature ranges, far from the phase change responses (e.g., in the solid and liquid stages), and shows a pronounced latent peak in the region close to the temperature where the phase change occurs (Tm = 24.5 °C for heating and Tm = 22.95 °C for cooling).

Thermal Conductivity of the PCM-RBA Mortars
The thermal conductivity of PCM-RBA mortar mixtures was determined using the Hot-Disk transient plane source method [28]. For this aim, three samples of 150 mm × 150 mm × 80 mm cuboids were tested. The measurements were done with a 9.9 mm diameter sensor, in three different specimen sides. Steady-state conditions with a temperature of 20 °C were considered. In Figure 2, it can be observed that all mixtures have quite comparable thermal conductivities, which range between 0.696 (min. value) and 0.846 (max. value) W/mK. As a general trend, it can be observed that the mixtures with PB RBAs deal with slightly lower conductivities than the SB RBA ones. This can be attributed to the higher porosity of the SB RBAs, which affect the overall conductivity of the composites. It could also be observed that for higher PCM contents, the thermal conductivity is lower. Moreover, it can be also observed that the thermal conductivities are only slightly affected by the different PCM volume fractions. A little increase of the thermal conductivity was measured for those mixtures with a higher pore filling ratio (i.e. PB-80 in comparison with PB-65). This could be the result of the allocation of paraffin wax into the capillary pore space of the composites, which leads to a slight increase in heat conductivity of the PCM-RBAs since the considered PCMs are more conductive than air. However, it can finally be concluded that immobilizing PCMs into the RBA porous structure does basically not (e.g. SB-80 vs. SB-65), or only slightly (see, PB-80 vs. PB-65), modify the overall thermal conductivity of RBA mortars.

Thermal-Energy Storage DKK Tests in Spherical Samples
Spherical-shaped specimens were used to monitor the time-dependent temperature evolution of the PCM-RBA mortars. The adopted and patented non-conventional testing technique (namely Dynamische Kugel Kalorimetrie), DKK [26]) was followed by the authors for the TES identification of the composite materials under investigation. For each considered mixture, three spherical samples were produced, and two thermocouples were positioned in the center of each sphere and at its outer surface, respectively. Heating tests were done by using an isothermal conditioned oven with a fixed temperature of 49 • C. Cooling tests were done with a climatic chamber fixing the temperature at ca. 9 • C. The graphs of Figure 3 show the average results (from three independent spherical specimens) of the measured temperature evolution, in the center of the spherical samples, versus time. Heating and cooling results are plotted for both PCM SB-brick and PB-brick mixtures. It can be seen that, for the mortar mixtures casted with either SB or PB bricks, a delay of the temperature development takes place when PCM-RBAs, with PCM filling degrees of 65 V.-% and 80 V.-%, are analyzed. The presence of PCM and their melting/solidification behavior actually shifts the temperature curve into the right/down direction for heating response (Figure 3a,c) and up/left for cooling (Figure 3b,d). Particularly, a quasi-horizontal plateau of the temperature evolution data can be appreciated during both temperature rise (heating) or temperature decrease (cooling), i.e. in the range between 21 and 26 • C the PCM and vice versa. By taking into consideration the effect of PCM volume fractions, it appeared that almost no thermal differences between the considered RBAs and PCM additions, i.e., 65 and 80%, exist. Comparison of the temperature evolutions SB-65 with SB-80 or PB-65 with PB-80 show an almost equal response. For each mixture and/or sample, the complete phase change of the paraffin was always fully occurring, meaning that all PCMs were in their final state at the end of each heating or cooling test. This is also shown in Figure 3 where the center of the spheres at the end of each cooling test has a temperature less than 10 • C and under heating close to 50 • C, while the melting points range between 21-26 • C. For a thorough discussion on the present experimental data, with emphasis on the DSC analyses, conductivity measurements and DKK tests, reference is made to Mankel et al. 2019 [28].

4.1.D Spherical-based Solution for the ACCM
Equation (10) is now transferred into spherical coordinates for predicting the thermal energy storage behavior and temperature evolution in the tested specimens presented in Section 3.3.
From this, the following relation can be derived: For a thorough discussion on the present experimental data, with emphasis on the DSC analyses, conductivity measurements and DKK tests, reference is made to Mankel et al. 2019 [28].

D Spherical-Based Solution for the ACCM
Equation (10) is now transferred into spherical coordinates for predicting the thermal energy storage behavior and temperature evolution in the tested specimens presented in Section 3.3.
From this, the following relation can be derived : where (r, θ, φ) are the radial distance and polar and azimuthal angles. Making use of the spherical symmetry, the previous relationship can be significantly simplified as follows:

Schematization and Discretization
The ACCM model, under the assumptions of spherical geometry and symmetry, was solved by means of the finite difference method. The heat-diffusion through the PCM-RBA mortar systems was calculated by solving the differential equation previously described in Equation (14), having a 1D spherical-based hypothesis and by adopting heat convection (Robin) boundary conditions for describing the environmental surface conditions of either a furnace or a climate chamber.
Thus, the boundary condition for the sample core at node "1" (Figure 4) was adiabatic, meaning that the heat flux, q 1 , is null. Thus, the boundary condition for the sample core at node "1" (Figure 4) was adiabatic, meaning that the heat flux, q1, is null.  Sample core: while the Robin boundary condition is applied at the outer surface (node "ns") Ambient condition: where Tf is the ambient temperature (fixed by the oven or climatic chamber), Tr is the surface temperature at r=R, while h is the heat transfer coefficient. In this context, the Finite Difference (FD) space domain was discretized into ns spaces (namely Finite Differences), leading to ns+1 nodes and nt time discretization steps. Furthermore, a fully implicit Sample core : Materials 2020, 13, 1064 9 of 16 while the Robin boundary condition is applied at the outer surface (node "n s ") Ambient condition : where T f is the ambient temperature (fixed by the oven or climatic chamber), T r is the surface temperature at r = R, while h is the heat transfer coefficient. In this context, the Finite Difference (FD) space domain was discretized into n s spaces (namely Finite Differences), leading to n s + 1 nodes and n t time discretization steps. Furthermore, a fully implicit Euler Method for the transient problem was applied.
The differential equation shown in Equation (14), within a space domain with a length of R between the node 1 to n s (see Figure 4a), can be solved as shown in the following scheme: By using the implicit backward Euler Method for time (j) and central for space (i) (see Figure 4b), Equation (17) can be discretized to The boundary condition, for the sample core at r = 0 (i = 1), has been further developed. The right-hand side of Equation (17) can be simplified as follows [32] lim Then, by adopting the central difference approximation of the adiabatic boundary condition expressed of Equation (15) and by using one ghost node (namely node "0"), the following expression can be stated: It follows that T j+1 2 = T j+1 0 and combining it into Equations (18) and (19), it can be easily achieved the following expression of the adiabatic boundary condition: For the implementation of the Robin boundary condition at the outer surface of the specimens, e.g. node r = R (i = n s ), Equation (15) can be discretized in the following way: In Equation (22), the temperature at the ghost node r = R + 1 (i = n s + 1) is known, T j+1 n s +1 = T f , thus the following expression can be rewritten: Thus, Equation (16) (24) and after some mathematical elaborations, the following expression can be achieved for the Robin natural boundary condition:

Numerical Results and Comparisons
This section reports the description of the numerical results and shows their comparisons against the experimental data, which has been briefly outlined in Section 3. The numerical simulations are based on the assumption that the PCM-RBA mortars can be considered as a continuum and homogenous media. In this context, homogenized (meso-scale based) parameters were considered in the selection of the input data for the numerical prediction.

Homogenized Macroscopic C eff Model for the PCM-RBA Mixtures
A homogenization technique was employed for evaluating the effective specific heat capacity C eff of the PCM-RBA mortars. It is based on the mixture theory by using the volume percentages of each individual component such as RBAs, cement paste and PCMs. More precisely, the model smears out the specific heat capacity of the RBA C RBA (T), cement paste C paste (T) and the apparent specific heat capacity of the PCM C app,PCM (T) through adopting the volume fraction of each component as the smeared out (weighting) factor χ.
The specific heat capacities of each component were experimentally determined with DSC measurements (see Section 3.1) and are shortly summarized in Table 2.   RBAs and cement pastes have only the sensible heat storage part (Figure 1a), while the PCMs have an apparent specific heat capacity C app,PCM (T) that incorporates the additional latent behavior within the temperature range of the phase change (during melting and solidification, as shown in Figure 1b).
The evaluation of the specific heat capacities was determined in two consecutive steps. First, at aggregate level, where the PCM-RBAs were considered as lumped components of RBAs (C RBA (T)) plus PCM (C app,PCM (T)) and weighting their volume fractions χ to achieve the smeared C eff,PCM-RBA (T):

PCM-RBAs
: C e f f ,PCM−RBA (T) = χ RBA × C RBA (T) + χ PCM × C app,PCM (T) (26) where χ RBA and χ PCM are the volume fraction of the recycled bricks and the filled PCMs, respectively. Then, at a composite level (i.e., PCM-RBA mortar) a homogenized overall system of the effective specific heat capacity, C eff , was determined by weighting the heat capacities C eff,PCM-RBA (θ), evaluated through Equation (26), and C paste (θ) of the individual material components by their volume fractions ψ: The exact volume fraction ratio between PCM-RBAs and cement paste, of the investigated mixtures, were investigated by performing µ3D-XCT-scans [33] of the spherical specimens. Nine 2D-slices were extracted from each 3D body of the scanned spherical specimens ( Figure 5) and mean volume fractions of PCM-RBAs and cement paste were accurately determined by image analyses. An average area ratio was evaluated for each slice by applying a recoloring of the surface area through white and red pictures ( Figure 5). The results of these analyses are shown in Figure 6 where the volume fractions of paste and PCM-RBA were determined in both PCM.RBA mortars SB and PB.
The effective specific heat capacity of the composite systems, C eff (T), were thus modelled using Equations (26) and (27) by taking into account the determined volume fractions between cement paste (including air voids) and PCM-RBAs (shown in Figure 6) and considering the specific heat capacities of each component as summarized in Table 2.   The evaluation of the specific heat capacities was determined in two consecutive steps. First, at aggregate level, where the PCM-RBAs were considered as lumped components of RBAs (CRBA(T)) plus PCM (Capp,PCM(T)) and weighting their volume fractions χ to achieve the smeared Ceff,PCM-RBA(T): where χRBA and χPCM are the volume fraction of the recycled bricks and the filled PCMs, respectively. Then, at a composite level (i.e., PCM-RBA mortar) a homogenized overall system of the effective specific heat capacity, Ceff, was determined by weighting the heat capacities Ceff,PCM-RBA(θ), evaluated through Equation (26), and Cpaste(θ) of the individual material components by their volume fractions ψ: The exact volume fraction ratio between PCM-RBAs and cement paste, of the investigated mixtures, were investigated by performing 3D-XCT-scans [33] of the spherical specimens. Nine 2Dslices were extracted from each 3D body of the scanned spherical specimens ( Figure 5) and mean volume fractions of PCM-RBAs and cement paste were accurately determined by image analyses. An average area ratio was evaluated for each slice by applying a recoloring of the surface area through white and red pictures ( Figure 5). The results of these analyses are shown in Figure 6 where the volume fractions of paste and PCM-RBA were determined in both PCM.RBA mortars SB and PB. The effective specific heat capacity of the composite systems, Ceff(T), were thus modelled using Equations (26) and (27) by taking into account the determined volume fractions between cement paste (including air voids) and PCM-RBAs (shown in Figure 6) and considering the specific heat capacities of each component as summarized in Table 2.    of each component as summarized in Table 2.   of each component as summarized in Table 2.

Numerical Prediction and Comparison
By implementing the ACCM procedure described in Sections 2 and 4, temperature evolutions were simulated and compared with the experimental data reported in Section 3. More in detail, the spherical samples made of RBA mortars (with and without PCMs) were simulated with the proposed heat flow model. The input values were obtained from the conducted experimental measurements, as well as from the C eff (T) curves described in Section 5.1. All input parameters employed in the aforementioned simulations are summarized in Tables 3 and 4. Table 3. Overview of the numerical parameters assumed for the SB mixtures. The thermal conductivities were assumed as temperature-independent and measured through Hot-Disk tests, as outlined in Section 3.2. Then, the effective specific heat capacity C eff (T) was modelled by using the homogenized model as described before in this section. The number of FD space discretization was chosen 100 while the number of time steps selected was 1000, in all simulations. Moreover, the calibrated Robin heat transfer coefficient (h), representing the heat transfer conditioning coefficient was the same for each mixture and reported in Tables 3 and 4.

Numerical Parameters
The numerical simulations, based on the input parameters above mentioned, have been compared against the experimental data of Section 3.3. In Figure 9, the six graphs show the temperature evolutions for the control mixtures (Figure 9a,b), and for that one having PCM: i.e., SB-65 (Figure 9c It can be observed that the modeling approach was able to simulate the experimental temperature evolutions very accurately. In particular, the simulations of the reference mixtures represented in Figure 9a,b show a very good agreement with the experimental results. Then, the simulations of the mixtures TES enhanced with PCM also show good comparisons and trends as the experimental values. A slightly overestimated latent effect can be observed, which results in a slightly amplified shoulder in the temperature evolution (see Figure 9c-f). A reason for this effect can be attributed to micro-and/or meso-structural effect, which can influence the effective thermal conductivity and the melting/solidification activations of the integrated PCMs. This could lead to a slight deviation of the latent effects.
In this context, it may be important to remark that almost all input parameters were chosen from the characterizations of the thermal tests as well as from the homogenized C eff model, which is actually represented by the experimentally determined heat capacities of the individual material components weighted by their volume fractions (see Section 5.1). Thus, with this set of input parameters, a sound numerical prediction, for all six PCM-RBA mortar systems, could be achieved, without the need of applying re-calibration and/or optimizations. The results also show that, with a unique heat transfer coefficient, described as the h (Robin) parameter, the simulations for the reference PCM-RBA mortar systems, REF-SB and REF-PB, are almost in perfect agreement with the corresponding experiments. This supports the assumption that both the spherical symmetry and the hypothesis that the composite PCM-RBA mortar systems could be considered by a homogeneous medium, were both effectively correct. In this context, it may be important to remark that almost all input parameters were chosen from the characterizations of the thermal tests as well as from the homogenized Ceff model, which is actually represented by the experimentally determined heat capacities of the individual material components weighted by their volume fractions (see section 5.1). Thus, with this set of input parameters, a sound numerical prediction, for all six PCM-RBA mortar systems, could be achieved, without the need of applying re-calibration and/or optimizations. The results also show that, with a unique heat transfer coefficient, described as the h (Robin) parameter, the simulations for the reference PCM-RBA mortar systems, REF-SB and REF-PB, are almost in perfect agreement with the corresponding experiments. This supports the assumption that both the spherical symmetry and the hypothesis that the composite PCM-RBA mortar systems could be considered by a homogeneous medium, were both effectively correct.

Conclusions
Based on the results shown in this paper, the following conclusions can be drawn: • An enthalpy-based model, formulated for spherical coordinates and symmetry, was proposed for predicting the thermal energy storage in the tested PCM-RBA specimens. • Thermal measurements, obtained from dynamic DSC and steady-state Hot Disk tests, were employed for calibrating the model of the numerical activities.
• A mixture theory, based on volume fractions deduced from 3D micro µX-ray computer tomography measurements, was used for generating the resulting meso-composite thermal parameters adopted in the numerical analysis. Particularly the C eff (T) curves have been based on this approach.

•
The numerical simulations for the temperature evolution, compared with the experimental DKK results, showed accurate and consistent agreement. It may be important to highlight that these numerical results were just based on input parameters obtained from the experimental characterizations of the mortar components.

•
No fitting adjustments, re-calibrations, or numerical adaptions were necessary for reaching good agreement between numerical to experimental data. This can confirm that on the one hand the experimental activities were performed in an accurate way, and on the other hand, that the numerical assumptions and procedures are very accurate to model TES responses in cementitious materials.
Future numerical developments, which follow this research, will include micro-to-meso scale analysis taking into considerations local effects like inclusions (PCM-RBAs), air bubbles, porosity, and interface effects. These further steps will lead to optimizing the "best" recipe for achieving the most performing energy-saving and sustainable cementitious composite.