赤藓糖醇-碳纳米管复合相变材料热特性模拟研究.docxVIP

? ? 赤藓糖醇/碳纳米管复合相变材料热特性模拟研究 ? ? 闫晓鑫，冯妍卉，邱 琳，张欣欣 1) 北京科技大学能源与环境工程学院，北京 100083 2) 北京科技大学冶金工业节能减排北京市重点实验室，北京 100083 人类社会正处于飞速发展时期，其背后离不开能源的支撑与消耗.伴随着一次能源的减少与二氧化碳排放的增加，全球能源与气候变暖状况均不容乐观.作为能源大国，我国提出“碳达峰、碳中和”的目标，在“双碳”背景下，能源结构从一次能源转向天然气和可再生能源对减少碳排放具有重要意义.目前可再生新能源发展势头强劲，然而风能、太阳能等新能源受到自然条件的制约，并且具有间歇性、不稳定的特点，因此有必要将这类新能源进行储存.相变储能技术借助相变材料在特定温度范围内物态转变进行储、释放能量，其原理简单、相变过程温度稳定，成为能源领域的研究热点.然而相变材料作为相变储能技术的主体材料，普遍具有热导率低的问题，这极大限制了能量利用效率，因此针对如何提高相变材料热导率进行了大量研究. 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模拟研究较少，目前的研究大多以实验方法为主，缺少系统的理论研究.针对其热导率低的问题，本文提出一种赤藓糖醇/单壁碳纳米管复合相变材料，基于分子动力学模拟手段，预测了其热导率，并探究了碳纳米管长度、质量分数、分布方式等因素对复合材料热导率的影响，并且从声子振动角度对其潜在机制进行了剖析. 1 模型与模拟方法 1.1 模型构建 本文采用Materials studio软件进行建模.赤藓糖醇单晶胞结构如图1（a）所示，其晶格常数a=b= 1.281 nm，c= 0.681 nm[