可压缩混合层流动的研究进展.docVIP

圆锥可压缩边界层流动国内最新研究进展1) 杨武兵 沈清 庄逢甘 中国航天空气动力技术研究院， 北京 100074 摘要： 本文将介绍近几年国内在圆锥可压缩边界层流动研究方面取得的部分最新进展，这些成果包括在分析方法、计算软件和实验技术上取得的突破，以及用于非对称转捩现象的研究进展，同时还包括转捩对气动力、气动热影响的一些研究。这些工作不仅在学术上有较高的价值，同时也是高速飞行器性能提升的需求，具有很强的工程应用背景。 引言： 圆锥是高速飞行器常采用的几何外形之一，对与之相关的外部流动，主要是边界层流动展开研究是此类飞行器性能提升的需求。边界层转捩是边界层研究中的重点，至今为止也是一个难点，高速使这一状况更为突出。首先，高速时气动热问题更为突出，对于马赫3以上的飞行器，热防护成为设计中的主要考虑，这实际上要求对转捩推迟技术给予更多的关注[1]；其次，高速时出现了一系列复杂气动现象，如激波/边界层干扰、激波/激波干扰等，这些现象对于气动力、气动热均有重要影响，对它们的正确预报一个很重要的前提是转捩预报是正确的[2]。 对转捩的控制首先需认识流动不稳定机制，因为确定流动的最主要不稳定机制将允许设计人员将精力集中在最有效的控制机制上[3]。对于边界层流动，线性不稳定性存在流向迁移不稳定性，交叉流不稳定性与G?rtler不稳定性。在流向迁移不稳定性中，Mack用无粘线性稳定性对平板的绝热边界层流动的计算表明在马赫数2.2以上同时存在具有不同波数和相速度的多种不稳定波，这些模态根据波数的增加编号，最低波数的模态为第一模态，通常称为涡模态，对应亚声速边界层中的T-S模态，更高模态通常称为Mack模态，常被描述为声模态。对于该不稳定性的控制存在较多研究，例如研究表明壁面冷却技术能稳定第一模态扰动，但是却对第二模态具有不稳定效应；又如被动多孔控制技术，该技术利用壁面对第二模态的频率范围的声学吸收特性，中断壁面和相对声学线之间的反馈机制，从而稳定第二模态扰动。交叉流不稳定性对于实际的三维飞行器而言非常重要，例如即便攻角很小，它也能支配流动稳定性并引起提前转捩。控制弹体形状以减小交叉流动是控制交叉流不稳定性的最实用的方法，但是飞行器几何形状的约束限制了可能造型的数目，而其它的控制选择也相当有限。例如研究已经表明抽吸控制技术在超声速时能推迟转捩，但是在高温流动中是不切实际的。G?rtler不稳定性因曲率而引起，对于高超声速引擎进气道和排气道处的转捩而言非常重要，但关于G?rtler不稳定性的控制仅有少量研究。 对于转捩预测，最常用的方法是eN方法，该方法在低速流动中得到了很好的验证，但是对于高速流动具体状况并不清楚。对高速时转捩位置及其长度的预测结果的改善同样依赖于对其物理机制的进一步了解。直接数值模拟与近来发展的PSE方法是研究转捩过程的有力工具；而在传统的计算方面，现有方法还不能考虑如三维平均流动、粗糙度以及大多bypass转捩等复杂效应；在实验方面，因传统风洞和激波风洞中的高噪声环境，使得大量地面实验数据并不明确，与此同时飞行实验也存在诸多限制[2]。 这种研究现状表明，转捩机制的研究，在高速时对相关方法提出了更加苛刻和迫切的需求，国内学者近几年对圆锥边界层问题展开了研究，在分析方法、计算软件、实验技术上取得了一系列成果，从不同的途径对某些问题如非对称转捩现象取得了新的认识，本文将对部分最新成果进行介绍。直接数值模拟与线性稳定性分析 近年，直接数值模拟在边界层和混合层等方面的工作日益增多，但计算结果的验证及其有效性是必须解决的一个问题。2006年，航天十一院的沈清等在之前解决了平板边界层线性稳定性问题的DNS与LST计算对比和相互验证问题的基础上，进一步在圆锥外形上解决了这一问题，如图1所示，该工作提升了DNS与LST在流动失稳机制上认识的能力，例如对非对称转捩机理的揭示。图2为圆锥边界层最不稳定波增长率，它揭示出，小攻角使迎风面趋于稳定，而背风面趋于不稳定，而且相对零攻角状态，背风面出现最不稳定波的起始位置前移，而迎风面出现最不稳定波的起始位置后移，从而形成非对称失稳。另外，背风面的最不稳定扰动波的增长率明显高于迎风面，这样不稳定扰动波在初始放大过程中将快于迎风面，表明了非对称转捩的出现趋势。据此可以推断，这是引起圆锥边界层非对称转捩的机理。 图1. 圆锥边界层T-S波特征函数DNS和LST结果的对比(左背风面，右迎风面) 图2. 圆锥边界层最不稳定波增长率（M=6） 在2006年，航天十一院还拓展了前期基于平板边界层发展的高精度DNS方法，开展了圆锥M=6.0边界层的转捩和湍流计算，初步给出了小攻角下出现的非对称转捩结果。图3、图4为计算结果，显示了圆锥边界层非对称转捩过程和湍流特性的初步计算结果，给出了非对称转捩的直