某氢氧发动机推力室氢喷嘴烧蚀问题仿真分析

氢氧火箭发动机推力室内的喷注均匀性不但影响燃烧效率,还有可能影响喷嘴、面板及内壁等结构的可靠性.针对某型氢氧火箭发动机推力室多次出现固定位置两个氢喷嘴的烧蚀问题,采用CFD方法模拟了此发动机氢头腔及喷嘴的内部流动.通过分析流动特性,并给出量化对比结果,得出了以下两个结论:一方面此发动机推力室喷注面氢流量分布不均,而多次产生烧蚀的喷嘴是所有喷嘴中氢流量最小的两个;另一方面氢喷嘴出口环形间隙内流量分布不均,在所有喷嘴中产生烧蚀的喷嘴出口流速分布不均匀度是最高的.这两个因素共同作用下导致喷嘴局部混合比过高,是造成固定位置喷嘴局部烧蚀的重要原因.     

轨控发动机真空流场计算

拦截弹轨控发动机在真空中的流场可为其设计提供重要的理论依据，同时其喷焰红外辐射特性也是防御上的重要研究对象。采用基于有限体积形式的LU格式离散N-S方程，通过时间推进法求解拦截弹轨控发动机喷管以及外场喷流区域在内的气相统一流场，同时考虑了各主要组分参与的化学反应，得到了轨控发动机喷管内外速度、温度、密度、组分浓度等参数的分布情况。研究表明：使用本文中的方法可以很好地计算出轨控发动机在真空中的内外流场。真空羽流膨胀迅速。     

飞行高度对固体火箭两相喷焰流动与辐射的影响

固体火箭发动机喷焰的辐射特性对火箭目标探测与跟踪识别具有重要工程意义。本文针对某模型固体火箭发动机,基于欧拉离散相模型描述气-固两相相互作用,采用详细化学反应机理计算喷焰复燃效应,基于逐线积分法+米散射理论求解气体及粒子辐射物性参数,通过视在光线法计算辐射传输,并利用测量数据验证了本文模型的适用性,仿真分析了不同飞行高度对气固两相喷焰流动及辐射特性的影响。研究结果表明：随着飞行高度增加,喷焰中燃气与固体颗粒间相互作用减弱,且粒径越大的粒子与燃气间差异越大;两相喷焰受不同飞行高度上掺混与复燃效应程度差异的影响,致其温度分布和不同组分浓度分布存在显著差异;各高度下两相喷焰光谱辐射呈现气体选择性发射谱结构,其中高温Al2O3对辐射贡献主要集中在短波,且高度越高,影响越小;不同谱带间辐射特性受不同组分发射带影响而存在差异,且2.7μm和4.3μm两个主要发射带的辐射峰值出现在不同飞行高度。     

机匣法兰螺栓连接建模和密封性能分析

由于燃烧室机匣法兰螺栓连接处存在漏油问题,影响发动机的使用寿命和工作性能,还会污染环境、浪费油料。基于悬臂梁挠度变形理论,采用解析法推导螺栓连接结构挠度变形方程,建立优化连接结构密封设计的方法;在 ANSYS 中建立机匣螺栓连接精细有限元模型,计算法兰结合面在不同载荷形式下的轴向间隙分布;采用加强连接密封性的实施方案,通过数值模拟分析不同厚度密封圈的封严效果。结果表明：本文推导的方程能计算出法兰连接结构受载变形量,适用于连接结构密封优化,采用的方案能够解决实际工程应用中机匣法兰螺栓连接处漏油问题。     

针栓喷注器中心推进剂偏转角模型分析研究

为了实现针栓喷注器中心推进剂偏转角的准确预测,基于流场分析建立了中心推进剂偏转角理论模型。从动量守恒方程推导了中心推进剂偏转角公式,通过数值仿真和试验结果对其进行验证,并分析了工况参数和结构参数对中心推进剂偏转角的影响规律。结果表明：理论模型预测值与数值仿真和试验结果很好地吻合,套筒遮挡喷注面积对偏转角影响最大,在变推力时偏转角随着套筒遮挡喷注面积增加而减小。喷注压降、中心筒壁厚和底部凹腔深度对中心偏转角影响很小,当套筒遮挡喷注面积一定时,中心筒底部有凹腔的偏转角比没有凹腔的偏转角约大6°,该模型为针栓喷注器工程设计和进一步精确计算变推力下的雾化角提供了重要参考。     

双螺杆固液混合燃料输送系统设计与实验研究

固液燃料混合推进是航天推进系统的重要环节,液体发动机燃料热熵低,为了提高火箭或导弹的续航能力,结合固体粉末燃料热熵高的特点,采取了固液混合燃烧的策略,进而提高发动机燃料的热熵,增加发动机的续航里程。本文根据齿轮啮合原理,对关键结构双螺杆齿形展开修形设计,建立了流道三维模型和仿真模型;在齿形设计的基础上,围绕固液混合均匀、变比例输送开展固液混输系统集成设计,完成整机的加工、装配、调试,搭建实验平台,建立液压测试系统,开展实验。结果表明,设计的两相流泵实现了固液两相物料的混合输送,实现了固液燃料变流量、变比例、变出口压力的均匀混合效果。该小型燃料泵能够实现出口流量30 ~ 600L/h、出口压力0 ~ 900kPa、固液混合比例0% ~ 10%的固液混输,且混输效果良好,研究结构可供设计变推力火箭发动机参考。     

CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes （RANS）方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会（DPW Ⅵ）。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明：迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。     

基于数值模拟的NSDBD等离子体激励器防冰特性

飞行器表面在一定气象条件下会产生积冰,积冰会使飞行器气动性能下降,是危害飞行安全的重要因素之一。常见的气热及电热防冰系统已经广泛运用于现有飞行器上。近些年,在纳秒脉冲阻挡介质放电（NSDBD）等离子体激励器的相关研究中发现NSDBD等离子体激励器可对周围流场进行快速加热,考虑到这种热效应可能作为飞机防冰的一种新方式。本文用数值方法对NSDBD等离子体激励器防冰特性开展了研究。首先,建立了基于Messinger模型的积冰模型,对典型积冰条件进行了验证计算;其次,耦合唯象学等离子体模型与非定常雷诺平均Navier-Stokes方程,计算等离子体对空气流场的影响;最后,将NSDBD等离子体激励器布置在NACA0012翼型前缘防冰区,结合积冰模型与唯象学等离子模型,对其防冰特性进行了研究。计算结果表明等离子体加热的热气流会覆盖在翼型表面防冰区。在相同的霜冰条件下,开启等离子体激励器时机翼前缘没有出现积冰,说明等离体子激励器应用于机翼防冰是有效的。针对不同的激励器参数对防冰特性的影响规律进行了研究,总体上防冰效果与峰值电压、激励器频率有关,从防冰效果和能耗方面考量,在给定计算条件下,存在最优电压值和最优激励器频率值。激励器分布方式对防冰特性的影响与其具体流场有关,需要具体分析。     

DDES方法在复杂旋翼流场计算中的应用

应用延迟脱体涡模拟（DDES）方法,进行了复杂分离流动旋翼流场的数值模拟研究。为控制计算规模,更好地模拟旋翼桨叶近壁面附近区域流动分离及远场尾迹发展,采用了贴体与背景自适应直角网格嵌套的网格系统;同时为提高计算速度,采用了基于分布式存储的数据并行模式,实现了贡献单元搜索、洞切割和流场计算的并行化。计算工作包含一个共轴双旋翼悬停和单旋翼下降状态,目的是分别讨论DDES方法在旋翼表面附近空间区域是否存在大范围分离时计算结果与雷诺平均Navier-Stokes结果之间的差异。计算表明,对于下降状态,DDES方法预测出了更强的尾迹区分离流动;对于桨叶表面基本为附着流动的双旋翼中等桨距角情况,两种方法预测的流场结构相似,仅在下旋翼桨叶气动力和旋翼桨根下方空间区域的诱导速度计算结果上存在微细差异。     

考虑榫槽/榫齿配合间隙的叶盘结构疲劳寿命稳健性优化研究

为了更好地满足榫齿结构的设计、加工和制造需求,同时提高叶盘结构的疲劳寿命稳健性,利用有限元方法对不同配合间隙下的应力应变进行计算,揭示并分析了叶盘结构疲劳寿命随配合间隙的变化规律,综合考虑载荷、材料参数、配合间隙不确定性的情况下,利用二次多项式响应面分别建立了随机变量与叶盘疲劳寿命的近似函数关系,并结合多目标规划理想点法建立了叶盘结构疲劳寿命多目标稳健性优化模型。优化结果表明:叶片、轮盘疲劳寿命均值分别增加了3.24%和1.93%,疲劳寿命概率区间分别降低了10.13%和8.16%,叶盘结构疲劳寿命对参数变化的敏感度降低,有利于更加准确地进行寿命评估。