轴向通流旋转盘腔流动换热的数值研究

为了探究旋转盘腔内的流动和换热规律,对轴向通流旋转盘腔进行了非稳态数值模拟,将计算结果与实验数据进行了对比,探究了流动不稳定性的发展过程,分析了盘腔内流动结构和盘面换热特性随旋转雷诺数的变化规律。结果表明:旋转引起的正旋涡从盘罩附近开始发展,随转速的增大而变大,挤压低半径区域的强迫对流区,最终扩展到整个盘腔,盘腔中轴面的涡对数与流动不稳定性的强度有关。上游盘和下游盘的高半径区域换热强度随转速的增大而增强,下游盘低半径区域的换热强度在低转速下由于冲击作用而较强,但该冲击作用随转速的增大而减弱,低半径区域的换热强度也就随之减弱。当旋转雷诺数增大到4.94×105时,下游盘低半径区域受到的冲击作用减小到可以忽略。     

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段;在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波;第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波;稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s;斜激波的存在使燃烧室出口平面流场产生了较大波动。      

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

航空紧固件生产车间流水线平衡量化分析及改善

为了量化地评价航空紧固件生产流水线的平衡性并对其进行改善,对一种典型的航空紧固件生产流水线进行了分析,结合流水线平衡的概念和航空紧固件生产流水线的特点,通过公式推导建立了一套可以从不同角度量化评价流水线平衡性的指标体系,阐述了其应用的方法和可行性;通过对比几种当前主要的研究流水线平衡的方法,指出航空紧固件生产流水线的平衡更适合采用工业工程的方法进行分析;研究了平衡性改进问题的解决思路,指出瓶颈工序对流水线的平衡性有重要影响,通过缩短瓶颈工序的加工时间可使整条流水线的平衡性得到有效改善,同时也使流水线上各工序的加工效率获得显著提升。最后通过实例计算对以上理论和方法进行了应用和验证。     

一种双余度非相似容错计算机软件设计

非相似余度设计是容错技术的重要设计方法,在很多领域都有广泛应用,在航空领域主要应用于飞行控制系统.在非相似双机架构下分析现有且适宜的操作系统,并在此基础上构建适合机载飞控系统的底层运行框架和任务调度算法,在此框架下研究双机余度管理的关键技术.研究结果表明,所设计的软件框架与算法满足飞控系统要求.     

飞机结冰冰形测量方法研究进展

不同气象环境下飞机部件的结冰外形是飞机结冰研究不可缺少的内容,它对开展结冰气动分析、结冰防护设计、结冰飞行操作和结冰适航取证等研究具有重要意义。结合国内外冰形测量技术发展现状,从接触测量方法和非接触测量方法两个方面分别介绍了现有冰形测量技术的操作流程及测量原理,系统分析总结了各自存在的优缺点。在此基础上,结合飞机结冰和结冰风洞试验的特点,归纳了冰形测量技术下一步发展面临的挑战,包括全类型结冰的精确测量、结冰全貌三维测量和结冰生长过程实时测量等,同时,从基础测量手段与数值计算融合发展的角度,展望了冰形测量方法未来的发展趋势。     

吸气式高超声速飞行器非均匀尾喷流试验

在中国空气动力研究与发展中心0.5m高超声速风洞中,开展了非均匀喷流条件下的吸气式高超声速飞行器后体尾喷流/外流干扰测压试验研究。采用非均匀内喷管,模拟飞行器尾喷管非均匀入流,测量了飞行器后体膨胀面及水平翼表面压力,采用高清纹影观测了喷流干扰区域的流场结构,获得了不同工况下非均匀入流对尾部及水平翼表面压力分布的影响规律。试验结果显示尾喷管非均匀入流对飞行器尾部壁面压力分布及流场结构有明显影响,喷管入流的非均匀特征在吸气式高超声速飞行器喷流模拟中不可忽视。非均匀喷流核心区压力分布明显高于均匀喷流时的结果;核心区域外,非均匀喷流的作用面积略小于均匀喷流,且非均匀喷流同外流交叉干扰区域的面积和强度要略小于均匀喷流;均匀喷流在喷管出口区域存在明显的膨胀波系,交叉干扰激波及剪切层的扩张角也大于非均匀入口条件时的结果。     

氧化层非均匀增长对热障涂层应力分布的影响

利用扫描电镜对不同氧化程度的热障涂层(TBCs)进行观察,分析发现热生长氧化层(TGO)各部位并非均匀增长。假设热生长氧化层的形状为正弦曲线,用统计方法处理扫描电镜图片可以得到热生长氧化层的基本几何参数。由此建立单波形二维轴对称有限元模型,计算并分析了非均匀增长模型的应力状况。结果表明:相同温度下,非均匀增长模型的轴向、剪切和Mises应力分量的变化趋势与均匀增长模型的对应趋势一致;但是非均匀增长模型的各个应力分量最大值均小于均匀增长模型的应力分量最大值。      

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。