超音速方腔流动的数值模拟研究.pdf

超音速方腔流动的数值模拟 王涛申义庆杨国伟 中国科学院力学研究所高温气体动力学实验室，北京，100080 摘要：超音速方腔流动普遍存在于航空航天飞行器中，本文首先利用K-G两方程湍流模型数值模拟了开腔状 态下的流场，计算结果与文献有关实验和计算的结果相吻合。进一步模拟了在方腔上不同位置增加盖板时的 流动状态，结果表明：有无盖板其流动有很大的不同，如腔底压力明显升高，后向台阶后的尾流涡内出现两 个反向涡等。 O引言 超音速方腔流动普遍存在于航空航天飞行器中，如现代隐身战斗机内置导弹的发射时仓 盖的打开／闭合过程，超音速发动机内的燃烧过程等。方腔流动存在边界层分离、旋涡、激波／ 膨胀波以及激波与边界层的相互干扰等现象，因此对超音速方腔流动的数值模拟不仅有重要 的学术意义，也具有较大的实际应用价值。 方腔流动由于其长深比的不同，具有开腔和闭腔两种流动状态，已有较多的文献报导了 有关的实验、计算和理论分析的结果【1‘副，其中Zhang等【41给出了不同长深比方腔底部的压力 实验数据以及比较了几种涡粘性模型，如代数湍流模型、七一倒模型等的数值计算结果。当长 深比三／D10时，剪切层跨越了整个方腔，将腔内流体与外部气流分开，方腔流动呈开腔 状态，参见图la；当长深比三／D13时，方腔流动一般为闭腔状态，分离尾流和压缩流能 较好分开，在它们之间存在再附边界层与主流流动，在两个激波之间形成一压缩扇形区域， 参见图lb。 图1a闭腔流动 图1b开腔流动 本文首先利用K。G两方程模型计算了长深比为15的闭腔流动状态，所得结果与Zhang 等‘41的计算与实验结果一致，进一步模拟了方腔打开腔盖在不同位置的流场，数值表明，当方 腔上方存在腔盖时，流场结构发生了较大的变化，如腔底压力升高，腔盖位于一定高度时， 后向台阶后的大涡内出现两个反向涡。 1控制方程和数值方法 本文采用完全NS方程和K-G两方程湍流模型”】 二+——+————=—上+—o—o+-4- ¨l 警十娑+挈掣：莩+孥孥+S4- (1) af。孤。却 az 舐 却 az 一 其中，Q=【P，pu，pv，jaw，pe，肚，昭】1， 第十一届全国激波与激波管学术会议 F=【肛，pu2+p，puv，puw，u(m+p)，puk，pug]r， G=【∥，puv，pv2+p，pvw，v(pe+p)，,ovk，印苫】r， H=【∥，puw，pvw,／升v2+p，w(pe+p)，pwk，p孵】r， 州o，o，0’o，o舡竿，口簪最一警 一(j,￡L+pT／Og)华(甜2矿2剃， 肚max(o．01旷∥w)，足=《善， a似,,。‘801， f孑=以c考+等，一；摹磊，一；肚毛，勺=以c等+警，一j2 2 ∥r=∥+触2，gJ P1)p。+r1V2 W2)+七 2一(彭三+rrJ-it)T，江(r 0I=09=2，D=5／9，芦=0．075，∥+=0．09。 对方程(1)的离散，本文采用LU．SGS隐式迭代法【61，粘性项及源项采用二阶中心差分离 散，对流项用Roe格式的改进形式HLLE格式【81来离散，由于增加的K、G两个方程可写为守 恒形式，因此在HLLE格式中特征值矩阵为七对角阵，女11 Jacobi矩阵彳=薏的特征值矩阵 求得各Jacobi矩阵的左、右特征向量矩阵。 2计算结果与分析 块网格划分，将计算区域划分为三块，在各个壁面附近进行了加密，图2为网格划分示意图。 304 图2网格划分示