A Unique Strategy for Polyethylene Glycol/Hybrid Carbon Foam Phase Change Materials: Morphologies, Thermal Properties, and Energy Storage Behavior

Polyethylene glycol (PEG)/hybrid carbon foam (CF) phase change materials (PCMs) were prepared by integrating PEG into CF via dynamic-vacuum impregnation. The hybrid CF was first synthesized by mixtures of graphene oxide (GO) and carbon nanotubes (CNTs) with different volume ratios. The morphologies, chemical structures, thermal conductivities, shape-stabilization levels, and photo-thermal energy conversion levels of these composite PCMs were characterized systematically. The prepared composite PCMs exhibited good shape-stabilization levels and showed their original shapes without any PEG leakage. It was found that the polyethylene glycol/carbon foam with multi-walled carbon nanotubes (PEG/MCF) composite PCMs had a better shape-stable performance below the temperature of 250 °C, and the thermal conductivity of the PEG/MCF composite PCMs reached as high as 1.535 W/(mK), which was obviously higher than that of polyethylene glycol/carbon foam with single-walled carbon nanotubes (PEG/SCF, 1.159 W/(mK)). The results of the photo-thermal simulation tests showed that the composite PCMs had the ability to absorb light energy and then convert it to thermal energy, and the maximum thermal energy storage efficiency of the PEG/MCF composite PCMs and the PEG/SCF composite PCMs was 92.1% and 90.6%, respectively. It was considered that a valuable technique to produce high-performance composite PCMs was developed.


Introduction
Nowadays, the grave environmental crisis and the shortage of fossil fuels drive the development and efficient utilization of various kinds of renewable energy. Thermal energy storage, a popular method of energy storage, plays an important role in improving energy utilization by coordinating the energy supply and demand. Three kinds of thermal energy storage (TES), including sensible heat storage, latent heat storage, and thermal-chemical heat storage, are widely used at present [1][2][3]. Among these, latent thermal energy storage using phase change materials (PCMs) is one of the most widely applied in various fields, including building energy, air-conditioning, solar thermal storage, smart textiles, heat pumps, waste heat recovery, and electronic devices [4][5][6][7][8][9]. All of these applications rely on the distinguished behaviors of PCMs, including their high latent heat storage, usability at larger scales, low cost, and quasi-constant phase change temperatures.
In accordance with the phase change process, PCMs are divided into four categories: solid-solid PCMs, solid-liquid PCMs, solid-gas PCMs, and liquid-gas PCMs [10]. At present, intensive research has been conducted on solid-liquid PCMs, where polyethylene glycol (PEG) is one of the most

Fabrication of Composite PCMs
The composite PCMs were prepared by the dynamic-vacuum impregnation method. This method utilizes a continuous and synchronous ultrasonic vibration and vacuum field, which can accelerate the dip rate of the PEG and enhance the impregnation effect. PEG was first melted at 90 • C in a beaker. After melting completely, the CFs with different volume ratios were added to the PEG melt with ultrasonic vibration and absolutely impregnated with the PEG melt under a vacuum condition with a pressure of 0.05 MPa. Finally, the samples were cooled at room temperature. The obtained composite PCMs were named PEG/MCF 0.3-1.1 (or PEG/SCF 0.3-1.1 ), and the mechanism is shown in Figure 1.

Fabrication of Composite PCMs
The composite PCMs were prepared by the dynamic-vacuum impregnation method. This method utilizes a continuous and synchronous ultrasonic vibration and vacuum field, which can accelerate the dip rate of the PEG and enhance the impregnation effect. PEG was first melted at 90 °C in a beaker. After melting completely, the CFs with different volume ratios were added to the PEG melt with ultrasonic vibration and absolutely impregnated with the PEG melt under a vacuum condition with a pressure of 0.05 MPa. Finally, the samples were cooled at room temperature. The obtained composite PCMs were named PEG/MCF0.3-1.1 (or PEG/SCF0.3-1.1), and the mechanism is shown in Figure 1.

Experimental Method
The morphologies and microstructures were characterized via a field-emission transmission electron microscopy (FE-TEM) instrument (JEM-F200, Japan electronic plant type association, Tokyo, Japan) and a field-emission scanning electron microscopy (FE-SEM) instrument (Gemini 500, Carl Zeiss AG, Brunswick, Germany). The fourier transform infrared (FT-IR) spectroscopy results for the composite PCMs were recorded on a Nicolet iS50 FT-IR spectrometer (Nicolet, Madison, WI, USA) in the wavenumber range of 400-4000 cm −1 . The X-ray diffraction (XRD) patterns of the pure PEG and composite PCMs were obtained using a D8 diffractometer (Bruker D8 Advance, Dresdon, Germany) in the range of diffraction angle 2θ = 5-50 • , with a scan speed of 3 • /min at room temperature. The Brunauer-Emmett-Teller (BET) surface area was measured using the nitrogen adsorption method at the temperature of liquid nitrogen on a Micromeritics ASAP 2020 Plus HD88 system (Micromeritics Instrument Corp, Atlantic, GA, USA). A thermogravimetry analysis (TG) was performed using a TGA-2 (Berne, Switzerland) thermal analyzer with a heating rate of 10 • C/min from 30 to 600 • C in a nitrogen atmosphere. The samples were maintained at 10 mg. The degree of crystallization was investigated using a polarizing microscope (Axioskop40, Olympus, Tokyo, Japan) with a magnification of 200× or 500×, and the samples were placed on a hot plate with a heating rate/cooling rate of 5 • C/min from 30 to 70 • C. The thermal energy storage of the pure PEG and the composite PCMs was tested using a differential scanning calorimeter (DSC, Q100, Berne, Switzerland) with a heating/cooling rate of 10 • C/min in a purified nitrogen. The samples were maintained at 10 mg to ensure instrument sensitivity. The thermal energy conversion was determined by a xenon lamp (CEL-HXUV300, Chengdu, China) and an optical power meter (CEL-NP2000, Chengdu, China). To characterize the shape-stabilization, samples of 3 × 1 cm (radius × height) in size were placed on a hot plate and digital photos of them were obtained at different temperatures. The thermal conductivity was detected with a thermal conductivity analyzer (DRH-300, Nanjing, China) using the double heat plate method. A sample of 3 × 1 cm (radius × height) in size was placed on the double hot plate, the temperature was adjusted to 30 • C, and the thermal conductivity value was recorded. Then, the temperature was adjusted to 40, 50, 60, and 70 • C, respectively, and the thermal conductivity value was recorded. The data was obtained from the average value of three measurements for each sample. Meanwhile, the shape of the samples after the measurement under different temperature was observed. The solar-to-thermal energy conversion experiment was carried out as follows. First, the intensity of a xenon lamp was adjusted to 100 mW/cm 2 as the simulated solar light source, and then samples of 2 × 1 × 0.1 cm in size were maintained for 600 s with light, after which the light was immediately removed. During this process, the temperature change of each sample was recorded.

Physical Properties and Micro-Morphology of Hybrid Carbon Foam
It is clear that the physical properties of a hybrid carbon foam, such as its pore size and adsorption quality, directly determine the PEG injection behavior in CF and then the thermal energy storage of composite PCMs. The N 2 adsorption analysis results for the hybrid carbon foam with different MWCNTs (or SWCNTs)/GO volume ratios are shown in Figure 2. All of adsorption hysteresis loops are between 0.8 and 1.0 P/P 0 , which can be classified as type I CF [22]. As the MWCNTs (or SWCNTs) content increases, the MCF absorption capacity decreases from 77.14 to 23.58 cm 3 /g, while that for SCF is found to decrease from 58.09 to 39.48 cm 3 /g, indicating that the introduction of CNTs affects the adsorption quantity of the carbon foam. The embedded figure in Figure 2 shows that the average nano-pores size of the carbon foam "cell walls" decreases from 12.79 to 5.71 nm or 11.28 to 8.48 nm, respectively, with an increase in the nanotube content. Meanwhile, the average micro-pores size of the carbon foams with different volume ratios is 19.43~27.19 µm. This indicates that nanotubes can contribute to the formation of the carbon foam framework. The abundant pores in carbon foam play the role of PEG adsorption. Materials 2018, 11, x FOR PEER REVIEW 5 of 16 The morphologies and microstructures of the hybrid carbon foams with the different MWCNTs (or SWCNTs)/GO volume ratios were observed using the FE-TEM equipment. As illustrated in Figure  3a,b,e,f, it was found that the MWCNTs (or SWCNTs), with lengths ranging from several hundred nanometers to tens of micrometers, were dispersed on or between the GO nanosheets. However, it can be obviously found that excessive introduction of CNTs resulted in its accumulation on GO nanosheets, as shown in Figure 3e,f. This could be harmful to the PEG dipping process.

Chemical Structures of PEG and Microstructures of Composite PCMs
FT-IR was used to detect the chemical structures of PEG both pure and in the composite PCMs. As shown in Figure 4, the characteristic peak of pure PEG is characterized by the stretching vibration of O-H at approximately 3450 cm −1 , the C-H stretching vibration at approximately 2900 cm −1 , the C=O stretching vibration at 1690 cm −1 , and the C-O stretching vibration at 1108 cm −1 [23]. Meanwhile, the absorption peaks of the main functional groups of PEG in the composite PCMs are almost the same as those of pure PEG, with only slight displacements. Furthermore, it was found that there were no obvious peak displacements of functional groups of PEG in the CF with a different composition and microstructure. The morphologies and microstructures of the hybrid carbon foams with the different MWCNTs (or SWCNTs)/GO volume ratios were observed using the FE-TEM equipment. As illustrated in Figure 3a,b,e,f, it was found that the MWCNTs (or SWCNTs), with lengths ranging from several hundred nanometers to tens of micrometers, were dispersed on or between the GO nanosheets. However, it can be obviously found that excessive introduction of CNTs resulted in its accumulation on GO nanosheets, as shown in Figure 3e,f. This could be harmful to the PEG dipping process. The morphologies and microstructures of the hybrid carbon foams with the different MWCNTs (or SWCNTs)/GO volume ratios were observed using the FE-TEM equipment. As illustrated in Figure  3a,b,e,f, it was found that the MWCNTs (or SWCNTs), with lengths ranging from several hundred nanometers to tens of micrometers, were dispersed on or between the GO nanosheets. However, it can be obviously found that excessive introduction of CNTs resulted in its accumulation on GO nanosheets, as shown in Figure 3e,f. This could be harmful to the PEG dipping process.

Chemical Structures of PEG and Microstructures of Composite PCMs
FT-IR was used to detect the chemical structures of PEG both pure and in the composite PCMs. As shown in Figure 4, the characteristic peak of pure PEG is characterized by the stretching vibration of O-H at approximately 3450 cm −1 , the C-H stretching vibration at approximately 2900 cm −1 , the C=O stretching vibration at 1690 cm −1 , and the C-O stretching vibration at 1108 cm −1 [23]. Meanwhile, the absorption peaks of the main functional groups of PEG in the composite PCMs are almost the same as those of pure PEG, with only slight displacements. Furthermore, it was found that there were no obvious peak displacements of functional groups of PEG in the CF with a different composition and microstructure.

Chemical Structures of PEG and Microstructures of Composite PCMs
FT-IR was used to detect the chemical structures of PEG both pure and in the composite PCMs. As shown in Figure 4, the characteristic peak of pure PEG is characterized by the stretching vibration of O-H at approximately 3450 cm −1 , the C-H stretching vibration at approximately 2900 cm −1 , the C=O stretching vibration at 1690 cm −1 , and the C-O stretching vibration at 1108 cm −1 [23]. Meanwhile, the absorption peaks of the main functional groups of PEG in the composite PCMs are almost the same as those of pure PEG, with only slight displacements. Furthermore, it was found that there were no obvious peak displacements of functional groups of PEG in the CF with a different composition and microstructure.  The microstructures of the hybrid carbon foams with different MWCNTs (or SWCNTs)/GO volume ratios and their composite PCMs were observed by FE-SEM.
As shown in Figure 5a,b,d,e, a porous structure of the carbon foams is found. The GO presents a sheet-like wrinkled surface, which can act as a platform to enhance the ligation of sheets with the CNTs and improve the structural stability of the hybrid carbon foam. The MWCNTs (or SWCNTs) were well-dispersed on or between the GO sheets and effectively combined with the GO sheets. The MWCNTs (or SWCNTs) played a bridge role for the entire hybrid carbon foam, effectively connecting the sheet-like structure of GO, which had a positive impact on the preparation of the composite PCMs. Therefore, the unique structure of the hybrid carbon foam was constructed by the interaction of the sheet-like GO and MWCNTs (or SWCNTs), which will improve the shape stability of the composite PCMs. As shown in Figure 5c,f, the hybrid carbon foam was filled with PEG after the impregnation. There was no terrible surface fracture or cracking in the cross-sections of the composite PCMs. This indicated that the PEG and hybrid carbon foam had good compatibility.

Macroscopic Phase Change of PEG and Composite PCMs
The macroscopic phase changes of the composite PCMs are very important for preventing the leakage of materials. As shown in Figure 5a,b,d,e, a porous structure of the carbon foams is found. The GO presents a sheet-like wrinkled surface, which can act as a platform to enhance the ligation of sheets with the CNTs and improve the structural stability of the hybrid carbon foam. The MWCNTs (or SWCNTs) were well-dispersed on or between the GO sheets and effectively combined with the GO sheets. The MWCNTs (or SWCNTs) played a bridge role for the entire hybrid carbon foam, effectively connecting the sheet-like structure of GO, which had a positive impact on the preparation of the composite PCMs. Therefore, the unique structure of the hybrid carbon foam was constructed by the interaction of the sheet-like GO and MWCNTs (or SWCNTs), which will improve the shape stability of the composite PCMs. As shown in Figure 5c,f, the hybrid carbon foam was filled with PEG after the impregnation. There was no terrible surface fracture or cracking in the cross-sections of the composite PCMs. This indicated that the PEG and hybrid carbon foam had good compatibility.  The microstructures of the hybrid carbon foams with different MWCNTs (or SWCNTs)/GO volume ratios and their composite PCMs were observed by FE-SEM.
As shown in Figure 5a,b,d,e, a porous structure of the carbon foams is found. The GO presents a sheet-like wrinkled surface, which can act as a platform to enhance the ligation of sheets with the CNTs and improve the structural stability of the hybrid carbon foam. The MWCNTs (or SWCNTs) were well-dispersed on or between the GO sheets and effectively combined with the GO sheets. The MWCNTs (or SWCNTs) played a bridge role for the entire hybrid carbon foam, effectively connecting the sheet-like structure of GO, which had a positive impact on the preparation of the composite PCMs. Therefore, the unique structure of the hybrid carbon foam was constructed by the interaction of the sheet-like GO and MWCNTs (or SWCNTs), which will improve the shape stability of the composite PCMs. As shown in Figure 5c,f, the hybrid carbon foam was filled with PEG after the impregnation. There was no terrible surface fracture or cracking in the cross-sections of the composite PCMs. This indicated that the PEG and hybrid carbon foam had good compatibility.

Macroscopic Phase Change of PEG and Composite PCMs
The macroscopic phase changes of the composite PCMs are very important for preventing the leakage of materials.

Macroscopic Phase Change of PEG and Composite PCMs
The macroscopic phase changes of the composite PCMs are very important for preventing the leakage of materials. The shape-stable properties of the pure PEG and the composite PCMs were observed in leakage tests using a hot plate and a digital camera. The pure PEG and the composite PCMs samples were heated from 70 • C to 100 • C and maintained for 10 min with a step increase of 10 • C.
As shown in Figure 6, the pure PEG began to melt gradually at 70 • C, and its solid-liquid phase was formed with an increase in temperature, which caused leakage in the pure PEG sample. In contrast, the shape-stable composite PCMs were maintained at a temperature above 90 • C. PEG/SCF-1.1 and PEG/MCF-1.1 started to melt, but the leakage could not wet the filter paper at 80 • C, while the other composite PCMs retained their original shapes. When heated to 90 • C and held there for 10 min, PEG/SCF-1.1 and PEG/MCF-1.1 presented small amounts of leaked liquid, while PEG/SCF-0.9 and PEG/MCF-0.9 remained stable, with only a little melting on the bottom. When the temperature reached 100 • C for 10 min, only a little melt was found for PEG/SCF-0.7 at such a high temperature. Additionally, PEG/SCF-0.5, PEG/SCF-0.3, PEG/MCF-0.5, PEG/MCF-0.3, and PEG/MCF-0.7 did not melt throughout the heating process as a result of the unique structure of the hybrid carbon foam, which provided enough micropores to facilitate the PEG impregnation and effectively prevent the leakage of PEG. It is worth noting that the shape of the PEG/MCF composite PCMs were more stable than those of the PEG/SCF composite PCMs, which could have resulted from the mechanical property of the multiple layers of MWCNTs being more stable than SWCNTs. Thus, the phase change materials could retain a steady morphology under a relatively high temperature and prevent the PEG from escaping through encapsulation. The shape-stable properties of the pure PEG and the composite PCMs were observed in leakage tests using a hot plate and a digital camera. The pure PEG and the composite PCMs samples were heated from 70 °C to 100 °C and maintained for 10 min with a step increase of 10 °C.
As shown in Figure 6, the pure PEG began to melt gradually at 70 °C, and its solid-liquid phase was formed with an increase in temperature, which caused leakage in the pure PEG sample. In contrast, the shape-stable composite PCMs were maintained at a temperature above 90 °C . PEG/SCF-1.1 and PEG/MCF-1.1 started to melt, but the leakage could not wet the filter paper at 80 °C , while the other composite PCMs retained their original shapes. When heated to 90 °C and held there for 10 min, PEG/SCF-1.1 and PEG/MCF-1.1 presented small amounts of leaked liquid, while PEG/SCF-0.9 and PEG/MCF-0.9 remained stable, with only a little melting on the bottom. When the temperature reached 100 °C for 10 min, only a little melt was found for PEG/SCF-0.7 at such a high temperature. Additionally, PEG/SCF-0.5, PEG/SCF-0.3, PEG/MCF-0.5, PEG/MCF-0.3, and PEG/MCF-0.7 did not melt throughout the heating process as a result of the unique structure of the hybrid carbon foam, which provided enough micropores to facilitate the PEG impregnation and effectively prevent the leakage of PEG. It is worth noting that the shape of the PEG/MCF composite PCMs were more stable than those of the PEG/SCF composite PCMs, which could have resulted from the mechanical property of the multiple layers of MWCNTs being more stable than SWCNTs. Thus, the phase change materials could retain a steady morphology under a relatively high temperature and prevent the PEG from escaping through encapsulation.

Crystallization Behavior of the PEG and the Composite PCMs
To discuss the effects of CF on the morphology of the composite PCMs, the corresponding characterization was carried out. The structures of the MCF (or SCF) with different MWCNTs (or SWCNTs)/GO volume ratios, pure PEG, and the composite PCMs were characterized by XRD, as shown in Figure 7.

Crystallization Behavior of the PEG and the Composite PCMs
To discuss the effects of CF on the morphology of the composite PCMs, the corresponding characterization was carried out. The structures of the MCF (or SCF) with different MWCNTs (or SWCNTs)/GO volume ratios, pure PEG, and the composite PCMs were characterized by XRD, as shown in Figure 7.
The typical diffraction peaks of the GO [24,25]. The characteristic diffraction peak of the PEG that was in the PCMs was presented, while the characteristic diffraction peaks of the GO and MWCNTs (or SWCNTs) disappeared. PEG can completely package these MWCNTs (or SWCNTs) in the MCF (SCF). When the PEG/MCF (PEG/SCF) composites are performed by XRD, only the PEG diffraction peaks can be seen. Notably, the PEG diffraction peaks at around 2θ = 31.50 • are visible, suggesting that MWCNTs (SWCNTs) can induce PEG crystallization. As can be seen from Figure 7a, an increase in the content of MWCNTs causes a less sharp characteristic diffraction peak of the PEG/MCF composite PCMs. This phenomenon can be explained by considering the MWCNTs as an impurity that has entered the PEG and hindered the growth of the PEG crystal itself. Compared with the MWCNTs, the characteristic diffraction peak of the PEG/SCF composite PCMs becomes more pointed with an increase in the loading of SWCNTs. This indicates that SWCNTs, as a nucleating agent, enhance the crystallization ability of PEG. The typical diffraction peaks of the GO and CNTs were observed at approximately 2θ = 10.30°, 21.03°, and 25.50°, respectively, while those of pure PEG appeared at 2θ = 19.30°, 23.30°, 26.20°, and 27.20°, respectively [24,25]. The characteristic diffraction peak of the PEG that was in the PCMs was presented, while the characteristic diffraction peaks of the GO and MWCNTs (or SWCNTs) disappeared. PEG can completely package these MWCNTs (or SWCNTs) in the MCF (SCF). When the PEG/MCF (PEG/SCF) composites are performed by XRD, only the PEG diffraction peaks can be seen. Notably, the PEG diffraction peaks at around 2θ = 31.50° are visible, suggesting that MWCNTs (SWCNTs) can induce PEG crystallization. As can be seen from Figure 7a, an increase in the content of MWCNTs causes a less sharp characteristic diffraction peak of the PEG/MCF composite PCMs. This phenomenon can be explained by considering the MWCNTs as an impurity that has entered the PEG and hindered the growth of the PEG crystal itself. Compared with the MWCNTs, the characteristic diffraction peak of the PEG/SCF composite PCMs becomes more pointed with an increase in the loading of SWCNTs. This indicates that SWCNTs, as a nucleating agent, enhance the crystallization ability of PEG.
A polarizing microscope (POM) was another method used to characterize the crystallization of the composite PCMs. POM photographs of the pure PEG and the composite PCMs are shown in Figure 8. A polarizing microscope (POM) was another method used to characterize the crystallization of the composite PCMs. POM photographs of the pure PEG and the composite PCMs are shown in Figure 8. The typical diffraction peaks of the GO and CNTs were observed at approximately 2θ = 10.30°, 21.03°, and 25.50°, respectively, while those of pure PEG appeared at 2θ = 19.30°, 23.30°, 26.20°, and 27.20°, respectively [24,25]. The characteristic diffraction peak of the PEG that was in the PCMs was presented, while the characteristic diffraction peaks of the GO and MWCNTs (or SWCNTs) disappeared. PEG can completely package these MWCNTs (or SWCNTs) in the MCF (SCF). When the PEG/MCF (PEG/SCF) composites are performed by XRD, only the PEG diffraction peaks can be seen. Notably, the PEG diffraction peaks at around 2θ = 31.50° are visible, suggesting that MWCNTs (SWCNTs) can induce PEG crystallization. As can be seen from Figure 7a, an increase in the content of MWCNTs causes a less sharp characteristic diffraction peak of the PEG/MCF composite PCMs. This phenomenon can be explained by considering the MWCNTs as an impurity that has entered the PEG and hindered the growth of the PEG crystal itself. Compared with the MWCNTs, the characteristic diffraction peak of the PEG/SCF composite PCMs becomes more pointed with an increase in the loading of SWCNTs. This indicates that SWCNTs, as a nucleating agent, enhance the crystallization ability of PEG.
A polarizing microscope (POM) was another method used to characterize the crystallization of the composite PCMs. POM photographs of the pure PEG and the composite PCMs are shown in Figure 8.  These images show that the change in the crystalline behavior of PEG is affected by the change in the volume ratio of carbon nanotubes and GO. It is observed that the spherulite morphology of the pure PEG has a unique black cross-light pattern and uniform distribution (Figure 8a). Compared with pure PEG, when the content of MWCNTs (or SWCNTs) increases, the crystallization shape of the PEG/MCF (or SCF) composite is smaller and more fragmented until it becomes a string of small spherulite crystals. This phenomenon shows that the CNTs and GO, as a kind of impurity and/or nucleating agent, destroy the spherulite morphological characteristics of PEG, reduce the PEG crystallization rate within the limited space, and limit the molecular diffusion motion and conformational change in the molecules during crystal growth.

Thermal Energy Storage of the PEG and the Composite PCMs
By means of the DSC test method, the phase change temperature and latent heat values of the pure PEG and the composite PCMs were tested, and their ability to store heat energy was investigated as shown in Figure 9. The results of a data analysis of the initial melting/ crystallization temperature (T ms /T cs ), end melting/crystallization temperature (T me /T ce ), peak melting/crystallization temperature (T mp /T cp ), and melting/crystallization enthalpy (∆H m /∆H c ) are listed in Table 1. The results of a data analysis of 150 cycles are listed in Table 2. The degree of crystallinity of the composite PCMs is listed in Table 3. These images show that the change in the crystalline behavior of PEG is affected by the change in the volume ratio of carbon nanotubes and GO. It is observed that the spherulite morphology of the pure PEG has a unique black cross-light pattern and uniform distribution (Figure 8a). Compared with pure PEG, when the content of MWCNTs (or SWCNTs) increases, the crystallization shape of the PEG/MCF (or SCF) composite is smaller and more fragmented until it becomes a string of small spherulite crystals. This phenomenon shows that the CNTs and GO, as a kind of impurity and/or nucleating agent, destroy the spherulite morphological characteristics of PEG, reduce the PEG crystallization rate within the limited space, and limit the molecular diffusion motion and conformational change in the molecules during crystal growth.

Thermal Energy Storage of the PEG and the Composite PCMs
By means of the DSC test method, the phase change temperature and latent heat values of the pure PEG and the composite PCMs were tested, and their ability to store heat energy was investigated as shown in Figure 9. The results of a data analysis of the initial melting/crystallization temperature (Tms/Tcs), end melting/crystallization temperature (Tme/Tce), peak melting/crystallization temperature (Tmp/Tcp), and melting/crystallization enthalpy (ΔHm/ΔHc) are listed in Table 1. The results of a data analysis of 150 cycles are listed in Table 2. The degree of crystallinity of the composite PCMs is listed in Table 3.     The endothermic peak of the composite PCMs is about 64 • C, and the exothermic peak is about 42 • C, which is very similar to the endothermic/exothermic peak range of PEG. As shown in Table 1, compared with pure PEG, the melting and crystallizing temperature regions of the composite PCMs are obviously expanded. It is considered that the GO and CNTs act as a nucleator to promote the crystallization of PEG and provide a thermal conduction pathway in the composite PCMs to accelerate the phase change process. The DSC test results show that the endothermic (∆H m = 135.16~155.09 J/g) and exothermic (∆H c = 131.38~149.83 J/g) function values of the composite PCMs are slightly lower than those of the pure PEG (∆H m = 173.60 J/g, ∆H c = 160.97 J/g). The ∆H m /∆H c values of the composite PCMs gradually decrease with an increase in the CNTs content because an excessive amount of nanotubes could destroy the molecular structure of PEG, which further impacts the re-crystallization and melting of the composite PCMs. It can be observed from Table 2 that the measured values are almost unchanged after 150 cycles of cooling and heating. This proves that the composite PCMs have good thermal cycling stability.

Melting Process Solidifying Process
As a key parameter of the composite PCMs, the degree of crystallinity is closely related to the thermal energy storage. The value of degree of crystallinity X can be obtained from Equation (1): where ∆H 0 is the melting enthalpy value of pure PEG when it is tested by DSC, and ∆H 0 = 173.6 J/g, w is the content of PEG in the composite PCMs. As can be seen from Table 3, the degree of crystallinity of the composite PCMs showed a gradual decline with the increase of the content of carbon nanotubes. This is because carbon nanotubes, as a kind of impurity and/or nucleating agent, limit the molecular diffusion motion and conformational change in the molecules during crystal growth, and reduce the PEG crystallization rate within the limited space.
The thermal stability levels of pure PEG and the composite PCMs were characterized using TG as shown in Figure 10.   As can be seen from Table 3, the degree of crystallinity of the composite PCMs showed a gradual decline with the increase of the content of carbon nanotubes. This is because carbon nanotubes, as a kind of impurity and/or nucleating agent, limit the molecular diffusion motion and conformational change in the molecules during crystal growth, and reduce the PEG crystallization rate within the limited space.

Samples Samples
The thermal stability levels of pure PEG and the composite PCMs were characterized using TG as shown in Figure 10. As can be seen from Figure 10a, the PEG/MCF composite PCMs have maintained a good thermal stability level with almost no mass loss under 250 °C, which shows that the PEG/MCF composite PCMs have high thermal stability levels in the range of the phase change temperature and can be used in energy storage. However, it can be seen from Figure 10b that the thermal stability levels of the PEG/SCF composite PCMs are worse than those of the PEG/MCF composite PCMs. This is because the unique interlayer structure of MWCNTs can enhance stability of CF.
The mass loss of the composite PCMs is over 90% from Table 4. It is illustrated that CF provides a huge network space that increases the amount of PEG injection. So, the leakage of PEG of the composite PCMs is prevented and the thermal stability of the PCMs is improved.  As can be seen from Figure 10a, the PEG/MCF composite PCMs have maintained a good thermal stability level with almost no mass loss under 250 • C, which shows that the PEG/MCF composite PCMs have high thermal stability levels in the range of the phase change temperature and can be used in energy storage. However, it can be seen from Figure 10b that the thermal stability levels of the PEG/SCF composite PCMs are worse than those of the PEG/MCF composite PCMs. This is because the unique interlayer structure of MWCNTs can enhance stability of CF.
The mass loss of the composite PCMs is over 90% from Table 4. It is illustrated that CF provides a huge network space that increases the amount of PEG injection. So, the leakage of PEG of the composite PCMs is prevented and the thermal stability of the PCMs is improved. Thermal conductivity is also relevant to energy storage applications because it directly affects the energy absorption and energy conversion.
As shown in Figure 11, the thermal conductivity of pure PEG is 0.289 W/(mK), while the thermal conductivities of the PEG/MCF composite PCMs and the PEG/SCF composite PCMs reached as high as 1.535 W/(mK) and 1.159 W/(mK), respectively, and the mean error of thermal conductivity measurements is 4.43% and 5.97%, respectively, which indicates that the introduction of a hybrid carbon foam is very effective in enhancing the thermal conductivity. It is clear that the network, by connecting the MWCNTs (or SWCNTs) and the GO nanosheets to each other, improves the thermal conductivity of the composite PCMs. The thermal conduction in the PEG/MCF and the PEG/SCF was five and four times that in the pure PEG, respectively. It can be understood that multiple layers of MWCNTs could more efficiently bridge the oxidized graphene, and this thermal conduction bridge could effectively promote the thermal conduction property of the material. The PEG/MCF-0.7 composite PCM was chosen to investigate the thermal conduction at various temperatures after 150 cycles because it has the highest thermal conductivity and the greatest shape stability among the PCMs prepared. It is interesting that there is almost no change that takes place on the thermal conduction property of the composite PCM, and the original shape of the PCMs remains unchanged without PEG leakage, indicating that the PCMs have steady thermal conduction performance, as shown in Figure 11b. Thermal conductivity is also relevant to energy storage applications because it directly affects the energy absorption and energy conversion.
As shown in Figure 11, the thermal conductivity of pure PEG is 0.289 W/(mK), while the thermal conductivities of the PEG/MCF composite PCMs and the PEG/SCF composite PCMs reached as high as 1.535 W/(mK) and 1.159 W/(mK), respectively, and the mean error of thermal conductivity measurements is 4.43% and 5.97%, respectively, which indicates that the introduction of a hybrid carbon foam is very effective in enhancing the thermal conductivity. It is clear that the network, by connecting the MWCNTs (or SWCNTs) and the GO nanosheets to each other, improves the thermal conductivity of the composite PCMs. The thermal conduction in the PEG/MCF and the PEG/SCF was five and four times that in the pure PEG, respectively. It can be understood that multiple layers of MWCNTs could more efficiently bridge the oxidized graphene, and this thermal conduction bridge could effectively promote the thermal conduction property of the material. The PEG/MCF-0.7 composite PCM was chosen to investigate the thermal conduction at various temperatures after 150 cycles because it has the highest thermal conductivity and the greatest shape stability among the PCMs prepared. It is interesting that there is almost no change that takes place on the thermal conduction property of the composite PCM, and the original shape of the PCMs remains unchanged without PEG leakage, indicating that the PCMs have steady thermal conduction performance, as shown in Figure 11b. The dimensional stability of the composite PCMs was tested by a thermomechanical analyzer. It can be seen from the Figure 12 that there is an inflection point in the pure PEG at 60.2 °C , while there is an inflection point in the PEG/MCF-0.3 and PEG/SCF-0.3 at 81.3 °C and 78.4 °C , respectively. It can be explained that carbon foam encapsulation of PEG can effectively improve the dimensional stability of the composite PCMs. However, as the content of carbon nanotubes increases, the deformation temperature of the phase change material decreases. This is because the excessive introduction of CNTs destroys the frame structure of CF and limits the adsorption of PEG. Thus, the deformation temperature of the composite PCMs will gradually decline. Subsequently, the samples were subjected The dimensional stability of the composite PCMs was tested by a thermomechanical analyzer. It can be seen from the Figure 12 that there is an inflection point in the pure PEG at 60.2 • C, while there is an inflection point in the PEG/MCF-0.3 and PEG/SCF-0.3 at 81.3 • C and 78.4 • C, respectively. It can be explained that carbon foam encapsulation of PEG can effectively improve the dimensional stability of the composite PCMs. However, as the content of carbon nanotubes increases, the deformation temperature of the phase change material decreases. This is because the excessive introduction of CNTs destroys the frame structure of CF and limits the adsorption of PEG. Thus, the deformation temperature of the composite PCMs will gradually decline. Subsequently, the samples were subjected to a thermal mechanical cycle analysis 20 times, which further proved that the composite PCMs have a good dimensional stability.

Thermal Energy Conversion
In order to detect whether the composite PCMs had energy storage and conversion behaviors, a simulation of a solar-to-thermal energy conversion experiment for the storage and release of the composite PCMs was conducted. Samples with a size of 2 × 1 × 0.1 cm were maintained for 600 s under a simulated solar light source with the intensity of the xenon lamp at 100 mW/cm 2 .
As shown in Figure 13, during the entire optical radiation process, the composite PCMs absorbed rays. This was followed by an increase in temperature until the inflection point appeared (at about 80~360 s), which indicated that the heat storage effects of the composite materials had already occurred. When the light source was removed, the temperature of the composite PCMs began to drop until a large platform emerged, indicating that the composite PCMs were releasing energy in the process of recrystallization. The thermal energy storage efficiency (η) was estimated by the light radiation energy that was absorbed during the phase change. The data results can be calculated by the following equation [26]: where m is the mass of the PCMs material, ΔH is the enthalpy value of phase change, P is the intensity of the simulated light source, S is the surface area of the sample, and ts and te are the starting and ending time of phase change, respectively.

Thermal Energy Conversion
In order to detect whether the composite PCMs had energy storage and conversion behaviors, a simulation of a solar-to-thermal energy conversion experiment for the storage and release of the composite PCMs was conducted. Samples with a size of 2 × 1 × 0.1 cm were maintained for 600 s under a simulated solar light source with the intensity of the xenon lamp at 100 mW/cm 2 .
As shown in Figure 13, during the entire optical radiation process, the composite PCMs absorbed rays. This was followed by an increase in temperature until the inflection point appeared (at about 80~360 s), which indicated that the heat storage effects of the composite materials had already occurred. When the light source was removed, the temperature of the composite PCMs began to drop until a large platform emerged, indicating that the composite PCMs were releasing energy in the process of recrystallization. The thermal energy storage efficiency (η) was estimated by the light radiation energy that was absorbed during the phase change. The data results can be calculated by the following equation [26]: where m is the mass of the PCMs material, ∆H is the enthalpy value of phase change, P is the intensity of the simulated light source, S is the surface area of the sample, and t s and t e are the starting and ending time of phase change, respectively. se () mH η PS t t    (2) where m is the mass of the PCMs material, ΔH is the enthalpy value of phase change, P is the intensity of the simulated light source, S is the surface area of the sample, and ts and te are the starting and ending time of phase change, respectively. By calculation, the maximum thermal energy storage efficiency of the PEG/MCF composite PCMs and the PEG/SCF composite PCMs was 91.2% and 90.6%, respectively. The results are significantly better than those of Yang et al. [27]. Moreover, the thermal energy storage efficiency of the PEG/MCF composite PCMs is higher than that of the PEG/SCF composite PCMs. The explanation for this is that the interlayer structure of MWCNTs in the composite PCMs affects its energy storage as well as its thermal conductivity and morphological stability.

Conclusions
In this study, composite PCMs, such as PEG/MCF (or PEG/SCF), with outstanding shape-stabilization, good thermal conductivity, and superior energy conversion, simultaneously, were prepared by introducing a hybrid CF. The hybrid carbon foams were fabricated by freeze-drying of a mixture of GO and MWCNTs (or SWCNTs), which showed a unique structure to enhance the impregnation of PEG and prevent its leakage in composite PCMs. It was found that there is no leakage of the melted PEG at 100 • C for a series of composite PCMs, including PEG/MCF-0.3, PEG/MCF-0.5, PEG/MCF-0.7, PEG/SCF-0.3, and PEG/SCF-0.5. It was indicated that the PEG/MCF composite PCMs had better thermal stability below 250 • C, and the thermal conductivities of the PEG/MCF composite PCMs reached as high as 1.535 W/Mk, which was obviously higher than that of PEG/SCF. The results of the photo-thermal simulation tests showed that the composite PCMs with different CNTs had the ability to absorb light energy and convert it to heat energy. The maximum thermal energy storage efficiency of the PEG/MCF composite PCMs and the PEG/SCF composite PCMs was 91.2% and 90.6%, respectively, which has potential application in energy storage.