碳管阵列湿润与滑移性质的分子动力学模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083403 碳管阵列湿润与滑移性质的分子动力学模拟 顾骁坤，陈民 （清华大学工程力学系，北京 100084） (Tel: 010 Email: mchen@tsinghua.edu.cn) 摘要：采用分子动力学方法研究了水在碳管阵列表面上的湿润和滑移现象。模拟获得了液滴在经过表 面处理的碳管阵列上的接触角和液滴展布特性。发现碳管阵列结构能改变纳米液滴在固体表面的展布 性质，使得疏水性表面的接触角大大增加。另外，通过模拟碳管阵列平板的 Couette 剪切流动，得到 了通道中流体的速度分布，观察到了速度滑移现象。模拟表明碳管阵列的排布能够影响滑移长度，进 而有可能改变和控制微纳米通道内的流体输运。 关键词：碳管 接触角 滑移长度 分子动力学 1 前言 在过去的 10 年中，超疏水表面在微流动、表面清洁等中的潜在应用引起了人们的 [1] 广泛兴趣。特别是在 Barthlott 等人系统地研究了“荷叶效应”的微观机理之后，超 疏水表面（接触角大于150°）的设计和制备技术得到了很大的发展。 碳纳米管由于其独特的结构使得它在力、热、光、电等多方面表现出优异的特性， 也使得它的制备和性能研究成为持久的热点。最近 Liu[2]等人报道了他们成功地制备出 可控管径、间距和层数，且高度规则排列的碳纳米管阵列。它天然的微观结构可能使得 经过表面处理的碳管阵列成为一种理想的超疏水表面和微流动通道壁面[3]。通常认为超 [4] 疏水表面能够减小流动阻力，原因是在疏水表面会产生较大的边界滑移。Truesdell 等 人的实验观察到牛顿流体在聚二甲硅氧烷表面上滑移长度可以达到250微米。此外分子 动力学模拟也被用来研究纳米尺度的流动滑移[5]，然而模拟所得到的滑移长度在纳米量 级，与实验结果之间有着3-4个数量级的差别，这种滑移长度的巨大差别主要来自于实 验和模拟中的纳米结构在尺度上有显著差别，但同时也对分子动力学模拟提出了挑战。 本文采用分子动力学方法研究了纳米液滴在碳管阵列上的湿润性质，以及水在碳管 阵列通道中的流动情况，目的在于：1）利用分子动力学预测液滴不同管径，不同间距碳 管阵列上的接触角；2 ）模拟微流动现象，了解影响滑移长度的因素。 2 分子动力学模拟 本文的分子动力学模拟分为两部分，液滴在碳管阵列上的湿润现象和水在碳管阵列 通道里的流动。模拟系统分别如图 1 和图 2 所示。模拟的壁面由两部分组成：壁面的上 部是带帽的碳纳米管阵列；壁面的底部是一层按照石墨六边形结构排列的碳原子。为了 使固壁保持静止或是以恒定速度运动，碳管的底层碳原子（“固定的碳原子”）所受的外 力被人为地设成 0 。 [6] 水分子之间的相互作用采用半经验三点势能函数 TIP3P 描述 ；因为碳管上的碳原 [7] 子不带电，所以碳管和水分子之间的相互作用只有碳-氧原子之间的范德华相互作用 ： σ σ CO 12 CO 6 U ( r) 4cε [( ) =−( ) ] CO