复合固体推进剂细观结构建模及脱黏过程数值模拟

为研究复合固体推进剂损伤演化规律,基于分子动力学颗粒填充算法构建了HTPB(hydroxyl terminated polybutadiene)推进剂细观结构模型,通过在AP(ammonium perchlorate)颗粒/HTPB基体界面处引入黏接接触替代传统的黏接单元,并基于Hooke Jeeves的参数优化算法反演得到颗粒/基体界面处内聚力模型参数,利用双线性和自定义指数型损伤内聚力模型模拟了AP颗粒和HTPB基体黏接界面处损伤的萌生、发展、聚合直至宏观裂纹破坏的过程。通过数值仿真与实验结果对比发现,指数型损伤内聚力模型比双线性模型能更准确描述推进剂单轴拉伸过程中颗粒与HTPB基体界面间脱黏过程。最后对比了多阶段加载实验结果与仿真结果曲线,发现两者变化趋势基本一致,最大偏差仅为10%,验证了所建细观模型的可靠性及反演所得界面参数的准确性。      

非轴对称端壁设计的高负荷涡轮气热性能研究进展

非轴对称端壁设计因能够有效地减少涡轮叶栅的二次流损失和提高气动性能而在高负荷涡轮设计中得到应用。简要回顾了涡轮叶栅二次流模型和非轴对称端壁造型方法,重点综述了非轴对称端壁设计的高负荷涡轮气动性能研究进展和抑制端壁二次流的作用机制,介绍了非轴对称端壁设计的高负荷涡轮端壁气热耦合作用的冷却特性研究进展,总结了高负荷涡轮非轴对称端壁设计技术的应用成果,展望了非轴对称端壁设计在高负荷涡轮的高效气动和冷却布局应用方面需要深入研究的内容。     

1 种副孔正交型超声速5 孔探针的设计与应用

针对超声速喷管出口流场Ma和方向角的测量,设计了直径为4 mm的副孔正交型超声速5孔探针。圆锥形测压头半角为20°,测压孔直径为0.4 mm,探针直段长140 mm。分别在亚声速风洞(196个校准点)和超声速风洞(294个校准点)中对5孔探针进行了校准,结果表明:5孔探针在不同Ma下的方向特性曲线具有几何相似性,有较高的角度灵敏度;总压特性系数和静压特性系数均具有良好的对称性。并通过数值计算研究了5孔探针的扰流特性。结合自动位移机构,完成了超声速喷管出口5个截面共315个坐标点的标定,获得了流场的Ma和方向角分布特性。     

融合式翼梢小翼减阻效应研究

有效减小诱导阻力对于飞机降低油耗、提高航程具有重要意义。针对某飞机翼身组合体构型,采用CFD数值模拟方法分析融合式翼梢小翼对飞机气动力特性的影响,尤其是其减阻效应;并给出翼梢小翼附近的空间流场。结果表明:带翼梢小翼后翼尖涡强度减弱,飞机阻力系数明显下降;固定升力系数0.5时,弯矩增加3.2%,阻力系数减小4.2%。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

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

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

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。     

航空发动机两级反向旋流燃烧室燃烧流场大涡模拟研究

为深入了解真实航空发动机内燃烧流场,采用动态亚网格模型结合单步快速化学反应(FC)、火焰面(FLM)及反应进度变量(FPV)等三种燃烧模型对径向两级反向旋流燃烧室单头部构型进行了大涡模拟。为避免模型简化误差,对燃烧室包括全部气膜冷却孔在内的精细结构进行了完全仿真。在已达到统计定常状态的冷态流场基础上首先模拟了燃料喷注和掺混过程,约2.6ms后燃料到达真实的点火位置,随后采用FPV模型在半径3mm的球形区域数值模拟了点火,展示了在主燃孔横向射流作用下初始火焰沿化学恰当比混合分数等值线传播并充满整个火焰筒的发展过程,结果显示火焰到达燃烧室出口的耗时约为6～7ms。不同模型算法预测的平均温度场与CARS测量结果作了对比,LES-FPV,RANS-FPV,LES-FLM以及参考文献中RANS-FLM计算平均误差分别为3.47%,4.17%,7.76%和11.22%,表明LES改进了模拟精度,且FPV模型显著优于FLM模型。RANS-FPV预测的出口存在严重热斑,导致其给出的出口温度分布因子(OTDF)及最大径向温度分布因子(RTDF)值分别达到0.593和0.313;LES-FPV结果均匀性最好,其预测值分别为0.284和0.193。     

电子回旋共振等离子体推力器放电机理数值模拟研究

电子回旋共振等离子体推力器(ECRPT)是一种高比冲、高效率且结构简单的新型电磁式推力器。为了研究推力器的放电原理和工作机制,采用漂移-扩散流体模拟方法,仿真模拟了微波等离子体放电过程。仿真结果表明,电子数密度达到10~(16)~10~(17)m~(-3)数量级,氙气的电子数密度比氩气高50%;电子数密度、碰撞功率损耗均随着计算域内压强的增大而增大,电子温度随压强的增大而减小;电子数密度、碰撞功率损耗随着入射微波功率的增大而增大。在未来ECRPT的实际应用中,可以通过使用氙气,适当增大推力器腔内压强以及入射微波功率,使其具有最佳的推力、比冲和工作效率。     

壁面粗糙度对高超声速进气道气动性能的影响

为探究壁面粗糙度对于高超声速进气道气动性能的影响,采用经过校验的数值仿真和理论分析相结合的方法进行研究。结果显示:同一飞行马赫数下,随着壁面粗糙度的增加,进气道的流量系数、总压恢复系数和出口马赫数逐渐降低,而静压比、压差阻力系数、摩擦阻力系数以及起动马赫数则逐渐增大;不同飞行马赫数/不同钝化半径下,进气道性能参数随壁面粗糙度的变化规律相似,均表现出较好的拟合规律,据此获得的拟合公式及光滑壁面进气道的气动性能可预估不同壁面粗糙度下进气道的气动性能;就本文研究的进气道而言,当壁面相对粗糙度从0增加至0.625%时,进气道起动马赫数从4.25增加至4.85。壁面粗糙度增加,导致进气道沿程附面层增厚是进气道气动性能参数出现上述变化规律的主要原因。