惯性静电约束推力器放电过程数值模拟研究

为了研究惯性静电约束推力器(Inertial Electrostatic Confinement Thruster,IECT)的放电原理和工作机制,采用漂移-扩散流体模拟方法,基于圆柱形惯性静电约束推力器的结构,研究不同栅网线直径、栅网个数、推力器尺寸条件下等离子体放电情况和阴极电压、背景气压对推力器放电的影响。结果表明:在所研究条件下,保证阴极必要的几何透过率的同时,适当增加推力器栅网个数可以提高喷射离子密度,减小羽流发散角;随着阴极电压和背景气压的增大,推力器喷射的等离子体密度增大。但是,压强继续增大会达到临界值,等离子体被约束在推力器内部无法喷出,即喷射模式无法运行,故阴极电压与背景气压对IECT均有较大影响。     

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

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

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

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

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

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

基于混合有限元解析法的永磁同步电机永磁体电涡流损耗估计

永磁同步电机永磁体受限于热约束,无法在温度较高的环境下运行,故需减少永磁体上的电涡流损耗,从而降低永磁体上的温度。针对使用有限元法对永磁体电涡流损耗估算时间较长,以及使用解析法估算时难以达到与有限元法相同的精度,采用混合有限元解析法估算永磁体上的电涡流损耗。结合电涡流的反作用,在模拟电机旋转时,无需重复划分三角形区域;使用MATLAB软件仿真模拟,将混合有限元解析法与Galerkin有限元法对比,减少三角形区域划分的个数。由此验证了永磁体上电涡流损耗符合端部效应以及集肤效应的特征,在保证精度的同时,减少了仿真的时间。     

基于同步电机电抗参数定义式的仿真计算方法

从同步电机各种稳态和瞬态电抗参数的定义、各种电抗参数与电机内部各绕组和磁路结构之间的电磁关系出发,提出了一种基于电抗参数定义式的仿真计算方法,给出了详细的仿真计算步骤。通过算例与样机实测结果的对比,表明仿真计算结果与实测值非常吻合,从而验证了所提计算方法的有效性和准确性。     

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

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

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构,层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟,研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型,基于结构化拼接网格,对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性,雷诺数对翼型气动特性具有显著影响。随雷诺数降低,时均分离泡外形增大、位置后移,平均阻力系数增大,特别是在较低雷诺数下,翼型升阻力系数随时间出现振荡现象。进一步研究表明,造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低,流动黏性增大,导致分离剪切层速度梯度减小,流动发生转捩及再附位置后移,直至翼型表面不再发生转捩和再附。     

一种内并联型内转进气道通道间干扰特性研究

为了研究工作马赫数范围Ma0～6的组合循环发动机进气道的工作特性,设计了一种三通道内并联型可调内转进气道,为了适应进气道能宽速域正常工作,在内转进气道的顶板压缩面上进行了变几何设计,可以调节进气道的内收缩比,兼顾了进气道低马赫数下的起动性能和高马赫数下气动性能。采用三维数值仿真的方法对进气道在过渡模态下的反压特性进行了分析研究。研究结果表明,可以通过改变隔离段反压或者涡轮通道反压的大小来调节两通道之间的流量分配和涡轮通道出口的稳态畸变;当增加涡轮通道反压时,涡轮通道的流量系数和稳态畸变会逐渐减小,而冲压通道的流量系数会逐渐增大,最大增加量约为原有流量的23.3%;当增加隔离段反压时,涡轮通道的流量系数和稳态畸变会逐渐增大,而冲压通道的流量系数会逐渐减小,最大减小量约为原有流量的29.3%。     

缝槽平板振动气膜冷却数值仿真研究

为了分析缝槽平板振动气膜冷却传热规律,采用基于Reynolds Averaged Navier-Stokes(RANS)和Detached-Eddy Simulation(DES)的CFD动网格仿真分析方法,研究了不同振动频率、不同振幅、不同吹风比及剪切涡结构振动的绝热壁面有效温比。结果表明:(1)在低频(f=100Hz)下,振动有利于提高换热面绝热壁面有效温比;但在高频(f=10kHz)下,振动降低了换热面绝热壁面有效温比。(2)相比频率,振幅对绝热壁面有效温比影响更大;在中频(f=1kHz)下,振动绝热壁面有效温比小于稳态绝热壁面有效温比;不同吹风比(M=0.7,1.2,1.9)对绝热壁面有效温比影响趋势基本一致。(3)相比基于RANS的振动仿真方法,基于DES的振动仿真方法能更详细地模拟混合层剪切涡结构流场信息,能更好地反映剪切涡结构对壁面换热效果的影响。