水在纳米管道中流动行为的分子动力学模拟.pdfVIP

Vol_33 高等 学校 化 学 学报 No．12 2012年 12月 CHEMICALJOURNALOFCHINESEUNIVERSITIES 2727～2733 水在纳米管道 中流动行为的分子动力学模拟 饶 欢 ，杨文泓 ，王思邈 ，孔 滨 ，刘文志 ，李晓霞 ，杨小震 (1．北京分子科学国家实验室，高分子物理与化学国家重点实验室， 中国科学院化学研究所，北京 100190； 2．多相复杂系统国家重点实验室，中国科学院过程工程研究所，北京 100190) 摘要 采用分子动力学方法模拟了纳米管道内水分子的流动行为．考察了压强差、管道直径和管道长度对 通量的影响，验证了流体在直径为2～2．8nm管道中的流动行为符合 Hagen．Poiseuille(HP)方程．研究发现， 末端效应具有长度依赖性 ，对于较短的管道其末端效应更为显著．为了深入了解真实非均匀管道的流动行 为，使用直径不同的2种管道以不同的连接顺序组合成4种非均匀管道模型，最终得到了不同的流量．提出 了二元管道模型，对非均匀管道内的流动行为进行了初步探讨． 关键词 分子动力学模拟；纳米管道；流动行为 ；Hagen—Poiseuille方程；末端效应 中图分类号 O641 文献标识码 A doi：10．7503／cjc 纳米尺度流体的流动行为研究是近年的热点之一．虽然纳米流体流动规律的研究刚刚开始 ¨．2]， 但是纳米流体技术已经有较广泛的应用．例如用多孔膜进行分离和纯化l3J、分子检测 和生物膜离 子通道 等．Hagen-Poiseuille(HP)方程一般用于预测宏观尺度流体通过圆柱形管道的体积流动通 量 ．然而，有些实验发现纳米尺度流体的流动行为偏离HP理论的预测 ，加J．因此，宏观上影响流体 流动行为的因素在微观尺度下是否起到同样的影响仍需要深入研究．Zhu和Whitby等_1卜”研究发现， 压强差与通过纳米管道流体的通量呈线性关系，这与HP理论预测相符．纳米管道直径和长度对通量 的影响较为复杂，而 目前的研究主要集中在管径小于 1．66nm的纳米管道，这种尺度的纳米管中流体 流动行为不符合 HP理论的预测 引．这是由于直径小于 1．66nm的纳米管道内的水分子会出现特殊 的结构 (如n元环状结构 。。和单列水分子结构 )．Thomas等 发现当纳米管道直径在 1．66～ 4．99B1TI之间时，通量随着直径的增加而增加，与HP理论预测一致．但该研究并没有深入讨论通量与 管道长度的定量关系．此外，纳米管道孔壁与流体之间的相互作用也对通量产生影响[15,23,24]． 在实验使用的多孔膜中，管道沿轴向的管径并不均匀．Goldsmith等 初步探讨了圆锥形管道内 的流体流动行为．因此研究这种不均匀性对流体在管道内的流动行为的影响有着十分现实的意义． 本文主要利用分子动力学方法在原子级别上模拟研究了压力驱动水分子通过管径介于2～2．8nm 之间的无定形碳纳米管道 ¨的流动行为．通过压强差、管道直径和长度对通量的影响，考察 HP理论 预测微观尺度行为的有效性．在管径固定的情况下，通过改变管道长度观察纳米管道的末端效应，发 现未端效应对长度有依赖性．考察了4种管径不均匀的纳米管道中流体的流动行为，并讨论了造成不 均匀管径的流体行为差别的顺序特征． 1 分子动力学模拟 1．1 模 型 实验所用的过滤膜是一种多孔介质．通过给活塞施加恒定的压力，驱动流体形成跨膜流动．模拟 实验水的过滤过程如图1所示模型．图1(A)是沿管道轴向的剖面图，体系由前后2个水箱、中间的无 收稿 日期：2012-03-22． 基金项 目：中国科学院创新项 目(批准号：KJCX2Yw-H．19)和国家 自然科学基金 (批准号资助． 联系人简介：杨小震，男，研究员 ，博士生导师，主要从事高分子计算机模拟研究．E—mail：yangx@iccas．ac．cn 孔 滨，男，博士，副研究员，主要从事微观及介观计算机模拟研究．E-mail：lbkong@iccas．ac．cn 饶 欢等：水在纳米管道中流