基于吸脫附理論的RDX驅(qū)酸洗滌動力學(xué)模型研究

謝 強,陳厚和,李 慧

(1.安徽理工大學(xué) 化學(xué)工程學(xué)院,安徽 淮南 232001;2.南京理工大學(xué) 化學(xué)與化工學(xué)院,江蘇 南京 210094)

引 言

黑索今(RDX)是一種應(yīng)用廣泛的單質(zhì)高能炸藥[1-3],其生產(chǎn)能力與本質(zhì)安全化程度直接影響著其應(yīng)用。目前常用的RDX制造方法為直接硝解法與醋酸酐法。在RDX的直接硝解法生產(chǎn)制造過程中,由于硝化液的氧化過程與RDX的結(jié)晶過程同時進(jìn)行,極易在結(jié)晶過程中于RDX晶體周圍包裹著含酸的氣泡[4-5]。此外,有學(xué)者[5-6]對RDX等一類硝胺炸藥的晶體形貌進(jìn)行表征發(fā)現(xiàn),其晶體存在著孔洞、凹坑麻點、溝槽、裂紋等缺陷,而且在這些缺陷中均充滿著硝酸的母液[7-8]。由此綜合導(dǎo)致了直接硝解法制造的RDX晶體含酸量較高。酸性粒子催化使得RDX貯存的危險性有所增加[9]。因此,研究RDX的驅(qū)酸洗滌動力學(xué)可為RDX驅(qū)酸洗滌的工藝參數(shù)優(yōu)化與選擇提供理論依據(jù)。

目前針對RDX驅(qū)酸洗滌過程的研究主要集中在機理、數(shù)值模擬[10]、動力學(xué)等方面。如張幺玄等[11]對超聲波驅(qū)酸洗滌RDX的機理進(jìn)行了探討,認(rèn)為超聲波清洗器內(nèi)的空化氣泡在RDX界面上發(fā)生了非對稱性潰滅,從而產(chǎn)生局部高速射流與渦流作用,強化了RDX驅(qū)酸清洗效果。2014年羅志龍等[12]基于RDX的超聲驅(qū)酸洗滌數(shù)據(jù),假設(shè)驅(qū)酸過程符合一級動力學(xué)規(guī)律,采用最小二乘法進(jìn)行了線性擬合,獲得了相關(guān)動力學(xué)參數(shù)。在此基礎(chǔ)上,2016年李婷婷等[13]對比分析了磁力常規(guī)清洗與超聲波清洗的動力學(xué)參數(shù),常規(guī)清洗的表觀活化能大于超聲波清洗的表觀活化能,表明超聲洗滌更容易驅(qū)酸。然而,對其驅(qū)酸洗滌的動力學(xué)研究僅停留于假設(shè)的一級反應(yīng),由擬合效果來看,動力學(xué)模型描述的準(zhǔn)確性不夠高。

基于此,本研究以吸附、擴散等理論為出發(fā)點,基于相應(yīng)理論的動力學(xué)模型進(jìn)行RDX驅(qū)酸洗滌的動力學(xué)探討,以獲得較為精確的動力學(xué)模型來較好地描述RDX驅(qū)酸洗滌過程,以期為RDX制造工藝參數(shù)優(yōu)化提供較為準(zhǔn)確的理論依據(jù),為相關(guān)的工藝改進(jìn)與連續(xù)化設(shè)計提供基礎(chǔ)。

1 驅(qū)酸洗滌動力學(xué)模型的建立

基于物理吸附與化學(xué)吸附的定義可知,由于RDX晶體表面存在著孔洞、凹坑麻點、溝槽、裂紋等缺陷且其晶體存在著非球形的菱與角,在含有硝酸的水溶液中RDX表面存在著不飽和力場。因此RDX表面自發(fā)地吸附硝酸與水。此外,因RDX與水、硝酸不存在化學(xué)反應(yīng),無任何形式的化學(xué)鍵,所以其吸附過程應(yīng)為物理吸附,吸附力為RDX表面分子與硝酸分子、水分子間的分子間作用力。此外,陳剛[14]在對RDX表面吸附溶劑分子的研究過程中發(fā)現(xiàn),溶劑分子與RDX晶面之間的吸附能主要是由靜電作用能和范德華作用能組成,且溶劑分子對RDX晶面的吸附是一個放熱過程。由于硝酸水溶液可作為RDX的溶劑,此研究結(jié)果同樣說明了RDX與硝酸水溶液之間的吸附不存在任何的化學(xué)變化。

由此可以推測,RDX表面酸脫除機制主要由兩個步驟組成:(1)RDX表面吸附的酸在吸收一定能量(能量大于吸附能)后從固體表面脫離;(2)脫離的酸由RDX晶體附近借助一定濃度梯度擴散遷移至溶液中。因此,在驅(qū)酸洗滌過程中,速率控制步驟有可能是其中一種或是兩種的綜合。

目前,國內(nèi)外對固液吸附動力學(xué)模型與機理的研究有很多。在探索與實踐的歷程中,逐漸形成了6種典型且常用的動力學(xué)模型[15-17],包括基于溶液濃度的擬一階與擬二階動力學(xué)模型、Elovich方程模型、Freundlich方程模型、Avrami方程模型、H-S分?jǐn)?shù)階動力學(xué)模型。結(jié)合脫附與吸附之間的可逆關(guān)系,此6種常用的脫附動力學(xué)模型如表1所示。

表1 常用的脫附動力學(xué)模型Table 1 Common desorption kinetic models

目前對質(zhì)量擴散的動力學(xué)機理及其函數(shù)已有很多,如表2所示[18]。其擴散動力學(xué)模型為:

表2 常用的質(zhì)量擴散動力學(xué)機理函數(shù)Table 2 Commonmass diffusion kinetic mechanism functions

(1)

(2)

式中:α為轉(zhuǎn)化率,即RDX表面酸的脫除率,見下式:

α=(qmax-qi)/(qmax-qe)

(3)

式中:qe為RDX表面吸附與脫附達(dá)到動態(tài)平衡時的酸度。

此外,不論是以脫附模型為主、還是以擴散模型為主,假設(shè)RDX的驅(qū)酸洗滌速率常數(shù)與溫度的關(guān)系服從阿倫尼烏斯方程,即:

k=Aexp(-E/RT)

(4)

式中:E為RDX表面脫酸的活化能。

在上述動力學(xué)模型進(jìn)行非線性擬合的同時,需建立評選準(zhǔn)則以選擇較優(yōu)的RDX表面脫附動力學(xué)模型。因此采用全局優(yōu)化的通用算法(UGO)將實驗數(shù)據(jù)與脫附動力學(xué)模型進(jìn)行擬合。相關(guān)系數(shù)是選擇的首要衡量標(biāo)準(zhǔn)。其次的衡量標(biāo)準(zhǔn)為均方根誤差(RMSE)、確定系數(shù)(DC)和殘差平方和(RSS)。相關(guān)系數(shù)與DC越大,RMSE與RSS越小,則說明模型擬合程度越高。這些參數(shù)可由式(5)～(7)獲得:

(5)

(6)

(7)

式中:qexp,i為第i個實測轉(zhuǎn)化率值;qpre,i為第i個預(yù)測的轉(zhuǎn)化率值;N為觀測的樣本數(shù)。

2 實 驗

2.1 實驗材料

RDX粗品,平均粒徑1.3mm,酸度約1.072%,自制;氫氧化鈉,分析純,國藥集團(tuán)化學(xué)試劑有限公司;酚酞,分析純,上海三愛思試劑有限公司。

2.2 實驗裝置

經(jīng)過實驗室直接硝解法合成RDX的步驟與抽濾后,得到RDX粗品。然后將裝有一定量的RDX粗品與去離子水的燒杯放置于集熱式磁力攪拌加熱器與超聲波清洗機內(nèi),在不同的洗滌溫度下分別進(jìn)行不同清洗時間的實驗。隨后將驅(qū)酸洗滌實驗后的RDX經(jīng)抽濾后得到含酸濾餅。最后,將濾餅均勻分散于去離子水中,滴加酚酞指示劑,并使用裝有氫氧化鈉標(biāo)準(zhǔn)溶液的微量滴定管進(jìn)行酸度測定。濾餅酸度計算公式為:

(8)

式中:q為濾餅酸度(%);c為氫氧化鈉標(biāo)準(zhǔn)溶液的濃度(mol/L);V1為空白實驗所需氫氧化鈉標(biāo)準(zhǔn)溶液的體積(L);V2為滴定濾餅所需氫氧化鈉標(biāo)準(zhǔn)溶液的體積(L);m為濾餅質(zhì)量(g)。

實驗流程與裝置如圖1所示。

圖1 RDX驅(qū)酸洗滌動力學(xué)實驗流程與裝置Fig.1 The kinetic test process and device of washing RDX

2.3 實驗過程

為了研究RDX驅(qū)酸洗滌動力學(xué),在20、25、30、35、40℃的洗滌溫度下分別進(jìn)行不同清洗時間的實驗(超聲清洗時間分別為10、20、40、60、80min,加熱攪拌清洗時間分別為10、40、70、100、130min)。為了與李婷婷等[13]研究的結(jié)果進(jìn)行對比分析,其他實驗條件均為一致:兩種清洗方式的液固比均為2∶1,超聲清洗頻率為25kHz,超聲功率為85%,攪拌速率為150r/min。然后基于不同溫度下不同清洗時間的酸度值及6個脫附模型和9個擴散模型,采用共享參數(shù)模式下非線性擬合方法進(jìn)行驅(qū)酸洗滌動力學(xué)參數(shù)的求解。其中,針對脫附模型中可轉(zhuǎn)換為積分形式的模型1、2、5、6來說,采用共享參數(shù)模式下非線性擬合方法進(jìn)行求解;針對不可轉(zhuǎn)換積分形式的模型3來說,可采用共享參數(shù)模式下常微分方程參數(shù)擬合方法進(jìn)行求解;針對不可轉(zhuǎn)換積分形式的模型4來說,微分方程求解器無法進(jìn)行計算。此外,針對9個擴散模型來說,基于不同機理函數(shù)的積分形式將實測數(shù)據(jù)進(jìn)行轉(zhuǎn)化,隨后采用共享參數(shù)模式下非線性擬合方法進(jìn)行求解。

共享參數(shù)模式下非線性擬合方法,即在具有共同的擬合參數(shù)且含有隨不同數(shù)據(jù)文件而變化的常數(shù)情況下采用的一種非線性擬合方法,可通過1stOpt軟件實現(xiàn)。以脫附模型6與磁力攪拌加熱清洗方式獲得的5種溫度條件下的數(shù)據(jù)為例進(jìn)行說明,其共有4個擬合參數(shù)(脫附活化能E、指前因子A、脫附時間指數(shù)n、脫附模型常數(shù)m),洗滌溫度(20、25、30、35、40℃)即為隨不同數(shù)據(jù)而變化的常數(shù)。

3 結(jié)果與討論

3.1 脫附模型的優(yōu)選

不同溫度下兩種清洗方式的脫附模型擬合狀態(tài)參數(shù)間的對比如圖2所示?；赗與DC參數(shù)均為最大且RMSE與RSS均為最小的選擇原則,黑色柱形則表示最優(yōu)的模型。擬合較優(yōu)模型的動力學(xué)參數(shù)見表3。

圖2 不同溫度下兩種清洗方式的脫附模型擬合狀態(tài)參數(shù)對比Fig.2 Comparison of fitting parameters of desorption models for two washing methods at different temperatures

表3 兩種清洗方式的較優(yōu)脫附模型擬合參數(shù)Table 3 The better fitting parameters of desorption models for two washing methods

由圖2可知,在20、25、30、35、40℃的清洗溫度下,磁力攪拌清洗的最優(yōu)脫附模型分別為模型6、6、5、6、6,而超聲波清洗的最優(yōu)脫附模型分別為模型1、1、1、6、6。此外,還可以看出,對磁力攪拌清洗來說,模型5、6的擬合狀態(tài)參數(shù)大都相近;而對超聲波清洗來說,模型1、5、6的擬合狀態(tài)參數(shù)大都相近,由表3可看出,模型5、6的脫附時間指數(shù)n均接近于1且模型6的脫附常數(shù)m也接近于1。綜上所述,模型6(即H-S分?jǐn)?shù)階脫附動力學(xué)模型)能夠更好地、統(tǒng)一地描述兩種清洗方式的驅(qū)酸洗滌過程。對磁力攪拌來說,擬合決定系數(shù)R與DC均在0.999以上、RMSE均小于0.012、RSS均小于0.0008;對超聲波清洗來說,擬合的R均在0.998以上、DC均在0.993以上、RMSE均小于0.031、RSS均小于0.006。

3.2 擴散模型的優(yōu)選

不同溫度下兩種清洗方式的擴散模型擬合狀態(tài)參數(shù)(R、DC、RMSE、RSS)分布云圖如圖3與圖4所示。其中模型9的機理函數(shù)積分形式將擬合的因變量進(jìn)一步擴大,基于RMSE與RSS的定義與計算公式,因變量的擴大直接導(dǎo)致RMSE與RSS的值擴大,使得RMSE與RSS的衡量準(zhǔn)則失效,故在圖3與圖4中,舍去模型9的RMSE與RSS值。此外,由于在擴散模型擬合中,因變量的值與機理函數(shù)的積分形式關(guān)系密切,因此擴散模型的優(yōu)選性僅能使用作為相對值的R與DC來衡量。

圖3 不同溫度下磁力攪拌驅(qū)酸洗滌擴散模型的擬合狀態(tài)參數(shù)云圖Fig.3 Cloud map of fitting parameters of diffusion models for the magnetic stirring cleaning at different temperatures

圖4 不同溫度下超聲驅(qū)酸洗滌擴散模型的擬合狀態(tài)參數(shù)云圖Fig.4 Cloud map of fitting parameters of diffusion models for the ultrasonic cleaning at different temperatures

如圖3所示,針對不同洗滌溫度下的磁力攪拌清洗來說,相關(guān)系數(shù)R的云圖中僅在模型6位置出現(xiàn)“紅色縱向貫穿”的現(xiàn)象,而在DC的云圖中出現(xiàn)“紅色縱向貫穿”的位置位于模型4、6、7,在RMSE、RSS的云圖中出現(xiàn)“藍(lán)色縱向貫穿”的位置位于模型6、7、8。因此用來描述磁力攪拌驅(qū)酸過程的最優(yōu)模型為模型6。

由圖4可知,在R云圖中,雖未出現(xiàn)“紅色縱向貫穿”,但在模型6的位置時縱向紅色區(qū)域較多,在DC的云圖中出現(xiàn)“紅色縱向貫穿”的位置位于模型4、6、7,在RMSE的云圖中藍(lán)色區(qū)域的現(xiàn)象與圖4相應(yīng)的云圖一致,在RSS云圖中,出現(xiàn)“藍(lán)色縱向貫穿”的位置位于模型4、6、7、8。因此,最優(yōu)的超聲波驅(qū)酸洗滌擴散模型為模型6。

綜上所述,擴散模型6,即Jander方程V類,能夠更好且統(tǒng)一地描述兩種清洗方式驅(qū)酸洗滌過程,其動力學(xué)相關(guān)參數(shù)為:當(dāng)洗滌方式為磁力攪拌時,活化能為13.49kJ/mol,指前因子為0.973min-1;當(dāng)洗滌方式為超聲波清洗時,活化能為11.02kJ/mol,指前因子為0.569min-1。

3.3 最優(yōu)驅(qū)酸洗滌模型的選擇

由上述對脫附模型與擴散模型的分析可得,兩種較優(yōu)的動力學(xué)方程如下式所示:

(9)

(10)

將上述兩式變化為:

(11)

qi=qmax{1-[Aexp(-E/RT)t]1/2}3

(12)

基于兩種清洗方式的實驗數(shù)據(jù),將以脫附模型6與擴散模型6得到的預(yù)測曲線進(jìn)行對比分析,其結(jié)果如圖5和圖6所示。

圖5 磁力攪拌清洗數(shù)據(jù)與脫附模型及擴散模型對比Fig.5 Comparison of data from experiments, the desorption model and the diffusion model for the magnetic stirring cleaning

圖6 超聲波清洗數(shù)據(jù)與脫附模型及擴散模型對比Fig.6 Comparison of data from experiments, the desorption model and the diffusion model for the ultrasonic cleaning

由圖5和圖6可知,脫附模型6的預(yù)測曲線均更接近于實驗數(shù)據(jù),且擬合的相關(guān)系數(shù)均在0.998以上,高于李婷婷等[13]的擬合結(jié)果。結(jié)果表明,用來描述的RDX驅(qū)酸洗滌過程的最優(yōu)動力學(xué)模型為脫附模型6,即H-S分?jǐn)?shù)階脫附動力學(xué)模型。此外,還可以看出,隨著溫度的升高,擴散模型越接近于脫附模型。這是因為洗滌溫度增加,由于脫附活化能高于擴散活化能,RDX表面脫附速率增長的幅度較大,從而導(dǎo)致其驅(qū)酸洗滌過程的速率控制步驟由脫附為主的綜合過程逐步轉(zhuǎn)化為酸的擴散過程。

綜上所述,RDX的驅(qū)酸洗滌動力學(xué)模型如式(13)所示:

(13)

式中:當(dāng)洗滌方式為磁力攪拌時,E=20.05kJ/mol、A=32.070min-1、n=1.4682、m=1.2984;當(dāng)洗滌方式為超聲波清洗時,E=21.42kJ/mol、A=329.516min-1、n=0.9856、m=1.0545。

此外,在式(13)中,脫附時間指數(shù)n體現(xiàn)了RDX驅(qū)酸洗滌過程中可能存在的復(fù)雜脫附機理綜合作用。

將超聲波清洗的動力學(xué)參數(shù)與磁力攪拌的情況相比(見表3),其表觀活化能的值僅增加了1.37kJ/mol,變化較小,而指前因子的值約提高了9倍。雖然超聲波清洗的脫附活化能稍高,但其指前因子的增加非常明顯,兩種清洗方式的酸脫除速率常數(shù)對比如圖7所示。

圖7 兩種清洗方式的酸脫除速率常數(shù)對比Fig.7 Comparison of rate constants of removal acid for two washing methods

由圖7可知,當(dāng)洗滌溫度為20～50℃時,與磁力攪拌清洗相比,超聲波清洗時的速率常數(shù)增加了4.86～5.17倍。結(jié)果表明,雖然超聲清洗的脫附能壘稍高,但即將脫附的酸活性位點數(shù)目較大,總體來說,反應(yīng)速率常數(shù)較大。此結(jié)果主要因為,超聲波的引入使其空化效應(yīng)在含RDX粗品的水溶液中進(jìn)行作用,從而引起溶液的湍動效應(yīng)與微擾效應(yīng),即空化效應(yīng)產(chǎn)生的聲沖擊波引起溶液宏觀湍流并于RDX晶體缺陷邊界發(fā)生氣泡破裂引起局部的微射流。而磁力攪拌僅是通過攪拌葉輪轉(zhuǎn)動引起溶液的湍流效應(yīng)。兩者相比,超聲波清洗較大幅度地提高了RDX表面吸附的酸脫效率。

4 結(jié) 論

(1)H-S分?jǐn)?shù)階脫附模型能夠較為準(zhǔn)確地描述RDX驅(qū)酸洗滌過程,其擬合相關(guān)系數(shù)均在0.998以上、擬合決定系數(shù)均在0.993以上。

(2)在本研究的實驗條件下,當(dāng)洗滌方式為磁力攪拌時,脫附活化能為20.05kJ/mol、指前因子為32.070min-1、脫附時間指數(shù)為1.4682、脫附模型常數(shù)為1.2984;當(dāng)洗滌方式為超聲波清洗時,脫附活化能為21.42kJ/mol、指前因子為329.516min-1、脫附時間指數(shù)為0.9856、脫附模型常數(shù)為1.0545。

(3)當(dāng)洗滌溫度為20～50℃時,與磁力攪拌清洗相比,超聲波清洗時RDX表面吸附的酸脫除速率約增大了4.86～5.17倍,且隨著洗滌溫度的升高,驅(qū)酸洗滌過程的速率控制步驟逐漸由脫附為主的綜合過程逐步轉(zhuǎn)化為酸的擴散過程。