Next Article in Journal
New Insights into Hybrid Materials Based on Conductive Polymers and Their Use in Energy-Related Applications
Next Article in Special Issue
Eutectic Fatty Acids Phase Change Materials Improved with Expanded Graphite
Previous Article in Journal
Numerical Study on Hydrodynamic Characteristics and Electrochemical Performance of Alkaline Water Electrolyzer by Micro-Nano Surface Electrode
Previous Article in Special Issue
Development of a Kinetic Model for the Redox Reactions of Co2.4Ni0.6O4 and SiO2/Co2.4Ni0.6O4 Oxides for Thermochemical Energy Storage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermomechanical Performance Analysis of Novel Cement-Based Building Envelopes with Enhanced Passive Insulation Properties

by
Jorge Marin-Montin
1,
Eduardo Roque
1,
Yading Xu
2,
Branko Šavija
2,
Juan Carlos Serrano-Ruiz
1 and
Francisco Montero-Chacón
1,*
1
Materials and Sustainability Group, Department of Engineering, Universidad Loyola Andalucía, Avenida de las Universidades, s/n, 41704 Sevilla, Spain
2
Microlab, Faculty of Civil Engineering and Geosciences, Delft University of Technology, 2628 CN Delft, The Netherlands
*
Author to whom correspondence should be addressed.
Materials 2022, 15(14), 4925; https://doi.org/10.3390/ma15144925
Submission received: 9 June 2022 / Revised: 11 July 2022 / Accepted: 13 July 2022 / Published: 15 July 2022
(This article belongs to the Special Issue Advanced Materials for Thermal Energy Storage)

Abstract

:
The design of new insulating envelopes is a direct route towards energy efficient buildings. The combinations of novel materials, such as phase-change (PCM), and advanced manufacturing techniques, such as additive manufacturing, may harness important changes in the designing of building envelopes. In this work we propose a novel methodology for the design of cement-based building envelopes. Namely, we combined the use of a multiscale, multiphysical simulation framework with advanced synthesis techniques, such as the use of phase-change materials and additive manufacturing for the design of concrete envelopes with enhanced insulation properties. At the material scale, microencapsulated PCMs are added to a cementitious matrix to increase heat storage. Next, at the component level, we create novel designs for the blocks, here defined as HEXCEM, by means of additive manufacturing. The material and component design process is strongly supported on heat transfer simulations with the use of the finite element method. Effective thermal properties of the mixes can be obtained and subsequently used in macroscale simulations to account for the effect of the volume fraction of PCMs. From the experimental and numerical tests, we report an increase in the the thermal inertia, which results in thermal comfort indoors.

1. Introduction

The building sector is responsible for 40% of the primary energy requirements in the EU and 36% of greenhouse gas emissions in Europe [1,2]. Therefore, there is a global need for improving energy management in buildings, which can be approached using different strategies [3] dealing with generation, energy storage, materials, or control systems that not only are more sustainable and environmentally-friendly, but also economical. One solution for providing thermal comfort while reducing the energy demand is the use of new insulating building envelopes [4,5].
Active/passive thermal energy management systems [6,7] based on renewable sources (e.g., solar radiation) are gaining much attention within the nearly zero energy buildings (NZEBs) paradigm [8]. While active systems are generally more difficult to implement in a building, passive systems can be easily integrated in building envelopes [9], with neither aesthetic nor functional drawbacks [7].
Concrete is one the most widely used materials worldwide [10,11]. It is mainly used in buildings for structural elements (i.e., columns and beams), foundations, slabs, and envelopes. Although concrete has caught attention for its mechanical properties, it can also be designed with targeted thermal insulation properties. For instance, higher internal porosity yields lower effective thermal conductivity [12]. Of course, this is a microstructural feature, but the same relationship applies when dealing with components. In this sense, if we place air gaps within a concrete block, the effective thermal conductivity of the component drops [13].
Another way to increase concrete’s thermal insulation properties is by means of so-called phase-change materials (PCMs). PCMs can absorb energy in the form of latent heat to complete a reversible solid-to-liquid phase transition [14,15]. In eutectic PCMs, this process takes place at a constant temperature. PCMs can be mixed with cementitious materials [16,17,18,19,20], and the resulting composite material has an increased heat capacity and better insulation properties, which are suitable for building envelopes [21,22]. In such cases, when the outdoor temperature increases, the material absorbs the heat, and if the melting point is surpassed, part of this energy is stored in the form of latent heat. This slows down the global heating process and the other way around.
PCMs themselves have low thermal conductivity [23]. Therefore, when embedded in the cement paste, the resulting thermal conductivity of the mix decreases [24,25]. This combination of low conductivity and high thermal energy storage capacity is very interesting in temperature management of buildings [26,27].
The integration of PCMs in cement-based materials is generally accomplished by encapsulation [28]. In the case of macrocapsules, they are generally placed in the mold before the casting [26] or mixed with the paste [14]. With the former method, a larger mass of PCM can be added and the risk of leakage is lower than that of the latter. One important drawback of micro-encapsulated PCMs is that it lowers the mechanical strength of the mortar in a drastic way [24]. However, the surface-area–volume ratio is much higher in microcapsules, which enhances the heat transfer process. Macrocapsules also may present partial phase transitions due to the low conductivity of the PCM and the larger sizes, which results in low performance of the energy storage feature [14].
Most of the works in the literature tackle the insulation design of building envelopes by means of PCMs [29]. These include experimental [30,31,32,33] and numerical approaches [34,35,36,37], mainly based on the use finite element method (FEM). Moreover, the FEM becomes a powerful method when dealing with coupled physical phenomena, such as heat transfer and solid mechanics [38,39]. However, there the coupled effect between PCMs and air voids have not been sufficiently covered. Wu et al. [40] presented a study on cellular cementitious composites with microencapsulated PCMs from a mechanical perspective. In this work, we focus on its thermomechanical behavior instead.
The design of building envelopes with PCMs and air cells is a multiscale problem. There exist different multiscale approaches [41] for solving thermomechanical problems. Homogenization-based schemes either by averaged properties [42] or so-called FE2 techniques [43,44,45] have been successfully used in these applications.
In this work we present a novel methodology for the design of cement-based building envelopes that combines the use of a multiscale, multiphysical simulation framework with advanced synthesis techniques, such as the use of phase-change materials and additive manufacturing. To improve the insulation properties, we tackled the design process at the material level, by means of microencapsulated PCMs, and the component level, via air cells. Moreover, additive printing (also commonly referred to as 3D printing) was used to prepare tailor-made specimens that were already analyzed numerically. Thus, we address the design process of building envelopes by combining advanced numerical and experimental techniques in a seamlessly fashion. In the first place, we designed the mix of the material and validated it numerically. Then, we used a homogenization approach to obtain macroscopic material properties that are used in FEM simulations of different envelope configurations. The numerical results are in perfect agreement with the experiments, fulfilling the validation process of the design and also validating a numerical framework to explore other designs that were not accounted for within the experimental campaign.
Although the design methodology presented in this paper focuses on building envelopes, it can be used without loss of generality in the design processes of other type of devices, such as thermal batteries [46], geothermal dissipation devices [47], freeze–thaw damage-resistant devices [48], or early-age cracking mitigation devices [16].

2. Materials and Methods

In this section we present the multiscale methodology followed for the design and validation of the building envelopes. In the first place, we designed a cementitious matrix with microencapsulated PCM at different volume fractions. These mixes were then characterized numerically using the FEM. The resulting material properties were subsequently upscaled. Then, we carried out the component design by considering different cell sizes. We used the FEM at the component scale for two purposes: (i) to validate the designs and (ii) to develop a multiscale simulation framework for new building envelopes.

2.1. Material Design

The material designed for the building envelopes consists of a CEM I 42.5 cementitious matrix (ENCI, The Netherlands) and microencapsulated PCM (Encapsys LLC, United States) with different volume fractions: V f = 0%, 10%, and 30%. The mix had a water-to-cement ratio of w/c = 0.45 and PCM capsules reach a random distribution within the cement paste, as shown in Figure 1. The PCMs used are organic paraffin encapsulated by a melamine-formaldehyde shell. The same PCMs have been used by the authors in reference [49]. Further details on the PCM characterization can be found in references [49,50]. Table 1 shows the composition of the three mixes.
Once ready, the mixes were poured into the molds of the bricks and hardened for 2 days. Then, these were demolded and placed in a cure room until 21 or 28 days. According to [18], the thermal conductivity of cement paste decreases during the first 7 days of curing, but remains practically constant afterwards.
The PCMs used are organic paraffin encapsulated by a melamine-formaldehyde shell, as in reference [16]. The heat of fusion, h f , is 143.5 J/g and the melting point, T f , is 19 °C [49]. The thermophysical properties of the constituents are shown in Table 2.

2.1.1. Unit Cell FEM Model

In our modeling approach, we considered 2D unit cells of the material including the PCM and the cementitious matrix, which were explicitly modeled at different volume fractions. For the FEM discretization, we considered isoparametric quadrilateral elements. Figure 2 shows three different meshes. The material properties considered for each phase were those from Table 2.
Heat transfer is described as an elliptical partial differential equation that reads:
ρ c p T t · ( k T ) = q
where ρ is the density, c p the specific heat, T is the temperature field, t is time, k is the thermal conductivity, and q is the volumetric heat source.
In this work, we account for the non-linear effect of the phase transition following the effective specific heat approach [51,52]:
c p ( T ) = c p , s for T < T f Δ T f 2 c p , s + h f Δ T f for T f Δ T f 2 T T f + Δ T f 2 c p , l for T > T f Δ T f 2
where c p , s and c p , l are the specific heat in the solid and liquid phases, respectively, h f is the enthalphy of fusion, T f is the phase-change temperature, and Δ T f the temperature window.
In order to solve an initial boundary value problem (IBVP), apart from the partial differential equations involved, we need to define the solution domain (i.e., the unit cell and time) and initial and boundary conditions. In our simulations we applied Dirichlet boundary conditions on the left side of the unit cell and evaluated the temperature evolution at different locations (Figure 3).
The heat transfer problem can be discretized using the FEM [53], yielding the following matrix form:
C t h T ˙ + K t h T = q
where C t h and K t h are the heat capacity and thermal conductivity matrices, respectively, T is the temperature vector, and q is the flux vector.
This is a time-dependent problem; therefore, we can discretize in time using Equation (3) with a time increment Δ t , such that at an instant of time t j , the system reads:
C t h T ˙ j + K t h T j = q j
We make use of the Crank–Nicolson time integration scheme [53] to determine the temperature at an instant of time t j + 1 :
T j + 1 = 1 Δ t C t h + 1 2 K t h 1 1 Δ t C t h 1 2 K t h T j + q
Since this problem is non-linear, we must update the capacity matrix at each time step, i.e., C t h = C t h ( T j ) , which gathers the phase-change phenomenon, as described above.

2.1.2. Upscaling Procedure

The unit cell model is heterogeneous, since it accounts for two different phases. However, for the component design, we require material properties of the mix. For this purpose, we use an upscaling homogenization approach [54] for determining the macroscopic thermophysical properties, namely, k and ρ c p .
The effective conductivity, k e f f , of a two-phase system with spherical inclusions can be obtained as [55]:
k e f f = 2 k m ( 1 ϕ ) + ( 1 + 2 ϕ ) k p ( 2 + ϕ ) + ( 1 ϕ ) k p k m
where k m and k p are the thermal conductivity of the matrix and particle, respectively, and ϕ is the volume fraction of particles. In our case, V f = ϕ , so we keep using ϕ for the mathematical description of the problem in the following.
The effective volumetric heat capacity, ρ c p e f f , is evaluated using mixing rule [52]:
ρ c p e f f ( T ) = ϕ ( ρ p c p p ( T ) ) + ( 1 ϕ ) ( ρ m c p m ( T ) )
where ρ m and ρ p are the densities of the matrix and particle, respectively. Similarly, c p m and c p p are the specific heats of the matrix and particle, respectively.
Using the definitions above, it is possible to set the homogenized phase-change (i.e., temperature-dependent) model as:
ρ c p e f f ( T ) = ρ c p , s e f f for T < T f Δ T f 2 ρ c p , s e f f + ϕ ρ p h f Δ T f for T f Δ T f 2 T T f + Δ T f 2 ρ c p , l e f f for T > T f Δ T f 2
where ρ c p , s e f f and ρ c p , l e f f are the effective volumetric heat capacity of the solid and liquid phases, respectively, ρ p is the density of the particle, and h f and T f are the heat of fusion and melting point of the PCM.

2.2. Component Design

The new building envelope proposed herein is made of a PCM-based cementitious matrix, described previously, with air cells, yielding a honeycomb-like structure. In the first place, tailor-made molds were fabricated using 3D printing technology. Then, different brick configurations were cast and tested experimentally using a hot plate setup. Finally, the designs were validated using a macroscopic FEM model. Moreover, the simulation framework was also validated experimentally so that it can be used in the analysis of different designs.

2.2.1. Experimental Methods

The building envelopes follow a honeycomb structure so that large amounts of air could be hosted within the blocks while preserving a compact and robust structure. Due to the complex geometry of the blocks, new tailor-made molds were fabricated using 3D printing technology [56] (Ultimaker 2 + 3D) with acylonitrile butadiene styrene (ABS) and polylactic acid (PLA). Namely, we used ABS for the large air cells design, and we used PLA for the other two. The advantage of PLA over ABS is that the former is biodegradable and allows a higher print speed [57]. In any case, the choice of the material does not affect the results of this research. Three blocks with size 100 × 52 × 20 mm3 and different air cell sizes were designed using CAD techniques (Figure 4).
Once the bricks were printed, these were filled with a mix of two components of PS8510 silicon rubber, which takes two hours to harden. In such a way, molds of the printed bricks were realized in a flexible material, which is better for the demolding process (Figure 5a). Another option is to directly print the molds in 3D, but this may result in lower precision [58].
Finally, three types of honeycomb blocks (i.e., small, medium, and large cells), as shown in Figure 5b, were cast using the material mix described previously. Ordinary cement paste (CEM I 42.5, w/c = 0.45) included encapsulated organic paraffin wax with different volume fractions ( V f = 0%, 10%, and 30%). Table 3 shows the different configurations considered in our analysis.
To analyze the thermal behavior of the new blocks, a hot plate setup (Figure 6) was used along with a FLIR A320 thermal camera and ThermaCAM Researcher Pro software. All the specimens were insulated on the free boundaries (i.e., back, top, and sides) except for the front, which was left open so as to evaluate the temperature evolution with the thermal camera. The insulation was made of styrofoam, sponges, and duct tape. The hot plate was set to 100 °C, and the room temperature was 22.5 ± 1 °C. The tests had a total duration of 3.5 h.
The evolution of the temperature was obtained for three sample heights, namely, 30%, 60%, and 90%, to see the spatial gradient. Temperature was evaluated every 9 min until the steady-state conditions were accomplished.

2.2.2. Brick FEM Model

At the component level, we used the homogenized phase-change FEM formulation presented in Section 2.1.2 for the non-linear unsteady heat transfer problem. For the FEM model, we used isoparametric quadrilateral elements, and the boundary conditions reflected those of the hot plate setup.
As shown in Figure 7, the contact of the bottom side of the brick with the hot plate was not perfect, and internal convective losses due to the surface roughness appeared. In order to model this, we used thermal contact resistance elements (i.e., interfacial conductance) that were calibrated with part of the experimental results. The temperature of the source was set to that of the hot plate, 100 °C. Air trapped in the cells would also generate convective losses that would affect the heat transfer process in the blocks. These were taken into account in our model by using film conditions:
q b o t t o m , j = h c A b o t t o m T p l a t e T b o t t o m , j q c e l l s , j = h a i r A c e l l s T a i r T c e l l s , j
where q b o t t o m , j is the imposed flux at the bottom nodes and instant t j , h c is the interfacial conductance, A b o t t o m is the contact surface, T p l a t e is the temperature of the hot plate, T b o t t o m , j is the temperature of the bottom surface of the brick, q c e l l s , j is the imposed flux on the inner surface of the cells, A c e l l s is the inner surface of the cells, T a i r is the air temperature at the cells, and T c e l l s , j is the temperature at the inner surface of the cells.
The time evolution of the temperature distribution was analyzed by placing three probes in the domain at the center of the horizontal axis, and three different heights: 30%, 60%, and 90%—those recorded in the experimental setup.
We extended our analysis of the component so as to evaluate the mechanical stresses promoted by the temperature gradients. Thus, according to the constitutive equation in solid mechanics:
σ = D : ε e l
where σ is the stress tensor, D the consitutive tensor, and ε the elastic strain tensor.
The elastic strain tensor can be evaluated as:
ε e l = ε ε t h
where ε is the total strain obtained from the displacement field gradient, and ε t h is the thermal strain tensor, which in isotropic solids is:
ε t h = α Δ T I
where α is the thermal expansion coefficient; Δ T is the temperature jump defined as the difference between the temperature field obtained from the heat transfer analysis, T, and the reference temperature, T r , i.e., Δ T = T T r ; and I is the identity tensor.
Regarding the FEM model, we used isoparametric quadrilateral elements for the space discretization, and the system was solved at every instant of time t j :
K u j = f j
where K is the stiffness matrix, u j the displacement vector, and f j the force vector, which is defined as:
f j = f m , j + f t h , j
with f m , j being the mechanical force vector and f t h , j the thermal force vector.
For an element el and isotropic expansion, the element thermal force vector, f t h , j e l , reads:
f t h , j e l = α Δ T B T D 1 1 0 T
B is the matrix containing first-derivatives of the shape functions, and D is the constitutive matrix (in our case, plane stress problem).
Finally, we solved the mechanical problem for each t j using the preconditioned conjugate gradient method.

3. Results and Discussion

3.1. Material Thermal Characterization

In the first place, we performed heat transfer transient simulations at the unit cell scale to characterize the thermal behavior of the two materials designs. For this purpose, we defined three unit cells with the following volume fractions of PCM: V f = 0 % , 10 % , and 30 % . The size of the unit cell was 10 µm, and the capsule sizes were 3.5 and 6 µm for the 10 and 30% contents, respectively, which are in the ranges of similar compositions [59,60,61]. The FEM mesh consisted of approximately 2000 elements.
The unit cells were initially set to T ( x , 0 ) = 15 °C, and the temperature at the left edge, x = 0 , identified with the label “a” in Figure 3, was suddenly set to T ( 0 , y , t ) = 30 °C. The non-linear heat transfer problem was solved for the first 200 s. As observed in Figure 8, at the microscale it took approximately 36 s to complete the phase transition, which in our paraffin wax takes place at T f = 19 °C. We have considered a temperature window of Δ T f = 2 °C and the thermophysical properties in Table 2.
Figure 9 shows the evolution of the temperature of the points labeled as a and e as per Figure 3, i.e., the left and right edges of the unit cell. We used two types of model: (i) a heterogeneous model accounting explicitly for the phase-change phenomenon in the PCM particles (results are shown as continuous lines in Figure 9), and (ii) a homogeneous model using the homogenization approach presented in the previous section (results are shown as dotted lines in Figure 9). We can observe how the temperature rises during the first 9 s from 15 °C to 18 °C, where the phase transition starts taking place in our approach (we use a temperature window of ± 1 °C with respect to the melting point ( T f = 19 °C). Thus, we observed a change in the slope of the curves when the phase transition, solid to liquid, started.
In the case of no PCMs ( V f = 0 % ), it can be seen that there was no change in slope, and the temperature at the right edge reached that of the one imposed at the left edge at 200 s. This means that steady-state conditions were reached. As we increased V f , we observed that the material required more time to achieve these conditions. Thus, the addition of PCM delays the rising of temperature, and we also confirmed that the higher the PCM amount, the higher the thermal inertia exhibited by the composite material. In fact, with the heterogeneous model, we measured 1.8× delay until reaching 95% of the steady-state conditions for V f = 10 % , and this number increased up to 3.4× in the case of V f = 30 % . We also observed some differences between the homogeneous model, in which one phase with homogenized thermophysical properties was modeled, and the heterogeneous one. These differences may have come from the fact that the homogenized expressions correspond to a 3D case (inclusions were treated as spheres), whereas the heterogeneous model considered inclusions as disks. In any case, the same meshes were used for both models. With the the homogenized model being more conservative in terms of thermal management capabilities, we used it in the brick simulation.

3.2. Brick Characterization

Once the new mix (i.e., PCM and cement) had been characterized, we proceeded with the characterization of the building envelopes. For this reason, we tested different designs (summarized in Table 3) using the hot plate setup. We also performed FEM simulations of the bricks to better understand the effects of V f and ϕ a i r on the global insulating properties of the bricks. As pointed out in Section 2.2.1, thermographic analysis was used to measure the temperature evolution within the brick. Figure 10 shows the temperature distribution obtained through thermographic analysis and the corresponding estimates using the FEM model with homogenized properties.
As a first step in our modeling tasks, we needed to calibrate the thermal conductance of the brick-plate interface. In Figure 11 we show the experimental results (in dotted lines) and the equivalent simulation results (in continuous lines) of the large air cell model, which was the one used for the calibration, since we have it for the three different grades of PCM. From the results, we can observe that a conductance of h c = 200 W/m2K provides a good fit. This value of hc is in agreement with that observed for concrete–steel interfaces [62]. We considered an air film coefficient of h a i r = 10 W/m2K.

3.2.1. Effect of Vf

In Figure 12 we show the temperature evolution during ca. 2.5 h, for the large air cell specimens with V f = 10 % (Figure 12a) and V f = 30 % (Figure 12b). In the first place, we can observe that the increase in PCM content provided higher thermal inertia and insulating properties. In fact, during the first 3000 s, the temperature rise in the mix with higher V f increased at a lower pace than that of V f = 10 % . The steady-state conditions were reached after close to 4000 s.
If we look at the 90% height, which is the most interesting line, for it is closest to the indoor one, we can observe again that after 1000 s, the difference in temperature was about 2.5 °C between both mixes. The steady-state conditions were also milder in the case of V f = 30 % .

3.2.2. Effect of Air Cells

We now study the effect of air cells in the insulating properties of the bricks. In Figure 13a it can be seen that at 60% of the height of the samples without PCM, the temperature difference after 6000 s was greatest between the samples with the small air cells and the large air cells. The sample with the large air cells heated up the slowest, as expected, since the volume of the air cells covered 44% of the total volume of the cementitious sample, whereas the small air cells and the medium air cells only covered 30% and 24% of the total volume, respectively. The fact that the sample with the medium-sized air cells heated up slower than the one with the small air cells does not seem to be in line with the literature. The same seems to happen in Figure 10, where the average temperature at 90% of the height of the samples is shown. The sample with small air cells heated up faster than the sample with medium air cells, which is not in proportion with the volume fraction of air within the different samples. Just like in the case of the PCM-containing samples, the graphs for 60% and 90% height are very similar in terms of shape and temperature difference between the samples.
In Figure 13b, it can be observed that larger hole specimens (L) presented a temperature decrease of at least 5 °C with respect to small hole specimens (S) when no PCM was added, and up to 7.5 °C with V f = 30 % .

3.2.3. Thermomechanical Assessment

Finally, we used the temperature field obtained from the heat transfer simulations to run a sequential FEM of the mechanical problem. We considered an elastic modulus E = 20 GPa, as observed by [49], a Poisson’s ratio ν = 0.2 , and thermal expansion coefficient α = 10 5 1/K. We ran the simulation for the first 3000 s, where the most important temperature gradients appeared.
In Figure 14, we show the total deformation of the brick, for large air cells and V f = 10 %. One can see the free expansion deformation mode with the highest values at the bottom tips.
Figure 15 shows the maximum (Figure 15a) and minimum (Figure 15b) principal stresses at the same instant of time. As expected from the deformation field, the highest maximum principal stress appeared at tips of the bottom surface. It must be remarked that we considered sliding boundary conditions for the bottom surface. Apart from the stress concentration due to boundary conditions, we observed well-developed stress bands below the tensile strength (ca. 2 MPa). Regarding the minimum principal stresses, these are well below the compressive strength as measured in [49]. This is in agreement with the experimental observations, where no cracks were identified.

4. Conclusions

The design of new insulating envelopes is a direct route towards energy efficient buildings. The use of advanced numerical and manufacturing techniques may lead to important advances in this sense. In this work we combined in a seamless fashion new materials (e.g., PCMs), manufacturing (e.g., 3D printing), and simulation techniques (e.g., finite element method) in the design of new insulating components with enhanced properties. Namely, we presented the so-called HEXCEM component, which has important benefits, such as improved thermal insulation, reduced energy consumption and GHG emissions, reduced specific weight, and applicability (it could be also used in new buildings as well as in the rehabilitation of older buildings).
For this purpose, a new synthesis approach based on additive manufacturing of tailor-made molds and microencapsulated PCM has been successfully implemented. As we showed in this work, complex geometries can be easily accomplished.
In order to reinforce the design process for both the material mix and the component, a two-scale (material and component) simulation strategy has been successfully developed. For that, the numerical simulations carried out with our methodology showed excellent agreement with the experimental tests. Moreover, it can be used to explore new designs for different applications before actually synthesizing them.
From the analyses, we found that microencapsulated PCM in cement, within the phase-transition temperature, may lead to 1.8× delay to 95% steady-state conditions in the case of V f = 10 % , and 3.4× for V f = 30 % . On the other hand, larger hole specimens (L) present a temperature decrease of at least 5 °C with respect to small hole specimens (S) when no PCM is added, and up to 7.5 °C with V f = 30%.
With respect to the thermomechanical behavior, the simulations of the bricks showed no thermal cracking. This is also in agreement with the experimental observations. This is important for ensuring mechanical integrity.

Author Contributions

Conceptualization, B.Š. and F.M.-C.; methodology, Y.X., B.Š. and F.M.-C.; numerical analysis, J.M.-M. and E.R.; resources, J.C.S.-R. and B.Š.; writing—original draft preparation, J.M.-M. and F.M.-C.; writing—review and editing, J.C.S.-R. and B.Š.; visualization, E.R. All authors have read and agreed to the published version of the manuscript.

Funding

Y. Xu and B. Šavija acknowledge the financial support of the European Research Council (ERC) within the framework of the ERC Starting Grant Project “Auxetic Cementitious Composites by 3D printing (ACC-3D),” grant agreement number 101041342.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

F.M.-C., J.M.-M., E.R. and J.C.S.-R. would like to acknowledge the partial financial support for this project from the Research Fund of the Universidad Loyola Andalucía. Y.X. and B.Š. would like to acknowledge Ties Verhagen for his support with the experimental setup.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Agency, I.E. Energy Technology Perspectives 2016: Towards Sustainable Urban Energy Systems; IEA: Paris, France, 2016. [Google Scholar]
  2. Zhou, D.; Zhao, C.Y.; Tian, Y. Review on thermal energy storage with phase change materials (PCMs) in building applications. Appl. Energy 2012, 92, 593–605. [Google Scholar] [CrossRef] [Green Version]
  3. Mariano-Hernández, D.; Hernández-Callejo, L.; Zorita-Lamadrid, A.; Duque-Pérez, O.; Santos García, F. A review of strategies for building energy management system: Model predictive control, demand side management, optimization, and fault detect & diagnosis. J. Build. Eng. 2021, 33, 101692. [Google Scholar] [CrossRef]
  4. Sarakinioti, M.V.; Turrin, M.; Konstantinou, T.; Tenpierik, M.; Knaack, U. Developing an integrated 3D-printed façade with complex geometries for active temperature control. Mater. Today Commun. 2018, 15, 275–279. [Google Scholar] [CrossRef] [Green Version]
  5. Yang, L.; Yan, H.; Lam, J.C. Thermal comfort and building energy consumption implications–a review. Appl. Energy 2014, 115, 164–173. [Google Scholar] [CrossRef]
  6. Navarro, L.; De Gracia, A.; Colclough, S.; Browne, M.; McCormack, S.J.; Griffiths, P.; Cabeza, L.F. Thermal energy storage in building integrated thermal systems: A review. Part 1. Active storage systems. Renew. Energy 2016, 88, 526–547. [Google Scholar] [CrossRef] [Green Version]
  7. Navarro, L.; De Gracia, A.; Niall, D.; Castell, A.; Browne, M.; McCormack, S.J.; Griffiths, P.; Cabeza, L.F. Thermal energy storage in building integrated thermal systems: A review. Part 2. Integration as passive system. Renew. Energy 2016, 85, 1334–1356. [Google Scholar] [CrossRef] [Green Version]
  8. D’Agostino, D.; Mazzarella, L. What is a Nearly zero energy building? Overview, implementation and comparison of definitions. J. Build. Eng. 2019, 21, 200–212. [Google Scholar] [CrossRef]
  9. Berardi, U.; Soudian, S. Experimental investigation of latent heat thermal energy storage using PCMs with different melting temperatures for building retrofit. Energy Build. 2019, 185, 180–195. [Google Scholar] [CrossRef]
  10. Mindess, S. Developments in the Formulation and Reinforcement of Concrete; Woodhead Publishing: Sawston, UK, 2019. [Google Scholar]
  11. Meyer, C. Concrete and sustainable development. ACI Spec. Publ. 2002, 206, 501–512. [Google Scholar]
  12. Chen, J.; Wang, H.; Xie, P.; Najm, H. Analysis of thermal conductivity of porous concrete using laboratory measurements and microstructure models. Constr. Build. Mater. 2019, 218, 90–98. [Google Scholar] [CrossRef]
  13. Al-Jabri, K.S.; Hago, A.; Al-Nuaimi, A.; Al-Saidy, A. Concrete blocks for thermal insulation in hot climate. Cem. Concr. Res. 2005, 35, 1472–1479. [Google Scholar] [CrossRef]
  14. Šavija, B. Smart crack control in concrete through use of phase change materials (PCMs): A review. Materials 2018, 11, 654. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Beltrán, J.I.; Wang, J.; Montero-Chacón, F.; Cui, Y. Thermodynamic modeling of nitrate materials for hybrid thermal energy storage: Using latent and sensible mechanisms. Sol. Energy 2017, 155, 154–166. [Google Scholar] [CrossRef]
  16. Šavija, B.; Schlangen, E. Use of phase change materials (PCMs) to mitigate early age thermal cracking in concrete: Theoretical considerations. Constr. Build. Mater. 2016, 126, 332–344. [Google Scholar] [CrossRef] [Green Version]
  17. Ling, T.C.; Poon, C.S. Use of phase change materials for thermal energy storage in concrete: An overview. Constr. Build. Mater. 2013, 46, 55–62. [Google Scholar] [CrossRef]
  18. Ricklefs, A.; Thiele, A.M.; Falzone, G.; Sant, G.; Pilon, L. Thermal conductivity of cementitious composites containing microencapsulated phase change materials. Int. J. Heat Mass Transf. 2017, 104, 71–82. [Google Scholar] [CrossRef] [Green Version]
  19. Drissi, S.; Ling, T.C.; Mo, K.H.; Eddhahak, A. A review of microencapsulated and composite phase change materials: Alteration of strength and thermal properties of cement-based materials. Renew. Sustain. Energy Rev. 2019, 110, 467–484. [Google Scholar] [CrossRef]
  20. Alkhazaleh, A.H.; Almanaseer, W.; Ismail, M.; Almashaqbeh, S.; Farid, M. Thermal and mechanical properties of cement based-composite phase change material of butyl stearate/isopropyl palmitate/expanded graphite for low temperature solar thermal applications. J. Energy Storage 2022, 50, 104547. [Google Scholar] [CrossRef]
  21. Royon, L.; Karim, L.; Bontemps, A. Thermal energy storage and release of a new component with PCM for integration in floors for thermal management of buildings. Energy Build. 2013, 63, 29–35. [Google Scholar] [CrossRef]
  22. Mankel, C.; Caggiano, A.; Koenders, E. Thermal energy storage characterization of cementitious composites made with recycled brick aggregates containing PCM. Energy Build. 2019, 202, 109395. [Google Scholar] [CrossRef]
  23. Khan, Z.; Khan, Z.; Ghafoor, A. A review of performance enhancement of PCM based latent heat storage system within the context of materials, thermal stability and compatibility. Energy Convers. Manag. 2016, 115, 132–158. [Google Scholar] [CrossRef]
  24. Lecompte, T.; Le Bideau, P.; Glouannec, P.; Nortershauser, D.; Le Masson, S. Mechanical and thermo-physical behaviour of concretes and mortars containing phase change material. Energy Build. 2015, 94, 52–60. [Google Scholar] [CrossRef]
  25. Kuznik, F.; David, D.; Johannes, K.; Roux, J.J. A review on phase change materials integrated in building walls. Renew. Sustain. Energy Rev. 2011, 15, 379–391. [Google Scholar] [CrossRef] [Green Version]
  26. Ying Kong, S.; Yang, X.; Chandra Paul, S.; Sing Wong, L.; Šavija, B. Thermal response of mortar panels with different forms of macro-encapsulated phase change materials: A finite element study. Energies 2019, 12, 2636. [Google Scholar] [CrossRef] [Green Version]
  27. Mankel, C.; Caggiano, A.; Ukrainczyk, N.; Koenders, E. Thermal energy storage characterization of cement-based systems containing microencapsulated-PCMs. Constr. Build. Mater. 2019, 199, 307–320. [Google Scholar] [CrossRef]
  28. Khudhair, A.M.; Farid, M.M. A review on energy conservation in building applications with thermal storage by latent heat using phase change materials. Energy Convers. Manag. 2004, 45, 263–275. [Google Scholar] [CrossRef]
  29. Cabeza, L.F.; Castell, A.; Barreneche, C.D.; De Gracia, A.; Fernández, A. Materials used as PCM in thermal energy storage in buildings: A review. Renew. Sustain. Energy Rev. 2011, 15, 1675–1695. [Google Scholar] [CrossRef]
  30. Bao, X.; Tian, Y.; Yuan, L.; Cui, H.; Tang, W.; Fung, W.; Qi, H. Development of high performance PCM cement composites for passive solar buildings. Energy Build. 2019, 194, 33–45. [Google Scholar] [CrossRef]
  31. Essid, N.; Eddhahak-Ouni, A.; Neji, J. Experimental and numerical thermal properties investigation of cement-based materials modified with PCM for building construction use. J. Archit. Eng. 2020, 26, 04020018. [Google Scholar] [CrossRef]
  32. Li, M.; Ma, T.; Liu, J.; Li, H.; Xu, Y.; Gu, W.; Shen, L. Numerical and experimental investigation of precast concrete facade integrated with solar photovoltaic panels. Appl. Energy 2019, 253, 113509. [Google Scholar] [CrossRef]
  33. Castell, A.; Martorell, I.; Medrano, M.; Pérez, G.; Cabeza, L.F. Experimental study of using PCM in brick constructive solutions for passive cooling. Energy Build. 2010, 42, 534–540. [Google Scholar] [CrossRef]
  34. Lagou, A.; Kylili, A.; Šadauskienė, J.; Fokaides, P.A. Numerical investigation of phase change materials (PCM) optimal melting properties and position in building elements under diverse conditions. Constr. Build. Mater. 2019, 225, 452–464. [Google Scholar] [CrossRef]
  35. Blanco, J.M.; Frómeta, Y.G.; Madrid, M.; Cuadrado, J. Thermal performance assessment of walls made of three types of sustainable concrete blocks by means of fem and validated through an extensive measurement campaign. Sustainability 2021, 13, 386. [Google Scholar] [CrossRef]
  36. Jayalath, A.; Aye, L.; Ngo, T.; Mendis, P. Multi-scale analysis on thermal properties of cement-based materials containing micro-encapsulated phase change materials. Constr. Build. Mater. 2020, 254, 119221. [Google Scholar] [CrossRef]
  37. Fachinotti, V.D.; Bre, F.; Mankel, C.; Koenders, E.A.; Caggiano, A. Optimization of multilayered walls for building envelopes including PCM-based composites. Materials 2020, 13, 2787. [Google Scholar] [CrossRef]
  38. Cifuentes, H.; Montero-Chacón, F.; Galán, J.; Cabezas, J.; Martinez-De la Concha, A. A finite element-based methodology for the thermo-mechanical analysis of early age behavior in concrete structures. Int. J. Concr. Struct. Mater. 2019, 13, 41. [Google Scholar] [CrossRef]
  39. Aguayo, M.; Das, S.; Maroli, A.; Kabay, N.; Mertens, J.C.; Rajan, S.D.; Sant, G.; Chawla, N.; Neithalath, N. The influence of microencapsulated phase change material (PCM) characteristics on the microstructure and strength of cementitious composites: Experiments and finite element simulations. Cem. Concr. Compos. 2016, 73, 29–41. [Google Scholar] [CrossRef] [Green Version]
  40. Wu, Z.; Xu, Y.; Šavija, B. Mechanical Properties of Lightweight Cementitious Cellular Composites Incorporating Micro-Encapsulated Phase Change Material. Materials 2021, 14, 7586. [Google Scholar] [CrossRef]
  41. Montero-Chacón, F.; Sanz-Herrera, J.A.; Doblaré, M. Computational multiscale solvers for continuum approaches. Materials 2019, 12, 691. [Google Scholar] [CrossRef] [Green Version]
  42. Eddhahak-Ouni, A.; Drissi, S.; Colin, J.; Neji, J.; Care, S. Experimental and multi-scale analysis of the thermal properties of Portland cement concretes embedded with microencapsulated Phase Change Materials (PCMs). Appl. Therm. Eng. 2014, 64, 32–39. [Google Scholar] [CrossRef]
  43. Özdemir, I.; Brekelmans, W.; Geers, M. Computational homogenization for heat conduction in heterogeneous solids. Int. J. Numer. Methods Eng. 2008, 73, 185–204. [Google Scholar] [CrossRef] [Green Version]
  44. Kodjo, K.M.; Yvonnet, J.; Karkri, M.; Sab, K. Multiscale modeling of the thermomechanical behavior in heterogeneous media embedding phase change materials particles. J. Comput. Phys. 2019, 378, 303–323. [Google Scholar] [CrossRef] [Green Version]
  45. Zhi, J.; Raju, K.; Tay, T.E.; Tan, V.B.C. Multiscale analysis of thermal problems in heterogeneous materials with Direct FE2 method. Int. J. Numer. Methods Eng. 2021, 122, 7482–7503. [Google Scholar] [CrossRef]
  46. Vigneshwaran, K.; Sodhi, G.S.; Muthukumar, P.; Subbiah, S. Concrete based high temperature thermal energy storage system: Experimental and numerical studies. Energy Convers. Manag. 2019, 198, 111905. [Google Scholar] [CrossRef]
  47. Santamouris, M.; Kolokotsa, D. Passive cooling dissipation techniques for buildings and other structures: The state of the art. Energy Build. 2013, 57, 74–94. [Google Scholar] [CrossRef]
  48. Rodríguez, C.R.; de Mendonça Filho, F.F.; Figueiredo, S.C.; Schlangen, E.; Šavija, B. Fundamental investigation on the frost resistance of mortar with microencapsulated phase change materials. Cem. Concr. Compos. 2020, 113, 103705. [Google Scholar] [CrossRef]
  49. Šavija, B.; Zhang, H.; Schlangen, E. Influence of microencapsulated phase change material (PCM) addition on (micro) mechanical properties of cement paste. Materials 2017, 10, 863. [Google Scholar] [CrossRef] [Green Version]
  50. Šavija, B.; Luković, M.; Kotteman, G.M.; Figuieredo, S.C.; Schlangen, E. Development of ductile cementitious composites incorporating microencapsulated phase change materials. Int. J. Adv. Eng. Sci. Appl. Math. 2017, 9, 169–180. [Google Scholar] [CrossRef] [Green Version]
  51. Lamberg, P.; Lehtiniemi, R.; Henell, A.M. Numerical and experimental investigation of melting and freezing processes in phase change material storage. Int. J. Therm. Sci. 2004, 43, 277–287. [Google Scholar] [CrossRef]
  52. Thiele, A.M.; Sant, G.; Pilon, L. Diurnal thermal analysis of microencapsulated PCM-concrete composite walls. Energy Convers. Manag. 2015, 93, 215–227. [Google Scholar] [CrossRef]
  53. Bergheau, J.M.; Fortunier, R. Finite Element Simulation of Heat Transfer; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  54. Montero-Chacón, F.; Zaghi, S.; Rossi, R.; García-Pérez, E.; Heras-Pérez, I.; Martínez, X.; Oller, S.; Doblaré, M. Multiscale thermo-mechanical analysis of multi-layered coatings in solar thermal applications. Finite Elem. Anal. Des. 2017, 127, 31–43. [Google Scholar] [CrossRef] [Green Version]
  55. Felske, J. Effective thermal conductivity of composite spheres in a continuous medium with contact resistance. Int. J. Heat Mass Transf. 2004, 47, 3453–3461. [Google Scholar] [CrossRef]
  56. Attaran, M. The rise of 3-D printing: The advantages of additive manufacturing over traditional manufacturing. Bus. Horizons 2017, 60, 677–688. [Google Scholar] [CrossRef]
  57. Kumar, R.; Singh, R.; Farina, I. On the 3D printing of recycled ABS, PLA and HIPS thermoplastics for structural applications. PSU Res. Rev. 2018, 2, 115–137. [Google Scholar] [CrossRef]
  58. Kazemian, A.; Yuan, X.; Cochran, E.; Khoshnevis, B. Cementitious materials for construction-scale 3D printing: Laboratory testing of fresh printing mixture. Constr. Build. Mater. 2017, 145, 639–647. [Google Scholar] [CrossRef]
  59. Shah, K.W. A review on enhancement of phase change materials-A nanomaterials perspective. Energy Build. 2018, 175, 57–68. [Google Scholar] [CrossRef]
  60. Inaba, H.; Kim, M.J.; Horibe, A. Melting heat transfer characteristics of microencapsulated phase change material slurries with plural microcapsules having different diameters. J. Heat Transf. 2004, 126, 558–565. [Google Scholar] [CrossRef]
  61. Wang, X.; Niu, J.; Li, Y.; Zhang, Y.; Wang, X.; Chen, B.; Zeng, R.; Song, Q. Heat transfer of microencapsulated PCM slurry flow in a circular tube. AIChE J. 2008, 54, 1110–1120. [Google Scholar] [CrossRef]
  62. Tao, Z.; Ghannam, M. Heat transfer in concrete-filled carbon and stainless steel tubes exposed to fire. Fire Saf. J. 2013, 61, 1–11. [Google Scholar] [CrossRef]
Figure 1. Micrograph of dispersed PCM microcapsules (reprinted from [49]). PCMs from the same batch were used in the current work.
Figure 1. Micrograph of dispersed PCM microcapsules (reprinted from [49]). PCMs from the same batch were used in the current work.
Materials 15 04925 g001
Figure 2. Mesh for different unit cells: (a) V f = 0%, (b) V f = 10%, and (c) V f = 30%.
Figure 2. Mesh for different unit cells: (a) V f = 0%, (b) V f = 10%, and (c) V f = 30%.
Materials 15 04925 g002
Figure 3. Boundary conditions of the unit cell problem. The centerline indicates different probe locations.
Figure 3. Boundary conditions of the unit cell problem. The centerline indicates different probe locations.
Materials 15 04925 g003
Figure 4. Summary of the specimens production based on CAD and additive manufacturing techniques.
Figure 4. Summary of the specimens production based on CAD and additive manufacturing techniques.
Materials 15 04925 g004
Figure 5. (a) Molds of the bricks. (b) Cast specimens.
Figure 5. (a) Molds of the bricks. (b) Cast specimens.
Materials 15 04925 g005
Figure 6. Hot plate setup. Labels: (A) insulation, (B) hot plate, (C) specimen.
Figure 6. Hot plate setup. Labels: (A) insulation, (B) hot plate, (C) specimen.
Materials 15 04925 g006
Figure 7. Boundary conditions of the component FEM model and probe locations.
Figure 7. Boundary conditions of the component FEM model and probe locations.
Materials 15 04925 g007
Figure 8. Phase transition at the unit cell problem ( V f = 30%. Top label (in parentheses): Average liquid fraction. Bottom label: time.
Figure 8. Phase transition at the unit cell problem ( V f = 30%. Top label (in parentheses): Average liquid fraction. Bottom label: time.
Materials 15 04925 g008
Figure 9. Temperature evolution of different compositions obtained for the heterogeneous model and compared to homogenized properties.
Figure 9. Temperature evolution of different compositions obtained for the heterogeneous model and compared to homogenized properties.
Materials 15 04925 g009
Figure 10. Temperature distribution within the brick during the hot plate test: (a) Thermographic analysis. (b) FEM simulations.
Figure 10. Temperature distribution within the brick during the hot plate test: (a) Thermographic analysis. (b) FEM simulations.
Materials 15 04925 g010
Figure 11. Numerical vs. experimental temperature evolution for V f = 0 % and large air cells at different heights of the brick.
Figure 11. Numerical vs. experimental temperature evolution for V f = 0 % and large air cells at different heights of the brick.
Materials 15 04925 g011
Figure 12. Numerical vs. experimental temperature evolution for large air cells at different heights of the brick: (a) V f = 10 % and (b) V f = 30 % .
Figure 12. Numerical vs. experimental temperature evolution for large air cells at different heights of the brick: (a) V f = 10 % and (b) V f = 30 % .
Materials 15 04925 g012
Figure 13. Experimental temperature evolution of: (a) V f = 0 % and (b) All specimens at 90%.
Figure 13. Experimental temperature evolution of: (a) V f = 0 % and (b) All specimens at 90%.
Materials 15 04925 g013
Figure 14. Total deformation in the brick after 3000 s.
Figure 14. Total deformation in the brick after 3000 s.
Materials 15 04925 g014
Figure 15. Principal stresses in the brick after 3000 s: (a) Maximum principal stress. (b) Minimum principal stress.
Figure 15. Principal stresses in the brick after 3000 s: (a) Maximum principal stress. (b) Minimum principal stress.
Materials 15 04925 g015
Table 1. Composition of the mixes (g/L) for different PCM volume fractions.
Table 1. Composition of the mixes (g/L) for different PCM volume fractions.
Constituent V f = 0% V f = 10% V f = 30%
CEM I 42.51354.81219.8948.6
water610.2548.8426.9
PCM090270
Table 2. Thermophysical properties of the phases [16,49].
Table 2. Thermophysical properties of the phases [16,49].
Material PropertyCementParaffin Wax
Density, ρ (kg/m3)1965900
Conductivity, k (W/mK)1.000.42
Specific heat, cp (J/kgK)15301900
Enthalpy of fusion, hf (J/kg)n/a143,500
Melting point, T f (°C)n/a19
Table 3. Building envelope configurations tested.
Table 3. Building envelope configurations tested.
ReferenceSize (mm) ϕ a i r ϕ Curing Time (Days)
Large (L) 100 × 52 × 20 0.440.0, 0.1, 0.321, 21, 28
Medium (M) 100 × 52 × 20 0.240.021
Small (S) 100 × 52 × 20 0.300.021
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marin-Montin, J.; Roque, E.; Xu, Y.; Šavija, B.; Serrano-Ruiz, J.C.; Montero-Chacón, F. Thermomechanical Performance Analysis of Novel Cement-Based Building Envelopes with Enhanced Passive Insulation Properties. Materials 2022, 15, 4925. https://doi.org/10.3390/ma15144925

AMA Style

Marin-Montin J, Roque E, Xu Y, Šavija B, Serrano-Ruiz JC, Montero-Chacón F. Thermomechanical Performance Analysis of Novel Cement-Based Building Envelopes with Enhanced Passive Insulation Properties. Materials. 2022; 15(14):4925. https://doi.org/10.3390/ma15144925

Chicago/Turabian Style

Marin-Montin, Jorge, Eduardo Roque, Yading Xu, Branko Šavija, Juan Carlos Serrano-Ruiz, and Francisco Montero-Chacón. 2022. "Thermomechanical Performance Analysis of Novel Cement-Based Building Envelopes with Enhanced Passive Insulation Properties" Materials 15, no. 14: 4925. https://doi.org/10.3390/ma15144925

APA Style

Marin-Montin, J., Roque, E., Xu, Y., Šavija, B., Serrano-Ruiz, J. C., & Montero-Chacón, F. (2022). Thermomechanical Performance Analysis of Novel Cement-Based Building Envelopes with Enhanced Passive Insulation Properties. Materials, 15(14), 4925. https://doi.org/10.3390/ma15144925

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