纳米变截面结构的导热性质.pdfVIP

中困工程热物理学会 传热传质学 2006年学术会议 编呼：063．．8 纳米变截面结构的导热性质 曹炳阳，张清光，过增元 (清华人学航天航空学院，传热与能源利用北京市重点实验室，北京100084) 摘要：采用分子动力学方法模拟了固态氩的纳米变截面导热性质，研究发现纳米变截面材料的热阻 和热流的大小与方向都相关：当纳米结构沿热流的方向为渐缩时，纳米结构的热阻随热流的增加而增 人，析当纳米结构沿热流的方向为渐扩时．纳米结构的热阻随热流的增加呈减小的趋势；当热流较人 时，热流沿渐缩方向时的热阻明显大于热流沿渐扩方向时的热阻，但当热流较小时纳米变截面结构的 热阻和热流方向的关系不人。最后依据热质的运动和传递理论的动能效应对该现象进行丫分析解释。 关键词：导热性质纳米变截面结构动能效应分子动力学模拟 1引言 近年来，微电子、光电子和微／纳机电系统(MEMS／NEMS)等微纳米科学与技术发 展迅速，其中涉及到重要的热设计、热管理和热优化问剧¨。微纳米结构材料，如微纳 米薄膜、纳米线、纳米管等的导热性质研究备受研究人员所关注，目前已揭示出一些特 别的导热规律，例如由于边界、晶界和缺陷等因素对载热粒子的散射作用而导致微纳薄 膜和纳米线的热导率较之体材料大大降低，而由于低维结构对载热粒子输运模式的影响 使得纳米管常具有非常高的热导率等。热流密度较高是纳米结构材料导热的一个重要特 点，过增元等基于热质概念对导热规律的研究表明热的动能效应在高热流密度时会显著 影响导热规律12J，并且，纳米变截面结构材料由于截面变化使得热量传递的速度变化增 大从而动能效应会更为显著。目前，关于变截面对纳米结构材料导热影响的研究还很少 见到报道。 ． 本文采用非平衡分子动力学模拟(NEMD)方法研究了氩模型二维纳米梯形变截面 结构材料的导热性质，意在研究热流密度和截面变化对纳米结构材料导热性质的影响规 律并对动能效应提供模拟验证。 2分子动力学模拟方法 本文模拟的系统示意图如图l所示。系统为固态氩原子按N[1111品格结构排列的周 期性梯形结构，最窄处的宽度为4．8nm，最宽处的宽度为14．4nm，模拟系统的总长度为 164．2nm，因此每一个梯形结构周期边界的梯度约为O．18。模拟中系统中间和两侧形成 不同的温度兀和死，沿热流方向实施周期性边界条件。根据乃和乃的相对人小，模拟 系统内会形成不同的热流方向以研究其影响。为消除边界条件带米的影响，模拟以图中 远离边界的A和B两个梯形结构单元为主要研究对象。 原子间的势能作用采用Lennard—Jonesl2-6势能函数，即 42 M=4cF咿2一(却 10引J，o=3．405x 对于氩原子之间的作用，取t=1．607x 10舢m。模拟计算中，能量和K 度分别以￡和盯为基本单位，温度的单位为r=120K。系统原予相空间的演化采用 K。计算前 Leapfrog—Verlet算法，系统的平均温度通过阻尼力约束方法控制为弘12 200000时间步为系统平衡时间，统计结果不少于500000时间步。 T． T 图1模拟系统示意图 t3l。在该方案中，系统在各边界可以采用周期性边界条件，将模拟系统沿热流方向划分 为刀层，每隔一定的时间选取第l层中速率最大(或最小)的原子和第力／2+1层中速率 最小(或最大)的原子交换速度矢量，在图l中相当于中间区域和两侧区域的能量流E， 也就意味着构造了高温端向低温端的热流，达到稳态时可以获得系统的温度分布。在计 算中可以通过调节速率交换的频率改变热流的大小，交换频率越大则热流越火。由于原 子间交换的是速度矢量，所以该方案具有自动保证系统总动量、总动能和总能量守恒的 优点。 3结果及讨论