田入園，張領(lǐng)科

(南京理工大學(xué) 能源與動(dòng)力工程學(xué)院，南京 210094)

【化學(xué)工程與材料科學(xué)】

環(huán)境壓力對(duì)鋁氧反應(yīng)點(diǎn)火過(guò)程影響的計(jì)算分析

田入園，張領(lǐng)科

(南京理工大學(xué) 能源與動(dòng)力工程學(xué)院，南京 210094)

為了研究環(huán)境的壓力對(duì)鋁氧反應(yīng)點(diǎn)火過(guò)程的影響，運(yùn)用零維均質(zhì)點(diǎn)火模型，利用CHEMKIN-PRO閉式均相反應(yīng)器進(jìn)行求解計(jì)算；分析了點(diǎn)火過(guò)程中主要基元反應(yīng)對(duì)點(diǎn)火過(guò)程的影響，計(jì)算了參與點(diǎn)火過(guò)程的主要物質(zhì)摩爾分?jǐn)?shù)隨時(shí)間的變化規(guī)律；重點(diǎn)計(jì)算了在2 300 K初始溫度不同環(huán)境壓力下(1～9 atm)的點(diǎn)火延遲時(shí)間，點(diǎn)火延遲時(shí)間與環(huán)境壓力的擬合關(guān)系式為t=3.36×10-7p0-0.99；獲得了不同環(huán)境壓力時(shí)組分O2、Al(l)、Al2O、Al2O2、AlO、Al2O3、O、AlO2、Al、Al2O(l)在點(diǎn)火過(guò)程中的摩爾分?jǐn)?shù)變化規(guī)律。

含能材料；鋁氧反應(yīng)；點(diǎn)火過(guò)程；環(huán)境壓力；數(shù)值分析

鋁作為金屬添加劑被廣泛地運(yùn)用于各種煙火、炸藥以及推進(jìn)劑中，它不僅可以抑制推進(jìn)劑的不穩(wěn)定燃燒，也能同氧化劑燃燒反應(yīng)放熱提高推進(jìn)劑的爆熱與比沖，其氧化產(chǎn)物Al2O3的氣化熱較高，冷凝時(shí)的放熱將進(jìn)一步增加推進(jìn)劑能量[1]。Alavi等[2]和Puri等[3]運(yùn)用分子動(dòng)力學(xué)分別研究了粒徑在1.0～3.4 nm以及2～9 nm鋁顆粒的熔點(diǎn)，結(jié)果表明，粒徑在1 nm的鋁顆粒的熔點(diǎn)為400 K，明顯低于大顆粒鋁的熔點(diǎn)(933 K)。Bazyn等[4]研究了初始溫度為2 640 K，環(huán)境壓力為3～30 atm粒徑為10 μm的鋁顆粒的燃燒過(guò)程，發(fā)現(xiàn)該粒徑鋁顆粒的燃燒時(shí)間受環(huán)境壓力影響非常顯著。Park等[5]和Trunov等[6]通過(guò)對(duì)鋁顆粒的點(diǎn)火燃燒實(shí)驗(yàn)研究，發(fā)現(xiàn)大顆粒鋁粒子在高壓下的燃燒主要受擴(kuò)散控制，小顆粒鋁粒子在低壓環(huán)境下主要受機(jī)理控制影響(即小顆粒燃燒時(shí)，顆粒表面反應(yīng)速率比氧化劑擴(kuò)散速率快得多)。Washburn等[7]提出納米鋁顆粒在點(diǎn)火過(guò)程中快速氣化后同氧化劑反應(yīng)，主要受反應(yīng)機(jī)理控制。Lynch等[8]通過(guò)納米鋁顆粒的燃燒實(shí)驗(yàn)證明了其反應(yīng)主要在氣相中進(jìn)行。

本文通過(guò)建立鋁氧零維均質(zhì)氣相點(diǎn)火模型，忽略顆粒尺寸大小，基于鋁與氧的反應(yīng)機(jī)理，主要研究鋁點(diǎn)火過(guò)程中參與反應(yīng)的物質(zhì)組分變化、2 300 K的初始環(huán)境溫度下不同環(huán)境壓力對(duì)反應(yīng)進(jìn)程的影響以及點(diǎn)火延遲時(shí)間與環(huán)境壓力的關(guān)系。

1 計(jì)算模型

為了簡(jiǎn)化計(jì)算，提出如下基本假設(shè)：鋁顆粒與空氣均勻混合在恒壓p0密閉環(huán)境中，在點(diǎn)火持續(xù)期內(nèi)試驗(yàn)環(huán)境為絕熱；為了排除輸運(yùn)(對(duì)流與擴(kuò)散)等因素的影響，在任意時(shí)刻該反應(yīng)環(huán)境溫度以及物質(zhì)組分的分布是均勻的；假設(shè)空氣中O2與N2的體積比VO2∶VN2=1∶3.76，忽略空氣中N2與其他氣體對(duì)反應(yīng)的影響，鋁與氧氣的總包反應(yīng)可以描述為2Al+1.5O2→Al2O3；只考慮氣相動(dòng)力學(xué)反應(yīng)，忽略表面動(dòng)力學(xué)反應(yīng)。

對(duì)于該均質(zhì)系統(tǒng)，各組分質(zhì)量分?jǐn)?shù)Yk(k=1，2，…，n′，n′為組分的數(shù)量)和溫度T的控制方程[9]：

(1)

(2)

式(1)中：ωk為k中組分的凈生成率；Wk為第k種組分的分子質(zhì)量；t為時(shí)間(s)；ρ為密度(kg/m3)；cp為定壓比熱容(J/(kg·K))；hk為k組分的焓值(J/kg)。

鋁點(diǎn)火反應(yīng)機(jī)理中動(dòng)力學(xué)數(shù)據(jù)由式(3)與式(4)確定，熱力學(xué)數(shù)據(jù)由式(5)確定：

(3)

kfi=AiTniexp(-Ei/RT)

(4)

(5)

其中i為化學(xué)反應(yīng)式的數(shù)量；k為組分的數(shù)量；vk為第k種組分的擴(kuò)散速度；上標(biāo)“ ′ ”為正方向的化學(xué)計(jì)量系數(shù)，上標(biāo)“ ″ ”為反方向的化學(xué)計(jì)量系數(shù)；χk為第k種組分的摩爾分?jǐn)?shù)；kfi為反應(yīng)速率；Ai為指前因子；ni為溫度系數(shù)；Ei為活化能；R為普適氣體常數(shù)。

本文采用Catoire等[10-11]提出的Al-O氣相反應(yīng)動(dòng)力學(xué)機(jī)理模擬分析納米鋁顆粒在空氣中點(diǎn)火燃燒的化學(xué)反應(yīng)過(guò)程，該機(jī)理由8種組分10個(gè)反應(yīng)式(R1～R10)組成，如表1所示。根據(jù)Starik等[12]的研究，在機(jī)理中添加了O2分解反應(yīng)R11。假設(shè)初始狀態(tài)下的納米鋁顆粒為液態(tài)，在反應(yīng)過(guò)程中需要考慮鋁的相變，所以在研究中加入了一個(gè)具有較高速率常數(shù)的局部平衡相變反應(yīng)過(guò)程式R12。由于Al2O3會(huì)迅速分解成低能氧化物和O等物質(zhì)，因此加入反應(yīng)式R13[13]。

表1 鋁與氧的反應(yīng)機(jī)理

2 點(diǎn)火模擬結(jié)果與分析

針對(duì)表1鋁與氧的反應(yīng)機(jī)理，采用CHEMKIN-PRO軟件包中封閉0維反應(yīng)器(閉式均相反應(yīng)器)進(jìn)行求解計(jì)算，反應(yīng)所涉及物質(zhì)的熱力學(xué)系數(shù)均通過(guò)JANNAF熱力學(xué)數(shù)據(jù)表[14]和Burcat熱力學(xué)數(shù)據(jù)庫(kù)[15]獲取。計(jì)算得到了環(huán)境壓力為1 atm，溫度為2 300 K時(shí)參與反應(yīng)的主要物質(zhì)摩爾分?jǐn)?shù)與環(huán)境溫度T隨時(shí)間變化曲線(xiàn)，如圖1(a)所示。

由圖1(a)可以看出，在相同工況下與文獻(xiàn)[13]所計(jì)算出的結(jié)果圖1(b)吻合較好；當(dāng)時(shí)間到達(dá)1×10-7s時(shí)各物質(zhì)摩爾分?jǐn)?shù)與溫度基本達(dá)到穩(wěn)定，此時(shí)認(rèn)為完成點(diǎn)火過(guò)程，溫度約為3 819 K，該條件下的點(diǎn)火反應(yīng)過(guò)程主要集中在1.0×10-8～1.0×10-7s時(shí)間。

分別計(jì)算了初始溫度為2 300 K，環(huán)境壓力p0分別為1 atm、3 atm、5 atm、7 atm和9 atm的鋁氧反應(yīng)點(diǎn)火過(guò)程，不同環(huán)境壓力下點(diǎn)火過(guò)程溫度變化曲線(xiàn)如圖2所示；環(huán)境壓力與點(diǎn)火延遲時(shí)間的關(guān)系如圖3所示；得到了點(diǎn)火過(guò)程中反應(yīng)物、中間產(chǎn)物和生成物摩爾分?jǐn)?shù)隨時(shí)間變化曲線(xiàn)，分別如圖4～圖6所示。

圖1 點(diǎn)火過(guò)程主要物質(zhì)摩爾分?jǐn)?shù)與溫度隨時(shí)間變化曲線(xiàn)

圖2 不同環(huán)境壓力下點(diǎn)火過(guò)程溫度變化曲線(xiàn)

圖3 環(huán)境壓力與點(diǎn)火延遲時(shí)間曲線(xiàn)

由圖2中可以看出：環(huán)境壓力越高，鋁點(diǎn)火延遲時(shí)間越短。這是因?yàn)閴毫ι呒涌炝虽X與氧氣的氧化反應(yīng)進(jìn)程；環(huán)境壓力p0的不同對(duì)點(diǎn)火完成時(shí)的溫度影響比較大，如當(dāng)環(huán)境壓力為1 atm時(shí)，點(diǎn)火溫度約為3 819 K，環(huán)境壓力為9 atm時(shí)，點(diǎn)火溫度約為4 141 K。由圖3可以看出點(diǎn)火延遲時(shí)間隨環(huán)境壓力的增加而減小，其變化趨勢(shì)與文獻(xiàn)[16]的實(shí)驗(yàn)結(jié)果吻合。通過(guò)擬合點(diǎn)火延遲時(shí)間與環(huán)境壓力的關(guān)系可得兩者關(guān)系近似成反比關(guān)系，擬合曲線(xiàn)關(guān)系式為t=3.36×10-7p0-0.99。

由圖4中可以看出O2和Al(l)變化趨勢(shì)基本相同，Al(l)的摩爾分?jǐn)?shù)首先減少。這是因?yàn)榉磻?yīng)是在氣相中進(jìn)行，Al(l)首先需要通過(guò)假定的反應(yīng)式R12蒸發(fā)成Al，隨后進(jìn)行反應(yīng)R2使得Al(l)和O2的摩爾分?jǐn)?shù)下降趨勢(shì)相同；不同的環(huán)境壓力改變了O2與Al(l)摩爾分?jǐn)?shù)達(dá)到平衡的時(shí)間，也改變了O2達(dá)到平衡時(shí)的摩爾分?jǐn)?shù)；環(huán)境壓力越大達(dá)到平衡時(shí)的O2摩爾分?jǐn)?shù)就越大。

圖4 不同環(huán)境壓力下反應(yīng)物摩爾分?jǐn)?shù)變化曲線(xiàn)

由圖5(a)可以看出在Al2O的摩爾分?jǐn)?shù)首先迅速增至峰值，隨后又快速下降達(dá)到平衡位置。這是因?yàn)殡S著R2反應(yīng)的進(jìn)行環(huán)境中不斷生成AlO，AlO與Al結(jié)合通過(guò)反應(yīng)R5生成Al2O，隨后Al2O的摩爾分?jǐn)?shù)下降的原因是因?yàn)锳l2O與O結(jié)合生成Al2O2(反應(yīng)式R8)，從而使Al2O2的摩爾分?jǐn)?shù)增加(圖5(b))，生成的Al2O2緊接著又和O結(jié)合生成Al2O3(反應(yīng)式R9)，所以Al2O2的摩爾分?jǐn)?shù)又下降，最終達(dá)到穩(wěn)定值；從圖5可以看出環(huán)境壓力的不同并沒(méi)有改變反應(yīng)中間產(chǎn)物的摩爾分?jǐn)?shù)變化趨勢(shì)，但改變了中間產(chǎn)物到達(dá)平衡時(shí)的摩爾分?jǐn)?shù)，改變了反應(yīng)歷程的時(shí)間；環(huán)境壓力越高反應(yīng)中間產(chǎn)物Al2O、Al2O2、AlO、Al2O3、AlO2達(dá)到平衡時(shí)的摩爾分?jǐn)?shù)也越大，中間產(chǎn)物O的摩爾分?jǐn)?shù)降低，使得燃燒更充分，點(diǎn)火完成時(shí)的溫度越高(圖2)。

圖5 不同環(huán)境壓力下中間產(chǎn)物摩爾分?jǐn)?shù)隨時(shí)間變化曲線(xiàn)

圖6 不同環(huán)境壓力下生成物摩爾分?jǐn)?shù)變化曲線(xiàn)

由圖6中可以看出：Al的摩爾分?jǐn)?shù)先增加再減少，這是因?yàn)辄c(diǎn)火初期Al(l)蒸發(fā)成Al，使Al的摩爾分?jǐn)?shù)上升至最大，隨后Al與O2和AlO反應(yīng)(反應(yīng)式R2和R5)消耗Al，使得Al的摩爾分?jǐn)?shù)下降，當(dāng)Al(l)的蒸發(fā)量與Al的消耗量達(dá)到平衡時(shí)，此時(shí)Al的摩爾分?jǐn)?shù)不變，點(diǎn)火完成；圖6(a)顯示反應(yīng)環(huán)境壓力越高，點(diǎn)火過(guò)程中Al的摩爾分?jǐn)?shù)上升過(guò)程中所到達(dá)的峰值越小，反應(yīng)平衡時(shí)Al的摩爾分?jǐn)?shù)也小，說(shuō)明環(huán)境壓力高不利于鋁的蒸發(fā)，但有利于鋁蒸氣的消耗；圖6(b)顯示環(huán)境壓力的變化并不影響Al2O3(l)平衡時(shí)的摩爾分?jǐn)?shù)，說(shuō)明環(huán)境壓力的變化并不會(huì)改變反應(yīng)平衡時(shí)Al2O3(l)的摩爾分?jǐn)?shù)。

3 結(jié)論

通過(guò)對(duì)鋁與氧氣反應(yīng)的機(jī)理研究，認(rèn)為鋁在點(diǎn)火過(guò)程中所經(jīng)歷的主要反應(yīng)過(guò)程為R2、R5、R8和R9；計(jì)算結(jié)果表明：在初始溫度為2 300 K，環(huán)境壓力為1 atm時(shí)的點(diǎn)火過(guò)程中，主要物質(zhì)摩爾分?jǐn)?shù)隨時(shí)間與溫度變化特征與文獻(xiàn)[13]吻合較好。通過(guò)環(huán)境壓力的變化對(duì)點(diǎn)火延遲時(shí)間的計(jì)算分析，發(fā)現(xiàn)點(diǎn)火環(huán)境壓力的升高明顯縮短了點(diǎn)火延遲時(shí)間，且點(diǎn)火延遲時(shí)間與環(huán)境壓力近似成反比關(guān)系，擬合曲線(xiàn)關(guān)系式為t=3.36×10-7p0-0.99。環(huán)境壓力越高，參與反應(yīng)的中間產(chǎn)物Al2O、Al2O2、AlO、Al2O3、AlO2達(dá)到平衡時(shí)的摩爾分?jǐn)?shù)也就越大，而中間產(chǎn)物O的摩爾分?jǐn)?shù)越低；Al的次級(jí)氧化物的摩爾分?jǐn)?shù)增加有利于點(diǎn)火過(guò)程，點(diǎn)火完成時(shí)的溫度越高。

[1] 李穎,宋武林,謝長(zhǎng)生,等.納米鋁粉在固體推進(jìn)劑中的應(yīng)用進(jìn)展[J].兵工學(xué)報(bào),2005,26(1):121-125.

[2] ALAVI S,THOMPSON D L.Molecular Dynamics Simulations of The Melting of Aluminum Nanoparticles[J]. The Journal of Physical Chemistry A，2006,110(4):1518-1523.

[3] PURI P,YANG V.Effect of Particle Size on Melting of Aluminum at Nano Scales [J].The Journal of Physical Chemistry C，2007,111(32):11776-11783.

[4] BAZYN T,KRIER H,GLUMAC N.Oxidizer And Pressure Effects on The Combustion of 10-Micron Aluminum Particles [J].Journal of Propulsion and Power，2005,21(4):577-582.

[5] PARK K,LEE D,RAI A,et al.Size-Resolved Kinetic Measurements of Aluminum Nanoparticle Oxidation With Single Particle Mass Spectrometry [J].The Journal of Physical Chemistry B，2005,109(15):7290-7299.

[6] TRUNOV M A,SCHOENITZ M,DREIZIN E L.Effect of Polymorphic Phase Transformations in Alumina Layer on Ignition of Aluminium Particles [J].Combustion Theory and Modelling，2006,10(4):603-623.

[7] WASHBURN E B,TRIVEDI J N,CATOIRE L,et al.The Simulation of The Combustion of Micrometer-Sized Aluminum Particles With Steam [J]. Combustion Science And Technology，2008,180(8):1502-1517.

[8] LYNCH P,FIORE G,KRIER H,et al.Gas-phase Reaction in Nanoaluminum Combustion[J].Combustion Science and Technology，2010,182(7):842-857.

[9] KEE R J,RUPLEY F M,MEEKS E,et al.CHEMKIN-III:A FORTRAN Chemical Kinetics Package for The Analysis of Gas-Phase Chemical And Plasma Kinetics [R].Sandia National Laboratories Report.1996,SAND96-8216.

[10]CATOIRE L,JEAN-FRAN-ATILDE,LEGENDRE O,et al.Kinetic Model For Aluminum-sensitized Ram Accelerator Combustion [J].Journal Of Propulsion And Power，2003,19(2):196-202.

[11]SWIHART M T,CATOIRE L,LEGRAND B,et al.Rate Constants For The Homogeneous Gas-Phase Al/HCl Combustion Chemistry[J].Combustion And Flame，2003,132(1):91-101.

[12]STARIK A M,KULESHOV P S,SHARIPOV A S,et al.Numerical Analysis of Nanoaluminum Combustion in Steam[J]. Combustion And Flame，2014,161(6):1659-1667.

[13]HUANG Y,RISHA G A,YANG V,et al.Effect of Particle Size on Combustion of Aluminum Particle Dust in Air [J]. Combustion And Flame，2009,156(1):5-13.

[14]CHASE M W.JANAF Thermochemical Tables [J].JANAF Thermochemical Tables,By Chase,MW Washington,DC:American Chemical Society; New York:American Institute of Physics for The National Bureau of Standards,C1986.United States.National Bureau of Standards.1986.

[15]BURCAT A,RUSCIC B.Third Millenium Ideal Gas And Condensed Phase Thermochemical Database For Combustion With Updates From Active Thermochemical Tables[M].Argonne National Laboratory Argonne,IL,2005.

[16]彭晶,肖進(jìn).壓力對(duì)納米鋁粉燃燒性能的影響[J].化學(xué)推進(jìn)劑與高分子材料,2015,13(3):53-58.

(責(zé)任編輯楊繼森)

Numerical Analysis of Aluminum/Oxygen Ignition Under Different Environmental Pressure

TIAN Ru-yuan, ZHANG Ling-ke

(School of Energy and Power Engineering, Nanjing University of Science and Technology, Nanjing 210094, China)

In order to study the influence of the environmental pressure to aluminum ignition in gas phase, the zero-dimensional model of homogeneous ignition was established with Al-O2reaction mechanism and calculated by using the CHEMKIN-PRO software package. The influences of main elementary reactions of the ignition process were illustrated. The mole fraction of main substances involved in the ignition process versus time was obtained and analyzed. The ignition delay time on different environmental pressure(1～9 atm)were analyzed, and their relationship is thatt=3.36×10-7p0-0.99. And the change rule of mole fraction influence of the environmental pressure to O2, Al(l), Al2O, Al2O2, AlO, Al2O3, O, AlO2, Al, Al2O(l) were obtained in the ignition process respectively.

energetic materials; Al/O reaction; ignition process; environmental pressure; numerical analysis

2016-07-22；

國(guó)家自然科學(xué)基金資助(51106075)

田入園(1992—)，男，碩士研究生，主要從事金屬添加劑燃燒特性研究。

10.11809/scbgxb2016.12.034

田入園,張領(lǐng)科.環(huán)境壓力對(duì)鋁氧反應(yīng)點(diǎn)火過(guò)程影響的計(jì)算分析[J].兵器裝備工程學(xué)報(bào)，2016(12):150-154.

format:TIAN Ru-yuan, ZHANG Ling-ke.Numerical Analysis of Aluminum/Oxygen Ignition Under Different Environmental Pressure[J].Journal of Ordnance Equipment Engineering,2016(12):150-154.

V311+.3

A

2096-2304(2016)12-0150-05

修回日期：2016-08-23