用可压缩性修正两方程湍流模型数值模拟激波%2f湍流附面层干扰流动.pdfVIP

中图航空学会推进系统气动热力学々业委员会 CSAA99_—02i 江苏(苏鲁皖)航空航天学会动力装置专业委员会 中国人民解放军空军工程学院 ’99学术会议 用可压缩性修正两方程湍流模型 数值模拟激波／湍流附面层干扰流动 王国庆梁德旺 (南京航空航天大学动力工程系，南京210016) 摘要本文采用可压缩性修正的两方程湍流模型，数值模拟了跨声速2D喷管管内三 种不同波前马赫数的激波／湍流附面层干扰流动．不仅对流场中时均参数的计算结果与 实验值进行了比较，而且对咏动参数也作了详尽地比较．结果表明，可压缩性修正的 两方程湍流模型不仅可以准确地模拟正激渡／湍流附面层干扰流动的时均参数．而且可 以较准确地模拟脉动参数：舾端动能女，雷诺应力等．它不仅可以准确地模拟无分离的 激波，湍流附面层干扰流动，而且能峙正确地模拟引起流动分离的激波，湍流附面层干 扰流动各参数分布的基本规律． 关键词可压流，激澍湍流附面层干扰，湍流模型，数值模拟 1引言 在许多工程实际中，如飞机机翼、压气机叶栅通道以及进气道等跨声速飞行器部件中 都可能遇到诸如激波，湍流附面层干扰等问题。在激波，湍流附面层干扰区域，沿着流动方 向流动参数变化剧烈，激波前后存在很大压力梯度。在激波根部很强的粘性作用使附面层 厚度迅速增长，当激波足够强时，甚至会引起附面层分离，粘性的影响此时更加突出。囡 此采用数值方法模拟这些流动时，激澍湍流附面层干扰是造成计算精度不高的关键所在。 就求解激波，湍流附面层干扰流动而言，人们现在一般采用全场求解NS方程和 对近壁雷诺应力的影响，从而大丈提高了正激波，湍流附面层干扰计算的准确度。尽管如此， B，L湍流模型是基于平板湍流附面层实验的经验模型．它对外流有比较好的计算精度，而 对内流则有较大的近似性。针对内流一般采用两方程湍流模型。 关于用两方程湍流模型计算正激波，湍流附面层干扰问题，LaunderE5l认为几乎不可能， 而粱德旺教授也曾采用两方程湍流模型数值模拟Bohning喷管fll内正激波，湍流附面层干扰 流动但未能成功。仔细分析现有的可压缩流体两方程湍流模型，可以发现这些模型都是将 不可压模型简单地扩展成为变密度的可压流体之形式，而本质上都承认Morkov／n假醺， 即湍流马赫数埘可以忽略不计，这种近似处理对高速及超高速流动来说是不适宜的Kl。 sarka一1研究发现，随着湍流马赫数的增大，流体参数(温度T，密度p，压力P等)的脉 动量也逐渐增大，此时流体的可压缩性带来的影响会引起湍流结构及湍流动力性能的变 化。所以忽略湍流马赫数赫的湍流模型，将难以准确估计可压缩性带来的影响．从而不 适用于高速及超高速湍流流动的横拟。DavidC．Wilcoxt”，SarkarS．pⅡ”zeman0．f9I等在流体 一106一 王国庆粱德旺 用可压缩性修正两方程湍流模型数值模拟激澍湍流附唾层干扰流动 的可压缩性对模拟高速及超高速湍流流动的影响方面进行了探索性研究，荠分别提出了各 自的关于膨胀耗散或压力膨胀的可压缩模型。在学习和研究前人成果的基础上，文献[10) 给出了同时包括膨胀耗散和压力膨胀的可压缩性修正的两方程湍流模型，并采用可压缩性 修正厢的两方程湍流模型数值模拟Bonhing喷管内正激波，湍流附面层干扰流动，比较了修 正后湍流流场中的时均参数(壁面压力分布和波后速度分布)计算结果和实验值。 由于对脉动量的分析能够较全面地反映揣流模型对湍流流场的模拟，所以本文作者将 采用可压缩性修正两方程湍流模型【tm对喷管内正激波，湍流附面层干扰流动进行数值模拟， 数据。文中不仅对流场中时均参数的计算结果与实验值进行比较，而且对脉动参数也进行 比较，同时对脉动量在未分离和出现分离等状态下的分布情况进行了分析。 2基本方程和湍流模型 2．1基本方程 昙(刃去(再)-0 妄(碱)+丢喊可+确)=旦dxj(勺) 丢(函+毒(西nt刃5毒(可～1) Prandtl数(晶=o．72)，晶为气体湍流Prandtl数(堵=0．90)，￡：i+V2萌+k。