Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.     

前后缘型线同时可控的乘波体设计

提出了一种前后缘型线同时可控的乘波体设计方法,在马赫数可控的外锥形曲面基准流场中,结合流线追踪技术和混合函数,实现了椭圆前缘转椭圆后缘的乘波体设计,并在设计点(Ma=6.0)和接力点(Ma=4.0)对其进行数值仿真研究。前后缘同时可控的乘波体在型面剧烈过渡处产生了较弱的激波,出口两侧存在高温高压区,后部对称面附近的激波形状由圆弧变为平直线且出口处流场基本均匀,非常有利于与进气道匹配设计。另外,该乘波体具有较高的容积率和预压缩效率,附面层修正后的容积率为0.24,设计点时乘波特性较好,接力点时前部完全乘波,具有较高的升阻比,有黏条件下设计点和接力点的升阻比分别为2.54和2.41。此外,与给定前缘的乘波体相比,其升力、阻力、俯仰力矩和出口增压比都有明显增加,但是升阻比和出口总压恢复系数有所降低,在设计点无黏升阻比由3.56降为3.00。以上研究表明,本文的设计方法可行且更加灵活,拓宽了乘波体的选择范围。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

Effect of blade tip winglet on the performance of a highly loaded transonic compressor rotor

The tip leakage flow has an important influence on the performance of transonic com-pressor. Blade tip winglet has been proved to be an effective method to control the tip leakage flow in compressor, while the physical mechanisms of blade tip winglet have been poorly understood. A numerical study for a highly loaded transonic compressor rotor has been conducted to understand the effect of varying the location of blade tip winglet on the performance of the rotor. Two kinds of tip winglet were designed and investigated. The effects of blade tip winglet on the compressor over-all performance, stability and tip flow structure were presented and discussed. It is found that the interaction of the tip winglet with the flow in the tip region is different when the winglet is located at suction-side or pressure-side of the blade tip. Results indicate that the suction-side winglet (SW) is ineffective to improve the performance of compressor rotor. In addition, a significant stall range extension equivalent to 33.74% with a very small penalty in efficiency can be obtained by the pressure-side winglet (PW). An attempt has been made to explain the fundamental mechanisms of blade tip winglet in detail.     

软式空中加油受油机头波数值仿真分析

针对软式空中加油对接过程中受油机头波加剧软管锥套飘摆、降低对接成功率的问题,建立了受油机三维模型,利用Delaunay方法自下向上生成非结构网格.通过计算ONERA M6三维翼型表面压力分布和压力系数,验证了三维Navier-Stokes方程和SST k-ω湍流模型的有效性.在此基础上对不同飞行条件下受油机头波进行静态数值计算,分析了头波随马赫数和迎角的变化规律,并通过仿真解释了头波影响下软管锥套的飘摆原因,预测了对接时锥套的运动趋势.     

烧蚀后的C-SiC隔离段内流场精细测量试验

为了获得极度粗糙内壁面对激波串流动特性的影响规律,通过基于纳米粒子示踪的平面激光散射技术和高频动态压力测量技术测量了Ma2来流下烧蚀后的C-SiC隔离段中激波串流场结构,获得了激波串初始激波形态、激波后附面层发展形态以及激波串动态特性。结果表明,烧蚀后的C-SiC隔离段中激波串结构与光滑不锈钢隔离段相似。但是极度粗糙的内壁面深刻影响了近壁区流动,附面层增厚效应非常明显。前者激波串内的附面层比后者厚约50%,前者激波分叉点比后者更接近唇口约30%。极度粗糙的内壁面也提高了附面层的分形维数,加剧了拟序结构的破碎程度。烧蚀后的C-SiC隔离段中附面层的分形维数在1.548～1.649,比光滑不锈钢隔离段高6.7%～8.9%。烧蚀后的C-SiC隔离段极度粗糙内壁面对激波串振荡频率几乎没有影响。激波串前传过程中的振荡频率约20Hz。     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

脉冲爆震燃烧室出口燃气能量分布特性

为了提高脉冲爆震燃烧室出口燃气能量转换效率,采用数值计算方法,研究了带气动阀和爆震增强结构的脉冲爆震燃烧室出口燃气能量分布特性,结果表明:脉冲爆震燃气的膨胀过程先后经历一次膨胀、二次膨胀和过度膨胀3个阶段,单个周期内二次膨胀阶段的时间和能量占比最大;燃烧室出口燃气压力势能、动能、内能和能流密度的分布主要受燃气压力的影响,且变化规律与燃气压力一致;一次膨胀阶段燃气增加的压力势能、动能、内能分别占单个周期的479%、259%、251%;二次膨胀阶段燃气增加的压力势能、动能、内能分别占单个周期的502%、576%、581%。