微尺度后台阶流动的Monte Carlo直接模拟.pdfVIP

计 算 物 理 第23卷第2期 V01．23．No．2 CHINESE OF 2006年3月 JOURNALCOMPUTATIONALPHYSICS Mar．．2006 [文章编号】1001．246X(2006}02—0171—06 微尺度后台阶流动的MonteCarlo直接模拟 秦丰华， 孙德军， 尹协远 (中国科技大学力学与机械工程系，安徽合肥230027) [摘 要] 对DSMC方法中压力边界条件的实现给出了一种改进的方法，并分区模拟了微尺度后台阶气体流动， 得到了其流动的一般典型结构．计算结果表明经典后台阶流动中的回流区在一定参数范围内仍然存在，但上壁 面的分离区消失．讨论了流动参数对流动结构的影响，并给出了再附长度与流动参数的关系． [关键词】DSMC；压力边界条件；分区；微尺度后台阶流动；再附长度 [中图分类号】035 [文献标识码】 A O 引言 ElectroMechanical MEMS(Micro 流动由于空间的细化而有其独特的性质，其中最重要的有两个方面：特征尺度的减小导致较大的Kn，使流动 处于过渡领域，稀薄气体效应明显…；较小的横向尺度导致了较强的粘性耗散以及气体密度的剧烈变化，使 流动在低马赫数下呈现出显著的可压缩效应心J． 一般认为，对于过渡领域的流动，N—S方程已不再适用，控制方程为高度非线性的微分一积分方程—— 展的蒙特卡罗直接模拟(DSMC)方法b’4。，避开对方程的求解，直接计算微观粒子的运动，通过统计平均得到宏 观物理量．该方法简单易行，只是由于计算量巨大，限制了其更广泛的应用．随着计算机技术的飞速发展，DSMC 方法不仅在大尺度稀薄气体流动的计算中得到成功，近年来在微流动中也得到广泛应用，实践证明DSMC方法 是研究过渡领域流动的有力工具．本文采用蒙特卡洛直接模拟方法，对微尺度后台阶流动进行研究． 后台阶流动边界几何形状(图1)简单，但流动结构复 杂．在实际工程中有着广泛应用，也得到了广大学者的系 ^f 列研究．Annaly等¨1对经典后台阶流动进行了系统研究， Ⅳ 200时的流动结 实验测量了不同Reynolds数110Re7 10。 构．结果表明在实验船下，台阶后存在回流区，其大小由 工l 再附长度石．表示，并据此可以将流动分为层流、过渡流 图1 后台阶流动的几何外形和一般流动结构 The andthe 和湍流3个区域；在一定m范围内，上壁面存在一个跨 Fig．1 geometry generalbackward—facing structure 再附点对应位置的分离区；当1200Re1500时，在回 flow 流区后部还出现一个很小的分离区． 当尺度降到微米量级，胁变小，‰变大，流动将表现出不同的特征．wu等㈨在研究压力边界条件时给 出了一个类似后台阶流动的算例，但是没有做进一步的分析．xue等nj采用DSMC方法模拟了微尺度后台阶 ER=H／h对流动的定性影响．遗憾的是没有给出类似Re一戈。的关系，也没有得到流动结构发生显著变化的 临界参数．本文试图在这方面进行进一步的研究． [收稿日期】2004—10—22；[修回日期]2005—04—18 [基金项目]中国科学院创新(KJCX2．SW—L2)及中国科技大学青年基金资助项目 [作者简介】泰丰华(1975一)，