连续爆轰发动机热力学性质的数值模拟

在二维圆柱腔中连续爆轰发动机数值模拟的基础上,对喷入流场中的粒子进行跟踪,得到经历不同燃烧过程的粒子轨迹图.分析影响粒子轨迹及其上物理参量变化的主要因素.数值模拟得到相应的示功图和示热图,计算连续爆轰发动机的机械功及热效率,将数值模拟结果与理想ZND(Zeldovich-Neumann-Doring)模型进行定性和定量的比较,所得数值模拟结果与理论结果符合 得很好,验证了连续爆轰发动机的优势.      

旋转爆轰发动机工作过程的数值模拟

为了深入研究旋转爆轰发动机工作过程,采用高精度高分辨率迎风通量分裂格式(AUSMPW+)、氢氧7组分8反应化学动力学模型,求解三维带化学反应的Euler方程。对旋转爆轰发动机工作过程进行了数值模拟,分析了采用预爆轰管点火过程的流场特征以及整个发动机推力、推力偏心距、侧向力随时间的变化规律。计算结果表明:切向预爆轰管紧贴喷注入口时,点火过程仅产生一道单向爆轰波,能成功点燃旋转爆轰发动机;在本文给定的计算条件下,旋转爆轰发动机平均推力约180N,旋转爆轰波传播频率约为14285Hz;旋转发动机正常工作过程中,推力偏心矩、侧向力随时间周期性变化。     

气云爆轰

气云爆轰实验研究在一根专门设计的管中进行，通过对压力和火焰阵面的测试，研究激波的成长与发展过程，气云爆轰的爆速及其影响因素以及激波扣的点火诱导时间。理论研究考虑了液滴在高温，高压和高速气流中的变形，剥离和破碎，液滴的点火和局部爆炸以及气云爆轰的Ｃ－Ｊ判据。     

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

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

双爆轰管脉冲爆轰发动机近场噪声特性的试验研究

为研究双管脉冲爆轰发动机近场爆轰噪声特性,设计了双管脉冲爆轰发动机爆轰噪声试验测试系统,对不同填充率、不同管间距下的爆轰噪声进行了测量。试验结果表明:随着传播距离的增加,爆轰噪声峰值声压逐渐减小。PDE管口附近,爆轰噪声峰值声压衰减速度最快;距离管口较远位置处,衰减速度逐渐减小。随着填充系数的增加,爆轰噪声逐渐增大,指向性逐渐变缓;距离管口400mm处,60%,100%和140%填充系数下峰值声压级依次为197.95,201.93和204.51d B。不同传播距离处,爆轰噪声指向性不同。当距离管口较近时,爆轰噪声最大值出现在0°方向;距离管口大于2m的区域内,最大值出现在30°方向。不同管间距情况下,爆轰噪声呈现出相同的指向性。管间距越小,爆轰噪声的指向性越明显。     

PIV中提取速度信息的一种新方法

通过对ＰＩＶ中粒子图像固有特点的研究，提出了一种新的速度信息提取方法。该方法先通过模拟人眼判断粒子图像中心点的过程来获得粒子的中心点，再利用这些中心点之间的相关关系求取速度。以内烯机缸内流场的ＰＩＶ实测图像为例，通过与直接ＦＦＴ变换求相关的方法进行比较，阐明了该方法的特点。     

喷管流场V3V试验示踪粒子特性研究

采用体3维测试系统(V3V)试验研究了拉伐尔喷管出口的速度场。为了选择适合V3V试验测量气流流场的示踪粒子,对3种不同物理属性的固体微粒进行测试分析,结果表明:二氧化钛粒径为5～10μm,对于激光散射性太差,粒子亮度较弱,在V3V处理软件识别时粒子匹配不成功;滑石粉和聚苯乙烯的粒径分别为20～25μm和30～35μm,在V3V软件处理4副图片时的粒子匹配率分别为30%和40%,试验得到的速度场与理论计算值误差分别约为2%和5%。     

复合固体推进剂危险性分析

一、引言 复合固体推进剂不单是一种能按一定规律燃烧的物质,而且潜在着爆炸或爆轰危险性,特别是高能复合推进剂,硝胺推进剂和改性双基推进剂的危险性就更大了。因为,就其化学本质来说,固体推进剂与炸药本质上很难区别,它的能量密度和反应热可与常规高效炸药相匹敌,甚至更高。因此,关于复合固体推进剂及其发动机的偶然性爆炸或爆轰问题,引起了人们的关注。     

氢氧爆轰直接起始的射流点火方法研究

为探讨氢氧爆轰驱动激波风洞中的高氢混合比氢氧混合气直接爆轰的点火方法，对射流点火管的点炎能力进行实验研究发现，点火管长度是进行射流点火能力的关键因素，射流点能力随点火管长度的增加而增强。点火管与爆轰管之间无隔膜，结构简单，操作方便。     

神经网络在立体成象测速技术中的应用

立体成象测速技术是一种用于测量三维流场中流体的三个速度分量的光学方法。该测量为非接触测量,其测量原理是基于由两个CCD传感器从不同的有利位置点获取的图象。采用这一方法时,如果流场中散布的粒子的密度较大的,由成象系统所获得的图象中就可能丢失单个粒子的部分象点或者是等效的数据点。该数据丢失大多数是发生在对重叠的粒子图象进行分解的过程中和对粒子进行跟踪阶段。为了尽可能恢复更多的数据点并提高测量精度,在立体成象测速技术的两个阶段中采用了神经网络技术。在对重叠的粒子进行分解的阶段采用了后向传播(Back Propagation)神经网络,这是由于此方法具有图形识别和分类的能力。在对粒子跟踪阶段,跳跃神经网络（Hopfield neural network)是获取相应粒子轨迹的一种公认的最佳求解方法。研究表明,神经网络对于提高测试性能具有潜在的作用,而且已被证明对于立体成象测速技术是非常有用的。     

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

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

狭缝内爆轰波传播模式的实验研究

为研究狭缝内爆轰波起爆及各种传播模式的特点,实验得到了不同初始压力下(p0=10~60k Pa)化学当量比的乙烯/氧气预混气体在高度1.0mm,宽度20mm狭缝内的爆轰性能。采用烟膜片记录爆轰波运行轨迹,高速摄影捕捉火焰面。结果表明:随着初始压力的降低,实验依次得到稳定爆轰模式、"结巴"爆轰模式、驰振爆轰模式和低速稳定传播模式,中间两者为不稳定模式。不稳定爆轰传播模式具有强烈的速度震荡,在一个周期内,烟膜板中周期性的出现过驱爆轰的胞格结构,随后"结巴"爆轰以单头形式与侧壁面发生碰撞并向前传播,驰振爆轰短暂熄灭。低速稳定传播速度为45%的C-J速度,烟膜板中无胞格结构出现。狭缝内爆轰波速度亏损值为10%~15%,高于大尺寸通道中的速度亏损。     

进气掺氢对爆轰波参数影响的试验研究

为研究氢气对气液两相脉冲爆轰发动机内爆轰参数的影响,对进气掺氢的脉冲爆轰发动机进行了试验研究,实现了发动机的稳定工作。研究发现,进气中掺混氢气对爆轰波参数有较大影响,氢气进气压力为0.1MPa时,掺氢后爆轰波的峰值压力和速度均明显提高,爆轰波速度由1045.9m/s升高至1846.2m/s,发动机尾部位置上爆轰峰值压力从1.9MPa上升至4.5MPa,爆轰波后平台压力提高了0.2~0.6MPa;当氢气进气压力由0.1MPa增加至0.22MPa时,爆轰波速度随之增大。进气掺氢后燃烧转爆轰位置提前,有利于减少燃烧转爆轰的距离,进气掺氢的脉冲爆轰发动机能在较少扰流片(5片扰流片)情况下起爆,但当扰流片少于5片时,发动机内未能正常起爆。     

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

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

连续旋转爆轰发动机流场三维数值模拟

对连续旋转爆轰发动机(continuous rotating detonation engine,简称CRDE)流场进行了三维数值模拟.数值计算获得同轴圆管腔内CRDE的多个爆轰循环过程.依据计算获得流场,分析了可燃气入射、提前燃烧、爆轰波结构等实现连续爆轰的几个关键机理问题.通过多个算例对不同来流总压时的旋转爆轰的推进性能进行了比较与分析,最后计算获得基于燃料的比冲约为2 200 s,燃料流量随总压变化呈线性增长.      

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

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

CE/SE方法在非定常爆轰计算中的应用

建立了非定常爆轰波的理论模型。根据CE/SE方法的基本思想，推导出适合于求解二维N-S方程的长方形网格的CE/SE方法。对于包含刚性反应源项的守恒方程组，将其分为两部分，其中刚性反应源项单独求解。本文以爆轰管内装填C3H8/O2预混气体为例，数值模拟管内爆轰波的传播过程以及爆轰波脱离爆轰管后流场的图谱。计算结果显示，爆轰波在管内传播时，其波系极其复杂。当爆轰波离开出口后，迅速熄灭为非反应的激波。管外不仅有激波、膨胀波的传播，而且还存在激波与涡的相互作用。数值实践表明，该CE／SE方法可以有效捕获爆轰波与激波等强间断。     

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

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

高温真实气体效应中催化效应对气动热影响的实验探索

主要介绍了在氢氧爆轰驱动高焓激波风洞中进行真实气体效应中催化效应对气动热影响的实验研究。首先从测试工作的角度，论述了测热模型、测热传感器及风洞改进等关键技术问题及其解决方法。其次描述了在氢氧爆轰驱动高焓激波风洞中，首次开展气动热风洞试验的过程及其初步结果。结果表明：热流数据随测点位置和迎角的改变呈有规律的变化；在同样条件下，完全催化表面比完全非催化表面热流数值有明显增加的趋向。     

当量比对连续旋转爆轰燃烧室特性的影响

为了探究当量比对甲烷-空气连续旋转爆轰燃烧室（CRDC）特性的影响,利用二维可压缩欧拉方程对CRDC进行了数值研究,分析了爆轰波的发展过程和贫燃熄火过程,对比了不同工况下CRDC特性参数的变化情况。结果表明:CRDC起爆后燃烧场在由不稳定状态到相对稳定状态的过程中发生了2次碰撞,当进气当量比较低时,CRDC未能完全发生2次碰撞过程就已经熄爆。随着进气当量比的降低,爆轰波传播速度、轴向平均速度、出口平均温度、出口平均总压均呈下降趋势;增压比随当量比降低而减小的根本原因在于旋转爆轰燃烧过程和等压燃烧过程的熵增差减小,使吉布斯自由能增量差减小。CRDC的燃料驻留时间处于亚毫秒量级,燃烧热效率保持在99%以上。