纳米孔内气体流动的MD模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术年会 编号：083356 纳米孔内气体流动的MD模拟研究* 刘其鑫 姜培学 向恒 热科学与动力工程教育部重点实验室，清华大学热能工程系， 北京 100084 联系电话： （010 E-mail: Jiangpx@mail.tsinghua.edu.cn 摘要：了解和掌握微纳尺度流动规律对 NEMS 等设计制造具有重要意义。当系统的特征尺寸与流体 的分子平均自由程相当时，气体流动表现出一些与宏观尺度下很不相同的特性。目前，实验研究只能 测得流动的宏观参数，对流动细节不能获得。本文对氩气在纳米孔内的流动进行分子动力学模拟研究， 得出气体在纳米孔内的速度分布和阻力系数。 关键词：纳米孔，努生数， 流动阻力系数 1. 引言 目前，新技术的发展正朝着微/纳尺度的方向发展，在此尺度下，了解和掌握微纳尺 度流动规律对NEMS等设计制造具有重要意义。因此，微/纳尺度流动的研究成为一大热 点。目前，已经有学者对纳米尺度下气体流动开展实验研究，但实验研究只能测得流动 的宏观参数，不能获得流动细节如速度分布，而分子动力学模拟能弥补实验这方面的不 足。Travis等用非平衡分子动力学模拟方法模拟了单原子流体及多原子流体[1,2]、L-J流体 [3]在间距为几个纳米两平板间的Poiseuille流，发现当平板间大于 5 个分子间距时，流体 在平板间流动时的速度分布除了在边界有速度滑移外，其它部分与N-S方程预测的相同。 Koplik [4]和Todd[5]采用非平衡分子动力学模拟了流体在两平板间的低雷诺数流动，发现 在流体与平板的界面存在速度滑移，其它部分的速度分布与N-S方程描述的一致。对于 气体在纳米孔内的Poiseuille流的非平衡分子动力学模拟，Bhattacharya等[6]用非平衡分子 动力学模拟了纳米尺度下努森数为 0.27 的气体流动，发现N-S方程在努森数为 0.27 时似 乎仍然有效。本文采用非平衡分子动力学模拟方法模拟氩气在纳米孔内的压差驱动下的 流动。 2.物理模型和模拟方法 本文选择氩气作为气体研究对象。由于关心的是气体的流动，在模拟中固体部分的 原子固定不动，选择氩固体代表固体通道，气固之间的作用通过势能参数调节。氩气分 子之间以及氩气分子和固体之间的势能作用均采用 LJ12-6 势能函数。 δ ⎡ δ 12 6 ⎤ ⎛ ⎞ ⎛ ⎞ U r () 4ε⎢⎜ ⎟ −⎜ ⎟ ⎥ （1） r r ⎢⎣⎝ ⎠ ⎝ ⎠ ⎦⎥ 其中能量参数ε和距离参数δ 取值为：氩气分子之间取ε=1.67 ×10-21 J，δ =3.405 -10 m ；对于氩气分子和固体之间的作用，δ 的取值与氩气分子之间相同，取ε=2.82 ×10 * 本文得到国家自然科学基金和清华大学基础研究基金(JCpy2005049)资助 -21 ×c ×10 J ，其中系数c在 0 到 1 之间，以调节气固之间作用强弱。模拟的几何模型如 图6.15 所示，图中外围的圆形代表纳米管的横截面壁面，气体分子在压差作用下在纳米 管内形成沿x方向的宏观流动。 图 1 气体在纳米管内流动的分子动力学模拟几何模型 本文采用非平衡分子动力学模拟方法模拟气体在纳米管内的Poiseuille流动，对纳米 dP dx g / −ρ x [1,2,6] 管中的氩气分子全部施加x方向的加速度场gx ，这等效于压力梯度