利用分子动力学模拟对纳米流体分散稳定性的研究初探.pdfVIP

中国工程热物理学会 工程热力学与能源利用 学术会议论文 编号：101200 利用分子动力学模拟对纳米流体分散稳定 性的研究初探 1 王楠 陈俊 安青松 史琳 ( 清华大学热科学与动力工程教育部重点实验室 北 京100084 ) Tel:010E-mail: anqingsong@ 摘 要 纳米流体的制备和分散稳定性是所有后续纳米流体研究工作的前提。本文利用平衡分子动力学 模拟方法，尝试从微观角度研究纳米流体的分散稳定性。以纳米颗粒 Cu 在液态 CO2 中的分散稳定为 算例，通过跟踪每个固体分子和液体分子的位置，分别观察了颗粒周围的液体吸附现象和颗粒的团聚 现象。通过改变固液分子间Lennard-Jones 势能函数中引力项与斥力项的系数，观察其对颗粒周围液体 分子的密度分布和颗粒团聚的影响。模拟结果显示，颗粒周围存在液体吸附层；固液间引力项增强， 可以提高颗粒周围吸附液体分子的数目，但不改变液体层厚度；固液间引力项增强还可以有效抑制或 延迟团聚的发生。这为纳米颗粒的分散稳定提供了新的思路。 关键词：分子动力学模拟，纳米流体，分散稳定，固液作用，团聚 1 引言 1995 年美国Argonne 国家实验室的Choi[1]等人首次提出了纳米流体这一概念，即将 纳米尺度的金属或金属氧化物粒子以一定的方式和比例悬浮于液体中，形成一类新的传 热冷却工质。很多实验研究结果表明，纳米流体的物性和传热性能相对基液会有所改变， 具有广阔的应用前景，成为了研究热点[2-6] 。 纳米颗粒在制冷剂中的分散稳定性研究是后续所有研究的基础性工作，目前大多数 学者从实验和宏观受力方面进行研究[7-10] ，包括不同因素对分散稳定的影响情况，分散 剂的作用等，但少有从微观角度研究颗粒在基液中的分散稳定性。本文以纳米颗粒 Cu 在液态CO2 中的分散稳定为对象，利用MD 模拟颗粒在基液中的运动情况，主要关注颗 粒表面吸附液体层的情况和颗粒的团聚现象，并且通过改变固液分子间 Lennard-Jones (L-J)势能函数中引力项和斥力项的系数，观察其对吸附层和颗粒团聚的影响。希望为纳 米流体稳定性的研究探索一种新的方法。 2 模拟方法 分子动力学模拟通过求解系统中的每个分子在其它所有分子作用下的牛顿运动方 程，获得每个分子的位置和速度，从而得到系统内粒子随时间的运动过程，最后统计出 宏观参量。本文模拟Cu 纳米颗粒悬浮于液态CO2 中形成的纳米流体。其中粒子的初始 排布满足面心立方晶格（FCC ），初始速度由波尔兹曼分布随机获得。分子间作用势采 1基金项目：国家973 资助项目（Grant No. 2010CB227305 ），国家自然科学基金资助项目（No） 用L-J 势能函数，其形式如下：  12 6      σ σ U( r ) 4ε −  (1) ij          r r  ij ij      ε σ 其中, r 为粒子i 和j 之间的距离， 和 分别为L-J 势能函数的能量参数和长度参 ij