Chinese Journal of Ship Research
后燃烧效应对约束空间内爆炸载荷的影响规律
1，2，吴卫国1，2徐维铮1 430063武汉理工大学高性能舰船技术教育部重点实验室，湖北武汉2 430063武汉理工大学交通学院，湖北 武汉
徐维铮，吴卫国
摘 要：［目的］为了研究后燃烧效应对约束空间内炸药爆炸载荷的影响规律，［方法］将反应率模型耦合到可WENO有限差分格式和三阶TVD-RK格压缩欧拉方程中，并以源项的方式进行后燃烧能量的添加。采用五阶式对耦合方程空间项和时间项进行离散求解，自主编写约束空间内炸药爆炸后燃烧过程的二维数值计算程序。基于自主程序开展不同工况下内爆炸后燃烧效应数值计算，探讨内爆炸过程中反应速率及后燃烧能量大小对爆炸载荷的影响规律。［结果］研究表明：在后燃烧能量大小一定的情况下，反应速率常数增大时，冲击波到达时间提前，冲击波峰值、冲量均增大，准静态超压峰值保持不变；在反应速率常数一定的情况下，随着后燃烧能量的增大，冲击波峰值、冲量及准静态超压峰值均增大；后燃烧能量的加入能显著增强爆炸载荷强度。［结论］研究成果可为抗爆结构设计及毁伤评估提供一定的参考和指导。WENO关键词：约束空间内爆炸；后燃烧效应；反应率模型；五阶 格式；爆炸载荷
Afterburning effect on blast load in confined space
Xu Weizheng1，2，Wu Weiguo1，2 1 Key Laboratory of High Performance Ship Technology of Ministry of Education，Wuhan University of Technology， Wuhan 430063，China 2 School of Transportation，Wuhan University of Technology，Wuhan 430063，China Abstract：［Objectives］ In order to study the afterburning effect on blast load in confined space， ［Methods］the reaction rate model is coupled to the compressible Euler equations，and the afterburning energy is added in the form of source term. The space term and time term of the coupled equations are solved discretely using the fifth-order WENO finite difference scheme and the third-order TVD-RK scheme. Then an in-house two-dimensional numerical calculation code is developed for afterburning effect in confined space，on the base of which the numerical calculations of the afterburning effects in different explosion cases in confined space are carried out，and the influences of reaction rate and afterburning energy on the blast load are discussed. ［Results］ The study shows that，1） in the case of constant afterburning energy，when the reaction rate increases，the shock wave arrives earlier，the shock wave peak and impulse increase，and the final quasi-static overpressure peak remains unchanged；2）in the case of constant reaction rate，when the afterburning energy increases，the shock wave peak，impulse and final quasi-static overpressure peak increase. The addition of the afterburning energy can significantly enhance the strength of the blast load.［Conclusions］The study in this paper can provide reference and guidance for the design and damage assessment of anti-explosion structures. Key words：explosion in confined space ；afterburning effect；reaction rate model；fifth-order WENO scheme；blast load
0引言
针对约束空间内的爆炸问题，爆炸后燃烧效应是需要重点关注的内爆炸物理现象。大多数炸药（TNT炸药、温压炸药、SDF 药［1］等）的混合型炸爆轰产物具有负氧性，高温高压爆轰产物在膨胀过程中会与周围空气中的氧气进行剧烈的燃烧反应并释放大量能量，这一物理现象为后燃烧效应［2］。后燃烧效应对爆炸冲击波的传播过程、冲击波壁面反射压力、准静态压力都会产生一定的影响。为了研究约束空间内后燃烧效应对爆炸载荷的影响规律，国内外学者开展了大量的实验和数等［3］采值模拟研究。金朋刚 用压力传感器测量了3在密闭罐体内分别存在氮气、空气和氧气 种条TNT件下 炸药在后燃烧过程中产生的准静态压力，结果表明，在氧气环境下，TNT炸药的后燃烧效应较明显，能够产生更大的准静态压力。李鸿500 L宾等［ 4 ］在容积为 的密闭爆炸罐中进行了TNT炸药爆炸实验，发现随着环境中氧气量的增大，准静态压力增大，说明提高环境中的氧气含量能够提高爆轰产物的反应率。李芝绒等［5］通过实验测量了在圆柱形密闭爆炸罐内分别存在空气和2氮气这 种条件下温压炸药爆炸的冲击波峰值以及罐体内的准静态压力，结果表明，在爆轰反应阶段，氧气与超细铝粉发生氧化反应，空气环境中的冲击波峰值和冲量比氮气环境的略高；在后燃烧阶段，氧气与铝粉混合产生后燃烧反应，释放出大量的热量，与氮气环境相比，空气环境中的准静态压力和热响应温度峰值显著增大。Kuhl等［6 ］提出采用多流体模型描述考虑后燃烧效应的爆炸场， 3将爆炸组分划分为燃料、空气和产物 种成分，分3别对这 种成分建立质量、动量、能量守恒方程，在爆炸过程中，采用狄拉克函数进行燃烧阵面的捕捉，释放后燃烧能量，并开发了网格自适应数值Togashi等［7-8］将自主开发的三维爆炸计算程序。波数值计算程序与劳伦斯利弗莫尔实验室开发的Cheetah Cheetah编码衔接，利用 编码计算后燃烧能量，再将后燃烧能量添加到自主开发的程序，实现爆炸过程中后燃烧效应的模拟。爆炸流场包含高密度比和高压力比强间断等复杂流场结构，对爆炸过程进行数值模拟时，需要高精度、强稳定性激波捕捉格式。Liu等［9］提出了Weighted Essentially NonOscillation加权本质无振荡格式（ Scheme，WENO），Jiang 等［10-11 ］发展了
该格式并扩展了其应用。目前，WENO格式作为一种典型的高精度激波捕捉格式，对流场内的激
波间断具有较高的分辨率，适用于求解包含激波、膨胀波以及接触间断等复杂结构的流场。由于后燃烧过程涉及复杂的多组分燃烧化学反应，如果考虑详细的化学反应过程，不仅程序编写复杂，且由于化学反应时间尺度与流场时间尺度存在差别，计算量大，难以应用于实际工程问题
计算中。为此，本文拟提出一种考虑后燃烧效应WENO的简化反应率模型，考虑到 格式精度高及Fortran稳定性较好的优势，基于 平台，采用五阶WENO有限差分格式，自主开发约束空间内考虑爆炸后燃烧效应的二维数值计算程序，并探讨反
应速率及后燃烧能量大小对约束空间内爆炸载荷的影响规律。
1 控制方程及数值求解
Miller为了近似考虑爆炸后燃烧效应，在 反应率模型［12］思想的启发下，提出一种考虑炸药爆
炸后燃烧效应的数值计算方法，初步探讨后燃烧
效应对约束空间内爆炸载荷的影响规律。Miller反应率模型最初提出的目的是近似描述含铝炸药爆轰过程中铝粒子在爆轰波阵面后的反应过程［12］。我们尝试将模型思想推广应用到约束空间内爆炸后燃烧效应的数值计算中。本文采
用的反应率模型如下：
1 1 dα =- aα p6 （ ） 1 2 dt式中：a 为反应速率常数；p为流体压力；α 为后燃烧过程中反应率（初始时 α = 1 ，反应完成后α = 0 ）。 不考虑后燃烧效应的可压缩欧拉方程为¶U ¶F ¶G （2） + + =0 ¶t ¶x ¶y
其中，
ρ ρu ρv ρu ρu2 +p ρvu 3 U= F= G= （ ） 2 ρv ρuv ρv +p E u ( E + p) v ( E + p)将反应率模型式（1）耦合到式（2）中，并以源项的形式进行后燃烧能量的添加，可得考虑后燃
烧效应的可压缩欧拉方程为
¶U ¶F ¶G =S （4）
+ + ¶t ¶x ¶y 其中，
ρu ρ ρu ρu2 +p U= ρv F= ρuv E u ( E + p) ρα ρuα
ρv ρvu （5） G= ρv 2 +p S= 1 1 ρaα 2 p 6 Q v ( E + p) af 1 1 ρvα -ρaα 2 p6
6 E = ρe + 1 ρu2 + 1 ρv2 （ ） 2 2 （7） p =( γ - 1) ρe以上式中：ρ为密度；u v 分别为 x y 方向上的速度分量； E 为单位体积流体的总能量；e 为比内能；Q 为爆炸后燃烧过程中单位质量释放的能af量，J/kg；γ 1.4。为气体的绝热指数，文中γ统一取为式（5）在每个方向上均可以看成是一个带有源项的双曲守恒律方程： ¶u ¶f ( u) 8 + = s ( u) （ ） ¶t ¶x方向，式（8）的半离散守恒型格式为例如，针对x dui =-¶f | 1 ( )+ + sx »- fˉ - fˉ sx dt ¶x Dx x = x = x i + 12 i - 12 =x i i i9
（ ）式中： fˉ 和 fˉ 分别为单元（ xi ， xi ） i - 12 i + 12 - 1/2 + 1/2左、右界面处的对流项数值通量； D x 为 x 方向的均匀网格间距， x和 y方向的网格间距相同。采WENO格式对式（9）中的对流项进行数值用五阶 TVD TVD-RK）对重构，采用三阶 龙格库塔法（式（9）中的时间项进行数值离散［13］。
2 约束空间内爆炸后燃烧效应数值计算
WENO采用五阶 有限差分格式，自主开发了约束空间内考虑爆炸后燃烧效应的二维数值计算程序。自主程序对于爆炸波的数值计算可靠性在文献［14-15］中已经得到验证。本节采用自主程序对约束空间内炸药爆炸过程进行数值模拟研究，主要探讨反应率演化过程及后燃烧能量的加入对爆炸载荷的影响规律。
2.1 舱室尺寸及测点布置
1（a长方形舱室的尺寸如图 ）所示，图中mm。在舱室壁面上设置了3数值单位为 个典型测点，用于监测爆炸载荷时间历程和反应率时间10 mm，如图1（b）所示。历程。均匀网格步长取为
2.2 初始条件和边界条件
基于瞬时爆轰假定，将方形炸药等效为均22.2 mm匀高压气团，具体参数为：边长 ，密度1 630 kg/m3，压力 3.057 9×109 Pa。炸药设置在舱2（a）中的红色区域，周围蓝色区域为室中间，即图 1.0 kg/m3，压 1.0×105 Pa。空气域，空气密度为 力为2（b初始时刻反应率分布见图 ）中的红色区域。初始时刻反应还没发生，炸药所在区域反应率1，舱室其他区域反应率为0。壁面边界条件设为置为刚性边界，由于冲击波与结构变形耦合效应十分复杂，这里暂时不考虑舱室结构的变形。
2.3 计算工况
为了探讨反应速率和后燃烧能量大小对爆炸
2过程的影响规律，选取 类典型工况进行计算。工1：后燃烧单位质量释放的能量 =4.69×106 J/kg，况 Q af =0，保持其他参数不变，反应速率常数分别取为a 10，40。工况2：反应速率常数为 =10，保持其他a参数不变，后燃烧单位质量释放的能量分别取为=0，3.0×106，4.69×106 J/kg。Q af
2.4 反应率演化过程分析
为了探讨爆炸后反应率在舱室内部的演化过2程，给出了 种工况在不同时刻的反应率分布图， 3 =4.69×106 J/kg；如图 所示，其中左图为 a = 0 Q af =4.69×106 J/kg。由图 3可知，2右图为 a = 10 Q af种工况下，反应率在不同时刻总体呈现相似的空=4.69×106 J/kg间分布形态。在爆炸工况a = 0 Q af中反应速率常数为零，即在整个爆炸过程中没有发生反应，没有后燃烧能量的加入。根据式（4）可知，该工况求解的是爆轰产物质量分数的演化过=4.69×106 J/kg程；由于爆炸工况a = 10 Q 中反应af速率常数不为零，反应速度快，从而导致舱室内部0。的反应率快速下降，最终趋近于 4给为定量显示舱室内部反应率的变化，图1（ =0，10，40；Q =4.69×106 J/kg）壁面出了工况 a af
1 2 4测点 和 的反应率时间历程曲线。由图 可见，随着反应速率常数 a 的增大，反应率迅速下降，反0。这里需要说明的是，当应完成后，反应率等于时，反应没有发生，根据式（4）反应速率常数 a = 0可知，反应率的演化过程就是爆轰产物质量分数
的演化过程，由于爆炸后爆轰产物与舱室内部的
空气进行了复杂的掺混过程，因此舱室内部爆轰0 1。产物的质量分数数值将趋近于大于 而小于 图5给出了工况2（Q =0，3.0×106，4.69×106 J/kg； af =10 1 2 a ）壁面测点 和 的反应率时间历程曲线。5由图 可以看出，在反应速率常数一定的情况下，后燃烧能量的加入对反应速率时间历程的影响显
著，然而随着后燃烧能量 Q 的增大，其大小对反af应速率时间历程的影响较小。
2.5 后燃烧能量的加入对爆炸载荷的影响
为了探讨不同反应速率常数a对舱室内部爆1炸 载 荷 的 影 响 规 律 ，给 出 了 工 况J/kg）壁面典型测点3 （a = 0 10 40; Q = 4.69 ´ 106 af 6 6的超压时间历程曲线（图 ）。由图 可以看出，后燃烧能量的加入明显增强了冲击波载荷和冲量；
随着反应速率常数 a 的增大，冲击波到达时间提
前，冲击波峰值增大，冲量增大。由于后燃烧能量=0，Q =4.69×106 J/kg的加入，与工况 a 相比，工af =10，40；Q =4.69×106 J/kg况 a 的准静态超压峰af =10，值增大。由于后燃烧能量相同，工况 a =4.69×106 J/kg =40，Q =4.69×106 J/kg Q 和工况 a af af具有相同的准静态超压峰值。
为了探讨不同后燃烧能量大小Q 对舱室内af 2（ =0，部爆炸载荷的影响规律，给出了工况 Q af 3.0×106，4.69×106；a =10）壁面典型测点3的超压7）。由图7时间历程曲线（图 可以明显看出，在反应速率常数一定的情况下，随着后燃烧能量 Q af的增大，载荷强度增大，最终的准静态超压峰值增大。为了初步验证本文后燃烧能量加入的可靠性，从准静态超压峰值的角度进行了对比分析。献［16］可知，封闭舱室内部准静态超压峰根据文值的计算公式为 10 Dp =( γ - m Q - （ ） 1)V -p0 ( γ 1) s tol ρE式中：m 为炸药质量；Q = Q + Q ，为爆炸过tol TNT af程中释放的总能量，包含后燃烧过程中释放的能=4.69×106 J/kg，为炸药的爆热；量 Q ，其中 Q af TNT =1 630 kg/m3， p 为初始大气压力；V 为舱室体积；ρ 0 E为炸药密度。1，2 10将工况 中的后燃烧能量Q 代入式（ ）， af得到准静态超压峰值的理论计算值，再将理论计6（a）算结果与数值模拟结果进行对比，结果如图7（a）所示。由图可以明显看出，理论计算结和图果与数值模拟结果吻合较好，证明了后燃烧能
量加入的可靠性。
3结论
通过研究得到如下主要结论： 1）在后燃烧能量大小一定的情况下，反应速率常数增大时，冲击波到达时间提前，冲击波峰
值、冲量均增大，准静态超压峰值保持不变。2）在反应速率常数一定的情况下，随着后燃烧能量的增大，冲击波峰值、冲量及准静态超压峰
值均增大，后燃烧能量的加入能显著增强爆炸载
荷强度。3）本文针对爆炸后燃烧过程的数值计算，尽管没有考虑复杂的多组分反应过程，然而采用一
种简化的反应率模型近似描述后燃烧能量对爆炸载荷的影响，不失为一种满足实际工程计算的有效方法。本文的研究方法及结果可为进一步研究内爆炸复杂多组分后燃烧效应及抗爆结构设计提供参考和借鉴。
参考文献：
［1］ Kuhl A L，Bell J B，Beckner V E． Turbulent combus⁃ tion in SDF explosions ［C］//Proceedings of the JAN⁃ NAF Conference. San Diego：Lawrence Livermore Na⁃ tional Lab，2009. ［2］ Kuhl A L，Reichenbach H． Combustion effects in con⁃ fined explosions［J］. Proceedings of the Combustion In⁃ stitute，2009，32（2）：2291-2298. 3］ 金朋刚，郭炜，王建灵，等. TNT ［ 密闭条件下 的爆炸压力特性［J］. 火炸药学报，2013，36（3）：39-41. Jin P G ，Guo W，Wang J L，et al. Explosion pressure characteristics of TNT under closed condition［J］. Chi⁃ nese Journal of Explosives & Propellants，2013，36 （3）：39-41（in Chinese）. 4］ 李鸿宾，金朋刚，郭炜，等. ［ 炸药在密闭空间中爆炸超压测试与分析［J］. 科学技术与工程，2013，13（28）： 8448-8451. Li H B Jin P G ，Guo W，et al. Overpressure test and ， analysis of TNT blast in confined chamber［J］. Science Technology and Engineering， 2013， 13 （28 ）： 8448-8451（in Chinese）.
5］ 李芝绒，王胜强，殷俊兰. ［ 不同气体环境中温压炸药
J］. 2013，36爆炸特性的试验研究［ 火炸药学报， （3）：59-61. Li Z R ，Wang S Q，Yin J L. Experiment study of blast performance of thermobaric-explosive under different gas environment［J］. Chinese Journal of Explosives ＆ Propellants，2013，36（3）：59-61（in Chinese）. ［6］ Kuhl A L，Oppenheim A K，Ferguson R E. Multi-flu⁃ id model of exothermic fields in explosions［C］//Pro⁃ ceedings of the 12nd All-Union Symposium on Com⁃ bustion and Explosion. Chernogolovka， Moscow Re⁃ gion：Lawrence Livermore National Lab，2000. ［7］ Togashi F，Baum J D，Mestreau E，et al． Numerical simulation of long-duration blast wave evolution in con⁃ fined facilities ［J］. Shock Waves，2010，20（5）： 409-424. ［8］ Togashi F，Baum J D，Soto O A，et al．Numerical simu⁃ lation of TNT-Al explosives in explosion chamber［C］// Proceedings of the 7th International Conference on Computational Fluid Dynamics. Hawaii： ICCFD7， 2012. ［9］ Liu X D ， Osher S， Chan T. Weighted essentially non-oscillatory schemes［J］. Journal of Computational Physics，1994，115（1）：200-212. ［10］ Jiang G S，Shu C W. Efficient implementation of weighted ENO schemes［J］. Journal of Computational Physics，1996，126（1）：202-228. ［11］ Shi J，Zhang Y T，Shu C W. Resolution of high order WENO schemes for complicated flow structures［J］. Journal of Computational Physics，2003，186（2）： 690-696. ［12］ Miller P J. A reactive flow model with coupled reac⁃ tion kinetics for detonation and combustion in non-ideal explosives［C］//Proceedings of the Materi⁃ als Research Society Symposium Proceedings. China Lake，CA：Naval Air Warfare Center，1995，418. 13］ 徐维铮，吴卫国. WENO ［ 格式精度对舱室内爆炸载究［J］. 舰船科学技术，2018，40（1）：荷影响规律研9-16. Xu W Z ，Wu W G. The influence of accuracy of WE⁃ NO schemes on the blast load in confined cabin［J］. Ship Science and Technology，2018，40（1）：9-16 （in Chinese）. 14］ 徐维铮，孔祥韶，郑成，等. ［ 约束空间爆炸冲击波高精度数值计算及实验验证［J］.哈尔滨工程大学学报，2017，38（7）：1151-1161. Xu W Z ，Kong X S，Zheng C，et al. High-resolution hydrocode and validation of explosive shock wavein a confined space［J］. Journal of Harbin Engineering University，2017，38（7）：1151-1161（in Chinese）. 15 徐维铮，吴卫国. ［ ］ 爆炸波高精度数值计算程序开发
及应用［J］. 中国舰船研究，2017，12（3）：64-74. Xu W Z ，Wu W G. Development of in-house highresolution hydrocode for assessment of blast waves and its application［J］. Chinese Journal of Ship Re⁃ search，2017，12（3）：64-74（in Chinese）. 16］ 徐维铮，吴卫国. ［ 泄压口大小对约束空间爆炸准静
律［J］. 高压物理学报，2017，态超压载荷的影响规31（5）：619-628. Xu W Z ，Wu W G. Effects of size of venting holes on the characteristics of quasi-static overpressure in con⁃ fined space［J］. Chinese Journal of High Pressure Physics，2017，31（5）：619-628（in Chinese）.