当量比和间隙尺寸对爆震波传播过程的影响

为探索微型脉冲爆震发动机推进系统的可行性,进行了微空间内爆震极限特性的研究.通过对预混氢气/氧气爆震火焰在平板狭缝中的传播过程进行测量,分析了爆震波压力、速度随当量比和间隙尺寸的变化趋势,发现在极限范围外,来自稳流段的爆震波都能在间隙内通过压力和速度的自我调整后达到稳定状态,调整过程所需的距离随间隙尺寸减小而增长.根据爆震波速度衰减的定义,给出了激波和火焰锋面速度沿间隙通道方向的变化特点,观察到了微尺寸下爆震波的4种传播模式:稳定爆震波、准稳定爆震波、低速爆震波和非爆震波.      

爆震管内爆燃到爆震转捩过程的实验研究

基于常温常压轴向进气条件下气动阀式两相爆震发动机模型开展爆燃到爆震转捩过程的研究，利用离子探针和高频压力传感器测量了管内火焰传播速度和燃烧压力波。研究结果表明，管内爆震波形成的初始位置在正常火焰面下游的一定位置处，如假设爆震管内波和火焰面是一维的，则在爆震波产生的最初一段时间内，爆震管内将出现三道火焰面：往出口方向传播的正常火焰面和爆震波火焰面及往进口方向传播的回传爆震波火焰面。     

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

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

点火方式对旋转爆震发动机工作特性的影响

为研究并改善旋转爆震发动机的点火-起爆性能,进一步深入了解点火方式对旋转爆震发动机工作特性的影响,采用小能量火花点火装置、高能量火花点火装置、爆震管三种点火方式,进行了一系列旋转爆震发动机起爆实验,发动机采用环缝-喷孔对撞式掺混方式,燃料为H2,氧化剂为空气。实验对比研究了不同点火方式下旋转爆震发动机的点火起爆性能及点火方式对发动机工作特性影响。实验结果表明三种点火方式均成功起爆旋转爆震波,并周向稳定传播。通过对旋转爆震波起爆过程详细分析发现,不同点火方式引燃的不同初始状态的火焰均需在环形燃烧室经历一个类似DDT的火焰发展过程才能成功建立爆震波,且火焰发展过程的时间间隔表现出很强的随机性,但总体来看爆震管点火时爆震波建立时间较其他两种点火方式短。此外,该工况条件下三种点火方式起爆发动机时其工作状况可重复性均可达100%,稳定工作过程中的传播特性与点火方式无明显关系,爆震波传播频率较为稳定,在5437~6440Hz波动。     

旋转爆震发动机爆震波建立过程实验研究

为了深入研究旋转爆震发动机爆震波建立过程及形成机理,采用小能量火花单次点火的方式进行了一系列旋转爆震发动机起爆实验。发动机采用环缝-喷孔对撞式掺混方式,燃料为H2,氧化剂为空气,实验成功起爆旋转爆震波,并连续旋转稳定传播,爆震波传播频率为5.09~6.45k Hz,传播速度为1286~1644.8m/s。在发动机稳定工作过程中,集气腔与燃烧室相互影响,二者处于平稳的动态平衡。其次,通过对旋转爆震波起爆过程详细分析发现,点火形成的初始火焰在环形燃烧室经历一个类似DDT的火焰发展过程,成功转变为爆震波,且从点火到爆震波建立之间的火焰发展传播过程和时间间隔均表现出很强的随机性。此外,为验证小能量火花点火的可靠性,还进行了小能量点火重复性实验,发现在稳定工况条件下采用小能量点火成功率最高可达100%,各组旋转爆震波传播速度在1440m/s附近波动。     

环形爆震波聚焦起爆

本文对三种带不同结构形式的环形预爆管的爆震发生器进行了数值模拟,来研究环形爆震波聚焦起爆现象及其气动特性。数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动的Navier-Stokes方程来模拟化学动力学和流体动力学过程。研究发现用低的点火能量对环形预爆管中的燃料和氧化剂点火产生层流火焰,层流火焰在狭窄管壁的作用下完成爆燃向爆震转捩,形成环形爆震波,环形爆震波在聚焦腔经过反射、汇聚,最终引爆主爆震室中的可燃混合物;并对稳定爆震的距离和强度进行了对比研究。通过分析对环形爆震波衍射、反射、聚焦的气体动力学特性及复杂波系发展规律有了更深的认识,为进一步试验提供参考。     

斜爆震波起爆特性及其波系结构研究

本文通过对斜爆震波起爆特性及其波系结构的研究进展进行综合评述,分析其研究现状与存在的不足,指导未来斜爆震波的基础研究。首先总结了斜爆震波起爆和驻定特性研究进展,重点回顾了斜爆震波起爆准则和驻定条件,接着对斜爆震波的诱导区结构和波面胞格结构研究进行了总结,随后对斜劈构型对斜爆震波的影响研究现状进行了概述。综述表明,斜爆震波的起爆、驻定和波系结构还需持续开展研究,特别需要关注的是平面斜爆震波的起爆准则和驻定条件,在真实流动条件下更深入研究斜爆震波宏观结构和波面结构的形成机理,探索受限空间内斜爆震波起爆、反射和驻定特性。     

圆盘结构下旋转爆震波的不稳定传播特性

为研究圆盘结构下旋转爆震波的不稳定传播特性,以2H2+O2+3.76N2为反应混合物,开展该结构下的二维数值研究。详细分析了旋转爆震波不稳定传播模态的流场特征,以及爆震波参数、出口流场参数和增压比的变化。结果表明,不稳定传播模态下,旋转爆震波重复进行"解耦—再起爆"过程,内圆扩张曲面的发散作用,使出口附近的爆震波首先解耦,解耦区域逐渐往流场内部扩张;反射激波与前导激波的碰撞产生局部热点,促使爆震波重新起爆;爆震波的压力、温度以及传播速度皆随爆震波的解耦及再起爆过程发生变化;旋转爆震波的不稳定传播对出口速度分量及马赫数的影响很小;出口增压比受不稳定传播的影响较大,随时间呈周期性变化,变化幅值较高且不稳定。出口增压比的循环周期与爆震波"解耦—再起爆"的周期一致,稳定后两者的循环频率约为21.6kHz。     

横向射流起爆爆震波数值研究

采用CFD软件对横向喷射的射流在预先充满丙烷/空气恰当比混合物的爆震室中起爆爆震波的过程和机理进行了数值模拟,并讨论了射流的压力、速度以及温度对爆震波起爆特性的影响.结果表明,对于计算的物理模型,爆震射流不能在爆震室中直接起爆爆震波;射流与壁面的碰撞形成激波反射,激波反射产生的横渡、热点和局部爆震引发起爆;能够成功起爆爆震波的横向射流存在最小的射流压力,该射流压力为0.9MPa;较小的射流压力、较小的射流速度或较低的射流温度都不利于横向射流在爆震室中起爆爆震波.     

激波绕射触发爆震波的试验研究

在单爆震试验平台上对激波绕射触发爆震波的特性进行了试验研究.针对不同的激波反射器结构尺寸,测量了爆震管内气体压力变化历程及火焰传播速度,获得了不同工况下激波绕射触发爆震波的特性.结果表明:激波反射器对爆震波的形成有促进作用,并取决于反射器中的反射锥及反射孔板的结构与尺寸;在反射锥堵塞比一定的情况下,存在最佳反射锥锥角使得爆震波的形成距离最小;而反射锥与反射孔板的间距增大,爆震波形成距离增大,一旦形成稳定的爆震波后,反射器对爆震波压力的影响不大.      

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

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

环截面管中的爆燃转爆震实验

为了实现在环截面爆震管封闭端同时点火,采用了先在预燃室中点火,再通过八路支管将火焰引入到爆震管封闭端,各个支管长度一样,因此可看成在环管封闭端多点同时点火.实验结果显示,在爆震管封闭端火焰汇聚形成激波点火,可在环管封闭端附近直接产生爆震波.在爆震管中安放扰流片时,流动阻力使爆震速度偏低,无扰流片时,爆震速度接近理论值.不加扰流片当量比在1.0左右,火焰速度会出现爆燃或爆震两种随机状态,当量比到1.1时可避免这种状况.      

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

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

热射流点火对爆震管内火焰加速及爆震波触发影响的数值研究

为研究热射流点火对爆震管内火焰加速及爆震波触发的影响,运用34步26组分丙烷基元反应进行了二维数值模拟,获得了5种不同的热射流发生器几何结构下的爆震管内火焰传播规律及缓燃向爆震转捩(deflagration to detonation transition,DDT)的时间与距离.结果表明:爆震管中湍流在火焰加速阶段起着重要作用,在爆震波触发阶段激波与火焰相互作用占有主导地位.根据DDT过程的定义,得到DDT时间在1.4~2.0ms之间,同时发现热射流发生器长度为150mm,热射流发生器孔径在8mm时DDT时间最短,热射流发生器长度及孔径对DDT距离的影响不大.      

基于自适应网格加密的超声速可燃气热射流起爆详细反应数值模拟

采用块结构自适应网格加密开源程序AMROC,在高性能计算集群中,进行氢气、氧气、氮气详细反应机理的二维超声速热射流爆震起爆自适应网格精细数值模拟,研究在热射流的持续喷射作用下超声速可燃气热射流爆震起爆,以及形成的精细的爆震胞格结构.结果表明:超声速可燃气流中热射流类似气动斜劈,形成局部激波诱导燃烧.超声速可燃气中持续的热射流喷射会形成过驱爆震,并导致不规则爆震胞格的生成.热射流的扰动压缩作用对过驱爆震的形成以及不规则爆震胞格的产起着关键的作用.热射流的扰动压缩波以当地声速通过爆震波后的亚声速区域作用于爆震波,使爆震波处于持续过驱状态,并形成不规则的爆震胞格结构.     

突跃型与平滑型斜爆震波起爆机制数值模拟

采用带有化学反应的Euler方程,对突跃型与平滑型2种形态的斜爆震波(ODW)在起爆机制和结构特点等方面的差异进行了研究.结果表明:斜激波(OSW)后是否存在亚声速区是判定斜爆震形态的依据.当波后火焰抬升斜激波面使得亚声速区存在时,横向激波与三波点结构就会出现并形成突跃型斜爆震.当波后所有区域皆为超声速区时,波后火焰燃烧无法影响到上游激波面,则会形成平滑型斜爆震.通过斜激波关系式给出了横向激波形成所需的临界斜激波角度,只有当斜激波角度大于此临界角度时才能形成突跃型斜爆震波,反之则形成平滑型斜爆震波.      

超声速预混气的热射流起爆过程数值模拟

采用带有化学反应的Euler方程,对超声速预混气的热射流起爆过程进行了研究。在来流速度大于CJ(Chapman-Jouguet)爆震速度的情况下,射流所导致的弓形激波在爆震波起爆过程中起到了重要作用,而爆震波不稳定性导致横波的产生是决定能否完全起爆的关键因素。射流角度的改变能够影响到起爆的过程与驻定性质。随着射流角度的逐渐减小,会出现起爆后驻定、起爆后不能驻定以及起爆失败三种结果。     

斜爆轰波的波角和法向速度-曲率关系初探

为研究斜劈诱导斜爆轰波的波阵面弯曲效应，以期为斜爆轰的不稳定性及其演化规律提供新的见解，基于加权本质无振荡（WENO）格式空间离散和附加Runge-Kutta方法时间离散的求解器，针对不同的化学反应参数（释热量、放热速率和化学反应区参考长度）条件，开展斜爆轰波的数值计算研究。结果表明斜爆轰波沿波阵面的波角变化可分为3个区域：区域I，波角平滑减小；区域II，波角跃升后衰减；区域III，波角有规律振荡。波阵面法向速度-曲率关系在区域I呈现准垂直直线变化趋势，并伴随着爆轰波强度的不断衰减；在区域III则呈现出"D"形曲线，即由极曲线段、光滑水平变化段和拟线性变化段组成，为类胞格结构的周期性演变；区域II可认为是以上两个区域特征的耦合。不同的化学反应参数对斜爆轰波波阵面的弯曲效应影响存在较大差别。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。