Journal of Energy Storage 25 (2019) 100874

Contents lists available at ScienceDirect

Journal of Energy Storage

journal homepage: www.elsevier.com/locate/est

Preparation, thermal characterization and examination of phase change T

materials (PCMs) enhanced by carbon-based nanoparticles for solar thermal
energy storage
⁎
Meizhi Hea,b, Luwei Yanga,b, Wenye Linc, , Jiaxiang Chena,b, Xiang Maoa,b, Zhenjun Mac
a
CAS Key Laboratory of Cryogenics, Technical Institute of Physics and Chemistry, Beijing, 100190, China
b
University of Chinese Academy of Sciences, Beijing, 100049, China
c
Sustainable Buildings Research Centre (SBRC), University of Wollongong, 2522, Australia

A R T I C LE I N FO A B S T R A C T

Keywords: This paper presents the preparation and thermal characterization of phase change materials (PCMs) enhanced by
Phase change material (PCM) carbon-based nanoparticles, including graphene nanoplatelets (GNPs), multi-walled carbon nanotubes
Graphene nanoplatelets (MWCNTs) and nano-graphite (NG). A systematic experimental framework, consisting of material selection and
Multi-walled carbon nanotubes preparation, material property characterization and thermal performance examination, was proposed and used
Nano-graphite
in this study to facilitate the development of nano-enhanced PCMs (NePCMs) for solar thermal energy storage
Thermal properties
Solar thermal storage
applications. By applying this framework, the characteristics and potential performance of PCM composites can
be comprehensively understood, and better assessed before practical applications. It was found that the thermal
conductivity of the myristic acid (MA) can be significantly enhanced by adding the nanoparticles in particular
GNPs as additives into the PCM. The thermal conductivity of the PCM composites can be improved by 176.26%,
47.30% and 44.01% respectively under the solid phase, by adding GNPs, MWCNTs and NG with a concentration
of 3 wt%. However, the concentration of the nanoparticles needs to be carefully determined to maximise the
benefit in thermal conductivity enhancement. Different from that under the solid phase, the thermal con-
ductivity enhancement of the NePCMs developed under the liquid phase followed linear increasing trends with
relatively low increasing rates, when increasing the concentration of the nanoparticles. Besides the thermal
conductivity enhancement, the adding of nanoparticles also modified the phase change process with a smaller
phase change temperature range and eliminated supercooling while maintaining the high latent heat capacity. A
further thermal performance examination demonstrated that the prepared NePCMs showed high thermal and
chemical stability, which can be used to substantially reduce the phase transition time, and therefore are good
potential candidates for solar thermal energy storage applications.

1. Introduction solution to rationalise the utilisation of solar thermal energy, compared

with sensible heat TES systems [7].
Due to increasingly excessive energy consumption and severe en- For proper functioning of a PCM TES system, rapid melting and
vironmental problems over the last two decades, the development and solidification of the PCMs are desirable to respond to fast load varia-
deployment of thermal energy storage (TES) systems using phase tions. However, most PCMs have the inherent disadvantage of low
change materials (PCMs) to facilitate energy savings, renewable energy thermal conductivity, which results in low charging and discharging
utilization, and efficiency enhancement of thermal systems have at- rates, and limits the overall energy storage capacity and efficiency of
tracted considerable attention [1–3]. PCMs are substances that can PCM TES systems [8,9]. To solve this problem, a number of heat
absorb, store and release a large amount of thermal energy at relatively transfer enhancement techniques have been proposed, such as the use
constant temperatures [4,5]. One of typical applications of PCM TES is of fins and heat pipes, embedment/dispersion or shape-stabilisation
for the storage of solar thermal energy which is considered to be the using high thermal conductivity materials, micro-encapsulation and
most abundant renewable energy [6]. The utilization of PCM TES can macro-encapsulation [10–14]. The dispersion of high thermal con-
overcome the intermittency of solar energy, thereby providing a better ductivity additives is one of the effective approaches to enhancing the

⁎
Corresponding author.
E-mail address: [email protected] (W. Lin).

https://doi.org/10.1016/j.est.2019.100874
Received 27 April 2019; Received in revised form 15 July 2019; Accepted 28 July 2019
Available online 10 August 2019
2352-152X/ © 2019 Elsevier Ltd. All rights reserved.
M. He, et al. Journal of Energy Storage 25 (2019) 100874

thermal conductivity of PCMs, however, normal millimetre-sized or 3.9 °C and 3.6 °C, respectively. Colla et al. [32] selected carbon black
micron-sized additives have the disadvantages of relatively low surface (CB) nanoparticles and paraffin to prepare nano-enhanced PCMs. It was
energy and relatively large density, which are more likely to cause se- found that CB based nano-PCMs showed high thermal conductivity
diments [15]. As an alternative solution, adding nanoparticles into enhancement of more than 25%. Wang et al. [33] dispersed treated
PCMs to prepare PCM composites with high thermal conductivity is carbon nanotubes (CNTs) into the palmitic acid (PA) for property
currently of great interest and receiving increasing attention [16–18]. modification. The thermal characterization showed that the thermal
For instance, Harikrishnan and Kalaiselvam [19] prepared a nano-en- conductivity of prepared treated CNTs/PA composite was 46.0% higher
hanced PCM composite using CuO nanoparticles as additives and oleic than that of the raw PA, when the concentration of treated CNTs was
acid as the base PCM. It was found that by adding 2 wt% of CuO na- 1 wt%.
noparticles, the thermal conductivity of the PCM composite can be This paper presents the preparation and thermal characterization of
improved significantly by 98.66%. As a result, the time needed for high-performance nano-enhanced PCM composites for solar thermal
solidification and melting processes can be reduced by 27.67% and applications using myristic acid (MA) and carbon-based nanoparticles,
28.57%, respectively. Zeng et al. [20] added copper nanowires into including GNPs, MWCNTs and NG. The main research idea of this work
tetradecyl alcohol to prepare nano-enhanced PCM composites. It was is to develop PCM composites with desired characteristics, especially
found that by adding 11.9 vol% copper nanowires, the thermal con- with high thermal conductivity, to facilitate the effective thermal en-
ductivity reached as high as 2.86 W·m−1· K−1, which was about 9 times ergy management of solar thermal systems. The novelty of this study is
that of the pure tetradecyl alcohol. Sahan and Paksoy [21] used syn- the development and utilisation of a systematic experimental frame-
thesized Fe3O4 nanoparticles as additives to prepare nano-Fe3O4/par- work to guide material preparation, properties characterisation and
affin PCM composites by a sol-gel method. It was interesting to find that performance examination of nano-enhanced PCMs (NePCMs).
the adding of nano-Fe3O4 can increase both thermal conductivity and Accordingly, the base PCM and different nanoparticles were selected to
thermal storage capacity. The Differential Scanning Calorimetry (DSC) develop various NePCM composites. The chemical and thermal prop-
test results showed that the latent heat storage capacity increased by erties of the prepared PCM composites were characterized by scanning
approximately 20% with the adding of 10 wt% nano-Fe3O4. Zeng et al. electronic microscope (SEM) analysis, Fourier transform infrared (FT-
[22] prepared nano-enhanced PCM composites by adding Ag nanowires IR) spectroscopy, thermal conductivity measurement, DSC test and
into 1-Tetradecanol. The thermal characterization showed that the Thermogravimetry (TG) analysis. The thermal performance of the
thermal conductivity of the PCM composite was more than 3.5 times NePCM composites was further examined by cooling bath and thermal
higher than that of the raw PCM, when adding 11.8 vol% of Ag nano- cycling stability analysis.
wires.
Besides metallic and metallic oxide nanoparticles, carbon-based 2. Research methodology
nanoparticles with extremely high thermal conductivity (e.g. 3000-
6600 W·m−1· K−1) [23–25] have been frequently used in many studies The overall research methodology used in this study was refined
[26–35] as additives to enhance the thermal conductivity of PCMs, based on a general framework proposed in a previous study [36], as
especially for organic PCMs which generally have lower thermal con- presented in Fig. 1. It mainly consists of three steps, including material
ductivity compared with inorganic PCMs. For instance, Ji et al. [26] selection and preparation (Step 1), material property characterization
added functionalized multi-walled carbon nanotubes (MWCNTs) into (Step 2) and thermal performance examination (Step 3).
palmitic acid to modify its thermal properties. It was found that the
well-dispersed MWCNTs in palmitic acid can effectively improve the 2.1. Materials selection and preparation of nano-enhanced PCM composites
thermal conductivity of the nano-enhanced PCM composites. Sahan
et al. [27] used activated carbon (AC) and MWCNTs as additives to 2.1.1. Materials selection
prepare AC/paraffin and MWCNTs/paraffin composites with better To facilitate solar thermal applications (e.g. solar hot water, and
thermal properties, respectively. The thermal characterization showed solar powered desiccant cooling etc.) of the nano-enhanced PCMs, a
that the thermal conductivity of the MWCNTs/paraffin and the AC/ commercial myristic acid (MA) with a nominal phase change tem-
paraffin composites increased by 38.5% and 17.2% respectively, by perature of 54–55 °C supplied from Aladdin biochemical technology Co.
adding 1 wt% and 10 wt% of the nanoparticles. The increase in the Ltd. was selected and used as based PCM. MA was selected as the base
latent heat (by 9.6%) was also observed when adding 1 wt% MWCNTs PCM mainly due to its suitable phase change temperature for low-grade
into the paraffin, compared with that of the raw paraffin. Choi et al. solar thermal storage applications. In addition, as one of the typical
[28] studied and compared the enhancement in thermal conductivity fatty acids, it also has the advantages of almost no supercooling and
when adding three different carbon nanoparticles to stearic acid. The long-term stability without phase segregation. Three different types of
comparison showed that with adding 0.1 vol% MWCNTs, graphite and carbon-based nanometre-sized materials, including GNPs, MWCNTs
graphene, the corresponding thermal conductivities of the mixtures and NG were used as high thermal conductive additives, respectively.
increased by 10.5%, 9.91% and 21.5%, respectively. Li [29] improved The properties of the base PCM and the nano-additives are summarized
the thermal conductivity of paraffin by adding nano-graphite (NG). in Table 1. It is worthwhile to mention that all chemicals used were
Random dispersion of the NG was observed in the paraffin-NG com- received from the suppliers without further functionalisation.
posites, and its thermal conductivity reached 0.9362 W·m−1· K−1 when
the NG concentration was 10 wt%, which was 7.41 times that of the raw 2.1.2. Preparation of NePCM composites
paraffin. Yuan et al. [30] used graphene nanoplatelets (GNPs) and ex- The two-step method which has been extensively used in the pre-
panded graphite (EG) as additives to modify the thermal conductivity of paration of nanofluids [37,38] was adopted to prepare the nano-en-
palmitic-stearic acid (PA-SA). It was found that although both GNPs and hanced PCM composites in this study. Firstly, MA was completely
EG can effectively improve the thermal conductivity of PA-SA linearly, melted above 70 °C for 10 min. The three types of carbon-based nano-
by adding 8 wt% of EG, the thermal conductivity was 15.8 folds higher additives were then added and mixed with the molten MA with dif-
than that of the raw PA-SA, and 2.7 folds higher than that using 8 wt% ferent weight concentrations individually following the below proce-
GNPs as additives. Seki et al. [31] modified the thermal properties of dure: magnetic stirring for 30 min first, and followed by the ultrasonic
adipic acid (AA) and sebacic acid (SA) eutectic by adding GNPs as vibration for 1 h to facilitate a homogeneous dispersion; lastly, the
additives. With different loading rates of 0.1 wt%, 0.3 wt% and 0.5 wt nano-enhanced PCM samples were cooled to the ambient temperature.
%, the thermal conductivity of AA-SA can be improved by 6%, 11% and The corresponding schematic diagram is shown in Fig. 2. The experi-
19% and the supercooling degree decreased from 17.7 °C to 6.8 °C, mental plan designed to prepare the nano-enhanced PCM composite

2
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Fig. 1. The NePCMs development framework used in this study.

samples is summarized in Table 2. (length) ×3.5 cm (width) ×0.2 cm (height). After cooled to the room
temperature and solidified, the cuboid-shaped samples were polished in
2.2. Material property characterization the upper and lower surfaces to provide two smooth surfaces for the
installation of the sensor. For the liquid phase, each sample was melted
The characterization of the PCM composites involved the mor- and poured into a 50 ml beaker and kept in a water bath at 70 °C, while
phology analysis by SEM, the chemical interaction examination by FT- the sensor was inserted into the centre of the beaker to test the thermal
IR spectroscopy, and the measurements of the thermal physical prop- conductivity under liquid phase. The measurement of each sample
erties including thermal conductivity, enthalpy-temperature relation- under solid or liquid phases was repeated for three times to obtain the
ship, and thermogravimetry feature. average values. For MA, 3 wt% GNPs/MA, 3 wt% MWCNTs/MA and
3 wt% NG/MA, their thermal conductivity was further measured under
2.2.1. SEM analysis and FT-IR spectroscopy the sample temperatures of 35 °C, 45 °C and 60 °C.
The microstructures of the carbon-based nano-enhanced materials, The enthalpy-temperature relationship of the prepared samples was
including 3 wt% GNPs/MA, 3 wt% MWCNTs/MA and 3 wt% NG/MA measured by using a DSC (DSC1 from Mettler Toledo, Switzerland) with
PCM composites were observed using a S-4800 SEM (from Hitachi, a temperature scanning rate of 10 °C min−1 under a constant nitrogen
Japan). The interactions between the carbon-based nano-additives and atmosphere of 50 ml·min−1. This DSC equipment used a so-called FRS
the based PCM were characterized by using Excalibur 3100 FT-IR 5 + DSC sensor, whose TAWN resolution and sensitivity were reported
spectrometer (from Varian, USA), with KBr pellets. The FT-IR spectrums as 0.12 and 11.9, respectively. It was calibrated using standard Indium
were recorded in the wavenumber range of 4000-400 cm−1. and Zinc samples provided by the supplier. In the DSC tests, aluminium
pans were used as the sample holders and the mass of each sample was
2.2.2. Measurement of thermal physical properties controlled within 5 mg. The resulted DSC curves within the temperature
The thermal conductivities of the prepared PCM composites, and the range of 0–100 °C were analysed in detail.
raw MA, were tested based on the transient hot-wire method by a TG curves of the prepared samples were recorded by a thermo-
thermal conductivity meter TC 3000E (from Xi’an XIATECH, China). gravimeter (TGA1 from Mettler Toledo, Switzerland) under nitrogen
The measurement accuracy of the thermal conductivity meter TC 3000E stream at a flow rate of 30 ml·min−1. The measurement accuracy of this
used was ± 3%. The thermal conductivity of these samples was mea- thermogravimeter was ± 0.5 °C. The rate of heating was 10 °C min−1 in
sured under both solid and liquid phases. For the solid phase, the the temperature range of 30–500 °C.
samples were first melted and filled into moulds with a size of 5 cm

Table 1
Main properties of MA and additives.
Substances Purity Density or tapped density Main physical parameters Manufacturers

−3
MA 98% 0.898 g·cm Melting point: 54-55 °C Aladdin biochemical technology Co. Ltd.
Latent heat of fusion: 194.90 J·g−1
−3
GNPs 99.5% 0.075 g·cm Thickness: 3-10 nm Nanjing XFNANO Materials Tech Co. Ltd
Diameter: 5-10 μm
MWCNTs ＞95% 0.270 g·cm−3 Length: 10-30 μm Nanjing XFNANO Materials Tech Co. Ltd
Diameter: ＜8 nm
NG 99.9% 0.252 g·cm−3 Thickness: < 40 nm Nanjing XFNANO Materials Tech Co. Ltd
Diameter: 3-6 μm

3
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Fig. 2. Schematic diagram of the sample preparation process.

Table 2 respectively.
Prepared NePCM composites with different weight concentrations of the ad-
ditives.
3. Results and discussion
Sample index Sample Additive concentration
3.1. Microscopic structure
0 MA –
1 GNPs/MA 1 wt%
2 GNPs/MA 2 wt% The microstructures of the PCM composites and carbon-based nano-
3 GNPs/MA 3 wt% additives are presented in Fig. 4. It can be seen from Fig. 4a that GNPs
4 MWCNTs/MA 1 wt%
characterized a glossy plane structure with the high surface-to-thick-
5 MWCNTs/MA 2 wt%
6 MWCNTs/MA 3 wt%
ness ratio. The GNPs seemed to be agglomerate, however, after mixed
7 NG/MA 1 wt% with MA through magnetic stirring and ultrasonic vibration, the GNPs
8 NG/MA 2 wt% seemed to be dispersed in the PCM, as can be seen in Fig. 4b. The plate-
9 NG/MA 3 wt% like structure of GNPs was maintained in the GNPs/MA composite,
which facilitated the formation of a functional GNPs thermal conduc-
tion network. The MA was introduced into the network and covered
2.3. Thermal performance examination
GNPs, making them look thicker (see Fig. 4b), similar to that reported
by Li [29]. Compared with that of the GNPs, the high aspect ratio of
2.3.1. Exothermic performance analysis for bulk samples
MWCNTs was achieved due to their long tubular structure, as can be
It is of great interest to investigate the thermal performance of bulk
seen in Fig. 4c. This long tubular structure with a high aspect ratio also
PCMs or PCM composites, due to the fact that it may be different from
enabled the MWCNTs to interlink with each other and resulted in a
that tested based on small samples. The exothermic performance of the
thermal conductive network after mixed with MA as shown in Fig. 4d.
bulk NePCM composites or raw MA was tested in this study using a
As MA penetrated and filled the network of MWCNTs quite well, the
water bath test apparatus, CC-K6 with Pilot ONE (from Huber,
composite looked homogeneous, indicating an acceptable dispersion of
German). It consists of a water jacketed glass vessel, a PT100 tem-
the MWCNTs, but some tips of MWCNTs still can be observed clearly.
perature sensor and a water bath thermostat, as shown in Fig. 3. The
Different from GNPs and MWCNTs, the aspect ratio of NG was relatively
melted MA or NePCM composites were filled in the water jacketed glass
lower (see Fig. 4e). As a consequence, the NG/MA composite showed a
vessel, while the glass vessel was supplied with the water flow of 30 °C
weak connection between the NG, as can be seen in Fig. 4f. The ag-
from the water bath thermostat. The accuracy of the PT100 tempera-
glomeration of NG can be clearly seen from Fig. 4e, but after prepared
ture probe used is ± 0.02 °C. It was placed in the centre of the vessel,
as NG/MA composite, a good dispersion of NG seemed to be achieved,
and the temperature measurement interval was 5 s.
as presented in Fig. 4f. It is worthwhile to note that even though ac-
ceptable dispersion of the nanoparticles in the PCM seemed to be
2.3.2. Thermal cycling test achieved through a rational composite preparation process (e.g. using
The thermal cycling stability, which is an important indicator to ultrasonic vibration), agglomerations of nanoparticles can still be found
assess the practical applicability of a PCM or PCM composite, was tested in the PCM composites, especially for GNPs/MA and MWCNTs/MA (see
through accelerated thermal cycling test in this study. In the test, the Appendix A for details).
prepared NePCM composite samples were encapsulated into 2 ml cen-
trifuge tubes, which were melted in the hot water bath at 70 °C and then
solidified at 30 °C periodically. The phase change characteristics of the 3.2. FT-IR analysis
PCM composites were measured by the DSC after per 100 thermal cy-
cles, while their crystalline and chemical structures before and after Fig. 5a presents the FT-IR spectrums of MA, GNPs and GNPs/MA
300 thermal cycles were characterized using an X-ray diffraction (XRD) composites. A series of characteristic absorption peaks can be observed
(D8 focus from Bruker, Germany), in a diffraction angle (2θ) ranging in the FT-IR spectrum of MA. For instance, the absorption peaks at
from 5°-90° and FT-IR in the wavenumber range of 4000-400 cm−1, around 2920 cm−1 and 2850 cm−1 were caused by the stretching vi-
brations of the groups −CH3 and −CH2 in MA, while the characteristic
absorption peak at around 1700 cm−1 represented the group C]O
stretching vibration. The peak of the bending vibration of −CH2, lo-
cated at around 1460 cm−1, and that of the C–H and CeC chemical
bonds laid at around 1300 cm−1. The absorption peaks of rocking vi-
bration and bending vibration can also be observed at around 933 cm−1
and 721 cm−1, respectively. Differently, no obvious absorption peak in
the FT-IR spectrum of GNPs was observed. It can also be seen from
Fig. 5a that the FT-IR spectrums of the GNPs/MA composites were
consistent with that of the raw MA, which contained the same char-
acteristic absorption peaks in MA without the emergence of any ap-
Fig. 3. Apparatus of exothermic experiments. parently new peak. It meant that no considerable chemical interaction

4
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Fig. 4. Microscopic structures for: a) GNPs; b) 3 wt% GNPs/MA; c) MWCNTs; d) 3 wt% MWCNTs/MA; and e) NG; f) 3 wt% NG/MA.

occurred between MA and GNPs, indicating that the GNPs/MA com- A similar phenomenon was reported and analysed by Harikrishnan and
posites prepared were chemically stable. Similar to that of GNPs and Kalaiselvam [19], Wu et al. [39] and Choi et al. [28]. Different from the
MA in the GNPs/MA composites, no considerable chemical interaction nonlinear increasing trend in the solid phase, the thermal conductivity
was observed between the nano-additives and MA according to the FT- of the GNPs/MA composite under the liquid phase experienced a linear
IR spectrums of MWCNTs and MWCNTs/MA, and NG and NG/MA increasing trend with the increasing of the GNPs concentration. How-
shown in Fig. 5b and c, respectively. ever, its corresponding values were much lower than that of the solid
phase. For instance, 3 wt% GNPs only contributed to the improvement
in the thermal conductivity of 59.70% for the GNPs/MA composite
3.3. Thermal conductivity analysis under the liquid phase, compared to that of 176.26% under the solid
phase.
The variation of the thermal conductivity of the nano-enhanced Similar to the GNPs/MA composites, the thermal conductivity of the
PCM composites with the change of the concentrations of different MWCNTs/MA composites increased with the increase of MWCNTs
carbon-based nano-additives is presented in Fig. 6, and the percentages concentration, following a nonlinear trend under the solid phase while
of enhancement in thermal conductivity are further summarized in a linear trend under the liquid phase, as can be seen in Fig. 6b. When
Table 3. The standard error in thermal conductivity measurements was adding 1 wt% MWCNTs into MA, the thermal conductivity of the
only up to 0.0024 W·m−1· K−1. It can be seen that the thermal con- MWCNTs/MA composite can be improved substantially by 33.71% in
ductivity of the GNPs/MA composites in both solid and liquid phases the solid phase and 3.63% in the liquid phase, respectively. By further
increased when increasing the weight concentration of GNPs (Fig. 6a). increasing the MWCNTs concentration to 3 wt%, the thermal con-
In the solid phase, the thermal conductivity of 1 wt% GNPs/MA com- ductivity reached 0.3220 W·m−1· K−1 (for solid phase) and
posite reached 0.4405 W·m−1· K−1, which was improved by 101.51% in 0.2083 W·m−1· K−1 (for liquid phase), which were 47.30% and 12.84%
comparison with that of the raw MA. By further increasing the GNPs higher than that of the raw MA, respectively.
concentration to 2 wt% and 3 wt%, the relative increments in thermal Even though similar trends were found, the thermal conductivity of
conductivity reached by 127.04% and 176.26%, respectively. However, the MWCNTs/MA composite was obviously lower than that of the
the increasing rate tended to reduce with the increase of the GNPs GNPs/MA composite. It can be explained by the heat transfer me-
concentration. It can be attributed to the agglomeration or even sedi- chanism within nano-enhanced PCMs. The enhancement in thermal
mentation phenomena which were likely to occur in the PCM compo- conductivity resulted from the high thermal conductive percolation
sites under high nanoparticle concentrations. As a consequence, the network established by the nanoparticles in the NePCM composites,
total amount of functional nanoparticles which were beneficial to the which enabled the efficient transport of phonons [40]. However, this
thermal conduction enhancement cannot be improved proportionally, enhancement was limited by the heat transfer across the interfaces
when further increasing the weight concentration of the nanoparticles.

5
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Fig. 5. FT-IR spectra for: a) MA, GNPs and GNPs/MA PCM composites; b) MA, MWCNTs and MWCNTs/MA PCM composites; and c) MA, NG and NG/MA PCM
composites.

Fig. 6. Effect of the weight concentration of nano-additives on thermal conductivity for: a) GNPs/MA PCM composites; b) MWCNTs/MA PCM composites; and c) NG/
MA PCM composites.

6
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Table 3 phase transition temperature range. When the temperature was above
The thermal conductivity of the samples. 60 °C in the liquid phase, the thermal conductivity became relatively
Samples Thermal conductivity Relative improvements [%] constant again, but it was substantially lower than that under the solid
[W·m−1·K−1] phase for each PCM composite. This phenomenon was due to the
random molecular motion in the liquid phase, which collapsed the
Solid Liquid Solid Liquid percolation network established by nanoparticles under the solid phase.
MA 0.2186 0.1846 – – This result was in agreement with the work of Qian et al. [47] and
1 wt% GNPs/MA 0.4405 0.2302 101.51 24.70 Babaei et al. [48].
2 wt% GNPs/MA 0.4963 0.2677 127.04 45.02 Base on the above discussion, it can be concluded that the carbon-
3 wt% GNPs/MA 0.6039 0.2948 176.26 59.70 based nanoparticles with inherent high thermal conductivity (reported
1 wt% MWCNTs/MA 0.2923 0.1913 33.71 3.63
to be as high as 3000-6600 W·m−1· K−1) [23–25] can effectively im-
2 wt% MWCNTs/MA 0.2988 0.2026 36.69 9.75
3 wt% MWCNTs/MA 0.3220 0.2083 47.30 12.84 prove the thermal conductivity of PCM composites. However, to max-
1 wt% NG/MA 0.2518 0.1899 15.19 2.87 imise the enhancement in the thermal conductivity of the PCM com-
2 wt% NG/MA 0.2746 0.1950 25.62 5.63 posites, the concentration rate, configuration of the nanoparticles, and
3 wt% NG/MA 0.3148 0.2059 44.01 11.54
dispersion of the nano-additives need to be carefully selected and
handled, and the attention needs to be paid on the solid/liquid state of
the PCM.
between nanoparticles within the percolation network, due to the high
thermal resistivity associated with low-frequency phonons transfer at
3.4. DSC analysis
the interface [22,41–44]. The special flake-like configuration for the
GNPs enabled the fewer interfaces in the percolation network, com-
The DSC curves of the raw MA and GNPs/MA PCM composites with
pared to the wire-like configuration for the MWCNTs, which therefore
different additive weight concentrations are presented in Fig. 8a, and
better facilitated the thermal conductivity enhancement.
the corresponding phase change temperature and latent heat are sum-
Different from the GNPs/MA and MWCNTs/MA composites, the
marized in Table 4. It can be seen that the GNPs/MA composites ex-
linear increasing trends in the thermal conductivity were found for the
perienced sharp phase change processes during both melting and
NG/MA under both solid and liquid phases. When adding 3 wt% NG
freezing processes (Fig. 8a), compared with that of the raw MA which
into MA, the thermal conductivity of the NG/MA composite was im-
had clearly wider phase change temperature ranges. The on-set tem-
proved linearly to 0.3148 W·m−1· K−1 and 0.2059 W·m−1· K−1 under
peratures of the melting process for GNPs/MA PCM composites with the
the solid and liquid phases, respectively. It was only enhanced by
GNPs concentration of 1 wt%, 2 wt%, and 3 wt% were the same of
44.01% under the solid phase, which was much lower than that of the
54.3 °C, while that of the freezing process were 54.2 °C, 54.4 °C and
GNPs/MA composite. This might be due to the existence of a threshold
54.4 °C, respectively. As a comparison, the raw MA had the on-set
for the nanoparticles, above which effective mutual interlinkage among
temperatures of 54.4 °C for the melting process, and 52.4 °C for the
the nanoparticles can be achieved to establish the percolating network
freezing process. It demonstrated the effectiveness of using GNPs to
[45,46]. This threshold was highly dependent on the aspect ratio of the
modify the phase change characteristics of MA, which approached a
nanoparticles. Since the aspect ratio of the NG particles was much
preferable isothermal phase change process. It can also be concluded
smaller, a higher weight concentration was required to establish the
that this effect was not sensitive to the concentration of the GNPs
percolation network and achieve the same level thermal conductivity
added. In the meanwhile, only a slight decrease in the latent heat was
enhancement effect as that of the GNPs.
observed for the GNPs/MA composites when compared with the raw
To gain an insight of the temperature-dependency of the thermal
MA. The latent heat for the 1 wt%, 2 wt% and 3 wt% GNPs/MA com-
conductivity enhancement, the thermal conductivity of MA, 3 wt%
posites were 189.18 J·g−1, 187.39 J·g−1 and 187.19 J·g−1 for melting,
GNPs/MA, 3 wt% MWCNTs/MA and 3 wt% NG/MA was further mea-
and 190.88 J·g−1, 189.70 J·g−1 and 188.58 J·g−1 for freezing, respec-
sured under different temperatures, as presented in Fig. 7. Similar
tively, which were only reduced by 2.93%, 3.85% and 3.96% (for
variation trends in the thermal conductivity with the variation of the
melting) and 2.19%, 2.79% and 3.37% (for freezing). It can also be seen
temperature were found for the three PCM composites. The thermal
that by adding the GNPs, the supercooling during the freezing process
conductivity was relatively constant in solid phase below 50 °C for the
of MA was effectively eliminated (Fig. 8a), due to the nucleating effect
individual PCM composites. However, a step decrease in the thermal
of the nano-additives [49] which served as nucleating agents to facil-
conductivity was observed when increasing the temperature across the
itate the crystallization of MA.
Similar to the GNPs/MA PCM composites, narrow endothermic
peaks and exothermic peaks were observed in the DSC curves for the
MWCNTs/MA composites, compared with the raw MA (Fig. 8b). The
adding of MWCNTs also showed a minor influence on the latent heat of
the PCM composites. By adding 1 wt%, 2 wt% and 3 wt% MWCNTs into
MA, the latent heat of fusion decreased from 194.90 J·g−1 to
191.99 J·g−1, 190.12 J·g−1 and 188.47 J·g−1, while the latent heat of
solidification dropped from 195.15 J·g−1 to 192.19 J·g−1, 190.99 J·g−1
and 188.21 J·g−1, respectively. The phenomenon with smaller phase
change temperature ranges can also be found when adding NG into MA,
together with the minor influence on the latent heat of the NG/MA
composites, as presented in Fig. 8c and Table 4. It seemed that the three
carbon-based nano-additives had a similar effect in terms of modifying
the phase change process of MA.

3.5. TGA analysis

Fig. 7. Thermal conductivity of MA, 3 wt% GNPs/MA, 3 wt% MWCNTs/MA, Fig. 9 presents the TGA profiles for PCM composites modified by
and 3 wt% NG/MA as a function of temperature. different carbon-based nano-additives, and the details are further

7
M. He, et al. Journal of Energy Storage 25 (2019) 100874

a) 20 b) 25
15 MA 20 MA
10 1 wt% GNPs/MA 1 wt% MWCNTs/MA
15 2 wt% MWCNTs/MA
2 wt% GNPs/MA

Heat flow (mW·mg -1)

Heat flow (mW·mg -1) 5 3 wt% GNPs/MA 10 3 wt% MWCNTs/MA

0 5
-5 0

-10 -5
-10
-15
-15
-20
-20
-10 0 10 20 30 40 50 60 70 80 90 100 110 -10 0 10 20 30 40 50 60 70 80 90 100 110
Temperature ( ) Temperature ( )

20
c)
15 MA
1 wt% NG/MA
10
2 wt% NG/MA
Heat flow (mW·mg -1)

5 3 wt% NG/MA

0
-5
-10
-15
-20
-10 0 10 20 30 40 50 60 70 80 90 100 110
Temperature ( )

Fig. 8. DSC melting and freezing curves for: a) GNPs/MA PCM composites; b) MWCNTs/MA PCM composites; and c) NG/MA PCM composites.

Table 4 composites. It can be seen that with the modification of 3 wt% GNPs,
Phase change temperature and phase change latent heat of GNPs/MA, the discharging time of GNPs/MA composite was significantly reduced
MWCNTs/MA and NG/MA PCM composites. to around 4900 s, compared to that of around 7450 s for the raw MA.
Samples On-set temperature (oC) Latent heat (J·g−1) Similarly, the adding of MWCNTs with a concentration of 3 wt% sub-
stantially decreased the discharging time of MWCNTs/MA composite to
Melting Freezing Melting Freezing around 5850 s. When using NG particles as additives with a con-
centration of 3 wt%, the discharging time of the PCM composite can
MA 54.4 52.4 194.90 195.15
1 wt% GNPs/MA 54.3 54.2 189.18 190.88
also be reduced to around 6500 s. The performance of the bulk PCM
2 wt% GNPs/MA 54.3 54.4 187.39 189.70 composites reflected the thermal conductivity enhancement effect due
3 wt% GNPs/MA 54.3 54.4 187.19 188.58 to the adding of different carbon-based nano-additives. The GNPs/MA
1 wt% MWCNTs/MA 54.6 53.2 191.99 192.19 composites whose thermal conductivity enhancement outperformed the
2 wt% MWCNTs/MA 54.4 53.1 190.12 190.99
other two PCM composites provided the best bulk PCM discharging
3 wt% MWCNTs/MA 54.4 53.0 188.47 188.21
1 wt% NG/MA 54.5 54.1 191.27 191.15 performance with the shortest discharging completion time.
2 wt% NG/MA 54.6 54.0 190.67 191.04
3 wt% NG/MA 54.6 54.1 188.90 188.85
3.7. Thermal cycling stability analysis

summarized in Table 5. By adding the three types of nano-additives into Fig. 11 compares the DSC curves of different PCM composite sam-
MA, lower initial decomposition temperatures were found, ranging ples after per 100 thermal cycles, and Table 6 further summarizes the
from 116.0 °C to 142.8 °C with different carbon-based nanoparticles corresponding phase change characteristics before and after the
concentrations, compared with that of 152.3 °C for the raw MA. The thermal cycles. It can be seen from Fig. 11a that the DSC curves of the
decrease in the initial decomposition temperature was acceptable, as it 3 wt% GNPs/MA PCM composite after 100, 200 and 300 thermal cycles
was far higher than the phase change temperature range of MA, even were consistent with the original DSC curve before thermal cycling, and
though it reduced the safe operation temperature limits for the appli- no obvious degradation can be observed. The corresponding latent heat
cation of the NePCM composites. Similar to the initial decomposition of fusion was 185.69 J·g−1 after 100 thermal cycles, 183.43 J·g−1 after
temperature, the lower maximal and final decomposition temperatures 200 thermal cycles, and 181.75 J·g−1 after 300 thermal cycles, re-
were observed for the PCM composites modified with the carbon-based spectively, which were only reduced by 0.80%, 2.01%, and 2.91% in
nanoparticles additives. comparison with the sample before the thermal cycling test. Slight
decrease was also observed for the latent heat of solidification with the
increasing number of thermal cycles, but it was still above 180 J·g−1
3.6. Exothermic process analysis after 300 thermal cycles (i.e. 180.44 J·g−1) with a degradation no more
than 5% (i.e. 4.32%). Similarly, negligible degradation in the latent
Fig. 10 presents the thermal exothermic process of the bulk MA, heat was found in 3 wt% MWCNTs/MA composite sample (see Fig. 11b)
3 wt% GNPs/MA, 3 wt% MWCNTs/MA and 3 wt% NG/MA PCM and 3 wt% NG/MA composite sample (see Fig. 11c) after the thermal

8
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Fig. 9. TGA profiles for: a) MA and GNPs/MA PCM composites; b) MA and MWCNTs/MA PCM composites; and c) MA and NG/MA PCM composites.

Table 5
TGA data for MA, GNPs/MA, MWCNTs/MA and NG/MA PCM composites.
Samples Initial decomposition temperature (oC) Maximal decomposition temperature (oC) Final decomposition temperature (oC)

MA 152.3 238.2 241.5

1 wt% GNPs/MA 130.2 214.5 219.5
2 wt% GNPs/MA 130.2 215.5 220.2
3 wt% GNPs/MA 116.0 210.8 215.7
1 wt% MWCNTs/MA 142.8 221.0 224.5
2 wt% MWCNTs/MA 133.3 219.0 223.3
3 wt% MWCNTs/MA 118.0 213.5 218.3
1 wt% NG/MA 130.0 214.5 220.2
2 wt% NG/MA 132.2 216.8 221.0
3 wt% NG/MA 132.3 217.7 223.3

cycling tests, even though a slight decreasing trend was found when the
number of thermal cycles increased. The above results indicated that
the nano-enhanced PCM composites prepared had excellent thermal
cycling stabilities.
To further obtain a clear insight into the thermal cycling stability,
the crystallization characteristics of the PCM composite samples were
studied before and after 300 thermal cycling tests, as presented in
Fig. 12. It can be seen that the two XRD curves were consistent with
each other for the GNPs/MA composite (Fig. 12a), with identical ab-
sorption peaks at the diffraction angles (2θ) of 5.6°, 8.3°, 14.0°, 21.5°
and 24.1°. It indicated that after 300 thermal cycles, the crystalline
structure of 3 wt% GNPs/MA composite sample remained the same.
Similar results can be concluded by comparing the XRD curves of 3 wt%
MWCNTs/MA (see Fig. 12b) and 3 wt% NG/MA (see Fig. 12c) before
and after the thermal cycling test. By further analysing the FT-IR
spectra before and after 300 thermal cycles of the PCM composites
Fig. 10. Cooling curves of MA, 3 wt% GNPs/MA, 3 wt% MWCNTs/MA and 3 wt
(Fig. 13), it was found that the chemical structures of the PCM com-
% NG/MA PCM composites.
posites also remained stable.

9
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Fig. 11. DSC melting and freezing curves after 0, 100, 200 and 300 thermal cycles for: a) 3 wt% GNPs/MA PCM composite; b) 3 wt% MWCNTs/MA PCM composite;
and c) 3 wt% NG/MA PCM composite.

Table 6 to characterize the NePCMs. The thermal cycling stability and bulk
Phase change characteristics of thermal cycles of 3 wt% GNPs/MA, 3 wt% thermal performance of the NePCMs prepared were also tested.
MWCNTs/MA and 3 wt% NG/MA PCM composites. It was found that the thermal conductivity of the MA can be sig-
Samples Number of On-set temperature (oC) Latent heat (J·g−1) nificantly enhanced by adding carbon-based nanoparticles as additives.
thermal Using GNPs as nano-additives showed the best performance for thermal
cycles conductivity enhancement, which can significantly enhance the
Melting Freezing Melting Freezing
thermal conductivity by 176.26% under the solid phase, compared to
3 wt% GNPs/ 0 54.3 54.4 187.19 188.58 that of 47.30% and 44.01% for MWCNTs and NG, under the weight
MA 100 54.3 54.3 185.69 183.00 concentration of 3 wt%. The high thermal conductivity in the GNPs/MA
200 54.5 54.1 183.43 181.83 composite was due to the thermal conductive percolation network es-
300 54.4 54.2 181.75 180.44 tablished by GNPs, which tended to have fewer nanoparticle interfaces,
3 wt% 0 54.4 53.0 188.47 188.21
MWCNTs/ 100 54.5 53.6 185.77 186.32
compared to the other carbon-based nanoparticles. Different increasing
MA 200 54.3 53.5 184.85 184.97 trends of the thermal conductivity variation under solid and liquid
300 54.3 53.8 181.21 182.86 phases were found for both GNPs/MA and MWCNTs/MA PCM com-
3 wt% NG/MA 0 54.6 54.1 188.90 188.85 posites, which were nonlinear and linear, respectively. However, due to
100 54.5 54.1 186.02 187.13
the existence of random molecular motion under the liquid phase, the
200 54.6 54.0 184.38 185.93
300 54.4 54.1 182.68 182.03 corresponding thermal conductivity enhancement was evidently lower
than that under the solid phase. To establish and maintain an effective
percolation network, the concentration of nano-additives needs to be
4. Conclusions carefully selected. It was also found that by adding the carbon-based
nanoparticles as additives into the MA, the phase transition can be
This paper presented the preparation, thermal characterization and modified to approach an isothermal process and the supercooling effect
thermal examination of phase change materials (PCMs) enhanced by can be eliminated, while still maintaining a high latent heat capacity.
carbon-based nano-particles for solar thermal energy storage applica- The thermal cycling and bulk thermal performance tests further de-
tions, by following a proposed experimental framework for the devel- monstrated the practical performance of the NePCMs developed, which
opment of nano-enhanced PCMs (NePCMs). The framework enabled the exhibited excellent thermal cycling stability, and great potential to re-
comprehensive assessment of the prepared NePCMs before practical duce the phase change time required. The NePCMs developed can be
applications. Myristic acid (MA) was selected as the base PCM, and used to develop solar thermal energy storage for various applications.
three different types of nanoparticles, i.e. graphene nanoplatelets
(GNPs), multi-walled carbon nanotubes (MWCNTs) and nano-graphite Acknowledgement
(NG), were used as nano-additives to modify the properties of the base
PCM. The properties of the nanoparticles/PCM composites, including This paper was financially supported by the National Key Research
morphology, chemical stability, and thermal properties were measured and Development Program of China (No. 2016YFD0400905).

10
M. He, et al. Journal of Energy Storage 25 (2019) 100874

Fig. 12. XRD spectra before and after cycling for: a) 3 wt% GNPs/MA PCM composite; b) 3 wt% MWCNTs/MA PCM composite; and c) 3 wt% NG/MA PCM composite.

Fig. 13. FT-IR spectra before and after cycling for: a) 3 wt% GNPs/MA PCM composite; b) 3 wt% MWCNTs/MA PCM composite; and c) 3 wt% NG/MA PCM
composite.

11
M. He, et al. Journal of Energy Storage 25 (2019) 100874

References [25] A.A. Balandin, S. Ghosh, W.Z. Bao, I. Calizo, D. Teweldebrhan, F. Miao, C.N. Lau,
Superior thermal conductivity of single-layer graphene, Nano Lett. 8 (2008)
902–907.
[1] W.Y. Lin, Z.J. Ma, P. Cooper, M.I. Sohel, L.W. Yang, Thermal performance in- [26] P.J. Ji, H.H. Sun, Y.X. Zhong, W. Feng, Improvement of the thermal conductivity of
vestigation and optimization of buildings with integrated phase change materials a phase change material by the functionalized carbon nanotubes, Chem. Eng. Sci. 81
and solar photovoltaic thermal collectors, Energy Build. 116 (2016) 562–573. (2012) 140–145.
[2] J. Giro-Paloma, Y. Konuklu, A.I. Fernandez, Preparation and exhaustive char- [27] N. Sahan, M. Fois, H. Paksoy, The effects of various carbon derivative additives on
acterization of paraffin or palmitic acid microcapsules as novel phase change ma- the thermal properties of paraffin as a phase change material, Int. J. Energy Res. 40
terial, Sol. Energy 112 (2015) 300–309. (2016) 198–206.
[3] B. Xu, P.W. Li, C. Chan, Application of phase change materials for thermal energy [28] D.H. Choi, J. Lee, H. Hong, Y.T. Kang, Thermal conductivity and heat transfer
storage in concentrated solar thermal power plants: a review to recent develop- performance enhancement of phase change materials (PCM) containing carbon
ments, Appl. Energy 160 (2015) 286–307. additives for heat storage application, Int. J. Refrig. 42 (2014) 112–120.
[4] P.B. Salunkhe, P.S. Shembekar, A review on effect of phase change material en- [29] M. Li, A nano-graphite/paraffin phase change material with high thermal con-
capsulation on the thermal performance of a system, Renew. Sust. Energ. Rev. 16 ductivity, Appl. Energy 106 (2013) 25–30.
(2012) 5603–5616. [30] Y.P. Yuan, N. Zhang, T.Y. Li, X.L. Cao, W.Y. Long, Thermal performance enhance-
[5] X.W. Fu, Z.M. Liu, B. Wu, J.L. Wang, J.X. Lei, Preparation and thermal properties of ment of palmitic-stearic acid by adding graphene nanoplatelets and expanded
stearic acid/diatomite composites as form-stable phase change materials for graphite for thermal energy storage: a comparative study, Energy 97 (2016)
thermal energy storage via direct impregnation method, J. Therm. Anal. Calorim. 488–497.
123 (2016) 1173–1181. [31] Y. Seki, S. Ince, M.A. Ezan, A. Turgut, A. Erek, Graphite nanoplates loading into
[6] N.L. Panwar, S.C. Kaushik, S. Kothari, Role of renewable energy sources in en- eutectic mixture of Adipic acid and Sebacic acid as phase change material, Sol.
vironmental protection: a review, Renew. Sust. Energ. Rev. 15 (2011) 1513–1524. Energy Mater. Sol. Cells 140 (2015) 457–463.
[7] M.M.A. Khan, N.I. Ibrahim, I.M. Mahbubul, H.M. Ali, R. Saidur, F.A. Al-Sulaiman, [32] L. Colla, L. Fedele, S. Mancin, L. Danza, O. Manca, Nano-PCMs for enhanced energy
Evaluation of solar collector designs with integrated latent heat thermal energy storage and passive cooling applications, Appl. Therm. Eng. 110 (2017) 584–589.
storage: a review, Sol. Energy 166 (2018) 334–350. [33] J.F. Wang, H.Q. Xie, Z. Xin, Y. Li, L.F. Chen, Enhancing thermal conductivity of
[8] J.M. Khodadadi, L.W. Fan, H. Babaei, Thermal conductivity enhancement of na- palmitic acid based phase change materials with carbon nanotubes as fillers, Sol.
nostructure-based colloidal suspensions utilized as phase change materials for Energy 84 (2010) 339–344.
thermal energy storage: a review, Renew. Sust. Energ. Rev. 24 (2013) 418–444. [34] A. Allahbakhsh, M. Arjmand, Graphene-based phase change composites for energy
[9] Z.A. Qureshi, H.M. Ali, S. Khushnood, Recent advances on thermal conductivity harvesting and storage: state of the art and future prospects, Carbon 148 (2019)
enhancement of phase change materials for energy storage system: a review, Int. J. 441–480.
Heat Mass Transf. 127 (2018) 838–856. [35] W.T. Wang, B.T. Tang, B.Z. Ju, Z.M. Gao, J.H. Xiu, S.F. Zhang, Fe3O4-functionalized
[10] C.Y. Zhao, W. Lu, Y. Tian, Heat transfer enhancement for thermal energy storage graphene nanosheet embedded phase change material composites: efficient mag-
using metal foams embedded within phase change materials (PCMs), Sol. Energy 84 netic- and sunlight-driven energy conversion and storage, J. Mater. Chem. A 5
(2010) 1402–1412. (2017) 958–968.
[11] S. Motahar, N. Nikkam, A.A. Alemrajabi, R. Khodabandeh, M.S. Toprak, [36] Z.J. Ma, W.Y. Lin, M.I. Sohel, Nano-enhanced phase change materials for improved
M. Muhammed, A novel phase change material containing mesoporous silica na- building performance, Renew. Sust. Energ. Rev. 58 (2016) 1256–1268.
noparticles for thermal storage: A study on thermal conductivity and viscosity, Int. [37] X.Q. Wang, A.S. Mujumdar, Heat transfer characteristics of nanofluids: a review,
Commun. Heat Mass Transf. 56 (2014) 114–120. Int. J. Therm. Sci. 46 (2007) 1–19.
[12] T. Khadiran, M.Z. Hussein, Z. Zainal, R. Rusli, Encapsulation techniques for organic [38] V. Kumaresan, P. Chandrasekaran, M. Nanda, A.K. Maini, R. Velraj, Role of PCM
phase change materials as thermal energy storage medium: a review, Sol. Energy based nanofluids for energy efficient cool thermal storage system, Int. J. Refrig. 36
Mater. Sol. Cells 143 (2015) 78–98. (2013) 1641–1647.
[13] H.M. Tauseef-ur-Rehman, M.M. Ali, U. Janjua, W.M. Sajjad, Yan, A critical review [39] S.Y. Wu, H. Wang, S. Xiao, D.S. Zhu, An investigation of melting/freezing char-
on heat transfer augmentation of phase change materials embedded with porous acteristics of nanoparticle-enhanced phase change materials, J. Therm. Anal.
materials/foams, Int. J. Heat Mass Transf. 135 (2019) 649–673. Calorim. 110 (2012) 1127–1131.
[14] M.M. Umair, Y. Zhang, K. Iqbal, S.F. Zhang, B.T. Tang, Novel strategies and sup- [40] F. Yavari, H.R. Fard, K. Pashayi, M.A. Rafiee, A. Zamiri, Z.Z. Yu, R. Ozisik, T. Borca-
porting materials applied to shape-stabilize organic phase change materials for Tasciuc, N. Koratkar, Enhanced thermal conductivity in a nanostructured phase
thermal energy storage–a review, Appl. Energy 235 (2019) 846–873. change composite due to low concentration graphene additives, J. Phys. Chem. C
[15] Z.W. Ge, Y.L. Li, D.C. Li, Z. Sun, Y. Jin, C.P. Liu, C. Li, G.H. Leng, Y.L. Long, Thermal 115 (2011) 8753–8758.
energy storage: challenges and the role of particle technology, Particuology 15 [41] S.T. Huxtable, D.G. Cahill, S. Shenogin, L.P. Xue, R. Ozisik, P. Barone, M. Usrey,
(2014) 2–8. M.S. Strano, G. Siddons, M. Shim, P. Keblinski, Interfacial heat flow in carbon na-
[16] M.A. Kibria, M.R. Anisur, M.H. Mahfuz, R. Saidur, I.H.S.C. Metselaar, A review on notube suspensions, Nat. Mater. 2 (2003) 731–734.
thermophysical properties of nanoparticle dispersed phase change materials, Energ. [42] J.N. Shi, M.D. Ger, Y.M. Liu, Y.C. Fan, N.T. Wen, C.K. Lin, N.W. Pu, Improving the
Convers. Manage. 95 (2015) 69–89. thermal conductivity and shape-stabilization of phase change materials using na-
[17] G.R. Dheep, A. Sreekumar, Influence of nanomaterials on properties of latent heat nographite additives, Carbon 51 (2013) 365–372.
solar thermal energy storage materials – a review, Energ. Convers. Manage. 83 [43] J.L. Zeng, Y.Y. Liu, Z.X. Cao, J. Zhang, Z.H. Zhang, L.X. Sun, F. Xu, Thermal con-
(2014) 133–148. ductivity enhancement of MWNTs on the PANI/Tetradecanol form-stable PCM, J.
[18] H. Nazir, M. Batool, F.J.B. Osorio, M. Isaza-Ruiz, X.H. Xu, K. Vignarooban, Therm. Anal. Calorim. 91 (2008) 443–446.
P. Phelan, A.M. Inamuddin, Kannan, Recent developments in phase change mate- [44] J.L. Zeng, Z. Cao, D.W. Yang, F. Xu, L.X. Sun, X.F. Zhang, L. Zhang, Effects of
rials for energy storage applications: a review, Int. J. Heat Mass Transf. 129 (2019) MWNTs on phase change enthalpy and thermal conductivity of a solid-liquid or-
491–523. ganic PCM, J. Therm. Anal. Calorim. 95 (2009) 507–512.
[19] S. Harikrishnan, S. Kalaiselvam, Preparation and thermal characteristics of CuO- [45] M. Shtein, R. Nadiv, M. Buzaglo, K. Kahil, O. Regev, Thermally conductive gra-
oleic acid nanofluids as a phase change material, Thermochim. Acta 533 (2012) phene-polymer composites: size, percolation, and synergy effects, Chem. Mater. 27
46–55. (2015) 2100–2106.
[20] J.L. Zeng, F.R. Zhu, S.B. Yu, L. Zhu, Z. Cao, L.X. Sun, G.R. Deng, W.P. Yan, L. Zhang, [46] Y.S. Liu, Y.Z. Yang, Use of nano-α-Al2O3 to improve binary eutectic hydrated salt as
Effects of copper nanowires on the properties of an organic phase change material, phase change material, Sol. Energy Mater. Sol. Cells 160 (2017) 18–25.
Sol. Energy Mater. Sol. Cells 105 (2012) 174–178. [47] T.T. Qian, J.H. Li, W.W. Feng, H.E. Nian, Single-walled carbon nanotube for shape
[21] N. Sahan, H.O. Paksoy, Thermal enhancement of paraffin as a phase change ma- stabilization and enhanced phase change heat transfer of polyethylene glycol phase
terial with nanomagnetite, Sol. Energy Mater. Sol. Cells 126 (2014) 56–61. change material, Energ. Convers. Manage. 143 (2017) 96–108.
[22] J.L. Zeng, Z. Cao, D.W. Yang, L.X. Sun, L. Zhang, Thermal conductivity enhance- [48] H. Babaei, P. Keblinski, J.M. Khodadadi, Improvement in thermal conductivity of
ment of Ag nanowires on an organic phase change material, J. Therm. Anal. paraffin by adding high aspect-ratio carbon-based nano-fillers, Phys. Lett. A 377
Calorim. 101 (2010) 385–389. (2013) 1358–1361.
[23] P. Kim, L. Shi, A. Majumdar, P.L. McEuen, Thermal transport measurements of [49] T.D. Dao, H.M. Jeong, Novel stearic acid/graphene core-shell composite micro-
individual multiwalled nanotubes, Phys. Rev. Lett. 87 (2001) 215502. capsule as a phase change material exhibiting high shape stability and performance,
[24] S. Berber, Y.K. Kwon, D. Tomanek, Unusually high thermal conductivity of carbon Sol. Energy Mater. Sol. Cells 137 (2015) 227–234.
nanotubes, Phys. Rev. Lett. 84 (2000) 4613–4616.

12