旋转爆轰燃烧室高度对流场影响

为了研究旋转爆轰燃烧室几何结构对旋转爆轰波的影响规律,针对旋转爆轰燃烧室高度对爆轰流场结构和爆轰波传播特性问题,采用2维数值模拟方法分析了爆轰波在4种高度的燃烧室内的参数变化和波系特征。结果表明:在4种高度的燃烧室内爆轰波均能自持传播,且爆轰流场结构相似;当燃烧室高度为4 mm时,爆轰波后的反射激波压力较高;在一定范围内,爆轰波传播速度随燃烧室高度的增加而增大,在爆轰波高度为10、16和20 mm时,爆轰波波锋面的压力分别为3.93、4.61、4.69 MPa。     

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

连续旋转爆轰波在无内柱圆筒内的数值模拟

在对同轴圆管内的连续旋转爆轰的已有研究成果进行总结分析后,提出无内柱的连续旋转爆轰的燃烧室模型.利用结合了一步化学反应模型的圆柱坐标系下的Euler方程三维数值模拟了在这种燃烧室内爆轰波的传播行为.数值结果表明:爆轰波由初始起爆的单向传播最终自动收敛稳定到若干个典型爆轰波头的流场结构.在当前进气条件下,流量可以达到450kg/m2s,其基于反应物的比冲约为1900s.与同以H₂/Air作为燃料时同轴连续旋转爆轰发动机(RDE)燃烧室模型的2000s相比,基于反应物的比冲十分接近,但由于内部有可燃气以爆燃形式燃烧等原因,有5%的损失.      

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

燃烧室轴向和周向长度对气液两相旋转爆轰特性的影响

为了研究旋转爆轰发动机燃烧室轴向和周向长度对气液两相旋转爆轰特性的影响,采用守恒元和求解元(CE/SE)方法对带化学反应的汽油/富氧空气两相旋转爆轰理论模型进行求解,获得了气液两相旋转爆轰流场结构,并分析了燃烧室轴向和周向长度对燃烧室流场、爆轰波传播特性以及发动机推力性能的影响。计算结果表明:轴向长度对燃烧室上游流场影响甚小,却对下游流场参数影响较显著。随着轴向长度增加,燃烧室出口压力、温度、密度以及周向速度均降低,轴向速度则逐渐增大,同时发动机平均推力密度和燃料比冲先增大后减小。当周向长度过短,燃烧室内难以形成自持传播的爆轰波,随着周向长度增加,上游爆轰波强度增加,对应的流场参数均有所增大,但发动机推力性能略有所降低。     

连续旋转爆轰燃烧室增压特性的数值研究

为了分析爆轰燃烧室几何尺寸对增压特性的影响,利用二维欧拉方程数值研究了当量H2/Air在连续旋转爆轰燃烧室中的燃烧流场情况。研究表明:在准稳定状态下爆轰燃烧流场特征参数随时间呈周期性振荡,且其振幅保持不变;横波及局部爆炸波的存在是爆轰波后侧附近总压沿x轴呈现波动变化的主要原因;影响连续旋转爆轰燃烧室增压特性的直接因素为斜激波的高度及爆轰波强度,减少轴向尺寸或增加周向尺寸可增强其增压特性。连续旋转爆轰燃烧室的增压比高达2.52,与传统的燃烧室相比,大大提高了燃气的做功能力。     

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段；在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波；第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波；稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s；斜激波的存在使燃烧室出口平面流场产生了较大波动。      

圆柱坐标系下连续旋转爆轰发动机的数值模拟

使用两步化学反应模型对连续旋转爆轰发动机(Continuous rotating detonation engine,CRDE)进行了二维数值模拟研究。数值计算获得了同轴圆管腔中间层曲面上连续爆轰的多个循环过程。分析了燃料入射、提前燃烧、爆轰波结构和波传播速度等几个关键问题。计算获得燃烧室内流场结构与之前实验研究结果定性符合。计算中燃料以最大1400m/s沿轴向入射,爆轰波可约以2400m/s沿周向连续循环传播。计算的燃烧室内爆轰波循环频率约为2300Hz。     

燃烧室宽度对液态燃料旋转爆轰发动机影响实验研究

为了研究燃烧室宽度对液态燃料旋转爆轰发动机工作特性的影响,搭建了气液两相旋转爆轰实验系统,以汽油/富氧空气为工质,氢气/氧气预爆轰管作为点火装置,在不同燃烧室宽度下开展了一系列实验研究,分析了爆轰波的起爆过程,以及燃烧室宽度对爆轰波传播特性与发动机推力性能的影响。实验结果表明：点火后,燃烧室内需要经过一个爆燃转爆轰过程才能形成自持传播的爆轰波;爆轰波在不同燃烧室宽度下均以双波对撞模态传播,对应的波速分布在850~1025m/s内,随着当量比增加,波速整体呈增加趋势;当燃烧室宽度减小,波速整体有所降低;不同燃烧室宽度下推力性能存在显著差异,其中燃烧室宽度在16.5mm下,发动机的推力和燃料比冲要明显低于11.5mm和9mm的;随着燃烧室宽度减小,内外壁面边界层在流场中的作用更为突出,降低了发动机推力的稳定性。     

旋转爆轰波与涡轮平面叶栅相互作用数值模拟

由于旋转爆轰燃烧室具有自增压特性，可提高热力循环效率，因此将旋转爆轰燃烧室应用于燃气轮机可进一步提高系统的性能。基于非稳态雷诺时均Navier-Stokes方法，采用剪切应力输运k-ω湍流模型，建立旋转爆轰燃烧室与涡轮平面叶栅耦合计算模型，研究旋转爆轰燃烧室内的复杂波系与涡轮叶片的相互作用，分析涡轮叶栅对高频爆轰压力振荡的抑制作用。结果表明：旋转爆轰燃烧室内的燃气在涡轮叶栅内加速，并且在斜激波后的局部区域马赫数的增加更为明显。斜激波与涡轮静转子叶片的前缘、压力面、吸力面以及尾缘相互作用，由于旋转爆轰波不同的传播方向，使得斜激波与静子叶片呈相互垂直或平行，进而形成两种不同的波系结构。涡轮叶栅对高频压力振荡存在明显的抑制作用，涡轮叶栅上下游高频压力振荡幅值的衰减率达到80%以上。研究结果展示了旋转爆轰波作用下涡轮叶栅内复杂波系结构特征，并对基于爆轰燃烧推进技术的应用提供了一定的理论基础。     

聚心火焰与诱导激波相互作用及爆燃转爆轰过程

基于带化学反应的二维轴对称Eu ler方程,利用带限制函数的波传播算法,在两端均敞开和一端封闭一端敞开的两种边界情况下,对圆柱型燃烧室中的氢气-空气预混气的聚心燃烧进行了数值模拟,探讨了聚心火焰引发爆轰过程中的波系结构以及火焰的变形。数值结果表明,在两端均敞开的情况下,燃烧诱导的激波在轴心、火焰和壁面之间来回反射,使火焰失稳,加速了火焰的燃烧,但没有转变成爆轰。而在一端封闭的情况下,封闭端反射产生的激波,不断穿越火焰,使火焰严重失稳,加剧了燃烧速度,最终导致爆轰的形成。稳定爆轰阵面的参数与CJ理论计算的结果进行了对比,两者符合的较好。同时,火焰在与激波的作用过程中,形状扭曲变形,对于两端敞开的情况呈对称的扁平头部蘑菇云,而对于一端封闭的情况,则呈封闭端小敞口端大的扁平头部蘑菇云。     

高总温空气与汽油燃料的旋转爆震验证试验

为验证高总温空气来流条件下汽油燃料旋转爆震的可行性,开展了气液两相旋转爆震发动机试验研究。旋转爆震发动机环形燃烧室外径和内径分别为202mm和166mm,长度为155mm。通过空气加热器模拟高总温空气来流环境,汽油和空气采用分开喷注的方式,分别通过高压喷嘴和环缝进入燃烧室。试验采用垂直安装的预爆震管成功起爆了旋转爆震波,并实现了旋转爆震波的连续稳定传播。试验结果表明：当空气质量流量为1110.0g/s,当量比为0.97,空气总温为713K时,旋转爆震波以双波对撞模态在燃烧室内连续传播,爆震波传播频率为1827.31Hz,与高频压力信号经快速傅里叶变换得到的主频一致,爆震波传播速度为1059.6m/s。在空气质量流量为1110.0g/s,当量比为0.84,空气总温为713K的工况下进行了3s的长程试验,验证了以高总温空气为氧化剂、汽油为燃料的旋转爆震发动机长时间连续稳定工作的可行性,获得的旋转爆震波传播频率为1907.5Hz。     

圆盘旋转爆震燃烧室中爆震波传播模态

本文的目的在于通过求解14组分19步CH4/O2反应的欧拉方程去研究燃料喷射温度、压力、燃烧室内外直径比值和燃烧室径向长度对圆盘结构下旋转爆震传播模态的影响。研究表明,只有预混气喷射温度范围为500-900K,喷射压力范围为0.5-3.5MPa的条件下能够在固定燃烧室内形成连续的旋转爆震波。受燃烧室结构影响,旋转爆震波的传播模态分为稳定模态和非稳定模态。整个非稳定传播模态根据爆震波对前周期的干涉又可以分为干涉阶段和不干涉阶段。在不干涉阶段,爆震波传播速度略高于Chapman-Jouguet (CJ)速度;在干涉阶段,爆震波传播速度低于CJ速度且爆震波被间断面分割成两部分。非稳定模态爆震波传播速度小于稳定模态,而爆震波夹角、燃烧室出口面积比则反之。     

Investigation of hot jet on active control of oblique detonation waves

The hot jet injection is utilized to actively control the oblique detonation wave, such as initiating and stabilizing an oblique detonation wave at a desired position that is shorter than the length of induction zone, and adjust the height of the oblique detonation wave at the exit of combustor when the oblique detonation wave engine is working on off-design flight conditions. The fifth order Weighted Essentially Non-Oscillatory (WENO) scheme and a two-step reversible reaction mechanism of the stoichiometric H₂/Air are adopted in the simulations. With the help of hot jet injection, the transition from inert oblique shock wave to the oblique detonation wave immediately occurs near the position of hot jet injection, and consequently the length of combustor can be reduced. The angle of oblique detonation wave also decreases as the hot jet injection approaches the nose of wedge. Additionally, the height of the oblique detonation wave at the exit of combustor can be flexibly adjusted, and also depends on the injection position and the strength of the hot jet. If the velocity of the hot jet is too weak to directly trigger the overall oblique detonation wave at the position of injection, increasing the injection pressure will improve the strength of the hot jet and results in a successful transition.     

脉冲爆震燃烧室与涡轮相互作用的试验

采用汽油和空气作为燃料和氧化剂,进行了脉冲爆震燃烧室(PDC)与涡轮相互作用的原理性模型试验装置的试验.试验结果表明:PDC工作平稳,在发动机爆震工作时压气机出口空气质量流量比用于PDC产生爆震的空气质量流量大100kg/h左右;在PDC出口前已形成充分发展的爆震波,压力波经过涡轮膨胀后峰值压力和波速明显降低;且充填系数越大,各位置处的平均峰值压力越高,压力波经过涡轮后的衰减越小;涡轮在经受累计40多分钟共12 000多次脉冲爆震波或压力波的冲击后仍然完好无损.     

多循环脉冲爆震发动机爆震室内部压力分布实验研究

利用基于多循环吸气式脉冲爆震发动机的实验获得了沿爆震室轴向9个位置的实时压力数据,经过对数据进行时域和频域的分析处理后,利用MATLAB软件的GRIDDATA函数拟合出了6个工作频率下某个循环时爆震室内压力随轴向位置及时间变化的三维图形。通过对拟合后的数据从不同角度进行分析,发现采用爆燃向爆震转捩（DDT）方式起爆会在爆震波形成前生成一个局部爆轰区、爆震波形成后平均速度维持在1000m/s左右,且峰值压力升高、随机性增大。该实验的研究结果为爆震室的结构强度研究提供了依据。     

预爆管爆震波衍射特性数值研究

为了实现预爆管出口爆震波向主爆管的成功传播,针对预爆管出口不同几何形状,开展了数值模拟研究,获得了爆震波从预爆管向主爆管传播的衍射特性.研究表明:当预爆管出口为直喷管时,不能实现向预爆管出口截面上游传播的爆震燃烧;当预爆管出口扩张角为45°时,可以实现爆震波在主爆管内向衍射面上、下游传播,研究结果获得的爆震波衍射特性得到了试验验证.