Alkane-based eutectic phase change materials doped with carbon nanomaterials †

Cold thermal energy storage is an issue of increasing importance on a global scale particularly in the format of passive thermal protection. This study presents three eutectic Phase Change Materials (ePCMs) composed of n -alkanes, which provide passive temperature control (their operation is automatically induced by exceeding the limit temperature without the need for a control system) around 4 1 C (277.2 (cid:2) 2 K) and are chemically neutral. The solid–liquid equilibrium (SLE) in the following binary systems was investigated: { n -tetradecane + n -heptadecane}, { n -tetradecane + n -nonadecane}, { n -tetradecane + n -heneicosane}, allowing the determination of two ePCMs with enthalpies close to 220 J g (cid:3) 1 and one significantly lower (155.5 J g (cid:3) 1 ). Moreover, two solid–liquid–liquid equilibrium (SLLE) phase diagrams were determined for systems: { n -tetradecane + 1,6-hexanediol} and { n -tetradecane + 1,12-dodecanediol}. In addition, the work provides a systematic analysis of the problem of designing ePCMs with specific properties and the aspects that need to be considered. The possibility of predicting the parameters of eutectic mixtures using the UNIFAC (Do) equation and the equation of ideal solubility was verified. A method for predicting the enthalpy of melting of eutectic was also proposed and confronted with the results of DSC analysis. Thermodynamic studies have been supplemented with measurements and correlation of experimental data of ePCMs density and dynamic viscosity as a function of temperature. The final issue is the improvement of the thermal conductivity of paraffins by the addition of nanomaterials such as Single Wall Carbon Nanotubes (SWCNTs), Expandable Graphite (GIC) or Expanded Graphite (EG). The possibility of forming a long-lasting composite material composed of ePCMs and 1 wt% of SWCNTs with a thermal conductivity significantly higher compared to pure ePCMs has been proven via stability testing under operating conditions.


Introduction
After previous energy crises in the 1970s and 1980s caused by turmoil in the oil markets, the time of the energy crisis related to gas supply constraints has arrived. 1 Under these circumstances, we note that countries are paying particular attention to diversifying their fuel suppliers and deepening their energy independence. Considering the climate changes that have been taking place in recent times, it is important to look for solutions that provide low temperature conditions at a low cost. The energy consumption of compressor refrigerators (which are the main refrigeration technology) is directly dependent on the temperature of the heat receiver, so care must be taken to ensure good heat exchange at the condenser. However, given the diurnal variation in temperature, energy savings can be achieved by using thermal energy stores that produce cold during the night and use it during the day. 2 Phase-change materials, PCMs, are either effective thermal energy stores (both heat and cold) or passive temperature controllers (as they stabilize the temperature without human intervention or control). When discharging PCMs, there is no threat of significant overcooling below the material's melting point when cooling products, nor is there any risk of burns when heating, as the temperature of the medium does not significantly exceed the crystallisation temperature. For this reason, PCMs, in comparison to other thermal energy storage (TES) systems, can be used without concerns for the protection of materials that are sensitive to minute changes in temperature; however, to ensure this protection, it is necessary to select a PCM with the desired conversion temperature. For this purpose, research is being conducted on eutectic phase change materials (ePCMs), which allow the use of compounds from a group with desirable properties while having an unfavourable melting temperature.
As paraffins in the role of PCMs are characterised by an almost 100% efficiency with simultaneously practically unlimited service time, they represent environmentally friendly and efficient accumulators of thermal energy. It is also important to note their high availability and low cost due to the large scale of production in petroleum refining processes. Companies currently supplying PCMs in bulk quantities to the Chinese market offer the reagents investigated in this paper at prices ranging from $2 to $15 per kilogram, which, considering the long lifetime, seems beneficial. The main advantages of long-chain alkanes include their chemical inactivity and ease of crystallisation resulting in negligible sub-cooling while having a very high latent heat value. In the case of long-chain n-alkanes, the energetic effect of crystallisation is related to a lowering of the system entropy as a result of the formation of a highly organised crystalline structure resulting from weak but multiplied van der Waals interactions and dihydrogen bonds (C-HÁ Á ÁH-C). 3 In order not to lose the high enthalpy of fusion at the ePCM formation stage, it is necessary to select components with positive interactions resulting in a positive mixing enthalpy value. Blending two alkanes together results in a negligibly small negative deviation from ideality due to entropic effects. In order to obtain positive deviations, an alkane can be combined with, for example, an alcohol or carboxylic acid, and such systems have been quite intensively studied. [4][5][6] In contrast, a poorly explored niche so far are diols and in particular, due to their lower propensity to undercool, a,o-diols. [7][8][9] The presence of two hydroxyl groups at the terminal positions of a,o-diols gives them the possibility of energetically favourable packing of the crystal lattice using strong hydrogen interactions resulting in a much higher energy effect of the phase transformation than for an alkane of similar size. Considering the favourable effect of positive deviations on latent heat, in the present work an ePCMs composed of an alkane and an a,o-diols was attempted in addition to the two-alkane systems.
The main problem observed with organic PCMs is their low thermal conductivity resulting in the inertia of the heat stores and temperature stabilizers made from them. The default solution is to increase the heat transfer surface area and therefore to use various types of metal ribbed heat exchangers 10,11 or to place PCM in porous materials with high thermal conductivity. [12][13][14][15] Among the most extensively studied additives with high thermal conductivity are carbon nanotubes (CNTs) whose small addition in the form of a suspension (nanofluid) repeatedly improves thermal properties with negligible degradation of the stored heat density. [16][17][18] Alternatively, the investigated nanomaterial is expanded graphite (EG), whose production by expansion (thermal decomposition) of exfoliated graphite (expandable graphite) is much simpler, resulting in its lower price. [19][20][21] In addition to thermal conductivity, the key properties of functional materials are the values of their physicochemical parameters such as density (which reflects energy storage density) and viscosity (crucial for hydraulic transport processes and reducing the inertia of phase transformations). The coefficients of thermal expansion, determined from the mathematical description, allow the design of suitable vessels resistant to unsealing.
Considering the practical nature of the undertaken research, it seems important to propose theoretical methods for selecting PCMs with the expected properties. For this reason, the appropriateness of using selected mathematical models for estimating the composition and melting enthalpy of eutectic with the expected melting point was experimentally verified.

Chemicals and materials
The reagents used in this study were applied without additional purification and the purity grades were assumed according to the specifications declared by the supplier (Table 1). Alkanes and diols exhibited a narrow melting point range (below 0.3 K) and a water content of no greater than 1000 ppm for diols and less than 20 ppm for alkanes. Expandable graphite (Graphite intercalation compound, GIC) was expanded by heating it at 520 K under nitrogen flow for 30 min resulting in expanded graphite, EG. 22

Solid-liquid phase equilibrium measurements
The course of the liquidus curves was determined using a highly accurate dynamic (synthetic) method, in which the sample is subjected to slow heating (slower than 2 K h À1 ) and the temperature at which the last fragment of the second phase (solid or liquid) disappears is recorded. The investigated samples were prepared by weighing the components on an analytical scale with an accuracy of 1 Â 10 À4 g (u(x) = 5 Â 10 À4 ) and then sealed in a vessel containing a stirring bar. The temperature of the system, measured using a calibrated P 750 (DOSTMANN electronic GmbH) thermometer, had an error of no more than u(T) = 0.05 K. A detailed description of the method can be found in our previous works. 6,23,24  Values for the heat effects of phase transitions (DH) and the temperatures of polymorphic transformations were determined using differential scanning calorimetry, DSC. A DSC 1 STARe (Mettler Toledo) calorimeter operating in a heat-flux mode with liquid nitrogen cooling system had been adjusted prior to measurements using a series of high-purity compounds according to the procedure described in our previous work. 6 The standard error of the onset temperature u(T) = 0.4 K and the relative error of the heat effects u r (DH) = 0.03 were determined for measurements carried out with a heating rate of 2 K min À1 or 5 K min À1 , as indicated by the thermograms. For the study of mixtures, samples of approximately 2 mL were first prepared by weighing the components with a composition accuracy of u(x) = 5 Â 10 À4 , followed by melting and thorough mixing. In liquid form, about 10 mg of the mixture was weighed out with an accuracy of 2 Â 10 À5 g and sealed in a tight aluminium vessel under ambient atmosphere.

Density and viscosity measurements
Density measurements of eutectic mixtures (EM) in the liquid phase were carried out using an Anton Paar GmbH DMA 4500 M oscillating densimeter (Graz, Austria) with two builtin Pt 100 thermometers providing a temperature uncertainty of u(T) = 0.01 K. The experiments were carried out at atmospheric pressure p = 100 kPa, with u(p) = 5 kPa. The precision of the instrument, taking into account correction for viscosity, is 1 Â 10 À5 g cm À3 , while considering the calibration performed using standards provided by the supplier and the purity of the reagents used in the tests, the relative error of the density is u r (r) = 0.001. 25 Viscosity was measured using a DVNext-LV cone/ plate rheometer from AMETEK Brookfield (Middleboro, USA) with a relative uncertainty of dynamic viscosity of u r (Z) = 0.03. Viscosity error was determined using standards: APN26E and APN75 from Paragon Scientific Ltd and n-tetradecane (Alfa Aesar, purity 499%). All measurements were carried out at an atmospheric pressure of p = 100 kPa, (u(p) = 5 kPa) and the temperature uncertainty is u(T) = 0.1 K. In both cases, samples prepared according to the procedure described for DSC were examined.

Composite stability examination
Composite systems consisting of solid carbon material and ePCM were prepared by weighing a presumed portion of the powder (SWCNT, GIC or EG) on an analytical scale and then filling with ePCM previously prepared according to the same procedure as described for DSC measurements. Samples containing a (C14 + C17) EM with SWCNTs (0.1%, 0.3%, 0.5% and    Standard uncertainties are: u(x) = 5 Â 10 À4 , u(T) = 0.05 K, u(p) = 5 kPa.
1% by weight), GICs (10% and 75% by weight) and EG (0.5%, 5%, 10% by weight) were prepared. Prior to testing, the samples were vacuumed to remove air from the porous materials and then, after sealing, homogenised in an ultrasonic cleaner at 298.2 K for one hour. The propensity of the melted systems to sediment under isothermal conditions was tested first. Samples were placed in a thermostat at 298.2 K for 50 days after which the photographs were taken. Then, after re-homogenisation, the vessels were subjected to a series of 1000 temperature cycles to assess the stability of the suspensions during PCM operation. Samples were cooled to 269.2 K, held at this temperature for 10 min and then heated to 283.2 K and also held for 10 min.
During the measurement the temperature of the water bath in which the samples were immersed was recorded. The full description of experiment with the apparatus scheme was published in our previous publication. 8 The appearance of the composites after 1000 cycles.

Thermal conductivity
The liquid-phase thermal conductivity (l) was measured by transient hot-wire method according to the ASTM D7896-19 standard using a THW-L2 thermal conductivity meter (Thermtest Inc., Canada). Samples of approximately 30 mL were prepared by weighing an ePCM and carbonaceous materials on an analytical scale (according to the same procedure as above). Prior to measurement, the samples were placed in a vacuum chamber at 338.2 K to remove air and then homogenised in an ultrasonic bath at 338.2 K for one hour. Composite systems containing SWCNTs at 1 wt% and EGs at 10 wt%, as well as ePCM without additives, were investigated to compare the change in thermal conductivity of the ePCM. A detailed description of the measurement methodology is included in our previous publication. 8

Solid-liquid phase equilibrium
The results of the solid-liquid phase equilibrium (SLE) measurements using the dynamic method are summarized in numerical form in Table 3 and presented graphically as dots in Fig. 1-3. SLE measurements for the system {C14 (1) + C17 (2)} gave a liquidus curve representing a simple eutectic with a single monotonicity change point corresponding to the eutectic point.   In the case of the {C14 (1) + C19 (2)} and {C14 (1) + C21 (2)} systems, in addition to the minimum, an inflection in the solubility line of the odd n-alkanes was observed. In the {C14 (1) + C21 (2)} system, the temperature of inflection (305.6 K) corresponds to the C21 polymorphic transition, suggesting that a pure alkane, rather than its solid solution, is present in equilibrium with the saturated solution. The {C14 (1) + C19 (2)} system, containing alkanes with a smaller chain length difference, despite the different crystal structure showed a slight decrease in the C19 solid-solid phase transition temperature, which is 296.1 K for pure C19 and is observed at 294.4 K on the liquidus. The dependence of the polymorphic transition temperature on the composition of the system shows that the solid phase present in equilibrium with the saturated solution is not pure C19, but its solid solution. The minor difference in the length of the n-alkanes leads, in the case of the {C14 (1) + C17 (2)} system, to a complete lack of inflection on the C17 liquidus despite the existing polymorphic transformation of pure C17 at 284.5 K and suggests the presence of a solid solution in this system over a significant range of concentration. In order to confirm the hypothesis, it was necessary to investigate the thermal effects accompanying the heating of samples with different ratios of C14 to odd n-alkanes.
The DSC thermograms plotted in blue on the phase diagrams confirm the composition dependency of the polymorphic transformation temperature and the melting temperature of the eutectic, while the small values of the thermal effects of these transformations do not allow a complete phase diagram to be explicitly determined using DSC. The continuous black lines represent relationships resulting from the DSC, while the dashed lines show the probable course of the boundaries of occurrence of the individual phases. The proposed phase diagram is consistent with studies carried out for the similar system {C16 + C17}, 26,27 in which all the phases: b T (triclinic of C14 and C16), a H (hexagonal of C17 occurring at higher temperatures) and b O (orthorhombic of C17 at temperatures below the polymorphic transition temperature) 28 form solid solutions in specific concentration ranges. The results of the DSC analysis of the remaining systems give no reason to suspect the presence of solid solutions in significant molar fraction ranges, and confirm the formation of simple eutectic systems by pairs of even alkanes with odd ones at a difference in aliphatic chain length of more than 3 carbon atoms.
The observed results are a consequence of very significant positive deviations from ideal behaviour leading to a miscibility gap in the liquid phase. Achieving the eutectic point requires the use of either a geminal, vicinal or very long a,o-diol as well as a monohydroxy alcohol, which in contrast have lower enthalpies of melting and a higher propensity to overcool.
All liquidus curves were correlated using the NRTL equation (Fig. 5) and compared with the prediction using the ideal solubility equation and the modified (Dortmund) UNIFAC equation (UNIFAC (Do)). The used values of the surface and volumetric parameters R and Q are included in Table S1 in ESI. † In each case, small positive deviations from ideality are visible and the curves predicted by both equations are equivalent.
The consistency of the UNIFAC (Do) results and ideal solubility is due to UNIFAC model's omission of interactions between CH 3 and CH 2 groups belonging to the same CH 2 main group. In a situation when the residual part of the calculated activity coefficient is equal to zero, deviations from ideality are influenced only by the combinatorial term. However, the geometric similarity of the CH 2 and CH 3 groups in most of the composition range (through g = 1) reduces the equation form to ideal solubility. The negligible effect of the difference in the surface parameters (Q) is noticeable only for g 2 values near the eutectic composition.
In all systems except {C14 (1) + C17 (2)}, where the error due to the formation of a solid solution is noticeable, it is unnecessary to determine the eutectic point experimentally in order to use as an ePCM (Table 5).

Latent heat examination
The results of the thermal analysis of n-alkanes are presented in Fig. S1-S4 in ESI. † The polymorphic transformation of C21, in comparison with the other investigated odd alkanes, occurs in a wider temperature range (the peak is broadened and uneven), which results in a less visible change in the monotonicity of the liquidus and blurred thermograms plotted on the phase diagram (Fig. 3). The temperature and transformation enthalpies of the pure components obtained in this work are compared with literature values in Table 2, while the properties of the obtained ePCMs are summarised in Table 6.
The eutectic compositions determined from the SLE experimental data were used to prepare the samples for DSC analysis, the thermograms of which, together with a full description, are presented in ESI, † in Fig. S5-S7. Among the ePCMs studied, the (C14 + C17) EM has the lowest latent heat value (Fig. 6), which is a direct result of the solid solutions presence in this system (Fig. 1).
According to the phase diagram presented in Fig. 1, in contrast to the other two ePCMs, the solid (C14 + C17) EM is not composed of crystals of the individual components, but of their solid solutions. As crystal lattice of the solid solution has a higher entropy than that of the pure alkane, its formation during crystallisation results in a minor entropy difference with respect to the liquid solution (D fus S). The difference in ordering in the solid phases implies the melting enthalpy value of the eutectic mixture according to the relation for the isobaric process: D fus H = TÁD fus S and results in a minor amount of heat that such an ePCM can store.
Knowing the composition and the melting temperature of the eutectic mixture (from experiment or predictive calculations), it   is possible to predict the thermal effect associated with its melting. The following equation, taking into account the temperature dependence of the enthalpy of melting of the components, is most commonly found in the literature. 29,30 where: T E -eutectic temperature, C L P i and C S P i -respectively, the specific heat capacity of the liquid and solid phases of component i. In order for the above equation to be applicable to compounds exhibiting a polymorphic transition, it should be supplemented with the expression D tr H i /T tr . In recent publications, the authors propose to perform an entropy balance starting with crystals of the two components with a eutectic ratio and ending with the molten mixture at the same, eutectic temperature. 30 This approach takes into account the entropy of mixing of the ideal solution, but requires knowledge of the C p of the liquid eutectic mixture (or the application of mixing rules for the regular solution) and is not an intuitive way to include imperfections in the system. According to the above, and bearing in mind that the total enthalpy change does not depend on the path of the transformation, we propose an enthalpy balance based on the thermodynamic cycle described in ESI, † and summarized as follows.
A prerequisite for the application of the above and previous calculations is that the investigated compounds form simple eutectic systems, i.e. systems in which no solid solution is present, and the liquid solution is in equilibrium with the pure solid component. In the case of a phase transition in a solid in the temperature range under study, the enthalpy of the polymorphic transition (D tr H i ) should be included in the above expression. On the other hand, due to the negligible difference in the heat capacity of the two solid phases of the same compound, it is justified to omit the distinction of the C S P i of the two polymorphic forms.
Only in the case of the {C14 + C17} system there is noticeable difference between the calculated and experimental values that exceeds the measurement error, and this is a further confirmation of the presence of solid solutions in this system (Fig. 7). In all cases, the actual values of melting enthalpy are lower than the calculated ones, indicating very small negative values of the excess enthalpy of mixing (DH E ).

Density and viscosity
The density (r) of the EMs has been measured in the temperature range from 283.15 to 353.15 K and for all samples it reaches similar values, typical of paraffins, i.e. 0.76 g cm À3 at 298.15 K. While the density of pure n-alkanes increases with an increase in molecular weight, the density of their eutectic mixtures drops with the difference in length of the aliphatic chains, i.e. the (C14 + C17) EM has the highest density and the (C14 + C21) EM exhibits the lowest. This relationship is a direct result of the C14 content of the investigated mixtures (Fig. 8).
The experimental points included in Table 7 were correlated using eqn (3), the parameters of which are included in Table 9.
The coefficient of isobaric thermal expansion (a P ) of the tested mixtures assumes values characteristic of paraffins (a P = 9.5 Â 10 À4 K À1 ). The relationships observed for the density measurements, are also fulfilled by the dynamic viscosity (Z) of the tested ePCMs presented in Fig. 9 and Table 8. The highest viscosity is exhibited by the (C14 + C17) EM and the lowest viscosity by (C14 + C21) EM, while the viscosity of pure C14 is lower than that of the above mixtures (black squares in the Fig. 9). The above differences are not significant from a practical point of view, but they allow the properties of other mixtures consisting of n-alkanes to be predicted.
The experimental data of dynamic viscosity were correlated using the following two-parameter De Guzman-Andrade equation. 31 ln Z= mPa s ð Þ ½ ¼A þ B T=K (4)
a lightweight material with a high thermal conductivity value is suspended in the PCM. Predominant among such are carbon nanomaterials having structured graphene planes. In order for the vibrations of the carbon network to be able to transmit energy over long distances, carbon nanotubes with a high length-to-diameter ratio were chosen. Nanofluids based on SWCNTs are frequently studied in the literature, while their dominant disadvantage is their price. A cheaper alternative seems to be EG, which can be easily produced on a large scale. Due to the very similar properties of all blends, (C14 + C17) EM was chosen to study composite systems. The macroscopic stability of the proposed nanofluids was first investigated. Mixtures containing different mass percentages of SWCNTs (0.1%, 0.3%, 0.5%, 1.0%), GIC (10% and 75%) and EG (0.5%, 5.0%, 10%) were maintained at 298.2 K after preparation according to the procedure described above. Table 10 contains the observations made before and after 50 days of testing. All nanofluids containing SWCNTs showed stability during the experiment. The GICs started to sediment almost immediately, only the use of a solid material as high as the liquid column guarantees the stability of the material. A more developed surface area and a lower density of EG allows the material to be stable under isothermal conditions using no less than 10 wt% EG, samples with a lower content are characterised by EG floating to the liquid surface. The appearance of the samples before testing (A), after 50 days at 298.2 K (B) and after 1000 cycles of melting and crystallisation (C) is illustrated in Fig. 10.
Another aspect investigated was the stability of the nanofluids under PCM operating conditions, i.e. at cyclically varying temperatures. The same samples after 1000 cycles of melting and crystallisation showed the stability of the SWCNTs nanofluids with a content above 0.5 wt%, and at lower concentrations solids were observed to float to the liquid surface. For the GICs, no differences were observed compared to isothermal conditions, and for all samples containing EG, complete separation of the carbonaceous material from the liquid was observed. The observations are only valid for the situation where free space is present above the mixture, while this is a necessary condition to protect against the negative effects of thermal expansion.
Considering the inhomogeneous nature of the suspensions, it is not possible to take a representative sample for DSC analysis after the composites stability testing. Therefore, only the liquid fraction of samples was taken to evaluate the stability of the thermal effect. On the example of the composite containing 1 wt% of SWCNTs, it was found that there was no influence of the series of 1000 melting and crystallization on the temperature and thermal effect of the transition. The differences in measured values do not exceed the error of the measurement method. No influence of variable parameters on the studied system results from the inert nature of alkanes, which are thermally stable in  the studied temperature range and do not react with the added carbon materials. In addition, both SWCNTs and EG show chemical passivity towards most materials. The suspensions containing 1 wt% SWCNTs and 10 wt% EG were selected out for thermal conductivity measurements and compared with (C14 + C17) EM. For all samples, the thermal conductivity at higher temperature takes on lower values and for EM these values are comparable to those of C14 (l C14 = 0.137 W (m K) À1 at 298.2 K and l C14 = 0.130 W (m K) À1 at 338.2 K). 32 An addition of SWCNTs at an amount as low as 1% increases the thermal conductivity by about 50%, while achieving the same effect with EG probably requires the use of around 6-8 wt%, since a 10% addition results in an increase of no more than 100% (Fig. 11). Apart from the matter of the costs of using these additives, an important aspect is the reduction of the heat storage density, which will be decreased to an imperceptible extent with a content of 1% SWCNTs (Table 11).

Conclusions
In the present work, three binary ePCMs containing a ntetradecane and an odd n-alkane with melting points around 4 1C (277.2 AE 2 K) are proposed. The melting temperatures are suitable for their use in the thermal stabilisation of sensitive materials such as tissues, foods and drugs. The heat storage density, and therefore the amount of heat that can be absorbed before the temperature of the material rises, exceeds 200 J g À1 in two cases. The formation of ePCMs allows the melting temperature to be adjusted while maintaining a high latent heat: (C14 + C19) EM has an enthalpy of melting of 217.4 J g À1 at T fus = 276.2 K, while in comparison, C14 with a similar melting point (T fus = 278.4 K) has a latent heat of 225.3 J g À1 . Among the three examined systems, only {C14 (1) + C17 (2)} forms a solid solution, which is due to the small difference in length of the aliphatic chains and therefore negatively affects the heat of fusion of such a mixture, which is 155.5 J g À1 . Thus, when designing an ePCM, it is necessary to take into account the possibility of various processes occurring in the solid phase causing modification of its physical properties and failure to fulfil the model assumptions.
While relying on the favourable effects of the existence of positive deviations from ideality in the designed ePCMs, it must be borne in mind that very significant deviations lead to the occurrence of a miscibility gap in the liquid phase, making it impossible to obtain the expected eutectic. Accordingly, the attempt to form ePCMs from n-alkanes and a,o-diols (both groups are distinguished by a very high latent heat value) failed due to the miscibility gap in the liquid phase occurring over a wide range of compositions.
The application of carbon nanomaterials like SWCNTs allows for a significant improvement in the conductivity of ePCMs, which is a prerequisite to ensure efficient response to temperature changes. SWCNTs suspensions in alkanes are stable under isothermal conditions at nanomaterial contents as low as 0.1 wt%, while the negative effect of cyclic melting and crystallisation on the stability of nanofluids ceases to be noticeable at SWCNTs contents above 0.5 wt%. The flame retardant properties of GICs unfortunately do not cooperate with the possibility of using expandable graphite as a heat transfer medium due to the large density difference. The stability of such a material is only possible after the formation of a paste in which the GIC is saturated with PCM. The behaviour of nanofluids with EG is similar, as the low density of porous EG causes it to flow onto the surface of the molten PCM. Of the additives tested, SWCNTs have the highest ability to increase thermal conductivity.

Conflicts of interest
There are no conflicts to declare.