同轴突扩燃烧室低频不稳定燃烧数值模拟

采用大涡模拟耦合预混燃烧方法,建立了液体冲压发动机不稳定燃烧计算模型。对模型发动机进行了不稳定燃烧数值研究,计算得到的发动机内压强振荡和火焰传播过程与文献中的实验数据吻合较好。结果分析表明,突扩构形发动机内大尺度的旋涡结构支配了燃烧室中的火焰传播过程,旋涡运动耦合非稳态的燃烧热释放是激发燃烧室低频压强振荡的重要原因。     

侧面突扩燃烧室冷态流场可视化研究

燃烧室流场中旋涡的不是稳定是造成整体式压式发动机侧面扩燃烧室燃烧振荡的重要原因。在透明矩形侧面突扩燃烧室模型上进行了水流模拟显示实验。显示出了燃烧室流场中的振荡涡系和稳定涡系。实验发现：燃烧室头部旋涡非常稳定；进口射流剪切层存在“马蹄涡”的周期性脱落；射流在燃烧室通道内卷绕形成涡强较大的二次流螺旋柱状涡对，涡对相撞又使螺旋涡失稳、振荡和破碎。旋涡不稳定性是侧面突扩燃烧室燃烧振荡的流体力学原因，分流     

加力燃烧室流场形态与振荡燃烧数值模拟

利用数值方法对发动机加力燃烧室内的流场形态和振荡燃烧进行模拟,获得流场平均参数及涡系结构,观察到燃烧室内火焰稳定器唇口附近剪切涡对选择模型的敏感性,发现火焰稳定器后的涡脱落和回流区的周期性变化是引起加力"蜂鸣"现象和振荡燃烧的重要因素,证明了Rayleigh准则在研究此类问题时的适用性.      

细长旋成体大迎角俯仰振荡的数值研究

为研究俯仰振荡对细长旋成体流场与气动特性的影响,利用雷诺平均方法与大涡模拟的混合方法相结合的脱体涡模拟方法对大迎角下细长旋成体俯仰运动进行了数值模拟研究。通过对背风面分离涡的强度和位置、细长体绕流形态、截面压强系数和侧向力系数的观察和分析,发现施加小振幅、特定频率的俯仰振荡对大迎角细长旋成体背风面流场有明显的控制能力,...     

细长旋成体大迎角非定常运动数值研究

利用DES方法对处于非对称流态下的细长旋成体进行了俯仰和偏航振荡的数值模拟,观察细长体气动特性尤其是背风面非对称分离涡的变化。计算结果表明在大迎角非对称多涡系情况下,固定频率和振幅的非定常运动可以改变流场结构和气动力,对细长体背风面流态有巨大的影响。俯仰振荡对非对称分离涡有明显的控制作用,抑制流场的非对称性,使分离涡趋于对称;而固定频率的偏航振荡则破坏背风面分离涡的稳定性,激励非对称背涡产生随时间变化的周期性脱落。所进行的非定常运动与细长体流场耦合作用研究在国内外研究尚少。     

模型固体火箭发动机内涡脱落过程的大涡模拟

为研究后向台阶结构引起的涡脱落对固体火箭发动机内流动稳定性的影响，采用大涡模拟技术，对模型发动机内的流动进行了数值模拟，获取了发动机燃烧室内压强随时间的变化曲线及其振荡频率，并和已有的计算结果及实验结果进行了对比分析。结果表明采用LES方法的计算结果明显优于采用k-ε湍流模型的计算结果；不考虑燃烧时模型发动机内由涡脱落引起的压强振荡的振幅较小，但频率较高，且随来流速度的增加而增大；在后向台阶这种结构中，涡一般会在台阶的边沿形成，并在其下游周期性地脱落；在保证y^＋较小的情况下，网格大小对计算结果的影响不明显。     

双燃烧室冲压发动机超声速混合层混合增强数值研究

研究了具有双燃烧室冲压发动机超声速燃烧室工程背景的特定超声速混合层的混合增强措施——低速流入口加入沿流向振荡.研究结果表明:它对二维混合层的混合增强效果明显,加振荡后流场燃烧强度加大,但仍属于较弱燃烧强度;振荡频率决定着涡量聚集区域的位置和大小,振幅决定着涡量分配和聚集的均匀程度.      

受限空间二维超声速混合层混合增强数值研究

本文研究了具有双燃烧室冲压发动机超声速燃烧室工程背景的受限空间超声速混合层的混合增强措施——低速流入口加入沿流向振荡。研究结果表明：它对二维混合层的混合增强效果明显,加振荡后流场燃烧强度加大,但仍属于较弱燃烧强度;振荡频率决定着涡量聚集区域的位置和大小,振幅决定着涡量分配和聚集的均匀程度。     

支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。     

预混燃烧大涡模拟和燃烧模型的检验

用代数二阶矩亚网格(ASOM-SGS)燃烧模型对文献中测量的钝体后方丙烷-空气预混燃烧进行了大涡模拟,模拟统计的时平均速度、速度脉动均方根值和温度分布与实验数据结果吻合很好,表明所采用的ASOM-SGS亚网格燃烧模型是合理的。模拟的瞬态结果显现了钝体后方湍流流动和火焰结构。将大涡模拟数据统计得到的反应率系数-浓度关联量的分布规律,与代数二阶矩RANS(ASOM-RANS)燃烧模型的模拟值进行对比,结果发现,大涡模拟统计值和ASOM-RANS模型的模拟值很接近,从而证明了湍流燃烧代数二阶矩RANS模型的合理性。     

发动机谐振与箭体模态耦合稳定性及机理研究

针对发动机谐振可能与箭体弹性振动之间发生耦合共振,从而引起姿态控制系统不稳定的问题,首先将发动机谐振特性简化为一个二阶环节,建立了包含发动机谐振方程的火箭姿态动力学模型;然后分析了发动机谐振对箭体模态极点频率和阻尼比的影响,并从理论上分析了出现负阻尼比时发动机和箭体模态间的耦合共振机理,在此基础上计算了导致姿态控制系统...     

基于全连接神经网络的分层旋流火焰燃烧振荡预报

燃气轮机贫油预混低排放燃烧室容易发生燃烧振荡,在宽工况极端条件下燃烧主动控制有潜力成为抑制燃烧振荡的有效方法。然而主动控制系统依赖于对燃烧振荡的快速监测和准确预报,因此有必要针对不同的燃烧振荡预报手段开展研究和验证。本文以甲烷预混同心分层旋流火焰的图像为基础,采用降低图像分辨率和提取火焰结构特征参数这两种不同的方式对火焰图像信息进行简化处理,并使用全连接神经网络对燃烧振荡进行预报研究。结果发现,两种方式都可以较为准确地预报燃烧振荡,精度均达到90%以上。预报精度随着图像分辨率的增加而升高,在极低的图像分辨率(3*3)下,预报精度也能达到90%以上。此外,对根据火焰平均图像提取的结构特征参数进行了敏感性分析,捕捉到了系统稳定性的转变,但参数变化范围受训练集限制。本文提出的基于数据驱动方法对燃烧振荡的预报时间小于2 ms,为实现燃烧振荡实时在线预报提供了支持。     

直升机空中振荡故障的试验研究与预防

 首先概述了故障的主要现象与特征。然后着重描述了直升机空中振荡故障的试验研究:包括研究途径;直升机模型的建立;两个回路振荡试验及其主要研究结果。据此提出了在设计,使用中的一些预防原则和措施。     

吸气式高超声速飞行器俯仰/滚转耦合运动特性

针对一种类似SR-72构型的吸气式高超声速飞机开展了进气道通流状态下俯仰/滚转耦合运动相关研究。通过数值模拟获得了滚转单自由度静稳定性、动稳定性以及强迫俯仰/自由滚转运动下的两自由度耦合动稳定性，研究了飞行器转动惯量以及俯仰运动频率对耦合运动的影响，简要分析了耦合运动的机理。研究发现虽然此飞行器具有滚转静稳定性和动稳定性，但是在强迫俯仰/自由滚转运动过程中，滚转通道却出现了小幅度振荡与大振幅振荡交替出现的情况，最大滚转角超过70°。小幅度振荡出现在正弦俯仰振荡的上半周期，其振荡频率随轴向转动惯量增加而降低，幅值随俯仰振荡频率增加而增大；大振幅振荡出现在下半周期，其幅值基本不变，而振荡频率与俯仰振荡一致。这种现象基本不受惯性耦合作用影响，可以认为是由气动力主导的。     

带抽吸二元进气道/隔离段激波串振荡特性

针对抽吸缝作用下激波串非定常振荡的复杂流动问题,采用高速纹影结合壁面动态压力测量的方法,在马赫数6的激波风洞中研究了高马赫数二元进气道/隔离段中激波串的自激振荡特性。隔离段出口不同堵塞度的实验结果表明:在低堵塞度下,隔离段内的分离激波无明显振荡;在中等堵塞度25.3%～32.3%和高堵塞度35.3%～38.2%工况下,隔离段内产生非定常激波串,受到隔离段内预先存在的背景波系以及抽吸缝泄流作用的影响,分别出现大幅度低频振荡和小幅度高频振荡;而当堵塞度超过临界值后,激波串被推出进气道,出现不起动。在大幅度振荡模式中,上壁面大分离区周期性地形成和消失,下壁面的激波串前沿分离激波在抽吸缝后缘到隔离段出口之间大幅度振荡,其振荡主频约为280Hz～480Hz,并且随着堵塞度升高而降低;在小幅度振荡模式中,上壁面始终存在大分离区,下壁面的激波串前沿分离激波在抽吸缝附近小幅振荡,其振荡主频约为900Hz～1800Hz。两种振荡模式均给隔离段壁面带来严酷的脉动压力载荷。     

Effects of pressure oscillation on aerodynamic characteristics in an aero-engine combustor

The effects of pressure oscillation on aerodynamic characteristics in an aero-engine combustor are investigated. A combustor test rig is designed to simulate the pressure drop characteristics of a practical annular combustor. The pressure drop characteristics are firstly measured under atmosphere condition with non-reacting flow (or cold flow), and the air mass flow proportion of each component (dome/liner) are obtained; these properties are base lines for comparison with combustion state. The combustion tests are then carried out under conditions of inlet temperature 340–450 K, fuel air ratio 0.010–0.028. The stability map and the oscillation frequencies are obtained in the tests, the results show that pressure oscillation amplitude increases with the increase of fuel air ratio. Phase trajectory reconstruction is applied to classify the pressure oscillation motion; there are three motions captured in the tests including: “disk”, “ring” and “cluster”. The pressure drops across the dome under strong pressure oscillation are distinctly divergent from the cold flow, and the changes of pressure drops are mainly affected by pressure oscillation amplitude, but is less influenced by pressure oscillation motion nor oscillation frequencies. Based on the mass flow conservation, the reduction of effective flow area of combustor under strong pressure oscillation is demonstrated. Liner wall temperatures are analyzed through Multiple Linear Regression (MLR) method to estimate the reduction of the air mass flow proportion of the liner cooling under strong pressure oscillation. Finally, the air mass flow proportions of each component under strong pressure oscillation are estimated, the results show that the pressure oscillation motion also has influence on air mass flow proportion.     

螺旋形强制谐波振荡对涡破碎位置的影响

本文用非线性波理论分析并计算了螺旋形振荡对涡破碎位置的影响,并在旋涡发生器中进行了施加不同频率的振荡对涡破碎位置影响的实验。实验与计算结果的趋势基本一致。     

浅谈微加速度计的整流误差

为明确微机械加速度计在振荡输入过程中直流输出的情况,着重分析了静电力反馈电容式微加速度计振荡频率低于微加速度计带宽和在微加速度计带宽附近两种情况下的失衡机理。为说明微加速度计与传统加速度计的整流误差差异,对两种整流误差的情况进行了半定量计算。本文的结论说明了微加速度计与石英挠性加速度计等传统惯性仪表在整流误差方面的区别;在此基础上,引申出由于整流误差决定不同种类的微加速度计有着不同标定要求和适用范围的结论。     

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。     

高速进气道低马赫数不起动特性及马赫数影响规律

为了探究进气道低马赫数不起动时的振荡特性,本文结合一体化前体/进气道构型,通过非定常仿真手段,对比研究了来流马赫数变化对进气道低马赫数不起动振荡流场以及飞行器气动力的影响规律。结果表明：低马赫数不起动时出现了稳定的振荡周期,且周期随着来流马赫数的增大而增长。由于拥塞发生在喉道处,其振荡流场单纯的表现为口部分离包的涨大和缩小,并且沿程压力的均值和幅值都呈现出喉道高两头低的分布趋势,而马赫数的增大会加剧此趋势。喘振周期中升力系数CL和阻力系数CD的变化趋势大致相反,升阻比曲线则表现为随分离包吐出而增大、吞入而缩小的趋势。CL和CD随着马赫数增大是整体下降的,但是脉动幅值变化不大,升阻比对马赫数的变化也并不敏感。此外,在进气道实现自起动过程中,当喉道瞬时流量高于起动时的流量一定程度,口部分离包将完全吞入。但定常仿真难以准确模拟该吞入过程,因此定常仿真得到的自起动马赫数偏高。