高马赫数来流条件下斜激波与流向涡对相互作用

针对斜激波入射流向涡对形成的复杂流动干扰现象,在来流马赫数6条件下,利用后掠斜坡式涡流发生器生成流向涡对,采用数值模拟结合激波风洞实验,研究了斜激波与流向涡对相互作用引起的激波变形和旋涡演变。结果表明,斜激波与流向涡对相互作用的剧烈程度大致分为强、中、弱三种。在强相互作用下,激波面向上游凸起的程度最剧烈,呈现类“T”字型特征,其对称面处接近于一道正激波,波后出现回流区和局部流动滞止,使得流向涡对发生破碎;在中等相互作用下,对称面处存在局部正激波,但没有出现回流区;在弱相互作用下,不再出现局部正激波。进一步地分析发现,斜激波与流向涡对相互作用产生的流向涡对破碎现象,不符合经典的激波作用下单个流向涡破碎理论的预测。采用流向涡对对称面处的气流参数,对经典的流向涡破碎理论进行修正后,能够预测本文流向涡对在斜激波作用下是否发生破碎。本文为掌握激波作用下流向涡对破碎现象的规律,提供了重要参考。     

宽速域翼尖涡及其与斜激波相互作用

针对宽速域飞行器产生的翼尖涡及其与斜激波相互作用问题,采用数值模拟方法,在来流马赫数0.2～6.0范围内,研究了翼尖涡的演化特性,并探讨了斜激波诱导翼尖涡破碎的机理。结果表明,在宽速域范围,翼尖涡的切向速度和环量分布仍然符合低速流动中的自相似关系式;在超声速至高超声速流动中,翼尖涡强度沿流向更快地减弱;当翼尖涡与斜激波相互作用时,涡核处压力越低的翼尖涡更容易破碎,而经典的涡破碎理论难以在宽速域范围内准确预测涡破碎现象;通过引入涡核处的压力亏损效应修正经典理论,大幅提升了预测斜激波诱导翼尖涡破碎的能力。     

圆柱跨声速绕流中的激波/湍流相互作用大涡模拟研究

采用大涡模拟方法研究了圆柱跨声速绕流中的激波/湍流相互作用问题,来流马赫数M∞取为0.75,基于圆柱直径D的雷诺数为2×105。计算结果表明,圆柱分离点处出现一道斜激波,并且以与涡脱泻Strouhal数一致的特征频率向上游传播。激波运动导致流场中出现反对称的流动模态,剪切层中压力信号的功率谱曲线中存在0.4、-1和-5次方的斜率关系,剪切层中的剪切应力角约为0°,脉动速度以流向脉动速度为主,并且沿剪切层的大尺度结构组织性减小。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

障碍物对激波触发爆震波影响的数值研究

为了研究不同几何形状障碍物对激波触发爆震波的影响,对氢气/空气混合物在不同马赫数的激波作用下的爆震波触发过程进行了数值模拟与分析。研究发现,对于不同几何结构的障碍物,触发爆震波的机理不同,热点的产生可以是激波与激波的相互作用,也可以是激波与壁面的相互作用。对于堵塞比较小的障碍物,合理的结构设计可以使成功触发爆震波的激波马赫数进一步减小。     

高超声速双模块内转式进气道的流动特性研究-Part I： 设计状态

本文通过特征线法设计了基于多道激波+等熵压缩波的基准流场,在此基础上通过流线追踪法设计了双模块内转式进气道。通过数值仿真和风洞试验相结合的方法,获得了内转式进气道的内外流特性。研究结果表明：在内转进气道最大半径对应的角区位置存在大量的边界层堆积,受第二道激波/边界层干扰,在激波根部卷起锥形旋涡;在内转式进气道内部,唇罩激波和管道边界层干扰显著,管道内存在自唇罩指向压缩面的强周向压力梯度,从而诱导管道内边界层均往一处汇聚,卷起大尺度流向涡。仿真和试验结果表明在来流马赫数5.74,攻角0度状态下,进气道气动性能优良,出口总压恢复系数系数达到0.58,最大抗反压为112倍。     

基于膨胀波效应的高超声速进气道肩部流动分离控制研究

为改善高超声速进气道唇口激波/附面层干扰诱导的肩部流动分离，从膨胀波及激波理论出发，推导出了膨胀波效应影响下的斜激波附面层干扰理论公式，获得了影响斜激波诱导分离的主要因素：膨胀角梯度、激波角及波前马赫数。在此基础上，开展了膨胀波效应影响下的流动分离控制研究，给出了膨胀波效应影响下斜激波诱导分离的判别及预测方法。结果表明：增大激波入射点处膨胀角梯度，可以显著减小甚至消除肩部流动分离；随着激波角增大，激波强度及逆压力梯度增加，分离区尺寸显著增大。而波前马赫数对分离区尺寸的影响不显著；在进口马赫数3.57～5.18，唇罩角度6°～10°范围内，当激波入射点处逆压比梯度小于250 m-1时，斜激波诱导的流动分离消失，可为改善超声速/高超声速进气道内流道流动分离提供技术支撑。     

马赫数2.5超燃室中超/亚声流流向涡掺混冷态流场数值研究

结合实验,对一种带亚声速预燃室和流向涡掺混器的超声速燃烧模型燃烧室,在其进口马赫数为2.5的来流条件下,进行了冷态流场的数值研究.计算与实验得到的燃烧室沿程压力分布相当一致,计算与实验得到的激波结构也基本吻合.计算结果表明:首先,在马赫数2.5的来流条件下,亚燃预燃室易于达到启动状态;其次,流向涡掺混器增强超/亚声速流之间的掺混的效果明显,但其掺混深度尚有限;最后,流向涡掺混的超燃室掺混段有着复杂的激波膨胀波波系,波涡干涉和激波附面层干涉结构.     

高超声速内转式进气道型面流场重构

内转式进气道流场参数分布不均，为改善该类进气道的气动性能，本文采用数值仿真方法开展了基于型面的内转式进气道流场重构研究。结果表明：流场重构型面中心线对进气道流场结构及流动特性影响较大，在给定偏距/长度与中心线末端斜率的约束条件下，选取合适的中心线起始角能够大幅提高进气道的气动性能，改善流场参数分布。与进气道原型方案相比，流场重构型面中心线10°起始角的进气道方案总压恢复系数、抗反压能力最大分别提升33.7%、26.4%，自起动马赫数下降1.1。随着流场重构型面中心线起始角增大，唇罩激波/侧壁边界层干扰诱发的流向涡减弱、流向涡传输轨迹向唇罩一侧偏移，低能流向唇罩两侧迁移趋势增强。在研究范围内，随着流场重构型面中心线起始角增大，隔离段出口总压恢复系数先增大后减小，自起动马赫数先下降后不变。     

斜激波与内凹半圆柱面无粘相互作用

针对三维内转式进气道流动中激波入射内凹壁面主导的复杂波系干扰问题，采用斜激波入射内凹半圆柱面的简化构型，在来流马赫数为6的条件下，通过无粘数值模拟结合理论分析，研究了激波角βi=14°～29°的斜激波在内凹半圆柱面反射形成的三维流场。结果表明，流场对称面均会出现显著高于二维情况的逆压梯度。当βi≤25°时，从侧壁到对称面，斜激波经历了从马赫反射（MR）到规则反射（RR）的转变，形成了MR-RR型流场，转变点处产生的桥激波向对称面延伸，桥激波在对称面反射后产生的压力峰值高于二维斜激波反射；当βi≥25°时，斜激波在侧壁和对称面均发生马赫反射，形成了MR-MR型流场，两种马赫反射分界点处产生的桥激波向侧壁发展，侧壁气流在对称面相撞后产生的压力峰值高于正激波后的压力；当βi=25°时，流场存在MR-RR型和MR-MR型双解现象。通过降维分析理论，揭示了两类流场中转变点和分界点的形成机制，并厘清了桥激波的产生原因和初期演化特征。当βi≥18°时，无粘激波干扰所主导的侧壁气流加剧向对称面汇聚，并在对称面附近产生流向涡对。无粘分析获得的认识，有助于揭示内转式进气道中流动汇聚和流向涡对等现象的形成机...     

斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。     

基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。     

耦合传热下激波对超声速气膜冷却影响

针对离散孔式超声速平板气膜冷却，在主流区引入楔角形成激波环境，以研究激波与超声速气膜之间的相互作用。通过计算楔角在0°、15°、20°和25°产生的四种激波强度下，超声速气膜与高温壁面的耦合传热。所得结果表明：适当强度的激波能够抑制气膜入射后产生的反向涡旋对，降低主流对气膜的卷吸，增大壁面平均H2摩尔分数并降低壁面温度。对金属层温度场的分析表明，壁面冷却效果随着激波角的增加而先增加后降低，其中楔角为20°时的流场结构最有利于壁面温度保护。小楔角生成的激波在低冷流马赫数下对冷却效果的改善更明显，大楔角则在高冷流马赫数下更明显，热障涂层（TBC）不影响这种变化趋势；激波的存在削弱了TBC的影响范围。可以揭示超声速气膜在耦合传热条件下的传热机理，为超声速气膜冷却的设计提供参考，或为现有超声速气膜冷却结构的优化提供依据。     

二维喷管的初始流动

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度混合WENO/TCD格式,对Ma=1.4的超声速平面射流初始流场进行了数值研究.数值结果清晰地描述了超声速平面射流初始流场的结构特征,包括主涡环与激波结构以及它们演变过程.因主涡环内存在涡导激波对,激波与涡相互作用加速射流剪切层失稳,使剪切层首次卷起形成小涡的位置出现在涡导激波对后,此与亚声速射流情况不同.小涡串卷起成后,相继与涡导激波对相互作用,使激波出现明显的变形并加速主涡环失稳.      

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.     

突跃型与平滑型斜爆震波起爆机制数值模拟

采用带有化学反应的Euler方程,对突跃型与平滑型2种形态的斜爆震波(ODW)在起爆机制和结构特点等方面的差异进行了研究.结果表明:斜激波(OSW)后是否存在亚声速区是判定斜爆震形态的依据.当波后火焰抬升斜激波面使得亚声速区存在时,横向激波与三波点结构就会出现并形成突跃型斜爆震.当波后所有区域皆为超声速区时,波后火焰燃烧无法影响到上游激波面,则会形成平滑型斜爆震.通过斜激波关系式给出了横向激波形成所需的临界斜激波角度,只有当斜激波角度大于此临界角度时才能形成突跃型斜爆震波,反之则形成平滑型斜爆震波.      

绕跨声速三角翼的激波/涡干扰流场数值模拟

 在绕三角翼的跨声速流动中,随着迎角的增加,三角翼上的涡破裂位置会出现突然前移的现象。针对这一与亚声速下不同的流动现象,采用带曲率修正的Spalart-Allmaras(SAR)湍流模型,求解定常雷诺平均Navier-Stokes(RANS)方程,对不同迎角下绕65°后掠尖前缘三角翼的跨声速流动进行数值模拟,并在此基础上,采用基于SAR湍流模型的脱体涡模拟(DES)方法,对由激波干扰导致的前缘涡破裂位置的运动规律进行了初步探讨。模拟结果与试验结果对比表明:SAR湍流模型能准确地模拟出三角翼上的激波系统和旋涡结构,并能准确模拟出由于激波干扰导致的涡破裂位置突然前移的现象。此外,对涡破裂后流场的非定常数值研究发现,支架前端正激波的干扰作用使得涡破裂位置向下游移动比较突然,而向上游移动则相对缓慢。     

Shock/shock interactions between bodies and wings

This paper examines the Shock/Shock Interactions (SSI) between the body and wing of aircraft in supersonic flows. The body is simplified to a flat wedge and the wing is assumed to be a sharp wing. The theoretical spatial dimension reduction method, which transforms the 3D problem into a 2D one, is used to analyze the SSI between the body and wing. The temperature and pressure behind the Mach stem induced by the wing and body are obtained, and the wave configurations in the corner are determined. Numerical validations are conducted by solving the inviscid Euler equations in 3D with a Non-oscillatory and Non-free-parameters Dissipative (NND) finite difference scheme. Good agreements between the theoretical and numerical results are obtained. Additionally, the effects of the wedge angle and sweep angle on wave configurations and flow field are considered numerically and theoretically. The influences of wedge angle are significant, whereas the effects of sweep angle on wave configurations are negligible. This paper provides useful information for the design and thermal protection of aircraft in supersonic and hypersonic flows.     

斜置机翼技术的研究与发展

介绍了斜置机翼技术的发展过程,并对国外相关研究进行了分析和整理;根据超声速机翼经典理论,说明了斜置机翼可有效减小高速飞行阻力的原因及其可变几何外形的能力;分析了斜置机翼独特的非对称布局带来的非线性高度耦合的气动特性和飞行动力学特性,并阐述了针对这些特性所需的各个方面的设计;最后对斜置机翼的应用前景进行了展望。