微通道中Poiseuille流动的分子动力学研究.pdf

微通道中Poiseuille 流动的分子动力学研究1 刘彬武，刘朝 重庆大学动力工程学院，重庆（400044 ） E-mail ：liuchao@ 摘 要：本文建立了外力作用下氩流体在平行平板间流动模型。以此模型为基础，模拟计算 了氩在不同性质壁面（亲水性和憎水性）间的流动过程，并对速度分布、温度分布等进行了 统计。模拟结果与N-S 方程和能量方程的解进行了比较,发现在憎水性壁面附近，易观察到 速度滑移和温度阶跃现象，而亲水性壁面上则反之。亲水性通道速度和温度的模拟结果与分 析解吻合较好,憎水性通道的模拟结果与分析解相差较大。同时发现，在系统的能量变化过 程中，憎水性系统的能量明显高于亲水性系统，这可以由亲水性壁面传热性能较好给出合理 解释。 关键词：分子动力学，Poiseuille 流动，亲水性，憎水性 1. 引言 随着精密加工技术和微电子技术的发展，出现了越来越多的微电子机械系统(MEMS)器 件。由于器件尺寸的减少，形成了许多微结构,这些微结构中流体的物理特性将会对器件的 性能造成显著的影响[1,2] 。所以对微结构中的微观非经典问题的研究尤为重要。目前，对于 微结构中流体性质的研究仍处于探索阶段。许多问题的研究在现有的试验环境下难以完成。 要对微尺度条件下的流体进行试验来得到相应流体的性质在目前的技术水平还是有一定的 困难性，这样我们必须借助其它的方法对微尺度下流体的性质进行科学的计算和分析。随着 计算机技术的发展，利用计算机来模拟流体的流动成为可能，为我们研究微观流动问题提供 新的研究手段。 2. 模拟细节 2.1 分子动力学研究方法 本文模拟了不可压缩流体氩在外力驱动下的Poiseuille 流动，固壁粒子为Pt 原子。如图 1 所示。在X,Y 方向采用周期边界条件，Z 方向为固壁边界条件。流体粒子和固体粒子均为 864 个。X 和Y 方向长度为无因次的38.669，流体在Z 方向的长度为19.334， 图1 微通道模型 壁厚为9.667 。外力作用于X 方向驱动流体流动。外力作用产生的功转化为两部分，一 部分为动能驱动流体流动，另一部分转化成热使流体升温。产生的部分热量必须通过壁面导 1本课题得到国家自然科学基金（No）的资助。 - 1 - 出，以阻止流体温度的无限制升高，处理方法是根据给定温度(T=1.0 ）调节壁面上反射粒子 的速度，相当于在壁面上固定一个恒温器[3] 。 流体粒子间的相互作用采用L-J 型双体势能函数： ⎡ 12 6 ⎤ u r 4=ε σ/ r − σ/ r ( ) ( ) ( ) ij ⎢⎣ ij ij ⎦⎥ 其中，rij 为分子间距，氩分子的基本参数σ =0.3410nm 、ε=1.6154×10-21J 。 流体粒子和固壁粒子间的相互作用采用如下势能函数： ⎡⎛ ⎞12 ⎛ ⎞6 ⎤ σ σ φ (r ) 4ε ⎢ sl =−β sl ⎥ sl ij sl ⎜ ⎟ ⎜ ⎟ ⎜