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

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了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-ω湍流模型的有效性.在此基础上对不同飞行条件下受油机头波进行静态数值计算,分析了头波随马赫数和迎角的变化规律,并通过仿真解释了头波影响下软管锥套的飘摆原因,预测了对接时锥套的运动趋势.     

水下超声速射流对上浮水雷受力特性影响研究

固体火箭发动机具有功率密度大、推力大等优势,常被用于上浮水雷的推进器。然而发动机在水下工作时燃气射流使流场压力发生剧烈脉动,进而影响发动机的推力性能及上浮水雷的受力特性。基于VOF多相流模型和理想气体模型,建立了上浮水雷在不同工作状态(欠膨胀、完全膨胀、过膨胀)下的数值模型,研究了水下燃气喷射流对上浮水雷的受力特性影响。结果显示,欠膨胀工况时,发动机推力平缓,大小为12.2KN,上浮水雷受力未出现负值;完全膨胀时,射流发生颈缩、胀鼓现象位置距离雷体较远,导致发动机推力及水雷壳体受力振荡不剧烈;过膨胀工况时,射流发生颈缩、胀鼓现象距离喷管较近,发动机推力发生剧烈脉动,产生21.37%的振荡幅度,胀鼓现象发生时,流场压力显著降低使得上浮水雷后体受力减小,壳体阻力增大,上浮水雷最大产生27KN的负推力。     

烧蚀后的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%。