Thermodynamic analysis and kinetics mechanism for direct nitridation reaction
-
摘要: 用TG-DSC热重同步热实验分析的方法,研究不同升温速率下,硅粉氮化机理及化学反应动力学.发现温度在1 000~1 300 ℃时:差示扫描量热曲线各出现一个吸热、一个放热峰,说明氮化机理已发生改变.在1 000~1 100 ℃温度范围内,氮化硅转化率显著增加,即温度是影响其转化率的主要因素.实验表明:氮化反应的限制性环节由反应开始阶段的界面化学反应控制和之后的界面化学反应与内扩散混合控制组成;通过动力学计算得到表观活化能E=404.5 kJ/mol,频率因子A=9.57×1015 m/s,反应级数n=0.95,最终得到反应的速率方程的数学表达式.Abstract: The mechanism and chemical kinetics of nitridation reaction at different heating rates was investigated by thermogravimetry. It is found that there is an exothermic peak and an endothermic peak in the differential scanning calorimetry curve at 1 000~1 300 ℃ centigrade, which indicates that the mechanism has changed. The conversion rate of silicon nitride is significantly increased within 1 000 to 1 100 ℃, and it indicates temperature is the dominant factor affecting the conversion rate. Experiments results show that the restrictive step of nitridation reaction is composed of chemical reaction in early stage of reaction and the mixed-control of chemical reaction and intraparticle diffusion in the subsequent reaction; the apparent activation energy is kinetically calculated to be 404.5 kJ/mol, the frequency factor A=9.57×1015 min/s, the reaction order n=0.95, and mathematical expression of reaction rate equation is finally obtained.
-
Keywords:
- silicon nitride /
- thermogravimetric analysis /
- thermodynamics /
- kinetic modeling /
- synthesis
-
0 矿井通风系统可靠性分析的特点
矿井通风系统可靠性的高低对矿井安全生产有着直接的影响.国内外自20世纪80年代以来,已取得了初步的研究成果[1-4].但由于问题本身的复杂性,还有许多问题尚未解决.矿井通风系统可靠性研究应在充分考虑以下几个特点的基础上进行.
(1)矿井通风系统是由风机、分支和通风构筑物三者组成的有机整体,在讨论矿井通风系统可靠性时,应充分考虑各单元的特点以及各单元之间的相互联系和相互制约;
(2)各通风构筑物单元在使用过程中的可靠性是相对的、动态变化的,其在某一时刻的可靠度不仅依赖于当时构筑物所处的状态,也与预计的剩余服役期有关;
(3)各分支的风量在系统结构不变时,不仅与风阻和风压有关,而且与风机和通风构筑物当时所处的状态有关.
1 分支的可靠性
分支可靠性定义为:在矿井通风网络结构不变的条件下,在规定的时间内某一分支i的风量值能够保持在一个合理的范围之内的概率,称为该分支的可靠度,记为Ri.并相应地称该分支在任意时刻能保持其风量值在一个合理的范围之内的概率为该分支的有效度,记为Ai.
以Qi(t)表示实际情况下第i分支在t时刻的风量值,Qi(0)表示在风机和通风构筑物都稳定可靠的前提下第i分支的风量值.根据安全规程规定或分支用风需求确定两个风量限定值Wi(1)和Wi(2),当第i分支满足
(1) 则第i分支正常,否则第i分支为不正常.因此第i分支的有效度Ai可表示为
(2) 通常,矿井通风网络的分支可靠度不仅是指风流的稳定程度或抗干扰能力,它与下列因素密切相关:①风流稳定性,角联网络分支造成的风流反向的可能性,风机的不稳定运行,井下提升运输设备因素等对风流运动的瞬时扰动现象;②调节设施工作的可靠性以及失效状态可能产生的危害程度和范围;③分支中的有毒性、窒息性和爆炸性气体与粉尘含量,一旦超限将使局部网络暂时失效;④漏风情况,风流的质量状况以及气候条件等.
由文献[5-8]可推导出,Qi(t)保持在闭区间[Qi(0)-Wi(1),Qi(0)+Wi(2)]内的概率(即第i分支的可靠度)为
(3) (4) (5) Ri值即为在综合考虑了分支风阻的随机扰动及通风构筑物的可靠性等因素后,第i条分支的风量保持在[Qi(0)-Wi(1),Qi(0)+Wi(2)]范围内的概率,即第i条分支的可靠度.
2 风机的可靠性
可用一般可修系统的可靠性计算方法来评价风机的可靠性[9-12].
对于单台风机若其寿命遵从指数分布1-e-λt,t≥0,故障后的修理时间也遵从指数分布1-e-μt,t≥0.当两个随机变量相互独立,且故障部件修复后的寿命分布与新部件相同时,单台风机的有效度为
(6) 对于主扇,根据矿山安全规程的要求,都需要设备用风机(一般只有电机备用).可将主扇单元看作用两个相同型号部件组成的冷储备系统,若这两台风机寿命和故障后修理时间均服从负指数分布,则主扇单元的有效度为
(7) 一般情况下,若某矿备有n台相同型号的主扇,正常情况下需要k台主扇工作,当工作的某台主扇发生故障后,则将剩余的n-k台中的某台启动,此时可将主扇单元视为k/n(G)表决系统,若每台主扇的寿命分布和故障后的修理时间分布均服从负指数分布,则主扇单元的有效度为
(8) 3 通风构筑物的可靠性
通风构筑物的类别不同,其可靠度的定义和计算方法亦不同.
3.1 截断风流的通风构筑物可靠性
对于风门、挡风墙、风帘截断风流的通风构筑物,其主要功能就是将风流截断,阻止风流通过.其可靠度定义为:对于截断风流的通风构筑物,在某一时刻t完成截断风流的能力,称为该构筑物在t时刻的可靠度.若Rs1(t)表示这类构筑物在t时刻的可靠度,则有
(9) 式(9)中,Qi(t)为截断风流的通风构筑物在t时刻的漏风量,m3/s;QSj(t)为截断风流的通风构筑物入风侧的设计风量,m3/s.其值为网络解算所得的风量值,此值不随构筑物的服役期t的变化而变化,而只与t时通风系统的网络结构有关.
新建截断风流的通风构筑物总有一定的漏风量,其值可通过实测而得,以Qi(0)表示.则新建截断风流的通风构筑物的可靠度为
(10) 该通风构筑物在使用过程中,随着其服役期t的增长,由于不断受到外部环境及构筑物本身的逐渐老化等因素的影响,在不维修的情况下,其漏风量会逐渐增加.在t时构筑物的漏风量Qi(t)可以按式(11)评估
(11) 式(11)中,Ci(t)是与构筑物在服役期内的经历有关的系数.
3.2 通过风流的通风构筑物可靠性
风桥等通过风流的通风构筑物的功能就是要使风流顺利地通过.可将其可靠度定义为:对于通过风流的通风构筑物在某一时刻t,其能使风流顺利通过的能力称为该构筑物在t时刻的可靠度.若以RS2(t)表示这类构筑物在t时刻的可靠度,则有
(12) 式(12)中,Qi(t)为通过风流的通风构筑物入风侧的风量,m3/s,即为网络解算所得之风量值;Q0(t)为该通风构筑物出风侧的风量值,m3/s.
新建构筑物在风流通过的过程中总有一定的损失,即有
(13) 即其风量损失为
(14) 在不维修的情况下,随着构筑物服役期t的增长,QS(t)逐渐增大.QS(t)可以实测,也可以按式(15)评估
(15) 式(15)中,CS(t)是与构筑物在使用期内的经历等因素有关的系数.
3.3 调节风流的通风构筑物可靠性
对于风窗等调节风流的通风构筑物,其可靠度定义为:对于调节风流的通风构筑物在某时刻t,其能通过的风量符合设计要求的程度,称为该构筑物的可靠度.并以RS3(t)表示.据此定义有
(16) 式(16)中,Qp(t)为t时刻该类构筑物的实际过风量,m3/s;Qd为设计构筑物的过风量,m3/s.Qp(t)值可以实测获得,也可利用与式(15)相似的公式评估而得.
4 矿井通风网络的可靠性
对于某一矿井通风网络,以m表示其节点个数,P表示其通路矩阵,A表示其关联矩阵,I为单位矩阵.若令C=A+I,取矩阵C中元素Cm1的余子式,并采用逻辑乘和逻辑加的运算方法对这一余子式展开,则得到从起始节点(进风井口)到终止节点(排风井口)的全部有向通路.
若以pij描述分支ej是否在第i通路上,并令
(17) 则可求得该矿井通风网络的通路矩阵为
(18) 式(18)中,k为该网络的通路数;n为该网络的分支数.
若以Rnet表示矿井通风网络的可靠度,则有
(19) 式(19)中,pi表示矿井通风网络的第i条通路.
式(19)是网络可靠度的计算式,它与网络中进风点和出风点间的最小通路密切相关.但由于最小通路间的不相斥性,具有大量公共边,使得网络可靠度的计算相当困难.为了改进算法,采用不交型布尔代数算法.通过不交型布尔代数的化简、展开吸收、归并等运算,使得最小通路间不含有公共边.这样对于通路,可按互不相容事件求概率,得到网络可靠度.
5 矿井通风系统可靠度优化分配
采用一个可行的状态S=(R,Q,H),使得S在满足按需分风和经济性的条件下,可靠度达到最大,这个问题可用下列非线性规划进行描述[13].
(20) 式(20)中,S1和S2分别为固定风量分支集和风机分支集;bij为该网络回路矩阵B中的元素;Ei为第i台风机的目标流线分支集;Hfi为第i台风机的压力上限值.
式(20)的目标函数比较复杂,直接求解困难较大,应设法把目标函数进行简化.由式(4)和式(5)不难看出,要使得每个Ri都尽可能的大,只有让均值μ*i尽可能地接近Qi(0),让方差σ*i尽可能小,而由μj(1)、μj(2)、σj(1)和σj(2)的不确定性,也只有让每个αij和βij的绝对值尽可能地小才行.根据已有研究结果[7-8],矿井通风网络的阻力定律、风量平衡定律和风压平衡定律可以综合写成下面的方程组.
(21) 由隐函数存在定理知,如矿井通风网络处于平衡状态S0的Jacobi行列式不为零,即
(22) 则式(21)在S0的某一邻域内存在唯一的连续偏导数解[14].
(23) 由此可知,当J的绝对值充分大时,可使每个αij和βij的绝对值尽可能地小.式(20)的目标函数可以改为对J进行极大化.各分支风量一定时,调节各分支的风阻极大化J的优化模型如下:设已知各分支的风量为Q1(0),Q2(0),…,Qn(0),至少存在一棵生成树,使其所有余树枝风量不为零.从经济方面考虑,约定风机的压力上限分别为Hf1,Hf2,…,HfF,则这时的优化模型为
(24) 式(24)中,ΔRj为分支j的风阻增量;J(ΔR1,ΔR2,…,ΔRn)为Jacobi行列式;S3为不可调分支[15]集;Ei和S2同前,M=106,且
(25) (26) (27) (28) (29) 6 应用
图 1为某矿井通风网络,10条分支,6个节点.风机工作方式为抽出式,型号为4-72-11№20,转速为315r/min,风流上行流动.各分支对应的风阻值(N·s2/m8)为式(30),该风机的H-Q特性方程为式(31),分支4、8、9为固定风量分支,该网络无火灾时的风量分配为式(32)
(30) (31) (32) 现分析火灾时各分支的风流稳定性.选一棵生成树,如图 2所示(实线部分).取Wi(1)=0.2Qi(0),Wi(2)=0.3Qi(0),μi(1)=μi(2)=0,σi(1)=0.04,σi(2)=1,i=1,2,…,10.上机计算求得各分支可靠度
(33) 图 2 图 1网络的生成树由可靠度知分支5风流最不稳定,分支9风流最稳定;固定风量分支中,分支4最不稳定,分支9最稳定.
当分支7发生火灾时,主干分支中的风量增加,旁侧分支中的风量减小甚至反向.与主干分支相交的分支4、6、8、9中,由可靠度可知分支4风流首先反向,分支6次之,分支8再次之,这与分支可靠度是对应的.解算结果证明由分支可靠度来判断风流稳定性是正确的.尽管分支5没有反向,但它的风流稳定程度受干扰最大,例如当Δt=1300℃时,分支5中风量变化率为401.40%,其他分支中风量变化率值均小于它.实际上,风流稳定性是风向和风量稳定性的总称.
-
表 1 硅粉的化学组成
Table 1 Chemical composition of silicon powder
成分 Si Al Ca Fe 含量/% ≥99.999 ≤0.001 ≤0.000 1 ≤0.000 1 表 2 高纯氮气的化学组成
Table 2 Chemical constituents of high-purity nitrogen
成分 N2 O2 H2O H2 总碳氢化合物 含量/% ≥99.999 ≤1.0×10-6 ≤1.0×10-6 ≤0.5×10-6 ≤1.0×10-6 表 3 不同升温速率下硅粉氮化反应的DSC 峰温数值
Table 3 DSC peak temperature values of nitriding reaction of silicon powder at different heating rates
β/(K·min-1) Tm /K 5 1 395.94 10 1 416.64 15 1 432.16 20 1 439.06 -
[1] SHUKLA P P, LAWRENCE J. Fracture toughness modification by using a fibre laser treatment of a silicon nitride engineering ceramic[J]. Journal of Materials Science, 2010, 45(23): 6540-6541. doi: 10.1007/s10853-010-4743-6
[2] 鲁元,杨建锋,李京龙. 碳热还原-反应烧结法制备多孔氮化硅陶瓷[J]. 无机材料学报, 2013, 28(5): 469-473. http://www.baidu.com/link?url=GZ4-5b11H224VJz2TEhTN8QlI9bjSEyAbeRy4oPH38SuYaD-vUYRRpEsQEsxzQSx55mvAlgKHeBXMAPqQ2fIL64NVDEqjydIpR-EfMOwSOWP0aqh-jtj-fpQupfkROAHnUgHv8LF996e4G3HWJHPNKtUF9nMdGrV8xF5xydv3yeB-sKUfjsi1ZuYN8NsFE0FVoFPaTJqIjnhxJShvkmo5ooKOknWO33IHhTCaOyf5cshb1s-r76B0BtTRPiORf7jzvHrFV6IDgGCsb4X6ziXetybGa02rXlEZP1mUQVuTKKHWmbhzn06VMI9vlGZPwI0G1p3WxWMooBkvOdFzCxHeU88iLqUtR7Q_0bYF60lowHwO5t3aJsVh9f6dgvuDGtC20RBR8tRW6EaqAMBC1uHi7PhZNOeMYH_t_WG90P4WbYOklQEZ76xx3znopLFcLEBH0rg0uuwzjSkbZ4mm9VpGK&wd=&eqid=86998a3500032b1a00000005581019e0 [3] 李金富,李康,李拥军,等. 工艺参数对自蔓燃制备氮化硅粉体的影响[J]. 硅酸盐通报,2007, 26(2): 252-255. http://www.cnki.com.cn/Article/CJFDTOTAL-GSYT200702007.htm [4] 李亚伟,张忻,田海兵,等. 硅粉直接氮化反应合成氮化硅研究[J]. 硅酸盐通报, 2003, 22(1): 30-34. http://www.cnki.com.cn/Article/CJFDTOTAL-GSYT200301007.htm [5] 王勇,沃银花,姚奎鸿,等. 流态床CVD 法纳米氮化硅粉体的制备[J]. 无机材料学报, 2006, 21(1): 41-45. http://www.cnki.com.cn/Article/CJFDTOTAL-WGCL200601007.htm [6] 李虹,黄莉萍,蒋薪. 碳热还原法制备氮化硅粉体的反应过程分析[J]. 无机材料学报, 1996, 11(2): 241-246. http://www.cnki.com.cn/Article/CJFDTOTAL-WGCL602.007.htm [7] 万小涵,张广清,JOHN SHARP,等. 高氮分压对碳热还原/氮化法合成氮化硅的影响[J]. 有色金属工程,2015, 5(4): 9-12. http://www.cnki.com.cn/Article/CJFDTOTAL-YOUS201504003.htm [8] 马啸尘,尹洪峰,张军战,等. 以木屑为碳源制备氮化硅粉体的研究[J]. 耐火材料, 2015, 49(1): 31-35. http://www.cnki.com.cn/Article/CJFDTOTAL-LOCL201501009.htm [9] 古亚军,曹迎楠,李发亮,等. 铁纳米颗粒催化氮化硅粉[J]. 硅酸盐学报, 2014, 42(12): 1585-1589. http://www.cnki.com.cn/Article/CJFDTOTAL-GXYB201412019.htm [10] 李勇,朱晓燕,王佳平,等. 反应烧结氮化硅–碳化硅复合材料的氮化机理[J]. 硅酸盐学报,2001, 39(3): 447-451. http://www.cnki.com.cn/Article/CJFDTOTAL-GXYB201103014.htm [11] 高梅,李勇,秦海霞,等. 闪速燃烧合成氮化硅铁的氮化机理[J]. 硅酸盐学报, 2015, 43(3): 358-362. http://www.cnki.com.cn/Article/CJFDTOTAL-GXYB201503019.htm [12] 李绍芬. 反应工程[M]. 3版. 北京: 化学工业出版社, 2013. [13] 郭汉杰. 冶金物理化学[M]. 北京: 高等教育出版社, 2006. [14] 杨福明,王立,尹少武,等. 硅粉常压直接氮化过程的非催化气固反应模型[J]. 北京科技大学学报,2013, 35(6): 785-792. http://www.cnki.com.cn/Article/CJFDTOTAL-BJKD201306012.htm [15] KISSINGER H E. Reaction kinetics in differential thermal analysis[J]. Analytic Chemistry, 1957, 2(11): 1702-1704.
[16] 沈兴. 差热、热重分析与非等温固相反应动力学[M]. 北京: 冶金工业出版社, 1995. [17] STAVA V, SETAK J. Computer calculation of the mechanism and associated kinetic data using a non-isothermal integral method[J]. Journal of Thermal Analysis and Calorimetry, 1975, 8(3): 477-489. doi: 10.1007/BF01910127