纳米通道内液态微流动密度分布特性数值模拟研究-物理学报.PDF

物理学报 Acta Phys. Sin. Vol. 62, No. 12 (2013) 124705 纳米通道内液态微流动密度分布特性数值模拟研究* † 胡海豹 鲍路瑶 黄苏和 ( 西北工业大学航海学院, 西安 710072 ) ( 2013年1月17 日收到; 2013年2月27 日收到修改稿) 微通道内流动因表面积/体积比值极大, 造成许多微尺度效应, 进而使微通道内出现完全不同于宏观流动的流体 密度分布特性. 本文以纳米通道内液态Poiseuille 流为对象, 采用非平衡分子动力学模拟方法研究了流体原子间相 互作用强度LL , 流体原子间平衡距离LL 以及壁面原子与流体原子间平衡距离LS 对通道内流体密度分布的影响 规律. 数值模拟中, 统计系综取微正则系综, 势能函数选用LJ/126 模型, 壁面设为Rigid-atom 壁面, 温度校正使用速 度定标法, 牛顿运动方程的求解则采用Verlet 算法. 模拟结果表明, 随LL 的减弱, 近壁面区密度分布的振荡幅度则 逐渐增大; 而LL 则同时影响流体原子的存在形态和密度分布, 较大的LL 会造成流体原子在整个通道内呈现面心 立方结构的类似固体排列, 较小的LL 会使得流体原子呈现不断变化的“ 团簇” 结构; 随LS 的变大, 近壁面区流体 密度振荡幅度增大, 且流体密度分布起点离壁面越远. 另外, 本文还从近壁面区流体原子的“俘获- 逃逸” 行为角度, 初步解释了原子间相互作用强度对密度分布的影响规律. 关键词: 纳米通道, 微流动, 密度分布, 分子动力学 PACS: 47.61.k, 83.50.Ha, 83.10.Rs DOI: 10.7498/aps.62.124705 法对纳米尺度下液体的Couette 流和Poiseuille 流 1 引言 进行了模拟, 发现了壁面“滑移” 现象和密度振荡 现象的存在; Travis 和Gubbins4 采用分子动力学 近十几年来, 自然科学和工程技术发展的一个 模拟方法对微尺度下的Poiseuille 流进行了系列模 重要趋势是朝微型化和纳米化迈进, 因而迫切需要 拟研究, 并对微流动模拟结果与基于N-S 方程解 弄清微、纳米尺度条件下的流体流动规律. 在这样 算结果的差别进行了深入分析; Song 和Wang5 研 的微系统中, 其通道的水力直径极小, 表面积/体积 究了不同晶格平面种类等条件下的Couette 流动和 的比值极大, 从而引起了许多微尺度效应, 如表面 Poiseuille 流动, 给出了不同条件下固液作用强度 张力效应、低雷诺数效应、双电层效应、分子力 对滑移长度和流体密度分布的影响规律. Song 和 效应等, 这些因素的影响使得用于宏观流体运动分 Wang 6 采用非平衡分子动力学方法模拟了不同壁 析的Navier-Stokes (N-S) 方程不再适用. 面润湿性时纳米通道内的Poiseuille 流, 给出了壁 为探究该类微流动的流动特性, 近年来越来 面润湿性对流体流量、滑移速度和滑移长度的影 越多的国内外学者开始关注该交叉学科领域. 如 响规律, 并且指出壁面润湿性对流体密度分布也有 Thompson 等 12 研究了纳米通道内的Couette 流 着重要的影响. Xin 和Zhang 7 采用分子动力学方 动, 给出了固液相互作用强度对流动系统滑移长度