工程科学学报 Chinese Journal of Engineering 赤藓糖醇碳纳米管复合相变材料热特性模拟研究 闫晓鑫冯妍卉邱琳张欣欣 Simulation of thermal properties of erythritol/carbon nanotube composite phase change materials YAN Xiao-xin,FENG Yan-hui.QIU Lin.ZHANG Xin-xin 引用本文： 闫晓鑫，冯妍卉，邱琳，张欣欣.赤藓糖醇/碳纳米管复合相变材料热特性模拟研究.工程科学学报，优先发表.d: 10.13374j.issn2095-9389.2021.09.14.002 YAN Xiao-xin,FENG Yan-hui,QIU Lin,ZHANG Xin-xin.Simulation of thermal properties of erythritol/carbon nanotube composite phase change materials[J].Chinese Journal of Engineering,In press.doi:10.13374/j.issn2095-9389.2021.09.14.002 在线阅读View online::htps:ldoi.org10.13374.issn2095-9389.2021.09.14.002 您可能感兴趣的其他文章 Articles you may be interested in 生物质多孔碳基复合相变材料制备及性能 Preparation and properties of biomass porous carbon composite phase change materials 工程科学学报.2020.42(1)：113 https::/doi.org/10.13374.issn2095-9389.2019.08.06.002 具有核壳结构的FS,微米球与碳纳米管原位复合介孔材料的构建及其在锂离子电池中的应用 Mesoporous composite of core-shell FeS,micron spheres with multi-walled CNTs and its application in lithium ion batteries 工程科学学报.2019,41(4：489htps:l1doi.org10.13374.issn2095-9389.2019.04.009 复合定形相变材料的封装及应用研究新进展 Research progress in encapsulation and application of shape-stabilized composite phase-change materials 工程科学学报.2020,42(11：1422htps:/1doi.org10.13374.issn2095-9389.2020.03.26.002 Zr-Cu薄膜生长及力学性能的分子动力学模拟 Molecular dynamic simulations of the growth and mechanical properties of ZrCu films 工程科学学报.2019,41(4：497 https:1doi.org/10.13374/j.issn2095-9389.2019.04.010 多孔基定形复合相变材料传热性能提升研究进展 Thermal conductivity enhancement of porous shape-stabilized composite phase change materials for thermal energy storage applications:a review 工程科学学报.2020,42(1)：26 https:/doi.org/10.13374.issn2095-9389.2019.07.19.001 电子皮肤新型材料与性能研究进展 Research progress on new materials and properties of electronic skin 工程科学学报.2020,42(6：704htps:/oi.org/10.13374.issn2095-9389.2019.07.18.001
赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 闫晓鑫 冯妍卉 邱琳 张欣欣 Simulation of thermal properties of erythritol/carbon nanotube composite phase change materials YAN Xiao-xin, FENG Yan-hui, QIU Lin, ZHANG Xin-xin 引用本文: 闫晓鑫, 冯妍卉, 邱琳, 张欣欣. 赤藓糖醇/碳纳米管复合相变材料热特性模拟研究[J]. 工程科学学报, 优先发表. doi: 10.13374/j.issn2095-9389.2021.09.14.002 YAN Xiao-xin, FENG Yan-hui, QIU Lin, ZHANG Xin-xin. Simulation of thermal properties of erythritol/carbon nanotube composite phase change materials[J]. Chinese Journal of Engineering, In press. doi: 10.13374/j.issn2095-9389.2021.09.14.002 在线阅读 View online: https://doi.org/10.13374/j.issn2095-9389.2021.09.14.002 您可能感兴趣的其他文章 Articles you may be interested in 生物质多孔碳基复合相变材料制备及性能 Preparation and properties of biomass porous carbon composite phase change materials 工程科学学报. 2020, 42(1): 113 https://doi.org/10.13374/j.issn2095-9389.2019.08.06.002 具有核壳结构的FeS2微米球与碳纳米管原位复合介孔材料的构建及其在锂离子电池中的应用 Mesoporous composite of core-shell FeS2 micron spheres with multi-walled CNTs and its application in lithium ion batteries 工程科学学报. 2019, 41(4): 489 https://doi.org/10.13374/j.issn2095-9389.2019.04.009 复合定形相变材料的封装及应用研究新进展 Research progress in encapsulation and application of shape-stabilized composite phase-change materials 工程科学学报. 2020, 42(11): 1422 https://doi.org/10.13374/j.issn2095-9389.2020.03.26.002 Zr-Cu薄膜生长及力学性能的分子动力学模拟 Molecular dynamic simulations of the growth and mechanical properties of Zr—Cu films 工程科学学报. 2019, 41(4): 497 https://doi.org/10.13374/j.issn2095-9389.2019.04.010 多孔基定形复合相变材料传热性能提升研究进展 Thermal conductivity enhancement of porous shape-stabilized composite phase change materials for thermal energy storage applications: a review 工程科学学报. 2020, 42(1): 26 https://doi.org/10.13374/j.issn2095-9389.2019.07.19.001 电子皮肤新型材料与性能研究进展 Research progress on new materials and properties of electronic skin 工程科学学报. 2020, 42(6): 704 https://doi.org/10.13374/j.issn2095-9389.2019.07.18.001
工程科学学报.第44卷，第X期：1-8.2022年X月 Chinese Journal of Engineering,Vol.44,No.X:1-8,X 2022 https://doi.org/10.13374/j.issn2095-9389.2021.09.14.002;http://cje.ustb.edu.cn 赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 闫晓鑫》，冯妍卉12)四，邱琳12)四，张欣欣2 1)北京科技大学能源与环境工程学院，北京1000832)北京科技大学冶金工业节能减排北京市重点实验室，北京100083 ☒通信作者，冯妍卉，E-mail:yhfeng@me.ustb.edu.cn,邱琳，E-mail:qiulin@ustb.edu.cn 摘要在“碳达峰、碳中和”的大背景下，能源结构从一次能源向新能源转变刻不容缓.由于新能源具有间歇性、波动性的 特点，储能技术可以有效解决上述问题而得到了广泛的关注.相变材料作为储能技术的关键，其热导率低的问题亟需解决 赤藓糖醇作为中低温区常用的高焓值相变材料，热导率仅为0.7WK,严重制约了实际应用中的能量利用效率.本文以 赤藓糖醇作为主要研究对象，采用具有超高热导系数的单壁碳纳米管作为导热增强材料，借助分子动力学模拟的方法探究了 碳纳米管长度、质量分数以及分布方式对赤藓糖醇/碳纳米管复合相变材料热导率的影响规律.当碳纳米管轴向长度小于其 声子平均自由程时，复合相变材料热导率随碳纳米管轴向长度增加而增大，同时随碳纳米管质量分数增加而增大，但表现出 显著的各向异性.由于引入赤藓糖醇-碳纳米管界面，复合相变材料径向热导率相比纯赤藓糖醇反而降低.当碳纳米管在赤 藓糖醇中随机分布时，热导率的各向异性得到了显著改善且各方向热导率均得到了提升.通过对比复合前后赤藓糖醇与碳 纳米管的声子振动态密度发现，由于两者间的相互作用.碳纳米管的声子振动受到抑制.而赤藓糖醇中声子热输运得到激发， 从而提高了热导率. 关键词赤藓糖醇；碳纳米管；复合相变材料：热导率；分子动力学模拟 分类号TK02 Simulation of thermal properties of erythritol/carbon nanotube composite phase change materials YAN Xiao-xin),FENG Yan-hui2,QIU Lin2,ZHANG Xin-xin 2) 1)School of Energy and Environmental Engineering,University of Science and Technology Beijing,Beijing 100083,China 2)Beijing Key Laboratory of Energy Conservation and Emission Reduction in Metallurgical Industry,University of Science and Technology Beijing, Beijing 100083,China Corresponding author,FENG Yan-hui,E-mail:yhfeng @me.ustb.edu.cn;QIU Lin,E-mail:qiulin@ustb.edu.cn ABSTRACT In the area of"carbon peaking and carbon neutralization,"changing energy structure from primary energy to new energy is an extremely important issue.Due to the intermittent and fluctuating characteristics of new energy,energy storage technology has proven a viable solution to this issue thus has attracted extensive attention.As a key to energy storage technology,the problem of the low thermal conductivity of phase change materials(PCMs)requires immediate attention.Erythritol is a high enthalpy phase change material commonly used in low-to-medium temperature processes.Its thermal conductivity of only 0.7 W.m-K-seriously hinders its energy utilization efficiency in practical application.In this paper,erythritol is the main research focus,and single-walled carbon nanotubes (CNTs)of ultra-high thermal conductivity are used as thermal conductivity reinforcements.The effects of length,mass fraction,and the distribution of CNTs on the thermal conductivity of erythritol/CNT composite PCMs were studied by means of molecular dynamics 收稿日期：2021-09-14 基金项目：国家重点研发计则资助项目(2018YFA0702300):国家自然科学基金资助项目(51876007,51876008)：北京市科技新星计划资助 项目(Z201100006820065):北京市自然科学基金资助项目(3202020)：北京科技大学青年教师学科交叉研究资助项目（中央高校基本科研 业务费专项资金)(FRF-IDRY-19-004)
赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 闫晓鑫1)，冯妍卉1,2) 苣，邱 琳1,2) 苣，张欣欣1,2) 1) 北京科技大学能源与环境工程学院，北京 100083 2) 北京科技大学冶金工业节能减排北京市重点实验室，北京 100083 苣通信作者， 冯妍卉，E-mail: yhfeng@me.ustb.edu.cn; 邱琳，E-mail: qiulin@ustb.edu.cn 摘 要 在“碳达峰、碳中和”的大背景下，能源结构从一次能源向新能源转变刻不容缓. 由于新能源具有间歇性、波动性的 特点，储能技术可以有效解决上述问题而得到了广泛的关注. 相变材料作为储能技术的关键，其热导率低的问题亟需解决. 赤藓糖醇作为中低温区常用的高焓值相变材料，热导率仅为 0.7 W·m–1·K–1，严重制约了实际应用中的能量利用效率. 本文以 赤藓糖醇作为主要研究对象，采用具有超高热导系数的单壁碳纳米管作为导热增强材料，借助分子动力学模拟的方法探究了 碳纳米管长度、质量分数以及分布方式对赤藓糖醇/碳纳米管复合相变材料热导率的影响规律. 当碳纳米管轴向长度小于其 声子平均自由程时，复合相变材料热导率随碳纳米管轴向长度增加而增大，同时随碳纳米管质量分数增加而增大，但表现出 显著的各向异性. 由于引入赤藓糖醇–碳纳米管界面，复合相变材料径向热导率相比纯赤藓糖醇反而降低. 当碳纳米管在赤 藓糖醇中随机分布时，热导率的各向异性得到了显著改善且各方向热导率均得到了提升. 通过对比复合前后赤藓糖醇与碳 纳米管的声子振动态密度发现，由于两者间的相互作用，碳纳米管的声子振动受到抑制，而赤藓糖醇中声子热输运得到激发， 从而提高了热导率. 关键词 赤藓糖醇；碳纳米管；复合相变材料；热导率；分子动力学模拟 分类号 TK02 Simulation of thermal properties of erythritol/carbon nanotube composite phase change materials YAN Xiao-xin1) ，FENG Yan-hui1,2) 苣 ，QIU Lin1,2) 苣 ，ZHANG Xin-xin1,2) 1) School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, China 2) Beijing Key Laboratory of Energy Conservation and Emission Reduction in Metallurgical Industry, University of Science and Technology Beijing, Beijing 100083, China 苣 Corresponding author, FENG Yan-hui, E-mail: yhfeng@me.ustb.edu.cn; QIU Lin, E-mail: qiulin@ustb.edu.cn ABSTRACT In the area of “carbon peaking and carbon neutralization,” changing energy structure from primary energy to new energy is an extremely important issue. Due to the intermittent and fluctuating characteristics of new energy, energy storage technology has proven a viable solution to this issue thus has attracted extensive attention. As a key to energy storage technology, the problem of the low thermal conductivity of phase change materials (PCMs) requires immediate attention. Erythritol is a high enthalpy phase change material commonly used in low-to-medium temperature processes. Its thermal conductivity of only 0.7 W·m–1·K–1 seriously hinders its energy utilization efficiency in practical application. In this paper, erythritol is the main research focus, and single-walled carbon nanotubes (CNTs) of ultra-high thermal conductivity are used as thermal conductivity reinforcements. The effects of length, mass fraction, and the distribution of CNTs on the thermal conductivity of erythritol/CNT composite PCMs were studied by means of molecular dynamics 收稿日期: 2021−09−14 基金项目: 国家重点研发计划资助项目（2018YFA0702300）；国家自然科学基金资助项目（51876007，51876008）；北京市科技新星计划资助 项目（Z201100006820065）；北京市自然科学基金资助项目（3202020）；北京科技大学青年教师学科交叉研究资助项目（中央高校基本科研 业务费专项资金）（FRF-IDRY-19-004） 工程科学学报，第 44 卷，第 X 期：1−8，2022 年 X 月 Chinese Journal of Engineering, Vol. 44, No. X: 1−8, X 2022 https://doi.org/10.13374/j.issn2095-9389.2021.09.14.002; http://cje.ustb.edu.cn
2 工程科学学报，第44卷，第X期 simulation.When the axial lengths of the CNTs were less than their phonon mean free paths,the thermal conductivity of the composite PCMs increased with increasing CNT axial length and mass fraction,although clear anisotropy was exhibited.Due to the introduction of erythritol-CNTs interfacial thermal resistance,the radial thermal conductivity of the composite materials was lower than that of pure erythritol.When CNTs were randomly distributed in erythritol,the anisotropy of thermal conductivity was significantly improved,as was thermal conductivity in all directions.Comparing the phonon vibration densities of the states of erythritol and CNT before and after recombination,it was found that,due to the interaction between the two,the phonon vibration of CNT was suppressed,and the phonon heat transport in erythritol was excited,thus improving thermal conductivity. KEY WORDS erythritol;carbon nanotubes;composite phase change materials;thermal conductivity;molecular dynamic simulation 人类社会正处于飞速发展时期，其背后离不 计算机计算能力的飞速发展，采用计算机模拟已 开能源的支撑与消耗.伴随着一次能源的减少与 经成为一种重要的研究手段.Li等例采用分子动 二氧化碳排放的增加，全球能源与气候变暖状况 力学方法(Molecular dynamics simulation,MD)探究 均不容乐观.作为能源大国，我国提出“碳达峰、 了碳纳米管与不同官能团的相变材料间的相互作 碳中和”的目标，在“双碳”背景下，能源结构从一 用以及热运动特征.Tafrishi等Io借助MD方法测 次能源转向天然气和可再生能源对减少碳排放具 量了石蜡及其复合材料的热性能，结果显示CNT 有重要意义.目前可再生新能源发展势头强劲，然 比石墨烯更好地提升了石蜡的热导率.Yu等叫设 而风能、太阳能等新能源受到自然条件的制约，并 计了一种单壁碳纳米管与NaCI的复合材料并采 且具有间歇性、不稳定的特点，因此有必要将这类 用MD预测了其密度、熔点、导热系数和熔化焓 新能源进行储存.相变储能技术借助相变材料在 等热性能，为熔盐相变材料的设计和应用提供了 特定温度范围内物态转变进行储、释放能量，其原 可能的指导.在采用MD方法预测材料热物性的 理简单、相变过程温度稳定，成为能源领域的研究 研究中，计算热导率通常有两种方法：平衡态分子 热点.然而相变材料作为相变储能技术的主体材 s动力学(Equilibrium molecular dynamics,EMD)方法 料，普遍具有热导率低的问题，这极大限制了能量 和非平衡态分子动力学(Non-equilibrium molecular 利用效率，因此针对如何提高相变材料热导率进 dynamics,.NEMD)方法.平衡态分子动力学方法基 行了大量研究 于Green-Kubo线性响应的涨落-耗散理论计算热 l991年碳纳米管(Carbon nanotubes.CNTs)被 导率，非平衡态分子动力学方法类似于实验方法 首次发现山此后纳米材料与纳米技术蓬勃发展 的测量机制，在模拟单元上施加温度梯度，由于计 纳米材料相比传统材料具有优异的电导率、热导 算速度快被广泛采用-M.NEMD的计算方法一 率、机械强度以及化学稳定性碳纳米管作为纳 般可分为三种，第一种方法是分别设置两个区域 米材料的代表性材料之一，其理论热导率高达 作为热源与热沉，控制热源与热沉的温度，形成温 6000WmlK-13常被用作相变材料的导热增强 度梯度，此时添加到热区的能量等于从冷区中移 材料.在不同种类的石蜡中添加质量分数15%的 除的能量，并且与两个区域间传递的热通量成正 碳纳米管可提高热导率323.2%，在散热器中表现 比，第二种方法是分别对热源和热沉施加/提取 出优异的蓄热降温作用啊.将脂肪酸与功能化的碳 固定量的热能（满足能量守恒定律），得到两个区域 纳米管枝接后，不仅热导率得到了提升，经过 间的温度梯度用于计算热导率6，第三种方法是 2000次循环后仍具有良好的相变稳定性和化学稳 反向扰动NEMD(Reverse non-equilibrium molecular 定性Cong等)在泡沫铜的基础上，使用碳纳米 dynamics,RNEMD)模拟，即Muller--Plathe7提出的 管进一步提高热导率，质量分数为5%的碳纳米管 动量交换算法，通过不同层内的原子速度交换构 可在泡沫铜的基础上将热导率进一步提升1.82 造热流，进而得到温度梯度 倍，达到棕桐酸热导率的32倍.与石墨纳米颗粒 赤藓糖醇(Erythritol))作为中低温区常用的相 相比，碳纳米管作为导热增强材料表现更加优异， 变材料具有熔化焓高的优点，然而其热导率仅有 在相同的添加量下，含有CNT的复合材料表现出 0.7WmlK-118对于赤藓糖醇/碳基复合相变材 更高的热导率圆 料热物性的MD模拟研究较少，目前得的研究大 尽管实验中观察到碳纳米管出色的导热增强 多以实验方法为主，缺少系统的理论研究.针对其 性能，但其背后的影响因素和机制仍需探索.随着 热导率低的问题，本文提出一种赤藓糖醇/单壁
simulation. When the axial lengths of the CNTs were less than their phonon mean free paths, the thermal conductivity of the composite PCMs increased with increasing CNT axial length and mass fraction, although clear anisotropy was exhibited. Due to the introduction of erythritol –CNTs interfacial thermal resistance, the radial thermal conductivity of the composite materials was lower than that of pure erythritol. When CNTs were randomly distributed in erythritol, the anisotropy of thermal conductivity was significantly improved, as was thermal conductivity in all directions. Comparing the phonon vibration densities of the states of erythritol and CNT before and after recombination, it was found that, due to the interaction between the two, the phonon vibration of CNT was suppressed, and the phonon heat transport in erythritol was excited, thus improving thermal conductivity. KEY WORDS erythritol；carbon nanotubes；composite phase change materials；thermal conductivity；molecular dynamic simulation 人类社会正处于飞速发展时期，其背后离不 开能源的支撑与消耗. 伴随着一次能源的减少与 二氧化碳排放的增加，全球能源与气候变暖状况 均不容乐观. 作为能源大国，我国提出“碳达峰、 碳中和”的目标，在“双碳”背景下，能源结构从一 次能源转向天然气和可再生能源对减少碳排放具 有重要意义. 目前可再生新能源发展势头强劲，然 而风能、太阳能等新能源受到自然条件的制约，并 且具有间歇性、不稳定的特点，因此有必要将这类 新能源进行储存. 相变储能技术借助相变材料在 特定温度范围内物态转变进行储、释放能量，其原 理简单、相变过程温度稳定，成为能源领域的研究 热点. 然而相变材料作为相变储能技术的主体材 料，普遍具有热导率低的问题，这极大限制了能量 利用效率，因此针对如何提高相变材料热导率进 行了大量研究. 1991 年碳纳米管（Carbon nanotubes, CNTs）被 首次发现[1] ，此后纳米材料与纳米技术蓬勃发展. 纳米材料相比传统材料具有优异的电导率、热导 率、机械强度以及化学稳定性[2] . 碳纳米管作为纳 米材料的代表性材料之一 ，其理论热导率高达 6000 W·m–1·K–1[3–4] ，常被用作相变材料的导热增强 材料. 在不同种类的石蜡中添加质量分数 15% 的 碳纳米管可提高热导率 323.2%，在散热器中表现 出优异的蓄热降温作用[5] . 将脂肪酸与功能化的碳 纳米管枝接后 ，不仅热导率得到了提升 ，经 过 2000 次循环后仍具有良好的相变稳定性和化学稳 定性[6] . Cong 等[7] 在泡沫铜的基础上，使用碳纳米 管进一步提高热导率，质量分数为 5% 的碳纳米管 可在泡沫铜的基础上将热导率进一步提升 1.82 倍，达到棕桐酸热导率的 32 倍. 与石墨纳米颗粒 相比，碳纳米管作为导热增强材料表现更加优异， 在相同的添加量下，含有 CNT 的复合材料表现出 更高的热导率[8] . 尽管实验中观察到碳纳米管出色的导热增强 性能，但其背后的影响因素和机制仍需探索. 随着 计算机计算能力的飞速发展，采用计算机模拟已 经成为一种重要的研究手段. Li 等[9] 采用分子动 力学方法（Molecular dynamics simulation, MD）探究 了碳纳米管与不同官能团的相变材料间的相互作 用以及热运动特征. Tafrishi 等[10] 借助 MD 方法测 量了石蜡及其复合材料的热性能，结果显示 CNT 比石墨烯更好地提升了石蜡的热导率. Yu 等[11] 设 计了一种单壁碳纳米管与 NaCl 的复合材料并采 用 MD 预测了其密度、熔点、导热系数和熔化焓 等热性能，为熔盐相变材料的设计和应用提供了 可能的指导. 在采用 MD 方法预测材料热物性的 研究中，计算热导率通常有两种方法：平衡态分子 动力学（Equilibrium molecular dynamics, EMD）方法 和非平衡态分子动力学（Non-equilibrium molecular dynamics, NEMD）方法. 平衡态分子动力学方法基 于 Green-Kubo 线性响应的涨落–耗散理论计算热 导率，非平衡态分子动力学方法类似于实验方法 的测量机制，在模拟单元上施加温度梯度，由于计 算速度快被广泛采用[11–14] . NEMD 的计算方法一 般可分为三种，第一种方法是分别设置两个区域 作为热源与热沉，控制热源与热沉的温度，形成温 度梯度，此时添加到热区的能量等于从冷区中移 除的能量，并且与两个区域间传递的热通量成正 比[15] ；第二种方法是分别对热源和热沉施加/提取 固定量的热能（满足能量守恒定律），得到两个区域 间的温度梯度用于计算热导率[16] ；第三种方法是 反向扰动 NEMD（Reverse non-equilibrium molecular dynamics, RNEMD）模拟，即 Müller-Plathe[17] 提出的 动量交换算法，通过不同层内的原子速度交换构 造热流，进而得到温度梯度. 赤藓糖醇（Erythritol）作为中低温区常用的相 变材料具有熔化焓高的优点，然而其热导率仅有 0.7 W·m–1·K–1 [18] . 对于赤藓糖醇/碳基复合相变材 料热物性的 MD 模拟研究较少，目前得的研究大 多以实验方法为主，缺少系统的理论研究. 针对其 热导率低的问题，本文提出一种赤藓糖醇/单壁 · 2 · 工程科学学报，第 44 卷，第 X 期
闫晓鑫等：赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 3 碳纳米管复合相变材料，基于分子动力学模拟手 的长度而定，随碳纳米管长度改变而发生变化.本 段，预测了其热导率，并探究了碳纳米管长度、质 文选择手性为(8-8)的单壁碳纳米管直径d= 量分数、分布方式等因素对复合材料热导率的 1.085nm,为了探究碳纳米管直径对复合材料热导 影响，并且从声子振动角度对其潜在机制进行了 率的影响，分别构建了长度为3.9、5.2、5.9、7.4和 剖析. 8.6nm的单壁碳纳米管，将其分别放入赤薛糖醇 1模型与模拟方法 中，如图1(b)所示.保证体系x、y方向尺寸为 5.124nm×2.562nm不变，体系：方向的尺寸与碳 1.1模型构建 纳米管长度相一致，此时5种模型中碳纳米管的 本文采用Materials studio软件进行建模.赤藓 质量分数均约为13%.为了探究碳纳米管添加量 糖醇单晶胞结构如图1(a)所示，其晶格常数a= 对复合材料热导率的影响，保持碳纳米管长度为 b=1.281nm,c=0.681nm9,对其进行4×2×n扩胞 8.6nm,通过改变碳纳米管周围赤藓糖醇分子的个 操作，得到尺寸为5.124nm×2.562nm×0.681nnm 数改变碳纳米管的添加量.图1(b)展示了质量分 的模拟体系，其中：方向扩胞倍数n依据碳纳米管 数为8.88%和13.35%的两种情况. a 6 图1赤藓糖醇/碳纳米管复合结构模型图.()赤藓糖醇单品胞结构：(b)碳纳米管质量分数分别为8.88%和13.35%的赤藓糖醇碳纳米管复合结 构模型 Fig.1 Erythritol/carbon nanotube composite structure model:(a)erythritol single cell structure;(b)erythritol/carbon nanotube composite structure model with a carbon nanotube mass fraction of 8.88%and 13.35% 1.2模拟方法 表1赤藓糖醇与碳纳米管不同种类原子间L势参数 本文所有的模拟计算均在开源软件Lammps!2o Table 1 L-J potential parameters between erythritol and different kinds 平台进行，OVITO软件包用于可视化处理.在 of carbon nanotubes 分子动力学模拟中，势函数的选择关系到模拟结 Atom type Ey/ev y/nm 果的准确性.本文采用GROMOS54A7力场2四1对 C2-Cp 0.0028 0.3000 赤藓糖醇进行描述，采用Tersoff势函数P1描述单壁 Cr-Cp 0.0028 0.3000 碳纳米管中碳原子间的相互作用，对于赤藓糖醇与 Oh-Cp 0.0033 0.3620 碳纳米管间的相互作用则采用Lennard-Jones(L-J) ho-Cp 0.0021 0.3025 势函数描述： h-Cp 0.0000 0.0000 E(rj)=4E1) (1) 递过程，是一种非平衡过程.因此，在NEMD方法中需 要构建温度梯度，等待体系达到稳态后，对体系热流 式中：E为原子相互作用势能：为原子位置；为 以及温度梯度进行统计，根据傅里叶定律计算热导率： 势阱深度；σ为相互作用势能为零时的平衡距离 (2) 根据Lorentz-Berthelot混合法则，e=√(e·s): oi+oj 式中：k为模拟体系热导率，WmlK;J为体系热 0i= ,具体参数见表1P 流，W:A为垂直于热流方向的横截面积，m:VT为 本文采用NEMD方法计算复合材料的热导率本 体系的温度梯度，Km.本文采用局部热浴法构 质上热传导是由于体系存在温度差而引起的能量传 建温度梯度，由于热浴内部不同的模态能级声子
碳纳米管复合相变材料，基于分子动力学模拟手 段，预测了其热导率，并探究了碳纳米管长度、质 量分数、分布方式等因素对复合材料热导率的 影响，并且从声子振动角度对其潜在机制进行了 剖析. 1 模型与模拟方法 1.1 模型构建 本文采用 Materials studio 软件进行建模. 赤藓 糖醇单晶胞结构如图 1（a）所示，其晶格常数 a = b = 1.281 nm，c = 0.681 nm[19] ，对其进行 4 × 2 × n 扩胞 操作，得到尺寸为 5.124 nm × 2.562 nm × 0.681n nm 的模拟体系，其中 z 方向扩胞倍数 n 依据碳纳米管 的长度而定，随碳纳米管长度改变而发生变化. 本 文选择手性为 （ 8-8）的单壁碳纳米管直 径 d = 1.085 nm，为了探究碳纳米管直径对复合材料热导 率的影响，分别构建了长度为 3.9、5.2、5.9、7.4 和 8.6 nm 的单壁碳纳米管，将其分别放入赤藓糖醇 中 ，如图 1（ b）所示. 保证体系 x、 y 方向尺寸为 5.124 nm × 2.562 nm 不变，体系 z 方向的尺寸与碳 纳米管长度相一致，此时 5 种模型中碳纳米管的 质量分数均约为 13%. 为了探究碳纳米管添加量 对复合材料热导率的影响，保持碳纳米管长度为 8.6 nm，通过改变碳纳米管周围赤藓糖醇分子的个 数改变碳纳米管的添加量. 图 1（b）展示了质量分 数为 8.88% 和 13.35% 的两种情况. (a) (b) 3.81 nm A y z x C B 5.1 nm 13.55 % 8.88 % 2.50 nm 图 1 赤藓糖醇/碳纳米管复合结构模型图. （a）赤藓糖醇单晶胞结构；（b）碳纳米管质量分数分别为 8.88% 和 13.35% 的赤藓糖醇/碳纳米管复合结 构模型 Fig.1 Erythritol/carbon nanotube composite structure model: (a) erythritol single cell structure; (b) erythritol/carbon nanotube composite structure model with a carbon nanotube mass fraction of 8.88% and 13.35% 1.2 模拟方法 本文所有的模拟计算均在开源软件 Lammps[20] 平台进行，OVITO 软件包[21] 用于可视化处理. 在 分子动力学模拟中，势函数的选择关系到模拟结 果的准确性. 本文采用 GROMOS 54A7 力场[22] 对 赤藓糖醇进行描述，采用 Tersoff 势函数[23] 描述单壁 碳纳米管中碳原子间的相互作用，对于赤藓糖醇与 碳纳米管间的相互作用则采用 Lennard–Jones (L–J) 势函数描述： E ( ri j) = 4εi j ( σi j ri j )12 − ( σi j ri j )6 （1） ri j εi j σi j εi j = √( εi · εj ) σi j = ( σi +σj ) 2 式中：E 为原子相互作用势能； 为原子位置； 为 势阱深度； 为相互作用势能为零时的平衡距离. 根据 Lorentz –Berthelot 混合法则 ， ， ，具体参数见表 1 [24] . 本文采用NEMD 方法计算复合材料的热导率. 本 质上热传导是由于体系存在温度差而引起的能量传 递过程，是一种非平衡过程. 因此，在NEMD 方法中需 要构建温度梯度，等待体系达到稳态后，对体系热流 以及温度梯度进行统计，根据傅里叶定律计算热导率： k = − J A∇T （2） ∇T 式中：k 为模拟体系热导率，W·m–1·K–1 ；J 为体系热 流，W；A 为垂直于热流方向的横截面积，m 2 ； 为 体系的温度梯度，K·m–1 . 本文采用局部热浴法构 建温度梯度，由于热浴内部不同的模态能级声子 表 1 赤藓糖醇与碳纳米管不同种类原子间 L-J 势参数 Table 1 L-J potential parameters between erythritol and different kinds of carbon nanotubes Atom type εij / eV σij / nm C2 -Cp 0.0028 0.3000 C1 -Cp 0.0028 0.3000 Oh-Cp 0.0033 0.3620 ho-Cp 0.0021 0.3025 h-Cp 0.0000 0.0000 闫晓鑫等： 赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 · 3 ·
工程科学学报，第44卷，第X期 激发，不同的热浴产生不同的非傅立叶热输运特 使用密度来验证模型以及所选力场的合理性1,27-2 性2.由于随机性和局域性，朗之万热浴(Langevin) 我们计算了293K温度下的赤藓糖醇在驰豫过程 优于Nose-Hoover热浴Pa,因此本文采用Langevin 中的密度，结果如图3所示.在100ps的时间内， 热浴分别对热源与热沉进行控温，计算300K温度 随着驰豫时间的延长，赤藓糖醇的密度逐渐趋于 下的热导率时，热源温度设置为310K,热沉温度 稳定，因此模拟过程中设置的250ps的弛豫时间 设置为290K,整体温差为20K体系热流J=dg 可以使系统达到较为稳定的状态.在293K的温 应满足能量守恒定律，其中E为提取/注入的能量， 度下，赤藓糖醇的密度模拟值为1422.44kgm3,与 J;1为模拟时间，S 文献值1440kgm3[29的误差为1.22%，这说明模 1.3模拟过程 型合理且所选GROMOS54A7力场具有较高的计 图2为NEMD方法计算热导率示意图.在进 算精度 行模拟前首先对体系进行能量最小化处理，使体 1600 系达到能量最低状态.采用Velocity Verlet算法求 Density 1550 Average density 解运动方程，时间步长为0.5s,三个方向均采用周 期性边界条件.模型沿热流方向均分为20层，为 C1500 了防止在模拟过程中体系整体发生漂移，将z轴两 1450 293K1422.44kgm3 端的一层原子设为固定层，分别将紧邻固定层的 两层原子设置为热源与热沉.首先将体系整体置 于NPT系综（等温等压系综，Constant pressure and 1350 constant temperature)下，Nose-Hoover热浴和压浴 分别用于调节系统的温度和压力.在目标温度以 1300 0102030405060708090100 Time/ps 及1个大气压下驰豫250ps,达到稳定状态后撤掉 NPT系综.之后在NVE系综（微正则系综，Micro- 因3293K温度下赤藓糖醇密度随驰豫时间变化 Fig.3 Density of erythritol changes with a relaxation time at 293 K canonical ensemble)下分别对热源与热沉进行 Langevin控温，250ps后热源与热沉温度达到稳定 2.2碳纳米管长度对复合相变材料热导率的影响 状态后开始进行计算各层原子的温度，计算2s, 本文采用NEMD方法计算热导率，向热源/热 取最后1s的数据进行统计，获得体系的温度梯 沉中施加提取的能量如图4(a)所示，满足能量守 度和热流，进而计算得到热导率，每个算例均进行 恒定律，其斜率即为热流.体系的温度分布如图4(b) 了三次独立计算 所示，热源温度为310K,热沉温度为290K,其线 性拟合的斜率即为体系的温度梯度.为了探究碳 纳米管长度对赤藓糖醇/碳纳米管复合材料热导率 的影响，分别计算了碳纳米管长度为3.9、5.2、5.9、 7.4和8.6nm的复合结构的热导率，结果如图5所 示.随着碳纳米管轴向长度的增加，复合相变材料 热导率随之提升.纯赤藓糖醇的热导率NEMD计 算值为(0.67±0.04)WmlK-12,在本文模拟范围 内，：方向（轴向）热导率达到(4.37±0.27)Wm1K-1 (纯赤藓糖醇热导率的6.5倍)，平均热导率达到 (1.66土0.10)Wm-K-1(纯赤藓糖醇热导率的2.5 国2NEMD方法计算热导率示意图 倍)，碳纳米管的加入显著提升了热导率.碳纳米 Fig.2 Schematic model of the NEMD method for thermal conductivity 管的声子平均自由程(Phonon mean free path,MFP) calculation 与其直径成正比30，室温下约为500~700nmB1-刘， 2 结果与讨论 由于本文模拟的碳纳米管尺寸(L)显著小于其声 子平均自由程，此时声子以弹道输运为主，而尺寸 2.1模型密度验证 对这种输运模式具有显著影响)，在弹道输运模 密度是相变材料的重要热物理参数，许多研究 式下，热导率几乎随几何尺寸线性增加，当几何尺
J = dE dt 激发，不同的热浴产生不同的非傅立叶热输运特 性[25] . 由于随机性和局域性，朗之万热浴（Langevin） 优于 Nosé–Hoover 热浴[26] ，因此本文采用 Langevin 热浴分别对热源与热沉进行控温，计算 300 K 温度 下的热导率时，热源温度设置为 310 K，热沉温度 设置为 290 K，整体温差为 20 K. 体系热流 ， 应满足能量守恒定律，其中 E 为提取/注入的能量， J；t 为模拟时间，s. 1.3 模拟过程 图 2 为 NEMD 方法计算热导率示意图. 在进 行模拟前首先对体系进行能量最小化处理，使体 系达到能量最低状态. 采用 Velocity Verlet 算法求 解运动方程，时间步长为 0.5 fs，三个方向均采用周 期性边界条件. 模型沿热流方向均分为 20 层，为 了防止在模拟过程中体系整体发生漂移，将 z 轴两 端的一层原子设为固定层，分别将紧邻固定层的 两层原子设置为热源与热沉. 首先将体系整体置 于 NPT 系综（等温等压系综，Constant pressure and constant temperature）下，Nosé–Hoover 热浴和压浴 分别用于调节系统的温度和压力. 在目标温度以 及 1 个大气压下驰豫 250 ps，达到稳定状态后撤掉 NPT 系综. 之后在 NVE 系综（微正则系综，Micro– canonical ensemble） 下 分 别 对 热 源 与 热 沉 进 行 Langevin 控温，250 ps 后热源与热沉温度达到稳定 状态后开始进行计算各层原子的温度，计算 2 ns， 取最后 1 ns 的数据进行统计，获得体系的温度梯 度和热流，进而计算得到热导率，每个算例均进行 了三次独立计算. x zFixedlayer Heat source Fixedlayer Heat flux Heat sink 图 2 NEMD 方法计算热导率示意图 Fig.2 Schematic model of the NEMD method for thermal conductivity calculation 2 结果与讨论 2.1 模型密度验证 密度是相变材料的重要热物理参数，许多研究 使用密度来验证模型以及所选力场的合理性[11,27–28] . 我们计算了 293 K 温度下的赤藓糖醇在驰豫过程 中的密度，结果如图 3 所示. 在 100 ps 的时间内， 随着驰豫时间的延长，赤藓糖醇的密度逐渐趋于 稳定，因此模拟过程中设置的 250 ps 的弛豫时间 可以使系统达到较为稳定的状态. 在 293 K 的温 度下，赤藓糖醇的密度模拟值为 1422.44 kg·m–3，与 文献值 1440 kg·m–3[29] 的误差为 1.22%，这说明模 型合理且所选 GROMOS 54 A7 力场具有较高的计 算精度. 1600 Density 1550 Average density 1500 1450 1400 1350 1300 0 10 20 30 40 50 Time/ps Density/(kg·m−3 ) 60 70 80 293 K 1422.44 kg·m−3 90 100 图 3 293 K 温度下赤藓糖醇密度随驰豫时间变化 Fig.3 Density of erythritol changes with a relaxation time at 293 K 2.2 碳纳米管长度对复合相变材料热导率的影响 本文采用 NEMD 方法计算热导率，向热源/热 沉中施加/提取的能量如图 4（a）所示，满足能量守 恒定律，其斜率即为热流. 体系的温度分布如图 4（b） 所示，热源温度为 310 K，热沉温度为 290 K，其线 性拟合的斜率即为体系的温度梯度. 为了探究碳 纳米管长度对赤藓糖醇/碳纳米管复合材料热导率 的影响，分别计算了碳纳米管长度为 3.9、5.2、5.9、 7.4 和 8.6 nm 的复合结构的热导率，结果如图 5 所 示. 随着碳纳米管轴向长度的增加，复合相变材料 热导率随之提升. 纯赤藓糖醇的热导率 NEMD 计 算值为（0.67±0.04） W·m–1·K–1 [28] ，在本文模拟范围 内，z 方向（轴向）热导率达到（4.37±0.27） W·m–1·K–1 （纯赤藓糖醇热导率的 6.5 倍），平均热导率达到 （ 1.66±0.10 ）W·m–1·K–1（纯赤藓糖醇热导率的 2.5 倍），碳纳米管的加入显著提升了热导率. 碳纳米 管的声子平均自由程（Phonon mean free path, MFP） 与其直径成正比[30] ，室温下约为 500～700 nm[31–32] ， 由于本文模拟的碳纳米管尺寸（Lz）显著小于其声 子平均自由程，此时声子以弹道输运为主，而尺寸 对这种输运模式具有显著影响[33] ，在弹道输运模 式下，热导率几乎随几何尺寸线性增加，当几何尺 · 4 · 工程科学学报，第 44 卷，第 X 期
闫晓鑫等：赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 5 寸大于声子平均自由程时，声子转为扩散热输运 在本文模拟范围内(L<MFP),复合材料的热导率 模式，此时热导率达到饱和值后趋于稳定因而 随碳纳米管轴向长度增加而显著提高 (a) (b) 2000 -Heat sink Temperature distribution of 310 1500 -Heat source erythritol/CNT (5.2 nm,13.06%) 1000 aLinear fitting curve J=0.97x-481.94 305 500 1=-4.11x+310.82 。-0.9954 -500 -1000 1=-0.97x+486.71 295 -1500 290 -2000 500 1000 1500 2000 2500 0 2 4 Time/ps Length/nm 图4NEMD方法热流与温度分布图.(a)体系施加/提取的能量（斜率为热流）：(b)体系温度分布图 Fig.4 Heat flow and temperature distribution diagram of the NEMD method:(a)energy added/subtracted by the system (the slope is heat flow),(b) temperature distribution of the system direction 赤藓糖醇分子的数量改变其添加量.如图6(a)所 x/y direction 示，随着碳纳米管添加量增加，复合材料整体热导 Mean thermal conductivity 率得到了提升，但是依旧存在显著的各向异性.尽 管碳纳米管长度相同，但是具有不同添加量的复 合材料：方向热导率却随着碳纳米管质量分数增 加而提高.相比于纯碳纳米管，在复合材料中赤藓 糖醇与其相互作用，这会导致碳纳米管自身的热 导率受影响而降低.当碳纳米管质量分数增加时， 其周围的赤藓糖醇数量减少，两者间的相互作用 39 5.2 5.9 7.4 8.6 Length of CNT/nm 减弱，如图6(b)所示（负值代表两者之间呈现相互 图5赤藓糖醇碳纳米管复合材料热导率随碳纳米管长度变化 吸引作用)，对碳纳米管热导率的抑制减弱，进而 Fig.5 Thermal conductivity of erythritol/CNT composites as a function 其轴向热导率得以提升 of the length of CNT 2.4碳纳米管分布方式对复合相变材料热导率的 然而，由于碳纳米管自身热导率具有较强的 影响 各向异性，使得复合材料热导率也表现出明显的 当碳纳米管在相变芯材中规则排布时，热导 各向异性.沿x、y方向(CNT径向)，热导率几乎不 率具有明显的各向异性，而通常实验中碳纳米管 随碳纳米管长度变化，约为0.3Wm1K.由于碳 在相变芯材中随机分布，为此，本文构建了5.1nm× 纳米管的引入，沿x、y方向热量传递时需要经过 5.lnm×5.4nm的模拟，盒子，采用Materials Studio 赤藓糖醇-CNT界面，界面热阻的存在造成了能量 软件中Amorphous Cell模块以接近真实密度向盒 损失，阻碍热传递，进而导致热导率降低B).甚至 子中随机填充赤藓糖醇分子和长度为5m的碳纳 低于纯赤藓糖醇.但整体而言，赤藓糖醇/CNT复 米管，碳纳米管个数分别为3、5、6根，模型结构如 合材料平均热导率仍然呈现一定程度的提升 图7所示.采用NEMD方法分别计算了沿三个方 2.3碳纳米管质量分数对复合相变材料热导率的 向的复合材料热导率结果如表2所示，碳纳米管 影响 随机分布方式下热导率的各向异性得到了明显改 碳纳米管作为导热增强材料，其添加的质量 善，沿各个方向的热导率均得到了提升，平均热导 分数对复合材料热导率具有显著影响.为排除碳 率随碳纳米管个数增加而升高，但增加趋势变缓. 纳米管长度对热导率的影响，本文在保证碳纳米 2.5声子振动态密度 管长度L=8.6nm的前提下，仅改变碳纳米管周围 对于非金属材料，声子作为热量传递的主要
寸大于声子平均自由程时，声子转为扩散热输运 模式，此时热导率达到饱和值后趋于稳定[34] . 因而 在本文模拟范围内（Lz < MFP），复合材料的热导率 随碳纳米管轴向长度增加而显著提高. 2000 (a) (b) 310 305 300 295 290 0 1 2 3 4 5 1500 1000 500 500 1000 Heat sink Heat source y=0.97x−481.94 y=−0.97x+486.71 y=−4.11x+310.82 R 2=0.9954 Length/nm Temperature/K Time/ps Energy/eV 1500 2000 2500 Temperature distribution of erythritol/CNT (5.2 nm, 13.06%) Linear fitting curve 0 −500 −1000 −1500 −2000 图 4 NEMD 方法热流与温度分布图. （a）体系施加/提取的能量（斜率为热流）；（b）体系温度分布图 Fig.4 Heat flow and temperature distribution diagram of the NEMD method: (a) energy added/subtracted by the system (the slope is heat flow); (b) temperature distribution of the system 5 Length of CNT/nm z direction x/y direction Mean thermal conductivity Thermal conductivity/(W·m−1K−1 ) 4 3 2 1 0 3.9 5.2 5.9 7.4 8.6 图 5 赤藓糖醇/碳纳米管复合材料热导率随碳纳米管长度变化 Fig.5 Thermal conductivity of erythritol/CNT composites as a function of the length of CNT 然而，由于碳纳米管自身热导率具有较强的 各向异性，使得复合材料热导率也表现出明显的 各向异性. 沿 x、y 方向（CNT 径向），热导率几乎不 随碳纳米管长度变化，约为 0.3 W·m–1·K–1 . 由于碳 纳米管的引入，沿 x、y 方向热量传递时需要经过 赤藓糖醇–CNT 界面，界面热阻的存在造成了能量 损失，阻碍热传递，进而导致热导率降低[35] ，甚至 低于纯赤藓糖醇. 但整体而言，赤藓糖醇/CNT 复 合材料平均热导率仍然呈现一定程度的提升. 2.3 碳纳米管质量分数对复合相变材料热导率的 影响 碳纳米管作为导热增强材料，其添加的质量 分数对复合材料热导率具有显著影响. 为排除碳 纳米管长度对热导率的影响，本文在保证碳纳米 管长度 L = 8.6 nm 的前提下，仅改变碳纳米管周围 赤藓糖醇分子的数量改变其添加量. 如图 6（a）所 示，随着碳纳米管添加量增加，复合材料整体热导 率得到了提升，但是依旧存在显著的各向异性. 尽 管碳纳米管长度相同，但是具有不同添加量的复 合材料 z 方向热导率却随着碳纳米管质量分数增 加而提高. 相比于纯碳纳米管，在复合材料中赤藓 糖醇与其相互作用，这会导致碳纳米管自身的热 导率受影响而降低. 当碳纳米管质量分数增加时， 其周围的赤藓糖醇数量减少，两者间的相互作用 减弱，如图 6（b）所示（负值代表两者之间呈现相互 吸引作用），对碳纳米管热导率的抑制减弱，进而 其轴向热导率得以提升. 2.4 碳纳米管分布方式对复合相变材料热导率的 影响 当碳纳米管在相变芯材中规则排布时，热导 率具有明显的各向异性，而通常实验中碳纳米管 在相变芯材中随机分布，为此，本文构建了 5.1 nm × 5.1 nm × 5.4 nm 的模拟盒子，采用 Materials Studio 软件中 Amorphous Cell 模块以接近真实密度向盒 子中随机填充赤藓糖醇分子和长度为 5 nm 的碳纳 米管，碳纳米管个数分别为 3、5、6 根，模型结构如 图 7 所示. 采用 NEMD 方法分别计算了沿三个方 向的复合材料热导率结果如表 2 所示，碳纳米管 随机分布方式下热导率的各向异性得到了明显改 善，沿各个方向的热导率均得到了提升，平均热导 率随碳纳米管个数增加而升高，但增加趋势变缓. 2.5 声子振动态密度 对于非金属材料，声子作为热量传递的主要 闫晓鑫等： 赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 · 5 ·
工程科学学报，第44卷，第X期 (a) (b)30 ￥ =direction x/y direction Mean thermal conductivity -29 2.95 4.256.658.8811.86 13.35 2.95 4.25 6.658.88 11.86 13.35 Mass fraction of CNT/% Mass fraction of CNT/% 图6()赤藓糖醇/碳纳米管复合材料热导率随碳纳米管质量分数变化：(6)赤藓糖醇与碳纳米管间相互作用能随碳纳米管质量分数变化 Fig.6 (a)Thermal conductivity of erythritol/CNT composites as a function of the mass fraction of CNT;(b)interaction energy between erythritol and CNT varies with the mass fraction of CNT (a) 6 (c) 图7赤藓糖醇碳纳米管随机填充复合模型.()赤藓糖醇/碳纳米管(3根)复合材料随机模型：(b)赤藓糖醇/碳纳米管(5根)复合材料随机模型： (c)赤藓糖醇碳纳米管(6根)复合材料随机模型 Fig.7 Random filling model of erythritol/CNT:(a)random filling model of erythritol/3 CNT composites;(b)random filling model of erythritol/5 CNT composites;(c)random filling model of erythritol/6 CNT composites 表2随机分布方式下赤藓糖醇/碳纳米管复合材料热导率 Table 2 Thermal conductivity of erythritol/CNT composites with random distribution Number of CNTs k/(Wm-K-)k/(W'm-K-) k./(W-m-K-) kn/(W'mK-) Enhancement ratio / 3 (mass fraction of 16.00%) 1.61±0.03 1.61±0.07 1.29±0.03 1.50±0.01 123.88 5 (mass fraction of 25.28%) 2.28±0.01 1.61±0.01 1.92±0.06 1.94±0.02 189.56 6 (mass fraction of 30.33%) 2.27±0.11 2.33±0.10 1.87±0.08 2.16±0.06 222.39 载体，承载着大量的能量.声子振动态密度 v为声子频率 (Vibration density of states,.VDOS)可以得到原子的 本文分别计算了赤藓糖醇以及碳纳米管在复 振动情况，利用声子振动态密度对材料的热导率 合前后的声子振动态密度，如图8所示.碳纳米管 变化机制进行分析是目前常用的手段之一B6 在复合材料中多处声子振动峰值降低，在低频区 声子振动态密度由原子的速度自相关函数 尤为明显，而对于纳米碳材料而言，低频区声子振 (Velocity autocorrelation function,VACF)傅里叶变 动具有较大的声子平均自由程，较长的能量传输 换得到[0， 距离可减少声子碰撞导致的热阻，对于热输运具 有重要作用的由于与赤藓糖醇相互作用，碳纳米 VACF(t)= (成(o)o+） (3) (14(o)-d(o》 管声子振动受到抑制，导致其在复合材料中的热 VDOSVACF 导率远低于独立的碳纳米管，从而复合材料的热 (4) 导率仅在10°~10Wm1K-1量级.在赤藓糖醇中 式中，N是原子个数，to)是第i个原子在to时刻 0~66THz范围内的声子对其热导率起主导作用28， 的速度矢量，<>代表初始时刻为的时间平均， 在此频率范围内，相比纯赤藓糖醇，复合材料中的
载体 ，承载着大量的能量 . 声子振动态密度 （Vibration density of states, VDOS）可以得到原子的 振动情况，利用声子振动态密度对材料的热导率 变化机制进行分析是目前常用的手段之一[36– 39] . 声 子 振 动 态 密 度 由 原 子 的 速 度 自 相 关 函 数 （Velocity autocorrelation function, VACF）傅里叶变 换得到[40] ： VACF(t) = ⟨∑N i=1 −→ui (t0)· −→ui (t0 +t) ⟩ ⟨∑N i=1 −→ui (t0)· −→ui (t0) ⟩ （3） VDOS = w + ∞ −∞ VACF(t) e −2π ivtdt （4） −→ui 式中，N 是原子个数， (t0) 是第 i 个原子在 t0 时刻 的速度矢量，代表初始时刻为 t0 的时间平均， v 为声子频率. 本文分别计算了赤藓糖醇以及碳纳米管在复 合前后的声子振动态密度，如图 8 所示. 碳纳米管 在复合材料中多处声子振动峰值降低，在低频区 尤为明显，而对于纳米碳材料而言，低频区声子振 动具有较大的声子平均自由程，较长的能量传输 距离可减少声子碰撞导致的热阻，对于热输运具 有重要作用[35] . 由于与赤藓糖醇相互作用，碳纳米 管声子振动受到抑制，导致其在复合材料中的热 导率远低于独立的碳纳米管，从而复合材料的热 导率仅在 100～101 W·m–1·K–1 量级. 在赤藓糖醇中 0～66 THz 范围内的声子对其热导率起主导作用[28] ， 在此频率范围内，相比纯赤藓糖醇，复合材料中的 表 2 随机分布方式下赤藓糖醇/碳纳米管复合材料热导率 Table 2 Thermal conductivity of erythritol/CNT composites with random distribution Number of CNTs kx / (W·m–1·K–1) ky / (W·m–1·K–1) kz / (W·m–1·K–1) kmean / (W·m–1·K–1) Enhancement ratio / % 3 (mass fraction of 16.00%) 1.61±0.03 1.61±0.07 1.29±0.03 1.50±0.01 123.88 5 (mass fraction of 25.28%) 2.28±0.01 1.61±0.01 1.92±0.06 1.94±0.02 189.56 6 (mass fraction of 30.33%) 2.27±0.11 2.33±0.10 1.87±0.08 2.16±0.06 222.39 5 (a) (b) −30 Interaction energy/eV−29 −28 −27 −26 −25 0 2.95 4.25 6.65 8.88 11.86 13.35 Thermal conductivity/(W·m−1K−1 ) 4 3 2 1 0 2.95 4.25 6.65 Mass fraction of CNT/% Mass fraction of CNT/% 8.88 11.86 13.35 z direction x/y direction Mean thermal conductivity 图 6 （a）赤藓糖醇/碳纳米管复合材料热导率随碳纳米管质量分数变化；（b）赤藓糖醇与碳纳米管间相互作用能随碳纳米管质量分数变化 Fig.6 (a) Thermal conductivity of erythritol/CNT composites as a function of the mass fraction of CNT; (b) interaction energy between erythritol and CNT varies with the mass fraction of CNT (a) (b) (c) 图 7 赤藓糖醇/碳纳米管随机填充复合模型. （a）赤藓糖醇/碳纳米管（3 根）复合材料随机模型；（b）赤藓糖醇/碳纳米管（5 根）复合材料随机模型； （c）赤藓糖醇/碳纳米管（6 根）复合材料随机模型 Fig.7 Random filling model of erythritol/CNT: (a) random filling model of erythritol/3 CNT composites; (b) random filling model of erythritol/5 CNT composites; (c) random filling model of erythritol/6 CNT composites · 6 · 工程科学学报，第 44 卷，第 X 期
闫晓鑫等：赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 7 赤藓糖醇声子振动在16、38和63THz处得到了 发了赤藓糖醇中的声子热输运，在一定程度上增 明显加强，振动峰值明显增高，碳纳米管的加入激 强导热 (a)0.040 (b)0.07 -CNT in composite -Erythritol in composite 0.035 ---Pristine CNT 0.06 -Pristine erythritol 0.030 0.05 0.025 A 0.04 0.020 0.03 0.015 0.02 0010 0005 0.01 0 20 40 60 80 100 20 40 60 80 100 Frequency/THz Frequency/THz 图8 赤藓糖醇碳纳米管复合材料声子振动态密度.()碳纳米管复合前后声子振动态密度：(b)赤藓糖醇复合前后声子振动态密度 Fig.8 VDOS oferythritol/CNT composite:(a)VDOS of CNTs before and after compound:(b)VDOS oferythritol before and after compound 3结论 4613 [5] Xu C,Xu S L,Eticha R D.Experimental investigation of thermal (1)碳纳米管作为导热增强材料可有效提高赤 performance for pulsating flow in a microchannel heat sink filled 藓糖醇的热导率，在其声子平均自由程范围内，赤 with PCM (paraffin/CNT composite).Energy Comvers Manag, 藓糖醇/碳纳米管复合相变材料热导率随碳纳米管 2021,236:114071 轴向长度增加而显著提高，同时随碳纳米管质量 [6] Al-Ahmed A,San A,Mazumder M A J,et al.Thermal energy 分数增大而提高.碳纳米管长度为8.6nm,质量分 storage and thermal conductivity properties of fatty acid/fatty acid- 数为1335%时，复合材料沿碳纳米管轴向热导率 grafted-CNTs and fatty acid/CNTs as novel composite phase 可达到纯赤藓糖醇的6.5倍，平均热导率可达到纯 change materials.Sci Rep,2020,10:15388 [7] Cong R S,Xu C L,Chen Y K,et al.Enhanced thermal 赤藓糖醇的2.5倍 conductivity of palmitic acid/copper foam composites with carbon (2)碳纳米管规则排列时复合相变材料热导率 nanotube as thermal energy storage materials.J Energy Storage, 表现出明显的各向异性，由于引入界面热阻导致 2021,40:102783 沿碳纳米管径向复合相变材料热导率反而低于赤 [8]Rabady RI,Malkawi D S.Thermal conductivity enhancement of 薛糖醇，碳纳米管随机分布方式可有效解决热导 sodium thiosulfate pentahydrate by adding carbon nano 率的各向异性，沿各个方向热导率均得到提升 tubes/graphite nano-particles.Energy Srorage,2020,27:101166 (3)由于碳纳米管与赤藓糖醇间的相互作用， [9] Li A.Hai G T.Cheng P,et al.Molecular insights into the interaction mechanism between C18 phase change materials and 碳纳米管的声子振动受到抑制，表现出远低于独 methyl-modified carbon nanotubes.Ceram Int,2021,47(16): 立碳纳米管的热导率，但赤藓糖醇中多处声子振 23564 动被激发，增强了赤藓糖醇中的声子热输运，从而 [10]Tafrishi H.Sadeghzadeh S,Ahmadi R.et al.Investigation of 促进了热量传递 tetracosane thermal transport in presence of graphene and carbon nanotube fillers-Amolecular dynamics study./Energy Storage, 参考文献 2020,29:101321 [1]lijima S.Helical microtubules of graphitic carbon.Nature,1991. [11]Yu YS,ZhaoC Y,Tao Y B,et al.Superior thermal energy storage 354(6348):56 performance of NaC-SWCNT composite phase change materials: [2]Qiu L,Zhu N,Feng Y H,et al.A review of recent advances in A molecular dynamics approach..4 ppl Energ☒，2021,290：116799 thermophysical properties at the nanoscale:From solid state to [12]Du Y P,Zhou T,Zhao C Y,et al.Molecular dynamics simulation colloids.PIrys Rep,2020,843:1 on thermal enhancement for carbon nano tubes (CNTs)based [3]Hone J,Whitney M,Piskoti C,et al.Thermal conductivity of phase change materials (PCMs).Int J Heat Mass Transf,2022 single-walled carbon nanotubes.Phys Rev B,1999,59(4):R2514 182:122017 [4]Berber S,Kwon Y K,Tomanek D.Unusually high thermal [13]Zhang H C,Rindt CC M,Smeulders D M J,et al.Nanoscale heat conductivity of carbon nanotubes.Phys Rev Lert,2000,84(20): transfer in carbon nanotubes-sugar alcohol composite as heat
赤藓糖醇声子振动在 16、38 和 63 THz 处得到了 明显加强，振动峰值明显增高，碳纳米管的加入激 发了赤藓糖醇中的声子热输运，在一定程度上增 强导热. 0.040 (a) (b) 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0 0.035 0.030 0.025 0.020 0.015 0.010 0.005 0 0 20 40 Frequency/THz 60 80 CNT in composite Pristine CNT 100 0 20 40 Frequency/THz VDOS/(a.u.) VDOS/(a.u.) 60 Erythritol in composite Pristine erythritol 80 100 图 8 赤藓糖醇/碳纳米管复合材料声子振动态密度. （a）碳纳米管复合前后声子振动态密度；（b）赤藓糖醇复合前后声子振动态密度 Fig.8 VDOS of erythritol/CNT composite: (a) VDOS of CNTs before and after compound; (b) VDOS of erythritol before and after compound 3 结论 (1) 碳纳米管作为导热增强材料可有效提高赤 藓糖醇的热导率，在其声子平均自由程范围内，赤 藓糖醇/碳纳米管复合相变材料热导率随碳纳米管 轴向长度增加而显著提高，同时随碳纳米管质量 分数增大而提高. 碳纳米管长度为 8.6 nm，质量分 数为 13.35% 时，复合材料沿碳纳米管轴向热导率 可达到纯赤藓糖醇的 6.5 倍，平均热导率可达到纯 赤藓糖醇的 2.5 倍. (2) 碳纳米管规则排列时复合相变材料热导率 表现出明显的各向异性，由于引入界面热阻导致 沿碳纳米管径向复合相变材料热导率反而低于赤 藓糖醇，碳纳米管随机分布方式可有效解决热导 率的各向异性，沿各个方向热导率均得到提升. (3) 由于碳纳米管与赤藓糖醇间的相互作用， 碳纳米管的声子振动受到抑制，表现出远低于独 立碳纳米管的热导率，但赤藓糖醇中多处声子振 动被激发，增强了赤藓糖醇中的声子热输运，从而 促进了热量传递. 参 考 文 献 Iijima S. Helical microtubules of graphitic carbon. Nature, 1991, 354（6348）: 56 [1] Qiu L, Zhu N, Feng Y H, et al. A review of recent advances in thermophysical properties at the nanoscale: From solid state to colloids. Phys Rep, 2020, 843: 1 [2] Hone J, Whitney M, Piskoti C, et al. Thermal conductivity of single-walled carbon nanotubes. Phys Rev B, 1999, 59（4）: R2514 [3] Berber S, Kwon Y K, Tománek D. Unusually high thermal conductivity of carbon nanotubes. Phys Rev Lett, 2000, 84（20）: [4] 4613 Xu C, Xu S L, Eticha R D. Experimental investigation of thermal performance for pulsating flow in a microchannel heat sink filled with PCM (paraffin/CNT composite). Energy Convers Manag, 2021, 236: 114071 [5] Al-Ahmed A, Sarı A, Mazumder M A J, et al. Thermal energy storage and thermal conductivity properties of fatty acid/fatty acidgrafted-CNTs and fatty acid/CNTs as novel composite phase change materials. Sci Rep, 2020, 10: 15388 [6] Cong R S, Xu C L, Chen Y K, et al. Enhanced thermal conductivity of palmitic acid/copper foam composites with carbon nanotube as thermal energy storage materials. J Energy Storage, 2021, 40: 102783 [7] Rabady R I, Malkawi D S. Thermal conductivity enhancement of sodium thiosulfate pentahydrate by adding carbon nanotubes/graphite nano-particles. J Energy Storage, 2020, 27: 101166 [8] Li A, Hai G T, Cheng P, et al. Molecular insights into the interaction mechanism between C18 phase change materials and methyl-modified carbon nanotubes. Ceram Int, 2021, 47（16）: 23564 [9] Tafrishi H, Sadeghzadeh S, Ahmadi R, et al. Investigation of tetracosane thermal transport in presence of graphene and carbon nanotube fillers—A molecular dynamics study. J Energy Storage, 2020, 29: 101321 [10] Yu Y S, Zhao C Y, Tao Y B, et al. Superior thermal energy storage performance of NaCl–SWCNT composite phase change materials: A molecular dynamics approach. Appl Energy, 2021, 290: 116799 [11] Du Y P, Zhou T, Zhao C Y, et al. Molecular dynamics simulation on thermal enhancement for carbon nano tubes (CNTs) based phase change materials (PCMs). Int J Heat Mass Transf, 2022, 182: 122017 [12] Zhang H C, Rindt C C M, Smeulders D M J, et al. Nanoscale heat transfer in carbon nanotubes - sugar alcohol composite as heat [13] 闫晓鑫等： 赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 · 7 ·
工程科学学报，第44卷，第X期 storage materials.J Plys Chem C,2016,120(38):21915 [28]Yan XX,Feng Y H,Qiu L,et al.Thermal conductivity and phase [14]Zou HY,Feng Y H,Qiu L,et al.Effect of the loading amount and change characteristics of hierarchical porous diamond/erythritol arrangement of iodine chains on the interfacial thermal transport of composite phase change materials.Energy,2021,233:121158 carbon nanotubes:A molecular dynamics study.RSC Adv,2020. [29]Hohlein S,Konig-Haagen A,Bruggemann D.Thermophysical 10(72):44196 characterization of MgCl,6H,O,xylitol and erythritol as phase [15]Ikeshoji T,Hafskjold B.Non-equilibrium molecular dynamics change materials (PCM)for latent heat thermal energy storage calculation of heat conduction in liquid and through liquid-gas (LHTES).Materials,2017,10(4):444 interface.Mol Phys,1994,81(2):251 [30]Chai Y,Xiao Z Y,Chan P C H.Horizontally aligned carbon [16]Jund P,Jullien R.Molecular-dynamics calculation of the thermal nanotube bundles for interconnect application:Diameter- conductivity of vitreous silica.Phrys Rev B,1999,59(21):13707 dependent contact resistance and mean free path.Nanotechnology, [17]Muller-Plathe F.A simple nonequilibrium molecular dynamics 2010,21(23):235705 method for calculating the thermal conductivity.J Chem Plnys [31]Yu C,Shi L,Yao Z,et al.Thermal conductance and thermopowe 1997,106(14):6082 of an individual single-wall carbon nanotube.Nano Lett,2005. [18]Zhang Q,Luo Z L.Guo QL.et al.Preparation and thermal 5(9):1842 properties of short carbon fibers/erythritol phase change materials. [32]Saaskilahti K,Oksanen J,Volz S,et al.Frequency-dependent Energy Convers Manag,2017,136:220 phonon mean free path in carbon nanotubes from nonequilibrium [19]Shimada A.Crystal and molecular structure of mesoerythritol. molecular dynamics.Phys Rev B,2015,91(11):115426 4 cta Cryst,1958,11(10:748 [33]Wang J,Li C,Li J,et al.A multiscale study of the filler-size and [20]Plimpton S.Fast parallel algorithms for short-range molecular temperature dependence of the thermal conductivity of graphene- dynamics.J Comput Phys,1995,117(1):1 polymer nanocomposites.Carbon,2021,175:259 [21]Stukowski A.Visualization and analysis of atomistic simulation [34]Bae M H,Li Z Y,Aksamija Z,et al.Ballistic to diffusive data with OVITO-the open visualization tool.Modelling Simul crossover of heat flow in graphene ribbons.Nat Commun,2013,4: Mater Sci Eng,2010,18(1:015012 1734 [22]Schmid N,Eichenberger A P.Choutko A,et al.Definition and [35]Qiu L,Zhang X H,Guo Z X,et al.Interfacial heat transport in testing of the GROMOS force-field versions 54A7 and 54B7.Eu nano-carbon assemblies.Carbon,2021,178:391 Biophys,J2011,40(7):843 [23]Tersoff J.New empirical approach for the structure and energy of [36]Feng D L,Feng Y H,Liu Y Z,et al.Thermal conductivity of a 2D covalent systems.Phys Rev B,1988,37(12):6991 covalent organic framework and its enhancement using fullerene [24]Qiu L,Zhu N,Feng Y H,et al.Interfacial thermal transport 3D self-assembly:A molecular dynamics simulation.Plys Chem properties of polyurethane/carbon nanotube hybrid composites.In C,2020,124(15):8386 J Heat Mass Transf,2020,152:119565 [37]Yousefi F,Khoeini F,Rajabpour A.Thermal conductivity and [25]Hu Y,Feng T L,Gu X K,et al.Unification of nonequilibrium thermal rectification of nanoporous graphene:A molecular molecular dynamics and the mode-resolved phonon Boltzmann dynamics simulation.Int J Heat Mass Transf,2020,146:118884 equation for thermal transport simulations.Phys Rev B,2020, [38]Yu Z P,Feng Y H,Feng D L,et al.Thermal conductance 101(15):155308 bottleneck of a three dimensional graphene-CNT hybrid structure: [26]LiZ,Xiong S Y,Sievers C,et al.Influence of thermostatting on A molecular dynamics simulation.Phys Chem Chem Phys,2019, nonequilibrium molecular dynamics simulations of heat 22(1):337 conduction in solids.JChem Phys,2019,151(23):234105 [39]Zhao C Y,Tao Y B,Yu Y S.Molecular dynamics simulation of [27]Feng B.Fan L W.Zeng Y,et al.Atomistic insights into the effects thermal and phonon transport characteristics of nanocomposite of hydrogen bonds on the melting process and heat conduction of phase change material.JMol Lig,2021,329:115448 erythritol as a promising latent heat storage material.IntJ Therm [40]Dickey J M,Paskin A.Computer simulation of the lattice Sci,2019,146:106103 dynamics of solids.Phys Rev,1969,188(3):1407
storage materials. J Phys Chem C, 2016, 120（38）: 21915 Zou H Y, Feng Y H, Qiu L, et al. Effect of the loading amount and arrangement of iodine chains on the interfacial thermal transport of carbon nanotubes: A molecular dynamics study. RSC Adv, 2020, 10（72）: 44196 [14] Ikeshoji T, Hafskjold B. Non-equilibrium molecular dynamics calculation of heat conduction in liquid and through liquid-gas interface. Mol Phys, 1994, 81（2）: 251 [15] Jund P, Jullien R. Molecular-dynamics calculation of the thermal conductivity of vitreous silica. Phys Rev B, 1999, 59（21）: 13707 [16] Müller-Plathe F. A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity. J Chem Phys, 1997, 106（14）: 6082 [17] Zhang Q, Luo Z L, Guo Q L, et al. Preparation and thermal properties of short carbon fibers/erythritol phase change materials. Energy Convers Manag, 2017, 136: 220 [18] Shimada A. Crystal and molecular structure of mesoerythritol. Acta Cryst, 1958, 11（10）: 748 [19] Plimpton S. Fast parallel algorithms for short-range molecular dynamics. J Comput Phys, 1995, 117（1）: 1 [20] Stukowski A. Visualization and analysis of atomistic simulation data with OVITO–the open visualization tool. Modelling Simul Mater Sci Eng, 2010, 18（1）: 015012 [21] Schmid N, Eichenberger A P, Choutko A, et al. Definition and testing of the GROMOS force-field versions 54A7 and 54B7. Eur Biophys J, 2011, 40（7）: 843 [22] Tersoff J. New empirical approach for the structure and energy of covalent systems. Phys Rev B, 1988, 37（12）: 6991 [23] Qiu L, Zhu N, Feng Y H, et al. Interfacial thermal transport properties of polyurethane/carbon nanotube hybrid composites. Int J Heat Mass Transf, 2020, 152: 119565 [24] Hu Y, Feng T L, Gu X K, et al. Unification of nonequilibrium molecular dynamics and the mode-resolved phonon Boltzmann equation for thermal transport simulations. Phys Rev B, 2020, 101（15）: 155308 [25] Li Z, Xiong S Y, Sievers C, et al. Influence of thermostatting on nonequilibrium molecular dynamics simulations of heat conduction in solids. J Chem Phys, 2019, 151（23）: 234105 [26] Feng B, Fan L W, Zeng Y, et al. Atomistic insights into the effects of hydrogen bonds on the melting process and heat conduction of erythritol as a promising latent heat storage material. Int J Therm Sci, 2019, 146: 106103 [27] Yan X X, Feng Y H, Qiu L, et al. Thermal conductivity and phase change characteristics of hierarchical porous diamond/erythritol composite phase change materials. Energy, 2021, 233: 121158 [28] Höhlein S, König-Haagen A, Brüggemann D. Thermophysical characterization of MgCl2 ·6H2O, xylitol and erythritol as phase change materials (PCM) for latent heat thermal energy storage (LHTES). Materials, 2017, 10（4）: 444 [29] Chai Y, Xiao Z Y, Chan P C H. Horizontally aligned carbon nanotube bundles for interconnect application: Diameterdependent contact resistance and mean free path. Nanotechnology, 2010, 21（23）: 235705 [30] Yu C, Shi L, Yao Z, et al. Thermal conductance and thermopower of an individual single-wall carbon nanotube. Nano Lett, 2005, 5（9）: 1842 [31] Sääskilahti K, Oksanen J, Volz S, et al. Frequency-dependent phonon mean free path in carbon nanotubes from nonequilibrium molecular dynamics. Phys Rev B, 2015, 91（11）: 115426 [32] Wang J, Li C, Li J, et al. A multiscale study of the filler-size and temperature dependence of the thermal conductivity of graphenepolymer nanocomposites. Carbon, 2021, 175: 259 [33] Bae M H, Li Z Y, Aksamija Z, et al. Ballistic to diffusive crossover of heat flow in graphene ribbons. Nat Commun, 2013, 4: 1734 [34] Qiu L, Zhang X H, Guo Z X, et al. Interfacial heat transport in nano-carbon assemblies. Carbon, 2021, 178: 391 [35] Feng D L, Feng Y H, Liu Y Z, et al. Thermal conductivity of a 2D covalent organic framework and its enhancement using fullerene 3D self-assembly: A molecular dynamics simulation. J Phys Chem C, 2020, 124（15）: 8386 [36] Yousefi F, Khoeini F, Rajabpour A. Thermal conductivity and thermal rectification of nanoporous graphene: A molecular dynamics simulation. Int J Heat Mass Transf, 2020, 146: 118884 [37] Yu Z P, Feng Y H, Feng D L, et al. Thermal conductance bottleneck of a three dimensional graphene–CNT hybrid structure: A molecular dynamics simulation. Phys Chem Chem Phys, 2019, 22（1）: 337 [38] Zhao C Y, Tao Y B, Yu Y S. Molecular dynamics simulation of thermal and phonon transport characteristics of nanocomposite phase change material. J Mol Liq, 2021, 329: 115448 [39] Dickey J M, Paskin A. Computer simulation of the lattice dynamics of solids. Phys Rev, 1969, 188（3）: 1407 [40] · 8 · 工程科学学报，第 44 卷，第 X 期