Elastic Properties of ReN2 under High Pressure
-
摘要: 超硬材料在工业上具有广泛的应用前景,如切割器具、研磨材料及耐磨涂层等。作为5d过渡金属双氮化合物之一的ReN2由共价键、离子键及金属键混合而成,因而具有诸多如高硬度、高熔点、耐腐蚀等优异的物理性质,进而具有潜在的研究价值。采用密度泛函理论中的平面波赝势法计算了零温零压下C2/m-ReN2的结构性质,并首次研究了高压下C2/m-ReN2的力学结构稳定性及弹性性质。研究得出了C2/m-ReN2的弹性常数、弹性模量、德拜温度、声速随压强的变化关系,除个别弹性常数,这些物理量皆随压强的增加而增加。还预测了C2/m-ReN2的韧脆性,并估算了C2/m-ReN2的维氏硬度。Abstract: Super hard materials have wide applications in industry, such as cutting tools, abrasive materials, wear resistant coatings. As one of the 5d transition metals double nitrogen compound, ReN2 contains both covalent bond, ionic bond and metallic bonding. In view of its many excellent physical properties, such as high hardness, high melting point and corrosion resistance, ReN2 earns much research interests. This article has calculated the structural properties of C2/m-ReN2 under zero temperature and zero pressure using the plane wave pseudo-potential method of density functional theory, and has studied the mechanical structure stability and elastic properties of C2/m-ReN2 under high pressure for the first time. The relations of the elastic constants, various modulus of elasticity, Debye temperature and the sound speed of C2/m-ReN2 with the pressure have obtained. In addition to the individual elastic constants, these quantities increase with the increase in pressure. In addition, we have also predicted the toughness and brittleness of C2/m-ReN2, and have estimated the Vickers hardness of C2/m-ReN2.
-
Key words:
- ReN2 /
- mechanical structure /
- elastic property /
- hardness
-
随着社会高速发展,特别是我国“双碳目标”的提出,对清洁能源的需求与日俱增,氢能作为一种环保、高效的新能源受到广泛关注。考虑到氢气管道建设成本高,利用现有的天然气管道,向甲烷中掺入一定量的氢气进而提升混合燃料热值,已成为一种经济可行的解决方案[1]。值得关注的是,近年来天然气输运事故屡有发生,若向其掺入氢气,有可能进一步加剧爆炸的危险性。倪靖等[2]讨论了不同掺氢比对甲烷-氧气爆轰特性的影响,发现掺氢后能够提高爆轰波的传播速度和爆轰敏感性。余明高等[3]通过实验研究了障碍物对甲烷-氢气爆炸特性的影响,发现最大爆炸超压和火焰传播速度随着障碍物阻塞率以及氢气体积分数的增大而增大。Yu等[4]研究了掺氢对甲烷-空气预混火焰传播特性的影响,发现随着氢气含量的增加,火焰前沿速度和爆炸超压显著升高。
鉴于可燃气体在运输过程中的危险性,大量学者对可燃气体的抑爆问题开展了研究。Li等[5]发现CO2对甲烷爆炸的抑制效果好于N2。Luo等[6]研究BC粉对氢气-甲烷-空气预混气的抑制后发现,BC粉对掺氢比低的预混气体有较好的抑制作用。陈鹏等[7]发现当金属丝网层数大于3时,甲烷-空气预混气爆炸火焰经过金属丝网时会发生淬熄。徐海顺等[8]通过研究发现铝镁合金泡沫对甲烷-空气预混气体爆炸有较好的抑制作用,材料对传播火焰的影响机制主要体现在湍流促进和冷却抑制两方面。袁必和等[9] 研究了新型多孔聚丙烯复合材料对瓦斯爆炸的影响,发现多孔材料的抑爆性能受填充位置、材料内径以及填充长度的影响。An等[10]通过对比球形非金属材料和铝合金材料对可燃气体的抑爆结果,发现由于球形非金属材料具有较高的结构强度,因而抑爆效果较好。邵继伟等[11]研究多孔材料对可燃气体的抑爆效果后发现,组合型多孔材料在密闭容器管道系统内的抑爆效果更为突出。综合以上研究发现,非金属多孔材料对抑制甲烷爆炸具有良好的效果,但甲烷掺氢的爆炸特性较甲烷出现明显改变,并且前人对非金属多孔材料在甲烷掺氢抑爆方面的研究较少。传统的非金属材料存在易燃、爆炸后熔融物易黏附管道内壁等缺陷,严重限制了多孔非金属材料在可燃气体抑爆中的实际应用。球形抑爆材料中空多孔结构能够有效增大比表面积,扩大热损失,并且耐高温耐火焰,具有良好的阻火隔爆性质[12-13]。
本研究通过实验探究掺氢比对甲烷掺氢爆炸特性的影响,比较单一球形多孔非金属材料和组合球形多孔非金属材料的抑爆效能,并通过改变材料填充长度,分析最佳阻火抑爆性能的搭配参数,为球形多孔非金属材料在混合燃料阻隔防爆领域的应用提供理论支撑与实验依据。
1. 实验装置与方案
1.1 实验装置与流程
为研究掺氢比对甲烷掺氢预混气体爆炸的影响并测试球形多孔材料的抑爆性能,自主搭建了气体管道爆炸实验平台,如图1所示。实验平台主要包括爆炸实验管道、点火系统、配气系统、同步控制、高速摄像系统(拍摄速度4 000帧每秒)以及数据采集系统。
实验管道由3节内径60 mm、管壁厚度10 mm、外径80 mm、长度2 m的水平圆形管道以及1节内径60 mm、长0.5 m、配备可视窗的方形管道组成,各个管道之间通过法兰-螺栓连接,实验管道为封闭状态。实验管道共安装2台压力传感器(CYG508微型高压传感器),分别距离管道左端1.2和6.3 m。点火系统(距离管道右端0.5 m处)由高压点火器以及2根钨丝组成,实验中点火能采用20 J。配气系统通过气瓶连接配气仪进行配气,预混气体中可燃气体的体积分数为10%。通过高速摄像系统采集火焰图像。实验过程如下:首先,清洁管道,连接管道和仪器;然后,检查管道气密性,采用抽真空法配气,向管道中通入预混气体;接着点火触发,数据采集系统对实验数据进行采集;最后,排出废气,重复实验,每组实验重复多次。
1.2 实验材料及实验方案
本实验的主要材料为球形多孔非金属材料和聚氨酯(polyurethane,PU)材料,如图2所示。球形多孔非金属材料的主要材质为聚偏氟乙烯,化学性质稳定,具有很强的抗冲击性和耐高温性,并且球形材料为多孔结构,孔隙率大,比表面积大。组合球形内部填充材料主要为聚氨酯,聚氨酯材料的可塑性较强,对爆燃火焰的传播具有一定的抑制作用[14]。
最大爆炸压力(pmax)是评价爆炸强度以及材料抑爆性能的重要指标。为了判断掺氢比对可燃气体爆炸的影响以及单一球形多孔非金属材料和组合球形非金属材料的抑爆性能,设计如下实验方案:首先研究不同掺氢比x对pmax的影响,分别对掺氢比为0%、5%、10%的甲烷掺氢预混气体进行测试,记录管道内爆炸压力并进行对比分析;在此基础上,研究多孔抑爆材料对pmax的影响,选取pmax最大的预混气体,分别对单一球形材料和组合球形材料进行抑爆实验,填充于距管道左端1.5 m的位置,填充长度分别为20、30、40 cm,对不同填充长度的同一材料以及相同填充长度的不同材料进行对比分析,判断其抑爆性能。
掺氢比x表示为
x=φ(H2)φ(CH4)+φ(H2)×100% (1) 式中:
φ(CH4) 、φ(H2) 分别为甲烷-氢气混合气体中甲烷和氢气的体积分数。2. 实验结果及分析
2.1 掺氢比对爆炸压力的影响
管道中不同掺氢比对预混气体爆炸压力随时间的变化曲线如图3所示。从图3可以看出,1号压力传感器所采集的压力峰值随着氢气体积分数的增大而增大。不同掺氢比对火焰图像的影响如图4所示。对比发现,随着掺氢比的提高,火焰的亮度明显变强,甲烷-空气预混气体掺氢后火焰传播特性及爆炸特性发生显著改变,氢气对甲烷爆炸具有促进作用[15]。
预混气体爆炸压力是反映爆炸破坏效应的重要指标,爆炸压力越大,爆炸造成的后果往往越严重。掺氢比对预混气体最大爆炸压力pmax和最大爆炸压力上升速率(dp/dt)max的影响如图5所示。当掺氢比为0%、5%、10%时,1号压力传感器处的pmax分别为200、235、245 kPa,(dp/dt)max分别为2 250、2 875、3 250 kPa/s。与未掺氢(掺氢比为0%)相比,掺氢比为10%时,1号压力传感器处的pmax提高22.50%,(dp/dt)max提高44.44%,2号传感器处的(dp/dt)max提高26.67%,说明掺氢比较低时,pmax和(dp/dt)max均随着氢气含量的增加而增大,爆炸强度增大。这是由于氢气的活性高于甲烷,反应比甲烷剧烈,发生爆炸反应的时间更短,随着掺氢比的增加,单位体积内氢气的量增多,氢气比例的升高增强了能量释放的集中程度,气体燃烧速率加快,缩短了热量损失时间,导致爆炸压力上升速率提高[16]。
2.2 单一球形材料对爆炸压力的影响
为研究球形非金属材料对甲烷掺氢后预混气体爆炸的抑制作用,采用爆炸压力较大的掺氢比为10%的预混气体进行对照。空管道与填充不同长度的抑爆材料对爆炸压力的影响如图6所示,其中0 cm代表无抑爆材料。由图6可知,填充抑爆材料后,各个工况下的pmax均小于无抑爆材料。1号、2号压力传感器测得无抑爆材料的pmax分别为245、200 kPa。当抑爆材料的填充长度分别为20、30、40 cm时,1号压力传感器处的pmax分别为205、185、175 kPa,对比无抑爆材料时的pmax分别降低了16.33%、24.49%、28.57%,爆炸压力得到有效抑制。这是由于球形抑爆材料是多孔中空塑料球形结构,有较大的比表面积,并且球形材料把管道空间切割成多个小空间,增大了反应面与球形材料表面的接触面积,促进了球形材料与火焰的热交换,增加了热损失,并且冲击波通过球形材料时,被切割的反应面与后面的球形材料发生碰撞导致能量损失,进而使爆炸压力降低[12]。2号压力传感器处各个工况下的pmax相较于无抑爆材料分别降低5.00%、5.00%、12.50%,抑制效果较弱。这是由于随着火焰传播,火焰通过球形抑爆材料后再无障碍物,燃烧速率加快,火焰到达2号压力传感器时已充分反应,致使降压效果不明显。
不同填充长度对(dp/dt)max的影响如图7所示。在低填充长度下,(dp/dt)max有所上升,当填充长度为20 cm时,1号压力传感器处的(dp/dt)max提高7.69%,2号压力传感器处提高36.84%。重复实验均显示(dp/dt)max提高,这是由于当球形材料的填充长度较低时,火焰在穿越过程中的湍流度增大,燃烧传质传热进程加快,燃烧强度增大,火焰传播速度进一步加快[8]。随着填充长度的增加,火焰被分离成多个离散的湍流火焰,最大爆炸压力上升速率逐渐降低,填充长度为40 cm时,相较于无抑爆材料,1号压力传感器处的(dp/dt)max降低33.85%。
2.3 组合球形材料对爆炸压力的影响
不同长度的组合球形材料对爆炸压力的影响曲线如图8所示。从图8可以看出:与无抑爆材料相比,填充组合球形材料后pmax明显降低,说明组合球形材料对管道内可燃气体爆炸有较好的抑制效果。当填充长度分别为20、30、40 cm时,1号压力传感器测得的pmax分别为160、155、120 kPa,与无抑爆材料相比,分别衰减了34.69%、36.73%、51.02%,相较于单一球形非金属材料,抑爆性能分别提升了112.43%、49.98%、78.58%;2号压力传感器测得的pmax为135、135、110 kPa,与无抑爆材料相比,分别衰减了32.50%、32.50%、45.00%。可以看到,当填充长度为20和30 cm时,抑爆性能差别不大,这是多孔材料冷却抑制与障碍物加压共同作用的结果: 当火焰经过球形多孔抑爆材料时,材料会吸收热量从而对爆炸产生抑制效果,材料的填充长度越长,抑制效果越好;球形材料在管道中相当于障碍物,火焰经过障碍物时层流会转变成湍流,增加障碍物数量能够明显增大火焰传播过程中的湍流强度,导致压力上升[3]。当填充长度较低时,材料冷却抑制效果与障碍物加压效果相互抵消,最终导致抑制作用较接近。
从图9中可以看出,(dp/dt)max也得到了较大程度的抑制,与无抑爆材料相比,填充长度为20、30、40 cm时,1号压力传感器处的(dp/dt)max分别衰减34.62%、46.15%、53.85%,2号压力传感器处的(dp/dt)max最高衰减了47.37%,抑制效果明显。与单一球形材料相比,组合球形材料的抑爆性能得到大幅增强。这是由于爆炸穿越填充材料时,被球形材料多孔结构离散后的湍流火焰和冲击波能量与球形多孔材料中的聚氨酯材料接触,聚氨酯材料粗糙的壁面消耗了链式反应中的自由基,阻碍燃烧链式反应的进行,并使部分火焰热量转移到聚氨酯材料中,在火焰与聚氨酯材料发生碰撞的过程中,聚氨酯材料与冲击波产生摩擦,导致部分冲击波能量转化为热量并消散掉[14]。
3. 结 论
基于自主搭建的气体爆炸管道平台,研究了不同掺氢比条件下甲烷掺氢的爆炸特性以及不同类型球形抑爆材料的抑爆性能,得到以下结论。
(1) 氢气对甲烷-空气爆炸具有一定的促进作用。在低掺氢比下,随着可燃气体中氢气体积分数的增加,火焰传播速率加快,最大爆炸压力和最大爆炸压力上升速率升高。
(2) 单一球形多孔非金属材料对可燃气体爆炸具有两方面作用:一方面为抑制作用,体现在填充材料后,球形材料会吸收能量,导致最大爆炸压力降低,并且填充长度越长,抑制效果越好;另一方面为促进作用,体现在填充长度较低时,会导致燃烧传质传热进程加快,火焰传播加速,导致最大爆炸压力上升速率上升。
(3) 组合球形多孔非金属材料与单一球形多孔非金属材料相比,抑爆性能大幅提升,降压效果明显,当填充长度为40 cm时,最大爆炸压力衰减51.02%。球形材料抑爆性能受填充长度影响,在冷却抑制与障碍物增压的共同作用下,填充20和30 cm的组合球形多孔非金属材料的抑爆效果相差不大。
-
表 1 在p=0 GPa和T=0 K下的平衡晶格参数a、b、c及
β ,平衡体积V0,体模量B0,体模量对压强的一阶导数B′0 及其他理论值[13, 27]Table 1. Equilibrium lattice parameters a, b, c and
β , equilibrium volume V0, bulk modulus B0, and its pressure derivationB′0 at p=0 GPa and T=0 K, together with other theoretical results [13, 27]Method a/nm b/nm c/nm β/(°) V0/nm3 B0/GPa B0′/GPa This work 0.682 0.282 0.939 142.38 0.027 59 361 4.78 Ref. [13] 0.682 0.284 0.936 142.30 0.027 60 Ref. [27] 0.683 0.284 0.939 0.027 77 表 2 在p=0 GPa和T=0 K下的体模量B、剪切模量G、杨氏模量E、泊松比
σ 、维氏硬度HV及其他理论值[13, 27]Table 2. Bulk modulus B, shear modulus G, Young’s modulus E, Poisson’s ratio
σ and Vickers hardness HV at p=0 GPa and T=0 K, together with other theoretical results[13, 27]Method B/GPa G/GPa E/GPa σ Hv /GPa This work 370 242 596 0.23 27.66 Ref. [13] 369 217 Ref. [27] 376 210 531 0.26 -
[1] OYAMA S T. Crystal structure and chemical reactivity of transition metal carbides and nitrides [J]. Journal of Solid State Chemistry, 1992, 96(2): 442–445. doi: 10.1016/S0022-4596(05)80279-8 [2] LÉVY F, HONES P, SCHMID P E, et al. Electronic states and mechanical properties in transition metal nitrides [J]. Surface and Coatings Technology, 1999, 120/121: 284–290. doi: 10.1016/S0257-8972(99)00498-3 [3] IVANOVSKII A L. Platinum group metal nitrides and carbides: synthesis, properties and simulation [J]. Russian Chemical Reviews, 2009, 78(4): 303–318. doi: 10.1070/RC2009v078n04ABEH004036 [4] GILMAN J J, CUMBERLAND R W, KANER R B. Design of hard crystals [J]. International Journal of Refractory Metals and Hard Materials, 2006, 24(1/2): 1–5. [5] HAINES J, LEGER J M, BOCQUILLON G. Synthesis and design of superhard materials [J]. Annual Review of Materials Research, 2001, 31(1): 1–23. doi: 10.1146/annurev.matsci.31.1.1 [6] CROWHURST J C, GONCHAROV A F, SADIGH B, et al. Synthesis and characterization of the nitrides of platinum and iridium [J]. Science, 2006, 311(5765): 1275–1278. doi: 10.1126/science.1121813 [7] GREGORYANZ E, SANLOUP C, SOMAYAZULU M, et al. Synthesis and characterization of a binary noble metal nitride [J]. Nature Materials, 2004, 3(5): 294–297. doi: 10.1038/nmat1115 [8] YOUNG A F, SANLOUP C, GREGORYANZ E, et al. Synthesis of novel transition metal nitrides IrN2 and OsN2 [J]. Physical Review Letters, 2006, 96(15): 155501. doi: 10.1103/PhysRevLett.96.155501 [9] ZHAO E, WU Z. Structural, electronic and mechanical properties of ReN2 from first principles [J]. Computational Materials Science, 2008, 44(2): 531–535. doi: 10.1016/j.commatsci.2008.04.016 [10] LI Y, ZENG Z. New potential super-incompressible phase of ReN2 [J]. Chemical Physics Letters, 2009, 474(1/2/3): 93–96. [11] DU X P, WANG Y X, LO V C. Investigation of tetragonal ReN2 and WN2 with high shear moduli from first-principles calculations [J]. Physics Letters A, 2010, 374(25): 2569–2574. doi: 10.1016/j.physleta.2010.04.020 [12] KAWAMURA F, YUSA H, TANIGUCHI T. Synthesis of rhenium nitride crystals with MoS2 structure [J]. Applied Physics Letters, 2012, 100(25): 251910. doi: 10.1063/1.4729586 [13] WANG Y, YAO T, YAO J L, et al. Does the real ReN2 have the MoS2 structure? [J]. Physical Chemistry Chemical Physics, 2013, 15(1): 183–187. doi: 10.1039/C2CP43010J [14] WANG Y, LV J, ZHU L, et al. Crystal structure prediction via particle-swarm optimization [J]. Physical Review B, 2010, 82(9): 094116. doi: 10.1103/PhysRevB.82.094116 [15] BOUHEMADOU A, KHENATA R. Pseudo-potential calculations of structural and elastic properties of spinel oxides ZnX2O4 (X= Al, Ga, In) under pressure effect [J]. Physics Letters A, 2006, 360(2): 339–343. doi: 10.1016/j.physleta.2006.08.008 [16] LOUAIL L, MAOUCHE D, ROUMILI A, et al. Calculation of elastic constants of 4d transition metals [J]. Materials Letters, 2004, 58(24): 2975–2978. doi: 10.1016/j.matlet.2004.04.033 [17] SEGALL M D, LINDAN P J D, PROBERT M J, et al. First-principles simulation:ideas, illustrations and the CASTEP code [J]. Journal of Physics: Condensed Matter, 2002, 14(11): 2717–2744. doi: 10.1088/0953-8984/14/11/301 [18] VANDERBILT D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism [J]. Physical Review B, 1990, 41(11): 7892–7895. doi: 10.1103/PhysRevB.41.7892 [19] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letters, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865 [20] PERDEW J P, CHEVARY J A, VOSKO S H, et al. Atoms, molecules, solids, and surfaces:Applications of the generalized gradient approximation for exchange and correlation [J]. Physical Review B, 1992, 46(11): 6671–6687. doi: 10.1103/PhysRevB.46.6671 [21] PFROMMER B G, COTE M, LOUIE S G, et al. Relaxation of crystals with the quasi-Newton method [J]. Journal of Computational Physics, 1997, 131(1): 233–240. doi: 10.1006/jcph.1996.5612 [22] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188–5192. doi: 10.1103/PhysRevB.13.5188 [23] WALLACE D C. Thermodynamics of crystals [M]. New York: Wiley, 1972: 582. [24] WANG J, LI J, YIP S, et al. Mechanical instabilities of homogeneous crystals [J]. Physical Review B, 1995, 52(17): 12627–12635. doi: 10.1103/PhysRevB.52.12627 [25] BARRON T H K, KLEIN M L. Second-order elastic constants of a solid under stress [J]. Proceedings of the Physical Society, 1965, 85(3): 523–532. doi: 10.1088/0370-1328/85/3/313 [26] BIRCH F. Finite elastic strain of cubic crystals [J]. Physical Review, 1947, 71(11): 809–824. doi: 10.1103/PhysRev.71.809 [27] ZHAO Z L, BAO K, LI D, et al. Nitrogen concentration driving the hardness of rhenium nitrides [J]. Scientific Reports, 2014, 4(1): 4797. [28] WATT J P. Hashin-Shtrikman bounds on the effective elastic moduli of polycrystals with orthorhombic symmetry [J]. Journal of Applied Physics, 1979, 50(10): 6290–6295. doi: 10.1063/1.325768 [29] HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society, Section A, 1952, 65(5): 349. doi: 10.1088/0370-1298/65/5/307 [30] PUGH S F.XCII Relations between the elastic moduli and the plastic properties of polycrystalline pure metals [J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1954, 45(367): 823–843. doi: 10.1080/14786440808520496 [31] FRANTSEVICH I N, VORONOV F F, BOKUTA S A. Elastic constants and elastic moduli of metals and insulators [M]. Kiev: Naukova Dumka, 1983: 60–180. [32] IVANOVSKII A L. Microhardness of compounds of rhenium with boron, carbon, and nitrogen [J]. Journal of Superhard Materials, 2012, 34(2): 75–80. doi: 10.3103/S1063457612020013 [33] TIAN Y J, XU B, ZHAO Z S. Microscopic theory of hardness and design of novel superhard crystals [J]. International Journal of Refractory Metals and Hard Materials, 2012, 33: 93–106. doi: 10.1016/j.ijrmhm.2012.02.021 -