模型燃烧室内不稳定燃烧发展过程的数值分析

为了对自燃推进剂燃烧室内出现的不稳定燃烧现象进行详细分析,采用欧拉-拉格朗日方法对该燃烧室内的气液两相非稳态燃烧过程进行了数值模拟,计算得到的压力振荡幅值和频谱特性与实验结果吻合较好,在此基础上对不同燃烧阶段内的压力和释热变化规律进行了分析。结果表明:在压力振荡幅值超过10%的不稳定燃烧阶段,压力振荡主频为9 200 Hz,燃烧室内横向压力分布与1阶切向振型一致,仿真中再现了1阶切向自激高频不稳定燃烧的产生及发展过程;稳定燃烧向不稳定燃烧转变早期,压力振荡从部分燃烧释热波动中获得能量,压力振荡幅值缓慢增长;随着燃烧进行,燃烧释热波动与压力振荡之间相位和频谱特性逐渐趋于一致,压力振荡幅值开始急剧增大;当二者完全耦合时,燃烧室内压力振荡幅值达到极限饱和状态,此时压力振荡幅值超过了平均室压的200%。      

转静干涉对转子叶片颤振特性的影响

在不同气动工况、不同几何模型下,采用自行开发程序对全环多排的高压压气机进行了流固耦合数值模拟,分析了进口导流叶片(IGV,以下简称导叶)对转子叶片颤振稳定性的影响。通过对气动弹性标准算例4的数值模拟验证了程序在颤振领域的有效性。针对导叶-转子模型和单转子模型,考察了近堵塞点、近设计点和近失速点3个工况下,节径变化对叶片颤振稳定性的影响,给出了气动弹性最不稳定状态对应的叶片振动形式。通过对比发现,导叶作用随工况而异,近堵塞点导叶使得转子1阶弯曲模态气动阻尼提高130.63%。研究表明:导叶引起的非定常压力波反射增强了转子叶片的非定常压力扰动幅值,使得弯曲振型的颤振稳定性增强。基于单转子模型的颤振分析给出了不准确的气动阻尼值。      

航空发动机一维数值仿真

分析了航空发动机一维气动热力数值仿真模型,此模型能够用于航空发动机整机数值仿真.进行了整机性能数值仿真,求解过程考虑了粘性效应以及气体流动损失对于整机性能的影响;采用了具有TVD性质的高阶精度Godunov格式,此格式具有良好的数值稳定性和激波捕捉能力;应用了时间推进法,交替使用显/隐格式.计算了某型带加力的混排涡扇发动机的速度特性以及节流特性并且与设计数据进行了比较分析,数值仿真结果与设计数据比较吻合.分析了典型工作点的计算过程,模型收敛速度较快.     

飞翼布局无人机进排气影响及机理分析

进气道尾喷管的存在不仅影响着飞机机身前后的流场特性与压力分布,而且影响着飞机升力特性、阻力特性以及力矩特性.为了研究进排气效应对飞机气动特性的影响,采用计算流体力学(CFD)数值模拟技术,建立飞翼布局无人机(UAV)堵锥整流模型和动力影响模型,其中动力影响模型是根据发动机不同的工作状态,在进气道以及尾喷管截面上设置不同的边界条件,对比分析这两种计算模型在流场特征和气动特性上的差异性,揭示内在的影响机理.进排气效应的研究有助于飞机/发动机的一体化气动综合设计.      

三维多组分羽流及其效应的数值模拟

考虑到实际卫星总体布局特点及姿控发动机燃气喷流为多组分的特点,采用差分求解N-S(Navier-Stokes)方程耦合DSMC(Direct Simulation Monte Carlo)方法建立三维多组分羽流场的数值模拟模型,对5 N无水肼单组元姿控发动机喷管产生的燃气羽流场及其效应进行数值模拟,计算得出羽流场的物理参数分布、化学组分分布,得到了羽流对太阳能帆板的羽流气动力和气动热效应.计算结果表明对于太阳能帆板的羽流效应防护应以N2,NH3为主,姿控发动机羽流会对太阳能帆板撞击产生气动力和气动热的影响,在太阳能帆板前沿中点处,存在气动压强峰值和气动热流峰值.计算结果为研究羽流效应提供了参考依据.     

面向气热耦合的涡轮叶片计算域模型建模方法

针对冷却结构复杂的涡轮叶片在气热耦合数值模拟中计算域模型建模效率低、模型质量不稳定以及对数值模拟适应性差等问题,分析面向制造的涡轮叶片模型及建模方法,结合气热耦合数值模拟对涡轮叶片数值模拟的实际需求,提出一种面向气热耦合的涡轮叶片计算域建模方法。通过外型冷却特征自动定位生成算法,创建涡轮叶片冷气域部分;通过自适应管道求交算法、边界自动匹配算法,创建能够适应不同叶型的涡轮叶片燃气域部分;此外,在创建流体域的过程中提取气热耦合数值模拟所需的关键几何特征及非几何信息,并与冷气域、燃气域、叶片实体集成,完成气热耦合计算域模型的生成。基于以上研究开发涡轮叶片气热耦合计算域模型快速建模系统,验证所提方法的有效性。      

非线性结构气动伺服弹性稳定性分析

真实的气动伺服弹性系统通常都含有各种非线性因素,它们会对系统的特性带来不可忽略的影响.为此采用一种频域方法对带有间隙的非线性结构进行气动伺服弹性稳定性分析.在气动伺服弹性运动方程中,非线性的结构刚度通过描述函数法转化为准线性的等效刚度,然后采用线性控制理论中的Nyquist方法来判断气动伺服弹性系统的稳定性并计算稳定裕度.以某一飞机的带有扑动和旋转间隙的受控翼面为例进行了数值计算.结果表明,旋转间隙对系统的稳定性影响是主要的;旋转间隙越大,系统稳定裕度越小.      

含间隙非线性机翼跨声速颤振时滞反馈控制

颤振主动控制会引入时滞,对气动弹性系统闭环稳定性具有显著影响。针对当前考虑时滞的机翼颤振主动控制多集中在亚、超声速域,采用线性气动力分析的研究现状,结合现代飞机大都以跨声速巡航、控制面偏转为作动器进行主动控制的应用特点,发展了考虑结构间隙非线性,基于气动力降阶模型的跨声速颤振时滞反馈主动控制方法。首先,以白噪声为激励信号,辨识得到跨声速下非定常气动力降阶模型,与间隙非线性结构模型耦合,构建被控对象状态空间模型;然后,通过一种含积分项的状态变换将输入信号存在时滞的被控系统转化为无时滞的系统;最后,采用最优控制理论设计最优时滞反馈控制。仿真结果表明:对于含时滞的系统,若施加不考虑时滞影响的控制方法,则无法抑制颤振,所提控制方法的有效性不受时滞大小的影响,可有效抑制颤振的发生。      

一种单喷嘴推力室燃烧内流场的方法

为了方便开展单喷嘴气气喷注器燃烧试验研究,设计了一种燃烧室壁面测温,并结合推力室燃烧流场的非稳态数值模拟的方法,以考察燃烧室内流场的发展和燃烧完成情况.应用试验数据对不同的湍流模型和燃烧模型进行了对比研究,得到与试验结果较为吻合的计算模型.并应用该方法开展了同轴剪切式喷注器氢/氧喷注动量比、燃烧室压力变化的试验研究.试验结果表明喷注器掺混燃烧效果随氢/氧喷注动量比增加而增强,而不随燃烧室压力的变化而变化.     

甲基肼/四氧化二氮发动机脉冲工况仿真与试验研究

基于简化化学反应模型,对甲基肼/四氧化二氮发动机脉冲工况燃烧特性进行了仿真分析,并结合高模热试车对仿真结果进行了对比验证.结果表明:燃烧室压力在3.7 ms时基本达到稳定;而温度场在8 ms左右达到稳定,室压先于温度场达到稳定;燃烧室壁面液膜热流密度先增大后降低,脉冲工况下发动机推力室壁面液膜冷却效果受脉冲宽度影响较大...     

探测器触地关机软着陆稳定性分析

深空探测器软着陆过程中发动机未能正常关闭将极大降低其软着陆稳定性,因此引入了触地关机着陆方案。针对触地关机方案,建立了腿式探测器动力学分析模型及发动机推力控制模型。考虑着陆环境与探测器着陆状态的不确定性,采用蒙特卡罗模拟分别对探测器在主发动机关闭失败时、触地关机着陆方案下、带有姿态控制的触地关机着陆模式下的软着陆稳定性进行了分析。采用描述性采样方法抽取样本点,基于均值估计相对误差建立模拟终止准则,在保证模拟精度的前提下提高蒙特卡罗模拟效率。计算并对比了3种着陆方式下探测器软着陆可靠度,触地关机方案可在主发动机关闭失败后将探测器稳定着陆可靠度提升11.6%,姿态控制的引入可进一步将安全着陆的可靠度提升9.7%。      

液体火箭发动机试验脉动压力数据分析方法研究

针对液体火箭发动机试验脉动压力数据的特点,对脉动压力数据和相关振动测点数据进行快速傅里叶变换（FFT）,频谱分析的结果能够有效判断发动机不稳定燃烧的声振类型。采用小波变换奇异点检测的方法,对启动段和关机段波动异常的脉动压力数据进行分析,结合相关振动测点的频谱信息能够在时域内准确定位奇异点位置,提取频域内的信息,有效判断推进剂燃烧情况和声振类型。分析结果表明上述方法在火箭发动机和组合件试验的脉动压力数据分析中具有重要的推广意义。     

关于太阳日冕磁环的稳定性

本文指出:在过去多数的日冕磁环稳定性的研究中,都存在一个未作说明、却又不容忽视的问题,即无形中都在磁环边缘处用了固壁边界条件,而它的采用与否对结果影响极大;另外在这些研究中只考虑m=1的模式也是不够的。这些问题都有待于进一步研究。      

跨音速风扇全环叶片颤振特性的流固耦合分析

发展了求解叶片颤振问题的流固耦合计算方法和全环叶片振动的气动弹性模型,在每一时间步同步求解流体运动方程和叶片振动方程并交换边界信息;流体域求解了非定常雷诺平均N-S方程,得到每一步由于叶片变形而引起的流场变化;叶片变形则由积分叶片表面受到的气动力并求解结构动力学方程得到.颤振分析是在全环叶片模型上进行的,并解除了预先设定叶片间相位角的限制.此方法的显著特征是在一次气动弹性计算过程中,可同时分析叶片多个固有模态、多个节径下的气动弹性稳定性,大大提高了使用时域法进行叶片排气弹分析的计算效率.考察了NASA rotor 67风扇全环模型在堵塞点、最高效率点和近喘点3个气动工况下,节径变化对叶片气动弹性稳定性的影响,给出了不同模态下气弹最不稳定状态对应的叶片振动节径形式.结果表明,振动形式对于叶片气动弹性稳定性的影响很大.     

涡扇发动机地面起动过程性能模拟

用积分方法建立了混合排气加力涡扇发动机起动过程模拟模型并用该模型进行了仿真计算,模型基于气动热力学方程、各部件特性以及经验关系式,根据给出的发动机起动过程初始条件、燃烧室供油规律等进行计算.为说明所建立的起动性能模型的模拟效果,将程序计算结果与试验结果进行了比较,结果表明该起动仿真模型能较准确地模拟混排加力涡扇发动机的地面起动过程.还分析了计算过程中对起动过程有比较大影响的因素,包括计算时间步长,起动机脱开转速,燃烧室总压恢复系数的经验系数,大气温度、压力以及海拔高度等.了解这些因素对起动过程及模拟结果的影响,便于通过修正的途径将该计算方法用于发动机地面起动性能的准确模拟.      

开放式微波结构中的时域模函数初探

为了探讨将Diakoptics算法应用到开放式微波结构时适宜的时域模函数形式,在对槽线、共面波导、微带线及耦合微带线的空间电磁场分布特性分析的基础上,提出零阶、一阶第一类贝塞尔函数是较适宜的时域模函数,并总结了选取规律.      

FDTD-Diakoptics时域模函数分析

时域模函数的选取是使用Diakoptics技术分析微波结构的一个重要步骤.二维贝塞尔函数系作为Diakoptics技术应用于开放式微波结构中的时域模函数是适宜的.将一维傅里叶-贝塞尔函数展开定理扩展到二维,对二维傅里叶-贝塞尔函数系的完备正交性进行了分析证明,从而阐述了贝塞尔函数系作为时域模函数的有效性.最后比较采用二维傅立叶－贝塞尔函数系做模函数的FDTD-Diakoptics计算结果与传统FDTD计算结果,二者吻合良好.