Fluid Dynamics Analysis of Early-Stage Shock Initiation ofPBX 9501 Using Improved CIM
-
摘要: 利用特征线法、简化的状态方程、炸药冲击阻抗与声阻抗的关系,可以得到凝聚炸药冲击起爆前期反应流场的解析模型(CIM2)。选择具有单一形式、带耗散项的反应速率函数,进而计算出不同深度炸药的粒子速度剖面。通过对不同深度反应速率的重新假设,改进了CIM2模型,并将其应用于PBX 9501炸药。结果表明,改进模型减缓了CIM2模型中的能量输送现象,从而能够更准确地计算冲击波阵面粒子速度、峰值粒子速度、峰值到达时间和粒子速度曲线形状。Abstract: In this work we obtained an analytic model CIM2 for the flow field of the reaction zone of condensed explosives in early-stage shock initiation using the characteristic line, the simplified equation of state and the simplified relationship between the material impedance and the shock impedance.By choosing a single form of the reaction rate with the dissipative term, we calculated the particle velocity at different depths of condensed explosives.The CIM2 model was improved by the new hypothesis of the reaction rates at different depths, and was applied to the calculation of PBX 9501.The results show that the improvement of CIM2 can slow down the energy transport after the shock, and that the improved CIM2 can give a more exact calculation of the particle velocity of the shock and crest, the time at which the crest arrives, and the curve of the particle velocity.
-
Key words:
- characteristic line /
- reaction rate /
- constant impedance model /
- PBX 9501
-
表 1 炸药及飞片参数
Table 1. Parameters of explosive and flyer
Material ρ0/(g/cm3) c/(km/s) ud/(km/s) σ U/(km/s) PBX 9501 1.842 5.28 3.80 0.25 Al flyer 2.785 5.33 0.75 表 2 不同深度采用的反应度-时间关系式
Table 2. Relationship of extent of reaction and time at different depths
d/(mm) Reaction rate function 0 λ=1-exp(-1.20τ2.50/2.50) 0.93 λ=1-exp(-1.08τ2.55/2.55) 1.94 λ=1-exp(-1.15τ2.75/2.75) 2.93 λ=1-exp(-1.20τ2.90/2.90) 3.92 λ=1-exp(-1.30τ3.15/3.15) 4.91 λ=1-exp(-4.00τ5.50/5.50) -
[1] LAMBOURN B D.An interpretation of particle velocity histories during growth to detonation [J].AIP Conf Proc, 2004, 706:367-370. doi: 10.1063/1.1780255 [2] JAMES H R, LAMBOURN B D.On the systematics of particle velocity histories in the shock-to-detonation transition regime [J].J Appl Phys, 2006, 100(8):084906. doi: 10.1063/1.2354416 [3] HILL L R, GUSTAVSEN R L, ALCON R R, et al.Shock initiation of new and aged PBX 9501 measured with embedded electromagnetic particle velocity gauges: LA-13634-MS [R].Los Alamos: Los Alamos National Laboratory, 1999. [4] 陈军, 田占东, 曾代朋, 等.PBX 9501炸药冲击起爆初期反应流体动力学特性研究[J].高压物理学报, 2013, 27(3):331-336. http://www.gywlxb.cn/CN/abstract/abstract1571.shtmlCHEN J, TIAN Z D, ZENG D P, et al.Reactive hydrokinetics in the early stage of shock initiation of PBX 9501 [J].Chinese Journal of High Pressure Physics, 2013, 27(3):331-336. http://www.gywlxb.cn/CN/abstract/abstract1571.shtml [5] 李述涛.固体炸药起爆早期反应速率研究[D].长沙: 国防科学技术大学, 2009: 25-32. http://cdmd.cnki.com.cn/Article/CDMD-90002-2010165540.htmLI S T.Analysis of the reactive rate for solid explosives during the early stage to detonation [D].Changsha: National University of Defense Technology, 2009: 25-32. http://cdmd.cnki.com.cn/Article/CDMD-90002-2010165540.htm