欠膨胀高速湍流射近场特征的数值模拟.pdfVIP

欠膨胀高速湍流射流近场特性的数值模拟 戚隆溪 马言宝 江先金 王柏懿 (中国科学院力学研究所，北京100080) 正摘要】本袁发展了分析欠膨胀高速湍流射流近场特性的敷值方法，计算了不同背压比务件下声速 喷管出口附近膨^k区的波系结构． 美t诩：欠膨胀射流，湍流，声速啧管，近场特性，胞格结构 一引言 近年来，高超声速推进器和超声速割焊喷枪的研制以及油气井喷发等灾害事故的控制等 众多应用引起了人们的重视，它们均与欠膨胀高速湍流射流有关，因此要求对其流场结构作 深入的研究。在各类射流问题当中，欠膨胀高速湍流射流由于激波的存在而导致了流场结构 相当复杂。本文采用数值方法分析这类射流的近场特性并给出射流喷嘴附近膨胀区内的波 系结构。 对于无粘射流，若干作者[卜4]基于Euler方程并分别应用有限体积方法和有限差分方 法(涉及到NND格式和自适应网格等)对欠膨胀高速射流进行了数值模拟，较好地捕获了激 波，给出了胞格结构。特别是第一个激波结构符合得很好。对欠膨胀射流来说．只要背压比 不接近l，第一个激波的位置很稳定，且与出口距离较近，湍流混合影响不大，所以采用Euler 方程计算时，第一个激波位置与实验值相比符合很好。但是，随着射流发展，湍流混合效应 增强，仍用Euler方程计算过渡区和远场区流场就会有较大误差。 对于粘性射流，Dash[5]等人基于抛物化NS方程(PNS方程)，采用空间推进方法发展了 一整套计算各种与射流性态有关的SCIPvIS程序，可计算定常欠膨胀高速射流流场。计算中 根据流动特性把流场划分为超声速区和亚声速区，超声速区叉可根据粘性影响大a吩为有粘 区和无粘区。这种方法把分析有粘和无粘的超声速区的激波捕捉法与分析亚声速区的压力分 裂法结合起来．并在声速线附近用特征耦合关系进行匹配。当马赫数在2．O附近且背压比不 大时，马赫盘较小，计算结果与实验结果[6]相比符合得较好，但当马赫敷远离2．O或在大 背压比情况下．计算结果与实验相比有较大误差。 二模型方程组 本文的理论模型为：(1)采用Ns方程描述流体运动；(2)计及射流内部压力场的非 均匀性；(3)采用零方程湍流模型。我们采用时间相关法计算定常解，这样控制方程为柱坐 标(x，r)中的F“re平均Navier-stokes方程： 型+丝+竺：堡+丝+墨二!!坚 … at ax西 dx ar r 、1 方程(1)中各个参量为： p“ l p” I 0 O U= E= p竺pI F： !”一f；矿： p“” l p“V+pl —f甜+p 引引 p“P+p“j p～e+p、) 0 0 I o I E。= o I； 只5 +c鲁+静罢l ～u ～v 。。。唯 + 旦丹 盟瞳 其中p，p，r和e为气体的压力、密厦、温魇平口比总能：“和V为远展分量(上划线一丰Ⅱ一分 量可表示为 铲扣巾罢一弘衅+争j 盘j r 口 rF=(∥+“)(罢+宰) ”扣“，雾一≯“x手+争 铲扣的写c…x罢+争 其中下标t表示湍流特性。气体的状态方程取下述形式 芦=(，一1)芦[虿一；(i2+彳2)】