不同Kn数下圆柱稀薄气体绕流规律地研究.pdfVIP

…+————．篓±三曼童鐾生簦鎏整墨兰垒婆笙苎塑窒：!!篓兰!璺!!壁=!!璧 不同Kn数下圆柱稀薄气体 绕流规律的研究 杜永乐，阎超 (1l京魅空靛天大学CFD嚣象实验室，100083) 摘要零文采用Ds∽方法，埘过渡区内不同Kn数下圆柱超蒋速绕流进行了数值计算。计算鳊果表明，本文采用的 os抵方毯谣算鞲发鞴，计算缩暴清楚熬攒拳了气体秘稀薄效盔对超音速流动娆律静澎摘。 关键词：DSMC：隧牲绕滤；稀薄气体 0引言 一般气体动力学从气体介质是连续豹这个假设出发来研究气体的流动规律，谯气体分予的平均自 IH程与物体f内特征尺度相比不是小蹙时，连续性假设即失效，这种稀薄气体将呈现出与连续气体显装 夺霹鞠流魂特注。在靛天飞季亍及擞橇逛镶域，这种滚动现象经常窭臻，嚣要对耱藩气体豹流动特整体 出精确预测…。然而，对稀薄气体流幼的实验研究由于实验方法尚不成熟，设备复杂昂贵，而建立襁 适续气体模型上的实验方法秘求解N喝方程获褥数值解的方法也已经失效，因此必须寻找薪的计算方 法。 孙!tzmafl方程靛一致瞧”“”。实践袭骧，Bs黔方法能够缀容易的模熬三缑舔薄气体涟动麴复杂渡场， 而且能真实地模拟包含热化学非平衡反应等复杂的稀薄气体流动问题“”””’。经过近40年的研究， DsMc方法在稀薄气体动力举研究领域获得了广泛的应用，DsMc方法也趋于成熟，成为数使求解稀薄 4i律囊力学蠲嚣瞧一获褥蠢犬戒璃豹方法。 本文将采用DSMC方法，研究不阀Kn数下圆柱稀薄气体绕流的流动规律。 i DSMC方法描述 扶分f运动谂燕疫攒遴气体运秘斌律的方程是Boltzman穷程： 鲁+手·轰小嚣=＆晤牙(，’fl-，址盯b，z)ene舌l 个仿真分子代替葵实气体，通过随机抽样跟踪仿真分予的运动轨迹，模拉仿真分子的碰撞，疆焉统计 网格内模拟分子的运动状态来实现对真实气体流动的模拟。 羹二楚闲为D潮e方法接述气终分子运动物理过程鲤姆性，入粕缀骞易提爨簿食物理搜簿携各种气 体分子横碰、分子之间碰擒传能计算模型及分子和物筒碰撞反射的计算模型，从黼达到缀简的计算精 J窖。 574 一 至旦堑墼工旦壁篓翌墨堡丝堕塑堡塑堕塞 2计算模型 的影响。 模型。碰撞抽样采用Bird提出的NTC方法，气体分子碰撞传能计算中考虑到转动能松弛。但由于来 流马赫数并不太大，所以不计及振动能的激发。 图1所示为计算域和计算网格，设定的边界条件为：cE为进口边界，取自由来流条件；ED为出 口边界，取真空条什；AB为物面边界；CA和BD为对称边界；由于本文编制的程序为三维计算程序， 冈此在z方向上生成～个网格，其两边界面设为对称边界。cE和ED边界距离物面的距离为5R(R为 圆柱的半径)。 ／7一r1’、、、、 、 ／／ ．b。 ’： 、、、、、11 图1计算域和计算网格 3计算结果 表1为三种来流状态下圆柱阻力系数的实验值与计算值的对比，实验值取自参考文献[9]。需要 指出的是，状态3和4的计算中均未考虑转动松弛。根据参考文献[4]，气体分子的转动特征温度很 低，而振动特征温度很高，因此，在本文计算条件下，不计及振动能激发不会引起较大的误差，但忽 略转动松弛将会带来很大的误差，表1的计算结果清楚的说明了这一点。表1的计算结果表明，本文 的计算结果真实可信。 图2为三种来流状态下的密度等值线图和压力等值线图。可以看出，在Kn=0．016的近连续过 渡流区，流动特征类似连续流，圆柱前方出现清脆的弓形脱体激波；当Kn=0．08时，激波依然明显， 但由于单位体积内气体分子数减少，气体分子由于要运动较长距离才能与别的气体分子发生碰撞，冈 此流场在空间的变化比较缓慢，激波厚度明显增加；当Kn=0．4时，激波已经变得模糊不消。 圈3为三种来流状态下的流线图。在Kn=0．016的近连续过渡流区，流动特征类似连续流，自 由来流还未碰到物面就在跨越脱体激波时改变流动方向；而在Kn=0．4时，单位体积内气体分子数 明显减少，气体分子间碰撞频率降低，分子的粘性输运能力减弱，流动呈现出明显的附体特征。