Research Survey on the Mechanism of the Synthesis of Ultrafine Diamond by Means of Explosives Detonation
-
摘要: 从实验探索、理论研究和数值模拟三方面对利用炸药分子中的碳合成超细金刚石的机理进行了系统的总结与分析,并介绍了目前这个领域的一些有待研究的问题。Abstract: In the article, the mechanism of ultrafine diamond synthesized from the carbon of explosives, are given a systematic summary and analysis from three aspects, experimental explores theoretical research and numerical modeling. Some problems expected to be researched at present are proposed.
-
弹药在储存、运输以及使用过程中,一旦遭受意外热刺激,则易发生燃烧、爆炸甚至爆轰,因此有必要对炸药的热安全性进行合理评估。在确保弹药安全的前提下提高炸药威力是目前亟待解决的问题。在混合炸药配方中添加氧化剂可改善其氧平衡,较大幅度地提升爆热。高氯酸铵(NH4ClO4,ammonium perchlorate,AP)具有密度较高、含氧量高、爆轰产物为气体、机械感度低等优点,是混合炸药配方的重要候选组分[1]。研究含有AP的浇注炸药的热响应特性对其推广应用以及新型钝感高聚物黏结炸药(plastic bonded explosives,PBX)研发具有重要意义。
国内外对含AP的含能材料热分解特性以及热安全性开展了较多研究。刘子如等[2]研究了AP-奥克托今(HMX)混合物的热分解特性,发现AP的存在使HMX的微商热重(derivative thermogravimetry,DTG)峰温大幅下降,证明AP与HMX在热分解中存在强烈的相互作用。Li等[3]研究了铝粉的加入对AP-HMX混合物热分解特性的影响,发现足量的铝粉可以消除AP与HMX之间的相互催化作用。陈中娥等[4]认为,AP热分解形成的孔隙是AP-HTPB推进剂慢速烤燃响应剧烈的原因。李苗苗等[5]对不同AP含量的HTPE推进剂进行了差示扫描量热法测试及慢速烤燃实验,发现随着AP含量的降低,推进剂的热稳定性和热安全性略有提升。Kim等[6]对含有AP的PBX炸药进行了快速和慢速烤燃实验,结果表明,在慢速烤燃中AP不参与点火,而在快速烤燃中点火由AP与硝胺类炸药共同参与。由此可见,AP在热刺激下的自热反应受到添加剂和温升速率等因素的影响,但是目前仍缺乏含AP的PBX炸药的热刺激响应特性相关研究[7]。
在炸药的热刺激响应特性实验中,温升速率对炸药的点火时间、点火位置以及装药的温度梯度等产生很大的影响。王沛等[8]对GHL炸药进行了慢速烤燃实验以及4种不同温升速率下的数值模拟,发现点火位置随着温升速率的增大从炸药内部移向炸药边缘。代晓淦等[9]进行了不同温升速率下PBX-2炸药的烤燃实验和数值模拟,发现温升速率越大,炸药中心温度在烤燃过程中的变化越小。刘静等[10]对模块装药进行了慢速烤燃实验和数值模拟,发现随着温升速率的提高,烤燃响应区域由1个中心环形区域变成关于中心对称的2个环形区域。由此可见,温升速率对炸药烤燃中的点火位置以及点火区域形状的影响较大,对获取热刺激下炸药的响应特性具有重要意义。
本研究将开展新型浇注炸药GOLA-1的慢速烤燃实验和数值模拟,分析GOLA-1炸药在不同条件下的热刺激响应特性及其在慢速烤燃中点火的主导因素,探讨温升速率对小型烤燃弹点火位置的影响。
1. 实验与数值模拟方法
1.1 小型烤燃弹实验
基于小型烤燃弹实验装置对GOLA-1炸药进行慢速烤燃实验。图1为小型烤燃弹实验测试系统示意图。该装置主要包括壳体、端盖、炸药、环形加热器、测温热电偶以及控温系统[11]。烤燃弹的外径为58 mm,长度为106 mm,壳体厚度为4 mm,药柱的直径为50 mm,长度为88 mm。药柱与端盖之间有高12 mm的空气域,用以缓解炸药受热膨胀对壳体的压力。端盖中心以及炸药上表面至炸药中心均开有直径为1 mm的小孔,用于放置K型热电偶,监测温升过程中炸药中心点的温度变化。壳体和端盖材质均为45钢,两者采用螺纹连接。装药壳体外侧包裹环形加热器,在慢速烤燃实验中温升速率为1.0 和1.5 K/min。加热器外部使用岩棉包裹进行保温,防止热量散失影响实验结果。在壳体与加热器之间放置两根K型热电偶,用来控制加热速率以及测量壳体外壁面的温度变化。
1.2 烤燃数值模拟
烤燃实验只能获得温度监测点的温度-时间关系,无法获取炸药内部的温度分布、点火位置、点火温度以及各组分的质量分数等信息。数值模拟可以较方便地获取以上信息,从而对炸药的热刺激响应特性进行较全面的评估。基于小型烤燃弹实验装置,建立如图2所示的数值模拟模型,其中端盖及壳体、空气域以及炸药分别用粉色、蓝色和红色表示。模型用六面体网格划分,网格尺寸为0.5 mm。
在直角坐标系下,热传导方程的形式为[12]
ρcp∂T∂t=λ∇2T+S (1) 式中:ρ为密度,kg·m−3;cp为比定压热容,J·kg−1·K−1;T为温度,K;λ为热导率,W·m−1·K−1;S为炸药的自热反应源项,W·m−3。
空气域与炸药之间的辐射热流qr采用P1辐射模型描述[13]
qr=13(a+σs)−Cσs∇G (2) 式中:a为吸收系数,m−1;σs为散射系数,m−1;C为线性各向异性相位函数系数;G 为入射辐射。
GOLA-1为HMX基浇注炸药,主要成分为HMX,其他成分包括Al、 AP、黏结剂(Kel-F)。组分HMX、Al、AP、Kel-F的质量比为44∶25∶20∶11,密度为1.837 g/cm3。对于炸药中的主要成分HMX,在慢速烤燃过程中,随着温度上升,HMX首先由β相转变为δ相,而后热分解生成最终的气态产物。基于过渡态理论,HMX在烤燃过程中的4步化学反应动力学模型[14]表示为:(1) β-HMX ↔ δ-HMX(反应1,单分子吸热反应);(2) β-HMX+δ-HMX ↔ δ-HMX(反应2,双分子吸热反应);(3) δ-HMX→产物(反应3,单分子吸热反应);(4) δ-HMX+产物→产物(反应4,双分子放热反应)。反应1和反应2描述β相HMX转变为δ相HMX的晶型转变过程,反应3和反应4描述最终产物生成的点火过程。各步反应速率方程为
r1=wβkThexp(TΔS1f−E1fRT)−wδkThexp(TΔS1r−E1rRT) (3) r2=wβwδkTh[exp(TΔS2f−E2fRT)−exp(TΔS2r−E2rRT)] (4) r3=wδZ3exp(−E3RT) (5) r4=wδwpZ4exp(−E4RT) (6) 式中:r1~r4为反应1~反应4的反应速率;wβ和wδ分别为β相和δ相HMX的质量分数,wp为产物的质量分数;k为玻尔兹曼常数;h为普朗克常数;R为普适气体常数;ΔS为活化熵,J∙mol−1∙K−1;E为活化能,J∙mol−1;下标f代表正反应,r代表逆反应;Z1~Z4为反应1~反应4的指前因子,s−1。
各组分质量分数的微分方程为
dwβdt=−r1−r2 (7) dwδdt=r1+r2−r3−r4 (8) dwpdt=r3+r4 (9) 单位时间内HMX热分解反应的反应源项SHMX为
SHMX=ρHMXQ1r1+ρHMXQ2r2+ρHMXQ3r3+ρHMXQ4r4 (10) 式中:ρHMX为HMX的密度,kg∙m−3;Q1~Q4为反应1~反应4的反应热,J∙kg−1。
Kel-F为含氟含氯黏结剂,分解时吸收热量,使用单步化学反应动力学模型描述其热分解过程。其热分解反应速率方程为[15]
r5=Z5exp(−Ea5RT) (11) 式中:r5、Z5、Ea5分别为Kel-F热分解反应的反应速率、指前因子和活化能。单位时间内Kel-F热分解反应的反应源项SKel-F为
SKel-F=ρKel-FQ5r5 (12) 式中:ρKel-F为Kel-F的密度,Q5为Kel-F热分解反应的反应热。
AP为氧化剂,其分解反应为放热反应,采用一步反应描述其热分解过程[16]
NH4ClO4→O2+H2O+HCl+N2 AP在240 ℃左右由斜方晶型(α-AP)转变为立方晶型(β-AP)[17]。在慢速烤燃的温升过程中,AP的温度几乎均在240 ℃以下,因此采用α-AP的热分解反应速率方程[18],即
r6=dαdt=Z6exp(−Ea6RT)αn(1−α)m (13) 式中:r6、Z6、Ea6分别为α-AP热分解反应速率、指前因子和活化能;α为AP的热分解转化率,初始阶段为零,完全分解为1;n和m为反应动力学参数,对于α-AP,n=1,m=1.38。根据文献[19],AP的热分解开始于150 ℃左右,因此在UDF用户子程序中设α在150 ℃以下为零,150 ℃以上为10−4,以启动对AP分解反应计算。设AP的密度为ρAP,则单位时间内AP热分解反应的反应源项SAP为
SAP=ρAPQ6r6 (14) 式中:Q6为AP热分解反应的反应热。
因为Al不参与点火前的热分解反应,所以不考虑铝粉的热分解反应源项,则炸药的总自热反应源项S为
S=0.44SHMX+0.11SKel-F+0.20SAP (15) 根据叠加原理,炸药的导热系数λ和比热容C分别表示为
λ=0.44λHMX+0.11λKel-F+0.20λAP (16) C=0.44CHMX+0.11CKel-F+0.20CAP (17) 表1和表2列出了炸药组分HMX、Kel-F和AP的反应动力学参数,表3列出了材料的物性参数。采用Fluent软件进行计算,设置温升边界条件为1.0和1.5 K/min,通过编写用户自定义函数UDF,将热导率、边界条件以及反应动力学参数加载到主程序中进行计算。
表 1 HMX的反应动力学参数Table 1. Reaction kinetic parameters of HMXReaction No. Z/s−1 ΔSf/
(J∙mol−1∙K−1)Ef/(kJ∙mol−1) ΔSr/ (J∙mol−1∙K−1) Er/(kJ∙mol−1) Q/(J∙g−1) 1 123.0 204 89.0 189 −25 2 −41.7 102 −75.2 865 −25 3 3.16×1016 200 −120 4 1.50×1014 173 3200 表 2 Kel-F和AP的反应动力学参数Table 2. Reaction kinetic parameters of Kel-F and APComponent Reaction No. Z/s−1 E/(kJ∙mol−1) Q/(J∙g−1) Kel-F 5 9.93×1017 272.0 −5870 AP 6 6.85×1012 146.3 297 表 3 材料物性参数Table 3. Parameters of material propertiesComponent Density/(kg∙m−3) C/(J·kg−1·K) λ/(W·m−1·K) HMX 1850 1004.26 0.5358 Kel-F 2020 1000.43 0.0527 AP 1950 1255 0.2760 Al 2719 871 1.39 Steel 8030 502.48 16.27 Air 1.225 1006.43 0.0242 2. 结果与讨论
2.1 烤燃实验
图3为不同温升速率下GOLA-1炸药点火后烤燃弹照片。图3(a)显示:点火后烤燃弹壳体膨胀变形,端盖被冲飞,烤燃弹底部断裂脱落,壳体螺纹完整;壳体表面有明显的烧灼痕迹,壳体内有少量未反应炸药。图3(b)显示了1.5 K/min温升速率下GOLA-1炸药点火前后烤燃弹的对比。结果显示,点火后烤燃弹壳体膨胀变形,端盖被冲飞,壳体底部与侧面有明显的烧灼痕迹,壳体内有少量未反应炸药。由实验现象判定,当温升速率为1.0和1.5 K/min时,GOLA-1炸药的热刺激响应等级均为燃烧,说明炸药在慢速烤燃条件下的热安全性较好。对比两种温升速率下炸药的响应程度,可以发现,温升速率为1.0 K/min时炸药的慢速烤燃响应更剧烈。
2.2 数值模拟
基于实验进行慢速烤燃数值计算,在1.0和1.5 K/min温升速率下实验与数值模拟得到的炸药中心测温点温度-时间曲线如图4和图5所示。可以看出,实验与数值模拟得到的温升曲线基本重合,点火温度和点火时间也基本吻合,说明数值模型及参数比较合理。图4显示,9900 s时实验测得的温升曲线的温升放缓,而后恢复正常温升速率。分析认为,此阶段的温升放缓是由于HMX在由β相转变为δ相的过程中吸收热量,温度约为451 K。在慢速烤燃实验中,炸药于11946 s时发生点火反应,点火时药柱中心点温度为488 K;计算得到的点火时间为11790 s,与实验结果的相对偏差为1.3%。图5显示:HMX相变造成的温升速率放缓出现在约6800 s,温度约为451 K;实验结果显示,炸药在约8207 s发生点火反应,点火时刻药柱中心温度约为480 K;计算得到的点火时间为8070 s,与实验结果的相对偏差为1.7%。由此可见,温升速率对点火时间的影响较大,随着温升速率的上升,点火时间缩短,然而温升速率对HMX的相变温度以及GOLA-1炸药的点火温度没有影响。
1.0和1.5 K/min温升速率下不同时刻炸药的温度分布如图6和图7所示。相比于炸药,壳体的比热容小而热导率大,因此壳体温度高于炸药温度,热量由壳体传递到装药内部。装药的部分热量散失至烤燃弹上部的空气域中,造成装药上部温度较低。图6(d)和图7(d)为点火时烤燃弹的温度分布,可见,点火区域位于装药底部的环形区域。这是由于药柱底部受到侧面壳体与底面的共同加热,使得装药底部温度相对较高。温度的升高使得炸药的自热分解反应速率上升,造成热量进一步累积,如此循环往复,致使点火最终发生在装药底部的环形区域。对比两者发现,在1.5 K/min的温升速率下装药的温度梯度较大,点火位置相对靠外。分析认为,这是由于炸药的导热系数远低于钢壳体,而温升速率相对较大,热量来不及传递到炸药内部造成的。
1.0 K/min温升速率下不同位置的温升曲线和质量分数曲线如图8、图9和图10所示。从图8可以看出,HMX在约9300 s发生相变,相变温度与中心点温度相同,均为451 K。对比中心测温点与点火点的温升曲线可以发现,点火点处HMX相变所造成的温度平台较小,且持续时间较短。点火点处的热量累积速率大于热量传递速率,温度上升较快,相同时刻点火点处的温度高于其他位置,从而使HMX相变时间早且速率快。临近点火时,由于HMX热分解反应释放大量的热量,造成点火点的温度急剧升高,点火温度约为517 K。图9和图10给出了中心测温点和点火点的HMX质量分数变化曲线。在点火时刻,β-HMX已经全部转变为δ-HMX,相变过程基本完成。点火时中心测温点仅有微量的HMX热分解产物,而点火位置则有质量分数为0.2的产物。说明在点火位置,第4步强放热反应已经开始进行,释放大量热量,导致温度急剧上升,主导了点火的发生。
1.5 K/min温升速率下不同位置的温升曲线和质量分数曲线如图11、图12和图13所示。图11中的点火点温度-时间曲线特征与图8基本相同,点火温度约为518 K。图12和图13显示:相比于温升速率为1.0 K/min的情况,点火时刻中心测温点处β-HMX的质量分数为0.04,δ-HMX的质量分数为0.96,表明部分β-HMX尚未转变为δ-HMX,第4步强放热反应尚未发生;而点火位置的HMX已完全相变,并且HMX分解的第4步强放热反应已经开始进行,有质量分数为0.19的产物生成,主导了点火的发生。
图14给出了温升速率为1.5 K/min时中心测温点处AP的转化率-时间曲线。转化率α与分解速率dα/dt均随着温度的升高而增大。在点火时刻,AP的转化率为1.06×10−4,表明此时中心点的AP基本不发生分解。7000 s左右,AP的分解速率放缓,而后快速上升,这是HMX在晶型转变过程中吸收热量而导致温升放缓造成的。图15为点火点处AP的转化率-时间曲线。可以看出:转化率在临近点火时刻迅速上升;点火时刻AP的转化率为2.26×10−4,相较于中心点处的转化率高约1倍。不论是中心测温点还是点火点,发生分解的AP的量均很少,表明在慢速烤燃过程中AP基本不参与反应,点火反应由HMX主导。这与Kim等[6]得到的含硝胺类及AP浇注炸药慢速烤燃实验结论一致。
2.3 温升速率对点火位置的影响
将实验与数值模拟相结合,分析温升速率对炸药热刺激响应特性的影响。在1.0和1.5 K/min两种温升速率的基础上,设定不同的温升速率,对GOLA-1炸药进行计算。当温升速率为0.5和0.4 K/min时,炸药点火位置的变化如图16所示。0.5 K/min时,点火位置仍为装药底部的环形区域,相比1.0和1.5 K/min的情况,环形区域缩小,并向装药的上部及中心移动。这是由于当温升速率下降时,炸药整体受热比较均匀,装药内部的温度梯度减小,相比于温升速率为1.0 K/min的情况,相同加热温度下炸药内部高温区域的体积较大,炸药的自热反应速率上升,炸药内部积累热量,导致点火发生在装药底部的环形区域。当温升速率降至0.4 K/min时,点火位置的形状发生改变,由环形变为扁的椭球形,并且位置移动至药柱下部中心。
为进一步分析温升速率对GOLA-1点火位置的影响,以0.05 K/min为间隔,进行0.1~0.4 K/min 7个不同温升速率的烤燃数值模拟。设温升速率为β,以壳体外侧底面的中心位置为坐标零点,沿烤燃弹中心轴线从壳体底面指向端盖为正方向,点火位置x为温度最高区域的中心,单位为m。温升速率对点火位置的影响如表4所示。随着温升速率的下降,点火位置向药柱上部移动。根据表4绘制点火位置随温升速率的散点图,如图17所示。观察图17发现,温升速率β与点火位置x之间存在线性关系,拟合关系式为
x=−0.05279β+0.04171 ,相关系数 R2=0.966。表 4 温升速率与点火位置的关系Table 4. Heating rate vs. ignition positionNo. β/(K·min−1) x/m 1 0.40 0.0212 2 0.35 0.0220 3 0.30 0.0273 4 0.25 0.0287 5 0.20 0.0298 6 0.15 0.0334 7 0.10 0.0374 3. 结 论
(1) 对新型HMX基含AP浇注炸药GOLA-1进行慢速烤燃实验,结果显示,当温升速率为1.0和1.5 K/min时炸药的烤燃响应等级均为燃烧,并且1.0 K/min下壳体的破坏更严重,响应更剧烈。从热刺激响应角度来看,温升速率越低,GOLA-1炸药的响应越剧烈。
(2) 开展了GOLA-1炸药烤燃数值模拟,得到的温升曲线与实验结果吻合较好。当温升速率为1.0 K/min时,GOLA-1炸药的点火时间与实验结果的相对偏差为1.3%,1.5 K/min时为1.7%,证明建立的烤燃模型及参数较合理。数值模拟结果显示,点火反应主要由HMX的第4步强放热反应主导,AP在慢速烤燃过程中基本不参与点火反应。
(3) 对不同温升速率下GOLA-1炸药的烤燃过程进行了数值模拟,结果表明:当温升速率为0.5 K/min时,点火仍发生在装药底部的环形区域,并向装药中心及上部移动;当温升速率下降至0.4 K/min时,点火区域的形状发生变化,由环形变为椭球形;随着温升速率继续下降,点火位置沿着药柱中心轴线向上部移动;当温升速率处在0.1~0.4 K/min之间时,温升速率与点火位置之间存在线性关系。
-
Urizar M J, James E Jr, Smith L C. Phys Fluids, 1960, 4: 262. Steinberg D J. Proc. of the 8th Symp. Int on Detonation. Maryland: NSWC White Oak, Silver Spring, 1985: 513. Staver A M, Gubareva N V, et al. ФГВ, 1984, 20(5): 100. Roy Greiner N, Robert Hermes. Proc of the 9th Symp Int on Detonation. Virginia: OCNR Arlington, 1989: 1170. 徐康, 金增寿, 等. 含能材料, 1993, 1(3): 19. Kuznetsov V L, Aleksandrov M N, Zagoruiko Z V, et al. Carbon, 1991, 29(4-5): 665-668. Vereschagin A L, Sakovich G V, Komarov V F, et al. Diamond and Related Materials, 1993, 3: 160-162. Roy Greiner N, Philips D S, Johnson J D, et al. Nature, 1988, 333: 440-442. Jhonson J D. Proc. of the 9th Sypm. Int. on Detonation. Virginia: OCNR Arlington, 1989: 417. Titov V M, Anisichkin V F. Proc. of the 9th Sypm. Int. on Detonation. Virginia: OCNR Arlington, 1989: 407. Staver A M, Ershov A P, Lyamkin A Z. ФГВ 1989, 20(3): 79. Anisichkin V F, Derndyaev B G, et al. ФГВ, 1988, 24(3): 121-122. Mcguire R R, Ornellas D L, Akst I B. Propellants Explosives Pyrotechnics, 1979, 4(2): 23-26. Shaw M S, Johnson J D. J Appl Phys, 1987, 62(5): 2080. Bdzil J B. Proc. of the 6th Symp. Int. on Detonation. Virginia: ONRN Arlington, 1976: 352. Tarver C M, Breithaupt R D, Kury J W. Proc. of Int. Symp. on Pyrotechnics and Explosives. Beijing, 1987. Mader C L, Craig B G. LA-5865, 1975. Badziag P, Verward W S, Ellis W P, et al. Nature, 1990, 343: 244-245. Guire R M, Ornellas D, Akst L. Symp HPD. Paris, 1978. Petrov Yu I. Clusters and Small Particles. Moscow: Nauka Publishers, 1986. Novikov N V. High Pressures in Sci, and Tech. Book of Abstracts of 11th AIRAPT Int Conf. Kiev, 1987. Fol' mer M. Kinetics of New Phase Formation. Moscow: Nauka Publishers, 1986. Landau L D, Lifshitz E M. Statistical Phys (2rd ed). Moscow: Nauka Publishers, 1964: 556. Babare L V, Dremin A N, Pershin S V, et al. ФГВ, 1969, 5(4): 528-539. Yakusheva O B, Yakushev V V, Dremin A N. ЖФХ, 1977, 51(7): 1657-1660. Yakusheev V V, Dremin A V, Nabatov S S, et al. ФГВ, 1970, 15(2): 132-140. Baute J, Chirat R. Proc. of the 8th Symp. Int. on Detonation. Maryland: NSWC White Oak, Silver Spring, 1985: 521-531. Mat Van Thiel, Francis H Ree. J Appl Phys, 1987, 62(5): 1761. Whittaker A G. Carbon, 1979, 17: 21. Hirano S, Shimono K, Naka S. J Mater Sci, 1984, 19: 1153. Pierk Tang. J Appl Phys, 1988, 63(4): 1041. Francis H Ree. J Chem Phys, 1984, 81(3): 1251. Francis H Ree. J Chem Phys, 1986, 84(10): 1. Francis H Ree, Mat Van Thiel. Proc. of the 8th Symp. Int. on Detonation. Maryland: NSWC White Oak, Silver Spring, 1985: 425. Yamada K, Sawaoka A B. Carbon, 1994, 32(4): 665. Franklin E Walker. Propellants Explosives Pyrotechnics, 1990, 15: 190-193. -

计量
- 文章访问数: 7697
- HTML全文浏览量: 565
- PDF下载量: 714