增强型喷射器对爆轰波DDT过程的影响

饶飞雄，雷知迪，丁 珏，翁培奋

(上海大学应用数学和力学研究所,上海 200072)

脉冲爆震发动机(pulse detonation engine, PDE)因其设计简单、热效率高等优点，未来有望替代传统发动机，因此越来越受到研究者们的重视[1]。

DDT混合增强技术是研究者们关注的焦点，是发展超声速飞行器推进系统的关键问题[6]。DDT中非常重要的2个参数是DDT时间和距离。Johansen等[7]通过研究激波与火焰相互作用产生压缩效应来加快DDT过程；Gamezo等[8]通过分析障碍物对爆燃转爆轰的作用得出障碍物最佳间距；张宝亮等[9]基于有限体积法分析障碍物与爆轰波的相互作用规律；王治武等[10]研究点火对两相无阀脉冲爆震发动机性能的影响，分析点火起爆对爆燃转爆轰时间的影响；Khokhlov等[11]分析出热点的形成可以促进火焰速度的加快，最终促进爆轰的形成；De Witt等[12]采用孔板与圆盘的组合装置通过激波发射形成热点成功触发爆震。但此试验中孔板的堵塞比较大，所形成的阻力较大，不太适合工程应用；马丹花等[13]采用二维黏性CE/SE方法研究扰流片对燃烧转爆轰的影响，得出扰流片最佳片数和间距；朱雨建等[14]采用孔栅干涉方法发现激波压缩效应和外界扰动对DDT起着十分重要的作用；刘云峰等[15]通过研究激波-爆震波转捩发现良好的点火位置和反射激波方向可以促进DDT过程。相关文献设计的扰流片可以促进DDT过程，但需要多个扰流片方能达到相应的效果，这对设计和制造相应的装置有着较高的要求，同时延长了装置的长度。

为了能在脉冲爆震发动机内高效点火并实现稳定爆轰，本文中设计一种新的喷射器，并开展爆燃转爆轰过程的数值模拟。讨论热点引起DDT的机理，分析喷射器表面阻断率及放置位置对DDT的影响，以期为进一步优化喷射器结构提供理论依据。

1 流场控制方程和数值方法

1.1 控制方程

对于爆震室可燃混合气体爆轰和爆轰波传播过程，采用Euler控制方程：

(1)

(2)

对于一个有N种组分和L个基元反应组成的系统，其反应机理可以表示为：

(3)

多数城市河流都有美化城市环境的功能，甚至以城市标志的地位存在，而一旦发生污染，会直接影响城市形象。为此，必须借助科学合理的方式对河水污染情况进行处理，促使其恢复原有功能，而在处理技术的选择上，需要结合城市河流污染实际情况。多元化处理就是综合多种治理技术对城市河流水体进行处理，包括水体净化功能恢复技术、人工生物处理技术等，但应用过程中需要充分考虑经济性原则。

(4)

1.2 数值方法

在计算对流通量时采用Roe格式：

(6)

式中：下标R、L代表单元面右、左的状态量。H为通量矢量，W为解向量，δW表示变量在单元面右左两侧值的差。A称为 Roe 平均矩阵。为了提高数值精度，对参量在单元面左右两侧的值通过Muscl插值得到。时间推进采取稳定性好、计算效率高的LU-SGS方法，时间步长采用10-8s来保证数值计算收敛与稳定性。

1.3 网格尺寸、初始条件和边界条件

爆轰波的化学反应速率是普通燃烧的103～108倍[16]。考虑二维爆震室中的H2、O2当量比为1[17]，采用二维直管模型，初始压力和初始温度分别为105Pa和300 K。采用二十步基元反应，并基于0.1 mm网格尺寸数值模拟爆轰波形成和发展过程。图1所示为41.7 μs时刻计算所得的爆轰波结构。表1给出了数值计算所得的爆轰参数与实验结果[18]、C-J理论值[18]的比较。结果显示：与C-J理论值相比，计算所得的爆速、温度、压力值与实验结果更接近，其与实验值的相对误差分别为0.20%、2.78%和2.69%，符合精度要求，说明采用的模型和数值方法是合理的。

图1 爆轰波结构Fig.1 Single detonation cell pattern

表1 验证算例中爆轰参数Table 1 Detonation parameters in verification example

2 喷射器增强燃烧转爆轰过程的数值研究

2.1 喷射器基本构型

在喷射器表面布置半径相同的孔。定义：表面阻断率φ为半球形喷射器表面上阻碍流体流动的面积与喷射器总面积之比。经过计算，喷射器各孔对应的表面阻断率分别为0.400、0.444、0.462和0.471。坐标原点位于左端壁面的中心，爆震室中喷射器沿着中轴线x方向放置。

爆震室模型的直径为2 cm，点火直径为1.6 cm，如图2所示[19]。划分网格时，在喷射器附近区域进行加密处理，最小网格尺寸为5×10-3mm。计算时，右端壁面为压力出口，其余为壁面边界条件。流场初始压力为105Pa，初始温度为300 K。

图2 计算区域和喷射器示意图Fig.2 Sketch of computation field and injector

充满气体C2H2、O2、Ar(初始化学当量比φ(C2H2)∶φ(O2)∶φ(Ar)=1∶2.5∶31.5的爆震室流场开展数值讨论。其中：单步反应的指前因子取3.66×1010kmol/(m3·s)，温度影响因子取0，活化能取1.26 ×108J/kmol[20]。

2.2 喷射器触发爆震波的形成

缩短燃烧转爆轰所需时间和转化的距离是优化喷射器的关键指标。表2给出了爆震室未安装喷射器以及喷射器安装在x=9.0 cm位置流场的参数值。可以看出：安装喷射器情况下，218 μs时刻爆震室形成了稳定爆轰波。由此可知，爆震管中安装喷射器可以快速发生燃烧转爆轰过程。

表2 有/无喷射器时爆轰波状态对比Table 2 Comparison of detonation wave state with or without injector

2.3 燃烧转爆轰DDT过程

图3所示为经过喷射器时火焰的发展历程。前导激波作用于装置的壁面并发生反射，大量燃料聚集在壁面附近，浓度升高，反应速率增大，喷射器附近流场的温度也随即升高，形成热点。

图3 喷射器影响下燃烧转爆轰的DDT过程Fig.3 DDT process affected by the injector

为了清晰显示喷射器附近流场密度、压力和温度的分布，对计算结果进行无量纲处理(选取未受扰动处流场参数作为无量纲基准尺度)。图4所示为101 μs时刻，无量纲化流场的压力、密度、压力和温度沿中轴线的分布。可以看出，101 μs时刻气体压力急剧升高。分析原因，主要是经过喷射器之后的火焰发生汇聚，发展的火焰与前导激波耦合在一起，使得气体的混合和燃烧得到增强，流场的压力快速升高。图5给出火焰沿x正方向传播速度随时间的发展历程。可以看到：经过喷射器之后，火焰速度在较短时间内升高至781.375 m/s，该值接近稳定爆轰速度的1/2，表明此时刻燃烧已受到喷射器的激励作用[21]，这符合文献[21]给出的流场DDT过程的触发机制。

图4 流场参数的分布Fig.4 Distribution of flow parameters

图5 沿x正方向的火焰速度Fig.5 Flame propagation velocity along x positive direction

2.4 燃烧转爆轰过程的影响因素

爆震波阵面是由前沿激波和紧跟其后的化学反应区构成。当前导激波和有限速率化学反应满足强耦合，爆轰波稳定传播。本文中将着重基于前导激波速度和反应区热量释放率2个重要因素，讨论喷射器放置位置和表面阻断率对爆震室燃烧转爆轰DDT的影响机制。

2.4.1基于前导激波速度分析喷射器放置位置的影响

对于喷射器表面阻断率为0.444，数值分析爆震室喷射器放置位置对爆轰流场的影响。表3所示为228 μs时刻，喷射器在不同位置下爆轰波形成情况。可以看到，喷射器位置为3.9 cm时，正好能成功触发爆轰波的形成，因此该处为触发的临界位置。同时，与其余情况相比，喷射器放置位置为4.0 cm时，爆轰波传播距离最大，达到35.74 cm。图6所示喷射器不同放置位置下，到达喷射器处前导激波的速度值。对于喷射器位置为3.9、4.0、5.0 cm情况，前导激波传播速度急剧增加，分别达到649.25、799.7、907.5 m/s 。喷射器位置引起前导激波速度的变化，进而影响流场的热点强度。因此，喷射器位置为3.9～9.0 cm时，燃烧转爆轰DDT过程所需时间较短，形成了稳定传播的爆轰波。

表3 喷射器位置对爆轰波的影响Table 3 Influence of injector positions on detonation

图6 前导激波速度与喷射器位置关系Fig.6 Relation between leading shock velocity and injector position

2.4.2基于波后反应区热量释放率分析喷射器表面阻断率的影响

在距离原点7.0 cm位置处，放置表面阻断率分别为0.400、0.444、0.462和0.471的喷射器，研究爆震室内燃烧转爆轰的过程。

表4给出表面阻断率为0.400、0.444、0.462和0.471喷射器流场的计算结果。得到：表面阻断率从0.444到0.471(增大6.08%)，对应228 μs时间内，爆轰波传播距离增大0.62%。

表4 喷射器表面阻断率对爆轰波形成和发展的影响Table 4 Influence of the blocking rate of injector on detonation

图7所示为74 μs时喷射器不同表面阻断率下流场的温度分布情况。可以看到：火焰经过喷射器后，火焰出现折曲、拉伸等变形现象。即燃料与氧化剂发生强烈的化学反应，燃烧速率急剧增大。热点附近无量纲化的流场的密度、压力和温度分布见图8所示。

图7 通过喷射器流场的火焰分布Fig.7 Distribution of the flame through the injector

结合表4，图7、8可以看出：较其他情况，对于阻断率为0.400的喷射器，前导激波经过后流场压力和温度峰值较低，因此在文中讨论的时间段内未形成稳定的爆轰过程。此外，随着喷射器表面阻断率越大，经过该装置后燃料混合更加充分，因此热点附近燃烧加剧、导致流场压力峰值急剧增大。如，对于表面阻断率为0.471的喷射器，74 μs时刻，压力峰值达到20.26 MPa。

图9所示为化学反应热释放率与时间之间的关系。得到：对于喷射器表面阻断率分别为0.444、0.462和0.471，燃烧热释放率随时间的历程具有相似的变化规律，且热释放率峰值均远大于喷射器阻断率0.400时的情况。如，对于表面阻断率为0.471的喷射器，85.5 μs时刻，燃料单位时间内反应释放的热量达到最大，为爆轰波的稳定传播提供能量。

图8 喷射器不同表面阻断率下热点位置附近流场的参数Fig.8 Flow parameters near the hot spots formed at different blocking rates

图9 反应阵面的热释放率时程曲线Fig.9 Histories of release heat rate at the reaction front

3 结 论

针对爆震室涉及的复杂流动、火焰和激波之间的相互作用，本文中基于二维反应欧拉方程，对具有喷射器的二维爆震室流场开展燃烧转爆轰过程及爆轰性质的数值研究，得到如下结论：

(1)喷射器放置于爆震室中，上游燃料在喷射器附近滞留时间增加，加速与周围流体介质的质量和热量交换，因此喷射器附近流场出现多个热点，为燃料燃烧提供活性反应中心，加快燃烧转爆轰的过程。

(2)基于前导激波速度数值分析喷射器位置对DDT的作用机制发现，喷射器位置为3.9 、4.0 、5.0 cm时，前导激波传播速度急剧增大，分别为649.25、799.7、907.5 m/s，均形成稳定传播的爆轰波。因此在一定范围内，前导激波速度增大，流场的热点较容易触发流场发生燃烧转爆轰现象。

(3)基于波后反应区热量释放率数值分析喷射器表面阻断率对DDT的作用机制发现，随着喷射器表面阻断率的增大，燃料混合充分，反应速率增大，燃烧热释放率也越大。对于表面阻断率为0.471的喷射器，85.5 μs时刻，燃料单位时间内反应释放的热量达到最大，为爆轰波的稳定传播提供能量。