滑移过渡区微喷管流动的Burnett方程研究.pdfVIP

滑移过渡区微喷管流动的Burnett 方程研究1) * *2) * + * 王晶鑫 ，包福兵 ，林建忠 ，余霞 * （中国计量学院流体检测与仿真研究所，杭州 310018 ） + （浙江大学流体传动与控制国家重点实验室，杭州 310027 ） 摘要：采用 Burnett 方程和滑移边界条件对滑移过渡区微喷管的二维流场进行了模拟研究。分析了不同努森数对 音速点位置的影响，发现随着喉部尺寸减小，微尺度效应的影响变大，音速点向着出口移动，已经不再固定在喉 部。分析了进口压强对推力和比冲的影响，随着入口压强的变大推力系数也变大，而比冲只有略微的增长。也分 析了不同收缩比和扩张比对推力和比冲的影响，并得到了高效率的微喷管结构。 关键字：滑移过渡区；Burnett 方程；边界滑移；推力；比冲 1 引言 近年来，微机电系统 （MEMS ）和纳机电系统 （NEMS ）迅猛发展，微纳机电系统中的微纳尺度流动 [1-3] 也日益受到国内外学者的关注 。微纳机电系统中气体流动和传热特性的预测，是其设计发展的重要部分。 在微纳尺度流动中，流场特征尺度小，接近于气体平均分子自由程，努森数，Kn  ，不再是个小数， L 流动达到滑移或者过渡流区。 微型航天器在航天和军事领域发挥着日益重要的作用，微喷管是微型航天器推进系统产生推力的主要 部件。尺度的微型化给喷管内部流场模拟、喷管外形设计、流场以及推进性能的测试带来新困难。近年来 [4] [5] 微喷管内部流动在国内外逐渐受到重视，许多学者对其进行了模拟研究 。孙建威等人 采用Navier-Stokes 方程对微型缩放喷管内的流场进行了二维数值模拟，在计算中引入滑移边界条件，通过比较有、无滑移边 界条件的数值结果，发现随着微喷管内Kn 数的增大，壁面上的边界滑移对流场和微喷管性能参数的影响 逐渐增大，当流场处于滑移过渡区时，滑移边界条件不仅对流场的影响已经较明显，对微喷管性能预测的 [6] 影响也不能忽略。张先锋等人 采用滑移边界条件的连续介质模型对喷管的二维流场进行数值模拟，分析 了不同进口温度及壁面换热条件对喷管性能的影响。其研究结果表明随着进口温度的增加，推力、推力效 率以及比冲效率减小，而比冲不断增大；在给定进口压力和温度下，绝热壁面的推力效率最小；壁面温度 越低，推力效率越大。但是对于滑移过渡区的流场，气体偏离了平衡态，Navier-Stokes 方程已经不再适用。 研究证明，采用 Chapman-Enskog 展开对 Boltzmann 方程取高阶近似所得到的Burnett 方程，能够描述气 体在轻微偏离平衡态时应力和应变之间非线性的本构关系，比Navier-Stokes 方程更适合于描述滑移过渡区 [7-8] 的稀薄气体运动 。 2 计算方法及模型 1) 浙江省重点科技创新团队 （2009R50024）和国家自然科学基金资助项目 2) E-mail ：dingobao@zju.edu.cn 2.1 Burnett 方程 [9] Burnett 方程是Boltzmann 方程的二阶近似，能够描述轻微偏离热力学平衡时的稀薄气体流动 。Zhong 为