激波聚焦火焰聚心起爆器的数值研究

针对短距离内的激波聚焦起爆问题,设计了一种激波聚焦火焰聚心起爆结构,通过数值模拟方法对该结构内部流场变化进行研究。结果表明:加速区是火焰失稳加速的重要区域;聚焦区能量汇聚造成局部爆炸是燃烧转爆轰加速的关键因素,且该区域内由于几何结构形成的卷吸涡能够增强聚焦并阻碍激波回传入加速区;稳燃区主要实现爆轰波从过驱状态到稳定爆轰状态的过渡区。确定中心轴线上温度、压力波动是由于聚焦腔聚焦作用造成的,而爆轰波锋面温度压力变化是由于三波点运动造成。      

受限空间内激波与火焰作用的三维计算

激波与火焰的相互作用常发生在超声速燃烧和燃烧转爆轰过程中,为深入了解这一现象,采用带化学反应的三维Navier-Stokes方程,对平面入射激波及其反射激波与火焰的作用进行了计算研究,其中燃烧过程使用单步反应模型描述。研究结果显示:在受限空间内模拟激波与火焰作用,能更好地符合实验结果,从而体现出受限空间的三维效应;火焰在入射激波的作用下主要经Richtmyer-Meshkov不稳定而发生变形,此时火焰变形以物理作用为主,燃烧膨胀效果相对不明显;当反射激波与变形火焰再次作用后,火焰迅速膨胀变形,放热率维持在较高水平,此时化学反应过程起主要作用;在反射激波的作用下,变形火焰复杂三维涡结构的形成能强化热量与质量的输送,提高燃烧速率。     

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

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

斜爆轰波系在受限空间内的演变及其临界条件的数值研究

作为一种面向未来高马赫数吸气推进技术的新概念发动机,斜爆轰发动机的运行依赖其稳定的宏观波系结构。本文针对斜爆轰发动机简化模型,采用多组分基元化学反应流动数值模拟技术,数值分析斜爆轰波系在受限空间内的宏观结构特征及其演变,并进一步分析爆轰波系结构转变的临界条件。研究结果发现,随着楔面角度的增加,依次出现四种结构：激波诱导燃烧,斜爆轰波双规则反射,回流区马赫反射,楔面燃烧。对于稳定的波系结构,楔面压缩角同时存在下临界、亚临界以及超临界三种极限条件,在波系演变过程还伴随激波规则反射和马赫反射的转变。     

旋转爆轰流场的数值模拟

基于带化学反应的二维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对充有当量比的氢气-空气预混气和空气的环形旋转爆轰流场,从点火燃烧到发展成旋转爆轰的过程进行了数值模拟。根据数值结果分析了波后流场中爆轰产物受激波、高温和离心力等作用而挤向外壁,形成有利于充入燃料,实现持续稳定旋转爆轰的流场特征。还讨论流场中爆轰波、激波与间断面和内外壁面反射或折射,从而形成多个激波相交的波系特征。为认识和理解旋转爆轰流场,开展旋转爆轰的实验研究等具有现实的指导意义。     

中心锥喷管喉道比参数对旋转爆轰燃烧影响的数值研究

旋转爆轰发动机的推进性能与尾喷管的设计有关。为探究中心锥喷管的喉道比参数如何影响旋转爆轰波的传播模态及旋转爆轰发动机的推进性能，采用三维欧拉方程结合氢气/空气基元反应模型，对不同喉道比参数下的旋转爆轰发动机进行了数值模拟与流动分析。数值结果显示，随着喷管喉道比的减小，旋转爆轰波依次呈现出稳定单波模态、单双波交替的混乱燃烧模态和爆燃燃烧模态。研究发现，喉道处产生的反射激波是影响RDW传播模态的关键因素。随着喉道比的减小，反射激波强度逐渐增强，并与新鲜来流混气作用引发局部热点，使单波模态转变为单双波交替的混乱燃烧模态；随着喉道比进一步减小，反射激波的高压作用使爆轰波熄爆，燃烧模态转变为爆燃燃烧。对推力性能的进一步分析表明，单波模态下发动机的推进性能随喉道比的减小而增强，且明显好于混乱模态和爆燃模态。     

旋转爆轰的三维数值模拟

旋转爆轰利用爆轰波在燃烧室的供气端传播,而爆轰产物从另一端排出。为了揭示其宏观特性,采用带化学反应的三维Euler方程,在贴体坐标系中,利用带限制函数的波传播算法,对圆环形燃烧室内的旋转爆轰进行了数值研究,讨论了爆轰传播过程中的波系结构。数值结果表明,即使仅在燃烧室内壁面的部分区域充入可燃气,爆轰波也能够以旋转方式在该区域稳定传播。     

曲面激波诱导斜爆轰的数值模拟

为了研究恰当比预混氢气-空气斜爆轰流场的波系结构和流动特征，基于带化学反应的Navier-Stokes方程，对弹头及楔穿越预混气体时诱导的斜爆轰进行了数值模拟。对流项的离散采用Steger-Warming格式，时间项采用二阶Runge-kutta方法。结果表明，对于弹头：（1）在亚爆轰条件下，能够模拟氢气-空气预混诱导爆轰流场的精细结构；（2）在超爆轰条件下，通过精细调整网格，能够很好地分辨强烈耦合的激波和燃烧波，且与Lehr实验吻合良好。对于楔结构：捕获到了清晰的三波点及其复杂精细的斜爆轰流场结构，预测的诱导燃烧距离、激波角和斜爆轰角与实验吻合良好；通过对流场波系结构变化过程的研究，获得了流场三波点随时间的演化过程。     

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

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

高温区碰并诱导气相爆轰二次起爆的数值研究

针对心形管内气相爆轰二次起爆过程进行了数值研究。采用二阶附加半隐的龙格-库塔法和5阶WENO格式求解二维欧拉方程,采用基元反应来描述爆轰化学反应过程。结果表明:气相爆轰在心形管中传播,由于壁面几何结构的作用,爆轰波熄灭后会发生二次起爆。该二次起爆由高温区碰并所诱导,且在管道不同位置诱导了两处二次起爆。爆轰流场中存在未反应气囊、射流、爆轰波马赫反射、激波与火焰作用及其所诱导的不稳定性等现象。计算得到的二次起爆过程胞格结构演变与实验基本一致。     

爆震管内波与火焰相互作用机理的试验

为了解波与火焰相互作用的机理,在60 mm×60 mm×2000 mm方爆震管内,用乙炔(C₂H₂)和空气混合物进行了单爆震性能研究.利用压力传感器与离子探针同时测得爆震管内的压力和火焰传播速度.根据波和火焰触发的不同时刻来分析爆震管内波与火焰相互作用的过程.结果表明:光滑爆震管内没有产生爆震波,压力波始终在火焰前面;加入阻塞比为0.4的扰流器后,爆震管内产生了爆震波.在扰流器内部,爆震波比较复杂,弱压缩波在火焰前面,但火焰在激波前面;在扰流器出口位置以后,激波在火焰前面.      

发射药在强约束条件下的燃烧转爆轰特性

以Ｄ－４５，Ｓ－１４，３／１三种发射药为研究对象，对钢管约束条件下的高密实火药颗粒床的燃烧转爆轰机理和敏感性进行了实验研究。分别用变片和离子地测量了药床中的压缩波和燃烧波的传播，通过对压缩波阵面和燃烧波阵面的分，发现压缩波形成火焰峰后的已燃区，它将不断加速与加强，最后赶上火焰峰形成强冲击波而诱发爆轰。     

汽油/空气两相脉冲爆震发动机触发爆震的研究

有效的触发爆震波是脉冲爆震发动机正常工作的关键所在。为了在多循环两相脉冲爆震发动机的缓燃到爆震转变过程中,有效的控制激波反射,成功的触发爆震。通过在爆震室内合理的排放障碍物来促进激波与火焰的相互作用,组织激波反射,适时的触发爆震,从而有效地缩短促发爆震的距离,并最终获得充分发展的爆震波。试验发现爆震管内同样的障碍物在不同的位置起到的作用不同。在直径58 mm,长度1275 mm的汽油/空气PDE爆震管内成功触发频率为50 Hz的爆震。      

Experimental research on rotating detonation with liquid hypergolic propellants

This paper describes experimental research into the initiation and propagation of rotating detonation for liquid Nitrogen TetrOxide (NTO) and liquid MonoMethylHydrazine (MMH). An annular rocket-type combustor without nozzle was designed to investigate detonation combustion. The propellants were injected through unlike impingement injectors. The combustion flame fronts and pressure waves were detected using optical diagnostics and dynamic pressure sensors, respectively. The propagation of rotating detonation was established spontaneously by increasing the mass flow rate of propellants. The velocity of propagation of the flame fronts and pressure waves was nearly equal and reaches supersonic speed. Two different detonation combustion patterns are present, single wave mode and double waves mode. And in double waves mode, the two detonation waves are always counter-rotating. The possibility of rotating detonation initiation in a combustor with nozzle was also checked. Stable rotating detonation can be initialized and sustained at similar operating conditions.     

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

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

高温来流下U型脉冲爆震燃烧室燃烧特性试验

为了获得高温来流下气/液两相多循环U型脉冲爆震燃烧室(U-PDC)的燃烧特性,使用汽油/空气为燃料和氧化剂开展了相关试验研究.结果表明:当来流温度为373 K,U-PDC能够在15～38 Hz范围内稳定工作,且随着工作频率提高爆燃向爆震转变(DDT)距离缩短;在室温条件下、工作频率为15～25 Hz时,未能形成充分发展...     

Experimental of combustion instability in NTO/MMH impinging combustion chambers

This paper presents an experimental study into dynamics of chamber pressure and heat release rate during self-excited spinning and standing azimuthal mode in NTO/MMH (nitrogen tetroxide/monomethylhydrazine) impinging combustion chambers. Nine cases including two combustion chamber configurations were conducted. The operating conditions of all unstable cases were located in the instability region according to Hewitt empirical correlation. The results show that chamber pressure oscillations keep pace with the corresponding OH* chemiluminescence intensity over the whole combustion region in the spinning and standing modes. It is indicated that the Rayleigh index is positive over the whole combustion area in all the unstable cases. A significant supersonic flame front structure of the first-order spinning mode was found in a cylindrical chamber, which means that a detonation wave could exist in the cylindrical chamber without a center body. The pressure and heat release rate oscillations at the pressure node are nonnegligible although their amplitudes are lower than those at the pressure antinode in the first-order standing mode with an annular chamber. Besides, the dominant frequency of pressure and heat release rate oscillations at the pressure node is twice as high as that at the pressure antinode.     

环形爆震室中火焰加速数值模拟及验证

对带不同数量孔板的环形爆震室进行数值模拟,并通过试验对数值计算进行验证,来研究火焰加速现象、爆燃向爆震转变过程和不同当量比下起爆距离.数值计算采用二维轴对称非定常Navier-Stokes方程来模拟流体动力学过程.研究发现用较低的点火能量对火焰混合区的可燃气点火产生低速火焰,低速火焰向环形爆震室射流并改变方向向出口传播,火焰在孔板的阻碍作用以及火焰诱导激波和反射波的加速作用下,由层流变为湍流,湍流火焰与其诱导激波相互加强,最终引爆未燃混气;还对爆震波在孔板区的传播过程进行了分析,对不同当量比下的火焰速度和起爆距离进行了模拟研究.