热冲击后气膜冷却叶片结构强度流固热耦合仿真

研究了冷气流量对气孔周围热应力的影响,为气膜冷却叶片可靠性设计提供参考。改变气孔的孔径,并建立有限元模型,结合有限元/边界元理论,通过流固热三场耦合技术获得热冲击后的叶片最大温度、温度不均衡程度及最大热应力。研究表明:增加冷气量有利于改善叶片冷却效率降低叶片温度,但也会使叶片温度不均衡程度增加,加剧尾缘气孔内的热应力载荷;增加前缘气孔直径可提升66%的平均冷却效率,有利于减缓气孔内的热应力,增加尾缘气孔的直径对冷却效率及热应力的影响均较小。此外,数值计算结果与试验及解析解较为吻合,对于气膜冷却叶片结构设计具有参考价值。      

高速气流冲击柱状泡沫多孔体的传热特性分析

采用流体/多孔区域一体化单区域算法,数值研究了高速绕流条件下前置于圆柱体前缘表面的柱状泡沫多孔体内部的传热特性。基于蒙特卡罗法考虑多孔域内的辐射热效应,分析了变化多孔区域长度和多孔阻力特性对模型激波阻力和前缘多孔区域气动热的影响。结果表明:在圆柱体前缘安置一定长度及带有适当阻力特性的泡沫多孔材料,可同时减小整体激波阻力并降低前缘表面的气动热效应。在模拟工况下,无量纲长度1.0、黏性阻力系数0.2×107m-2及惯性阻力系数200m-1的前缘泡沫多孔可减小激波阻力13.5%,降低约75%的前缘表面的平均气动热流密度。保持无量纲长度不变,减小泡沫多孔区域惯性阻力系数会降低激波阻力,但会略微增加前缘壁面气动热流密度。      

三维并联式TBCC发动机排气系统设计与实验

为实现三维并联式TBCC排气系统的设计,利用基于二维最大推力理论的密切法设计了圆转矩变截面非对称冲压发动机喷管,并在冲压发动机喷管上壁面为三维曲面的条件下,利用转动加滑动共用面的调节方式实现了涡喷发动机喷管的内型面设计及其喉道面积的调节,从而完成了排气系统设计。随后对设计结果进行了冷流实验验证。结果表明,这种排气系统设计方法具有可行性。在冲压单独工作状态,壁面压力分布的数值计算与实验结果吻合较好,相对方均根偏差低于10.8%;在涡喷单独工作状态及共同工作状态,上壁面压力在激波附面层干扰区域内偏差较大,最大相对偏差可达30.6%。但壁面压力分布的数值计算与实验结果的整体变化趋势一致。      

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符;反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。      

大涵道比风扇/增压级进气畸变数值模拟研究

为探究进口总压畸变下大涵道比风扇/增压级内部流场的主要特征,基于三维彻体力模型的思想,开发了一套能够实现风扇/增压级内外涵联算功能的三维数值计算程序。利用该程序模拟了某大涵道比风扇/增压级在周向总压畸变进气下的三维流场。计算结果表明：大涵道比风扇单转子不同叶高处的畸变传递特征差异较大,转子出口总压畸变强度由叶根到叶尖逐渐降低,在叶尖处衰减为最小值1.5%;在转子出口相应诱导出的总温畸变强度由叶根到叶尖逐渐升高,在叶尖处达到最大值1.4%;进口周向总压畸变导致风扇转子总压比下降0.5%,而风扇转子出口形成的总压总温复合畸变导致增压级总压比下降2%;总压畸变在增压级中呈逐级衰减趋势,而高温畸变区的周向范围在逐级增加。     

超临界环境下煤油替代物液滴燃烧特性数值研究

为研究液氧/煤油火箭发动机燃烧室内经喷注形成的煤油液滴的燃烧过程,基于实际气体状态方程、高压热物性修正、高压气液平衡和详细化学反应动力学,建立一维的全瞬态液滴燃烧模型,对超临界环境下两组分煤油替代物液滴的燃烧特性及液滴初始直径的影响进行仿真研究。结果表明,在超临界环境下,相比于煤油液滴纯蒸发过程,煤油液滴燃烧过程的迁移时刻大大提前;煤油液滴着火之后很快进入超临界燃烧阶段,此时液滴燃烧过程可以看成中心附近的燃料高浓度区与外侧氧气高浓度区之间的扩散燃烧过程;煤油液滴的火焰半径先增大,达到最大值之后开始减小,并减小为零,火焰温度在着火之后快速上升至最大值,并基本保持不变,在火焰半径减小为零之后开始降低;随着液滴初始直径的增大,火焰特性以及液滴中心参数变化曲线趋势不变、整体延迟,着火时间、迁移时间和液滴寿命增大。     

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机,针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响,对全机带动力风洞试验模型进行了大规模并行非定常数值计算,再现了风洞试验现象,通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂,这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速,为保证飞行安全,继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示,在飞行雷诺数和螺旋桨转速下,相同侧风范围内风洞试验显示的流动不稳定因素基本消失,偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析,建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系,进行了偏航稳定性的虚拟试飞评估,研究成果也为AG600飞机的首飞和飞行试验所验证。     

涡桨飞机螺旋桨滑流气动干扰效应及流动机理

螺旋桨滑流对飞机各气动部件的干扰是涡桨飞机气动设计中面临的难点之一。以某双发涡桨支线客机为对象,采用数值模拟手段分析了滑流对机翼、平尾、垂尾的影响及其流动机理。计算采用基于动态面搭接网格技术的非定常方法,气动力与表面压力分布的计算结果均与实验值吻合良好。研究结果表明：在滑流影响下,全机升力和阻力有所增加,升阻比和纵向静稳定度有所降低,并在无侧滑状态下产生了滚转力矩和偏航力矩;在不同展向位置,滑流对机翼表面流动分离的影响存在显著的差异。在螺旋桨下行运动一侧,滑流的旋转减小了当地迎角,同时桨后气流速度较高,翼面流动分离被有效抑制。在螺旋桨上行运动一侧,滑流的旋转增大了当地迎角,而且桨后气流速度较低,因而翼面流动分离并未得到明显改善;在中小迎角下,滑流对背景飞机平尾当地的动压没有产生影响,但增加了下洗角变化率,进而导致平尾效率降低;垂尾的侧力与偏航力矩是由滑流的侧洗引起的,而滑流的侧洗又与左右两侧机翼升力分布的不对称性有关。     

空间电动力绳系电流峰值点变化规律研究

为研究电流峰值点在几种影响因素下的变化规律,建立一种能够捕捉绳系电子发射与收集自洽平衡过程的新算法——电路空间耦合算法。为验证该方法的计算精度,以1.35kW霍尔推力器为等离子体源,在真空舱内开展绳系的电荷收集试验,在电路参数方面,计算精度约为12.8%,在场参数方面,计算精度约为3.6%。在此基础上,针对不同偏置电压、绳系长度以及绳系直径,对绳系的绳上电流分布、电势分布以及空间电势分布等参数进行数值计算。结果显示:电动力绳系的电流峰值点会随着偏置电压升高、绳系长度增加及绳系直径增大而发生比例上的向阳极端漂移,揭示了电子轨道运动限制的机制在各类壁面电荷输运机制中占优,导致电流峰值点漂移的产生。