非傅立叶导热的分子动力学模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术年会 073278 非傅立叶导热的分子动力学模拟研究木 刘其鑫姜培学 向恒 热科学与动力工程教育部重点实验室，清华大学热能工程系， 北京100084 联系电话： (010 摘要：当系统的尺寸为纳米量级时，其导热系数呈现强烈的尺寸效应，而纳米尺度固体中的导热过程 也可能偏离傅立叶定律。在纳米尺度下，非傅立叶现象的实验研究比较困难。本文采用分子动力学模 拟的方法研究了纳米氩薄膜一端受热后的非傅立叶导热过程。通过对比MD模拟结果与cFD的计算结 果发现纳米氩薄膜的导热过程呈现非傅立叶导热的特点。 关键词：氩薄膜，分子动力学模拟， 非傅立叶导热 一 1．引言 、 目前，新技术的发展正朝着微彪内尺度的方向发展，在此尺度下，能量的传递发生在 原子量级，理解此尺度下能量传递的物理规律具有十分重要的意义。因此，最近一二十 年内，微／纳尺度传热成为了国际传热界的一大热点f】】【21。当系统的尺寸为纳米量级时， 固体的导热过程将偏离傅立叶定律‘”。傅立叶定律为： 2一兄·V丁 gF (1) and 此定律隐含了热量传递速度为无限大的前提。Cananeo 能量传递的速度有限性的模型： dq—qF—q dt tr (2) 这里，乃是非傅立叶过程的松驰时间。将(2)式代入能量守恒方程就得到了(3)式： a丁 a2丁 a2丁 百竹”瓦卸瓦C)f G}‘f ci‘x f气) ．(3) a是热扩散率。Tzou【61修改傅立叶定律，修改后的傅立叶定律为： g(，+丁。，r)=一尼V丁(r+丁，，，) (4) 式中f。和f，分别为热流和温度相延迟时间。在此基础上结合能掌守恒提出了双相滞后 模型： ’本文得到国家自然科学基金和清华大学基础研究塔金(JCpy2005049)资助 50 瓦+o矿硝l旷H一‘一l (5) 孙窘斗n一掣] 当f，．为零时，(5)式就变为(3)式了。 Y 分子动力学方法已经广泛应用于纳米尺度热物理的研究。分子动力学模拟(~∞)方 法将组成系统的微观粒子(原子／分子)视为经典粒子，将所研究的系统视为经典多体系 统，即满足牛顿运动定律。通过数值求解每个微观粒子的牛顿运动定律方程。记录下在 各时刻系统中原子／分子的位置和动量，从而模拟系统随时间演进的微观过程，用统计方 法计算系统的各种参数和输运性质。文献【7H9】报道了用分子动力学方法模拟研究非傅立 叶导热现象，但这些文献都集中在碳纳米管、硅等材料，这些材料中能量的传递不全是 靠声子进行的，本文将采用分子动力学方法研究固体氩热传导的非傅立叶现象。选取氩 为研究对象是因为在固体氩中，热能的传递全靠声子，并且固体氩的势能模型比较准确。 2．计算模型