Thermodynamic insights into n-alkanes phase change materials for thermal energy storage

Huimin Yn,Huning Yng,Jipeng Luo,Nn Yin,Zhicheng Tn,Qun Shi,∗

a Thermochemistry Laboratory,Liaoning Province Key Laboratory of Thermochemistry for Energy and Materials,Dalian National Laboratory for Clean Energy,Dalian Institute of Chemical Physics,Chinese Academy of Sciences,Dalian 116023,China

b Beijing Institute of Spacecraft System Engineering,Beijing 100094,China

c University of Chinese Academy of Sciences,Beijing 100049,China

Keywords:n-Alkanes Phase change materials Thermal energy storage Heat capacity Thermodynamic functions PPMS Adiabatic calorimetry

ABSTRACT n-Alkanes have been widely used as phase change materials (PCMs) for thermal energy storage applications because of their exceptional phase transition performance,high chemical stability,long term cyclic stability and non-toxicity.However,the thermodynamic properties,especially heat capacity,of n-alkanes have rarely been comprehensively investigated in a wide temperature range,which would be insufficient for design and utilization of n-alkanes-based thermal energy storage techniques.In this study,the thermal properties of n-alkanes (C18H38-C22H46),such as thermal stability,thermal conductivity,phase transition temperature and enthalpy were systematically studied by different thermal analysis and calorimetry methods,and compared with previous results.Thermodynamic property of these n-alkanes was studied in a wide temperature range from 1.9 K to 370 K using a combined relaxation (Physical Property Measurement System,PPMS),differential scanning and adiabatic calorimetry method,and the corresponding thermodynamic functions,such as entropy and enthalpy,were calculated based on the heat capacity curve fitting.Most importantly,the heat capacities and related thermodynamic functions of n-heneicosane and n-docosane were reported for the first time in this work,as far as we know.This research work would provide accurate and reliable thermodynamic properties for further study of n-alkanes-based PCMs for thermal energy storage applications.

The development and application of energy storage materials have attracted wide attention,due to the contradiction between energy shortage and increasing energy demand caused by the massive consumption of fossil energy.As a commonly used thermal energy storage material,phase change materials (PCMs) can store and release a large amount of latent heat in the phase change process and keep almost constant temperature,which provides favorable conditions for thermal energy recovery and secondary utilization,as well as temperature control.In addition to relatively high energy storage density,PCMs are inexpensive,non-toxic and readily available [1],making them used in solar energy storage,electronic chip cooling,battery thermal management,industrial waste heat recovery,energy saving building,establishment of thermal energy storage (TES) system,intelligent temperature control clothing and other fields [2–7].PCMs are generally divided into inorganic phase change materials and organic phase change materials.Although inorganic PCMs have higher enthalpy of phase change,they also have obvious shortcomings,such as undercooling,corrosion,phase separation and poor thermal stability [8].Compared with inorganic PCMs,organic PCMs are always better choice of PCMs in TES system because of their non-corrosion,almost nonundercooling and self-nucleation [4,9].n-Alkanes,one of the most important organic PCMs [1],are unbranched saturated open-chain hydrocarbons,which not only have the advantages of organic PCMs mentioned above,but also have stable chemical properties,good cycle stability and low vapor pressure.Although the current study onn-alkanes is mainly focused on constructing composite PCMs for preventing liquid phase leakage and how to improve the thermal conductivity of PCMs [10–19],it is necessary to understand the basic physical and chemical properties of pure substances for practical applications ofn-alkanes in thermal energy storage field.The thermodynamic properties of puren-alkanes,such as heat capacity,enthalpy and entropy,and other thermal properties of thermal stability and thermal conductivity,are the prerequisites for the reasonable and effective preparation of PCMs that meet the needs of practical production,and the design of TES system also needs to understand the phase transition behavior ofn-alkanes and the thermophysical parameters of condensed matter.

Vélezet al.[20,21]studied the thermal properties of somenalkanes from C15H32to C20H42,and measured their thermal conductivities and heat capacities in the temperature range of 258–350 K by the transient hot wire method and differential scanning (DSC) method,respectively.The reliable liquid phase heat capacity data were obtained in their work,but the heat capacity data in the solid phase region deviated from the literature values of more than 8%,as mentioned in their paper.Duruptet al.[22]measured the liquid phase heat capacities ofn-octadecane,n-nonadecane andn-docosane in the temperature range of 373–473 K using a DSC method.However,it was found that the DSC instrument used in their work had a systematic deviation of 4%between the two sets of data,and only 50 determinations were presented from the heat capacity data of each sample.Fadenet al.[23]reviewed the thermophysical property data ofn-octadecane.According to their summary of heat capacity data from the literature,most researchers used DSC to only measure the heat capacity of liquid phasen-octadecane,and also relatively few observations were presented.For the heat capacity of solidn-octadecane,except that there were about 50 heat capacity data provided by Vélezet al.[20]and Messerlyet al.[24],the solid phase heat capacity data in other literature were not comprehensive.Generally,regarding the research on heat capacity ofn-alkanes,most of the reports so far are mainly focused on the measurement of liquidnalkanes,while most of solid phase heat capacity measurements are reported several decades ago and no related heat capacity studies are performed currently.On the other hand,these reported heat capacity studies onn-alkanes are generally performed using the DSC method,the accuracy of which may be insufficient for further investigate their thermodynamic properties.

In our previous research work,we reported heat capacity studies on some sugar alcohols and polyethylene glycols (PEG) [25,26],which are two of hottest PCMs currently concerned for new type PCMs synthesis and phase transition performance investigation.However,rather surprisingly,the heat capacity and related thermodynamic properties ofn-alkanes,commonly used PCMs in thermal energy storage application field,have not been systematically studied yet until now.In this work,we have studied thermodynamic properties especially heat capacity of fiven-alkanes ofn-octadecane,n-nonadecane,n-eicosane,n-heneicosane andndocosane using various thermal analysis and calorimetry methods.The crystal structure,thermal stability and thermal conductivity of thesen-alkanes have been determined using X-ray diffractometer (XRD),thermogravimetric analyzer (TGA) and thermal constants analyser,respectively.The phase transition process ofn-alkanes has been evaluated by the differential scanning calorimeter.The solid and liquid phase heat capacity of then-alkanes have been measured in a wide temperature range from 1.9 K to 370 K using a combined relaxation (physical property measurement aystem,PPMS),DSC and adiabatic calorimetry (AC) method,and the corresponding thermodynamic functions have been consequently obtained by the heat capacity curve fitting.

Then-octadecane,n-nonadecane,n-eicosane andn-docosane samples were commercially obtained from Sigma Aldrich,and then-heneicosane sample was from Supelco.The corresponding information of these samples is listed in Table 1.The detailed information about the sample characterization experiment and instrument correction is provided in Supporting information [27–31].

Fig.1.TG curves of n-alkanes samples.

n-Alkane (CnH2n+2) is a special kind of organic matter,and its low-temperature crystal may have different crystal structure because of the different number of carbon atoms:triclinic,12 ≤n(even) ≤26;monoclinic,2836 andn(odd) [32].The X-ray pattern of these fiven-alkanes are shown in Fig.S1 (Supporting information).Forn-octadecane,n-eicosane andn-docosane of evenn-alkanes,the XRD patterns presented a series of characteristic diffraction peaks of (010),(011),(100) and (111) reflections of the triclinic crystal [33–35].The XRD patterns ofn-nonadecane andn-heneicosane revealed two main diffraction peaks assigned to the (110) and (020) planes of orthorhombic crystal [36–38].The difference in crystal structure between odd-numbered and even-numberedn-alkanes is due to the different ways of terminal methyl packing.The end methyl groups of oddn-alkanes pack in the same direction,while the evenn-alkanes protrude in opposite directions.For short chains,the terminal methyl-methyl interactions occupy the dominant position of solid-state packing which makes the formation of orthogonal cell difficult and leads to the formation of a unimolecular triclinic lattice.With the increase of chain length,the influence of parallel chain packing interactions is higher than that of interlamellar packing.For odd-numberednalkanes,the mirror symmetry of molecules tends to prevent the best close packing and form orthogonal crystals [39,40].As the temperature increases,somen-alkanes still have a rotator state between the crystal and the liquid state,which is also called a disordered state.The transition from crystal to rotator state is called order-disorder transition (o-d transition).Among the fiven-alkanes studied,n-nonadecane andn-heneicosane andn-docosane all exist o-d transition,whilen-octadecane andn-eicosane only have a crystal-liquid phase transition process [41].The o-d transition may be related to the thermal motion of particles.When the o-d transition temperature (To-d) is reached,the molecules in the crystal begin to vibrate randomly at the equilibrium position.Then with the increase of temperature,the thermal motion becomes more intense and the molecular chains behave thermal torsional motion around their axes,losing the fixed orientation in the structure [42].Until the temperature reaches the melting point,n-alkanes tends to lose their solid structure and begin to transform into the liquid phase.

The thermal stability evaluation ofn-alkanes was performed by the TGA technique.It can be seen from the measured curves in Fig.1 thatn-alkanes of C18-C22began to decompose in the temperature range of 373–413 K.And as the number of carbon atoms increases,the initial and end decomposition temperature ofn-alkanes of C18-C22increases.The weight loss percentage ofn-octadecane,n-nonadecane,n-eicosane,n-heneicosane andn-docosane was 94.76%,95.16%,94.66%,94.95% and 95.33%,respectively.According to these results,it can be determined thatn-alkanes are thermally stable in the temperature region below 373 K.

The thermal conductivity ofn-alkanes was measured by the transient plane source (TPS) technique,and the obtained data are listed in Table S1 (Supporting information).The thermal conductivity ofn-octadecane,n-eicosane andn-docosane at room temperature is close to each other and significantly larger than that ofnnonadecane andn-heneicosane.After measuring the thermal conductivity ofn-undecane,n-dodecane andn-tridecane by transient hot wire method,Hans Forsmanet al.[43]found that the thermal conductivity of even-numberedn-dodecane was higher than that of adjacentn-undecane andn-dodecane with odd number of carbon atoms.In this study,the same conclusion is reached,that is,the thermal conductivity of even-numberedn-alkanes is higher than that of adjacent odd-numberedn-alkanes.The previous results for thesen-alkanes are shown in Table 2 [6,45–50],and it can be seen that our measured thermal conductivities ofnoctadecane,n-nonadecane andn-eicosane in the solid phase are in good agreement with these results.However,the thermal conductivity of crystallinen-heneicosane is rarely reported.In the term ofn-docosane,the thermal conductivity values obtained by different measurement techniques are quite different in a relatively close temperature range.Felipeet al.[44]obtained that the thermal conductivity ofn-docosane was 0.668 W m-1K-1by a C-Therm Technologies TCi Analyzer based on a modified transient plane source(MTPS) method,while Sari and Karaipekli [45]measured the thermal conductivity ofn-docosane to be 0.22 W m-1K-1by transient hot-wire method.In this work,the thermal conductivity ofn-docosane was determined to be 0.4091 W m-1K-1by a thermal constants analyser (Hot Disk TPS2500S) based on the TPS method.

Table 1 Sample information in this work.a

Table 2 Comparison of thermal conductivity data of pure n-alkanes in the solid (crystalline) state.

The DSC curve can reflect the phase transition behavior of the materials,including phase transition enthalpy and phase transition temperature.Differentn-alkanes have different DSC curve characteristics due to the existence of the rotator phase.For oddnumberedn-alkanes ofn-nonadecane andn-heneicosane,there are two endothermic peaks in the DSC curve of the heating run at a rate of 10 K/min corresponding respectively to the o-d transition and solid-liquid transition as shown in Fig.2a.For evennumberedn-alkanes ofn-octadecane,n-eicosane andn-docosane,only one obvious endothermic peak can be obtained in the DSC curve of the heating run at a rate of 10 K/min,moreover,there is a shoulder in the DSC curve ofn-docosane as shown in Fig.2b.It should be noted that the endothermic peak represents the melting peak ofn-octadecane andn-eicosane while the endothermic peak includes both the o-d transition and solid-liquid transition ofndocosane.Therefore,the calculated enthalpy value ofn-docosane,which agreed well with the sum of phase transition enthalpy from literature,derived from the DSC curve was the total enthalpy value including the o-d transition enthalpy and the melting enthalpy.n-Docosane undergoes the o-d transition from the triclinic crystal to the rotator phase,and then the solid-liquid transition from the rotator phase to the liquid phase.Wanget al.[51]obtained a perfect DSC curve at a heating rate of 0.5 mK/s using a high resolution and super-sensitive DSC in which there were conspicuous o-d transition peak and melting peak.It was obvious from the curve that the onset/end temperature of the o-d transition and the onset melting temperature of the solid-liquid transition was 316.22 K,316.86 K and 317.25 K,respectively.According to the previous work of Wanget al.,n-docosane exists in the rotator phase only within the temperature range of less than 0.5 K,which means the heating rate of 10 K/min may be too fast to obtain the independent o-d transition peak and melting peak forn-docosane.The phase transition properties ofn-octadecane,n-eicosane andn-docosane were determined again with a heating rate of 0.5 mK/s and the DSC curves are shown in Figs.2c and d.It can be seen that there was only one peak forn-octadecane andn-eicosane while two peaks were presented forn-docosane at the slow heating rate of 0.5 mK/s,and these DSC results are consistent with those reported in the literature [41,51].

For the AC experiment ofn-octadecane,the enthalpy of fusion is calculated by subtracting the contribution of the heat capacity of the empty container,and the hypothetic normal heat capacity of crystal and liquid from the total energy input by the heater during the phase change in the temperature range of (T1toT2) (Eq.1) [52],

whereQTis the total energy input during the fusion region,QEis the heat absorbed by the empty container,QCandQLare the heat absorbed by crystal and liquid which are calculated by extending the heat capacity baselines derived from heat capacities outside of the phase transition region in the temperature region of (T1toTm)and (TmtoT2),respectively,Tmis the melting point,andnis the molar substance amount of a sample.The detailed continuous energy input during the phase transition,the heat capacity contribution from the empty container,crystal and liquid ofn-octadecane and the calculated fusion enthalpy were listed in Table S2 (Supporting information).Meanwhile,the triple point temperature can be determined by AC fractional melting experiments based on the Van’t Hoff equation (Eq.2) [53,54],

Fig.2.DSC curves of (a) the odd-numbered n-alkanes at heating rate of 10 K/min;(b) the even-numbered n-alkanes at heating rate of 10 K/min;(c) n-octadecane and n-eicosane at heating rate of 0.5 mK/s;(d) n-docosane at heating rate of 0.5 mK/s.

Fig.3.Melting curve (Ti vs. 1/F) of n-octadecane.

whereT0is the melting point of purity sample;Tiis the temperature during fusion of the measured sample;Ris the gas constant;ΔfusHmis the molar enthalpy of fusion;xis the molar fraction of impurity;Fis the fraction melted,which means the ratio of the amount of melted sample to the total amount of tested sample.Fcan also be expressed as Eq.3,

whereqiis the heat absorbed by the sample during the melting process at the temperatureTi.

For a certain tested sample,the values ofT0,R,ΔfusHmandxare constant,so thatTiis proportional to 1/F,and when 1/Fequals 1 and 0,Tmeans the melting point of sample and purity substance,respectively.The results of AC fractional melting experiment are listed in Table S3 (Supporting information) and plotted in Fig.3.The coefficient of determinationR2is 0.9982,which indicated the good linearity of the melting curve.

Fig.4.Plot of the experimental heat capacities of n-alkanes measured with AC and PPMS.

The detailed data of phase transition properties ofn-alkanes from AC experiment and DSC experiments are listed in Table S4(Supporting information) and the results of comparison with the literature are listed in Table 3 [24,51,55-67].It is concluded that the o-d transition temperature and enthalpy increase with the increasing number of carbon atoms.However,the enthalpy of solid-liquid phase change shows a different changing law as the solid-liquid transition temperature rises with the increase of relative molecular weight.Same as the thermal conductivity rules ofn-alkanes,the melting enthalpy of even-numberedn-alkanes is larger than that of adjacent odd-numberedn-alkanes.Additionally,the average total phase transition enthalpy ofn-docosane measured at the heating rate of 10 K/min was 78.44 kJ/mol,which agreed well with the sum of o-d transition enthalpy and solid-liquid transition enthalpy from different literature [51,55–58].The average total phase transition enthalpy ofn-docosane measured at the heating rate of 0.5 mK/s was 73.47 kJ/mol,which was a little lower than the value of literature.Although then-docosane DSC curves with two separated peaks were obtained in Fig.2d,the phase transition enthalpy obtained was different from those in the literature likely due to the measurement accuracy of the instrument.Wanget al.[51,57]performed the phase change properties with a high resolution and super-sensitive DSC constructed by their research group,which could measure at a very slow rate scan rate like 10 μK/s with a baseline stability of ±3 nW.For our DSC with the slowest scan rate of 0.01 K/min,the same heating rate of 0.5 mK/s as Wanget al.is relatively difficult to control because the temperature dropped on the baseline during the heating measurement.Moreover,the baseline fluctuation was nearly 6 μW,much larger than the baseline stability of the DSC used by Wang.Based on the above two reasons,there may be a difference between the measured phase transition enthalpy ofn-docosane and the literature value.

The molar heat capacities ofn-alkanes measured by AC and PPMS are listed in Tables S5 and S6 (Supporting information),respectively,and plotted against the temperature from 0 to 300 K in Fig.4.It can be seen that at the temperature below 300 K,the heat capacities of fiven-alkanes increase with the temperature increasing,and no obvious thermal anomaly phenomenon or phase transition is observed.The heat capacity ofn-octadecane at 227.70 K deviated a little from the normal heat capacity increase trend.The reason may be caused by the measurement error,and consequently this abnormal value is not used in the fitting heat capacity curve and the calculation of the corresponding thermodynamic functions.At the temperature of about 260 K,the heat capacity curve ofnnonadecane crosses with that ofn-eicosane.In the temperature range of 260–300 K,the heat capacity ofn-nonadecane is larger than that ofn-eicosane,which is consistent with the experimental results of Van Miltenburget al.[68].This cross phenomenon may be due to the entry ofn-nonadecane into the pre-melting temperature region.

Table 3 Comparison of transition properties of n-alkanes between this work and the previous results.

Fig.5.Plot of the experimental heat capacities of n-alkanes from AC,PPMS and DSC measurements.

In order to study the thermal properties and calculate the thermodynamic functions above 300 K,the heat capacities ofn-alkanes were measured using DSC in the temperature range of 190–370 K.The heat capacity results from DSC,PPMS and AC are shown in Fig.5.It can be seen from the heat capacity curves that the molar heat capacities ofn-alkanes increase with the increasing molecular weight ofn-alkanes in the liquid phase.Moreover,the phase transition behavior such as o-d transition and solid-liquid transition in Fig.5 corresponds well to the heat flow curves in the DSC measurements.

For the calculation of thermodynamic functions,a set of theoretical and empirical models were used to fit the experimental heat capacity data.In the low temperature region (T <10 K),a harmonic lattice model was used to fit the low temperature heat capacities (Eq.4) [25,69],

where the odd power terms in temperature represent the contribution from the lattice vibration,and this model may extrapolate the heat capacity to 0 K to calculate the third law entropy [26].In the temperature region above 10 K without phase transition,the heat capacities of then-alkanes were fitted to the following orthogonal polynomial function,while the data in the phase transition region were fitted using a spline fitting method (Eq.5) [25],

The fitting parameters,temperature regions and the corresponding percent mean-square fitting deviations (%RMS) are listed in Tables S7 and S8 (Supporting information).All curve fitting %RMS are much smaller than the heat capacity uncertainties of ±3% fromT=1.9–20 K and ±1% fromT=20–300 K,which indicates that these fits can represent the heat capacity data well.Based on these heat capacity fitting parameters,the thermodynamic functions ofn-alkanes have been derived using the thermodynamic Eqs.6 and 7 in the temperature range from 0 K to 370 K,and the calculated results are presented in Table S9 (Supporting information).Sincen-octadecane andnhexadecane begin to enter the phase transition region at near 298.15 K and 0.1 MPa,the standard molar heat capacity,enthalpy and entropy under this condition are extremely large due to the latent heat of phase transition.Therefore,the standard molar heat capacity,enthalpy and entropy at 273.15 K and 0.1 MPa were provided for alln-alkanes samples.The standard molar heat capacity,entropy and enthalpy at 273.15 K and 0.1 MPa were determined to beCop,m=(428.066 ± 4.28) J K-1mol-1,Som=(324.412 ± 3.24) J K-1mol-1andHom=(55.8122 ± 0.558)kJ/mol forn-octadecane;Cop,m=(463.925± 4.64) J K-1mol-1,Som=(467.791 ± 4.68) J K-1mol-1andHom=(65.8221 ± 0.658)kJ/mol forn-nonadecane;Cop,m=(469.718 ± 4.70) J K-1mol-1,Som=(480.386 ± 4.80) J K-1mol-1andHom=(67.2148 ± 0.672)kJ/mol forn-eicosane;Cop,m=(488.907 ± 4.89) J K-1mol-1,Som=(500.013 ± 5.00) J K-1mol-1andHom=(70.2356 ± 0.702)kJ/mol forn-heneicosane;Cop,m=(503.581 ± 5.04) J K-1mol-1,Som=(523.975 ± 5.24) J K-1mol-1andHom=(73.3775 ± 0.734)kJ/mol forn-docosane.

Fig.6.Comparison/deviations of measured AC,PPMS and DSC heat capacities with/from those reported in literatures [21,22,24,60,68,70–72].●denotes the reference data obtained from experiments,▲denotes the reference data obtained from fitting heat capacity experiments data to a polynomial function,♦denotes the reference data calculated from a polynomial function obtained by group contribution method.

The data on the heat capacity of previously reported alkanes were compared with that of this study and plotted in Fig.6.For solid-phasen-octadecane,the deviations of our AC data from those of Messerlyet al.[24]and Parkset al.[60]are respectively within-2.5% to 2.2% in the temperature range of 80–280 K and-0.8% to 2.7% in the temperature range of 80–270 K except for the deviation at 230 K;the deviations of our DSC data from those of Messerlyet al.[24]and Parkset al.[60]are respectively within-1.3% to 1.9%in the temperature range of 193–278 K and-1.4% to 2.9% in the temperature range of 190–270 K except for the deviation at 230 K.The reason for the large deviation between the data of AC/DSC and the value of Parkset al.is that the heat capacity data at 230 K from Parkset al.is obviously larger than the overall trend of heat capacity.For liquid-phasen-octadecane,the deviations of our AC data from those of Messerlyet al.[24]and Miltenburg [70]are respectively within 0.1% to 0.8% in the temperature range of 304–388 K and-0.2% to 0.4% in the temperature range of 307–371 K,and the deviations of our DSC data from those of Messerlyet al.[24]and Miltenburg [70]are respectively within-0.3% to 0.2% in the temperature range of 318–379 K and-0.5% to 0.7% in the temperature range of 323–371 K.

For solid-phasen-nonadecane,the deviations of our PPMS data from those of Miltenburget al.[68]and Vélezet al.[21]are respectively within ±2% in the temperature range of 14–281 K and-0.2% to 3.3% in the temperature range of 265–280 K,especially the deviations of our PPMS data from those of Miltenburget al.[68]in the temperature range of 20–260 K are within ±0.5% the deviations of our DSC data from those of Miltenburget al.[68]and Vélezet al.[21]are respectively within-0.5% to 1.1% in the temperature range of 191–281 K and-3.1% to 4.0% in the temperature range of 265–280 K.For liquid-phasen-nonadecane,the deviations of our DSC data from those of Miltenburget al.[68]and Vélezet al.[21]are respectively within-0.2% to 1.1% in the temperature range of 321–372 K and-0.5% to 0.4% in the temperature range of 316–334 K.The deviation of our DSC data from that of Grigor’ev and Andolenko [71]is 3.5% at 313.15 K,due to the sample still in phase transition at that temperature.Moreover,Atkinsonet al.[72]measured a series ofn-alkanes includen-nonadecane andndocosane by AC and then obtained liquidn-alkanes heat capacity fitting equations.The deviations of our DSC data from the fitting equation ofn-nonadecane are within-0.9% to 0.4% in the temperature range of 321–372 K.

For solid-phasen-eicosane,the deviations of our PPMS and DSC data from those of Miltenburget al.[68]are respectively within-1.2% to 1.1% in the temperature range of 16–281 K and 1.7% to 2.5% in the temperature range of 192–285 K.For liquid-phaseneicosane,the deviations of our DSC data from those of Miltenburget al.[68]are within-0.5% to 0.7% in the temperature range of 326–370 K.

Forn-heneicosane andn-docosane,few experimental data in the range of 0–370 K can be found for direct comparison.For liquid-phasen-heneicosane,the heat capacity at 315.93 K determined by Grigor’ev and Andolenko [72]was 666.4 J K-1mol-1,while the sample was still in the solid-liquid transition process at the same temperature.For liquid-phasen-docosane,the deviations of our DSC data from the fitting equation ofn-docosane [71]are within-2.4% to 0.3% in the temperature range of 330–370 K.Taking previous publications including the work of Atkinsonet al.into account,Miltenburg [70]fitted a heat capacity function applicable to liquidn-alkanes from C5H12to C26H54with a 0.18% average percentage deviation,which was used to calculate the deviation between our DSC data and fitted data.It can be found that the deviations of our DSC data for liquidn-heneicosane andn-docosane from those of Miltenburg [70]are respectively within 0 to 2.3%in the temperature range of 329–370 K and-2.5% to 0.3% in the temperature range of 333–370 K.For bothn-heneicosane andndocosane,the deviation decreases with the temperature rising and then tends to flatten to around-2.5% gradually.In short,except for the pre-melting temperature region and individual temperature points,the heat capacities of these fiven-alkanes have a deviation of no more than 2.5% from the heat capacity data in the literature,and the DSC heat capacity curves are in good agreement with the AC or PPMS heat capacity data.Therefore,it can be pointed out that our heat capacity data ofn-alkanes are accurate,and it is recommended to use our measured data for the related PCMs study ofn-alkanes.

In summary,the solid structure,thermal stability,thermal conductivity,phase transition temperature,phase transition enthalpy and heat capacity ofn-alkanes have been studied using different analysis method,and these newly measured thermal properties have been compared with those reported in literature.In addition,the heat capacities were measured using a combination of AC/PPMS and DSC calorimetric methods in the temperature range of 1.9–370 K and the corresponding thermodynamic functions ofn-alkanes were derived from the heat capacity curve fitting.The standard molar heat capacity,entropy and enthalpy were determined.Most importantly,the heat capacities and related thermodynamic functions ofn-heneicosane andn-docosane have been reported for the first time in this work,as far as we know.This research work would provide accurate and reliable thermal properties for further study ofn-alkanes as PCMs for the material development and related thermal energy storage system design and construction.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Nature Science Foundation of China (No.22003065),Liaoning Provincial Natural Science Foundation of China (No.2019-MS-318),Science and Technology Major Project of Liaoning Province (No.2019JH1/10300002),the Scientific Instrument Developing Project of the Chinese Academy of Sciences (No.YJKYYQ20190046) and Dalian Institute of Chemical Physics (No.DICP I202036).Q.Shi would like to thank Dalian Outstanding Young Scientific Talent Program (No.2019RJ10).

Supplementary materials

Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.cclet.2021.05.017.