乙醇-空气稀液雾值班火焰的大涡模拟研究

为加深对气液两相湍流燃烧现象的认知,检测火焰面模型在液雾燃烧中的适用性,本文在欧拉-拉格朗日架构下使用火焰面/反应进度变量模型(FPV)数值模拟了乙醇-空气稀液雾值班火焰。欧拉坐标系下的气相湍流场使用大涡模拟方法模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑了相间质量、动量、能量交换。燃烧模型中采用碳元素定义混合物分数,在混合物分数方程源项中体现液相对燃烧模型的影响。模拟得到的气相温度分布和液相统计值均和实验数据较好吻合,验证了该燃烧模型对稀液雾扩散类型火焰的适用性。分析瞬时图发现,该稀液雾火焰的最高燃烧温度往往位于当量混合物分数附近,在出口下游20倍直径处火焰完全点着,此处上游FPV模型能给出局部点火熄火现象。蒸发作用在剪切层和点火区域较强,而液滴和火焰作用较弱,单一液滴很少被火焰包围。     

高超声速复杂气动问题数值方法研究进展

高超声速流场具有复杂流动特征,其中真实气体效应、磁流体干扰效应和力热结构耦合效应等对气动力分析产生了重要影响。将流体力学研究扩展到分子动力学、电磁流体力学以及流固耦合等交叉学科领域,这给数值模拟方法带来了巨大挑战。针对高超声速气动力/热分析的热点问题,重点关注高温效应与低密度流动效应、磁流体干扰效应和力热结构耦合效应等,结合算例分析了相应的数值求解技术;在气动热方面主要比较了3类求解方法(纯工程方法、纯数值方法和基于Prandtl边界层理论的方法),并给出了相应算例;对于气动力/热/结构耦合问题,从耦合模型及耦合计算方法两方面开展了分析。最后指出了高超声速复杂气动问题数值求解技术未来需重点关注的几个方面。     

高超声速飞行器表面温度分布与气动热耦合数值研究

针对高超声速飞行器热防护设计中的高温气体非平衡效应问题和气动热环境精确预测问题,基于流场的非平衡Navier-Stokes方程、表面的能量守恒方程和内部的热传导方程,考虑流场的非平衡效应、表面的热辐射效应、催化效应和烧蚀效应以及热防护层内部的热传导效应,建立了初步的表面温度分布与气动热的耦合计算方法,完善了高超声速飞行器气动物理流场计算软件(AEROPH_Flow)。在表面材料为碳-碳(C-C)条件下,对飞行高度为65km和飞行速度为8,10km/s的半球以及飞行高度为50km和飞行速度为8km/s的球锥模型,开展了表面温度分布与气动热的耦合计算,验证了计算方法和计算软件,分析了表面温度分布对气动热环境的影响。研究结果表明:表面温度分布对气动热的计算结果有较大影响,在气动热环境的预测中,不仅要考虑热化学非平衡效应和表面催化效应的影响,还要考虑表面温度分布的影响,最好是采用表面温度分布与气动热耦合计算的方法,以减小表面温度分布对气动热计算结果的影响。为此,需要发展完善非平衡流场/表面催化和烧蚀/热传导温度场(气/表/固)的计算模型、耦合求解技术和计算软件,实现对高超声速飞行器的真实飞行条件下高温气体非平衡效应和气动热环境的精确模拟。     

高超声速气动热预测技术及发展趋势

高超声速气动热预测技术是高超声速飞行器发展的关键技术之一。对高超声速气动热预测技术的发展情况进行了分析探讨。首先,简要回顾了国内外高超声速气动热理论预测及地面实验技术的发展历程;在此基础上,结合典型外形的计算与风洞试验结果的比较,重点介绍分析了气动热工程计算方法、数值模拟方法、气动热风洞试验设备的模拟能力及目前实验测试技术的研究水平;最后,对气动热预测技术的发展趋势进行了讨论,提出了气动热预测技术应研究解决的问题。     

发动机动力模拟风洞试验中的空气桥技术

介绍了发动机动力模拟风洞试验中的空气桥设计技术和影响修正方法.通过自由度分析,选择了合理空气桥布局,采用有限元方法对空气桥关键受力梁进行了优化,建立了空气桥天平一体化数值模拟技术,使得空气桥和天平刚度更加匹配.通过这些设计及优化使得空气桥的作用力最小,并具有较强的克服压力影响、温度影响的能力.通过试验建立了空气桥附加刚度影响、内部压力影响、温度影响和内部流量影响的修正方法,进一步减小了空气桥对天平的影响.在8m×6m低速风洞进行了某大展弦比飞机全模涡扇动力模拟器(TPS)短舱动力模拟试验.试验结果重复性好,阻力系数精度达到0.0003,和相关文献吻合.这表明空气桥技术是成功的,满足了TPS短舱动力模拟试验要求.      

7150铝合金热变形行为及微观组织

采用Gleeble热模拟机进行热压缩实验,研究7150铝合金在变形温度为300~450℃、应变速率为0.01~10s-1条件下的变形行为,采用Zener-Hollomon参数法构建合金高温塑性变形本构方程,并对变形后的微观组织进行分析。研究表明:7150铝合金的流变应力随应变速率增大而增大,随变形温度增大而降低。该合金热压缩变形的流变应力行为可用双曲正弦形式的本构方程描述,其参数A为4.161×1014s-1,α为0.01956 MPa-1,n为5.14336,热变形激活能Q为229.7531k J/mol。随着温度升高和应变速率降低,动态再结晶逐渐取代动态回复成为合金的主要软化机制。     

基于微观结构分析的CMDB推进剂动态断裂韧性加载率敏感性研究

为研究改性双基推进剂(CMDB)的动态断裂特性,利用霍普金森实验技术(SHPB)对CMDB推进剂进行了冲击实验,并对实验数据的可靠性进行了校验;通过数值仿真的方法进一步验证了实验数据的有效性,获得了CMDB推进剂的动态断裂韧性;利用电镜扫描技术对试件断面进行了微观表征和分析。结果表明:CMDB推进剂在加载率58906~100056MPa·m1/2·s-1范围内,表现出明显的加载率敏感性;加载率达到105000 MPa·m1/2·s-1左右时,断裂韧性不再继续增大,反而出现下降趋势;电镜分析结果表明,CMDB推进剂加载率敏感性与其内部高氯酸铵(AP)颗粒的微观破坏机理直接相关。     

基于DDES的有／无舱门腔体气动噪声仿真分析

采用延迟分离涡仿真（ DDES）的方法研究了空腔和带有舱门的腔体的气动噪声特性。仿真在自由来流M＝0．85下,长宽深比值为5：1：1的矩形开放腔体展开。空腔仿真的声压级幅值结果与试验数据相比误差在4％以内,频率位置的预测较好。 DDES方法对于气动噪声问题具有很好的仿真准确性。空腔仿真结果和带有舱门的结果对比后表明,舱门迫使剪切层气流更多的撞击后壁进入腔体内部,使腔体内部尤其是中部压强波动增大,声压级在低频部分增量显著,达到了5 dB左右。     

复合材料副翼典型结构件的VARI工艺模拟及试验验证

在复合材料液体成型工艺过程中,通过计算机软件的模拟研究,可以预先判断树脂在增强体中的流动形式和获得树脂的注模时间以及注模压力,进而为实际成型工艺提供参考设计。本文利用PAM-RTM模拟软件对复合材料副翼典型结构件的树脂流动过程进行了模拟计算。根据达西定律推导出树脂在增强纤维中一维流动渗透率计算公式,分别对增强纤维、缝线纤维进行了渗透率测试,在此基础上对渗透率数值进行了修正。利用渗透率修正值采用不同的流道设计方案对树脂注模过程进行了模拟研究,根据模拟结果最终选择合理的注射方式和流道布局。对所选方案进行了工艺试验验证,并对结构件进行了超声无损检测。结果显示:理论充模时间和试验充模时间基本一致,结构件没有内部缺陷,树脂填充完全,模拟工艺计算有效。     

二元双喉道射流推力矢量喷管流动参数影响的数值研究

采用数值模拟方法研究了不同流动参数对二元双喉道射流推力矢量喷管(Dual-throat fluidic Thrust-vectoring Nozzle,DTN)内流特性和推力矢量控制效果的影响。结果表明,DTN在非推力矢量情况下,NPR在3~4范围时,推力系数较大,达到0.968,而流量系数较小,仅为0.93;NPR再增大,推力系数迅速下降。在推力矢量情况下,落压比一定时,随着次流流量比的增加,推力矢量角增加,而流量系数、推力系数、推力矢量效率减小;次流流量比一定时,随着落压比的增加,推力矢量角减小,系统推力系数先增加后减小,流量系数略微增加。