纳米流体中液体分子吸附与分子动力学模拟.pdfVIP

中国：I：程热物理学会 ：I：程热力学与能源利用 学术会议论文 编号：071123 纳米流体中液体分子吸附的 分子动力学模拟 刘娟芳，高虹，蔡治勇，曾丹苓 重庆大学动力工程学院(400044) ll2l 84 Vacolj鼠孕yahoo．com．cn，023巧5 摘要：本文用分予动力学模拟方法在分析纳米颗粒表面外流体的密度分布和径向分布函数的基础上， 对不同尺寸的铜纳米颗粒对液体氮分子的吸附作用进行研究。模拟结果表明，液氮分子与铜纳米颗 粒之问的相互作用能的绝对值随纳米粒子尺寸的增大而减小，使得在较小的纳米颗粒与液体分子之 间存在着较强的吸附作用。因此，吸附量随纳米颗粒减小而增大。 关键词：纳米流体，分子动力学模拟，吸附常量，径向分布函数 1、引言 近年来，由于纳米技术的迅速发展，纳米流体被用于越来越广泛的领域。所谓的 纳米流体是将直径小于100nm的纳米级金属或非金属氧化物粒子以一定的方式和比例 均匀地分散于流体中所得到的不同于纯液体性质的流体。最近有关实验研究结果表明： 在液体中加入纳米粒子可明显地提高液体的导热系数和对流换热系数。另外的一些研 究结果表明：纳米颗粒对流体分子的表面吸附不仅影响纳米颗粒的动力学行为，而且 影响纳米颗粒表面所吸附的液体分子的结构。但用实验方法对吸附机理进行直接的研 究非常困难，而数值方法一般可对吸附现象进行较为细致的研究。随着计算机和计算 技术的迅猛发展，分子动力学模拟方法成为一种从原子尺度上研究吸附现象及机理的 有力工具，它通过使用一些经验或半经验势函数来描述原子间的相互作用。本文主要 是利用分子动力学模拟方法对纳米流体中液体分子的吸附机理进行研究。 2、模拟方法 在用原子描述相互作用的基础上，利用标准的分子动力学(MD)模拟技术，对单 个纳米铜同体颗粒悬浮于液体氮中所形成的纳米流体的动力学性质进行研究。系统中 原子间相互作用势均采川Lennard-Jones势函数进行描述，其形式如一卜．： ㈨=4％例2一删 ㈩ 式中，r为相互作川两个原子间的距离；占鲫为最小势能，对应于原子间吸引的强度： 109 仃筇为口、∥原子排斥核重叠时的距离。对于不同种类的原子间相互作j_}j的LJ使能 i+-一1 ㈤一 口筇=等产．j 而占鲫则以JL何平均近似 s。，=际 (3) 文中各原子的LJ参数见表1。本研究中，对金属纳米铜粒子并未采用嵌入势函数、 或紧束势函数，是因为LJ势函数较其简单可靠。另外也是因为本文只是对吸附机理进 行定性的研究。 表1分子动力学模拟中符原予的LJ势能参数 ‘图1为铜纳米颗粒与液体氮所形成 的纳米流体示意图。在模拟盒的三个方 向均使用了周期性边界条件，为了避免 模拟盒和像盒子中的铜纳米颗粒作用， 模拟盒的边K需设定得足够大。液体氮 分子的密度取为O．70039／cm3。在模拟 中，采J4j速度修正方法将液体氮分子利 俐体颗粒原子的温度修正到平衡温度。 牛顿运动方程采用Beeman’s Leapfrog算 法进行积分，其中的时间步艮取为r|毛 秒。初始时刻，同体颗粒位于模拟盒的 中心，总动量为0。 图l模拟模型的示意图 3、模拟结果与讨论 3．1密度分布 图2示出了液体氮分子在铜纳米颗粒周同的局域密度分布，其中横坐标为液体氮 分子离铜纳米颗粒表面的距离，而纵坐标为