颗粒无序分布的堆积床内部流动与传热分析

为了详细分析单组元发动机催化床内部的流动与传热,针对颗粒无序分布的堆积床,提出了一种网格生成途径,适合于对大量颗粒自然堆积的管道流动进行数值模拟。根据该方法,对不同颗粒直径、不同进口速度的管道堆积床进行了液体介质的流动与非定常传热模拟,并与实验结果和Ergun方程计算进行了对比。结果表明,本文提出的方法可以有效模拟颗粒无序分布的堆积床内部流动与传热,适合于大量颗粒的堆积床仿真;网格尺度小于1/20颗粒直径时,模拟结果与实验结果符合很好;在相同雷诺数下,摩擦系数随颗粒直径增大而减小,Ergun方程在高雷诺数下计算的流阻偏大;非定常传热时,壁面效应会使壁面附近温度先接近流体温度。     

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

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

基于正激波位置计算的进气道/发动机系统实时仿真模型

研究了超声速条件下,基于进气道内正激波位置计算的进气道/发动机系统综合仿真问题。提出了将超声速进气道分为内、外两部分独立计算的方案,外部通过斜激波计算公式获得进气道内部的边界条件,对进气道内部采用准一维CFD模型模拟背压对正激波位置的影响。相比二维CFD计算,在保证参数精度要求的前提下有效缩短了计算时间,提高了仿真实时性。进一步将进气道模型与发动机部件级模型匹配运行,构建了进气道与发动机综合实时仿真模型。进气道/发动机综合模型的仿真表明,其能够准确模拟进气道正激波位置对涡扇发动机的动稳态影响,当来流发生扰动时,通过喷管喉道面积和压缩部件导叶角的快速调节,能够抑制正激波位置的变化,并使推力下降量减少50%。      

叶片丢失后发动机整机响应模拟试验与仿真

依据大涵道比涡扇发动机结构设计了叶片丢失试验台,开展了一系列模拟叶片丢失试验,并采用显式有限元方法进行了数值仿真,研究了发动机叶片丢失后整机结构响应与载荷传递规律。结果表明:丢失叶片与机匣存在叶尖与叶身两次撞击,对应的加速度曲线存在两个响应峰值;转子转速越高,加速度响应幅值越大。叶片飞断后转子不平衡载荷传递路径为前轴承支承-中轴承支承-中介支板-机匣结构;叶片撞击机匣导致的冲击载荷则由风扇机匣向后传递,最终传给吊装结构;发动机承受的载荷是由不平衡和冲击影响耦合得到,其中冲击载荷为主要部分。该研究为掌握真实发动机叶片丢失下整机响应规律提供了试验模拟方法与数值仿真分析工具。      

平流层电动机控制器热特性仿真分析和实验

为得到某电动机控制器在平流层低温低气压环境下工作时的热特性,开展了基于软件Icepak的热仿真分析,并通过地面环境模拟实验和平流层飞行实验对仿真模型进行了验证,仿真结果与实验结果基本吻合。结果表明:控制器的热设计可以满足保温和散热要求;在平流层低气压条件下,内部风扇可有效地保持控制器内部温度的均衡;平流层高度的太阳辐射对控制器内部温度分布有较大影响,在进行热设计时应予以考虑。研究结果为平流层电动机控制器的工程应用提供重要的参考依据。      

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

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

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）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。