Introduction

The world energy consumption is expected to rise nearly 50% by 2050 [1]. This quick and high growing energy demand can generate a depletion of resources, supply difficulties and destructive environmental impact. 50% of the EU´s energy consumption is used for heating and cooling, from which 80% is consumed in buildings [2, 3]. Consequently, the 2015 Paris Agreement established a series of objectives to renovate the housing stock, prioritizing energy efficiency by improving the thermal insulation or applying new building designs, promoting the development and use of sustainable energy resources [3, 4]. Buildings present a large surface exposed to the sun light, being possible to use this renewable energy as an alternative source to satisfy their thermal energy demands [5]. However, this energy has a low density and is intermittent, which hinders its collection and storing.

The sensible heat of materials for the sun energy storage has been widely used in buildings. However, in recent years, the usage of the latent heat has gained importance being now one the most studied alternatives because of the high thermal capacity per mass unit and the absence of temperature change during the phase transition make this alternative very attractive for thermal energy storage (TES) [6]. Phase change materials (PCMs) are substances that are able to absorb and store large amounts of thermal energy in the form of latent heat. The mechanism of PCMs for energy storage relies on the increased energy need of some materials to undergo phase transition. A PCM absorbs or releases energy by taking advantage of its latent heat when temperature undergoes or overpasses its melting point [7].

PCMs can be applied in two kind of TES systems, the passive energy storage systems and active ones. In the passive energy storage, the PCM is incorporated into building materials. The active energy storage strategies cover heating, ventilation and air conditioning (HVAC) systems, where a new kind of fluid appears, the named Phase Change Slurry (PCS). The PCS is a fluid composed of a continuous liquid phase, such as water, and a dispersed solid one (PCM) which is stable and behaves as a fluid, permitting its transport, pumping an implementation in active energy storage systems [8, 9].

The passive systems have been implemented and studied since 1980 [10] by incorporating the PCM in the walls, shutters, ceilings and floors, for which gypsums, concrete and polyurethane foams are frequently used [4, 11, 12].

However, some authors have reported undesirable interactions of the PCMs and the building material or their leakage when they are incorporated directly [13, 14]. These problems can be solved by encapsulating the PCMs before mixing with the building material. One of the first and better studied polymer shells for PCM microencapsulation is polystyrene [15,16,17,18]. However, these microcapsules presented low mechanical resistance [19], resulting in a bad choice for long-lasting applications such as the enhancement of building thermal comfort. To overcome this drawback, polystyrene has been copolymerized with methyl-methacrylate [20,21,22] and with divinyl-benzene [23], resulting in similar thermal energy storage capacities. Polymethylmethacrylate [24,25,26] and melamine–formaldehyde [27,28,29,30] polymer shells have been also widely investigated for PCMs microencapsulation. However, their application in buildings seems to reduce the mechanical properties in a far extension than those with poly(styrene-co-divinylbenzene) shells [31, 32]. Besides, the methyl-methacrylate derived polymers tend to present related swelling problems. On the other hand, low density polyethylene and ethyl vinyl acetate copolymers [19] or polymers derived from chitosan (CS) and carboxymethyl cellulose [33] have been also employed but they are produced by methods incompatible with the production of a stable slurry. Thus, a PCS from poly(styrene-co-divinylbenzene) copolymer was selected as the most suitable thermal liquid for manufacturing gypsum composites.

On the other hand, although there is extensive research on the use of microcapsules (1–1000 µm) (MPCM) in TES [11, 16, 34], the current trend is the submicron- (0.1–1 µm) or nanoencapsulation (\(\le\) 100 nm) of the PCM [35, 36], since they presented better long-term stability, higher heat transfer area and better structural stability [37, 38]. Borreguero, A. M. et al.[12] manufactured smart gypsum composites by adding different kinds of MPCMs in order to develop building materials with high TES capacity. The incorporated MPCMs had an average particle size between 3 and 117 µm. They achieved up to 15 mass% of MPCMs with respect to the hemihydrate (Hem), obtaining a TES capacity of 10.57 kW m−1 K−1 and a decrease in thermal conductivity, with a minimum value of 0.20 W m−1 K−1. However, they observed a non-homogeneous distribution of MPCMs throughout the gypsum matrix, as well as a large amount of particle agglomerations, and this effect worsened with increasing capsule size and the amount of MPCMs in the composites. There are further investigation of MPCMs addition in gypsums [39,40,41] with these same problems, which may result in areas with worse thermal properties than others, as well as a worsening of mechanical properties.

On the other hand, Mohammadi, B. et al. [20] incorporated nanocapsules of butyl palmitate as core and polystyrene-co-methyl methacrylate as shell in a gypsum wall. They demonstrated that compared with the unmodified gypsum, the gypsum wall containing NPCM took longer time to get to the thermal equilibrium, 600 and 8400 s, respectively. Maleki, B. et al. [21] synthesized submicrocapsules (SPCMs) from polymethylmethacrylate modified with CuO nanoparticles containing n-dodecanol as phase change material. They manufactured gypsum-SPCM composites at reduced-scale considering a small room with PCM and another one without it. The results indicated that the PCM system exhibited narrow indoor air temperature fluctuations and maintained the indoor thermal comfort at most of the time of the year. They achieved an increase in heat capacity and a decrease in thermal conductivity of ~ 50 and ~ 47%, respectively, for 50% by mass of SPCMs with respect to the hemihydrate. However, these previous investigations did not study the NPCM distribution in the composite material and their influence on the physical and mechanical properties of the new building material. Smart gypsum composite boards up to 30 mass% n-octadecane NPCMs were fabricated by Khadiran, T. et al. [22]. They observed a reduction in the density, both bulk and real, of the building materials with the increase of the NPCMs load, as well as an increase in porosity. They detected, as well, a non-homogenous dispersion of the NPCMs in the gypsum matrix by FESEM. Nevertheless, the smart gypsum had good thermal properties, obtaining a reduction in indoor temperature of 4.6 K respect to the gypsum without NPCMs. As in the previous mentioned work, the mechanical properties of the smart gypsum were not studied.

In a previous work, it was found that the particle size of the capsules and the agglomeration played an important role in the final composites properties [12]. To the best of our knowledge, the use of a concentrated slurry containing PCM nanocapsules with particles size lower than 100 nm to hydrate the hemihydrate, allowing the development of gypsum composites, presenting homogeneous filler distribution, large thermal energy accumulation and low subcooling, without losing their insulating capacities has not been reported yet. Some articles described the production of gypsum boards containing PCMs but the used fillers do not satisfy the European Commission and ISO/TS 80,004–1 classification as nanomaterials, where most of the half of nanoparticles (in number) must be below 100 nm and also those gypsum blocks did not keep the thermal conductivity of the original gypsum.

In this research, a NPCS has been synthesized following our patent [42], containing NPCMs from poly(styrene–divinylbenzene) with paraffin wax Rubitherm®RT27 as PCM. Hence, thermally enhanced gypsums have been produced, employing such NPCS as additive. It was possible to reach a filler content up to 41% of NPCM with respect to the gypsum hemihydrate (Hem). The physical, morphological, thermal and mechanical properties of the developed composites have been studied.

Materials and methods

Materials

Black gypsum, supplied by BigMat S.A. and milli-Q water with a conductivity value lower than 5 µS cm−1 were used for fabricating the gypsum blocks composites. The nanoencapsulated PCM slurry (NPCS) was synthesized by the chemical engineering department of the University of Castilla-La Mancha according to the patent number ES2804063 [42]. The characteristics of this NPCS, the NPCM concentration (in mass) in the NPCS ([NPCM]) measured by gravimetry, the zeta potential (ζ) by Doppler laser microelectrophoresis technique, the particle size in volume and number (dv0.5 and dn0.5, respectively) by Dynamic Light Scattering (DLS) and the latent heat of the NPCM (ΔHNPCM) by Differential Scanning Calorimetry (DSC) are listed in Table 1.

Table 1 Characteristics of the NPCS added to the gypsums

Gypsum composites fabrication

Gypsum blocks were synthesized by weighting firstly the known mass of nanoencapsulated PCM slurry (NPCS) for obtaining the desired mass ratio of PCM-based nanocapsules regarding to the hemihydrate (NPCM/Hem). Then, when an extra amount of water was required (in addition to the water from NPCS), it was added to the mixer container to keep a constant water/Hem mass ratio at 0.6. After that, the required mass of Hem was included together with the mixture NPCS-water using a mixer that satisfies the European regulation UNE-EN 196–1 [43]. The mixing process is carried out at slow speed for 1.5 min. Finally, the mixture was poured into a mold of 40 \(\times\) 40 \(\times\) 160 mm for the mechanical test, to a mold of 30 \(\times\) 60 \(\times\) 100 mm for the thermal analysis and to a mold of 20 \(\times\) 110 \(\times\) 110 mm for the studying the thermal conductivity in lambdometer. The composites will be named as C-NPCM/Hemx, where “x” is the PCM nanocapsule content in mass regarding hemihydrate (NPCM/Hem mass ratio). The recipes used for the manufactured C-NPCM/Hem are shown in Table 2.

Table 2 Amount of Hem, NPCS and water in mass for the manufactured composites

Characterization

Fourier transform infrared (FT-IR) spectroscopy

Fourier Transform Infrared (FT-IR) spectroscopy analysis was carried out on a Spectrum Two spectrometer (Perkin Elmer, Inc) equipped with a Universal Attenuated Total Reflectance (UATR) accessory. Three samples in different positions of the same specimen were carried out, and an average value has been taken as a result. All infrared spectra were collected using 16 scans and 16 cm−1 resolution in the wavelength range of 4000 to 500 cm−1.

Scanning electron microscopy (SEM) analysis

Manufactured gypsum composites were depicted by means of Scanning Electron Microscopy (SEM) by using a Quanta 250 (FEI Company) with a tungsten filament operating at a working potential of 15 to 20 kV. The Backscattered Electron Detector (BSED) was applied for imaging.

X-Ray diffraction (XRD) analysis

The crystallography of the obtained materials was analyzed by using an X-ray diffractometer (Philips) model X, Pert MPD. The equipment has a radiation CuKa1, automatic divergence slit, graphite monochromator and xenon gas sealed detector.

Thermogravimetric analyses (TGA) of gypsum composites

The thermal stability of the synthesized composites was analyzed by thermogravimetric analysis (TGA) using a SDT Q600 Simultaneous DSC-TGA (TA Instruments). Three different samples in different positions from the same gypsum composite specimen were carried out, in order to study the homogeneity of the composites. Those samples were taken with a spatula from the mechanical specimen gypsums blocks at 0, 200 and 400 mm. The total sample mass was around 2 g. Temperature was tuned from room temperature to 973 K, under a nitrogen atmosphere and using a heating rate of 10 K min−1.

Density and porosity

The bulk density of the gypsum block (\(\rho_{{{\text{bulk}}}}\)) was determined by weighing and sizing the specimens. The real density (\(\rho_{{{\text{real}}}}\)) was determined by helium pycnometry (Micromeritics Accupy II 1340). The porosity of the gypsum composites (\(\varepsilon\)) were estimated with Eq. 1 from the \(\rho_{{{\text{bulk}}}}\), \(\rho_{{{\text{real}}}}\) and assuming that the gypsum pores are filled of air which remains trapped once the gypsum has solidified. The air density (\(\rho_{{{\text{air}}}}\)) was assumed to be 0.001 g cm−3, which corresponds to the \(\rho_{{{\text{air}}}}\) at normal conditions.

$${\upvarepsilon } = { }\frac{{\rho_{{{\text{real}}}} - \rho_{{{\text{bulk}}}} }}{{\rho_{{{\text{real}}}} - \rho_{{{\text{air}}}} }}$$
(1)

Procedure and experimental equipment thermal behavior tests of the composite gypsum blocks

The thermal behavior of gypsum composites additive with NPCS has been analyzed using a homemade, where the different temperature profiles of the specimens are registered with time when they are subjected to a hot or cold focus (only in one of its faces) and always starting from a steady state to a new one as temperature function. The equipment consisted of an aluminum box having a thickness of 1 mm through which a liquid coming from a thermostatic bath was fluid by means of a peristaltic pump. The temperature of the aluminum cell is controlling with the temperature of this liquid. The box has 30 \(\times\) 60 \(\times\) 100 mm, being the same dimensions of the tested specimen. Each gypsum block was place on the upper surface of the cell and further insulated with foams boards of 4 cm in thickness. This set is put inside of a thermostatic chamber which allows to keep the internal temperature at 298 K. A more exhaustive description of the equipment and of its performance for the thermal characterization of materials were described and proved in previous works [11, 44,45,46].

Tests were carried out applying a thermostatic bath set-point steps change from 291 to 318 ± 0.1 K and from 318 to 291 ± 0.1 K. Six thermocouples of K-type were used to measure temperatures: two were put in the external test tube surface (\(T_{{{\text{up}}}} )\), other two were placed at the middle of the test tube (\(T_{{{\text{middle}}}} )\) and the last ones in the cell (\(T_{{{\text{down}}}} )\). The liquid flow was a flow of 9 L h−1 to ensure the absence of thermal profile in liquid direction. Six heat flow sensors gSKIN®-XI and gSKIN®-XP were putted on the six faces of the test tubes in order to monitor online the inlet and outlet heat fluxes. The signals of the thermocouples and heat flow sensors were registered continuously using the NOKEVAL program and recorded by means of a computer. Using these signals, it is possible to quantify the, the \({\text{TES}}\) (kWh m−3) capacity per cubic meter (Eq. 3) of the gypsum composite:

$${\text{TES capacity}} = { }\frac{{q_{{{\text{acc}}}} }}{{m_{{{\text{sample}}}} }} \times \frac{{\rho_{{{\text{bulk}}}} }}{{3.6 \times 10^{6} }}$$
(2)

where \(q_{{{\text{acc}}}}\) is the amount of accumulated heat in the sample during the experiment (J), \(m_{{{\text{sample}}}}\) is the specimen mass (kg).

The thermal characterizations were performed three times.

Differential scanning calorimetry (DSC)

The DSC model 214 Polyma from Netzsch (equipped with a cooling system and nitrogen at 40 mL min−1 as the purge gas) was used to measure the \(c_\text{p}\) of the gypsums block composites employing Proteus Analysis software to analyze data. In particular, the software extension “Ratio Methat” allowed to obtain the \(c_\text{p}\) value of the material. DSC was carried out in three samples taken at different positions (0, 20 and 40 mm) of the same gypsums composite specimen to demonstrate the homogenous distribution of the nanocapsules contained in it. Sapphire was used as calibration or reference material before measuring the sample in order to obtain the real heat capacity of composites at each temperature (not only the differential heat flux regarding air), using the abovementioned proteus Analysis software extension to analyze data. So, in order to carry out the measure of the sample \(c_\text{p}\), the following methodology was followed:

  • 1. Blank measurement (empty crucible).

  • 2. Sapphire measurement as standard (reference material).

  • 3. Sample measurement of the gypsum composite to evaluate the material \(c_\text{p}\).

For the three steps, the same temperature ramp was used (253 to 323 K, up and down, at 3 K min−1).The latent heat capacities were determined by integrating the \(c_\text{p}\) curve regarding to the baseline.

For thermal cycling (heating/cooling) of the thermoregulating gypsum with the maximum load in NPCM (41%, in mass), such temperature cycle was performed between 263 and 313 K, at a heating/cooling rate of 3 K min-1 to ensure the applicability of the gypsum composite in realistic operating cycles.

Lambdometer

A lambdometer model HFM 446 Lambda medium (NETZSCH) was used, in order to quantify the thermal conductivity. In that way, thermal conductivity of different specimens with dimensions of 110 \(\times\) 110 \(\times\) 20 mm were measured at 281, 289, 297, 301, 305, 313 and 318 K with a temperature gradient of 20 K between upper and bottom plates.

Mechanical tests

The mechanical properties of the gypsum composites were studied by using a Servohydraulic Test System (MTS Landmark). The flexion and compressive strengths were calculated according to standard UNE-EN 13,279–2 [47]. The flexion test was done using specimens samples of 40 \(\times\) 40 \(\times\) 160 mm applying a load rate of 50 N s−1 until fracture. The compression test was done using the broken flexion test specimens. A load rate of 1 N mm−2 s−1 was applied over 40 \(\times\) 40 mm surface load application until fracture. From these results, the Young´s modulus (E) is determined as the slope value of the initial part of the compression curves. The mechanical tests were performed three times at room temperature (293 K), three times at PCM melting point (308 K) and three times when the whole PCM is melted (318 K).

Results and discussion

Fourier transform infrared spectroscopy

To confirm the presence of the NPCM in the gypsum blocks, FT-IR spectroscopy analyses were performed in three different positions of the same specimens for each NPCM/Hem mass ratio. As an example, the media of the FT-IR spectra of the C-NPCM/Hem0, C-NPCM/Hem0.2 and C-NPCM/Hem0.41 and the FT-IR spectra of the pure NPCM are shown in Fig. 1. In the FT-IR spectra of the NPCM appear the three characteristic absorption bands for paraffin wax: rocking vibration of –CH2 at 721 cm−1; deformation vibration of –CH2 and –CH3 at 1466 cm−1; and three intensive peaks between 3000—2840 cm−1, which correspond to alkyl stretching vibrations of –CH2 and –CH3 [48, 49]. This spectrum presents also absorption bands at 3020, 750 and 695 cm−1, attributed to aromatic C–H stretching vibration, corresponding to the polymeric shell of poly(styrene–divinylbenzene). In addition, peaks located between 1600 and 1500 cm−1 confirmed the presence of C=C phenyl stretching [50,51,52]. In the spectrum of gypsum without NPCM, the band within at 3200–3800 cm−1 denotes asymmetric stretching band of crystalline water. The band observed at 1150 cm−1 represented the symmetric stretching vibration of O–H group in the water. Besides, the band recorded at 1619 cm−1 was attributed to the bending vibration of S–O group in CaSO4 [53, 54]. The spectra of the composite gypsums with NPCM present the characteristic bands of the NPCM and the pure gypsum, demonstrating, undoubtedly, that the NPCMs are part of the gypsum matrix. Further, in Fig. 1, it can be observed that the characteristic bands from NPCM increase with its content in capsules.

Fig. 1
figure 1

FT-IR analyses of the gypsum composites C-NPCM/Hem0, C-NPCM/Hem0.2 and C-NPCM/Hem0.41 and the pure NPCM

Scanning electron microscopy (SEM) analysis

The presence of fillers can modify the morphology of calcium sulphate dihydrate crystal [47]. As an example, the structures of the composites C-NPCM/Hem0, C-NPCM/Hem0.2, C-NPCM/Hem0.3 and C-NPCM/Hem0.41 are shown in Fig. 2.

Fig. 2
figure 2

SEM images of gypsums with different NPCM/Hem mass ratio

According to Fig. 2, the crystal size of the gypsum increases with the presence of the NPCM; as an example, comparing the width of the biggest crystals found for C-NPCM/Hem0 and C-NPCM/Hem0.41, they increased from 2.0 to 7.3 µm, respectively. Besides, the shape of the crystal is affected by them. The unmodified gypsum shows a needle-shape crystal structure (C-NPCM/Hem0). However, with the increase of NPCM tends to be hexagonal in shape. The use of slurries containing NPCMs and surfactants that are often crystal morphology modifiers to hydrate the hemihydrate, can affect the morphology in two ways; (1) the setting behavior with respect to time, and (2) the shape formation. Hence, according to the morphology, either PCM nanocapsules joined with the stabilizers promote Ostwald ripening during later stages of crystallization, or the agglomeration of nanoparticles takes place reducing the nucleation points during the setting stage. The nanocapsules agglomeration can be confirmed by SEM, since the areas indicated by white circles correspond to pores occupied by NPCMs. The location of the lower density NPCM (dark areas) can be distinguished from the higher density gypsum crystals (bright areas) in the SEM images taken in backscattered mode.

X-ray diffraction (XRD) analysis

Figure 3.a shows the XRD analysis of the gypsums manufactured with different amounts of NPCM.

Fig. 3
figure 3

a X-Ray Diffraction patterns of the manufactured gypsums. b Full width at half maximum (FWHM) and corresponding mean size of the crystallite (t) obtained for the main peak located at 2θ = 11.6º using the Debye–Scherrer formula for the gypsums prepared with different NPCM/Hem mass ratio

In Fig. 3.a, it is observed that all the gypsums are mainly constituted by crystalline monoclinic calcium sulfate dihydrate (CaSO4·2H2O), with no important fingerprints from unreacted hemihydrate in XRD measurements, only revealing certain minor marks of unreacted orthorhombic hemihydrate (CaSO4·½H2O), in the peaks located at  = 25.4º, 29.5º and 32.1º. These vestiges of unreacted CaSO4·½H2O were a bit more pronounced for the reference gypsum without PCM (C-NPCM/Hem0) than for the gypsums including NPCM, exposing the dihydrate formation is slightly enhanced by the addition of NPCM. In contrast, all the gypsums showed very clear and narrow peaks related with monoclinic dihydrate (CaSO4·2H2O), principally located at  = 11.6º, 20.7º, 23.4º, 29.1º, 31.1º, 33.3º, 35.9º, 40.6º, 43.5º, 47.9º, 50.3º and 51.3º.

In Fig. 3.b, it is represented the decrease of FWHM as the NPCM/Hem ratio is augmented, revealing an increase of the crystallinity and domain size as the dose of NPCM is increased. From the reference gypsum (C-NPCM/Hem0) up to the gypsum C-NPCM/Hem0.3, the FWHM loss is moderate (from 0.113 to 0.108); however, for the gypsum containing the maximum amount of NPCM (C-NPCM/Hem0.41), such decline is more significant (FWHM = 0.091). In order to determine the mean size of the ordered crystalline domains in the manufactured gypsums, the Debye–Scherrer formula has been employed:

$$t = \frac{K \times \lambda }{{\beta \times {\text{cos}}\theta }}$$
(3)

where t is the mean size of the crystallite, K is the Scherrer constant and taken equal to 0.9, λ is the wavelength of the Cu-Kα radiation (1.54056 Å), β is the Full Width at Half Maximum (FWHM) of the X-ray peak in radians and θ is the Bragg angle [55]. Thus, the mean size of the crystallite varies from 3427.0 up to 6491.8 nm, confirming the trend of crystallinity enhancement as the NPCM/Hem ratio is increased. Regarding t, the increase of crystallite size is modest between NPCM/Hem ratios 0.00–0.30, whereas a sudden and prominent augmentation is produced for the maximum dose of NPCM. These XRD findings are in good agreement with the evidence observed in SEM images, since the mean crystallite size increases as the NPCM amount is augmented, exactly the same behavior as the observed for the growth of the gypsums grains in SEM images when the NPCM/Hem ratio was raised.

Density and porosity

Density (Fig. 4.a) and porosity (Fig. 4.b) of the manufactured gypsum composites are shown in Fig. 4.

Fig. 4
figure 4

a Density and b porosity of the composite gypsum blocks. Points were the mean of three specimens with whiskers representing the standard deviation interval

As expected, the higher the NPCM content, the lower the composite bulk and real density (Fig. 4.a). This tendency can be explained attending to the lower density of the NPCM (0.96 g cm−3) respect to that of the gypsum (2.38 g cm−3). Although the bulk densities of the composites decrease with the NPCM content, all of them satisfy the minimum density value of 0.60 g cm−3 established by the European regulation UNE-EN 13,279–1 for gypsum binders and gypsum plasters [56]. Figure 4.b shows a decrease in the composite porosity with the NPCM content. In the SEM images (Fig. 2) can be confirmed that the presence of NPCMs modifies the shape structure of the gypsum crystals, but also their small size allows them to occupy the composite pores.

Thermogravimetric analyses (TGAs) of gypsum composites

TGA of the gypsum composites with different NPCM/Hem mass ratios and the used NPCS are shown in Fig. 5.

Fig. 5
figure 5

a Average of the TGA curves of the gypsum composites with different NPCM/Hem mass ratios and b) TGA curves of the different position of the same specimen with the higher NPCM/Hem mass ratio

In Fig. 5.a, the TGA curve of the C-NPCM/Hem0 presented three regions of mass loss. The first one, observed between 353 and 433 K, is attributed to absorbed water (hygroscopic water) and the chemically bound water of the hydrated salt. The second one, that occurs between 783 and 823 K, corresponds to the water bound to hydraulic compounds and the last one, the smallest one, refers to the CO2 formed during the decomposition of carbonates [57, 58]. The NPCS TGA curve presented four steps of mass loss. The first one, between room temperature and 393 K, is due to the water evaporation; the second one, between 413 and 483 K, corresponds to the RT27 evaporation; the third step, between 483 and 533 K, can be attributed to the surfactant degradation, and the last one, between 653 and 723 K is appropriate to the polymer degradation. Finally, a 1.5% remains none degraded at 973 K due to the inorganic content of the surfactants. In the TGA curves of the gypsum composites can be observed four clear mass losses due to the absorbed water, the PCM, the polymer and the water bound to hydraulic compounds. The PCM and the polymer content present in each composite according to the TGA results is shown in Table 3. It is worthy to point out that the biggest deviation, between the three TGA in the different points of the specimens, was only a 0.7% in the case of PCM, for C-NPCM/Hem0.3. So, it is possible to say that the NPCMs have a good distribution in the composite when it is added in the form of NPCS, as it was confirmed by the TGA analyses of different points of the block from C-NPCM/Hem0.41 (Fig. 5.b).

Table 3 Average and standard deviation of mass loss of the PCM and polymer in the gypsum composites with different NPCM/Hem mass ratios analyzed three times

Thermal behavior of the composite gypsum blocks

Figure 6 shows the heat capacity of the synthesized composites using the NPCS as additive. The heat capacity has been calculated by DSC from two kind of experiments, one with a heating ramp and the other with cooling ramp.

Fig. 6
figure 6

Heat capacity versus temperature of the different composite gypsums measured by DSC

As expected, the heat capacity (\(c_{p}\)) of the gypsum blocks increases with the increase in the NPCM/Hem mass ratio (Fig. 6). It was calculated the average \(c_{p}\) (Table 4) in the PCM fusion interval (between 292 and 305 K) and the increase in the \(c_{p}\) with respect to the unmodified gypsum (C-NPCM-Hem0) was determined, with values of around 80, 140, 165, 208 and 269%, for C-NPCM/Hem0.1, C-NPCM/Hem0.15, C-NPCM/Hem0.2, C-NPCM/Hem0.3 and C-NPCM/Hem0.41, respectively. The latent heat of each composite was calculated by integrating the \(c_{p}\) versus temperature curve in the paraffin melting temperature range, obtaining values of 0, 9.7, 16.5, 18.0, 25.7 and 30.2 J g−1 for C-NPCM/Hem0, C-NPCM/Hem0.1, C-NPCM/Hem0.15, C-NPCM/Hem0.2, C-NPCM/Hem0.3 and C-NPCM/Hem0.41, respectively. The latent heat relationship with the NPCM/Hem mass ratio was found to be linear, as can be observed in Fig. 7.

Table 4 Summary composite gypsums thermal properties. Represented the mean values of three replicated measurements together with the standard deviation
Fig. 7
figure 7

Latent heat versus temperature of different composites measured by DSC. Points were the mean of three specimens with whiskers representing the standard deviation interval

Figure 8 shows the temperature profiles of the specimens containing NPCM when they were subjected to the heating process.

Fig. 8
figure 8

Temperature profiles of the external surface of the composite blocks containing NPCS additive with the NPCS, a Heating and b Cooling

The results shown in Fig. 8 demonstrate that the higher the NPCM content in construction material, the lower the slope of the temperature profiles when the temperature reaches the region in which the PCM starts to melt or solidify (298 K). On contrary to what happens for the case of gypsum without PCM, when PCM is present, an inflection point appears, since the thermal energy increases proportionally with respect to the amount of PCM present. This inflection point indicates that the PCM is absorbing the thermal energy leading its melting instead of increasing the gypsum block temperature (Fig. 8.a). The same happens when a temperature drop ramp is applied to the specimen (Fig. 8.b), but in this case the decrease in the temperature change is due to the fact that the PCM is releasing the thermal energy that was previously absorbed. In that way, the specimen C-NPCM/Hem0.41 had the biggest buffering behavior allowing to delay in 2.36 h the heating steady state and in 3.89 h the cooling steady state. This was also confirmed by the increase of the absorbed and released with the addition of the NPCM, determined from the experimental equipment measurements (Fig. 9). Furthermore, it was observed a reduction of the temperature in the steady state for each composite respect to the pure gypsum, which could indicate a decrease of the thermal conductivity.

Fig. 9
figure 9

Instant power and energy a absorbed and b released for the gypsum composites synthesized with the NPCS

The difference between the heating or cooling sources respect to the temperature of upper surfaces is motivated by the thermal conductivity, the specimen thickness and the chamber temperature (298 K). The thermal conductivity and the thickness are practically the same, then the gap is mainly influenced by the difference between the steady state with the temperature of the environmental chamber. In that way, the narrower temperatures in the cooling process respect to the heating process is due to the gap between 291 and 318 K respect to 298 K, respectively.

So, it can be also said that the insulating capacity of the composites has been improved, it could be explained attending to the low thermal conductivity of the NPCM which is within 0.1—0.2 W m−1 K−1, which is lower than that of pure gypsum (0.2–0.8 W m−1 K−1) [46, 59,60,61].

From the thermal analyses results, TES capacities per cubic meter in the range of the specimen temperature change were determined by using Eq. [2] (Table 4). These values increase with the increase of the NPCM content due to the latent heat of fusion of the PCMs and the polymer heat capacity. The TES capacities per cubic meter of the composite C-NPCM-Hem0.1, C-NPCM-Hem0.15, C-NPCM-Hem0.2, C-NPCM-Hem0.3 and C-NPCM-Hem0.41 increased by around 57, 89, 100, 133 and 140%, respectively, compared with C-NPCM-Hem0. The TES capacity values are represented in Table 4. On the other hand, when temperature is below the melting point, the manufactured composite gypsums release all the thermal energy that they had absorbed. Hence, with these materials it is also possible to increase the accumulation of energy and also reduce the thermal loss because they have a thermal conductivity lower than the standard gypsum.

Figure 10 presents the heat flux of the external surface of the specimens, which would correspond with the heat flow throughout the indoor surface of a wall, it means, the incoming or outgoing heat in a room. In that way, the negative heat flow values correspond to the heat transfer from the composite to the inside of a room.

Fig. 10
figure 10

Heat flux of the external surface of the composite blocks C-NPCM-Hem0 and C-NPCM-Hem0.41 for heating and cooling step

It is plain to see that, in the heating step at the transient state the C-NPCM-Hem0 exhibited a high incoming heat flux. On the contrary, the uptake of energy by the PCM in the composite (gray area) softened the incoming heat until temperature reached 298 K. Once that temperature was achieved, the greater portion of the composite had overpassed its melting temperature (Fig. 6), having no capacity to uptake more energy and increasing in that way the incoming heat flux into the room. Thus, the incorporation of NPCM greatly reduced the incoming heat from the outside during the transient states.

During the cooling step, when temperature was below 298 K, the heat flow reaches a plateau region that finishes at 296 K. This plateau area represents the solidification of the PCM, process in which the composite would transfer its previously stored heat to the room.

Finally, the thermal conductivity of all the manufactured specimens at different temperatures were measured by lambdometer (Fig. 11).

Fig. 11
figure 11

Thermal conductivity the composite gypsum blocks at different temperatures determined by using a lambdometer. Points were the mean of three specimens with whiskers representing the standard deviation interval

Figure 11 shows the thermal conductivity (\(k\)) decrease with the mass ratio of NPCM/Hem increase. The percentages of thermal conductivity reduction were around 17, 17, 30, 35 and 50% at 305 K for C-NPCM-Hem0.1, C-NPCM-Hem0.15, C-NPCM-Hem0.2, C-NPCM-Hem0.3 and C-NPCM-Hem0.41, respectively, compared with the pure gypsum (C-NPCM-Hem0). The measured \(k\) values at 305 K are gathered in Table 4. This reaffirms what was seen with the experimental setup (Fig. 8), where the maximum temperature reached in the steady state decreases when the mass ratio NPCM/Hem increased. Therefore, it can be concluded that the use of NPCS as additive to produce composite gypsums improves their insulating properties.

As a summary, Table 4 shows the thermal properties of the developer composites.

Another important point related to the energy savings is the reduction of the CO2 emissions to the atmosphere. The use of 1 m3 of panels in the correct orientation of gypsum composites containing a 15 and 41 mass% of NPCM respect to the hemihydrate would allow to reduce 2.2 and 3.4 kg of CO2 emissions per absorption/releasing operating cycle within 292–306 K, respectively. This is obtained considering that 0.42 kg of CO2 are released to the atmosphere per one electric kWh of produced energy [62]. In Fig. 12, a DSC thermal cycling for the thermoregulating gypsums of maximum load is shown.

Fig. 12
figure 12

Thermal heating/cooling cycling (− 263 to 313 K) for the C-NPCM/Hem0.41 by DSC

In Fig. 12, it can be perceived that the gypsum composite with maximum NPCS content revealed an excellent thermal cyclability, even after 100 heating/cooling cycles, so leading to a gypsum composite material that can resist storage-release cycles without losing its thermal energy storage capacity, neither presenting direct PCM leakage nor capsule breaking (which can also promote PCM losses) during this process. This large cyclability guarantee the long-term thermal stability of this kind of composites. Additionally, the endothermic and exothermic peaks are very close in temperature, indicating that gypsum reduces the nanocapsule subcooling.

Mechanical tests

Once confirmed the improvement of the thermal properties of the gypsum by using the NPCS as additive, the mechanical properties of the synthesized composites were studied according to the European Standard UNE-EN 13,279–2 [47].

Figure 13 shows the variation of the maximum flexion and compression strengths of the gypsum composites, respectively, as function of the NPCM/Hem mass ratio and temperature, being they tested at room temperature (293 K), PCM melting point (300 K) and once the PCM is melted (318 K).

Fig. 13
figure 13

Mechanical testing of the gypsum composites as function of mass ratio of NPCM/Hem, a flexural strengths and b compression strengths. Points were the mean of three specimens with whiskers representing the standard deviation interval

The behavior of the maximum flexural and compressive strengths of the composites was quite similar, being possible to be explained in the same way. This is, the reduction in the mechanical strengths of the composites occurs either by temperature or NPCM content. The larger the NPCM/Hem mass ratio the lower the maximum flexion and compression strength of the composite. It is due to the reduction in density and the lower strength of the NPCM respect to the pure gypsum. The maximum flexural and compression strengths are higher for the composites tested at 293 K, when the PCM is solid, respect to the values obtained at 300 and 318 K, when the PCM is liquid. It could be explained attending to the PCM changes from solid to liquid, decreasing its strength and, therefore, its contribution to the strength of the composite. Besides, the mechanical strength decreases with the NPCM content. Taking into account the glass transition temperature of PS-DVB, which is higher than 95 ºC, it does not contribute to the losing strength. Then, it is possible that the calcine of gypsum at 45 ºC has a minimum contribution to this effect but it is more dependent on the melting temperature of the PCM because it was not evidenced for the PCM-free gypsum. The maximum flexion and compression strengths reduction was observed for the composite C-NPCM-Hem0.41, reaching up to 85 and 87%, respectively, when the tested temperature is shifted from 293 to 318 K. However, there is no significant variation between the measurements at 300 and 318 K, indicating that the composite thermal conductivity is high enough for ensuring that practically the total PCM melts at its melting point. Nevertheless, it is important to point out that the application of these composites is limited by mechanical strength, being only possible to consider these composites as gypsum binders or gypsum plasters those obtained from a NPCM/Hem mass ratio up to 0.15 which satisfy the European mechanical standard EN 13,279–1 [56].

Considering the physical, thermal and mechanical properties of developed composites, some of these new materials can be used as construction gypsums for special applications, such as, gypsum plasterboard, drywall, plasterboard for ceilings and fiber-reinforced plasterboard [56].

Conclusions

Smart gypsum composites were manufactured using a nanoencapsulated phase change material slurry containing nanocapsules from poly(styrene-co-divinylbenzene) as shell and paraffin wax as core materials. A uniform distribution of the nanocapsules in the composite was observed, being possible to use the slurry to hydrate the hemihydrate. A linear relationship between the latent heat of the gypsum composites and the nanocapsules content was verified. The heat capacity, the thermal energy storage capacity and thermal conductivity of the new building material were improved in 269, 140 and 50%, respectively, when a 0.41 mass ratio of nanocapsules/hemihydrate was used. Composites from a nanocapsules/hemihydrate mass ratio up to 0.15 satisfied European mechanical standard EN 13,279–1 for gypsum binders and gypsum plasters and all of them, presented a bulk density higher than 0.60 g cm−3. The presence of nanocapsules and surfactant changed the morphology of gypsum crystals during the setting stage. Thanks to the incorporation of such high amounts of nanocapsules, it is possible to reduce 3.4 kg of CO2 emissions per operating cycle by using 1 m3 of panels of this lightweight gypsum.