非定常平板激波%2f湍流边界层干扰.pdfVIP

第十四届全国计算流体力学会议论文集 非定常平板激波／湍流边界层干扰 潘宏禄，马汉东，沈清 (中国航天空气动力技术研究院，北京100074) 摘 要：以高超声速发动机进气道湍流分离控制为应用背景，采用大涡模拟方法进行Mach数为3．o (唇口附近Mach约为3)激波／湍流边界层干扰流场机理研究。当前含激波的非定常湍流模拟是一个难点， 本文利用扰动循环引入的方法，先得到充分发展湍流场，后根据斜激波关系式引入激波的方法进行激波／ 湍流干扰模拟，结果显示：充分发展湍流场在激波作用下产生逆压梯度并发生分离，摩阻分布与实验一致， 湍流可以有效减小分离区长度和强度，这为进气道分离区湍流控制研究提供理论依据。反射激波形成过程 中其对下部湍流边界层产生严重抑止和压缩作用，随着反射激波抬升抑止作用减弱，同时边界层对激波形 成反作用，边界层附近激波发生扭曲。 关键词：湍流：转捩；超声速边界层；激波／边界层干扰． 0引言 超燃发动机研究囊括了复杂的气动力和气动热问题，其前端进气道包含了多种激波干扰(侧板 压缩使其更为复杂)，激波／湍流边界层干扰(SWTBLI)便是其中一类重要的气动热、力问题。 激波／湍流边界层干扰是超／高超声速流场中典型的流动问题之一，如自由干扰模型、Hump模型、 压缩拐角、前(后)台阶流动等(图1)，包含了复杂的非定常 动力学特征(如强压缩性、多尺度干扰等)，给高速飞行器气 动力、气动热研究带来了重大挑战，是目前最为困难的模拟 和预测问题之一。 SwTBLI问题研究始于上世纪中期，到目前为止虽然发展 了众多新的实验技术和数值方法，但SWTBLI问题仍然没有得 到很好的解决，正如Georgeu。指出的那样经典的边界层湍流 理论(如壁面律、速度拐点理论、过渡层对数律等)已经能 够很好的融合到当前的数值(如RANS、大涡模拟(LargeEddy 图1slrrBLl问题示意图 Simulation，LES)等)研究中，所建立的数值方法也已经能 够对弱压力梯度流场给出很好的模拟，但对于涉及强压力梯度及更复杂的流动情况，目前的研究手 直接模拟蒙特卡洛方法(DSMC)的黏性边界层求解已经较为成熟，目前需要侧重于复杂流动(如存 在粗糙元、吹／吸气、激波干扰)下的流场物理现象和作用特征(如流场内及壁面附近的化学反应、 边界层或自由剪切层转捩、湍流等)的研究。 Numerical Grasso【5】针对SWTBLI中的自由干扰问题，采用局部自干扰方法进行研究，分析了流场结构、噪声等 对不稳定模态形成机理的影响，同年，Zheltovodov[6】给出了压缩拐角直接数值模拟与实验的对比， 两者分离区大小及摩阻分布相当吻合，充分展示了高级数值模拟方法在SWTBLI问题研究上的有效 性。相似工作还有Wu8d4artinu1、Morgan[8】和Postl【9】等。 潘宏禄，马汉东等：非定常平板激波／湍流边界层干扰 国内在SWTBLI研究中，主要还是以RANS为手段进行预测和分析。 本文以平面激波湍流边界层干扰为抽象对象，采用基于结构网格的LES方法进行三维非定常数 值模拟分析，目的是结合实验来开展高超进气道湍流控制研究。 1 计算模型及数值方法 闭处理，方程如下： 望+堕：o 8t ％ 亟+望堕：一至+旦 al a)cj axl ax} 攀+毒(厕+刃乃=丢(ej+Fe舀i)Ot 出。 。 观 其中： 万=础于， rF=岛+巧，oJ=吼+彩 岛=衍)(善+鲁一；岛挈，乃叫扔善 勺=一声(珥～一玩乃)，彰=一砩(“，