密闭多舱室内柱形气云爆炸过程二维数值模拟
|
摘要
为了研究密闭多舱室空间内爆炸波传播过程及舱室壁面爆炸载荷分布特性,采用三阶WENO有限差分格式,编写了多舱室内爆炸冲击波二维数值计算程序。利用平面激波与矩形方腔相互作用实验验证了本文数值程序的可靠性与正确性。采用验证的程序开展了密闭多舱室内柱形高压气云爆炸波数值计算研究,探讨了多舱室内爆炸波传播路径,不同舱室内部爆炸载荷特性及爆源位置对爆炸载荷特性的影响规律。本文研究可为工程抗爆结构设计及爆炸毁伤评估提供一定的参考和指导。
In order to study on the blast wave propagation progress and the distribution characteristics of the blast load in multi-cabin, third order WENO finite difference scheme is adopted to develop an in-house two-dimensional numerical calculation code for simulating the blast in closed multi-cabin. The reliability and correctness of the numerical program are verified by the interaction between the plane shock wave and the rectangular cavity. Then, the validated code is used to conduct numerical simulation researches on the blast of cylindrical high pressure gas. The blast wave propagation path, blast load distribution characteristics in different cabins and the influence of explosion source position on blast load are investigated. The research in the present work can provide some references and guidance for the design of anti-explosion structure and evaluation of explosion damage.
In order to study on the blast wave propagation progress and the distribution characteristics of the blast load in multi-cabin, third order WENO finite difference scheme is adopted to develop an in-house two-dimensional numerical calculation code for simulating the blast in closed multi-cabin. The reliability and correctness of the numerical program are verified by the interaction between the plane shock wave and the rectangular cavity. Then, the validated code is used to conduct numerical simulation researches on the blast of cylindrical high pressure gas. The blast wave propagation path, blast load distribution characteristics in different cabins and the influence of explosion source position on blast load are investigated. The research in the present work can provide some references and guidance for the design of anti-explosion structure and evaluation of explosion damage.
引文
[1]徐胜利,汤明钧.近地空中气云爆炸波遇地面反射的研究[J].爆炸与冲击,1996(4):298-304.XU Sheng-li,TANG Ming-jun.Studies on reflection of blast waves for symmetric cloud explosion close to the ground[J].Explosion and Shock Waves,1996(4):298-304.
[2]徐胜利,岳鹏涛,彭金华.多爆源云雾爆炸波相互作用的三维数值研究[J].爆炸与冲击,2000,20(1):1-6.XU Sheng-li,YUE Peng-tao,PENG Jing-hua.Three dimensional computaion on the interaction of blast waves generated by multi-sources of FAE[J].Explosion and Shock Waves,2000,20(1):1-6.
[3]岳鹏涛,徐胜利,彭金华.FAE爆炸波对地面目标作用的三维数值研究[J].爆炸与冲击,2000,20(2):97-102.YUE Peng-tao,XU Sheng-li,PENG Jing-hua.3D numerical simulations on the interaction between FAE blast waves and ground targets[J].Explosion and Shock Waves,2000,20(2):97-102.
[4]陈明生,李建平,白春华.非圆截面云雾爆炸超压场数值模拟[J].含能材料,2015,23(5):484-489.CHEN Ming-sheng,LI Jian-ping,BAI Chun-hua.Simulation of explosion overpressure distribution for non-circular crosssection cloud[J].Chinese Journal of Energetic Materials,2015,23(5):484-489.
[5]LIU X D,OSHER S,CHAN T.Weighted essentially nonoscillatory schemes[J].J Comput Phys,1994,115:200-212.
[6]JIANG G S,SHU C W.Efficient implementation of weighted ENO schemes[J].J Comput Phys,1995,126:202-228.
[7]SHI J,ZHANG Y T,SHU C W.Resolution of high order WENO schemes for complicated flow structures[J].J Comput Phys,2003,186:690-696.
[8]IGRA O,FALCOVITZ J,REICHENBACH H,et al.Experimental and numerical study of the interaction between a planar shock wave and a square cavity[J].Journal of Fluid Mechanics,1996,313:105-130.
[2]徐胜利,岳鹏涛,彭金华.多爆源云雾爆炸波相互作用的三维数值研究[J].爆炸与冲击,2000,20(1):1-6.XU Sheng-li,YUE Peng-tao,PENG Jing-hua.Three dimensional computaion on the interaction of blast waves generated by multi-sources of FAE[J].Explosion and Shock Waves,2000,20(1):1-6.
[3]岳鹏涛,徐胜利,彭金华.FAE爆炸波对地面目标作用的三维数值研究[J].爆炸与冲击,2000,20(2):97-102.YUE Peng-tao,XU Sheng-li,PENG Jing-hua.3D numerical simulations on the interaction between FAE blast waves and ground targets[J].Explosion and Shock Waves,2000,20(2):97-102.
[4]陈明生,李建平,白春华.非圆截面云雾爆炸超压场数值模拟[J].含能材料,2015,23(5):484-489.CHEN Ming-sheng,LI Jian-ping,BAI Chun-hua.Simulation of explosion overpressure distribution for non-circular crosssection cloud[J].Chinese Journal of Energetic Materials,2015,23(5):484-489.
[5]LIU X D,OSHER S,CHAN T.Weighted essentially nonoscillatory schemes[J].J Comput Phys,1994,115:200-212.
[6]JIANG G S,SHU C W.Efficient implementation of weighted ENO schemes[J].J Comput Phys,1995,126:202-228.
[7]SHI J,ZHANG Y T,SHU C W.Resolution of high order WENO schemes for complicated flow structures[J].J Comput Phys,2003,186:690-696.
[8]IGRA O,FALCOVITZ J,REICHENBACH H,et al.Experimental and numerical study of the interaction between a planar shock wave and a square cavity[J].Journal of Fluid Mechanics,1996,313:105-130.