过渡区气体微尺度流动的格子Boltzmann模拟.pdfVIP

第 41卷 第 6 期 力 学 学 报 Vo1．41，NO．6 2009年 11月 ChineseJournalofTheoreticalandApplied M echanics NOV．2009 过渡区气体微尺度流动的格子 Boltzmann模拟 田智威 ) 郑楚光 t王小明 (中国地质大学 (武汉)资源学院，武汉 430074) 十(华中科技大学煤燃烧国家重点实验室，武汉 430074) 摘要 采用格子 Boltzmann方法 (LBM)对过渡区内的微尺度气体流动进行了模拟研究．在已有滑移区微流 动 LBM 模型中引入 Knudsen层速度修正，选取合适的修正函数表达式并依据动理论确定了可调参数的合理 取值．在边界条件的处理格式上，采用了适合过渡区模拟的高阶滑移边界的替代格式来捕捉过渡区微流动的滑 移速度，避免了直接求解高阶速度导数项的数值困难．通过对两类不同的微流动进行模拟的结果表明：与数值 解吻合得较好，尤其是对 KnO．5微流动滑移速度的预测，与已有 LBM 的模拟结果相比有明显的提高． 关键词 格子 Boltzmann方法，微尺度流动，过渡区，速度滑移，非线性压力分布 中图分类号： TK121，0354 文献标识码： A 文章编号： 0459—1879(2009)06—0828—07 引 言 方法模拟结果进行对 比，得出了 格“子Boltzmann模 型并不适合模拟过渡区微流动”的结论．但是需要指 近年来，微尺度下的流动特性研究已经成为了 出的是，其采用 的是 2002年 Nie提 出不完善的微流 国内外研究热点问题之一 _l12J．由于相关实验研究的 动LBM模型 8【】，该模型即使是对滑移区的模拟结果 限制、困难以及测量结果的不确定性，本文重点放 也很难令人满意． 在 了可信度更高的数值模拟的研究上．常规模拟方 随着模型的改进和发展，尤其是代表性文献 [6，7] 法包括了从连续介质模型方法到量子分子动力学方 的发表，极大地推动了格子Boltzmann模型应用于模 法一系列的有效方法．基于连续介质假设的模拟方 拟微流动的研究．在 2006年后，有少数研究者采用 法，在引入了速度滑移边界条件后，可以对微尺度流 LBM 来模拟过渡区的微流动：其中一类 [14,15】是基 动进行基本的模拟研究，但仅限于滑移区的模拟，对 于分子运动论，加入高阶的平衡态分布函数和高精度 于高 Kn数的微流动则不再适用．分子动力学方法 的Gauss—Hermite积分，构造出满足高阶矩条件的新 (MD)更适合于模拟液体或稠密气体流动，而不适合 模型，从理论上克服了格子Boltzmann模型的局限， 模拟稀薄气体流动和微尺度流动，而且其计算代价 以更精确地描述高 Kn数流动的特征．但是由于高阶 十分巨大，仅用于理论上的研究．直接 Monte—Carlo 矩的引入，使得模型的构造过于复杂，在一定程度上 模拟方法 (DSMC)在一定程度上解决了MD方法计 丧失了LBM 的一些固有优势，且不适合于实际模拟 算量大的问题，但却存在较大的统计噪声，尤其是 和应用， 目前仅仅停留在理论的分析和推导上． 对于一般的低速微流动． 另一类 [16,17J则是在 已有 LBM 模型中引入壁 而随着格子 Boltzmann方法 (1atticeBoltzmann 面修正函数，对 Knudsen层内的分子平均 自由程进 method，LBM)[ J这～新兴的计算流体力学方法的 行相应的修正，并提出了 有“效平均 自由程”的概 迅速发展和完善，已有相关的文献实现了LBM模拟