微尺度流动的扩展二阶滑移边界条件研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083200 微尺度流动的扩展二阶滑移边界条件 1* 1 2 郭照立 ，郑楚光 ，施保昌 1 2 ( 华中科技大学煤燃烧国家重点实验室，华中科技大学数学系，武汉 430074) 联系电话：027 E-mail:zlguo@hust.edu.cn 摘要：本文提出了求解微流动Navier-Stokes 方程的扩展二阶滑移边界条件。该边界条件依赖一个有效的分 子平均自由程，并考虑了气体分子与壁面碰撞的影响，能够刻画Knudsen 层内的流动。基于该边界条件， 对微通道内的气体流动进行了研究，计算结果与实验数据和其他理论结果相吻合。 关键词：微尺度流动、边界条件、分子平均自由程 0 前言 当气体流经固体壁面时，在壁面附近存在一个分子边界层，称为 Knudsen 层(Knudsen Layer, KL) 。KL 不同与一般流动的粘性边界层，其厚度为分子平均自由程(mean-free-path, MFP) 的量级。在 KL 内，气体分子之间的碰撞非常稀少，局部热力学平衡假设不成立， Navier-Stokes 本构关系不再适用。对一般的大尺度连续流动，由于KL 层的大小可以忽略， Navier-Stokes 方程及相应的无滑移边界条件可以用于描述整个区域内的流动。但是，对微尺 度流动，KL 层在整个流场中的比例较大，因此对其中的流动进行准确描述起着关键的作用。 为了描述KL 层内的流动，人们提出了多种超越Navier-Stokes 的流动模型，如Burnett 方程、super-Burnett 方程、Grad 13 矩方程、正则化Grad 矩方程等[1-7] 。但是，这些模型用 于微流动时存在各自的困难。例如，如何为这些模型配置适定的边界条件至今都是尚未解决 的问题，并且这些方程往往由于内在的不稳定性而难以计算求解。 由于这些困难，近年人们又开始探索在 Navier-Stokes 模型的框架内研究非平衡流动的 新途径[8-15] 。最近文献[15]提出了一个考虑KL 影响的扩展Navier-Stokes 本构关系，初步 计算结果表明该模型可以捕捉KL 流动。但是，通常Navier-Stokes 方程使用的边界条件不适 用于这一模型，需要研究更合理的边界条件。本文将就这一问题开展研究。 1 依赖流场几何形状的有效粘性 周知，流体的粘性是表征流体分子碰撞过程中动量交换的物理特性。分子连续发生两次 碰撞之间的分子自由飞行距离与分子的动量交换有重要关系。从分子动理学理论可知，无壁 面约束的分子平均自由程  与粘性系数  有如下关系[16]：    RT  , (1) p 2 其中T 为温度，R 为气体常数，p 为压力。当有壁面存在时，在远壁面区域上述关系成立； 但在KL 内，由于气体分子的飞行距离受壁面的影响非常大，上述关系不再成立。 壁面对平均自由程的影响可以通过考察分子自由程的概率分布来进行分析。由分子动理 r /  学理论可知，分子自由程满足一个指数分布，即p (r) e  /  。对无壁面约束的系统，  资助项目：国家自然科学基金 （No.50606012