1. Introduction
The interest in phase change materials (PCM) grows continuously due to its wide application in sectors as the storage of solar thermal energy, industrial heat recovery, construction, agricultural greenhouses, aerospace, health, refrigeration, among others [
1,
2,
3]. An adequate characterization is a mandatory requirement for numerical modeling of processes involving PCM, equipment design, and the development of this technology and its applications, since thermal characterization provides information about the amount of energy that a material can store. There are multiple thermal properties that must be defined for an adequate characterization of the PCM. Most of the existing standards were not developed specifically for thermal storage applications. However, they are currently used as the basis for determining the thermal properties of these kind of materials [
4].
Among the main standards applied to the characterization of PCMs using DSC, the following technical standards are distinguished: ASTM D87, which is the standard test method for melting point of petroleum waxes; ASTM D4419 for measurement of transition temperatures of petroleum waxes; ASTM E793, which is a standard method to determine enthalpies of fusion and crystallization with DSC. Similarly, there are more standards that are used for the determination of thermal conductivity, which depend on the measurement devices, as well as other standards for determination of density, viscosity, among others. These standards are applicable since there is currently no exclusive standard for PCMs, and they are applied more as a general use of the measuring device.
In the last thirteen years, several authors have worked in the development of standards for PCM, being notable the contributions especially in the determination of thermophysical properties by DSC. Differential scanning calorimetry is the most robust technique for determining the phase change enthalpy, there are several criteria that must be considered when characterizing a PCM with this technique [
5], for example, (i) the type of calorimetry, which can be made in dynamic or step mode. In this sense, Barreneche et al. [
6] did not observe significant differences between both methods for paraffins. (ii) the heating/cooling rate. As is proposed by Mehling et al. [
7], a thermal equilibrium must be achieved, which only can occurs at low heating/cooling rates, since high ramps of velocity generate artificial effects such as hysteresis and (iii) the mass of the sample must be constant. Günther et al. [
8] found that the signal peak shifts towards higher temperatures as the mass or heating rate increases.
Other procedures have been developed in recent years to standardize DSC measurements such as equipment calibration, determination of the required cooling/heating rate, experimental procedure, and representation of the obtained data [
9]. Some of the recommendations for traceability in measurements suggest very extensive tests, which involve testing at different cooling/heating rates until, for example, in the case of the enthalpy-temperature curve in the DSC the maximum value in temperature differences cannot exceed the limit of 0.2 K. The quality criteria and measurement methods for PCM were established by the ZAE Bayern research center and the Fraunhofer ISE, and since 2008 the RAL label has been awarded to products that have been tested under these exhaustive monitoring and testing time criteria [
10].
Specific heat is also other important property in the characterization of PCM, especially the values before and after the phase change (solid and liquid phase). In this case, the ASTM E1269 standard is used for the determination of the specific heat by differential scanning calorimetry, however Ferrer et al. [
11] compared the dynamic, isostep and area methods, finding that the results obtained by the area method had an error lower than 3%, while the dynamic and the isostep errors up to 6% and 16%, respectively. Thermal conductivity can be measured by different methods, but the three main measurement methods used are: (i) hot disk, (ii) hot wire and (iii) laser flash, each of these techniques require specialized equipment and are mainly based on ISO 22007, ASTM D7896, ASTM E1461.
Other devices for measuring thermal properties are especially required when evaluating large or non-homogeneous samples, as material encapsulated in polymeric matrices or expanded graphite, likewise composite materials for the construction industry, textile, among others. Cabeza et al. [
12] presented a series of devices that were built for the determination of properties, among which stand out for the measurement of enthalpy-temperature, diffusivity, and thermal conductivity curves.
The T-history method was proposed by Yinping & Yi [
13], this requires three measurements during the cooling process of a PCM sample with unknown properties and a reference substance with well-known thermal properties, which are cooled in ambient air, that is monitored too. The cooling curve of the sample and the reference are recorded, from an initial temperature value, higher than the melting temperature of the sample until both reach an ambient temperature. Some improvements have been proposed to the original method over the years, in which the values of specific heat (cp) and enthalpy (h) can be calculated as temperature dependent values using the same experimental procedure [
14]. Reference [
15] applied the T-history method to a wide variety of PCMs such as paraffins and lauric acid having no or low degree of supercooling. While some authors propose assemblies with the test tubes arranged in a vertical position according to the original method [
15], others propose the horizontal configuration to avoid the temperature difference in the tube [
16]. With T-history also is possible estimate of the “effective thermal capacity”
(T), which contains more information than any other phase change property [
17].
Comparative results between measurements of dynamic DSC, step mode DSC and T-history, show that the precision of the latter is similar, but allows the investigation of large samples, which is crucial for heterogeneous materials like composite PCM or salt hydrates [
8]. There are many investigations that have been carried out taking as reference the T-history for the determination of the thermophysical properties [
18], some salt hydrates like disodium hydrogen phosphate dodecahydrate, sodium acetate trihydrate and STL-47 have been evaluated by this method [
19]. The thermophysical properties of other PCMs such as paraffin waxes, salt hydrates and mixtures of fatty acids have also been evaluated by T-history and compared by DSC results [
20]. Thermophysical properties for alcohol, WFI, Greek paraffin, RT20, RT27, RT58, SP25A8 are reported in the literature [
21]. Other non-homogeneous materials have been evaluated by t-history method, this is the case of the commercial sodium acetate trihydrate-based PCM (Climsel C58), which was blending with additives for their characterization [
22]. Five paraffin waxes and wood resin (all available in Indian market) were investigated to determine their suitability as a thermal energy reservoir in solar drying applications [
23]. Also, this method has been adapted to be able to determine the latent heat in some granular phase change composites like RUBITHERM GR27 and GR41 EPCM [
24]. PCM for air conditioning applications, such as semiclathrate hydrates of tetra-n-butylammonium fluoride (TBAF) have been evaluated by this method, confirming that the T-history method is a practical method to determine the phase change enthalpy of hydrates [
25]. Recently, the use of T-history has been reported in nano-enhanced (NEPCM) applications, in which the thermophysical parameters of multiwall carbon nanotubes (MWCNT) in capric acid were obtained [
26].
Some setups of the T-history have been made using insulated sample holders. In some cases the effects of the thermal mass of the insulation can be neglected using sample holders with a high ratio between the thermal mass of the PCM to the insulated sample holder [
27], The evaluation method for insulated T-History measurements was reported taking into account the heat capacity of the insulation material [
28]. The T-history has critically assessed based on their mathematical formulation and experimental configuration [
29]. The effect of radial thermal gradients inside T-history samples on the enthalpy temperature curve measurement have been analysed [
30], and has been and has been proposed a differential formulation of the T-History method in order to achieve thermodynamic consistency with phase transitions theory, including the experimental effect of the speed of the thermal process in calorimetric calculations [
31].
In the present work, the T-history method was used to determine the enthalpy-temperature curves, on which the specific heat in solid and liquid phase are also obtained, this method was used due to its simplicity and low cost, results were comparing with measurements in DSC. However, both methods have their limitations, which generate convenience in the use depending of the application case, for instance, in the case of samples with subcooling phenomenon, the use of T-history is recommended, since with the DSC is not possible to obtain a representative measure of the enthalpy variation with temperature [
32]. If the sample of material to be tested is not pure and/or non-homogeneous, it is more useful the use T-history, since the sample sizes used in DSC is only a few milligrams, while for the T-history reach several grams, it is a greater challenge to obtain a representative sample in non-homogeneous materials if the sample size is to small [
12].
One of the major limitations of the original T-history lies in temperature control, which is performed at room temperature [
13]. However, researchers at ZAE-Bayern [
32] proposed configurations where tests between −20 to 65
C could be performed, Kravvaritis et al. [
17] added a cooling and heating system in which tests could be performed between −30 to 120
C, Stankovic [
33] used an environmental chamber with advanced temperature control, which operated in a temperature range between −10 up to 100
C with an accuracy of ±0.2
C. All these modifications contribute to the reliability of the measurements, but they make the method much more expensive than the simplicity that characterized the original method.
The conventional method DSC have a small chamber in which samples are evaluated in much wider temperature ranges, typically between −90 C up to 550 C, depending on the heating and cooling system they have incorporated. Therefore, this is a versatile and robust equipment with the capability to evaluate almost all the organic and inorganic substances currently used in thermal storage applications.
Another proposed method for the determination of several thermophysical properties simultaneously was proposed in 2018 by Yang & Liu [
34], in which Stefan’s problem is solved in 1D with constant heat flow boundary condition, using a quasi-static fusion model, the method is known as “T-melting CHF”.
This paper presents a general evaluation of two low-cost experimental methods (T-history and T-melting) proposed in the literature for the characterization of PCM, with a comparison of the results obtained by a specialized method (DSC) and the data provided by the manufacturer. Considering the simplicity and cost of the methods, these results can serve as the basis for determining properties in laboratories were specialized equipment for this purpose are not available. Likewise, the aim of this work is to establish a consensus between the difficulties and advantages presented by each method on the determination of the main thermophysical properties of PCM. In the first part, the experimental methods developed are presented, then the conditions and results of each of the tests developed are defined, and finally a comparison is made between the calculated properties and the effectiveness of each of the methods used.
4. Discussion
To compare the results obtained between the T-history method and the DSC, the cooling rate of the PCM samples is calculated from the T-history. It is represented in
Figure 16, which shows that for both RT45 and RT55, the cooling rate vary in some sections, but in general terms they are lower than 0.2
C/min on average, which represents cooling rates slightly lower than those evaluated by DSC.
Figure 17 shows one of the greater advantages of T-history method over DSC and its use in thermal energy storage applications: It allows to evaluate a representative sample of non-pure PCM in conditions close to those usually presented in actual practice. Other advantage is that when insulated PCM samples are evaluated, lower cooling speeds are obtained, and in the general case of DSC equipment that is cooled by nitrogen, noise signals are induced at low cooling speeds, which complicates the estimation of the specific heat.
Results are compared in an Enthalpy-Temperature diagram in
Figure 17 When comparing the amount of heat stored in the range given by the manufacturer, an average
value of 159.74 kJ/kg was obtained, which shows the consistency between the results obtained by the T-history with respect to those obtained by DSC.
Figure 18 presents the superimposed curves between the results obtained by the DSC, those calculated with the T-history method, and the equations proposed by Kravvaritis et al. [
17] for the effective heat capacity. Correspondence is evident between the results; it was expected that at much lower cooling rates in the T-history method than in the DSC, a higher peak would be obtained.
When comparing the conventional, unconventional experimental results and the values reported in the manufacturer’s data sheet, discrepancies are evident (
Table 8), among which it is highlighted that: (i) values of the phase change temperatures and enthalpies of fusion/solidification can be obtained by the T-history and T-melting methods. (ii) For the calculation of the conductivity, results obtained by the T-melting method show good correspondence with the values presented by the manufacturer. Colla et al. [
37] experimentally determined the conductivities of the RT45 and RT55 paraffins using an instrument based on the hot disk technique, (ThermTest TPS 2500s) to measure thermal conductivity. The instrument operates in transient plane source mode using a Kapton insulated disk sensor in order to be thermally neutral. Values obtained experimentally by these authors for the RT45 paraffin was 0.2415 W/m·K while for the RT55 it was 0.3336 W/m·K. Therefore, when results obtained and provided by the manufacturer of 0.2 W/m·K and those calculated by the T-melting CHF are compared, greater differences are found. Nevertheless, the T-melting CHF method allows obtaining a general estimate of the property. (iii) The stored energy can be determined by the T-history in the evaluated range, while with the T-melting CHF method an estimated phase change value is obtained that is consistent with the results obtained by the DSC. (iv) the specific heat in liquid phase calculated by the CHF method differs from that reported by the manufacturer, which can be justified not only by a poor insulation of the experimental module, as mentioned by [
34], but rather by the simplification of the model, which does not consider the convective effects in the liquid phase. so it is not recommended for the estimation of this property. However, the calculated solid phase specific heat is a good overall estimate of the property. Although the manufacturer provides the same value for the specific heat in both phases for its possible proximity and ease of handling, it is clear that the value of this property is lower in the liquid phase than in the solid phase, this is in agreement with the results obtained from calorimetry, which is consistent with the results found with the T-history method.
Finally, it is evident that the data provided by the manufacturer serves as a general guide for its use in thermal energy storage applications. However, it is advisable to determine the thermal properties under the same experimental conditions in which they will be used. Although the need of standards for the determination of each of the properties of phase change materials is evident, it is understandable that this remains a technical challenge since they can be affected by the inherent characteristics of the process in which it will be used. However, when it is required to guarantee the equilibrium values, it is required to proceed in accordance with what was reported by Gschwander et al. [
9], in which a very slow velocities for heating and cooling are recommended.