有变高度圆柱干扰的平板表面流谱拓扑结构-计算物理.PDF

[文章簟号]1001—246x(2000)04—0355．05 有变高度圆柱干扰的平板表面 流谱拓扑结构 马汉东， 周伟江 (北京空气莉力研究所．北京100074) [摘要]采用分区方法和R0e三阶流通量差分分裂格式数值求解雷诺平均N—s方程，湍流附加 扑学理论给出了有变高度圆柱干扰的平板表面流谱拓扑结构，并结台对称截面流谱进行了简要的 分析与探讨，指出了需要进一步深^研究的问题。 f关t词]圆柱；湍流；数值模拟；拓扑分析 f中圈分类号】0354 [文献标识码]A O 引言 在高超声速流动中，平板表面局部凸起的圆柱将使平板边界层流动发生根本性的变化。 同时，随着圆柱高度效应的引人，流场更趋复杂化，对应的物面流谱往往也会有所不同¨J。 的流场控制参数，相应地可将干扰流场划分为三种不同的类型，即非半无限干扰或称有尺度干 扰、过渡态和半无限干扰三类，这说明变高度圆柱将会导致不同的流场参数分布，而且干扰流 场的变化必然要反映到物面流谱上来。受实验条件以及流场复杂性的限制，实验油流谱很难 捕捉到细微的流动结构，如二次分离线。本文试图通过数值模拟来探讨物面流谱的变化。 由于数值模拟的流场分辨率往往受计算方法和网格布局的限制，这类问题的数值研究实 阶有限差分方法，数值模拟了变高度圆柱诱导的激渡边界层干扰流动，得到了与实验结果基本 一致的物面气动载荷分布和流场波系结构，本文将根据拓扑学理论对平板表面流谱进行分析 研究，进一步探讨需要研究的问题。 1数值模型及计算方法 ． 面传热的影响。基本计算条件【2]为：来流马赫数MⅡ=5，总温To=400 5．3×107m，未扰动模型当地湍流边界层厚度6=2lmm，壁温Tw=288．16 25mm．柱高分别为^／D=0．4、1．0、2．0、2．4和一 0922 [收稿日期]1999．07．22；【謦订日期】1999 [作者简开]马汉求(1961～)．男，安徽宿县，研究员，博士，从事计算流体动力学方面的研究 }} 葬 鐾 淫 筹17卷 港诗及举霹鏊霜露静糖互影噙。分区透嚣采霜拼萋鑫弼掩蓑拳，匿避对瘟耪壤量可矗接浃菇，最 慧数约为弼万，壶对阚裙美法遥截计算，其体细节译觅文：3]。 2平板表衙流谱拓棒结构分祈 其它主赘特征上看，它们与裳验袖流谱“a究全一致。下面将分析其拓扑结构。 根攥拓扑努析理沧，一个曲爵向量场援满足妁拓扑规律为}{1 《∑N+专∑v)(∑s+专∑s’)=z “ t ， 、 一 ， 7 为半鞍点数总襁，?为蛙面蛉黢拉拳性数，对本铡情况越霞可考虑为一环辅，匿此y=一】． 黼b归一e4 封^弗一l、0 ∞^∞=2，0 (d}^妒一。 隔l牛板表面极限流犊 1 sur如ce Fjg Limtting剥r钮mlin髓onplate 关子分离戥有如下结论：5]：在物津表两，着分离续藏奇点始，菝奇点盛为鞍点；若分离线 雀奇最终，该奇点必为结点：蒋分离缝上有很多奇点，其鞍点和焙点是交替分布的。 0．4)瘫谱奇点分布如熙2(a)，奇点计数如下 ∑s=4，∑s’=7，∑N=6，∑N+=l 第4期