长脉宽爆炸波与圆柱壳作用的波系演化与载荷分布特征

张军; 杨茂; 毛勇建; 穆云飞; 张焕好; 陈志华; 冯晓伟

doi:10.11883/bzycj-2023-0470

长脉宽爆炸波与圆柱壳作用的波系演化与载荷分布特征

doi: 10.11883/bzycj-2023-0470

1.
中国工程物理研究院总体工程研究所，四川绵阳 621999
2.
南京理工大学瞬态物理重点实验室，江苏南京 210094

基金项目: 国家自然科学基金（12072334, 11602256）

详细信息

作者简介:
张　军（1988－　），男，博士，副研究员，aabuzhangjun@foxmail.com

通讯作者:
毛勇建（1976－　），男，博士，研究员，maoyj@caep.cn

中图分类号: O381
计量
- 文章访问数: 376
- HTML全文浏览量: 97
- PDF下载量: 108
- 被引次数: 6
出版历程
- 收稿日期: 2023-12-29
- 修回日期: 2024-02-09
- 网络出版日期: 2024-03-11
- 刊出日期: 2024-08-05

Wave evolution and pressure distribution characteristics of the interaction between long-duration blast load and cylindrical structure

1.
Institute of Systems Engineering, China Academy of Engineering Physics, Mianyang 621999, Sichuan, China
2.
Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China

摘要

摘要: 长脉宽爆炸波与结构作用的波传播及其载荷分布规律是大型爆炸防护设计与安全评估的重要基础。为了掌握长脉宽爆炸波与圆柱壳的相互作用机制及其作用下圆柱壳表面的载荷分布规律，开展了100 ms级冲击波圆柱绕流的激波管试验，并采用大涡模拟方法和高阶WENO-TCD（weighted essentially non-oscillatory-tuned centered difference）混合格式进一步对150 ms长脉宽冲击波与圆柱作用过程的波系演化与压力分布进行了数值分析。结果表明：数值仿真与压力实测结果吻合较好，长脉宽作用下圆柱壳载荷分布呈现出明显的角度和高度相关性，背面压力高于侧面甚至与迎爆面相当，具有不同于传统短脉宽冲击波传播中的压力衰减模式。流场的压力云纹和波系三维结构演化揭示了侧端面的突然扩张是压力初期震荡及相比于正面和背面压力更低的主要原因；系列绕射激波在壳体背面碰撞与反射，以及序列减速激波在135°相位附近的驻定与叠加作用，是引起背面压力呈现出载荷整体提升的主要机制。此外，背风面上尾涡结构的形成与演化过程受边界效应的影响是导致长脉宽圆柱壳载荷分布沿高度方向出现差异的关键因素。
- 长脉宽 /
- 爆炸冲击波 /
- 圆柱壳 /
- 激波管试验
Abstract: The wave propagation and pressure distribution during the interaction between long duration blast waves and structures are important foundations for the large scale explosion protection design and safety assessment. In order to understand the interaction mechanism between long duration blast waves and cylindrical shells, as well as the distribution law of the surface load on the cylindrical shells under their action, the overpressure histories on the cylindrical structure surface were obtained through the 150 ms long duration blast wave shock tube experiment, and the shock wave evolution and the pressure load distribution were investigated numerically using the large eddy simulation and hybrid WENO-TCD (weighted essentially non-oscillatory-tuned centered difference) method. The results show that the overpressure load of the calculation results is in good agreement with the experimental results, and the overpressure load on the cylindrical shell appears a clear angle and height correlation. The pressure on the back shell is higher than that on the side surface or even comparable to the blast on the facing surface, which exhibit different pressure attenuation modes from the traditional short duration blast wave propagation. The sudden expansion on the side surface is the main reason for the initial oscillation of pressure, and has a lower pressure than that at the windward and leeward sides. On the other hand, a series of diffracted shock waves collides and reflects on the symmetry plane of the shell leeward, as well as the stationary and superimposed effects of the series of decelerating shock waves near the 135° phase, which are the main mechanisms that cause the overall increase of the pressure load on the cylindrical shell. In addition, the formation and evolution of wake vortex structures on the leeward side due to the boundary effects is a key factor leading to differences in the load distribution along the height direction. The above analysis methods and related results lay the foundation for the subsequent study of load distribution models for typical structural components under long duration blast waves.
- long-duration /
- blast waves /
- cylindrical shell structures /
- shock waves experiment

HTML全文

爆轰波是一种以高超声速传播的燃烧波。经过前导激波压缩的可燃气体，在非常短的时间和距离内完成自点火过程，通过化学反应释放大量的热，推动爆轰波继续前行，完成爆轰波的自维持传播。因此，爆轰波研究中最本质、最基础的研究是关于爆轰波自维持传播机理的研究^[1-5]。

现有研究多从爆轰波复杂的波系结构以及激波动力学角度出发，定性分析爆轰波传播过程中不同物理因素的作用。V.N.Gamezo等^[6]采用单步反应模型对爆轰波进行数值模拟，研究指出，波后未反应区域的爆炸推动了爆轰波的前进，其形成取决于活化能，而横波的相互碰撞使得中间部分的未反应区域燃烧更快。G.J.Sharpe^[7]研究了爆轰波传播过程中横波的发展和运动情况，使用不同的网格分辨率分析横波的碰撞以及波后未反应袋中燃烧速率对网格的敏感度，研究指出，在分析爆轰波的波系结构时，除高效优化的网格以外，还要处理好反应中的时间、空间尺度之间的匹配关系，尤其在化学反应时间很短的情况，准确捕捉波面结构更为重要。J.Y.Choi等^[8]总结了不同的化学反应模型对爆轰波胞格的影响，讨论了边界条件、计算域长度、时间步长、网格尺度等因素对胞格宽度数值模拟结果的影响，并且定性分析了活化能和反应的指前系数对稳定与不稳定爆轰波传播的影响。

爆轰波的传播涉及到非线性与多尺度的三维复杂流动过程。实验中，通过烟熏法以及瞬态流场捕捉技术可以得到三维爆轰波面传播所留下的鱼鳞状胞格结构。爆轰波传播过程中，前导激波由多个间隔排列的马赫杆、入射激波和横波组成。马赫杆、入射激波和横波相交于三波点，并形成三波结构，三波点的运动轨迹即为胞格结构。爆轰现象包含复杂的激波动力学与化学反应过程，胞格结构是这一过程的外在表现，最能够反映爆轰波传播的重要特征^[9]。近年来，已经利用理论、实验以及数值模拟方法开展了大量关于胞格结构的研究工作。实验研究是通过爆轰波在激波管中的传播，通过不同的测量仪器和手段记录三波点的轨迹以及瞬时的波系结构；数值模拟和理论分析则多是考虑多种因素，如活化能、指前系数、边界条件以及网格尺度等，对胞格大小的定性影响。对于胞格的定量化研究，现有理论均从反应区长度入手，将ZND模型理论计算所得半反应区长度与实验测得的胞格大小做线性拟合，针对特定配比的反应物得到特定的拟合公式。反应区的特征长度能和胞格的大小联系在一起，其中一定存在重要的物理机理。但是，拟合公式却没有普适性，不能说明线性关系存在的原因。

现有理论表明，爆轰波传播过程的主要机制在于非线性波与化学反应带的耦合^[10-12]。燃烧理论中，在高温热化学环境下，自点火机制扮演着重要的角色^[13]。火焰的稳定性主要由化学反应模型控制^[13]。本研究从胞格的大小出发，通过对不同化学反应模型的数值模拟结果的比较，分析自点火机制在爆轰波传播过程中的作用，讨论非线性波与化学反应带的关键耦合参数，从爆轰波传播过程中的耦合机理出发，探讨其中的关键物理机制，指出影响胞格大小的关键因素及其根本原因。

1. 控制方程与模型

爆轰波在二维无限长的直管道中传播，上、下边界满足反射壁面边界条件，计算域见图 1。左侧为高温高压点火区，爆轰波从左向右传播；初始爆轰气体为满足化学当量比的H₂-air混合气体，初始压力p₀=0.1 MPa，初始温度T₀=300 K，当量比ϕ=1.0。为保证爆轰波传播足够长的时间，在数值模拟过程中，保证声速点以及稀疏波完整的前提下，计算区域在计算过程中不断右移，使得计算区域足够大，即爆轰波能够在足够长的管道中传播直至达到稳定状态。由于流场中含有化学反应，流动方程组需要与化学反应动力学模型耦合起来数值求解，计算过程中的“刚性问题”和计算效率问题都会影响到计算所需的计算机内存和CPU时间，本研究在程序中加入并行计算OpenMP(Open Multiprocessing)，在保证网格尺度和网格数目的前提下，使得计算效率大大提高。每个数值算例中的爆轰波的总传播时间都在毫秒量级，以保证胞格尺度统计的可靠性。

图 1 计算模型示意图

Figure 1. Schematic illustration of the problem considered

下载: 全尺寸图片幻灯片

控制方程为二维Euler方程组，忽略了爆轰波传播过程中的扩散效应、粘性效应和热传导过程。化学反应采用一个修正的总包单步反应模型^{[12, 14]}和一个共包含9组元、19个反应的基元反应模型^[15]。式(1)~(7)为单步反应模型控制方程，基元反应模型的详细反应参数参见文献[15]。

$\frac{{\partial \mathit{\boldsymbol{U}}}}{{\partial t}} + \frac{{\partial \mathit{\boldsymbol{F}}}}{{\partial x}} + \frac{{\partial \mathit{\boldsymbol{G}}}}{{\partial y}} = \mathit{\boldsymbol{S}}$

(1)

$\mathit{\boldsymbol{U}} = \left( {\begin{array}{*{20}{c}} \rho \\ {\rho u} \\ {\rho v} \\ {\rho e} \\ {\rho Z} \end{array}} \right), \;\;\mathit{\boldsymbol{F}} = \left( {\begin{array}{*{20}{c}} {\rho u} \\ {\rho {u^2} + p} \\ {\rho uv} \\ {(\rho e + p)u} \\ {\rho uZ} \end{array}} \right), \;\;\mathit{\boldsymbol{G}} = \left( {\begin{array}{*{20}{c}} {\rho v} \\ {\rho uv} \\ {\rho {v^2} + p} \\ {(\rho e + p)v} \\ {\rho vZ} \end{array}} \right), \mathit{\boldsymbol{S}} = \left( {\begin{array}{*{20}{c}} 0 \\ 0 \\ 0 \\ 0 \\ {\dot \omega } \end{array}} \right)$

(2)

$p = (\mathit{\gamma} - 1)\rho [e - ({u^2} + {v^2})/2 - Zq]$

(3)

$T = p/(\rho R)$

(4)

$\dot \omega = - K\rho \;Z\;{\text{exp}}\left( { - \frac{{{E_{\text{a}}}(\mathit{\gamma} (Z) - 1)}}{{R(Z)T({\mathit{\gamma} _{\text{B}}} - 1)}}} \right)$

(5)

$\mathit{\gamma} (Z) = \frac{{\frac{{{\mathit{\gamma} _{\text{U}}}{R_{\text{U}}}Z}}{{{\mathit{\gamma} _{\text{U}}} - 1}} + \frac{{{\mathit{\gamma} _{\text{B}}}{R_{\text{B}}}(1 - Z)}}{{{\mathit{\gamma} _{\text{B}}} - 1}}}}{{\frac{{{R_{\text{U}}}Z}}{{{\mathit{\gamma} _{\text{U}}} - 1}} + \frac{{{R_{\text{B}}}(1 - Z)}}{{{\mathit{\gamma} _{\text{B}}} - 1}}}}$

(6)

$R(Z) = {R_{\text{U}}}Z + {R_{\text{B}}}(1 - Z)$

(7)

式中：p、ρ、u、v、e、γ、q和Z分别表示压力、密度、轴向速度、横向速度、比内能、比热比、混合气体的比放热量和化学反应进行度参数；源项中的 $\dot \omega$ 表示爆轰产物的质量生成率；在化学反应速率的阿伦尼乌斯(Arrhenius)公式中，K、T、E_a和R分别为指前系数、温度、活化能和气体常数；式(5)~(7)中的下标“U”代表反应物，下标“B”代表生成物。单步反应模型的参数取值如表 1所示。

表 1 单步反应模型的参数取值

Table 1. Parameters of one-step reaction model

Z_U	Z_B	γ_U	γ_B	R_U/(J·kg^-1·K^-1)	R_B/(J·kg^-1·K^-1)	E_a/(J·kg^-1)	K/s^-1	q/(J·kg^-1)
1.0	0	1.40	1.24	398.5	368.9	4.794×10⁶	7.5×10⁶	3.5×10⁶

下载: 导出CSV

| 显示表格

在数值模拟中，在空间方向上采用三阶的ENO格式离散，对矢通量进行了Steger-Warming分解。在时间方向上采用三阶TVD Runge-Kutta法。在下面的结果与讨论中，首先研究两个模型的点火延迟时间，然后对二维爆轰波的胞格结构进行模拟。通过对三波点的形成与发展以及自点火过程的分析，给出激波动力学和化学反应耦合的关键机理。

2. 结果与讨论

2.1 点火延迟时间

作为自点火过程中的重要参数，点火延迟时间表征了化学反应的特征时间，是化学反应模型的重要参数。点火延迟时间的定义有多种方法，本文基于温度变化来定义点火延迟时间：对于特定的化学反应模型，利用等容燃烧所得温度T随时间t的变化曲线，定义最大温度变化率(dT/dt)_max所对应的时间为点火延迟时间。对于稳定传播的爆轰波，诱导区内的气体热力学状态对自点火过程更为重要。准确预测与诱导区内气体温度T^*和压力p^*相对应的点火延迟时间对爆轰波的数值模拟是至关重要的。根据一维ZND模型，对于H₂-air爆轰波，当初始压力p₀=0.1 MPa和初始温度T₀=300 K时，对应的诱导区内的气体温度和压力分别为T^*=1544.1 K和p^*=2.77 MPa。但是，在二维爆轰波的数值模拟中，爆轰波面更为复杂，波后诱导区内的气体状态分布是不均匀的，T^*和p^*分别在1 300~1 500 K和1.5~3.0 MPa范围内变化。

图 2给出了总包单步反应模型(model-1)和基元反应模型(model-2)预测的初始压力为0.1 MPa下的点火延迟时间随温度的变化情况，同时还给出了部分实验结果作为对比。可以看出，在对数坐标下，在1 000~1 500 K(对应的横坐标取值范围为10.0~6.6 K^-1)的高温下，爆轰模型预测的点火延迟时间随温度近似呈线性变化，实验结果也基本反映了该规律。点火延迟时间与温度成反比，温度越高，点火延迟时间越短。相同温度下，两个模型计算得到的点火延迟时间并不相同。两个模型预测的点火延迟时间曲线斜率不同，在T=1 200 K (横坐标取值为8.3)左右相交。在1 300~1 500 K的温度范围内(对应的横坐标取值范围为7.7~6.6 K^-1)，同一温度下，利用model-2计算得到的点火延迟时间比model-1要长，并且温度越高，差别越大，但是两者的差值在18 μs以内。曲线的斜率表示点火延迟时间对温度的敏感程度，斜率越大，意味着点火延迟时间对温度变化越敏感。需要说明的是，model-1的计算结果是与压力无关的，即改变初始压力，对利用model-1计算得到的点火延迟时间没有影响。这是单步反应模型的重要特性，也是利用这类模型进行爆轰波数值研究时需要注意的事项。

图 2 模型和实验预测的点火延迟时间对比

Figure 2. Comparison of ignition delay times predicted by two models and the experimental data

下载: 全尺寸图片幻灯片

与单步反应模型不同的是，基元反应模型的点火延迟时间是与压力相关的。图 3给出了不同压力下利用Model-2计算得到的点火延迟时间，并将其与model-1的预测结果作比较。首先，比较在初始压力为0.1、1.0和2.0 MPa时model-2预测的点火延迟时间可知：在1 200~1 500 K的高温段(对应的横坐标取值范围为8.3~6.6)，随着压力的增加，相同温度下的点火延迟时间逐渐减小，并且与压力基本呈线性反比变化，这是基元反应模型的重要特性。以T=1 300 K(横坐标取值为7.7)为例，model-1预测的点火延迟时间(τ_ig)为19.82 μs，且该数值不随压力变化。但是，对于model-2，当压力为0.1 MPa时，τ_ig=30.90 μs；当压力为1.0 MPa时，τ_ig=3.00 μs；当压力为2.0 MPa时，τ_ig=1.00 μs。压力越高，model-2预测的点火延迟时间就越短，相比于model-1预测结果，两者的差值可以达到一个数量级甚至更大。高温高压下的实验结果比较缺乏，因此无法与实验结果进行比较。

图 3 不同压力下model-1和model-2的点火延迟时间

Figure 3. Comparison of ignition delay times predicted by model-1 and model-2 at different pressures

下载: 全尺寸图片幻灯片

对于初始压力0.1 MPa、初始温度300 K、满足化学当量比的H₂-air爆轰波，诱导区内温度和压力分别在1 300~1 500 K和1.5~3.0 MPa范围内，研究这一范围内的点火延迟时间更为重要。图 4比较了初始压力为1.0 MPa时不同模型预测的点火延迟时间(除model-1和model-2的计算结果外，其他结果参见文献[16-17])。可以看出，对于基元反应模型来说，点火延迟时间在不同的温度区间对温度的敏感性是不同的，不同的反应模型在相同温度下的点火延迟时间都不相同，其差异甚至在一个数量级以上。在1 300~1 500 K的范围内，model-1的点火延迟时间比其他基元反应模型的结果都要长。

图 4 不同化学反应模型预测的点火延迟时间

Figure 4. Comparison of ignition delay times predicted by different chemical reaction models

下载: 全尺寸图片幻灯片

2.2 三波点的运动

利用model-1和model-2两个爆轰模型分别模拟了H₂-air混合气体中稳定传播的爆轰波及其胞格结构。为使数值模拟结果具有可靠性，在不同网格分辨率和计算域大小的条件下开展了研究。网格尺度满足Δx=Δy，其取值为10或20 μm；计算域宽度(w)取值为2、3、4、8或16 mm。计算区域足够长，以避免出口边界对爆轰波面的影响。爆轰波的数值模拟中，首先采用直接点火的方式，形成过驱爆轰波，过驱爆轰波逐渐衰减为稳定传播的CJ爆轰波。当爆轰波传播足够长的时间，三波点个数不再发生改变时，爆轰波即达到稳定传播的状态。本文对爆轰波传播过程的分析均针对达到稳定传播状态的爆轰波。

图 5给出了利用model-1模拟所得稳定传播的爆轰波的流场压力分布，观察图 5可得爆轰波波面三波点的个数：当w=4 mm时，波面上存在一个三波点；当w=8 mm时，有一对三波点；当w增加到16 mm时，有两对三波点。3种情况下，平均每8 mm的计算宽度，波面有一对三波点的存在，说明计算结果与网格数以及网格分辨率均无关，计算结果可靠。

图 5 采用model-1数值模拟二维爆轰波传播，流场的压力等值线

Figure 5. Instantaneous contours of pressure for a two-dimensional detonation propagation simulated by model-1

下载: 全尺寸图片幻灯片

一对三波点从一次碰撞开始到下一次碰撞为止，完成一个运动周期。经过两个运动周期的一对三波点，其运动轨迹形成一个完整的胞格。对于稳定传播的爆轰波，三波点运动周期是固定不变的。图 6给出利用model-1模拟所得稳定传播的爆轰波三波点的运动情况，以计算域为w=8 mm为例说明了一对三波点在一个周期内的运动规律，其中每两幅图之间的时间间隔是2 μs。从图 6中可以看出，一对三波点完成一个完整胞格需耗时6 μs，因此这一对三波点的运动周期则约为3 μs。在诱导区内的气体的平均压力为1.9 MPa，平均温度为1 500 K左右。在该状态下，单步反应模型model-1预测的点火延迟时间为3.46 μs(见图 2)。由此可见，三波点的运动周期与诱导区内气体对应的点火延迟时间一致，其他算例也给出了同样的结论。

图 6 采用model-1数值模拟的一对三波点的压力等值线分布

Figure 6. Movement of triple-wave points shown with the pressure contour maps in the numerical simulation by model-1

下载: 全尺寸图片幻灯片

图 7给出了利用model-2模拟所得稳定传播的爆轰波的流场压力分布图。当w=2 mm时，波面上存在2对三波点；当w=3 mm时，为3对三波点；当w增加到4 mm时；波面上则有4对三波点。3种情况下，平均每1 mm的计算宽度，波面有一对三波点的存在，保证了计算结果的网格无关性。

图 7 采用model-2数值模拟二维爆轰波传播，流场的压力等值线

Figure 7. Instantaneous contours of pressure for a two-dimensional detonation propagation simulated by model-2

下载: 全尺寸图片幻灯片

图 8为利用model-2模拟所得稳定传播的爆轰波三波点的运动情况。在计算域为w=3 mm的情况下，得到了3对三波点，其中每两幅图之间的时间间隔是0.1 μs。从图 8中可以看出，一对三波点运动完成半个胞格用时0.9 μs，即这一对三波点的运动周期约为0.9 μs。而在1.9 MPa、1 500 K时，基元反应模型model-2预测的点火延迟时间为0.94 μs(见图 2)。由此可见，三波点的运动周期与诱导区内气体对应的点火延迟时间一致。其他算例也给出了同样的结论。

图 8 采用model-2数值模拟的一对三波点运动的胞格结构

Figure 8. Movement of triple-wave points shown with the cellular structure simulated by model-2

下载: 全尺寸图片幻灯片

2.3 胞格结构

对于稳定传播的爆轰波的胞格结构，胞格尺度是指所有胞格的平均宽度。图 9为利用model-1模拟得到的稳定的胞格结构。当w=4 mm时，垂直于爆轰波的传播方向上形成了半个胞格；当w=8 mm时，形成1个胞格；当w=16 mm时，形成2个胞格。所有3种计算条件下，平均胞格尺度均为8 mm，即计算结果的网格无关。实验测得的胞格尺度在8~15 mm之间^[18]，因此利用model-1模拟得到的胞格尺度与实验值完全一致。图 10给出了利用model-2模拟得到的稳定胞格结构。从图 10中可以看出：在w=2 mm的计算域中，垂直于爆轰波传播方向得到2个完整的胞格；当w=3 mm时，得到3个胞格；当w=4 mm时，胞格数增加到4个。胞格宽度同样与网格无关，利用model-2数值模拟得到的平均胞格宽度为1 mm，是model-1和实验所得胞格宽度的1/8。

图 9 采用model-1数值模拟二维H₂-air爆轰波的胞格结构

Figure 9. Cellular structures for a two-dimensional H₂-air detonation simulated by model-1

下载: 全尺寸图片幻灯片

图 10 采用model-2数值模拟二维H₂-air爆轰波的胞格结构

Figure 10. Cellular structures for a two-dimensional H₂-air detonation simulated by model-2

下载: 全尺寸图片幻灯片

通过以上分析可知，在诱导区温度为1 300~1 500 K、压力为1.5~3.0 MPa的范围内，model-1给出的点火延迟时间比model-2的结果要长，并且采用model-1计算得到的胞格比model-2的胞格更宽，说明点火延迟时间和胞格宽度成正比。另外，利用model-1所得点火延迟时间的斜率比model-2要大(见图 2)，其对应的自点火过程对温度变化就更为敏感；在相同的温度和压力变化下，斜率大的点火延迟时间变化范围更大，波面上各点的点火延迟时间相差也比较大，自点火过程相对更加不稳定，胞格更不规则。

现有模型预测的胞格尺度都比实验值至少小两倍^[19]，主要原因在于高温高压下的点火机制上存在问题。现在化学反应模型的参数设置是根据层流火焰速度与点火延迟时间的实验数据来确定的，而实验中不仅缺乏低温的结果，在高温高压的实验结果也是微乎其微，使得模型在预测高温高压下的点火机制时，有某些重要的物理机制被忽略或表达有误。当然，研究中使用的单步反应模型之所以能够准确的模拟胞格结构，主要在于该模型能够准确模拟当量比的H₂-air混合气体在1.9 MPa和1 500 K下的点火延迟时间。实际上，model-1的点火延迟时间与压力无关，其胞格宽度也与压力无关，无论在什么初始压力条件下，模拟得到的胞格大小都是不变的，这一点与物理规律不符，需要引起注意。

上述研究结果表明，爆轰波波后诱导区的气体的压力、温度对应的点火延迟时间与三波点的运动周期一致。运动周期反映的是爆轰波的动力学特性，而点火延迟时间反映的是爆轰波的化学反应动力学特性，二者是通过同一个时间尺度耦合在一起的。从爆轰波传播的物理机制来看，首先，前导激波压缩混合气体，使其压力和温度升高达到可燃条件，经过该状态对应的点火延迟时间后，这部分气体开始燃烧，完成自点火过程。自点火中释放的大量的热量，推动爆轰波继续向前，完成激波和燃烧的耦合。这两个特征时间尺度的匹配就尤为重要。如果点火延迟时间过长，那么化学反应的能量释放太慢，就无法维持爆轰波的自维持传播；如果点火延迟时间过短，又会导致在三波点到达之前，不应该燃烧的气体提前燃烧，产生新的三波点，那么得到的胞格数量就会变多，胞格尺寸就会变小。只有二者相互匹配，才能够得到稳定的胞格结构，同时保证激波和化学反应带的耦合。前导激波压缩过的爆轰气体可以实现自点火，自点火过程释放的能量得以推动爆轰波继续传播，这就是气相爆轰波流动与化学反应的耦合机理。对于化学反应模型来说，为了准确地模拟爆轰波的传播过程，就要保证模型能够准确模拟关键状态下的点火延迟时间，再现自点火过程，并确保化学反应动力学和气体动力学的特征时间尺度完全匹配。

3. 结论

基于二维欧拉方程数值研究的爆轰波，采用基元反应模型和单步反应模型描述化学反应过程，对不同模型的点火延迟时间进行比较，同时统计了爆轰波在稳定传播的状态下，两个模型得到的胞格的宽度以及三波点的运动周期，研究结果表明：(1)爆轰波的胞格宽度与点火延迟时间成正比；(2)化学反应带中诱导区内气体对应的点火延迟时间等于三波点的运动周期，点火延迟时间是非线性波的特征时间，三波点的运动周期是化学反应带的特征时间，爆轰波传播过程中，非线性波和化学反应带通过特征时间这一关键参数耦合在一起，爆轰波依靠爆轰气体的自点火过程得以自维持传播。

图 1 激波管试验系统

Figure 1. Shock tube experiments system

下载: 全尺寸图片幻灯片

图 2 圆柱壳与测点布局

Figure 2. Cylindrical shell and measuring points layout

下载: 全尺寸图片幻灯片

图 3 计算模型

Figure 3. Schematic of the computational model

下载: 全尺寸图片幻灯片

图 4 流场计算域网格分布

Figure 4. Grid distribution in computational domain

下载: 全尺寸图片幻灯片

图 5 不同网格尺度下圆柱壳体表面压力载荷分布曲线

Figure 5. Grid scales effect on the surface pressure load distribution curves of cylindrical shell

下载: 全尺寸图片幻灯片

图 6 H = 0.9 m时，表面压力监测点计算结果与测试结果的对比

Figure 6. Comparison of the overpressure with the corresponding experimental results at H = 0.9 m

下载: 全尺寸图片幻灯片

图 7 H = 0 m时圆柱壳体表面压力载荷变化曲线

Figure 7. Overpressure histories on the surface of cylindrical structure at H = 0 m

下载: 全尺寸图片幻灯片

图 8 H = 0 m时压力载荷峰值沿角度的分布曲线

Figure 8. Distribution of the peak overpressure along phase at H = 0 m

下载: 全尺寸图片幻灯片

图 9 长脉宽平面冲击波与圆柱壳体作用过程的密度纹影图（H = 0 m）

Figure 9. Numerical schlieren image of the interaction between long-duration blast and cylindrical structure at H = 0 m

下载: 全尺寸图片幻灯片

图 10 单脉冲激波与圆柱壳体作用过程的激波结构^[8]

Figure 10. Shock structure of the interaction between a pulse shock wave and a cylindrical structure^[8]

下载: 全尺寸图片幻灯片

图 11 H = 0 m时入射激波I₁～I₄与壳体作用过程中壳体表面4个测点压力载荷变化曲线

Figure 11. Overpressure on the surface of cylindrical structure at H = 0 m

下载: 全尺寸图片幻灯片

图 12 长脉宽平面冲击波与圆柱壳体作用过程中壳体表面压力分布图（左图为云纹图；中间为迎爆面压力分布；右图为背部压力分布）

Figure 12. Pressure distribution on the surface of a cylindrical structure (The left image is a moiré pattern; the middle is the pressure distribution on the blasting surface; the right image is the pressure distribution on the back)

下载: 全尺寸图片幻灯片

图 13 不同方位处圆柱壳体表面压力载荷变化曲线

Figure 13. Overpressure load on the surface of cylindrical structure at different orientations

下载: 全尺寸图片幻灯片

图 14 测点压力载荷峰值沿角度和高度的分布规律

Figure 14. Distribution of the peak overpressure along the angle and height

下载: 全尺寸图片幻灯片

图 15 圆柱壳体下游流场的三维尾涡结构

Figure 15. Three-dimensional wake structures at the leeward of a cylindrical structure

下载: 全尺寸图片幻灯片

参考文献(13)

[1]	吴昊, 方秦, 龚自明, 等. 冲击爆炸作用对核电站安全壳毁伤效应研究的进展 [J]. 防灾减灾工程学报, 2012, 32(3): 384–392. DOI: 10.3969/j.issn.1672-2132.2012.03.020. WU H, FANG Q, GONG Z M, et al. State of arts of impact and blast effects on the NPPC [J]. Journal of Disaster Prevention and Mitigation Engineering, 2012, 32(3): 384–392. DOI: 10.3969/j.issn.1672-2132.2012.03.020.
[2]	CLUBLEY S K. Long duration blast loading of cylindrical shell structures with variable fill level [J]. Thin-Walled Structures, 2014, 85: 234–249. DOI: 10.1016/j.tws.2014.08.021.
[3]	CLUBLEY S K. Non-linear long duration blast loading of cylindrical shell structures [J]. Engineering Structures, 2014, 59: 113–126. DOI: 10.1016/j.engstruct.2013.10.030.
[4]	沙莎, 陈志华, 姜孝海. 激波与障碍物作用加速衰减的数值研究 [J]. 工程力学, 2014, 31(9): 239–244. DOI: 10.6052/j.issn.1000-4750.2012.11.0914. SHA S, CHEN Z H, JIANG X H. Investigations into the accelerating attenuation of shock waves by obstacles [J]. Engineering Mechanics, 2014, 31(9): 239–244. DOI: 10.6052/j.issn.1000-4750.2012.11.0914.
[5]	高康华, 李斌, 刘宇都, 等. 外爆条件下冲击波与组合壳结构的相互作用 [J]. 爆炸与冲击, 2019, 39(5): 052101. DOI: 10.11883/bzycj-2017-0370. GAO K H, LI B, LIU Y D, et al. Interaction of a shock wave with a composite shell structure under an external explosion [J]. Explosion and Shock Waves, 2019, 39(5): 052101. DOI: 10.11883/bzycj-2017-0370.
[6]	杨科之, 刘盛, 周布奎. 光滑坑道中基于波形参数的冲击波衰减规律 [J]. 防护工程, 2022, 44(4): 30–35. DOI: 10.3969/j.issn.1674-1854.2022.04.005. YANG K Z, LIU S, ZHOU B K. The law of shock wave attenuation based on waveform parameters in smooth tunnels [J]. Protective Engineering, 2022, 44(4): 30–35. DOI: 10.3969/j.issn.1674-1854.2022.04.005.
[7]	GAUTIER A, SOCHET I, COURTIAUD S, et al. Evaluation of shock wave properties interacting with a parallelepipedal obstacle [J]. International Journal of Protective Structures, 2022, 13(4): 716–731. DOI: 10.1177/20414196221114908.
[8]	GAUTIER A, SOCHET I, COURTIAUD S. Analysis of shock wave interaction with an obstacle by coupling pressure measurements and visualization [J]. Sensors, 2022, 22(9): 3325. DOI: 10.3390/s22093325.
[9]	ZHANG H H, CHEN Z H, GUO Z Q, et al. Characteristic behavior of shock pattern and primary vortex loop of a supersonic square jet [J]. International Journal of Heat and Mass Transfer, 2017, 115: 347–363. DOI: 10.1016/j.ijheatmasstransfer.2017.07.002.
[10]	ZHANG H H, AUBRY N, CHEN Z H, et al. The evolution of the initial flow structures of a highly under-expanded circular jet [J]. Journal of Fluid Mechanics, 2019, 871: 305–331. DOI: 10.1017/jfm.2019.285.
[11]	PANTANO C, DEITERDING R, HILL D J, et al. A low numerical dissipation patch-based adaptive mesh refinement method for large-eddy simulation of compressible flows [J]. Journal of Computational Physics, 2007, 221(1): 63–87. DOI: 10.1016/j.jcp.2006.06.011.
[12]	张焕好, 肖焱丰, 郑纯, 等. 超声速横向射流三维流场结构特征 [J]. 推进技术, 2023, 44(1): 210730. DOI: 10.13675/j.cnki.tjjs.210730. ZHANG H H, XIAO Y F, ZHENG C, et al. Three-dimensional flow characterizes of transverse jet injection into supersonic crossflow [J]. Journal of Propulsion Technology, 2023, 44(1): 210730. DOI: 10.13675/j.cnki.tjjs.210730.
[13]	VITI V, NEEL R, SCHETZ J A. Detailed flow physics of the supersonic jet interaction flow field [J]. Physics of Fluids, 2009, 21(4): 046101. DOI: 10.1063/1.3112736.

施引文献

期刊类型引用(4)

1.	赵焕娟，汪乐溪，张英华，林敏，周正青，蒋慧灵. 孔隙比对氢气爆轰波结构的定量影响. 工程热物理学报. 2024(01): 303-309 . 百度学术
2.	刘思远，汪洋，李真珍，滕宏辉，田诚. 中心锥喷管喉道比参数对旋转爆轰燃烧影响的数值研究. 推进技术. 2023(09): 120-131 . 百度学术
3.	王硕，王中，余杰皓. 喷射气流灭弧防雷器的爆轰弹丸设计研究. 电瓷避雷器. 2020(02): 20-25 . 百度学术
4.	武郁文，褚驰，翁春生，郑权. 孔板扰动对爆轰波胞格结构特性影响的实验研究. 爆炸与冲击. 2019(11): 15-23 . 本站查看