壁面催化效应对高超声速气动热影响研究

高超声速飞行器壁面催化效应会导致激波层中原子在壁面处复合释热,加剧周围气动热环境。针对高超声速流动壁面催化特性,选择不同飞行马赫数及高度条件,采用完全催化和非催化两种条件对球锥模型壁面热流率进行数值模拟计算,研究壁面催化效应对气动热的影响规律。结果表明,固定飞行高度时,壁面催化效应对气动热的影响随马赫数增加而加强,Ma=25条件下驻点处完全催化与非催化热流比值高达1.92;固定飞行马赫数时,在50km高空以上壁面催化效应对气动热的影响随高度增加而减弱;壁面催化效应不仅会影响壁面附近的流场特性及组分分布状态,而且对整个激波层都有一定的影响作用。     

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

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

变循环发动机过渡态性能直接模拟方法

为了实现多变量可调的变循环发动机（VCE）过渡态性能模拟和控制规律设计,在稳态逆算法和过渡态性能隐式格式计算方法的基础上开发了变循环发动机过渡态直接模拟方法。在稳态逆算法模型中加入容积效应和转子动力学方程,实现了同时给定加速率、涡轮前温度和压缩部件工作点的条件下直接模拟计算过渡态过程中变几何参数和燃油流量的调节规律,验证了方法的精度和可行性,并将该方法用于变循环发动机转模态性能模拟和控制规律设计。计算结果表明,过渡态直接模拟方法的误差在0.58%以内,超声速巡航状态下由单外涵模态转换到双外涵模态的时间约为1 s,海平面静止状态下双外涵转单外涵的时间约为3 s,且推力、转速、涡轮前温度、喘振裕度等参数过渡平稳。该方法可简化转模态控制规律设计流程,并提高设计精度。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

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

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

横向矩形微槽对高超边界层失稳的控制作用

模拟了马赫数为6的空间发展平板边界层,通过在平板表面添加二维横向微槽研究了其对基本流及第二模态扰动波的影响。结果表明:所研究的二维微槽构形与常规多孔涂层相比具有较大的尺寸(100 μm以上),微槽对基本流的流向速度影响很小,空间采样点上的频谱并未出现新的不稳定模态。微槽有利于减小摩擦阻力,槽内流动导致的压差阻力比摩擦阻力低1~2个数量级,总阻力随着开槽率的增加而减小,当开槽率为06时,减阻效果超过40%。在入口以两种方式添加了第二模态扰动波,包括单个扰动和多个叠加扰动,下游的演化结果说明横向微槽能够在一个宽带频率范围内对第二模态的增长起到明显的抑制作用,且控制效果随着开槽率的增大而增强。     

基于实战使用的涡轴发动机空中起动飞行试验

以涡轴发动机实战使用中空中紧急起动相关技术指标和特情处置措施验证为目的,设计了贴近实战使用特征的发动机空中起动试飞方案,并以某涡轴发动机设计定型试飞为依托,进行了不同飞行工况下的飞行试验验证和数据分析研究。结果表明:该空中起动试飞方案合理可行,能够最大程度的贴近部队实战使用的技术特征,满足对发动机起动高度包线、起动时间等研制总要求规定的技术指标验证需求;试飞结果表明该型涡轴发动机4000m以上高度停车至再起动期间直升机平均下滑高度约为670m,平均停车时间约为185s。该试飞方法和数据结果可为部队实战使用过程中空中遭遇停车的紧急处置提供支持和参考。     

跨声速压气机转子的二次流旋涡结构

为了明确跨声速轴流压气机内部流场结构,数值模拟了NASA Rotor37转子,结合λ2准则分析流场参数,探索流动的规律和旋涡结构。研究发现,压气机转子的旋涡模型主要由马蹄涡、壁角涡、径向涡、脱落涡、泄漏涡、诱导涡和分离涡等7个旋涡组成。马蹄涡吸力面分支耗散,压力面分支向相邻的吸力面发展。壁角涡与脱落涡位于叶根角区,引起流动损失和角区失速。径向涡位于激波后吸力面的分离区内,它扩大吸力面分离、引起低能流体向叶顶堆积。激波与叶尖泄漏在叶顶通道中形成3涡:泄漏涡、诱导涡和分离涡,而叶栅通道出口存在分离涡和由泄漏涡与诱导涡合成的叶顶通道涡。泄漏涡与诱导涡破碎在流道中间产生的堵塞区,分离涡造成吸力面尾缘的低速区,共同触发跨声速压气机的失稳。     

大涵道比涡扇发动机TPS短舱低速气动特性分析

为了评估民机低速带动力试验时进排气效应的影响,选取大涵道比发动机涡轮动力模拟器(TPS)短舱和真实发动机短舱作为研究对象,采用数值模拟方法对其起飞、进近状态的低速气动特性进行对比分析。结果表明:由于TPS流量低于真实发动机需求,其唇口、外罩流场特征和真实发动机短舱有所不同,阻力特性也有差别;在进气道处于亚临界状态时,TPS短舱阻力系数比真实短舱大了约1.7个阻力单位,又由于唇口当地气流攻角更大,使得TPS短舱失速攻角相对降低了约1.0°;当进气道工作于超临界状态时,TPS短舱虽然也可以反映真实短舱的流动特性,但由于捕获流管收缩情况和气流驻点随攻角的变化,使得在0°~20°攻角时TPS短舱的阻力系数高于真实短舱,而在20°~30°攻角时其阻力系数略低,差量最大约为1.8个阻力单位。对于研究的大涵道比发动机,未经唇口及外罩修正的TPS短舱其低速气动特性基本可以反映真实进排气效应的影响,但在气动特性分析中可以考虑进一步修正进气效应的影响。