可压缩流向涡与激波轴对称干扰的数值模拟力学学报.PDF

第 31 卷 第 3 期 力 　学 　学 　报 Vol. 31 , No. 3 1999 年 5 月 ACTA MECHANICA SINICA May ,1999 可压缩流向涡与激波轴对称干扰的数值模拟1) 刘延宁 　尹协振 　韩肇元 ( 中国科学技术大学力学与机械工程系,合肥 230026) 摘要 　用 NS 方程数值模拟了可压缩流向涡和激波轴对称相互作用现象. 数值模拟包括定常 和非定常两种情况 ,计算结果分别与相应的实验进行了比较. 结果表明数值模拟成功地捕捉到 了激波和旋涡相互作用过程中发生的激波波面变形 ,激波振荡 ,涡核变大以及激波波后出现驻 点、回流区等流场特征. 提出了判断流向涡与运动激波相互作用中旋涡破碎的准则. 关键词 　可压缩流向涡 ,激波 ,数值模拟 ,激波和旋涡相互作用 引　言 激波与可压缩流向涡相互作用现象具有很强的工程实用背景 ,避免和控制激波与旋涡相 互作用后旋涡的破碎有重要的实际应用价值. 60 年代以来 ,飞行器设计中广泛采用了边条翼 和近距耦合鸭翼的布局方案. 当飞行器在跨、超声速区飞行时 ,边条翼和鸭翼产生的流向涡系 会与后翼面上产生的激波发生相互作用 ,从而产生非定常载荷 ,影响飞行器的机动性能. 同样 , 当这些旋涡与进气道产生的激波波系相遇时 , 由于激波和旋涡相互作用发生的旋涡破碎可能 会降低进气效率. 在这些情形中 ,旋涡的破碎是应当尽量避免的. 但是 ,在最近兴起的超燃发动 机的研究中 ,有人提出 ,可以利用旋涡与激波相互作用后旋涡的破碎现象来提高燃料与空气的 混合效率. 基于此 ,近年来 ,各国学者围绕可压缩流向涡与激波的相互作用 ,在实验与计算两方 面进行了大量的研究工作. 在实验研究中 ,有两种产生流向涡的方法 :一种由旋涡发生器产生[1 ,2 ] ,另一种利用机翼 产生的翼尖涡作为所研究的流向涡[3 ,4 ] . 这些实验都是在超声速风洞中进行的. 不同实验方案 下的结果都表明 ,激波与流向涡相互作用后 ,流场呈现出非定常性 ,激波波面发生程度不同的 变形 ;激波位置在平衡位置附近振荡 ;旋涡涡核尺度变大 ;一定条件下 ,激波波后出现驻点、回 流区 ,显示出旋涡破碎的流场特征. 激波与流向涡的相互作用可以分为强相互作用和弱相互作用两种. 一般认为 ,强相互作用 ( ) 时 ,激波波后出现驻点、回流泡 这也是旋涡破碎的标志 ,激波波面变形尺度与涡核尺度相比 大很多 ,有量级差别. 弱相互作用时 ,激波波面变形程度与涡核尺度在一个量级上 ,旋涡在通过 激波波面后 ,仍保持着紧致的涡核. 数值计算方面 , Copening 和 Anderson [5 ]用 Euler 方程模拟了流向涡与驻定斜激波的相互 作用 ,他们未捕捉到旋涡破碎现象. Kandil[6 ]用层流全 NS 方程分别模拟了准轴对称和三维两 种情况下的激波与旋涡相互作用 ,其中包括内流和外流两类问题. 其结果显示了激波与旋涡相 　1998 - 03 - 03 收到第一稿 , 1998 - 06 - 29 收到修改稿. 1) 国防预研基金资助项目. 　　　　　　　 力 　　　学 　　　学 　　　报 1999 年 第 31 卷 258 互作用后流场的非定常性 ,并捕捉到了旋涡的螺旋型和泡型破碎. Rizzeta[7 ] 则直接对翼尖涡 的产生过程以及旋涡与激波的相互作用进行了数值模拟. 以上报道的工作主要集中研究了驻定激波与流向涡相互作用的现象,未见有运动激波与 流向涡相互作用方面的报道. 我们曾经在激波管内开展了运动激波与流向涡的相互作用的实 验[4 ] ,并用轴对称 Euler 方程进行了相应的数值模拟[8 ] . 在本文中 ,我们用 NS 方程计算了两种 激波和流向涡轴对称相互作用的现象. 第一种情况 ,是定常激波与流向涡相互作用 ,如图 1 (a) 所示. 在一直通管的管口放置一旋涡发