局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。     

内吹式襟翼环量控制翼型升力响应特性

传统尖尾缘翼型通过控制迎角,综合利用襟翼、缝翼来改变升力,升力对迎角变化的时间响应历程可以用Wagner函数来描述,而内吹式襟翼（IBF）主要通过控制分离来拓展最大升力,并在一定范围内通过调节射流强度改变驻点位置和环量来对升力进行有效控制,其升力随吹气动量变化的时间响应尺度是否与传统尖尾缘翼型相同还不是很清楚。本文主要研究内吹式襟翼升力响应过程,并将其与传统尖后缘翼型升力响应特性进行对比。首先通过某襟翼偏角为30°的双圆弧环量控制翼型对数值方法进行验证,再对某最大厚度为18%弦长的亚声速翼型内吹式襟翼定常吹气控制下的流场进行非定常数值模拟,并分析了其中的瞬态特征。结果表明内吹式襟翼环量控制翼型对激励响应的时间依赖特征与Wagner函数有很好的相互关系,并可以用该函数来描述。     

电弧风洞加热条件下聚四氟乙烯烧蚀热响应特性

高温加热条件下,由于聚四氟乙烯的热解,对烧蚀温度场计算结果有较大影响,为了提高聚四氟乙烯烧蚀温度场计算精度,建立了聚四氟乙烯烧蚀温度场计算方法,对电弧风洞加热条件下聚四氟乙烯表面烧蚀热响应特性进行了验证研究。理论计算和试验测量结果对比表明:230～323℃升温区间内,随时间增长,温度逐渐升高,理论计算与试验测量结果变化趋势一致;323～680℃升温区间内,随时间增长,试验测量温度逐渐升高,理论计算温度为定值,理论计算与试验测量结果存在一定偏差;680～390℃降温区间内,随时间增长,温度降低,理论计算高于试验测量值,这与理论计算烧蚀量存在偏差有关。采用聚四氟乙烯材料烧蚀温度场计算方法,可以有效模拟高温加热条件下聚四氟乙烯热响应特性,从而为产品设计提供参考。     

刮削作用对涡轮转子机匣通道涡形成的影响

采用数值方法求解耦合剪切应力输运(SST)湍流模型的雷诺平均Navier-Stokes方程组,研究了不同间隙尺寸下GE-E3高压涡轮第一级转子内刮削作用对机匣通道涡形成和发展的影响。通过与轮毂通道涡结构的对比,发现机匣通道涡层次结构与经典二次流理论存在明显的差异,并对该差异形成的原因进行了深入探索。结果表明:叶顶对机匣边界层的刮削作用在机匣通道涡的形成过程中占主导作用,刮削作用使得流向叶片吸力面的来流机匣边界层在交汇点区域从内层向外层卷起,形成层次结构相反的机匣通道涡;叶顶浸入比值是影响转子机匣通道涡形成的重要参数,随着比值的增大,机匣通道涡损失先增大后减小;只有在间隙尺寸较大情况下,叶尖间隙的"抽吸作用"才能抑制机匣通道涡的发展。     

高速气流冲击柱状泡沫多孔体的传热特性分析

采用流体/多孔区域一体化单区域算法,数值研究了高速绕流条件下前置于圆柱体前缘表面的柱状泡沫多孔体内部的传热特性。基于蒙特卡罗法考虑多孔域内的辐射热效应,分析了变化多孔区域长度和多孔阻力特性对模型激波阻力和前缘多孔区域气动热的影响。结果表明:在圆柱体前缘安置一定长度及带有适当阻力特性的泡沫多孔材料,可同时减小整体激波阻力并降低前缘表面的气动热效应。在模拟工况下,无量纲长度1.0、黏性阻力系数0.2×107m-2及惯性阻力系数200m-1的前缘泡沫多孔可减小激波阻力13.5%,降低约75%的前缘表面的平均气动热流密度。保持无量纲长度不变,减小泡沫多孔区域惯性阻力系数会降低激波阻力,但会略微增加前缘壁面气动热流密度。      

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机,针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响,对全机带动力风洞试验模型进行了大规模并行非定常数值计算,再现了风洞试验现象,通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂,这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速,为保证飞行安全,继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示,在飞行雷诺数和螺旋桨转速下,相同侧风范围内风洞试验显示的流动不稳定因素基本消失,偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析,建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系,进行了偏航稳定性的虚拟试飞评估,研究成果也为AG600飞机的首飞和飞行试验所验证。     

液体火箭发动机推进剂泵诱导轮与离心轮的匹配

为获得诱导轮离心轮周向匹配的时序效应对离心泵外特性以及压力脉动的影响规律,阐释相关作用机制,采用基于分离涡仿真（DES）的离心泵三维全流道数值仿真方法,引入熵产理论以及压力脉动强度系数等先进分析方法对不同匹配角度下离心泵内能量损失机制及压力脉动特性进行了研究。结果表明：离心轮诱导轮的时序效应对泵外特性有一定的影响,随着匹配角度的增加,扬程和效率均呈现先减小后缓慢增大的趋势,扬程变化为0.8%,效率变化为1.2%,其影响机制由不同匹配角度下叶轮通道分离涡、叶轮叶片尾迹以及靠近隔舌处扩压器通道回流涡变化决定;时序效应对离心轮扩压器动静干涉效应影响显著,当诱导轮叶片尾缘位于离心轮相邻主叶片中间位置时,能够有效消除3倍频成分,显著降低泵内压力脉动水平,其中动静干涉区域以及隔舌处扩压器叶片表面压力脉动平均降幅分别达到14.5%和16.7%。     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

含有SMA弹簧驱动器的可变倾斜角翼梢小翼研究

 针对传统翼梢小翼在非设计状态减阻效果不佳的缺点,提出一种含有形状记忆合金(SMA)弹簧驱动器的变体翼梢小翼结构,它能根据飞行状态主动调整小翼的倾斜角,实时优化飞机的阻力特性.采用力-热-应变耦合法设计了所需的SMA弹簧驱动器,并通过有限元仿真与风洞试验验证了变体翼梢小翼的变形能力,最后初步研究了变体翼梢小翼的闭环控制方法.研究结果表明,在飞机的起飞阶段(自由来流流速为26 m/s,迎角为3°),变体翼梢小翼的倾斜角能在1 min内自主完成预定变化过程,倾斜角的最大变化量为23°,控制精度的最大误差为12%,各项指标均符合设计要求.     

大攻角地效翼非定常流动实验研究

大攻角地效翼非定常流动研究对大攻角地效翼非定常空气动力特性和流动控制具有重要的意义。对攻角为18°的NACA0012地效翼在模型风洞中进行实验,利用三维热线风速仪对地效翼后缘附近尾迹区内测点的速度脉动进行采样,分析不同飞行高度下大攻角地效翼的非定常流动特性。结果表明：地面效应下,大攻角地效翼后缘涡系形成位置后移;随着地效翼高度的降低,非定常流动的脉动频域宽度减小,速度脉动特征频率减小。