纳米流体导热系数的分子动力学研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093359 纳米流体导热系数的分子动力学研究 李凌 杨茉 卢玫 余敏 上海理工大学能动学院 上海市军工路516 号 200093 021 liling@usst.edu.cn 摘要：论文主要应用平衡分子动力学方法研究了纳米流体中纳米颗粒表面附近的液体吸附层的特性， 并分析了该吸附层对纳米流体导热性能的影响。通过对固体纳米颗粒以及液氩分子运动轨迹的跟踪， 发现在纳米颗粒表面和液体界面处由于纳米颗粒较强的吸引力形成了一个液体吸附薄层，这个薄层内 分子的排列结构不同于基础流体，呈现出类似于固体的特性，因而使得纳米流体的导热物性比基础流 体有了显著提高。 关键词：平衡分子动力学 纳米流体 液体薄层 1．前言 很早之前研究人员就发现相对于构成纳米流体的基础流体而言，纳米流体的导热性 能得到了很大程度的提高。例如：Xie 等[1]在研究泵油中悬浮铝纳米颗粒构成的油－铝 纳米流体的特性时发现，当铝纳米颗粒的体积百分比为 5％，单位换热面积为 25 m2/g 时，纳米流体的导热系数比泵油提高了38％。Das 等[2]发现当水中悬浮有三氧化二铝或 者一氧化铜纳米颗粒时，在很小的一个温度范围内纳米流体的导热性能提高了两到四 倍。然而让科学家们感到困惑的却是纳米流体导热性能的显著提高无法用现有的理论进 行解释，虽然有大量的文献[3-4]对此问题进行了实验和数值模拟研究，并且都提出了一 些可能的原因，但依然没有一个确定的解释。 有些研究人员认为由于纳米颗粒分子之间较强的吸引作用，基础流体的分子会在纳米颗 粒的表面上形成一层薄层，这个薄层里面的原子结构比宏观液体里面的原子的结构要更 加的有序。由于固体内部的导热主要是靠声子的传播，那原子排列更加有序的液体薄层 的导热系数就会比较高，因而这个液体薄层对于纳米流体导热性能的显著提高具有很大 的贡献。Yu 和Choi [5] 发展了一个改良之后的迈克斯韦模型来说明纳米流体的导热机理， 他们认定虽然这个液体薄层可能只有几个纳米厚，但是这个液体薄层的存在却非常可观 地提高了纳米流体的有效体积比从而提高了纳米流体的导热性能，尤其是当这个薄层的 厚度小于10 纳米的时候。Leong 等[6]提出了一个模型来模拟纳米流体的导热系数。他们 认为纳米流体中固体颗粒和液体分界面之间的液体层以及纳米颗粒的尺寸大小是纳米 流体导热性能显著提高的主要机理原因，他们在研究的时候假定了液体层 的厚度是 1 纳米。Xie 等[7]也研究了这个液体薄层对纳米流体导热性能的作用，他们分别研究了液 体层的厚度、纳米颗粒的尺寸、纳米颗粒的体积百分比以及纳米颗粒和基础流体导热系 数之比对纳米流体导热性能的作用。 资助项目：国家自然科学基金，国家自然基金海外合作项目，上海市自然科 学基金（09ZR1422400），上海理工大学博士启动基金 上述这些模型从不同的角度证明了液体层对纳米流体导热性能的提高作用，但是它 们的缺点就是无法测量出液体层的厚度，因而只能是通过寻找适合实验数据的方法来获 得液体层厚度。论文将通过跟踪纳米颗粒及其附近液体分子的运动状况来分析液体薄层 的形成、特性，并对其厚度将进行定量计算，进而探讨液体薄层对纳米流体的导热物性 显著提高的机理。 2．基本模型 应用分子动力学方法进行模拟研究的时候分子之间势函数的正确确定是关键，多原 子分子之间的原子间作用复杂并且计算时间长，因此目前多采用单元子工质进行研究。 氩由于适用于普遍应用的较简单的 Lennard-Jones 12-6 势函数而被广泛采用研究对象 来应用分子动力学方法研究工质的特性[8]。本章也将应用液氩作为基础流体研究纳米流 体导热物性的提高。在进行模拟的过程中，每一个原子都要受到它周围的每一个原子的 Lennard-Jones 作用力，其受到的合外力应该满足牛顿定理： d 2 r mi 2i ∑Fij dt j ≠i 其中m 和r 分别是第 i 个原子的质量和位置，F 是第 i 个原