Effect of Pressure Environment on Graphitization of Nano-Crystalline Diamond
-
摘要: 高压下纳米金刚石石墨化对多晶金刚石烧结体的性能具有非常重要的影响。为此,在5~9 GPa、600~1500 ℃的压力温度范围,分别对平均晶粒尺寸为50 nm的纯纳米金刚石粉末和NaCl-纳米金刚石混合粉末的石墨化温度进行了研究。结合粉末X射线衍射的物相定量分析方法,分析了非静水压(纯金刚石粉末)和准静水压(NaCl-纳米金刚石混合粉末)下纳米金刚石在不同压力和温度下的石墨化程度。结果表明:5 GPa时,纯纳米金刚石粉末石墨化的起始温度在800 ℃以上,9 GPa时则在1000~1300 ℃之间;在约7 GPa的压力下,较短的保温时间内纳米金刚石的石墨化温度由非静水压环境中的1000 ℃提高到准静水压环境中的1500 ℃以上。Abstract: The graphitization process of nano-crystalline diamond (NCD) under high pressure significantly influences the performance of sintered polycrystalline diamond bulks. Here, we investigated the graphitization temperature of both pure nano-diamond, with an average grain size of 50 nm, and a mixture of NaCl and nano-diamond powder in a pressure and temperature range of 5−9 GPa and 600−1500 ℃, respectively. With a quantitative analysis method employing powder X-ray diffraction, we analyzed the graphitization degree of NCD under different pressures and temperatures, examining both non-hydrostatic pressure conditions (pure NCD powder) and quasi-hydrostatic pressure conditions (NaCl-NCD mixed powder). Our findings indicate that the initial graphitization temperature of pure NCD powder exceeds 800 ℃ at 5 GPa, and ranges between 1000 and 1300 ℃ at 9 GPa. Notably, under quasi-hydrostatic pressure conditions at about 7 GPa, the graphitization temperature of NCD increases from 1000 ℃ in non-hydrostatic pressure conditions to 1500 ℃ or higher within a short holding time.
-
海床已成为一个全新的冲突领域,近海底爆炸可以对海底光缆、海底管道等设施造成严重的破坏,其破坏过程涉及近海底反射、多相流掺混、结构与流体的耦合作用,开展近海底爆炸冲击波时空演化规律研究对于水中兵器研制及海底设施抗爆抗冲击设计具有重要的意义。
目前,国内外对近海底水下爆炸研究已经取得了一定成果。基于Kelvin的冲击理论,Blake等[1]得出了判断自由面附近气泡射流方向和气泡运动方向的方法,该方法对近壁面也同样适用。张永坤[2]对沉底水雷在水下爆炸作用下的毁伤情况进行了研究,总结了水雷在不同情况下被毁伤的特点,为提高水下灭雷武器的作战能力奠定了理论基础。邵建军等[3]基于相似理论,研究了海底对爆炸的影响,发现海底底质变化显著影响水下爆炸的冲击波能和气泡能。杨莉等[4–7]在泥底、砂底及石底3种条件下进行了沉底爆炸试验,发现在3种底质条件下,气泡后期溃灭形态存在差异,最大峰值压力通常出现在装药靠近水底面一侧的斜上方,对于爆距较小的水底面测点,反射波追上入射波会形成马赫波。黄潇等[8]采用镜像法模拟海底边界的影响,发现自由场爆炸气泡比近海底爆炸气泡对潜艇施加的总纵弯矩更大。姚熊亮等[9]基于光滑粒子流体动力学(smooth particle hydrodynamics,SPH)方法模拟沉底水下爆炸,研究了水底底质厚度和炸药当量对冲击波压力的影响。邵宗战等[10]提出了沉底水雷海上爆炸威力的测量方法,给出了冲击波峰值压力的拟合方法及爆炸能量计算方法。前人在近海底水下爆炸冲击波载荷分析方面已取得了一些成果,然而,针对不同底质条件下近海底水下爆炸冲击波载荷时空分布规律的定量研究尚不充分。
本研究拟采用耦合欧拉-拉格朗日(coupled Eulerian-Lagrangian,CEL)方法建立近海底水下爆炸模型,探讨不同底质条件下近海底爆炸冲击波的时空分布规律,以期为近海底设施的抗爆抗冲击结构设计提供支撑。
1. 数值模型与数值验证
1.1 状态方程
本研究采用TNT炸药,爆轰产物通过JWL方程[11]描述
p=A(1−ωR1v)e−R1v+B(1−ωR2v)e−R2v+ωev (1) 式中:A、B、R1、R2及
ω 为与炸药状态有关的常数,v为爆轰产物的相对比容,e为炸药单位质量所蕴含的内能。TNT炸药的材料参数及JWL状态方程参数如表1[12]所示,其中:ρ为密度,D为爆速。水采用Mie-Grüneisen状态方程描述,表达式[11]为
p=pH(1−Γ0η/2)+Γ0ρ0Em (2) 式中:pH为Hugoniot压力;ρ0为初始密度;Em为单位质量内能;Γ0为Grüneisen常数;η为名义体积压缩应变,
η=1−ρ0/ρ 。pH的表达式[11]为pH=ρ0C20η(1−Sη)2 (3) us=C0+Sup (4) 式中:us为冲击速度,up为粒子速度,S和C0为us-up曲线参数。水的状态方程参数列于表2[12]。
在数值模型中,设空气为理想气体,其状态方程[11]为
p+pa=(γ−1)ρEm (5) Em=cV(θ−θz) (6) 式中:γ为绝热指数,pa为外界压力,cV为比定容热容,
θ 为当前温度,θz 为绝对零度。空气的状态方程参数如表3[12]所示。1.2 数值模型
为研究近海底水下爆炸冲击波的传播规律,建立了数值模型,欧拉域尺寸为15.0 m×6.5 m×15.0 m,模拟50 kg TNT在50 m水深处的近海底爆炸,TNT底部距海底0.05 m。将整个欧拉域的边界条件设置为流出无反射,海底底质采用模型1或模型2描述。模型1是一种较软的土壤模型,参考了Ambrosini等[13]和Luccioni等[14]的研究结果,具体参数见表4[13–14],其中:E为杨氏模量,ν为泊松比,φ为内摩擦角,c为黏聚力。模型2为刚性固壁,通过在欧拉域中的海底底质区域放置拉格朗日刚体实现。
为了更准确地获得同一位置的冲击波压力,在模型的不同位置设置测点,测点布局如图1所示。设装药半径为r,测点到爆心的距离为R,测点与水平方向的夹角为测点角度α。在距离爆心分别为7r、12r、17r、22r、27r,测点角度分别为0°、10°、20°、30°、40°、50°、60°、70°、80°、90°处布置测点,共计50个测点。为方便后续分析,用测点距爆心的距离(爆距)以及测点角度命名测点,爆距7r、12r、17r、22r、27r用数字2~6表示,测点角度0°~90°用数字0~9表示,例如:将爆距为7r、测点角度为0°的测点命名为2-0。
为验证结果的可靠性,对测点3-0、4-0、5-0、6-0测得的冲击波峰值压力进行收敛性分析,得到冲击波峰值压力pmax随网格数量k的变化曲线,如图2所示。当网格数量为3080988时,测点测得的冲击波峰值压力基本收敛。为此,选取网格数量为5324000来研究近海底水下爆炸冲击波的载荷特性。
1.3 数值验证
为验证模型的准确性,选取测点2-0、3-0、4-0、5-0、6-0在无海底底质时测得的峰值压力与Zamyshlyayev经验公式计算结果进行对比,如图3和表5所示。可以看出,数值模拟结果与经验公式拟合结果符合得较好,说明该模型能够用于近海底水下爆炸冲击波研究。
表 5 数值模拟与经验公式结果对比Table 5. Comparison of numerical simulation and empirical formula resultsR/r Peak pressure/MPa Error/% Simulation Empirical formula 7 205.03 193.72 5.84 12 96.32 86.90 10.84 17 60.65 58.62 3.45 22 43.04 43.81 1.75 27 32.58 37.91 14.07 1.4 工况设置
为研究水深(H)及海底底质对近海底水下爆炸冲击波载荷特性的影响,仅改变海底底质及静水压力,其他条件保持不变,进行数值模拟,计算工况如表6所示,其中:工况1和工况2用于研究自由场与近海底冲击波压力的差异,工况2和工况3用于研究底质条件对近海底爆炸冲击波载荷特性的影响,工况2、工况4和工况5用于研究水深对近海底水下爆炸冲击波载荷特性的影响。
表 6 数值模拟工况设置Table 6. Settings of simulation casesCase H/m Seafloor sediment Explosive environment 1 50 Nothing Free field 2 50 Model 1 Near the seabed 3 50 Model 2 Near the seabed 4 100 Model 1 Near the seabed 5 150 Model 1 Near the seabed 2. 结果讨论与分析
2.1 自由场和近海底爆炸冲击波特性对比
图4为工况1和工况2下不同时刻的冲击波压力云图,图5为测点5-3在工况1和工况2下测得的爆炸冲击波压力时程曲线。可以看出,自由场水下爆炸与近海底水下爆炸存在明显差别,受近海底条件的影响,近海底水下爆炸冲击波形状不是规则的圆形,冲击波波面上的压力峰值并非处处相等,而是越靠近海底区域,冲击波压力越小,测点5-3在近海底条件下测得的冲击波峰值压力及后续的脉动压力均小于自由场条件下测得的压力。由于水的冲击阻抗大于海底底质模型1的冲击阻抗,因此,水中冲击波在近海底的反射波为稀疏波,稀疏波追上冲击波波面将会导致冲击波压力降低;另外,海底底质模型1是一种较软且易变形的底质,冲击波作用在海底底质上,不仅会发生反射,还会发生透射,部分能量传递到海底底质中,使海底底质发生形变,冲击波能量也会发生损耗。因此,在近海底条件下,冲击波峰值压力和气泡脉动压力均小于自由场下的压力,理论与数值模拟结果相符。
近海底水下爆炸时,靠近海底区域处的冲击波压力明显偏小,即不同测点角度(α)条件下的冲击波峰值压力存在差异。定义测点测得的冲击波峰值压力与自由场时的冲击波峰值压力之比为该测点的反射系数。计算所有测点的反射系数,结果如图6所示。当爆距比一定时,随着α的增大,反射系数也逐渐增大,即海底吸能现象和稀疏波的影响逐渐减小;α在0°~10°区间时,海底的吸能现象过于剧烈,故不对该角度范围进行深入研究;当α处于20°~30°区间时,测点的反射系数随爆距比增加呈减小趋势,即近海底反射的影响随爆距比的增大而增强;当α > 40°时,海底吸能现象不明显。
为了更细致地研究近海底反射规律,将自由场与近海底工况下测点2-2~测点2-9、测点3-2~测点3-9、测点4-2~测点4-9的冲击波峰值压力进行对比,结果如表7所示,可以看出:当α处于20°~30°时,近海底工况下,冲击波峰值压力是自由场工况下冲击波峰值压力的81%~91%;随着α的增大,测点距海底越来越远,近海底反射现象也越来越弱,近海底反射影响逐渐消失。
表 7 自由场与近海底测得的冲击波峰值压力对比Table 7. Comparison of peak pressures of free-field and near-seabed underwater explosion shock waveR/r α/(°) pmax/MPa Reflection
coefficientR/r α/(°) pmax/MPa Reflection
coefficientNear seabed Free field Near seabed Free field 7 20 179.29 201.44 0.89 12 60 92.73 95.67 0.97 7 30 186.21 203.80 0.91 12 70 95.00 97.86 0.97 7 40 193.43 204.80 0.94 12 80 98.33 93.98 1.05 7 50 194.98 204.80 0.95 12 90 92.42 96.30 0.96 7 60 192.78 203.80 0.95 17 20 47.78 58.76 0.81 7 70 195.10 201.43 0.97 17 30 53.24 59.25 0.90 7 80 202.46 201.74 1.00 17 40 56.52 58.58 0.96 7 90 196.41 195.79 1.00 17 50 58.40 58.58 1.00 12 20 82.82 97.87 0.85 17 60 59.44 59.27 1.00 12 30 86.62 95.67 0.91 17 70 61.07 60.10 1.02 12 40 93.96 94.66 0.99 17 80 60.00 61.17 0.98 12 50 96.23 94.65 1.02 17 90 57.77 59.91 0.96 2.2 海底底质对近海底爆炸冲击波载荷特性的影响
图7显示了工况2和工况3下相同测点测得的近海底冲击波反射系数对比。可以发现,2种底质条件下冲击波传播规律存在明显差异:当底质为模型1时,由于近海底反射波为稀疏波,稀疏波追上冲击波波面将导致冲击波峰值压力变小,且海底底质模型1较软,吸能作用较强,因此,绝大部分测点的反射系数小于1;当海底底质为模型2时,受刚固边界的影响,近海底反射波为压缩波,压缩波追上冲击波波面会导致冲击波峰值压力增大,且海底底质模型2的吸能作用较弱,绝大部分测点的反射系数大于1。无论底质是模型1还是模型2,反射系数异常均发生在爆距比为7~15、α在60°~90°范围内。
为更清楚地研究不同底质条件对反射系数的影响,将不同底质条件下反射系数随测点角度和爆距比的变化绘制成图8。由图8(a)可知,当海底底质为模型1、爆距比保持不变时,反射系数随测点角度的增大呈增大趋势。由图8(c)可知,当海底底质为模型2、爆距比一定时,反射系数随测点角度的增大呈减小趋势。
当海底底质为模型1时,海底底质较软,海底反射稀疏波,对TNT爆炸冲击波的吸收较强,因此,近海底冲击波峰值压力减小,反射系数小于1是常态。图8(a)显示,当爆距比一定时,测点角度越大,反射系数就越大,因而可以推断,海底底质对不同测点角度处反射系数的影响随测点角度的增大而减弱。当海底底质为模型2时,海底底质较硬,海底反射压缩波,对冲击波的吸收很弱,因此,近海底冲击波峰值压力更大,反射系数大于1是常态。图8(c)显示,当爆距比一定时,测点角度越大,反射系数越小,同样可以得到海底底质对不同测点角度处反射系数的影响随测点角度增大而减弱的结论。
从图8(b)可以看出,当测点角度处于20°~40°区间时,不同测点测得的近海底冲击波反射系数随着爆距比的增大而减小;而图8(d)显示,当测点角度处于20°~60°区间时,不同测点测得的近海底冲击波反射系数随着爆距比的增大而增大。由于2种底质对近海底冲击波反射系数的影响截然相反,因此,综合图8(b)和图8(d)可以看出,测点角度在一定范围内时,海底底质对反射系数的影响随着测点爆距比的增大而增强。当测点角度超出该范围时,该现象减弱甚至消失。综合图8(a)~图8(d)可以发现,虽然海底底质发生了变化,不同海底底质对冲击波的影响效果不同,但是测点反射系数受显著影响的区域均集中在20°~50°范围。
为研究近海底水下爆炸反射波对冲击波的影响是否与起爆深度有关,将底质条件为模型1时不同水深处各测点的反射系数进行对比,结果如图9所示。从图9可以看出,随着水深的增加,同一测点处的反射系数基本一致,说明水深并不能显著影响反射系数,即水深的变化并不会对近海底反射现象造成显著影响,也不会对海底底质的吸能作用造成显著影响,近海底反射稀疏波对冲击波的削减作用并不会随着静水压力的变化而发生显著变化。
3. 结 论
基于CEL方法建立了近海底水下爆炸数值模型,对近海底爆炸冲击波的时空分布规律进行了研究,探究了海底底质和水深对近海底爆炸冲击波载荷特性的影响规律,得到如下结论。
(1) 近海底水下爆炸冲击波载荷特性与自由场明显不同。当测点角度为20°~30°时,近海底反射系数为0.81~0.91;随着测点角度进一步增大,近海底反射的影响逐渐减弱;当测点角度达到80°~90°时,近海底反射的影响基本消失。
(2) 在一定的角度范围内,近海底反射的影响随着爆距比的增大而增强,超出该角度范围后,该现象基本消失。改变底质时,近海底反射的影响使冲击波峰值压力增强或减弱,并且影响区域基本一致。水深对近海底水下爆炸反射系数无显著影响。
(3) 海底底质材料属性不同时,其对爆炸冲击波的吸收作用也存在差异,近海底水下爆炸反射波的种类也不一致。当海底底质较硬时,海底底质对冲击波的吸收作用相对较弱,反射波为压缩波,致使测点角度在20°~50°范围内的冲击波峰值压力增大,当测点角度在20°~50°范围内且测点角度一定时,反射系数随爆距比的增大而增大;当海底底质较软时,海底底质对冲击波的吸收作用较强,近海底反射波为稀疏波,致使20°~50°角度范围内冲击波峰值压力减小,当测点角度在20°~50°范围内且测点角度一定时,反射系数随爆距比的增大而减小。海底底质对反射系数的影响区域主要集中在20°~50°测点范围内,超出该范围时,海底底质对反射系数的影响随着测点角度的增大而逐渐消失,海底底质对反射系数的影响随爆距比的增大而增强的现象也随测点角度的增大而逐渐消失。
-
表 1 高温高压处理后回收样品中石墨的质量分数
Table 1. Mass fraction of graphite in the recovered samples after high-temperature and high-pressure treatment
Raw material Pressure environment p/GPa T/℃ wg/% Pure NCD
powderNon-hydrostatic
pressure5 1000 2.6 1300 19.1 1500 22.3 7 1000 3.3 1300 8.0 1500 12.1 9 1000 0 1300 6.8 1500 7.5 NaCl-NCD
powderQuasi-hydrostatic
pressure5 1000 0 1300 3.8 1500 17.6 7 1000 0 1300 0 1500 0 表 2 不同保温时间处理后回收样品中石墨的质量分数
Table 2. Mass fraction of graphite in the recovered samples after different holding time treatment
Raw
materialPressure
environmentp/GPa T/℃ t/min wg/% Raw
materialPressure
environmentp/GPa T/℃ t/min wg/% Pure NCD
powderNon-
hydrostatic
pressure5 1500 0.5 22.3 NaCl-NCD
powderQuasi-
hydrostatic
pressure5 1500 0.5 17.6 1500 5.0 93.8 1500 5.0 58.9 7 1300 0.5 8.0 7 1300 0.5 0 1300 5.0 15.6 1300 5.0 0 1500 0.5 12.1 1500 0.5 0 1500 5.0 17.1 1500 5.0 34.2 -
[1] BOVENKERK H P, BUNDY F P, HALL H T, et al. Preparation of diamond [J]. Nature, 1959, 184(4693): 1094–1098. doi: 10.1038/1841094a0 [2] AKAISHI M, YAMAOKA S, UEDA F, et al. Synthesis of polycrystalline diamond compact with magnesium carbonate and its physical properties [J]. Diamond and Related Materials, 1996, 5(1): 2–7. doi: 10.1016/0925-9635(95)00332-0 [3] BUNDY F P, BOVENKERK H P, STRONG H M, et al. Diamond-graphite equilibrium line from growth and graphitization of diamond [J]. The Journal of Chemical Physics, 1961, 35(2): 383–391. doi: 10.1063/1.1731938 [4] QIAN J, VORONIN G, ZERDA T W, et al. High-pressure, high-temperature sintering of diamond-SiC composites by ball-milled diamond-Si mixtures [J]. Journal of Materials Research, 2002, 17(8): 2153–2160. doi: 10.1557/JMR.2002.0317 [5] HANSEN N. Hall-Petch relation and boundary strengthening [J]. Scripta Materialia, 2004, 51(8): 801–806. doi: 10.1016/j.scriptamat.2004.06.002 [6] ARMSTRONG R W. The influence of polycrystal grain size on several mechanical properties of materials [J]. Metallurgical Transactions, 1970, 1(5): 1169–1176. doi: 10.1007/BF02900227 [7] IRIFUNE T, KURIO A, SAKAMOTO S, et al. Ultrahard polycrystalline diamond from graphite [J]. Nature, 2003, 421(6923): 599–600. [8] BUNDY F P, HALL H T, STRONG H M, et al. Man-made diamonds [J]. Nature, 1955, 176(4471): 51–55. doi: 10.1038/176051a0 [9] NAGENDRA NATH N S. The dynamical theory of the diamond lattice [J]. Proceedings of the Indian Academy of Sciences: Section A, 1935, 2(2): 143–152. doi: 10.1007/BF03035794 [10] PHINNEY F S. Graphitization of diamond [J]. Science, 1954, 120(3114): 393–394. doi: 10.1126/science.120.3114.393 [11] SEAL M. Graphitization of diamond [J]. Nature, 1960, 185(4712): 522–523. doi: 10.1038/185522a0 [12] DAVIES G, EVANS T, DITCHBURN R W. Graphitization of diamond at zero pressure and at a high pressure [J]. Proceedings of the Royal Society of London A: Mathematical and Physical Sciences, 1972, 328(1574): 413–427. doi: 10.1098/rspa.1972.0086 [13] XU N S, CHEN J, DENG S Z. Effect of heat treatment on the properties of nano-diamond under oxygen and argon ambient [J]. Diamond and Related Materials, 2002, 11(2): 249–256. doi: 10.1016/S0925-9635(01)00680-X [14] QIAN J, PANTEA C, HUANG J, et al. Graphitization of diamond powders of different sizes at high pressure-high temperature [J]. Carbon, 2004, 42(12): 2691–2697. [15] GUAN S X, PENG F, LIANG H, et al. Fragmentation and stress diversification in diamond powder under high pressure [J]. Journal of Applied Physics, 2018, 124(21): 215902. doi: 10.1063/1.5051749 [16] LI Q, ZHANG J W, LIU J, et al. Effect of stress state on graphitization behavior of diamond under high pressure and high temperature [J]. Diamond and Related Materials, 2022, 128: 109241. doi: 10.1016/j.diamond.2022.109241 [17] QIAN J, PANTEA C, VORONIN G, et al. Partial graphitization of diamond crystals under high-pressure and high-temperature conditions [J]. Journal of Applied Physics, 2001, 90(3): 1632–1637. doi: 10.1063/1.1382832 [18] 吕世杰, 罗建太, 苏磊, 等. 滑块式六含八超高压实验装置及其压力温度标定 [J]. 物理学报, 2009, 58(10): 6852–6857. doi: 10.3321/j.issn:1000-3290.2009.10.031LYU S J, LUO J T, SU L, et al. A slide-type multianvil ultrahigh pressure apparatus and calibrations of its pressure and temperature [J]. Acta Physica Sinica, 2009, 58(10): 6852–6857. doi: 10.3321/j.issn:1000-3290.2009.10.031 [19] 戴逸, 王文丹, 法志湘, 等. 八面腔压机中一定尺寸的二级压砧上运行的最大组装 [J]. 物理学报, 2021, 70(14): 144702. doi: 10.7498/aps.70.20210006DAI Y, WANG W D, FA Z X, et al. Maximum multianvil cell assembly running on WC anvils with certain size [J]. Acta Physica Sinica, 2021, 70(14): 144702. doi: 10.7498/aps.70.20210006 [20] 黄继武, 李周. 多晶材料X射线衍射 [M]. 北京: 冶金工业出版社, 2012: 95–116. [21] SEAL M. Graphitization and plastic deformation of diamond [J]. Nature, 1958, 182(4645): 1264–1267. doi: 10.1038/1821264a0 [22] LIU J, ZHAN G, WANG Q, et al. Superstrong micro-grained polycrystalline diamond compact through work hardening under high pressure [J]. Applied Physics Letters, 2018, 112(6): 061901. doi: 10.1063/1.5016110 [23] EDMOND J M, PATERSON M S. Strength of solid pressure media and implications for high pressure apparatus [J]. Contributions to Mineralogy and Petrology, 1971, 30(2): 141–160. doi: 10.1007/BF00372255 -