Open access peer-reviewed chapter

Experimental and Numerical Studies on Phase Change Materials

Written By

Cheng Wang and Ye Zhu

Submitted: October 30th, 2017 Reviewed: March 27th, 2018 Published: November 5th, 2018

DOI: 10.5772/intechopen.76807

Chapter metrics overview

1,548 Chapter Downloads

View Full Metrics


Phase change materials (PCMs) are attracting significant attentions in research and application, categorized into mainly three types, that is, organic (O), inorganic (IO) and eutectic (E). This section introduces the experimental and numerical investigations conducted in recent decades, mainly focused on the properties enhancement of PCMs and the performance improvement of its application in latent heat storage (LHS) units, as well as the evaluation and optimization of LHS units. It was concluded that lots of contribution have been made to PCMs and LHS units analysis. However, there is still some weakness in research, such as the lackness of detailed and systematic research on properties, the non-uniform standard on testing method as well as the contradictory conclusions. The most evaluation of LHS units is based on energy, instead of exergy, entropy and entransy. There is another issue that most of the research is based on numerical analysis, while less experimental research is conducted, especially in the case of LHS unit.


  • phase change materials (PCMs)
  • latent heat storage (LHS)
  • numerical
  • experimental

1. Introduction

Energy is the basis of modern society and is important for the survival of mankind as well as the development of civilization. Non-renewable and renewable sources are two kinds of energy source. Since the non-renewable energy source will be run out someday, the utilization of renewable energy source has been paid more and more attention in research. However, in most cases of renewable energy sources, such as solar and wind, intermittent nature is found. What is more, there is always a gap between energy supply and energy demand, as far as power, space and time in concern. Thus, energy storage technologies are proposed to solve or diminish this issue.

Energy is usually stored in energy storage (ES) system, in the form of mechanical, chemical, biological, magnetic and thermal. These energy storage forms can also be subsidized further in details. For instance, mechanical energy can be stored in compressed air, flywheel or hydro-pool, etc. and chemical energy can be stored in battery, reversible-reaction or hydrogen, etc. Among these energy storage forms, the most commonly used is the thermal energy storage (TES) with phase change materials (PCMs), due to its merits of low-cost, environmental-friend, easy-to-operate and abundant sources of storage facilities.

As a matter of fact, we human being has used renewable energy and conducted thermal energy storage, since quite a long time ago. For example, the ancient people utilized wind or hydro power to drive wheels for irrigation and they collected ice or snow in winter for cooling in summer. In modern society, we try to fully utilize the clean energy source, to deeply understand the process involved in TES process and to seek nature materials or manufacture artificial materials for TES. For the performance improvement of TES, the thermophysical properties are important limitations. For instance, the limited thermal conductivity of PCMs strongly constrains the conductive heat transfer process. The viscosity of PCMs at liquid phase also constrains the convective heat transfer process. In this section, we will introduce some progress of the research on PCMs and TES and discuss on the weak points.


2. Experimental studies of PCMs

Experiments studies were conducted on the properties of PCMs and the performance in LHS, as well as the enhancement.

2.1. Types and properties of PCMs

The materials involved in LHS are called as phase change materials (PCMs). There are varieties of PCMs under development, categorized as organic (O), inorganic (IO) and eutectic (E) materials, available in a wide range of melting/solidification temperature (Figure 1).

Figure 1.

Categories of PCMs.

Desired property of PCMs includes thermodynamic, kinetic, physical and chemical properties, as well as economic availability, as shown in Table 1.

Thermal propertiesPhysical propertiesKinetic propertiesChemical propertiesEconomic availability
High latent heat of transitionSmall volume changeSufficient crystallization rateLong-term chemical stabilityAbundant
High thermal conductivityLow vapor pressureNo supercoolingNo toxicityCost effective
Suitable melting/solidification temperatureHigh densityNon-flammableAvailable
Non-corrosiveCommercially viable

Table 1.

Properties of PCMs in demand.

Unfortunately, there is not a single kind of PCMs that satisfies all the properties listed above. The most undesired character of PCMs property is its thermal conductivity, since it will limit the heat transfer during energy storage/release process and correspondingly lead to deteriorated performance of LHS unit. This is often the case, except for some metallic PCMs. For instance, several measures have been taken. The major technical method is to composite with materials of high thermal conductivity.

2.2. Types of TES

TES can be categorized as three types, that is, sensible heat storage (SHS), latent heat storage (LHS) and thermo-chemical storage (TCS). In the first type (SHS), thermal energy is stored as the temperature increase of certain matters, usually with large thermal capacity. So, the amount of energy storage Q can be easily estimated as the product of mass m, thermal capacity C and temperature uplift ΔT, as shown:


In the second type (LHS), thermal energy is stored as the phase change process of certain materials, including the transformation of phase between solid and liquid (S-L) in melting/solidification process, between solid and gas (S-G) in sublimation/desublimation process, from liquid to gas (L-G) in condensation/evaporation process, as well as the transition from one solid phase to another (S-S). The amount of energy storage Q is the sum of the sensible heat stored in both phases and the latent heat involved in phase-transformation, which is the main portion of energy storage amount.


In the third type (TCS), thermal energy is stored in similar way as LHS. The major difference is that thermal energy is mainly stored as the enthalpy change in thermo-chemical reaction, instead of phase-transformation process.

Comparing with LHS and TCS, SHS technology often requires larger vessels. Comparing with LHS, TCS is associated with larger energy storage density, but is still at pre-mature state in terms of research and development. Therefore, latent heat storage (LHS) attracts the most attention in research and is believed as the most promising technology.

The performance of heat storage/release in PCMs is realized in LHS unit. The major favored characteristics of LHS unit includes faster rate of heat charging/discharging, higher efficiency of heat release, based on thermodynamic evaluations, including the basis of energy, exergy, entropy and entransy as well. Another important research field is the optimization of LHS unit.

It is widely accepted that performance of LHS unit is mainly constrained by heat transfer process. Therefore, heat transfer enhancement is a major task for LHS unit performance improvement in most research. Since heat transfer is generously expressed as:


where K represents the heat transfer coefficient in conduction or convection process, A represents the surface area for heat transfer and ΔTm represents the temperature difference between PCMs and HTF. The enhancement of heat transfer implies the increase of Q. So, it is obvious that there are three major methods for heat transfer enhancement of LHS unit, that is, increase of K, A and ΔTm.

The heat transfer coefficient should include both conductive and convective. The increase of conductive heat transfer coefficient can be mainly traced back to the thermal conductivity enhancement on PCMs. The only exception would be the encapsulation of PCMs. As far as the increase of convective heat transfer coefficient, the theoretical basis is convective heat transfer. Therefore, the progress on convective heat transfer can be applied directly in the performance improvement of LHS unit, such as the influence of passage size and shape, the effects of faster flow velocity and the disturbed flow pattern. More effective increase of convective heat transfer coefficient should be attributed to the application of heat pipe (HP) technology.

2.3. Heat transfer enhancement techniques

Heat transfer for thermal energy storage applications with phase change materials is reviewed by Ref. [1]. The major measures include conductive heat transfer enhancement and convective heat transfer enhancement.

2.3.1. Composites with porous materials

Composite with porous materials is an effective method for thermal conductivity enhancement of PCMs. Impregnation is a fast-growing technology. The porous material offers space for PCMs and the high thermal conductivity of porous materials supports more effective heat transfer in composite. Expanded graphite and metal foam are the mostly adopted porous materials.

Zhao and Wu [2] reported the considerable improvement of thermal conductivity of sodium nitrate with porous expanded graphite matrix and metal foams. The experimental results of embedding non-metallic PCMs in porous graphite showed tens or hundreds times of thermal conductivity improvement. Siaphush et al. [3] reported the effective thermal conductivity increased from 0.423 W/m/K to 3.06 W/m/K, when 95% porosity copper foam is adopted in PCMs of eicosane. It is also reported that critical porosity value exists for enhancement with porous carbon graphite foams. Yin et al. [4] reported in experiments that the critical mass fraction of porous graphite is 6.25%. Exceeding this value, the reinforcement effect decreases. Similar phenomenon is also found in our research on the effect of adoption expanded graphite in octadecane for performance improvement of LHS unit. However, the critical value is found at around 20%. Gao et al. [5] investigated the thermal performance of a direct contact thermal energy storage container with erythritol and expanded graphite. The thermal conductivity is reported as increased by about 2.5 times, with 4% mass fraction of EG, and the melting time is reduced by 16.7%.

Comparing with the amount of porous materials in composite, it is also reported that the pore structure is more important for composite. Lafdi et al. [6] investigated experimentally on the effects of foam porosity and pore size on the melting rate of PCMs. Zhong et al. [7] reported that small pore size and thick ligament in graphite foam leads to higher thermal diffusivity, while large pore size and thin ligament leads to larger latent heat storage capacity. Since the thermal diffusivity and latent heat storage density are both important factors, the pore size and ligament should be optimized in the design of LHS unit. Wu and Zhao [8] reported that mixed porous base is more effective than single porous base. Zhang et al. [9] studied the performance of metal foam (copper) and paraffin composite. Gulfam et al. [10] investigated the enhancement of thermal conductivity with expanded graphite in paraffin wax. Teppei et al. [11] reported high thermal conductivity of erythritol enhanced with porous nickel. Similar works are also reported by Nomura et al. with metal-stabilized carbon-fiber network [12] as well as expanded perlite, diatom earth and gamma-alumina [13].

Besides expanded graphite and metal foams, ceramic is also adopted recently as enhancement medium, as Li et al. [14] reported. With the development of material science, there should be more materials with porous structure, such as graphene, aerogel, etc., under consideration for the enhancement of thermal conductivity as well as other properties.

2.3.2. Dispersion of high conductive particles

The effect of addition of metal particles, especially with nano-size, on the enhancement of thermal conductivity of PCMs is also reported widely. Different from the obvious and established structure of compressed expanded graphite and metal foam, the distribution of particles are more like to expanded graphite in composite. However, the effects are usually better in the case of nano-particles addition. The reinforcement effect by the dispersion of nano-particles in continuous PCMs should be attributed to the unique phenomenon at microscopic scale, for example, reduce the internal resistance for heat transfer, which is also reported for the thermal conductivity enhancement of heat transfer fluid (HTF) in literatures.

Fan et al. [15], Qi et al. [16], Tao et al. [17], Kim et al. [18] and Shi et al. [19] investigated the effects of several types of carbon-based nano-particles. It is found that disk carbon nano-particles can improve the thermal conductivity by 10 times. Besides graphite or graphene, carbon nanotubes (CNTs) are more typical nano-particles in composite. Zhang et al. [20] investigated on nanoscale heat transfer in composite of sugar alcohol and carbon nanotubes. It is reported that specific improvement of heat transfer depends on the material and the diameter of CNT. Carbon nano-fibers (CNFs) are another common nano-particles in composite. Fereshteh et al. [21] analyses the application of phase change material with carbon fibers for thermal management of a Li-ion battery cell. It is concluded that the application of carbon fibers influences the temperature distribution in cells. Higher concentration of carbon fibers leads to more uniform temperature distribution. The maximum thermal conductivity enhancement degree is reported as 115% and averaged at 105%. Nomura et al. [22] reported a significant reinforcement degree of thermal conductivity for erythritol. Different CNFs groups and its mixture are adopted. It is found that with the mixture of CNFs at different length, thermal conductivity is more enhanced, comparing with the addition of single CNFs. To further construct the network for heat transfer inside composite, low-melting metal, such as indium is added to help bridge CNFs nearby. In recent studies, other carbon materials are under research.

Besides the carbon materials, the addition of nano-metal-particles, such as Cu, Ti and other magnetic metals, in composite are also conducted by researchers Kibria et al. [23], Zhang et al. [24], Luo et al. [25], etc. It is concluded that the thermal conductivity is enhanced, and sometimes the thermal capacity is also enhanced. Mettawee et al. [26] reported the effect of Al powder on thermal conductivity enhancement of paraffin wax. Motahar et al. [27] reported non-monotonic behavior of thermal conductivity, and optimum value of nano-particles occurs in composite. Wang et al. [28] reported the increase of thermal conductivity with mass fraction of CNTs. Similar result of MWCNTs is also reported by Zeng et al. [29] for palmitic acid. Oya et al. [11] studied the thermal conductivity enhancement of erythritol with graphite and nickel particles. The largest enhancement is reported as 6.4 times higher, comparing with the thermal conductivity of pure phase change material, at 15% volume fraction of expanded graphite. Khyad et al. [30] adopted 1% mass fraction of aluminum or copper to enhance thermal conductivity of paraffin.

Since the nano-particle can enhance the thermal conductivity of PCMs with the similar mechanism, more research is expected on this scope, with the development of materials science on materials as well as the manufacture method.

2.3.3. Using extended surface

Surface area for heat transfer is the most common method applied for the heat transfer enhancement of LHS unit, mainly in the form of fin-structure. The adoption of fins increase the contact surface between HTF and PCMs. Research is mainly focused on the selection of fin materials as well as the configuration and number of fins in LHS unit. As far as fin materials are concerned, thermal conductivity, density, cost and corrosion as well as mechanical performance are major concerns. Recently, mostly metal, such as copper, bronze, steel, stainless steel, aluminum alloy, etc., and sometimes graphite or ceramic are used as fin materials.

The core of fin-structure is its configuration, including shapes and orientation. The performance of single structured-fin will lead to the number of fins in demand is influential to the configuration of fins in LHS unit. The investigation of fin configuration is similar to heat exchanger (HE) with almost constant temperature boundary. The typical structure of LHS unit is tube-and-shell. So, there will be two forms of PCMs arrangement, that is, inside of tube and outside of tube as well as annual space between tubes.

As far as the orientation of fins is concerned, there are two mainly forms, that is, alongside and perpendicular to the axial direction. The fins can be arranged inside and outside of tubes (Figure 2).

Figure 2.

Sketch of fin configurations. (a) Perpendicular to the axial direction and (b) alongside the axial direction.

Sparrow et al. [31] experimentally investigated the solidification process of PCMs in a finned vertical tube. It was concluded that conduction controls the process, when liquid temperature is lower and at melting temperature. While convection is the controlling mode for temperature above melting temperature. Tao et al. [32] investigated numerically with the performance of LHS unit in a photo-thermal (PT) application. Velraj et al. [33] reported with numerical and experimental analysis on vertical finned tube. The results show the reversal decrease of solidification period with number of fins. Zhao and Tan [34] investigated the effects of HTF temperature and flowrate as well as fin height on the charging rate of LHS unit. It is concluded that with the increase of HTF inlet temperature and mass flowrate, as well as the increase of fin height, the charging period is shortened, implying enhanced heat transfer process. Erek et al. [35] analyzed the effects of fin parameters, such as fin size and fin space, as well as the effects of HTF on the dimensionless energy storage value. Liu et al. [36] conducted similar experimental research on the melting of stearic acid in annual space. It is concluded that heat conduction and natural convection are both the factors for heat transfer enhancement in LHS unit. Hosseini et al. [37] concluded that with the increase of fin’s height, the reduction of melting time exhibits a descending trend. Comparing with the melting process, effects of increasing fin’s height is more significant in solidification process.

Besides the shell-and-tube configuration, numerical investigation on plate-type LHS is also conducted by Gharebaghi and Sezai [38] for rectangular heat sink. Sharifi et al. [39] developed an analytic model to predict the melting period of PCMs. Mahmoud et al. [40] conducted different arrangements for heat sink with PCMs at different melting temperature. It is concluded that increasing fin number is good for heat distribution in LHS unit and leading to lower the peak temperature of heat sink. Arshad et al. [41] studied the effects of pin thickness as well as the volume fraction of PCMs on the cooling performance of heat sinks for electronic devices. The volume fraction of PCMs is kept at 9%. The heating boundary is assumed as uniform heat flux. Heat sinks are finned or not finned. The thickness of fins is ranging from 1 mm to 3 mm, with the interval of 1 mm. The volume fraction of PCMs is ranging from 0 to 1, with the interval of 0.33. Rahimi et al. [42] investigated with the charging and discharging processes of PCMs in finned-tube heat exchanger in experiments. The utilization of fins reduces the melting and solidification periods. It is also reported that the increase of inlet temperature is more effective for melting time reduction in the bare tube heat exchangers. The variation of flow rate of HTF is also more intensely influential on the solidification time for the bare tube heat exchangers.

As reported in the review paper contributed by Nasiru et al. [1], the presence of fins improve the heat transfer during the phase change process, regardless of the make-up and geometry of the LHTES systems. However, few studies considered the effects of fin numbers on thermal response of the LHS unit. Although the trend is easy to find, the quantitative analysis will help the optimal design in practice.

2.3.4. Using multi-PCMs

The increase of ΔTm should be expressed more precisely as the uniform distribution of the temperature difference between PCMs and HTF in LHS unit, during charging and discharging process. The benefits can be not only evaluated with the energy basis, but also with the exergy/entropy basis as well as the entransy basis, which is proposed in the recent decade.

According to the demand of uniform temperature difference between PCMs and HTF, the melting temperature of PCMs should decrease alongside the flow direction of HTF in charging process and increase alongside the flow direction of HTF in discharging process. This is usually realized by the transverse of flow direction of HTF in two processes.

Fang and Chen [43] numerically investigated the effects of multiple PCMs on the performance of LHS unit. It is concluded that difference of melting temperature between multiple PCMs is crucial for performance improvement. Wang et al. [44] proposed a new concept of homogeneous phase change process using multiple PCMs to significantly decrease the melting/solidification periods. Cui et al. [45] numerically analyzed the structure with three types of PCMs for solar receiver. It is reported that the fluctuation of HTF outlet temperature is better controlled, comparing with single PCMs. More energy flowrate is also expected. Hu et al. [46] proposed a thermal storage system with frustum-shape. Along the flow direction of HTF, volume of PCMs change. Maximum five types of PCMs are adopted in the LHS unit. It is found effective even at small temperature difference. However, there is also some report about the asynchronous effects of multiple PCMs on the charging process and the discharging process, by Kurnia et al. [47]. It is concluded that the arrangement of PCMs with high melting temperature at the inlet of HTF would improve the heat transfer in discharging process, but may slightly worsen the charging process.

Thus, the major factor for multiple PCMs design is the match of melting temperature of PCMs in LHS unit. To better understand this issue and to provide guidance for the design of LHS unit, optimization of multi-stage LHS unit with multiple PCMs is conducted. Since there is no heat-and-work conversion during the operation of LHS unit, entransy theory is also adopted for optimization, besides the traditional energy and exergy/entropy analysis.

Tao et al. [48] reported the melting temperature match for double-stage LHS unit in charging process. Zhao et al. [49] reported the melting temperature match for multi-stage LHS unit in charging process. Wang et al. [50, 51, 52] reported the optimized match of melting temperature and surface for heat transfer of double- and multi-stage LHS unit in charging and cycle period.

However, less attention has been paid to the transient process optimization as well as other factors influencing the operation of LHS unit. Moreover, the comparison between entransy analysis and exergy analysis is important. Works are undergoing in our group. It is found that the difference between optimum melting temperature of nearby PCMs is constant in entransy analysis, while the ratio is constant in entropy analysis. The detailed discussion will be made. However, still the reason lies there, not so easy to answer, although we know that is superficially due to the difference between the optimization goals.

2.3.5. Encapsulation

Encapsulation of PCMs is also an effective method for heat transfer improvement in PCMs region. The mechanism may be explained as the reduce of heat transfer path as well as the increase of surface conducting heat transfer. Encapsulation of PCMs is to disperse PCMs in LHS unit into groups of small-sized particles closed and surrounded by other materials or the derivatives of PCMs itself after procedure of treatment. So, the direct property of PCMs is actually not changed, and the benefits are mainly contributed to the performance improvement of LHS unit as discussed later. The main research lies on the selection of raw-material and the method of encapsulation, as reviewed by Jacob and Bruno [53], Liu et al. [54], Saman et al. [55], Liu et al. [56], etc.

Jamekhorshid et al. [57] and Su et al. [58] reviewed the microencapsulation methods of PCMs. Milian et al. [59] reviewed on specific encapsulation techniques for inorganic phase change materials and the thermophysical properties. Sketch of encapsulated PCMs is expressed in Figure 3. The shell can be single layered or multiple layered or linked matrix, and the core can be single zone or several isolated zones. The shape could be regular, such as spherical, tubular or oval and irregular.

Figure 3.

Sketch of encapsulated PCMs shapes. (a) Mononuclear; (b) polynuclear; (c) multi-wall; (d) matrix.

The methods of encapsulation are summarized in Table 2.

Physical methodsChemical methodsPhysic-chemical methods
Pan coatingInterfacial polymerizationIonic gelation
Air-suspension coatingSuspension polymerizationCoacervation
Centrifugal extrusionEmulsion polymerizationSol-gel
Vibration nozzle
Spray drying
Solvent evaporation

Table 2.

Methods of encapsulation of PCMs.

As far as the shell is concerned, Jacob and Bruno [53] reviewed on the shell materials in the encapsulation for high temperature thermal energy storage. Steel, nickel, sodium silicate, silicon dioxide, calcium carbonate and titanium dioxide are identified as shell materials. It is better to further consider the corrosion between shell materials and PCMs encapsulated, which is important for long-term stability as well as cost reduction. Ma et al. [60] reported an encapsulated metallic phase change materials. The shell material is iron and the core material is copper. The preparation is based on aerodynamic levitation method. It is concluded that the morphology evolution is attributed to the combined effects of liquid phase fraction of two not-miscible liquids, Stokes and Marangoni velocities of droplets, as well as the rotation direction of particles in the solidification process. Chen et al. [61] reported the preparation of nanocapsules. The core PCMs is n-dodecanol and the encapsulation method is miniemulsion polymerization with polymerizable emulsifier. The diameter is measured as 150 nm and the phase change temperature is 18.2°C. Yang et al. [62] proposed a hybrid elastomeric spherical structure. It is composed of foam core and solid shell. The performance is predicted with numerical investigation.

2.3.6. Application of heat pipe

Heat pipe is a thermal carrier to transfer heat from hot medium to cold medium spatially separated. Heat pipe has its own working fluid, flowing inside at closed or open mode. At the end of hot medium, the liquid phase working fluid evaporates and flows to the end of cold medium, where the gas phase working fluid condensates and flows back to the end of hot medium and makes a cycle. Since phase change is involved, heat pipe usually can supply better performance of heat transfer between hot medium and cold medium.

There are two operation modes of HP in LHS unit. The first is simultaneous heating for discharging and cooling for charging, and the other is intermittent heating and cooling. To some extent, the latter mode is easier to understand, and the former mode is better for the power match and good for continuous operation. In the intermittent mode, PCMs operates as the hot end of HP in discharging process of LHS unit and as the cold end of HP in charging process.

Shabgard et al. [63] developed a thermal network model to investigate the performance of cascaded PCMs and conducted exergy analysis. Shabgard et al. [64] considered the transient response of HP-assisted LHTES system with a 2D model. It is concluded that HP spacing is the key parameter for LHS unit design and controls the dynamic response of the system. Robak et al. [65] experimentally investigated the performance of HP-assisted LHTES system. It is concluded that with the assistance of HP, heat transfer during discharging process is almost twice improved. Bergman et al. [66] numerically investigated the performance of LHTES with HP in solar thermal power plant and reported increased charging/discharging rate of PCMs for two kinds of HTF flow pattern. Nithyanandam and Pitchumani [67] conducted a similar computational analysis on 3D physical model. In another work of Nithyanandam and Pitchumani [68], dynamic performance behavior of HP-assisted LHTES system is investigation with the consideration of cyclic operation.

However, most of the research is based on gravity-assisted HP. With the development of HP technology, other kinds of HP should also be considered for application in LHS unit. Moreover, it is found that most of the HP-assisted LHTES system is analyzed numerically and less attention has been paid to the experimental analysis. This will be an open field for research in the future.

2.3.7. Combined heat transfer enhancement techniques

With two or more techniques, such as the combination of fin-structure and heat pipe, or the combination of multiple PCMs and heat pipe, it is expected possible to further improve performance of LHS unit. Jung and Boo [69] numerically investigated the transient behavior of a LHS unit with fin-structured heat pipe. They used a row-by-row analysis method to estimate the layer necessary for system design. It is concluded that the increase of pitch would help increase heat transfer rate. Khalifa et al. [70] compared the performance of bare heat pipe and finned heat pipe. It is concluded that with fin-structure, energy efficiency is improved significantly. Nithyanandam and Pitchumani [71] conducted numerical analysis on LHSTES system with metal foam and embedded heat pipe. It is reported that the augmentation in heat transfer rate during charging decreased with pore-density of metal foam, due to the restriction in the formation of buoyancy-induced convection.


3. Numerical studies of PCMs

Although the results of numerical analysis are not always the case in practice, it offers an important way to investigate the process as well as the performance of phase change materials as well as latent heat storage units, characterized with less cost and short time occupation as well as convenience of parameter adjustment. The focus and the core in numerical analysis rely on the model used as well as the verification and modification of numerical models with experimental results.

3.1. Numerical models

Esen et al. [72] applied two models to describe the diurnal transient behavior of energy storage tanks. In the first model, HTF is flowing outside of pipe, the inside of which is filled with PCMs. In the second model, HTF is flowing inside of pipe, the outside of which is surrounded with PCMs. Two-dimensional analysis is conducted with enthalpy-based method, coupled with convective heat transfer between HTF and PCMs. The effects of properties of PCMs, parameters of geometry, such as the radius and height of cylinder or pipe, and characters of HTF, such as velocity and inlet temperature, on the melting time are discussed (Figure 4).

Figure 4.

Configuration of HTF passage. (a) Outside of tube; (b) inside of tube; (c) annual between tubes.

The expression of Nu is listed as:

For heat transfer inside tube:


For heat transfer outside tube:


Xia et al. [73] analyzed the heat transfer of latent thermal energy storage (LTES) system based on the effective packed bed model. The flow field is simplified as the flow through voids of a bed packed with PCMs particles. The random packing model is proposed for better simulation. The material properties and the thickness of encapsulation are two major factors for the heat transfer performance of a LTES bed.

The porosity is listed as:


Izquierdo-Barrientos et al. [74] presented a dimensionless numerical model for the evolution of enthalpy with temperature, instead of constant phase change temperature assumption.

The dimensionless parameters include:


Modi and Perez-segarra [75] developed a one-dimensional numerical model for a single-tank thermocline thermal storage system in the concentrated solar power plant. The influence of types of heat transfer fluid, the temperature difference stored in HTF as well as the cycle cut-off on system performance is investigated. Two aspects are taken as criterion for comparison, that is, cyclic behavior of the system and the time required for equilibrium state attainment.

The heat transfer coefficient is expressed as:


Opitz and Treffinger [76] developed a general heterogeneous model of heat transfer in packed beds. Lumped element formulation is implemented. The results are verified with two different experiments cited from references. No necessary to calibrate theoretical model with experiment results is reported.

The pressure drop for one layer of the packed bed is expressed as:


The heat transfer is summarized as:

With correlation of Coutier and Farber:


With correlation of Gnielinski:


Amin et al. [77] developed an effectiveness-NTU model of a thermal storage system with packed bed of encapsulated PCMs with the sphere shape. The two-dimensional representation is proposed to predict the heat transfer during phase change, comparing with one-dimensional phase change assumption in other configurations. A new definition of thermal resistance between HTF and PCMs is developed, taking the phase change process into consideration.

The heat transfer is expressed as:


Karthikeyan and Velraj [78] compared several mathematical models for numerical investigation of packed bed with encapsulated spherical PCMs. The enthalpy formulation technique is used to accommodate the phase change behavior of paraffin. Fully explicit finite difference method is adopted for solving numerical models. It is reported that the validity of model depends on the kind of HTF.

The governing equations are listed as:

For model 1:


For model 2:


For model 3:


3.2. Numerical simulation

Pakrouh et al. [79] present a numerical investigation on geometric optimization of heat sinks. Paraffin is selected as PCMs and aluminum is adopted as materials for heat sink and fins. The optimization parameters include the number of fins, the height of fins and the thickness of fins as well as the base. Natural convection is also taken into consideration. It reported a complex relation between PCMs and the volume percentage of thermal conductivity enhancers (TCEs). Shmueli et al. [80] investigated numerically with melting of PCMs in a vertical cylindrical tube and compared with experiments. The model is based on enthalpy-porosity formulation. The effects of parameters, such as the term describing the mushy zone in the momentum equation and the pressure-velocity coupling as well as pressure discretization schemes, are examined. No difference between PISO and SIMPLE schemes is found, while there is considerable difference between PRESTO and Body-Force-Weighted schemes. Local heat transfer and melting are compared and verified for numerical results. It is concluded that at the beginning of melting process, the heat transfer is mainly in the form of conduction in solid phase; while at the end of this process, the heat transfer is dominated by convection in liquid phase. Cascetta et al. [81] utilized FLUENT software to simulate the flow and heat transfer in an axisymmetric tank of cylindrical shape. Incompressible turbulent flow and fully developed forced convection is adopted in two-phase transient (LTNE-local thermal non-uniform) model to calculate the temperatures of fluid and solid phases. The porosity of filled bed is also considered variable in the radial direction and the thermal properties of both phases are related to temperature. The results agree well with experiments. Sciacovelli et al. [82] used enthalpy method to analyze the phase change phenomenon. Natural convection is neglected, due to the fully resolved fluid flow in the liquid phase. The evaluation of melting front as well as the temperature and velocity fields is studied in details. However, it is concluded that natural convection significantly affects the phase change process. Also in this paper, the effects of enhancement of thermal conductivity with the adoption of highly conductive nano-particles in PCMs are considered. Augment of thermal performance is found, due to the application of nano-particles. The melting time is reduced by 15% with 4% volume fraction of nano-particles. Similar results are also found for the heat transfer performance.


4. Conclusions

For the properties of PCMs under research, besides thermal conductivity and phase equilibrium, others such as supercooling [83], corrosion [84] and transportation [85] are also characterized and discussed. However, less attention has been paid on the systematic discussion. This is partly due to the diverse results in different research groups, and sometimes the conclusions are contradictory. For instance, Teng et al. [86] reported the advantage of multi-wall carbon nanotubes (MWCNTs) over graphite for effective enhancement of thermal conductivity. However, Choi et al. [87] reported the contrary conclusion. Another case is the reported results of the same method and the same materials at different ages or by different groups are sometimes at significant variations. For instance, the heat of fusion for Paraffin Wax is reported as 173.6 kJ/kg [88] and 266 kJ/kg [89]; the melting temperature of myristic acid is reported as 49–51°C [90] and 58°C [91].

One of the reasons lies on the lackness of uniform standard or detailed information of preparation, manufacture and raw materials as well as the diversified methods of properties measurements. As far as thermal conductivity is concerned, researchers can utilize stationary and non-stationary methods to measure. Even in detailed non-stationary measurement, point-, linear- or surface- heating source is available for choice. Therefore, it seems difficult to collect the results in reference to obtain the regular of physical properties for theoretical estimation or analysis.

In this section, we mainly introduce the progress of property enhancement of PCMs and performance improvement of LHS unit. The detailed information of reported results is referred to cited references. It is found that lots of work has been done in the past decades and great progress has been made. However, there is still some weakness in research. For PCMs, most research is based on experimental measurement of properties, and less attention has been paid on the regular summary for theoretical estimation in the future or optimal design of composite material as well as energy storage unit. As far as LHS unit is concerned, it is the opposite condition, where most research is based on numerical analysis and less experimental research is conducted. This may lead to the deviation of the performance of LHS unit in application from the designed values, especially when the properties of PCMs as well as its composite are still not clear in details. What is more, the lackness of uniform standard and detailed report on information of preparation, manufacture and raw materials makes it difficult to collect the results of different groups and different ages all together.



This work is supported by Natural Science Foundation of China (51306023) and Advanced Catalysis and Green Manufacturing Collaborative Innovation Center.


Conflict of interest

None declared.


  1. 1. Ibrahim NI, Al-Sulaiman FA, Rahman S, Yilbas BS, Sahin AZ. Heat transfer enhancement of phase change materials for thermal energy storage applications: A critical review. Renewable and Sustainable Energy Reviews. 2017;74:26-50
  2. 2. Zhao CY, Wu ZG. Heat transfer enhancement of high temperature thermal energy storage using metal foams and expanded graphite. Solar Energy Materials and Solar Cells. 2011;95:636-643
  3. 3. Siaphush A, O’Brien J, Crepeau J. Phase change heat transfer enhancement using copper porous foam. Journal of Heat Transfer. 2008;130:082301
  4. 4. Yin H, Gao X, Ding J, Zhang Z. Experimental research on heat transfer mechanism of heat sink with composite phase change materials. Energy Conversion and Management. 2008;49:1740-1746
  5. 5. Gao L, Zhao J, An Q, Zhao D, Meng F, Liu X. Experiments on thermal performance of erythritol/expanded graphite in a direct contact thermal energy storage container. Applied Thermal Engineering. 2017;113:858-866
  6. 6. Lafdi K, Mesalhy O, Shaikh S. Experimental study on the influence of foam porosity and pore size on the melting of phase change materials. Journal of Applied Physics. 2007;102:083549
  7. 7. Zhong Y, Guo Q, Li S, Shi J, Liu L. Heat transfer enhancement of paraffin wax using graphite foam for thermal energy storage. Solar Energy Materaisl and Solar Cells. 2010;94:1011-1014
  8. 8. Wu ZG, Zhao CY. Experimental investigations of porous materials in high temperature thermal energy storage systems. Solar Energy. 2011;85:1371-1380
  9. 9. Zhang P, Meng ZN, Zhu H, Wang YL, Peng SP. Melting heat transfer characteristics of a composite phase change material fabricated by paraffin and metal foam. Applied Energy. 2017;185:1971-1983
  10. 10. Gulfam R, Zhu W, Li X, Cheema II, Sheng P, Zhao G, Deng Y. Design, fabrication and numerical analysis of compact thermal management system integrated with composite phase change material and thermal bridge. Energy Conversion and Management. 2018;156:25-33
  11. 11. Oya T, Nomura T, Okinaka N, Akiyama T. Phase change composite based on porous nickel and erythritol. Applied Thermal Engineering. 2012;40:373-377
  12. 12. Nomura T, Zhu C, Nan S, Tabuchi K, Wang S, Akiyama T. High thermal conductivity phase change composite with a metal-stabilized carbon-fiber network. Applied Energy. 2016;179:1-6
  13. 13. Nomura T, Okinaka N, Akiyama T. Impregnation of porous material with phase change material for thermal energy storage. Materials Chemistry and Physics. 2009;115:846-850
  14. 14. Li Y, Guo B, Huang G, Kubo S, Shu P. Characterization and thermal performance of nitrate mixture/SiC ceramic honeycomb composite phase change materials for thermal energy storage. Applied Thermal Engineering. 2015;81:193-197
  15. 15. Fan L, Zhu Z, Liu M. A similarity solution to unidirectional solidification of nano-enhanced phase change materials (NePCM) considering the mushy region effect. International Journal of Heat and Mass Transfer. 2015;86:478-481
  16. 16. Qi G, Yang J, Bao R, Liu Z, Yang W, Xie B, Yang M. Enhanced comprehensive performance of polyethylene glycol based phase change material with hybrid graphene nanomaterials for thermal energy storage. Carbon. 2015;88:196-205
  17. 17. Tao YB, Lin CH, He YL. Preparation and thermal properties characterization of carbonate salt/carbon nanomaterial composite phase change material. Energy Conversion and Management. 2015;97:103-110
  18. 18. Kim S, Drzal LT. High latent heat storage and high thermal conductive phase change materials using exfoliated expanded nanoplatelets. Solar Energy Materials and Solar Cells. 2009;93:136-142
  19. 19. Shi J, Ger M, Liu Y, Fan Y, Wen N, Li C, Nenwen P. Improving the thermal conductivity and shape-stabilization of phase change materials using nanographite additives. Carbon. 2013;51:365-372
  20. 20. Zhang H, Rindt CCM, Smeulders DMJ, Nedea SV. Nanoscale heat transfer in carbon nanotubes - sugar alcohol composite as heat storage materials. Journal of Physical Chemistry C. 2016;120:21915-21924
  21. 21. Samimi F, Babapoor A, Azizi M, Karimi G. Thermal management analysis of a Li-ion battery cell using phase change material loaded with carbon fibers. Energy. 2016;96:355-371
  22. 22. Nomura T, Tabuchi K, Zhu C, Sheng N, Wang S, Akiyama T. High thermal conductivity phase change composite with percolating carbon fiber network. Applied Energy. 2015;154:678-685
  23. 23. Kibria MA, Anisur MR, Mahfuz MH, Saidur R, Metselaar IHSC. A review on thermophysical properties of nanoparticle dispersed phase change materials. Energy Conversion and Management. 2015;95:69-89
  24. 24. Zhang XL, Chen XD, Zhao QZ, Ding L. The research on the dispersion effect improvement for nano-copper in erythritol. Materials Research Innovations. 2015;19(S1):9-13
  25. 25. Zhichao L, Qiang Z, Gaohui W. Preparation and enhanced heat capacity of nano-titania doped erythritol as phase change material. International Journal of Heat and Mass Transfer. 2015;80:653-659
  26. 26. Mettawee ES, Assassa GMR. Thermal conductivity enhancement in a latent heat storage system. Solar Energy. 2007;81:839-845
  27. 27. Motahar S, Nikkam N, Alemrajabi AA, Khodabandeh R, Toprak MS, Muhammed M. Experimental investigation on thermal and rheological properties of n-octadecane with dispersed TiO2 nanoparticles. International Communications in Heat and Mass Transfer. 2014;59:68-74
  28. 28. Wang J, Xie H, Xin Z, Li Y, Chen L. Enhancing thermal conductivity of palmitic acid based phase change materials with carbon nanotubes as fillers. Solar Energy. 2010;84:339-344
  29. 29. Zeng JL, Cao Z, Yang DW, Xu F, Sun LX, Zhang XF, Zhang L. Effects of MWNTs on phase change enthalpy and thermal conductivity of a solid-liquid organic PCM. Journal of Thermal Analysis and Calorimetry. 2009;95:507-512
  30. 30. Khyad A, Samrani H, Bargach MN, Tadili R. Energy storage with PCMs: Experimental analysis of paraffin’s phase change phenomenon & improvement of its properties. Journal of Material & Environment Science. 2016;7:2551-2560
  31. 31. Sparrow EM, Larson ED, Ramsey JW. Freezing on a finned tubed for either conduction-controlled or natural-convection-controlled heat transfer. International Journal of Heat and Mass Transfer. 1981;24:273-284
  32. 32. Tao YB, He YL, Qu ZG. Numerical study on performance of molten salt phase change thermal energy storage system with enhanced tubes. Solar Energy. 2012;86:1155-1163
  33. 33. Velraj R, Seeniraj RV, Hafner B, Faber C, Schwarzer K. Experimental analysis and numerical modelling of inward solidification on a finned vertical tube for a latent heat storage unit. Solar Energy. 1997;60:281-290
  34. 34. Zhao D, Tan G. Numerical analysis of a shell-and-tube latent heat storage unit with fins for air-conditioning application. Applied Energy. 2015;138:381-392
  35. 35. Erek A, Ilken Z, Acar MA. Experimental and numerical investigation of thermal energy storage with a finned tube. International Journal of Energy Research. 2005;29:283-301
  36. 36. Liu Z, Sun X, Ma C. Experimental investigations on the characteristics of melting processes of stearic acid in an annulus and its thermal conductivity enhancement by fins. Energy Conversion and Management. 2005;46:959-969
  37. 37. Hosseini MJ, Rahimi M, Bahrampoury R. Thermal analysis of PCM containing heat exchanger enhanced with normal annular fines. Mechanical Science. 2015;6:221-234
  38. 38. Gharebaghi M, Sezai I. Enhancement of heat transfer in latent heat storage modules with internal fins. Numerical Heat Transfer, Part A: Application. 2007;53:749-765
  39. 39. Sharifi N, Bergman TL, Faghri A. Enhancement of PCM melting in enclosures with horizontally-finned internal surfaces. International Journal of Heat and Mass Transfer. 2011;54:4182-4192
  40. 40. Mahmoud S, Tang A, Toh C, Al-dadah R, Soo SL. Experimental investigation of inserts configurations and PCM type on the thermal performance of PCM based heat sinks. Applied Energy. 2013;112:1349-1356
  41. 41. Arshad A, Ali HM, Ali M, Manzoor S. Thermal performance of phase change material (PCM) based pin-finned heat sinks for electronics devices: Effect of pin thickness and PCM volume fraction. Applied Thermal Engineering. 2017;112:143-155
  42. 42. Rahimi M, Ranjbar AA, Ganji DD, Sedighi K, Hosseini MJ. Experimental investigation of phase change inside a finned-tube heat exchanger. Journal of Engineering, 2014. 641954
  43. 43. Fang M, Chen G. Effects of different multiple PCMs on the performance of a latent thermal energy storage system. Applied Thermal Engineering. 2007;27:994-1000
  44. 44. Wang J, Chen G, Jiang H. Theoretical study on a novel phase change process. International Journal of Energy Research. 1999;23:287-294
  45. 45. Cui H, Yuan X, Hou X. Thermal performance analysis for a heat receiver using multiple phase change materials. Applied Thermal Engineering. 2003;23:2353-2361
  46. 46. Zhipei H, Li A, Gao R, Yin H. Enhanced heat transfer for PCM melting in the frustum-shaped unit with multiple PCMs. Journal of Thermal Analysis and Calorimetry. 2015;120:1407-1416
  47. 47. Kurnia JC, Sasmito AP, Jangam SV, Mujumdar AS. Improved design for heat transfer performance of a novel phase change material (PCM) thermal energy storage (TES). Applied Thermal Engineering. 2013;50:896-907
  48. 48. Tao YB, He YL, Liu YK, Tao WQ. Performance optimization of two-stage latent heat storage unit based on entransy theory. International Journal of Heat and Mass Transfer. 2014;77:695-703
  49. 49. Xu HJ, Zhao CY. Thermodynamic analysis and optimization of cascaded latent heat storage system for energy efficient utilization. Energy. 2015;90:1662-1673
  50. 50. Wang C, Zhu Y. Optimization of double-stage latent heat storage unit in whole cycle with entransy analysis. International Journal of Heat and Mass Transfer. 2017;114:1013-1024
  51. 51. Wang C, Zhu Y. Entransy analysis on boiler air pre-heater with multi-stage LHS unit. Applied Thermal Engineering. 2018;130:1139-1146
  52. 52. Wang C, Zhu Y. Entransy analysis on optimization of a double-stage latent heat storage unit with the consideration of an unequal separation. Energy. 2018;148:386-396
  53. 53. Jacob R, Bruno F. Review on shell materials used in the encapsulation of phase change materials for high temperature thermal energy storage. Renewable and Sustainable Energy Reviews. 2015;48:79-87
  54. 54. Liu C, Rao Z, Zhao J, Huo Y, Li Y. Review on nanoencapsulated phase change materials: Preparation, characterization and heat transfer enhancement. Nano Energy. 2015;13:814-826
  55. 55. Gunasekara SN, Martin V, Chiu JN. Phase equilibrium in the design of phase change materials for thermal energy storage: State-of-the-art. Renewable and Sustainable Energy Reviews. 2017;73:558-581
  56. 56. Liu S, Li Y, Zhang Y. Review on heat transfer mechanisms and characteristics in encapsulated PCMs. Heat Transfer Engineering. 2015;36:880-901
  57. 57. Jamekhorshid A, Sadrameli SM, Farid M. A review of microencapsulation methods of phase change materials (PCMs) as a thermal energy storage (TES) medium. Renewable and Sustainable Energy Reviews. 2014;31:531-542
  58. 58. Weiguang S, Jo D, Kokogiannakis G. Review of solid-liquid phase change materials and their encapsulation technologies. Renewable and Sustainable Energy Reviews. 2015;48:373-391
  59. 59. Milian YE, Gutierrez A, Grageda M, Ushak S. A review on encapsulation techniques for inorganic phase change materials and the influence on their thermophysical properties. Renewable and Sustainable Energy Reviews. 2017;73:983-999
  60. 60. Ma B, Li J, Zhe X, Peng Z. Fe-shell/Cu-core encapsulated metallic phase change materials prepared by aerodynamic levitation method. Applied Energy. 2014;132:568-574
  61. 61. Chen Z, Fei Y, Zeng X, Zhang Z. Preparation, characterization and thermal properties of nanocapsules containing phase change materials n-dodecanol by miniemulsion polymerization with polymerizable emulsifier. Applied Energy. 2012;91:7-12
  62. 62. Yang W, Yue Z, Baoxing X. A hybrid elastomeric foam-core/solid-shell spherical structure of enhanced energy absorption performance. International Journal of Solids and Structures. 2016;92-93:17-28
  63. 63. Shabgard H, Robak CW, Bergman TL, Faghri A. Heat transfer and exergy analysis of cascaded latent heat storage with gravity-assisted heat pipes for concentrating solar power applications. Solar Energy. 2012;86:816-830
  64. 64. Shabgard H, Faghri A, Bergman TL, Andraka CE. Numerical simulation of heat pipe-assisted latent heat thermal energy storage unit for dish-stirling systems. Journal of Solar Energy Engineering. 2013;136:021025
  65. 65. Robak CW, Bergman TL, Faghri A. Enhancement of latent heat storage using embedded heat pipes. International Journal of Heat and Mass Transfer. 2011;54:3476-3484
  66. 66. Shabgard H, Bergman TL, Sharifi N, Faghri A. High temperature latent heat thermal energy storage using heat pipes. International Journal of Heat and Mass Transfer. 2010;53:2979-2988
  67. 67. Nithyanandam K, Pitchumani R. Computational studies on a latent thermal energy storage system with integral heat pipes for concentrating solar power. Applied Energy. 2013;103:400-415
  68. 68. Nithyanandam K, Pitchumani R. Computational modeling of dynamic response of a latent thermal energy storage system with embedded heat pipes. Journal of Solar Energy Engineering. 2013;136:011010
  69. 69. Jung EG, Boo JH. Thermal analytical model of latent thermal storage with heat pipe exchanger for concentrated solar power. Solar Energy. 2014;102:318-332
  70. 70. Khalifa A, Tan L, Date A, Akbarzadeh A. A numerical and experimental study of solidification around axially finned heat pipes for high temperature latent heat thermal energy storage units. Applied Thermal Engineering. 2014;70:609-619
  71. 71. Nithyanandam K, Pitchumani R. Computational studies on metal foam and heat pipe enhanced latent thermal energy storage. Journal of Heat Transfer. 2014;136:051503
  72. 72. Esen M, Durmus A, Durmus A. Geometric design of solar-aided latent heat store depending on various parameters and phase change materials. Solar Energy. 1998;62:19-28
  73. 73. Xia L, Zhang P, Wang RZ. Numerical heat transfer analysis of the packed bed latent heat storage system based on an effective packed bed model. Energy. 2010;35:2022-2032
  74. 74. Izquierdo-Barrientos MA, Sobrino C, Almendros-Ibanez JA. Modeling and experiments of energy storage in a packed bed with PCM. International Journal of Multiphase Flow. 2016;89:1-9
  75. 75. Modi A, Perez-segarra CD. Thermocline thermal storage systems for concentrated solar power plants: one-dimensional numerical model and comparative analysis. Solar Energy. 2014;100:84-93
  76. 76. Opitz F, Treffinger P. Packed bed thermal energy storage model—Generalized approach and experimental validation. Applied Thermal Engineering. 2014;73:245-252
  77. 77. Amin NAM, Belusko M, Bruno F. An effectiveness-NTU model of a packed bed PCM thermal storage system. Applied Energy. 2014;134:356-362
  78. 78. Karthikeyan S, Velraj R. Numerical investigation of packed bed storage unit filled with PCM encapsulated spherical containers—A comparison between various mathematical models. International Journal of Thermal Sciences. 2012;60:153-160
  79. 79. Pakrouh R, Hosseini MJ, Ranjbar AA, Bahrampoury R. A numerical method for PCM-based pin fin heat sinks optimization. Energy Conversion and Management. 2015;103:542-552
  80. 80. Shnueli H, Ziskind G, Letan R. Melting in a vertical cylindrical tube: Numerical investigation and comparison with experiments. International Journal of Heat and Mass Transfer. 2010;53:4082-4091
  81. 81. Cascetta M, Cau G, Puddu P, Serra F. A comparison between CFD simulation and experimental investigation of a packed-bed thermal energy storage system. Applied Thermal Engineering. 2016;98:1263-1272
  82. 82. Sciacovelli A, Colella F, Verda V. Melting of PCM in a thermal energy storage unit: Numerical investigation and effect of nanoparticle enhancement. International Journal of Energy Research. 2013;37:1610-1623
  83. 83. Safari A, Saidur R, Sulaiman FA, Yan X, Dong J. A review on supercooling of phase change materials in thermal energy storage systems. Renewable and Sustainable Energy Reviews. 2017;70:905-919
  84. 84. Vasu A, Hagos FY, Noor MM, Mamat R, Azmi WH, Abdullah AA, Ibrahim TK. Corrosion effect of phase change materials in solar thermal energy storage application. Renewable and Sustainable Energy Reviews. 2017;76:19-33
  85. 85. Tay NHS, Liu M, Belusko M, Bruno F. Review on transportable phase change material in thermal energy storage systems. Renewable and Sustainable Energy Reviews. 2017;75:264-277
  86. 86. Teng T, Cheng C, Cheng C. Performance assessment of heat storage by phase change materials containing MWCNTs and graphite. Applied Thermal Engineering. 2013;50:637-644
  87. 87. Da HC, Lee J, Hong H, Kang YT. Thermal conductivity and heat transfer performance enhancement of phase change materials (PCM) containing carbon additives for heat storage application. International Journal of Refrigeration. 2014;42:112-120
  88. 88. Dincer I, Rosen MA. Thermal energy storage: systems and applications. 2nd ed. John Wiley & Sons; 2011
  89. 89. Heckenkamp J. Baumann H. Sonderdruck aus Nachrichten: Latentwaermespeicher; 1997
  90. 90. Sari A, Kaygusuz K. Thermal performance of myristic acid as a phase change material for energy storage application. Renewable Energy. 2001;24:303-317
  91. 91. Lane GA. Low temperature heat storage with phase change materials. International Journal of Ambient Energy. 1980;1:155-168

Written By

Cheng Wang and Ye Zhu

Submitted: October 30th, 2017 Reviewed: March 27th, 2018 Published: November 5th, 2018