驻定斜爆轰波起爆条件与结构的数值模拟研究

为了探究驻定斜爆轰波的起爆条件以及来流马赫数和楔面角对斜爆轰波形态结构的共同影响,基于二维欧拉方程和单步化学反应模型,对驻定斜爆轰波进行数值模拟研究。结果发现,斜爆轰波的起爆和结构受到来流马赫数和楔面角的显著影响,来流马赫数越小,形成斜爆轰波的楔面角条件越苛刻。随着来流马赫数和楔面角增加,波后混合气体温度升高,爆轰反应变得更加剧烈,使得斜激波诱导区域变短,未反应气团减少,接触间断面更靠近楔面,下游出现的涡结构也减少。     

预混超声速气流斜激波诱导脱体爆轰研究

为了研究高静温预混超声速气流中爆轰的机理问题,利用自主设计的连续式加热器来产生马赫数为2.6,静温大于700 K的氢/空气预混超声速气流,并采用高速纹影技术来观测斜激波诱导脱体爆轰的直接起爆和发展的动态过程.研究表明,实验得到的起爆前斜激波和起爆后脱体爆轰波的角度与理论分析结果非常一致.并发现由激波诱导爆燃向激波诱导脱体爆轰的转变过程非常快;而且对于当前初温较高当量比较低的混合物,会出现爆轰波突然熄灭与重新起爆的现象.当前研究为爆轰燃烧在高超声速推进中的应用提供了重要参考.      

气相爆轰基元反应模型数值模拟

采用基元反应模型和高精度的ENO格式对气相爆轰进行数值研究，对H2／O2／Ar混合气体起爆和爆轰波传播过程的数值模拟结果表明，计算的爆轰波阵面参数和实验相当符合。数值研究的结果还表明，爆轰波反应区中参与反应的不同组分具有不同类型的变化特征。计算结果的精度随着网格尺寸增加而增加，并能保持较好的收敛性。     

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

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

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

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

激光支持爆轰波的仿真研究

采用了爆轰波的C-J起爆模型,计算了激光支持爆轰波初始参数,得到激光功率密度与爆轰波波面初始参数的关系.利用流体力学的方法建立爆轰渡在靶与约束层双重约束作用下的二维辐射膨胀模型,用Fluent软件进行仿真研究,得到了激光冲击强化后不同时间对应的爆轰波压力和轴向速度分布,仿真结果与激光冲击强化GH742试件后表面残余应力分布的实验结果一致.     

旋转爆轰燃烧室结构及稳定性研究

基于二维可压欧拉方程,对充有当量比为1的氢气/空气预混气的不同结构燃烧室进行数值计算。研究了爆轰波在不同坡度燃烧室内的传播过程,分析了燃烧室内不同波头数目对其性能的影响。研究表明:在燃烧室入口端采用适当的楔形结构,有助于提高沿单一方向传播的爆轰波的强度,提高燃烧室做功能力;当燃烧室周向尺寸较大时,同时起爆双波头有助于提高其工作的稳定性。     

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

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

倾斜环缝喷孔式连续旋转爆轰发动机试验

为了验证连续旋转爆轰发动机起爆过程和传播机理的理论研究,采用了1.0mm,1.6mm以及2.0mm三种尺寸的倾斜环缝喷孔结构,以H_2/O_2切向喷注的预爆轰管进行起爆,进行了H_2/Air组合方式的连续旋转爆轰发动机试验研究。试验结果表明,1.6mm环缝获得了最佳的试验结果,爆轰波频率为5.0~5.3kHz,对应的传播速度为1460.1~1547.7m/s。通过对旋转爆轰波传播机理的分析,验证了爆轰波存在三种传播方式:正转、反转和双波头对撞。     

基元化学反应一维爆轰波的数值模拟

应用基元化学反应模型对一维氢氧混合气体爆轰波进行数值模拟。采用一端高温高压点火起爆的方式,模拟爆轰波的传播。对一维爆轰波的结构、传播特征以及各个化学组分的变化趋势进行了研究。通过将模拟的结果同国内外的研究结果进行对比,发现了一些规律,为多维爆轰波的数值模拟和脉冲爆轰发动机的研究奠定一定的基础。     

脉冲爆震发动机爆震管壁换热及冷却的实验

利用Pt100热电阻对脉冲爆震发动机(PDE)在有、无水冷状态下爆震管壁沿轴向的温度分布进行了实验研究.实验结果发现爆震管同一轴向位置处,爆震外管壁温度随爆震频率的增加而增加.同一爆震频率时,距推力壁位置越远,爆震外管壁温度相对越高.向爆震管换热系统中通入冷却水之后,爆震外管壁温度有明显的降低,并且爆震外管壁温度梯度也有一定的下降,20Hz时加水冷之后,PDE工作16s时,爆震外管壁温度最高可降低106.5℃.      

爆震管同一截面处压力测量试验

为得到爆震管内同一截面上压力波的分布特性,在60mm×60mm×2000mm方爆震管内,用乙炔(C₂H₂)和空气混合物进行了单爆震研究.在距离封闭端1400,1600mm两个截面的不同位置处分别安装4个高频响压力传感器测得压力值.实验在光滑爆震管和安装扰流器的爆震管内分别进行,结果表明:①同一截面的压力波大小并不一样,安装扰流器的爆震管内得到充分发展的爆震波,且爆震波不是一个平面波;②壁面处压力波的压力大于爆震管道中间压力.      

球形点火腔对脉冲爆震DDT过程影响的数值研究

为了减小脉冲爆震发动机PDE（pulse detonation engine）单次循环时间,缩短爆燃向爆震转变DDT（deflagration to detonation transition）距离是关键.为此,提出了一种在头部添加球形点火腔的新型爆震管道结构设计方案,并采用丙烷和空气为可爆混合物,通过对新型爆震管DDT过程的二维数值模拟,研究其对DDT距离以及DDT时间的影响.数值模拟结果表明,头部采用球形点火腔后,点火腔中的压缩波经过多次反射后,能够在爆震管中更快地促使爆震波形成,DDT距离大大缩短;当球形点火腔直径为1.5倍等直爆震管直径时,相对于常规等直管爆震管结构,其DDT距离和过程时间分别减少了14%和16.26%.      

预混气爆震波在不同直径管间的传播特性

在管径20mm与60mm的两管间研究了预混气爆震波的传播特性,着重测量了爆震波传递过程中的压力与速度变化,并分析了预混气压力、浓度、当量比以及连接段尺寸等参数对爆震波传播的影响与规律。研究发现:小管中没有爆震时大管中仍然可以出现爆震波,小张角的连接段对爆震波的再触发有利,再触发后的爆震波其压力、速度受到预混气压力、浓度及当量比等参数的影响。      

爆震室压力测量可靠性试验

为了提高脉冲爆震发动机(Pulse Detonation Engine,简称PDE)内爆震波压力测量值的可靠性,解决多循环PDE爆震室压力测量问题.通过动态压力测试系统,对单次爆震和多循环PDE爆震室内的压力进行测量.单循环爆震波压力测量的比较性试验研究发现:当主爆震室的旁支套管长度小于100mm左右,可以测到爆震波的峰值;当套管长度大于150mm左右时,火焰不能在旁支小管内传播,传感器测到的只是爆震波后的平台的压力.当通过冷却套管来测量多循环PDE爆震室压力时,采用较短长度的套管,可以得到爆震波压力.研究结果解决了多循环爆震波的压力测量问题.      

脉冲爆轰发动机预爆轰点火数值模拟

本文采用二维数值模拟对脉冲爆轰发动机点火方式中的预爆轰点火进行了分析研究。模拟了预爆轰管在主爆轰管上的两种最基本的布置方式：平行布置方式和正交布置方式。讨论了预爆轰管的布置方式和几何形状对点火的影响。结果表明：正交布置方式优于水平布置方式，因为反射激波在点火过程中起到了非常重要的作用。预爆轰管和主爆轰管的内径比是影响点火的重要参数，而长度的影响则较小。     

弯曲方管内爆燃向爆震转变特性

为缩短爆震管轴向尺寸,采用曲管爆震管替代直管爆震管.基于弯曲方管实验器,进行了二维数值仿真计算,研究了点火位置、点火能量及弯曲段曲率半径对爆震波形成及传播的影响规律.数值模拟结果表明:爆燃向爆震转变(DDT)时间随弯曲段曲率半径增大而减小,DDT距离在弯曲段曲率半径最大时达到最小值;点火能量越大,爆震管内激波强度越大,DDT时间越小;在距离爆震管进口300mm处点火时,爆震管内爆震波压力峰值最大,DDT时间最短.     

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

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

狭缝内乙烯/氧气预混气体爆轰几何极限的实验

为研究狭缝内爆轰波传播极限,实验得到了不同初始压力(0.004~0.04MPa)下化学恰当比的乙烯/氧气预混气体在狭缝高度为1.0~4.0mm狭缝内的爆轰性能.采用烟膜板记录爆轰波运行轨迹,高速摄影捕捉火焰面.结果表明:狭缝高度越小,爆轰极限对应的临界初始压力越高.对于近极限条件内不稳定爆轰传播模式,包括"结巴"爆轰和驰振爆轰,其对应的初始压力范围随着狭缝高度的降低而变宽.考虑初始和边界条件,将水力直径与胞格宽度之比作为合理的爆轰敏感性参数描述预混气体爆轰特性,得出在不同狭缝通道内爆轰极限范围为,即该比值在0.326~0.403之间.      

共用喷管多管爆震发动机工作过程数值模拟

为研究共用喷管多爆震管间爆震流场的相互作用,利用CFD商业软件,选用通用有限速率模型,以氢气和氧气混合物为例,对不同喷管喉部面积和不同爆震管数的共用喷管多管爆震发动机工作过程进行了数值模拟,并对其性能进行了比较与分析。计算结果表明,喷管喉部面积和爆震管数目对爆震管间的相互影响产生很大影响,模拟结果可为多爆震管脉冲爆震发动机的设计提供理论基础。