振动Couette流稀薄效应的直接模拟Monte Carlo方法.pdfVIP

第 32卷第4期 南 昌大学学报(理科版) Vo1．32NO．4 2008句：8月 JournalofNanchangUniversity(NaturalScience) Aug．2008 文章编号：1006一o464(2008)04—0353一O3 振动 Couette流稀薄效应的直接模拟 MonteCarlo方法 谷 黎，姜水生 (南昌大学 机 电工程学院，江西 南昌 330029) 摘 要：以不稳定直接模拟MonteCarlo方法 (DSMC)对 以时间周期振动的Couetle流作出特征参数方面的研究 汁 算以Knudseu(Kn)andStokes(fi)数为特征函数。另外，本文以适用于振动Couette流 自由分子流动区的线性尤碰撞 Boltzmann方程导出解析方法 ，并且利用它验证 了在 自由分子流动区的DSMC结果 关键词 ：振动 Couette流 ；DSMC；Bohzmann 中图分类号：TBI26 文献标识码 ：A 随着没备零件的小 化与在低于大气压下装配 1 DSMC方法 零件出现问题，人们玎始需要考虑稀薄效应的影Ⅱ向。 稀薄效应的特征参数定义为 Knudsen数，Kn= 我们假设稀薄气体流处于两_尤限大平行平板之 A／L，这里，A是气体分子平均 自ff{ ，L是流动特征 间，平板间距 L，顶板在横 向方向 卜做有规律振动， 长度。流动可以分为以下儿个区域：连续统一体流 底板静止不动，如图1。两板保持在相等温度 T = (Kn0．O1)，滑移区(0．O1Kn0．1)，过渡 273K。气体介质初始状态静止，分子数密度为 2／，介 (0．1Knl0)， 自由分子区(Kn10) 质温度 =T 振动稀薄Couette流的流动特征以 卣接模拟蒙特 卜洛方法(DSMC)足一种基于粒 特征参数 Knudsen(K2／)和 Stokes数 ()表示。 子的方法 ，并}_I．适川于所有流动 域。它应用一种·随 Stokes数定义为 机算法估算碰撞可能性 j散射分布。尽管有更简单 u m~singot) 一 — — 一 一 的理论模型，例如带有滑移修iE的Navier—Stokes ’L_-___lI_l-__。■ll_-_l· -·__I·l一 方程或无碰撞 Bohzmann方程 ，但它们只适用于滑 移区和 自由分子区 Bird 的DSMC方法使之成为 可能。从此方法引入以来，它被大量应用于稀薄气体 的流动模拟。由于MEMS系统存在大缱吖移动配件， 不稳定DSMC计算就成为必需 现有文献对稳定的稀薄状态 Couetle流有所研 图 1 振动 Couette流图 究，但对以时间周期的振动Couette流没有很深的动 力学研究理论 。在一【程应用方面，不稳定流由于数值 = = (LZ／v] 计算的复杂 Ej存储的高要求，一直是个难题：振动 Couette流在各种 MEMS系统叶1都有