连续旋转爆轰发动机外流场三维数值模拟

为研究连续旋转爆轰发动机(CRDE)内外流场的变化特性,采用氢气-空气单步有限速率化学反应模型,对内径为40 mm、外径为60 mm、长度为50 mm的连续旋转爆轰发动机进行三维数值模拟,获得了CRDE内外流场结构特征和旋转爆轰波相关参数的变化特性,分析了不同进气总压条件对流场结构和发动机性能的影响。结果表明:爆轰产物在燃烧室出口附近膨胀加速,压力和温度大幅降低,在流场下游产生激波使压力回升,且随进气总压的升高,激波距燃烧室出口距离增加;出口附近羽流中心形成低压高温区域,中心平面上的平均压力低于环境压力,给发动机推力带来了副作用;羽流外围的空气受出口处斜激波的扰动,压力呈现出周期性变化;发动机推力随进气总压的升高而呈线性增加,进气总压为0.55MPa时,发动机推力达到了1160 N。计算仿真结果对掌握连续旋转爆轰发动机外流场特性具有一定的参考价值。     

径向入射激波聚焦实验和数值模拟研究

为得到不同入射激波马赫数对凹腔内径向入射激波聚焦过程的影响,对马赫数为1.25和2.18的径向入射激波在实验段内反射聚焦的过程进行了实验和数值模拟研究。用高速CCD拍摄了凹腔中气流流场的纹影图片,并采用保持强稳定性的龙格—库塔格式、加权基本无振荡算法和块结构动态自适应网格加密对激波反射聚焦的过程进行了数值模拟,实验和数值模拟结果吻合较好。复杂的激波形状及其之间的相互作用(如激波反射,激波聚焦和不同激波的相互作用)在实验和数值模拟中都能很好地观测到。通过比较不同马赫数入射激波的反射聚焦过程,发现较高的马赫数下聚焦后能产生射流,而较低的马赫数下却没有。同时发现,入射激波马赫数能影响二次激波出现的位置和强度,进而影响聚焦过程流场的紊乱程度。     

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

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

三维超声速不等密度平面混合层的小激波形成与演变

基于大涡模拟方法与高精度混合格式,对两种密度比(1.0与0.138)条件下的高对流马赫数(Mac=1.1)三维空间发展超声速平面混合层进行了数值模拟。数值结果清晰描述了扰动激励下超声速平面混合层的演变过程,讨论了三维混合层中小激波的结构及其形成过程。分析了混合层两侧密度变化对混合层失稳及小激波形成的影响,发现当上、下层密度不等时,混合层失稳加速,小激波强度减弱,而密度较低侧的小激波会更弱。     

CARDC 激波风洞 TSP 技术研究进展

从基本原理、关键技术和验证应用三个方面总结了近两年在中国空气动力研究与发展中心激波风洞中开展的温敏涂层(TSP)技术相关研究工作。通过解决快速响应温敏发光材料研制、模型研制、数据处理等一系列关键技术,完成图像采集系统、光学系统及标定系统的配套和系统集成,建立了一套适于激波风洞试验的高速 TSP 测量及标定系统。该技术可在激波风洞试验中获取模型被测面温敏涂层的发光图像,基于该图像可以直接观察模型表面热流分布和捕捉峰值热流的准确位置。结合温敏发光材料的物性参数标定数据,能够实现对模型表面热流的定量测量。不同于传统的传感器点热流测量技术只能得到模型表面有限数量的离散点的热流值,TSP 技术能够以高空间分辨率得到较大面积区域的详细热流分布信息,可更加全面的测量模型外表面的热环境,并且可以据此进一步分析和辨别边界层流态以及确定边界层转捩位置。试验对比表明,TSP 技术的测量结果与点热流传感器的测量结果具有良好的一致性。目前该技术已趋于成熟,在Φ2 m 和Φ0.6 m 激波风洞上成功应用于边界层转捩研究、局部干扰区热环境研究和复杂外形飞行器热环境研究等领域,已成为激波风洞除点测热技术之外又一重要测热技术。     

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

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

后掠翼三维鼓包串激波控制参数的影响

针对跨声速后掠翼,三维鼓包串作为一种有效的减阻方式具有结构简单、高效及鲁棒性好等优点.利用全局优化算法探索了鼓包设计参数空间的整体特性,并对鼓包长度、三维鼓包展向设计参数对鼓包减阻效果的影响进行了研究,发现鼓包顶点位置和高度对阻力系数最敏感,三维鼓包的展向设计参数则对阻力系数不敏感,而鼓包长度和鼓包相对展长越长越有利于减阻.在此基础上开展了小后掠角自然层流机翼加3种不同类型鼓包串的优化研究,通过优化结果发现,增加优化后的三维鼓包串,可将小后掠角自然层流机翼阻力发散马赫数向后推移,并且鼓包平均长度和控制区越大,效果越好.三维鼓包串具有良好的局部控制特性,可用于局部较强激波的抑制.三维鼓包串对常规后掠翼波阻具有良好的控制效果,同时能够抑制激波诱导的机翼后缘气流分离.      

导流块深度对抛物形凹面腔中激波反射聚焦的影响

为研究导流块深度对抛物形凹面腔内径向入射激波聚焦过程的影响,对马赫数为1.41的径向入射激波在导流块深度分别为0,5,10,15mm的凹面腔内反射聚焦过程进行了实验研究。结合高速CCD拍摄到的凹面腔中气流流场纹影照片和动态压力传感器测得的聚焦过程中流场的压力变化,对径向入射激波在凹面腔内的反射聚焦过程进行了描述。通过比较不同导流块深度下激波反射聚焦过程,发现随着导流块深度的加深,前导激波聚焦和反射激波聚焦的时间差会减小,使激波聚焦的强度增大,当导流块深度从0mm增加到15mm时,激波聚焦所致的峰值压力从0.39MPa增加到0.51MPa。但是,随着导流块深度的加深,压力增益下降,并使排气过程的难度加大,因此导流块深度为10mm左右时能取得较好的聚焦效果。     

预设堵块法检测进气道自起动能力的数值研究

在进气道内预先设置轻质堵块迫使进气道不起动,堵块被气流吹出后流道恢复畅通,为激波风洞提供了一种检测进气道自起动能力的方法。为了深入认识预设堵块方法的检测过程,将堵块简化为一个自由度的刚体运动,采用k-ωSST湍流模型,结合铺层动网格技术,对该检测方法进行了二维非定常数值模拟。在Ma_∞=5.9条件下,采用改变堵块质量的方式,获得了三种典型的检测过程。通过分析堵块在气动力作用下的运动与进气道非定常波系演化相互耦合的过程,揭示了堵塞作用产生和消除的流动机理。结果表明,堵塞时间没有显著地改变进气道的自起动能力,预设堵块方法成功检测到进气道在Ma_∞=4.9自起动。此外,建立了预估堵块运动的匀加速模型,预估的堵塞时间与数值模拟结果较为符合。在Ma_∞=5.9条件下,堵块的临界质量约为3.8g。考虑到激波风洞实验时间短暂,应选择合适的堵块。     

新型翅片式减涡器减阻特性数值研究

为了降低压气机径向引气过程中的压力损失,在设计出新型翅片单元结构的基础上,研究了新型翅片单元结构对径向引气压力损失的影响规律,对不同转速、新型翅片结构的去旋系统开展了数值研究,得到了不同工况下压气机共转盘腔径向引气的流场结构及压力损失分布曲线。研究结构表明:新型翅片单元结构能够抑制盘腔内气流旋流比,降低引气压力损失;翅片单元通道宽度和高度均存在最佳值使得减涡器减阻效果较好,在优选结构翅片单元通道宽度L=0.78,通道高度R3=0.97的条件下,其减阻效果较简单盘腔模型提高86.5%。高低翅片结构能起到较好的减阻效果,随着单侧翅片高度的升高减阻效果逐渐增强,在本文结构下增加单侧翅片高度L1=0.3时减阻效果最优,且A侧或B侧翅片增加带来的减阻效益相同。一方面,最优高低翅片结构其减阻性能相比于简单盘腔模型、典型翅片式减涡器模型以及翅片单元通道宽度L=0.78,通道高度R3=0.97的结构模型分别提高87.5%,29%,7.8%;另一方面,最优高低翅片结构能够减轻翅片单元的质量,具有较高的工程应用价值。