The Impact of Additives on the Main Properties of Phase Change Materials

: The main drawback of phase change materials (PCMs) is their low thermal conductivity, which limits the possibilities of a wide range of implementations. Therefore, the researchers, as found in the literature, proposed several methods to improve the thermal conductivity of PCMs, including inserting high thermal conductivity materials in nano-, micro-, and macro-scales, as well as encapsulation of PCMs. However, these inserts impact the other properties of PCMs like latent heat, melting temperature, thermal stability, and cycling stability. Hence, this paper aims to review the available in the open literature research on the main properties of enhanced PCMs that undergo solid–liquid transition. It is found that inserting high thermal conductivity materials and encapsulation results in improved thermal conductivity of PCMs, but it decreases their latent heat. Moreover, the insertions can act as nucleating agents, and the supercooling degree can be reduced. Some of the thermal conductivity enhancers (TCEs) may prevent PCMs from leakage. However, some test results are inconsistent and some seem to be questionable. Therefore, this review indicates these discrepancies and gaps in knowledge and points out possible directions for further research.


Introduction
Globally, in 2017, the total primary energy supply (TPES) was 13,972 Mtoe, of which about 81% was from fossil fuels like coal, oil, and natural gas [1]. Burning fossil fuels results in the release of harmful substances like greenhouse gases, including CO2, which contribute to climate change. Moreover, fossil fuel resources are limited, and they must also be converted to a usable form of energy, which is usually related to heat generation. On the other hand, according to the International Energy Agency [1], the share of renewable energy sources (RES) in the TPES is growing. The main disadvantage of RES is their unpredictability, resulting in a mismatch between demand and supply. Due to the abovementioned reasons, fossil fuels should be saved through increasing the efficiency of processes in the energy supply chain, and reduction the abovementioned mismatch. One of the solutions, or at least the minimization of these problems, is to implement thermal energy storage (TES).
There are three main types of TES, depending on the occurred phenomenon: sensible, latent, and thermochemical [2]. In sensible thermal energy storage (STES), heat is stored by increasing the temperature of a material. In STES, solid materials like rocks, metals, concrete, sand, bricks, and liquid materials like water or ethanol are usually used [3]. Latent heat thermal energy storage (LHTES) comes from isothermal or near isothermal heat absorption and release during the phase change of a material [4]. Materials that store latent heat are called phase change materials (PCMs). The third type Solid-liquid transition is the most commonly used type because of its large heat storage density and relatively small volume changes compared to solid-gas or liquid-gas transitions in which occurring gas occupies a large volume [13]. Solid-solid PCMs had not gained much attention yet mainly due to their low phase transition enthalpies [13], so this paper focuses only on solid-liquid PCMs and in the following parts of paper the term PCM will stand for only solid-liquid PCMs. The second classification of PCMs, based on chemical composition, divided them mainly into organic and inorganic due to their different chemical structure and various properties, which can be found elsewhere in the literature [10,12]. Required properties of PCMs include for example high latent and specific heat, high thermal conductivity, high density, small volume changes, suitable phase change temperature, low supercooling, good thermal stability, and good cycling stability [10][11][12]14]. Good cycling stability means no change in the PCM properties after many melting and solidification cycles. Detailed data on PCMs properties can be found elsewhere in the literature [7,8,10,12,15].
Due to the possibility of holding a large amount of heat, phase change materials may be used in many different applications such as maintaining thermal comfort in buildings by incorporating PCMs into floors, walls or ceilings [16], thermal energy storage of waste heat [17], thermal energy storage in domestic hot water systems [18], solar applications [19][20][21], thermal management of batteries [22], and electronic devices [23], incorporating PCMs into textiles to maintain constant body temperature [24], increasing cooling performance of air-conditioning units [25], the food industry [26], medicine [27] and many others which can be found in [28][29][30][31]. Although PCMs have many advantages over other TES technologies, their main drawback is low thermal conductivity, which in the solid state varies from approximately 0.6 W/(m•K) to 1.3 W/(m•K) for inorganic PCMs, and from approximately 0.2 W/(m•K) to 0.7 W/(m•K) for organic PCMs [8,12]. The low thermal conductivity of PCMs limits in practice the possibilities of a wide range of implementations. However, it is worth mentioning that there is a group of promising PCMs, i.e., metal alloys, with relatively high thermal conductivities [32], but due to their high thermal conductivity, they are not discussed in this paper. Melting temperatures, latent heat of fusion, and thermal conductivities of selected PCMs are presented in Table 1. Table 1. Properties of selected PCMs [8,12].

PCM
Melting Temperature, Due to the low thermal conductivity of PCMs, several methods of heat transfer enhancement in LHTES were proposed. The equation of heat transfer rate in the simplest form can be written as: where Q̇, U, A, and ΔT are the heat transfer rate (W), overall heat transfer coefficient (W/(m 2 •K)), surface area (m 2 ), and temperature difference (°C), respectively. According to Equation (1), the heat transfer rate can be improved by increasing the overall heat transfer coefficient, surface area, or temperature difference. Overall heat transfer coefficient depends, among others, on the thermal conductivity of a material. Thus, heat transfer enhancement techniques in LHTES can be divided into four main groups: 1. improving thermal properties, especially the thermal conductivity of a material, 2. increasing the heat transfer surface area, 3. improving the heat transfer process, 4. combined (hybrid) techniques.
Increasing the thermal conductivity of PCMs can be achieved by inserting high thermal conductivity materials into PCMs. Increasing the heat transfer surface area can be achieved by using fins or encapsulation PCMs. The third method includes using heat pipes, multiple PCMs, modifications of heat exchanger constructions, or changing process conditions, for example, by increasing heat transfer fluid (HTF) temperature and/or flow rate. Heat pipes can improve the heat transfer between HTF and PCM due to the evaporation and condensation of a working fluid, which results in a high heat transfer coefficient, thus augmenting the heat transfer between the HTF and PCM. Like was mentioned, the heat transfer rate between the HTF and PCM depends, among others, on the temperature difference between them. However, as the HTF flows through the heat exchanger, its temperature decreases, thus the temperature difference between the HTF and PCM also decreases. Therefore, to overcome the problem of decreasing temperature difference, multiple PCMs, i.e., PCMs with decreasing melting temperatures located along the HTF flow path, which ensures the constant temperature difference between the PCM and HTF. Constructions of heat exchangers can be modified in different ways, for example, Kadivar et al. [33] designed the double-tube heat exchanger with a non-concentric inner tube, which enhanced the natural convection and reduced the melting time of PCM. Combined techniques consist of a combination of at least two different techniques listed above.
The proposed classification of heat transfer enhancement techniques in LHTES is presented in Figure  2. In many cases, it is difficult to assign a heat transfer enhancement method in just one category since some techniques influence at least two parameters. For example, embedding PCMs into high thermal conductivity metal foam leads to both an increase in the effective thermal conductivity of composite material and an increase of the heat transfer surface area. In this paper, embedding PCMs into metal foam is included in improving the heat transfer conductivity group because most research on metal foam focuses only on thermal conductivity enhancement and neglects the fact of an increased heat transfer surface. A similar situation occurs in PCM encapsulation because encapsulated PCMs have a significant heat transfer surface area, and the shell of a capsule can be made of a high thermal conductivity material. Encapsulation is qualified for increasing the heat transfer surface method because shells are not always made of a high thermal conductivity material, but they always have an increased heat transfer surface area. Nevertheless, there is some research on encapsulation, which focuses on improving the thermal conductivity of shell materials, and such research is also discussed.
A detailed analysis of the mechanism of thermal conductivity enhancement of the enhanced PCMs can be found elsewhere [34]. Additionally, Liu et al. [35] and Rehman et al. [36] gathered the available correlations used to calculate the thermal conductivity of PCMs with thermal conductivity enhancers (TCEs).
There are several papers on heat transfer enhancement methods available in the literature [36][37][38][39][40]. Qureshi et al. [37] and Lin et al. [38] reviewed the addition of substances with high thermal conductivity and encapsulation as thermal conductivity enhancement techniques. Rehman et al. [36], in their work, focused only on using metallic and carbon-based porous materials to enhance thermal conductivity. Eanest Jebasingh and Arasu [39] and Leong et al. [40] reviewed research on nanoparticle insertion into PCMs. The abovementioned papers focus on either only one heat transfer enhancement technique or only thermal conductivity enhancement without considering how the applied method impacts other PCM properties. Therefore, there is a necessity to conduct a more comprehensive analysis of the impact of TCEs on the main properties of PCMs, which has been carried out in this paper.
This paper aims to review, discuss, and compare the research carried out, mainly between 2010 and 2020, on the main properties of enhanced PCMs, including thermal conductivity, latent heat, melting temperature, specific heat, thermal and cycling stability. The main purpose of this paper is to review a possibly wide range of thermal conductivity enhancers types and their influence on the main properties of the most commonly used PCMs, however, without emphasizing the types, i.e., organic or inorganic, of the PCMs. The effect of encapsulation and addition of the high thermal conductivity materials on the nano-, micro-, and macro-scale into PCMs to improve the PCMs' properties is discussed. In each section, the influence of the TCEs on the thermal conductivity, latent heat, specific heat, melting temperature, and supercooling degree of the PCMs, as well as the other features, e.g., phase change time, or additional benefits of enhanced PCMs are discussed. Despite the numerous papers on heat transfer enhancement in LHTES, to the best of the authors' knowledge, none of them considers the impact of different TCEs on properties other than the thermal conductivity. Moreover, some missing information in the knowledge, as well as possible directions for further research, is pointed out. This is expected to give a more comprehensive insight into the impact of TCEs on PCMs' properties.

Nanomaterials
Nanomaterials, i.e., materials that possess at least one dimension less than 1 μm, have gained much attention in the scientific community since the past few decades. They have found an application in LHTES as well. Nanomaterials are added to PCMs to form composite nano-enhanced phase change materials (NEPCMs) with improved thermal conductivity. The following types of nanomaterials are used to prepare NEPCMs: metals, metal oxides, and carbon-based materials. Generally, their thermal conductivity is higher than the PCMs, and therefore they are treated as TCEs. Although nanomaterials can increase the thermal conductivity of PCMs, they also affect the other properties of the PCMs, which is overlooked in the vast majority of the papers. Therefore, there is a necessity to conduct a more complex analysis of the properties of NEPCMs, which has been undertaken in this paper.
Before discussing the properties of NEPCMs (or the other composite materials in Sections 3 and 4), it should be noted that the addition of nanomaterials (or the other thermal conductivity enhancers) into PCMs does not change the PCMs' intrinsic properties, e.g., thermal conductivity, latent heat, etc. because the chemical composition of the PCMs remains unchanged. However, composite PCMs can be treated as homogenous materials with effective properties, which means that the properties of such composite materials are somehow a combination of the properties of pure PCMs and additives. Therefore, each time properties of composite materials are presented in this paper, it is referred to as their effective properties.
Thermal conductivity, latent heat, and melting temperature of NEPCMs are usually investigated by the researchers. A summary of such results is shown in Table 2. Thermal conductivity enhancement was calculated according to Equation (2): where n, keff, and kPCM are the thermal conductivity enhancement, the thermal conductivity of NEPCMs, and the thermal conductivity of pure PCM, respectively. Generally, NEPCMs possess higher thermal conductivity than pure PCMs. Maximum improvement, up to 3.36 times, was obtained by Wang et al. [41], who added expanded graphite (EG) into the octanoic acid/myristic acid mixture. In the case of metal-based and metal oxide-based nanomaterials maximal improvement was 1.86 times for paraffin wax [42] and 1.79 times for paraffin wax [43], respectively. The improvement of thermal conductivity depends mainly on the mass or volume fraction of the nanomaterial added into a PCM (detailed data about the PCMs and their thermal conductivity enhancers discussed in this section are shown in Table 2). However, the results of the research on this issue are inconsistent. Some researchers indicated that the thermal conductivity increases linearly with increasing nanoparticles fraction [44][45][46][47][48]. While others showed that the progression is nonlinear and the improvement starts to drop after reaching a certain value [42,[49][50][51][52][53]. The results obtained by Soni et al. [44] are in the first group, but they were achieved by using numerical calculations. Hence, there is a necessity to conduct experimental studies to confirm the results. He et al. [54] found in their experimental investigation that the thermal conductivity of NEPCMs increased linearly in a liquid phase. In a solid state the improvement for nano-graphite (NG)/PCM was also linear, but for graphene nanoplatelets (GNP)/PCM and multi-walled carbon nanotubes (MWCNT)/PCM there was saturation tendency, caused by agglomeration or sedimentation of the nanoparticles. The agglomeration of the nanoparticles was also reported in [48,49,51,54]. As a result, the thermal conductivity of NEPCMs with high fraction of the nanoparticles started to decrease. Therefore, dispersion stabilizers are used, like: sodium stearoyl lactylate [43], sodium dodecyl benzene sulfonate [46], polyvinyl pyrrolidone [52,55], and oleoyl sarcosine [56]. Some studies [41,44,49,51,52,54,57] investigated the insertion of different types of the nanomaterials into one PCM. The results of [44,52,57] indicated that the thermal conductivity of NEPCMs increased with increasing thermal conductivity of the nanomaterials. On the other hand, the results obtained in [41,49,51,54] do not confirm the abovementioned one, because the thermal conductivity improvement was not the highest for the highest-thermal conductivity nanomaterial. This indicates, that the improvement depends not only on the thermal conductivity of the nanoparticles, but also on their effective concentration [49], the nature of fractal structures, aggregation dynamics, and variation in Kapitza resistance [51], as well as the method of preparation [56].
The second important parameter of the PCMs is latent heat, which might be influenced by nanoparticles addition. As a result of nanoparticles addition, a part of the PCM is replaced by the non-phase-changing material and less latent heat per NEPCM mass can be stored. Such phenomenon was confirmed by many researchers [41,42,44,[46][47][48][54][55][56]58]. Based on Table 2, it can be concluded, that the decrease of the latent heat varied from about 2% to 15%, depending on the nanomaterial type and its mass fraction. The largest reduction of latent heat, by about 68%, was obtained by Ghossein et al. [56] who added silver nanoparticles (10 wt.%) into the eicosane. Theoretically, the latent heat of NEPCMs should decrease linearly with an increasing volume fraction of nanomaterials [41,44,46,55,56], which was confirmed experimentally by some researchers [41,47]. However, Masoumi et al. [46], Liu et al. [55], and Ghossein et al. [56] reported that the measured value of latent heat decreased more than predicted from the theory. According to Masoumi et al. [46], it was caused by the morphological, dimensional, and structural properties of the nanomaterials, while Ghossein et al. [56] attributed that finding to the colligative properties of the NEPCM.
On the other hand, the results on the latent heat of NEPCMs reported by some researchers [43,45,53,59] are significantly inconsistent with the abovementioned one, due to the fact that the NEPCMs' latent heat was higher than the latent heat of pure PCMs. Martin et al. [45] found that the latent heat of capric acid/SiO2 and capric-myristic acid/SiO2 increased by 8.0-10.7% and 6.8%, respectively, compared to pure PCM. Moreover, in the case of capric acid/SiO2 the improvement was proportional to the SiO2 mass fraction. Similar results were achieved by Zhu et al. [59], who prepared paraffin with aligned carbon nanotubes (a-CNT) and disordered carbon nanotubes (r-CNT). In the case of PCM/r-CNT the latent heat of NEPCM decreased, but for PCM/a-CNT an improvement by 11.2% was obtained, compared to pure PCM. The improvement of latent heat was proportional to the mass fraction of the nanoparticles. Vivekananthan and Amirtham [53] showed that the latent heat of melting increased from 311 kJ/kg to 338.6 kJ/kg, but the latent heat of solidification decreased from 308.2 kJ/kg to 252.8 kJ/kg for erythritol/GNP. Sami and Etesami [43] revealed that when the mass fraction of the TiO2 increased to 3%, the latent heat of NEPCM also increased. However, further increasing the fraction of the TiO2 resulted in decreased latent heat, compared to pure PCM. It was reported that increased latent heat of the NEPCM was caused by some interactions between nanoparticles and PCM molecules [43,59], however, it is not further discussed in this paper, and the readers are referred to [43,59] for details. Such abnormal results on latent heat like the abovementioned ones are in contrast with the findings reported previously in the literature and with theoretical predictions. Therefore, it seems that they should be verified by retesting, and then possible explanations of such results could be given.
A more practical approach to investigate the thermal properties of NEPCMs is an analysis of their phase change time, which was reported in [41,44,52,57,60]. Phase change time of NEPCMs depends mainly on their thermal conductivity and latent heat. In the same operating conditions, the phase change time of NEPCMs will decrease with increasing thermal conductivity and decreasing latent heat, because of the two reasons. Firstly, as thermal conductivity increases, the heat is transferred into the PCM faster. Secondly, if the latent heat of material decreases, less heat is required to melt the PCM, thus, for constant heat flux, the PCM melting rate increases with decreasing latent heat. Soni et al. [44] numerically investigated the solidification time of the NEPCMs. In the first case of their study, they assumed that solidification is a conduction-dominated process, whereas in the second case, both conduction and convection were included. The solidification time for pure PCM in the first case was 18,150 s, while in the second case, it was 19,000 s. This shows that in the case of numerical investigations, the results obtained depend on the assumptions and the models used during calculations (e.g., the conduction or conduction-convection model of heat transfer). Therefore, each result obtained by numerical calculations should be checked experimentally. Nevertheless, in research carried out by Soni et al. [44], the solidification time of PCM with Cu, Al, SiO2, and TiO2 nanoparticles, including both convection and conduction, was shorter by 9.7%, 6.8%, 2.1%, and 5.2% respectively, compared to pure PCM. Prabakaran et al. [60] reported that the solidification time of PCM/GNP composite was reduced by 39.21%, compared to pure PCM. Similar results were obtained experimentally by Choi et al. [52], who concluded that the heat transfer rate during solidification increased with increasing concentration of graphite (average thickness of 5 μm), and graphene (average thickness of 7 nm) in the NEPCMs. However, in the case of MWCNT/PCM, the heat transfer rate started to decrease when the volume fraction of MWCNT was greater than 1%, due to the high viscosity of the NEPCM and slowed natural convection. Wang et al. [41] found experimentally that the solidification time of EG/PCM composite was 9.41 times shorter compared to pure PCM, but the melting time of NEPCM was longer than the pure PCM's, due to slowed natural convection. The opposite results were achieved by Gupta et al. [57], who investigated the melting and solidification time of Fe/PCM and Cu/PCM nanocomposites. Compared to pure PCM the melting time was reduced by 7.8% and 5.6% for Fe/PCM and Cu/PCM, respectively. The solidification time of Fe/PCM and Cu/PCM was reduced by 35% and 30%, respectively, compared to pure PCM. It is worth noting that despite the Cu/PCM had slightly higher thermal conductivity, its phase change time was longer than the Fe/PCM, although the nanoparticles fraction was equal in both cases.
Another property of NEPCMs, which influences the possibility of energy storage, is their specific heat. It should be noted that there are limited studies on this property. Soni et al. [44] found that the specific heat of NEPCM decreased linearly with increasing nanoparticle volume fraction. Moreover, the greatest decrease was noticed for copper nanoparticles (by 4.3%), while the smallest decrease was obtained for aluminum nanoparticles (by 3.3%) compared to the other nanoparticles analyzed by Soni et al. [44]. In contrast, Martin et al. [45] showed in experimental research that the specific heat of capric acid/SiO2 and capric-myristic acid/SiO2 increased by 19.0-22.0% and 20.4-23.5%, respectively, depending on the SiO2 mass fraction. Similar findings were achieved by Aslfattahi et al. [58]. The specific heat of paraffin/Ti3C2 increased by 43% (compared to pure PCM) with an increasing mass fraction of the nanomaterial. According to the researchers it was caused by the interactions between the nanoparticles and the PCM.
The melting temperature of NEPCMs is generally only slightly different than the melting temperature of pure PCMs. Based on the information in Table 2, the abovementioned difference is less than 3 °C. In most cases, the melting temperature of NEPCM is a little lower than the pure PCM [41][42][43]46,47,55,56,59], but the opposite results can also be found [53,58]. Moreover, the nanoparticles can act as nucleating agents and decrease [42,48,53,55] or even eliminate [54] the supercooling degree, which is an advantage.
The next features of NEPCM are density and viscosity. These features are very rarely investigated by the researchers. Like was said in Section 1, a high density of PCMs is required, especially when the volume of the LHTES system is limited [40]. Based on the numerical model Soni et al. [44] showed that the density of NEPCMs increased linearly with an increasing volume fraction of the nanoparticles. Compared to the pure PCM, from the four investigated NEPCMs the greatest increment of density (by 30%), and the smallest increment (by 2.5%) was obtained for Cu/PCM, and for Al/PCM, respectively. Viscosity influences the behavior of NEPCMs, for example, high viscosity might suppress natural convection during the melting process. The viscosity of NEPCMs increases as the mass fraction of the nanoparticles increases [45,50,60]. In studies carried out by Á guila et al. [50] the viscosity of NEPCM increased by 61%, compared to pure PCM.
NEPCMs are expected to change their phase periodically in TES systems. Therefore, they should possess good cycling stability, which means that their properties should not change after many melting/solidification cycles. A slight decrease of the thermal conductivity [41,48,49,51], melting temperature [41,43,45,46,48,55], and latent heat [41,43,45,46,54,55,59] were observed after phase change cycles among the various literature reviewed. Martin et al. [45] showed that after 2000 cycles the latent heat of NEPCM decreased by 2.03%. Sami and Etesami [43] found that after 80 cycles the dispersion of the nanoparticles decreased. Nevertheless, some studies [48,53] confirmed that the changes in the properties of pure PCMs were higher than the changes of NEPCMs' properties. This confirms that the cycling stability of the NEPCMs can be better than pure PCMs. Ranjbar et al. [61] showed that after 500 cycles, the chemical structure of NEPCM did not change. It should be noted that among the reviewed literature, there are no chemical reactions in NEPCMs between the nanoparticles and the molecules of PCM [46,48,49,51,[53][54][55]57,59,61].
Another benefit of NEPCMs can be their higher decomposition temperature, compared to pure PCM, which was proven by some researchers [42,43,45,46,58,59]. The decomposition temperature can be improved by 52 °C, which was obtained by Zhu et al. [59]. This indicates better thermal stability of the NEPCM than the pure PCMs. However, the opposite findings were also achieved [48,54], and the reduction of decomposition temperature might restrict potential applications of NEPCMs. On the other hand, some NEPCMs possess the ability to prevent PCM from leakage [59,61]. The lack of leakage means that the PCM after melting is still kept in the composite material and no loss of composite material's mass is observed. Such ability allows applying PCMs directly, for example, the NEPCMs without leakage can be placed directly on walls to maintain the constant temperature in the room, without the necessity of placing them into tanks or containers.
Nanomaterials are promising candidates for improving the properties of PCMs because they might provide many advantages. From the discussed literature in Section 2 and the data collected in Table 2, it can be concluded that nanomaterials can increase the thermal conductivity of PCMs about 3 times, but they decrease the latent heat, and generally, the latent heat of NEPCMs decreases with increasing thermal conductivity. Nevertheless, nanomaterials can decrease or even eliminate the supercooling degree, improve thermal stability, and cycling stability. However, there are still some issues to be investigated. Firstly, there is a lack of complex studies that investigate all the properties of NEPCMs. Most of them focus only on one aspect, mainly the thermal conductivity, and the other properties are neglected, while the fact of decreased latent heat (usually up to 15%) of the NEPCMs, which was found in the majority of research, seems to be crucial in LHTES systems. Latent heat determines heat storage density, and therefore the optimum between thermal conductivity enhancement and latent heat reduction has to be found. However, the cases in which the latent and sensible heat of NEPCMs was improved require further studies to explain such results. Most of the researchers do not investigate the phase change time of NEPCMs. Like it was said before, increased thermal conductivity not always leads to the faster melting process. In many cases, researchers presented the value of thermal conductivity in one state only and/or they did not specify if the given value refers to the liquid or solid state of the material. Furthermore, sometimes researchers reported only the value of thermal conductivity enhancement, but they did not explain how they calculated it. For example, Mishra et al. [51] defined the thermal conductivity enhancement as: where n, keff, and kPCM are the thermal conductivity enhancement, the thermal conductivity of NEPCMs, and the thermal conductivity of pure PCMs, respectively, but in [45] the thermal conductivity enhancement was calculated as: With n, keff, and kPCM as before. There are also some problems with agglomeration and/or sedimentation of the nanoparticles. Solving this issue could be the subject of further studies. Additionally, the cost of the NEPCMs is overlooked. Only Mishra et al. [51] suggested that carbon black nano-powder (CBNP) might be one of the best candidates to be a TCE, due to its relatively low price of $1/g, compared to MWCNT, single-walled carbon nanotubes, and EG price of $5/g, $75/g, and $100/g, respectively. Martin et al. [45] suggested that the price of NEPCMs is high and there is a necessity to develop low-cost NEPCMs production methods on a large scale.
The results of such research are shown in Table 3, in which the melting and/or solidification time improvement was calculated as: where timp, tPCM, and tCPMC are melting (or solidification) time improvement, melting (or solidification) time of pure PCMs, and melting (or solidification) time of CPCMs, respectively. The research that investigated the thermal properties of CPCMs is summarized in Table 4, in which the thermal conductivity improvement was calculated according to Equation (2). Some experimental research showed that the thermal conductivity of paraffin, KNO3-LiNO3-Ca(NO3)2, and paraffin wax can be increased about 25 times [77], 50 times [80], and 45 times [63], when the carbon fibers, EG, and copper foam were added into the abovementioned PCMs, respectively. Wang et al. [64] estimated that the thermal conductivity of PCM/aluminum foam was 220 times higher than the pure PCM's. However, this value was calculated theoretically, and this prediction was not confirmed experimentally. Similarly, Fleming et al. [68] used a semi-analytical approach to estimate the thermal conductivity of PCM/aluminum foam composite.
The thermal conductivity of a CPCM depends on a few factors, mainly on the type and a mass or volume fraction of an additive. Various research showed that the thermal conductivity of CPCMs increases linearly with an increasing fraction of additive [62,78,[80][81][82][83]89]. However, numerous studies concluded that the abovementioned relation was nonlinear [69,78,79], and the thermal conductivity enhancement showed a saturation tendency [84][85][86], which means that the increase in the thermal conductivity was getting smaller as the fraction of additives increased. Nevertheless, based on the literature reviewed, it seems that the linear relation between the thermal conductivity enhancement and the fraction of TCE was observed mainly for carbon-based porous materials, although the non-linear relation for these materials was also reported. Therefore, it seems that the available results of investigations are insufficient to make general conclusions on this issue.
Secondly, the features of an additive may play an important role in thermal conductivity enhancement. The influence of carbon fibers' length was studied in [62,77,79]. Fukai et al. [77] showed that the influence of carbon fibers' length on the thermal conductivity improvement was insignificant, but the arrangement of carbon fibers was crucial. The fibers in the form of brush performed about 4 times better in terms of thermal conductivity enhancement than the fibers in the random arrangement. Frusteri et al. [62] found that the shortest carbon fibers ensured the highest thermal conductivity enhancement because their distribution in the PCM was the most homogenous. On the other hand, Zhang et al. [79] showed that longer carbon fibers resulted in greater thermal conductivity enhancement than the shorter ones. However, it should be noted that in [62,77] the length of the fibers was in the order of millimeters (the most advantageous length of carbon fibers was 0.2 mm [62]), while in [79] in the order of micrometers (the most advantageous length of carbon fibers was 0.225 mm [79]). Therefore, according to these studies, the carbon fibers' length of 0.2 mm seems to be optimal for thermal conductivity enhancement.
When porous materials are used to improve the properties of PCMs, porosity and pore density might affect the properties of CPCMs. According to [63,64,[67][68][69] porosity is the ratio of the volume of pores to the total volume of a material, usually given in percentages. Pore density, measured in pores per inch (PPI), can be defined as the number of pores in one linear inch. Usually, thermal conductivity improvement is inversely proportional to the porosity of porous materials [63,67,68,70,90]. In research conducted by Li et al. [67], the thermal conductivity of CPCM increased from 3.3 W/(m·K) to 6.8 W/(m·K) when the porosity of copper foam decreased from 97.3% to 92.4%. In research conducted by Zheng and Wang [70], the thermal conductivity of PCM/copper foam was 1.22 W/(m·K) and 7.07 W/(m·K) for the porosity of copper foam 98% and 93%, respectively. Xiao et al. [63] found that when the porosity of copper foam decreased from 96.95% to 88.89%, the thermal conductivity of CPCM increased from 5.04 W/(m·K) to 16.01 W/(m·K), and when the porosity of nickel foam decreased from 97.45% to 90.61%, the thermal conductivity of CPCM increased from 1.24 W/(m·K) to 2.33 W/(m·K). Li et al. [90] showed that reduction of the porosity of porous ceramics from 90% to 67% resulted in the improvement of thermal conductivity from 0.30 W/(m·K) to 0.51 W/(m·K). However, the results of research on the influence of pore density on the thermal conductivity enhancement are more inconsistent. Huang et al. [66] concluded that the thermal conductivity of CPCM decreased with increasing pore density of the copper, as well as nickel foam. On the contrary, Zhu et al. [71] found that as the pore density of iron foam increased, the thermal conductivity of CPCM also increased. Xiao et al. [63] concluded that the pore density did not affect significantly thermal conductivity improvement. Due to the inconsistency of the abovementioned results of experimental investigations, it is difficult to conclude unambiguously what effect the pore density has on the thermal conductivity enhancement. However, it seems that there should be an optimal pore density for which the thermal conductivity improvement is the greatest. Moreover, the measured value of the thermal conductivity of composite materials depends, among others, on the homogeneity of a tested sample, and the homogeneity of composite materials may vary significantly depending on porosity or pore density. Therefore, this issue could be a subject for further studies. Other features of materials that might affect the properties of CPCMs are the particle size of the material and its density. Two research teams found that thermal conductivity enhancement was greater when the size of TCE particles was smaller [84,91]. Additionally, Li et al. [91] concluded that thermal conductivity improvement was inversely proportional to the density of MgO, and the maximum thermal conductivity of CPCM was 5.81 W/(m·K). On the other hand, the opposite relation was found by Xiao et al. [86].
In addition to the abovementioned features of additives, the properties of CPCMs depend as well on the method of preparation [78] and a form of the prepared composite material [87]. Cheng et al. [78] used tetradecanol as PCM, expanded perlite (EP) as shape-stabilizing material, and copper powder, or carbon fibers as TCEs. From these materials, they prepared CPCMs by two different methods: mixing and implanting. In the mixing method, the TD and copper powder (or carbon fibers) were physically mixed, and then impregnated into the EP, whereas in the implanting method, the copper powder (or carbon fibers) were firstly implanted (added) into the EP and then the TD was impregnated into the previously prepared composite (copper powder (or carbon fibers)-EP). The results showed that the mixing method gave a better performance in terms of controllability and thermal conductivity enhancement compared to the implanting method. Ren et al. [87] made two CPCMs from Ca(NO3)2-NaNO3 and EG: one in granular form, and second in a shape-stable form prepared in cold compression and sintering process. The thermal conductivity of CPCM in the shapestable form was much higher (5.02-6.27 W/(m·K)) than in the granular form (1.56-1.90 W/(m·K)). Such phenomenon was the result of better connections between the particles of materials caused by the shape-stabilization process. To further improve the thermal conductivity of PCMs, some scientists proposed adding simultaneously two (or more) types of TCEs into one PCM. For example: iron foam and GNP [71], MgO and graphite flakes [91], high-density polyethylene (HDPE) and EG and MWCNT (or carbon nanofibers (CNF)) [88,94], and mullite and graphite flakes [92]. The combination of two or more thermal conductivity enhancers provides a synergistic thermal conductivity enhancement effect [88,94].
The main goal of inserting high thermal conductivity materials into PCMs is to increase the thermal conductivity of PCMs and hence to improve heat transfer rate and reduce charging (melting) or discharging (solidification) time. However, it should be noted that different researchers might define the charging and/or discharging time differently. For example, Yang et al. [72] defined the charging time as the time after which all PCM mass will be melted, but Liu et al. [73] and Righetti et al. [69] defined the charging time as the time after which temperature inside the LHTES unit will achieve a set value, i.e., the temperature set above the melting (or below the solidification in case of discharging) temperature of PCM. Moreover, the range of initial and end temperature of the PCM can be different in various research. Such discrepancies might make it difficult to compare the results obtained by different scientists. Therefore, the results given below in this paragraph are not compared to each other. Nevertheless, it has been proven in many experimental studies that inserting high thermal conductivity materials into PCMs results in shorter melting [63,66,[68][69][70]76,78,84,[86][87][88][89], and/or solidification time [69,73,74,77,80,81,83,88,89,92]. For example, the addition of copper, and alumina foam can decrease the melting time by 76.8% [74], and 93.6% [69], respectively, compared to pure PCM. The solidification time can be decreased by 72.2% [74] and 88.3% [69] for copper and alumina foam, respectively. Xu et al. [75] found that the melting time of PCM/copper foam can be reduced by 86.2%, but this result was obtained by numerical calculations. Hence, that finding should be checked experimentally. When TiO2 was added into n-octadecane, the melting and solidification time decreased by 45% and 26%, respectively [89]. Qu et al. [88] found that adding TCEs resulted in reduced melting and solidification time by 8.9-54.6% and 8.9-42.7% respectively, compared to pure PCM. The addition of EG can reduce the melting time by 73% [82]. These results indicate that the addition of the abovementioned inserts influences more the melting time than solidification time. However, Fu et al. [81] showed that the solidification time of the EG/PCM composite was reduced more than the melting time. A possible explanation for these discrepancies could be the different nature of the melting and solidification process, in which the convection and conduction are the dominant heat transfer mechanism, respectively. In the works discussed above in this paragraph, the different PCMs and TCEs were used, thus their thermophysical properties were different, which could affect the heat transfer, for example, it could suppress natural convection (as in the [81]).
Like was mentioned before, in the case of the metal foams application, thermal conductivity enhancement depends on their porosity and pore density. Righetti et al. [69] concluded that the pore density did not affect phase change time, but the melting and solidification time decreased with decreasing porosity of aluminum foam, which is consistent with the statement of the inversely proportional relationship between the porosity and thermal conductivity enhancement. The same result was reported by Chen et al. [74]. In contrast, Xu et al. [75] showed that the melting time of copper foam/PCM decreased with increasing porosity of the foam.
Although in the vast majority of research, phase change time of the CPCM decreases as its thermal conductivity increases, Fukai et al. [77] obtained different results. They found that the melting time of the CPCM with the volume fraction of carbon fibers less than 1% was longer compared to pure PCM. The reason for such a result was suppressed natural convection by the carbon fibers. Only when the volume fraction of the fibers increased to 2% the melting time of CPCM was the same as the pure PCM, because conduction heat transfer increased, and it compensated suppressed natural convection. On the other hand, the solidification time decreased with increasing volume fraction of the fibers. In addition to the already mentioned materials, the addition of other structures, like metallic meshes [95], aluminum honeycomb structures [96], aluminum or carbon fins [97] into PCMs can also be found in the literature. Another benefit of inserting high thermal conductivity materials into PCMs is a more uniform temperature distribution in the LHTES units [64,67,69,72,73].
The next desired property of LHTES systems is large heat storage capacity, which is directly related to latent heat and specific heat of the CPCMs. Based on data in Table 4 it can be concluded that the latent heat of CPCM is lower than the latent heat of pure PCMs. The decrease of a latent heat varies in a range of 1.4% [84] to 83.5% [89], depending on the type and fraction of additive. The latent heat of CPCMs decreases with increasing mass fraction of the additive [61,66,[67][68][69][70][71]74,75,81,84,86]. The specific heat of CPCMs is very rarely investigated by researchers. Xiao et al. [68] found that the specific heat of CPCMs decreased in both solid and liquid phases by 14.8-23.5% and 13.5-21.7%, respectively, compared to pure PCM. This was because the specific heat of additive was lower than pure PCM. On the other hand, Qu et al. [88] found that the specific heat of CPCMs in the liquid state decreased by 4.7-12.6%, but the specific heat in the solid state increased by 4.5%, compared to pure PCM.
Regarding the latent heat and the specific heat, some researchers investigated a more practical aspect, i.e., the heat storage capacity of LHTES units. Frusteri et al. [62] showed that the total heat storage capacity in the temperature range of 20-70 °C for pure PCM and CPCM was 292 kJ/L and 270 kJ/L, respectively. According to Liu et al. [73], the total heat storage capacity of the LHTES unit with CPCM increased by 3.1%, compared to the unit with pure PCM. Such improvement was caused by sensible heat stored in the metal foam. A similar result was obtained by Yang et al. [72]. It should be noted, that the results reported by Liu et al. [73] and Yang et al. [72] are possible only when the sensible heat stored in TCE is larger than the latent heat lost due to reduced mass of PCM.
The second practical aspect investigated by the researchers is the heat transfer rate, which can be improved by using CPCMs instead of pure PCMs [73,75,76]. In the LHTES system investigated by Liu et al. [73], the heat transfer rate increased by 49.3% and 10% during melting and solidification, respectively. Comprehensive numerical investigation of the LHTES with PCM locally enhanced with copper foam (Figure 3) was conducted by Xu et al. [75]. It was concluded that to achieve the lowest melting time simultaneously maintaining the heat storage capacity value of the LHTES system as high as possible (by reducing the mass of TCE), foam inserts should be gathered in the lower part of the annular region (Figure 3a). Compared to a case in (Figure 3d), the average heat transfer rate increased to 5.1, 6.8, and 6.9 times for cases (  (Figure 3a) and (Figure 3b)) might reduce the costs of metal foam. This is consistent with the conclusion made by Gasia et al. [76], who investigated four heat exchangers with: pure PCM, 17 aluminum fins, metallic wool distributed in a finned shape around the HTF tubes, and metallic wool distributed randomly around the HTF tubes. It was concluded that although the heat exchanger with fins gave the best performance, the advantage of metallic wool is its availability and low price. Lu et al. [80] estimated that the cost of CPCM consisted of KNO3-LiNO3-Ca(NO3)2 and expanded graphite (EG) was lower than the pure PCM, because the price of EG was lower than the price of PCM. Nevertheless, a majority of researchers neglect the cost of CPCMs. Phase change temperature is another important property of composite PCMs with inserts at micro/macro-scale, which defines their application suitability. Based on data in Table 4, it can be concluded that in most cases the melting temperature of CPCMs decreases slightly, by less than 4 °C [78,80,81,[83][84][85]89,90], which is an acceptable change. Nevertheless, much more reduction of melting temperature, by 18 °C, was also observed [82]. In some studies [66,68,71,79,88] a slight increase, from 0.96 °C [68] to 1.9 °C [79], of the melting temperature was detected. Usually, the changes in melting temperatures were independent of the mass or volume fraction of additive. However, some researchers found that the changes in melting temperatures were proportional to the fraction of additives [78], especially when the EG was used as TCE [80,82,83]. Different types of MgO (light and heavy, with densities of 0.2 g/mL and 0.5 g/mL, respectively) were used by Li et al. [91], who found that the melting temperature of CPCM with light and heavy MgO decreased by 3-5 °C and 6-9 °C, respectively, compared to pure PCM. Kenisarin et al. [84] used two types of EG with different sizes: 200-1200 μm and 50-200 μm. The melting temperature of CPCMs with large EG size, and with small EG size decreased by 1 °C, and 7-8 °C, respectively, compared to pure PCM. According to Kenisarin et al. [84], that difference was caused by some intermolecular interactions between PCM and the smaller EG. It should be noted, that only Kenisarin et al. [84] discovered some intermolecular interactions between the PCM and the additive. According to the authors, the nature of these interactions should be further investigated. The vast majority of the research found that PCMs and additives are physically mixed and no chemical reactions occur among them [57,61,[66][67][68][69][70][71][72]74,75,79,81,83,86].
Another advantage of CPCMs is their reduced supercooling degree, compared to pure PCMs [85,86,89,92]. For example, in research conducted by Xiao et al. [86] the supercooling degree was reduced from 13 °C (pure PCM) to 2.4 °C (composite PCM). The supercooling degree of CPCMs prepared by Qu et al. [88] was only 0.1 °C, but the researchers did not give the value of the supercooling degree for pure PCM. Fu et al. [81] and Yuan et al. [82] found that a larger reduction of supercooling degrees can be obtained by increasing the mass fraction of additive. On the other hand, Huang et al. [66] showed that the supercooling degree increased from 0.51 °C for pure PCM, to 2.75-3.50 °C for CPCMs with copper, and nickel foam. A similar finding was obtained by Xiao et al. [68].
The properties of CPCMs should not change after many melting-solidification cycles, which means that they should possess good cycling stability. Generally, after numerous melting/solidification cycles of the CPCMs listed in Table 4, their latent heat decreased [59,78,81,[84][85][86]89,90,92] in the range from 0.7% [78] to 19% [84] after 200 phase change cycles. The changes in melting temperature were usually referred to as negligibly small [78,81,[84][85][86]92]. However, Xiao et al. [86] found that the supercooling degree of CPCM increased after 400 melting-solidification cycles, but its value was still below 5 °C. Additionally, Li et al. [90] found that after 25 cycles the weight loss of CPCM was less than 10%, which indicated leakage of the PCM. A small leakage, i.e., the CPCM's mass decreased by 1.38%, of the PCM after 40 cycles was reported also by Yuan et al. [82]. Furthermore, Li et al. [91] observed a breakage of graphite flakes in the CPCM after 100 cycles.
Additives at micro/macro-scales should not impact the thermal stability of CPCMs, especially they should not reduce a decomposition temperature. Some studies confirmed that there were no noticeable changes in the decomposition temperature of CPCMs, compared to pure PCMs [66,71,83,85,89,91,92]. For example, the decomposition temperature of CPCM prepared by Cheng et al. [78] decreased by 2.9 °C, compared to pure PCM. It could be concluded, that such small change does not influence the working temperature range of CPCMs.
Adding materials at micro/macro-scale into PCMs has also an operational advantage. It can prevent PCMs from leakage due to the possibility of preparing shape-stable composite materials [57,66,[68][69][70]74,79,[82][83][84]86,95]. Depending on the type of additive, different mass fractions of additives must be used to prevent PCM from leakage. Yuan et al. [82], Zhang et al. [85], and Zhu et al. [59] prepared the erythritol/EG, n-eicosane/EG/SiO2, and paraffin wax/expanded vermiculite composite materials, respectively, and the maximum mass fractions of PCMs in these composite materials were 90%, 70%, and 60%, respectively. If the mass fractions of PCMs in the CPCMs exceeded the abovementioned values, a leakage of PCMs occurred. Another interesting property that CPCMs might possess is flexibility, which was reported by Wu et al. [83]. They prepared the flexible CPCM using paraffin wax, olefin block copolymer (OBC) and EG. Such flexible CPCM may considerably extend the scope of CPCMs use. On the other hand, the use of porous ceramics as TCE allowed preparing CPCMs with high mechanical strength [90,91].
Summarizing, the application of additives at micro/micro-scales can provide many benefits. Firstly, it can significantly (up to 50 times) improve the thermal conductivity of PCMs, and as a result, shorten phase change time (up to 93.6%). However, some studies proved that the additives may suppress natural convection during melting and therefore extend phase change time. It should be noted that the phase change time was not always investigated by researchers. Moreover, the charging/discharging time was defined differently by different researchers, which makes it difficult to compare. Another advantage of CPCMs is that they provide more uniform temperature distribution in LHTES units, and they possess good thermal and cycling stability. Generally, melting temperatures of CPCMs are very similar to those of pure PCMs, while the supercooling degree of CPCMs might be reduced and reach value as low as 0.1 °C. Additionally, the additives can prevent PCMs from leakage. On the other hand, the main negative effect of using additives at micro/macroscales is decreased latent heat of CPCMs (up to 83.5%). Hence, the heat storage density of a LHTES unit might be reduced. However, some studies showed that the heat storage density of a TES unit with CPCM might be higher compared to a unit with pure PCM, because of specific heat stored in the additive. The next important feature of CPCMs is their cost, which is overlooked by most of the researchers. Therefore, it can be generally concluded that there is a lack of studies that investigate every aspect of preparing CPCMs and all their most important properties.

Encapsulation
Heat storage in PCMs is connected with the melting process, and therefore, the presence of materials in a liquid state is inevitable. Hence, there is a necessity of locating PCMs in enclosures to prevent the material from leakage. One possible solution to the problem of leakage is the encapsulation of PCMs. In the literature, there are different definitions of encapsulation-for some researchers, encapsulation means a tank of any shape and size filled with PCMs, but for others, a capsule is a spherical shell with a PCM core, and this approach is the most commonly used [98]. Depending on the size, PCM capsules have been divided into three groups: macrocapsules, microcapsules, and nanocapsules, wherein dimensions of these capsules are more than 1 mm (but usually less than 10 cm), 1 µ m-1 mm, and less than 1 µ m, respectively [98]. PCM capsules, besides the fact that they prevent material from leakage, they have another advantage-a possibility to disperse them in liquid, e.g., water, to create PCM slurries, which can accumulate the heat and transfer the heat [99]. Usually, the heat capacity of such heat transfer fluids (HTFs), i.e., PCM slurries, is much greater compared to HTFs that use only sensible heat; furthermore, the thermal conductivity of PCM slurries may be higher than the thermal conductivity of HTFs that uses only sensible heat [100].
As mentioned in Section 1, encapsulation has been classified as one of the heat transfer enhancement techniques, mainly due to the large heat transfer surface area of capsules [101,102]. Nevertheless, the capsules' shells can be made of high thermal conductivity materials, which lead to the greater thermal conductivity of the whole capsule. Studies on encapsulation that focus on the high thermal conductivity shells are summarized in Table 5. The thermal conductivity of capsules can be improved up to 11 times, which was reported by Yu et al. [103], who prepared the capsule made of n-octadecane (core) and CaCO3 (shell). The improvement of thermal conductivity depends, among others, on the mass fraction of shell. If the thermal conductivity of shell is higher than the thermal conductivity of PCM, then the thermal conductivity of capsules increases with an increasing shell to core mass ratio [103,104]. The second factor that affects the thermal conductivity of the PCM capsule is a type of shell material and some of the researchers proposed to create hybrid shells [100,[105][106][107]. Such hybrid shells consist of at least two different materials. Xia et al. [106] fabricated the capsule with a melamine-formaldehyde (MF)/boron nitrate (BN) hybrid shell. The thermal conductivity of capsule increased almost linearly with increasing mass fraction of BN. Similar findings were obtained by Wang et al. [107], who created the capsules with the shell made of poly(melamine-formaldehyde) (PMF) and different mass fraction of SiC. Although the thermal conductivity of capsules increased with increasing fraction of SiC, some aggregates in the products occurred when SiC content increased to 7 wt.%. However, the application of SiC resulted in additional benefit, i.e., increased solar irradiation absorbance. As a result, the temperature of the hybrid-shell capsule was higher by 8 °C compared to the single-shell capsule after exposure to sunlight due to photo-thermal conversion. Such property could be utilized in solar energy applications. The capsules prepared by Zhu et al. [100] possessed the same property (photo-thermal conversion). The capsules were made of SiO2/graphene shells with different graphene mass fractions. Moreover, the thermal conductivity of capsules increased to about 1.5 W/(m·K) as graphene mass fraction increased to 2.5%. However, when graphene concentration further increased to 3.75 wt.%, and 5.0 wt.% the thermal conductivity of capsules decreased to 0.99 W/(m·K), and 0.97 W/(m·K), respectively. According to the authors, this abnormal phenomenon was caused by differences in the capsules' sizes, i.e., the size of the capsules with graphene mass fraction of 2.5%, 3.75%, and 5.0% was 390 nm, 257 nm, and 256 nm, respectively. Hybrid shells were also prepared by Zhu et al. [105]. The shells consisted of polystyrene (PS) and SiO2 or poly (hydroxyethyl methacrylate) (PHEMA) and SiO2. The capsules with PS/SiO2, and PHEMA/SiO2 shells possessed higher thermal conductivity than the capsules with SiO2 shell, which, according to Zhu et al. [105], was inconsistent with their predictions because the organic materials (PS and PHEMA) possessed lower thermal conductivities than inorganic material (SiO2). The researchers explained that result by a unique, more compact structure of the hybrid shells. Moreover, the mechanical properties of capsules with hybrid shells were superior compared to capsules' with SiO2 shell. Another way of capsules thermal conductivity enhancement was proposed by Praveen et al. [108], who did not create a hybrid shell, but mixed GNP with capsules. As a result, voids between the microcapsules were filled with GNP, which improved the heat transfer.
Increased thermal conductivity of capsules should reduce the phase change time of PCMs, which was examined by some researchers [98,101,104,105,107]. The charging and discharging time of capsules prepared by Zhang et al. [102] was reduced by 6.4%, and 29.4%, respectively, compared to pure PCM. Zhu et al. [105] found that the charging time of PCM/PS/SiO2, and PCM/PHEMA/SiO2 capsules was reduced by 23.4%, and 18.0%, respectively, compared to PCM/SiO2 capsule. The charging time reduction obtained by Wang et al. [107] was 11.9%, while the greatest charging time reduction, by 54.4%, was achieved by Zhu et al. [100].
The other method of phase change time reduction was proposed by Jia et al. [109], who investigated numerically the solidification process in the capsule with circular pin-fins, as shown in Figure 4. It was found that increasing the number of fins, their length, and diameter led to a reduction of solidification time by 60%. However, the reduction was insignificant after the number of fins, their length, and diameter exceeded 6, 15 mm, and 2 mm, respectively. Moreover, such operations resulted in decreased PCM mass by about 1.5%. Puertas et al. [110] used numerical simulations to examine the influence of the capsule's thermal conductivity on the melting time of PCM. They found that the melting time decreased by 20% when the thermal conductivity of the capsule increased 100 times. Similarly, Feng et al. [111] compared the phase change time of PCMs placed in tubular capsules with different thermal conductivities. It was concluded that the thermal conductivity of shell had a minor impact on the phase change time than the thermal conductivity of PCMs. Therefore, increasing the thermal conductivity of PCMs seems to be more advantageous than increasing the thermal conductivity of shells. Praveen et al. [108] investigated the thermal performance of heat sink without PCM (case I), with pure PCM (case II), and with capsules (case III). The heat was supplied to the heat sink until its (heat sink) temperature rose to 80 °C and it was assumed that the longer time to reach this temperature, the better was the thermal performance of the heat sink. It was found that the time to reach 80 °C was about 50 min, 88 min, and 62 min for cases I, II, and III, respectively. Therefore, despite the PCM capsules had higher thermal conductivity than pure PCM, their thermal performance was worse, because of lower latent heat. The low latent heat of PCM capsules is their main drawback. Based on information from Table 5, it can be concluded that in most cases the latent heat of PCM capsules decreased by approximately 30-50% [100,102,[105][106][107][108]112], compared to pure PCMs. However, the smallest and the largest reduction were 15.5% [104] and 77.6% [103], respectively. The latent heat reduction depends mainly on the core to shell mass ratio, and the latent heat decreases with decreasing value of that ratio [102][103][104]105,112]. However, Zhu et al. [100], Wang et al. [107], and Xia et al. [106] found that there is no relationship between the mass fraction of shell and the latent heat of capsule. Moreover, Xia et al. [106] showed that the increase of boron nitrate (BN) mass fraction in shell resulted in higher latent heat, compared to the capsule with shell made of only melamineformaldehyde (MF). According to the authors, that phenomenon was a result of increased thermal conductivity of the capsule which allowed the entire mass of PCM to change its phase.
Another property of encapsulated PCMs is their melting temperature. Generally, changes in capsules' melting temperatures are negligibly small, less than 1 °C [102][103][104][106][107][108]112]. However, in two studies [100,105] the melting temperature of capsules decreased by about 3 °C, independently from the shell to the core mass ratio. Additionally, the shell material can act as a nucleating agent, hence reduce the supercooling degree [100,103,107,112]. In contrast, Zhu et al. [105] obtained different results that indicated extended supercooling for capsules.
Some of the researchers investigated the thermal stability of PCM capsules and changes in their properties after melting-solidification cycles i.e., cycling stability. Zhu et al. [100], Yu et al. [103], and Cheng et al. [104] did not observe significant changes in the melting temperature and latent heat of capsules after 500, 200, and 20 melting/solidification cycles, respectively. Only slight fluctuations in the melting temperature and the latent heat of capsule were detected by Xia et al. [106], after 40 cycles. However, some of the studies showed that the latent heat of capsules decreased by 9% [105], and 19.7% [102] after 100 and 500 cycles, respectively. Regarding the thermal stability, it was determined, that encapsulation did not influence the thermal stability and decomposition temperature [102,105,108,112]. Moreover, in some studies [100,103,106,107], the decomposition temperature increased, because shell acted as a protective barrier. It should be noted that encapsulation does not change the crystal structure of the encapsulated PCM [100,102,105,112], and no chemical reactions between the core and shell materials occur [106,107].
Summarizing, PCMs' encapsulation might bring considerable benefits. It prevents PCMs from leakage and gives the possibility to prepare PCMs slurries, which can store, as well as transport heat. Generally, capsules filled with PCMs possess good thermal and cycling stability. Their melting temperature is comparable to that of pure PCMs, while a supercooling degree can be reduced. Their thermal conductivity can be 11 times higher than PCMs. Nevertheless, their main disadvantage is low latent heat, reduced by about 80% in comparison with pure PCM. On the other hand, the specific heat of capsules, as well as the cost of capsules preparation was overlooked by the researchers. Moreover, n-octadecane was a core material in the majority of the studies. Hence, there is an opportunity for researchers to create and investigate the properties of capsules with other types of PCMs as core materials.

Conclusions
In this paper, the research available in the open literature from 2010 to 2020 on the influence of TCEs on the properties of PCMs has been reviewed. The impact of various nanomaterials like: metals, metal oxides, and carbon-based nanoparticles has been studied. The addition of materials at the micro-and macro-scale such as carbon fibers, metal foams, expanded graphite, and the other materials have been also investigated. Additionally, the properties of encapsulated PCMs have been presented.
All of the TCEs might impact significantly the properties of PCMs. The properties depend on the type of TCE, its mass or volume fraction, as well as its physical properties like shape or size. A preparation method of enhanced PCMs is also important. Nevertheless, the TCEs can improve the thermal conductivity of PCMs. The greatest improvement (up to 50 times) can be achieved by using materials at the micro/macro-scale. PCMs encapsulation can improve the thermal conductivity up to 11 times. The lowest improvement is provided by nanomaterials (about 3 times). However, to achieve the same thermal conductivity improvement, up to 3 times, a much smaller mass fraction of the nanoparticles is required than the materials at micro/macro-scale. The reduction of the TCE amount leads to a lower decrease in latent heat. Therefore, the latent heat reduction of NEPCMs (up to 15%) is much lower than CPCMs (up to 83.5%) and capsules (up to 77.6%). Moreover, in some cases, the latent heat and the specific heat of NEPCMs was higher than pure PCMs which should be further investigated, because the reason for this phenomenon has not been thoroughly explained.
The addition of TCEs leads to higher thermal conductivity values, but some studies showed that it might be not enough to shorten the phase change time of PCMs, because TCEs may suppress natural convection. Therefore, the mass or volume fraction of the added materials should be chosen very carefully. However, in the majority of studies phase change time of enhanced PCMs was reduced, compared to pure PCMs. Scientists attributed it to the improved thermal conductivity, but it should be noted that decreased latent heat might also contribute to reduced phase change time.
Generally, melting temperatures of NEPCMs, CPCMs, and capsules undergo minor changes, while the TCEs or shells can act like nucleating agents, which reduces the supercooling degree. Enhanced PCMs possess good thermal and cycling stability, which means that their properties do not change after many thermal cycles. Moreover, the thermal stability of enhanced PCMs can be even improved. Another benefit of capsules is preventing PCMs from leakage and the possibility of preparing PCMs slurries, which possess the ability to store and transport heat. CPCMs can also prevent PCMs from leakage, and they allow to create shape-stable or flexible materials. Furthermore, the temperature distribution in TES units with CPCMs is more uniform than in units with pure PCMs.
Although researchers dedicate much effort to improve the thermal conductivity of PCMs, knowledge is still incomplete. Many issues should be investigated and could be the subject of future work. Many researchers focus only on thermal conductivity enhancement, and they do not investigate the influence of additives on other PCM properties. Furthermore, most of the research concerns organic PCMs with low melting temperatures (below 130 °C). Therefore, more inorganic PCMs and PCMs with higher melting temperatures should be investigated, and further research should include the influence of additives on all of the properties of enhanced PCMs. The economic aspect of preparing enhanced PCMs is also neglected. The influence of shape, size, and geometry of additives on properties of enhanced PCMs is not fully explored, and therefore requires more research. In most research the thermal conductivity of enhanced PCMs is measured. However, the time of phase change should also be investigated because some additives can suppress the natural convection, and despite increased thermal conductivity, the phase change time may be increased. Moreover, there is not enough comprehensive research that would cover the impact of TCEs on all of the properties of PCMs such as thermal conductivity, latent heat, specific heat, melting temperature, thermal stability, cycling stability, and its influence on the performance of LHTES units.
Based on the presented work, it can be concluded that although the greatest improvement of the thermal conductivity can be achieved by using TCEs at the micro/macro-scale and PCMs encapsulation, these additives and encapsulation may reduce the latent heat such significantly, that the CPCMs or PCM capsules will not fulfil their role as thermal energy storage materials. Moreover, the use of such composite materials could be associated with application problems, i.e., the necessity of creating thermal energy storage tanks with large volume and mass, which is unfavorable. Therefore, the addition of nanoparticles seems to be the most promising technique to improve PCMs' properties. A low nanoparticle mass fraction can considerably improve thermal conductivity with a relatively small decrease of latent heat. On the other hand, preparing shape-stable PCMs also seems to be a prospective direction of development because they do not require special tanks or containers, they can be made in any shape, and there are even attempts to create flexible ones, which can significantly expand the possibilities of their applications. Nevertheless, research on the preparing NEPCMs and shape-stable PCMs is relatively new, but this area of research is developing, and the number of studies on the NEPCMs and shape-stable PCMs is increasing every year. Therefore, this direction of research seems to be the most promising.