纳米通道内液体流动的分子动力学研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号；。73280 纳米通道内液体流动的分子动力学研究 向恒姜培学刘其鑫毛志方 热科学与动力工程教育部重点实验室，清华大学热能工程系，北京100084 联系电话：(010 摘要：本文采用分子动力学方法模拟液体在定截面及变截面纳米通道内的三维Poiseuille流动，研究液体 浸润性及通道截面变化对纳米通道内液体流动的影响。’研究结果表明液体对壁面不湿润时，壁面处有速度滑 移存在，并且速度滑移随接触角的增大瓶增火，间时，摩擦阻力系数减小；液体对壁面湿润，则不存在速度 滑移，同时，摩擦阻力系数增大。通道截面形状的变化对流场的影响随着远离变截面位置而迅速减弱，对流 体通道平直处的流体速度相对值分布影响不大，对摩擦阻力系数影响较大。 关键词：纳米通道液体Poiseuille流动速度滑移变截面分子动力学 1 引言 高度集成的需要和微加工技术的发展，将各领域里所使用材料的特征尺寸推进到了纳米 尺度。在微型器件封装的内部加入微型流道，将冷却流体直接引入到器件内部进行换热，是 解决器件发热密度日益增长问题的一个探索方向，这也使得微通道内流体流动及换热成为研 究热点。微通道内流体流动及换热的实验研究因为实验技术及测量上的难度，很难将所研究 通道的尺寸做到非常小。分子动力学模拟从原子／分子的尺度来研究问题，成为研究纳米尺度 通道内流体流动和换热的有力手段。 面粗糙度对气体在壁面处的速度滑移的影响。徐超等15】模拟了三维Couetce流，研究了壁面温 度及壁面作用势强弱对滑移速度的影响。本文采用非平衡分子动力学方法，通过模拟液体在 定截面及变截面绱米通道内的三维PoiseuiJle流动，研究液体湿滑陛及通道截面变化对纳米通 道内液体流动的影响。 2 模型及模拟方法 为了对比验证，本文先对液体在平直通道内PoiSeuille流动进行了模拟，模拟采用已有模 型n如图1所示。在由两平行平板相距一定距离构成的通道内，液体在x方向的压力驱动 下流动。模型在x和z方向上采用周期性边界条件来模拟无限大的情况。在已有模型的基础 上，本文还构造了变截面通道，对液体在变截面通道内的流动进行模拟。相关模型如图2所 示(均为模型在妒)，平面上的投影)。图中的灰色部分为通道的壁面，黑色部分为液体。与平直 通道模型一样，变截面通道模型在x和z方向上也使用了周期性边界条件。两个模型是在壁 本文得到国家自然科学疑会和清华火学皋础研究基金(JCpy2005049)资助 854 面上构造菱形凹陷或突起以实现变截面通道。 图1平直通道Poiseuille流动的模型 图2变截西通道P0jseuille流动的模型 本文关注的是纳米通道壁面对液体流动的影响，没有考虑换热，因此将流道壁面的分子 设置成固定状态，与液体分子有力的作用但没有热量交换。壁面分子与液体分子之间以及液 体分子之间的作用势均采用u12．6势能模型来描述，即 ‘ ≯(，)：{三c占Kcr7，)12一(盯厂厂)6】，，≤t (1) 【u，，．I 式中，势能常数占=1．67×10-21J，平衡常数仃=3．4×l仃20 m，截断半径t=2．5仃。式中的参 数c是为了调整分子之间的作用强弱，液体分子之间的作用势取％=1；而液体和壁面分子之 间作用势随c，，取值的不同而表现液体与壁面不同的浸润性质。因为在本文所有的模拟中c，， 均取1，为了简化起见，在本文后面部分中将c。，表述成c。 对于U12．6势能函数，液体与壁面的接触角可表示为【7J： cos侈=2占。，／qf一1 (2) 式中s。，为壁面与液体分子之间的势能参数，占，，为液体分子之间的势能参数。在本文中，此 式转化为：