微圆柱气体绕流流体力学分析.PDF

LHD 2011 年度夏季学术研讨会 401 微圆柱气体绕流的流体力学分析 胡远，孙泉华，樊菁 中国科学院力学研究所高温气体动力学国家重点实验室 （筹），北京海淀区 100190 摘要 本文采用基于分子动理论的耦合算法和基于连续流假设的传统CFD 方法对空气绕微尺度圆柱的 流动进行了数值模拟。Re 数的变化从0.01 到40，Ma 数则从不可压缩变化到 1。通过包含不同物理 机理的多层次的数值模拟，分析了气体绕微尺度圆柱流动的流体力学行为，特别是流动的可压缩效 应和稀薄效应对流动的影响。研究发现，流动可压缩性受低 Re 数效应的作用而增强，而稀薄效应 在一定程度上削弱流动的可压缩性。可压缩效应的影响使分离涡的长度随 Ma 数的增长呈现先增后 减的非单调变化趋势；而稀薄效应使分离显著推迟，分离涡缩小甚至消失。由于稀薄效应的影响， 在圆柱表面可以观察到速度滑移以及非平衡正应力，减小了表面所受的剪切应力。总的来说，可压 缩效应增大圆柱的阻力系数，而稀薄效应则相反，两种效应对圆柱阻力的变化贡献存在竞争关系。 关键词 圆柱绕流，低Re 数，可压缩性，稀薄效应 试 [4]。在蠕流流动中，惯性力的作用远小于粘 0B 引 言 性力，因而Stokes 将Naiver-Stokes （NS ）方程 圆柱绕流是流体力学研究的一个经典问题， 中的非线性惯性项全部忽略，得到了线性化的 人们对该问题的关注已经有一个多世纪了 [1]。 Stokes 近似。虽然 Stokes 给出了该近似的三维 在最近几十年微机电系统（MEMS ）迅速发展 球对称解，但遗憾的是，在求解二维圆柱问题 和广泛应用的刺激下，这些经典的流体力学问 时，Stokes 发现不存在同时满足物面条件和自 题受到了新的关注。这是因为在微尺度下，流 由来流条件的自洽解。这也是著名的 Stokes 佯 动表现出许多非同一般的特性 [2]。正因如此， 谬。Oseen 在 1910 年指出被 Stokes 近似忽略的 在微机电设备中的流体力学问题也成为了新的 惯性项会在远场起主导，并在近似中部分考虑 研究热点[3]。 惯性项的作用，得到所谓的 Oseen 近似，从而 微流动因其流动特征长度非常小，通常属 使 Stokes 佯谬得以解决。Oseen 近似方程的求 于低雷诺（Re ）数流动的范畴。此外，由分子 解比较复杂，Lamb 第一个给出其二维圆柱问 平均自由程（ ）与流动特征长度的比值定义 题的近似解 [5] 。实际上，蠕流问题是一种奇异 了努森（Kn ）数，它表征了流动非连续或者稀 摄动问题。对奇异摄动问题，一种比较有意义 薄程度的大小。标准状况下，空气的 大约为 的求解方法是匹配渐进展开法。Proudman 和 65nm 。对于一般的流动来说，Kn 数的影响可 Pearson [6]，以及Kaplun [7] 是较早采用匹配渐 以忽略，流动被认为是连续的。但是在微尺度 进展开方法来求解极低 Re 数圆柱绕流问题的 流动中，Kn 数的值并不小到可以被忽略，流动 学者。其中，Kaplun [7]给出了在壁面附近的三 开始表现出一些稀薄效应。事实上，Kn 数可以 阶近似解。当Re 大于1 以后，蠕流流动的假设 表示为两个熟悉的无量纲数——马赫数（Ma ） 不再成立，这些分析解也随之失效，实验研究 和 Re 数的函数，即 。Re ， 和数值求解 NS 方程开始成为主要研究手段。 Ma 和Kn 数分别代表着粘性效应，可压缩效应， 在实验研究方面，Tritton 开展