曲踐，李保衛(wèi)，鄭坤燦，武文斐

?

O2/CO2氣氛下多種碳基隨機孔模型的建立

曲踐1,2，李保衛(wèi)1，鄭坤燦1,2，武文斐1,2

（1內(nèi)蒙古自治區(qū)白云鄂博礦多金屬資源綜合利用國家重點實驗室，內(nèi)蒙古包頭 014010；2內(nèi)蒙古科技大學能源與環(huán)境學院，內(nèi)蒙古包頭 014010）

傳統(tǒng)隨機孔模型基于簡單一步反應(yīng)建立，不適用于處理O2/CO2氣氛下焦炭顆粒復(fù)雜氣固反應(yīng)。針對此問題，基于焦炭本身具有多種碳基的特點，以及焦炭顆粒在O2/CO2氣氛下燃燒的特性，建立復(fù)雜氣固反應(yīng)下的多種碳基隨機孔模型和孔隙結(jié)構(gòu)模型。模擬直徑為100 μm的焦炭顆粒在O2/CO2氣氛下燃燒的過程，使用FORTRAN語言自主編程計算并分析結(jié)果。研究表明，燃燒初期顆粒呈現(xiàn)競爭效應(yīng)，孔隙內(nèi)部氣體濃度產(chǎn)生劇烈波動。波動的生成原因是化學反應(yīng)與物理擴散之間的競爭，可以通過增加環(huán)境氧濃度和減小焦炭顆粒粒徑來改善。所提出的多種碳基隨機孔模型對于表征O2/CO2氣氛下焦炭顆粒的燃燒特性有著良好的適應(yīng)性。

O2/CO2氣氛；多孔介質(zhì)；隨機孔模型；單顆粒焦炭燃燒；傳質(zhì)；介尺度；數(shù)值模擬

引 言

隨著全球氣候變暖等環(huán)境問題的加劇，富氧燃燒技術(shù)在節(jié)能減排方面體現(xiàn)了巨大的優(yōu)勢，因此越來越受到全球范圍研究者們的關(guān)注[1-5]。與此同時，富氧燃燒技術(shù)在工程應(yīng)用中也暴露了煤焦顆粒燃燒溫度低、燃燒速率慢等缺點，因此，對富氧燃燒技術(shù)進行深入研究，對于進一步節(jié)約能源、降低CO2等溫室氣體的排放有著重要意義。

目前研究煤焦顆粒燃燒特性的手段包括實驗和單顆粒數(shù)值模擬，其中針對O2/CO2氣氛的研究大多以實驗為主。Struis等[6]在動力學控制條件下采用等溫熱重分析儀研究煤焦在O2/CO2氣氛下氣化反應(yīng)。Bejarano等[7]采用非嵌入式光學彩色測溫儀來測量單顆粒燃燒條件下?lián)]發(fā)分火焰及煤焦燃燒情況。實驗盡管能夠從宏觀角度描述燃燒現(xiàn)象，但是無法分析燃燒的機理。

單顆粒模型的數(shù)值計算因為能夠分析顆粒燃燒的介觀機理而受到研究者們的關(guān)注，諸多研究者提出了表征煤焦顆粒燃燒的數(shù)學模型[8-10]。其中Bhatia等[11-12]所提出的隨機孔模型由于能夠解釋最大反應(yīng)速率現(xiàn)象而被國內(nèi)外研究者們廣泛采用與修 正[6,13-14]。然而基于一步反應(yīng)的傳統(tǒng)隨機孔模型無法適應(yīng)O2/CO2氣氛下的復(fù)雜氣固反應(yīng)，Lu等[15-16]利用隨機孔模型表征煤焦與H2O和CO2的反應(yīng)，發(fā)現(xiàn)模擬結(jié)果與實驗數(shù)據(jù)并不符合。研究者們目前仍然難以將隨機孔模型應(yīng)用于O2/CO2氣氛下模擬焦炭顆粒復(fù)雜燃燒過程的研究。

本工作在前人研究的基礎(chǔ)上，建立了一個基于O2/CO2氣氛下多步復(fù)雜氣固化學反應(yīng)的多種碳基隨機孔模型（various char-RPM）。對Maffei等[17]研究者提出的詳細化學反應(yīng)方程加以整理改進，提出八步異相反應(yīng)加一步同相反應(yīng)表征焦炭顆粒燃燒，整理出3種氣體及3種碳基的反應(yīng)速率分別進行計算。新提出的模型能夠適應(yīng)O2/CO2氣氛下多組分多步復(fù)雜氣固反應(yīng)，并通過實驗數(shù)據(jù)驗證其適 應(yīng)性。

1 物理模型

本工作提出一個直徑為100 μm的焦炭顆粒在氧氣濃度為30%的O2/CO2氣氛下燃燒的物理模型，如圖1所示。焦炭是煤粉經(jīng)熱解析出揮發(fā)分形成的，經(jīng)歷了軟化、塑形及再凝固階段，產(chǎn)生了強度較 高的骨架結(jié)構(gòu)，并在內(nèi)部形成了發(fā)達的孔道。根據(jù)焦炭本身的特性及文獻[6,11-15,17-20]中的研究，對焦炭顆粒做出以下合理假設(shè)：①焦炭顆粒可近似為具有各向同性的球體；②顆粒內(nèi)部包括大孔和小孔，大孔主要為反應(yīng)氣體的運輸通道，小孔主要提供化學反應(yīng)面積；③大孔小孔均為迂曲度相同的曲折孔，均勻分布在焦炭顆粒中，氣體在孔隙中擴散包括容積擴散和努森擴散，不考慮外部擴散；④焦炭顆粒經(jīng)高溫熱解，具有強度較高的骨架結(jié)構(gòu)，因此視為反應(yīng)過程中不會發(fā)生坍塌與破碎，且燃燒反應(yīng)發(fā)生在內(nèi)孔表面，因此視為整體粒徑不發(fā)生變化；⑤顆粒粒徑較小，產(chǎn)生的化學熱很快釋放到環(huán)境中，使得顆粒溫度與環(huán)境溫度相差不大，因此不考慮輻射換熱。

圖1 焦炭顆粒物理模型

2 數(shù)學模型

2.1 化學反應(yīng)方程及反應(yīng)速率

單顆粒焦炭經(jīng)高溫熱處理，其中的碳基構(gòu)成會發(fā)生變化，部分無序結(jié)構(gòu)部分轉(zhuǎn)化為有序結(jié)構(gòu)[19-20]。因此在焦炭顆粒燃燒過程中，碳基有3種存在形式，分別為碳氫比為2:1的碳氫化合物（CHARH）,焦炭中的無序結(jié)構(gòu)（CHARC）和有序結(jié)構(gòu)（CHARG），3種形式碳基的初始相對數(shù)量取決于煤種。焦炭燃燒的化學反應(yīng)即這3種存在形式分別與反應(yīng)氣體發(fā)生氧化反應(yīng)，焦炭顆粒燃燒過程中的化學反應(yīng)方程及動力學參數(shù)見表1[21-26]。

表1 焦炭燃燒過程化學反應(yīng)方程及動力學參數(shù)

以上9步反應(yīng)得到了廣泛的驗證[27]，通過熱重分析儀、夾帶流反應(yīng)器等不同的實驗裝置獲取大量的實驗數(shù)據(jù)，并與其他研究者的研究相比較，保證整個動力學方案具有高度的準確性和廣泛的適用范圍。以上反應(yīng)的反應(yīng)溫度近似，各個反應(yīng)發(fā)生時間有很高的重合度，因此本工作假設(shè)以上反應(yīng)在焦炭顆粒燃燒過程中同時進行，將各氣體的消耗速率與生成速率疊加，計算得各反應(yīng)氣體的化學反應(yīng)速率如下

在焦炭顆粒燃燒過程中，O2作為主要的反應(yīng)氣體，始終處于被消耗的狀態(tài)，因此的值為負。CO2涉及到的反應(yīng)較復(fù)雜，消耗與生成同時發(fā)生。CO為主要的生成氣體，只在第9步中作為反應(yīng)氣體被消耗。表1中的9步反應(yīng)不僅描述了3種氣體的生成情況，同時也描述了3種碳基的消耗情況，將各碳基的消耗速率與生成速率疊加，計算得各存在形式碳基的化學反應(yīng)速率如下

以上反應(yīng)中，碳基CHARH在與O2反應(yīng)生成CO的同時，也有部分轉(zhuǎn)化為CHARC。

2.2 熱質(zhì)傳遞方程

2.2.1 一維球坐標傳熱方程 焦炭燃燒為放熱反應(yīng)，燃燒過程中產(chǎn)生化學熱使顆粒整體溫度升高，同時將顆粒的溫度在計算過程中與化學反應(yīng)動力學參數(shù)相耦合。焦炭顆粒為典型多孔介質(zhì)，由于顆粒孔隙比表面積大，氣固傳熱快，在同一時間下顆粒內(nèi)部氣固溫度相同，將9步化學反應(yīng)過程中的放熱吸熱量疊加作為源相，忽略輻射放熱。建立一維球坐標傳熱方程

2.2.2 多組分氣體傳質(zhì)方程 焦炭顆粒內(nèi)部有著復(fù)雜且發(fā)達的孔隙結(jié)構(gòu)，這使得顆粒在燃燒過程中，環(huán)境中的氣體通過擴散作用進入孔隙內(nèi)部，同時內(nèi)部化學反應(yīng)產(chǎn)生的氣體也通過擴散作用離開顆粒。在焦炭顆粒燃燒過程中，顆?？椎纼?nèi)部同時存在O2、CO2、CO 3種氣體。其中O2來源于環(huán)境氣氛，是主要的反應(yīng)氣體，CO2主要來源于環(huán)境氣氛，少部分來源于化學反應(yīng)生成，CO來源于化學反應(yīng)。以表示氣體種類。建立一維球坐標多組分傳質(zhì)方程

式中的ej為顆粒孔隙內(nèi)部氣體的有效擴散系數(shù)。由于焦炭顆粒燃燒過程中，孔隙內(nèi)部同時存在3種氣體，因此需要基于多組分特性對傳統(tǒng)的二元計算方式進行修正。在計算過程中，將孔隙內(nèi)的氣體簡化為計算氣體和另外兩種氣體的混合氣體′。顆粒大孔內(nèi)部的擴散為容積擴散，主要作用為氣體運輸，小孔擴散為努森擴散，分子自由程大于孔隙直徑，分子與壁面碰撞頻率大，是主要的化學反應(yīng)場所，反應(yīng)氣體的消耗和產(chǎn)物的生成使孔隙內(nèi)部形成了氣體濃度差，進一步促進了擴散作用。本工作分別計算大孔內(nèi)的容積擴散系數(shù)bj，小孔內(nèi)的努森擴散系數(shù)kj，然后根據(jù)焦炭顆?？紫堵始翱紫队厍扔嬎愕玫娇紫秲?nèi)部氣體有效擴散系數(shù)ej。建立擴散系數(shù)方程如下

2.2.3 初始條件和邊界條件 反應(yīng)初始，顆粒各部分溫度一致，為1073 K。此時環(huán)境中的O2還沒有進入顆?？紫秲?nèi)部，CO尚未生成，因此O2初始濃度為0。

2.3 多種碳基隨機孔模型

焦炭顆粒由3種不同形式的碳基構(gòu)成，3種碳基對不同溫度和不同氣體的敏感度不同，在實際反應(yīng)過程中的轉(zhuǎn)化特性也不相同，為保證結(jié)果更加真實準確，應(yīng)當對3種碳基分別進行計算。在此基礎(chǔ)上進一步求出焦炭顆粒總轉(zhuǎn)化率。焦炭顆粒在燃燒過程中，孔隙率、迂曲度等結(jié)構(gòu)參數(shù)始終發(fā)生變化，有研究表明該變化與焦炭顆粒轉(zhuǎn)化率和初始結(jié)構(gòu)參數(shù)線性相關(guān)[17]。本工作將計算得出的焦炭顆粒轉(zhuǎn)化率耦合進結(jié)構(gòu)參數(shù)的變化方程中，建立多種碳基隨機孔模型及結(jié)構(gòu)參數(shù)變化方程如下

3 數(shù)值計算及模型的實驗驗證

3.1 數(shù)值計算方法及初始參數(shù)

本工作所用程序，采用有限體積法（finite volume method）對焦炭顆粒進行離散；為了保證輸出值穩(wěn)定且始終具有物理意義，計算格式選擇全隱式格式；采用附加源項法處理邊界條件；采用TDMA算法計算結(jié)果。整個程序基于Developer Studio環(huán)境，采用FORTRAN語言自主編程。顆粒各項初始參數(shù)見表2。

表2 計算參數(shù)

3.2 模型的實驗驗證

本工作采用Maffei等[17]于2013年所做的實驗進行模型適用性的驗證。實驗采用壁溫恒定1400 K的滴管爐，在O2/CO2氣氛下進行了粒徑為75～90 μm的煤焦顆粒燃燒實驗，模型及實驗所采用的煤焦工業(yè)分析及元素分析參數(shù)見表3。

表3 煤焦的工業(yè)分析與元素分析

在相同的條件下，采用本工作提出的多種碳基隨機孔模型（various char-RPM），Everson等[28]于2011年提出的考慮擴散傳輸?shù)男拚S機孔模型（modified RPM），Bhatia等[11-12]提出的傳統(tǒng)隨機孔模型（RPM），分別計算焦炭顆粒在不同環(huán)境O2濃度下轉(zhuǎn)化率達到90%所需要的時間，并與兩種煤樣在氧氣濃度為20%、30%、40%、60%、80%和100%的情況下燃盡時間的實驗數(shù)據(jù)進行對比，結(jié)果見圖2。

圖2 焦炭顆粒的燃盡時間實驗結(jié)果與數(shù)值模擬結(jié)果比較

通過對比結(jié)果可以看出，傳統(tǒng)隨機孔模型只考慮了簡單的一步反應(yīng)，通過O2的反應(yīng)速率表征焦炭燃燒，因此通過隨機孔模型計算得出結(jié)果過于依賴環(huán)境中的O2濃度，在O2濃度較高的情況下與實驗數(shù)據(jù)比較吻合，但是當O2濃度較低時，模型計算所得的燃盡時間明顯高于實驗數(shù)據(jù)。修正隨機孔模型考慮了擴散對焦炭顆粒燃燒的影響，并且對結(jié)構(gòu)參數(shù)進行了優(yōu)化，整體吻合度要優(yōu)于傳統(tǒng)隨機孔模型，但是并沒有針對O2/CO2氣氛進行調(diào)整，因此在低O2濃度時吻合度仍然較差。

多種碳基隨機孔模型在Everson等[28]提出的修正隨機孔模型的基礎(chǔ)上，采用9步反應(yīng)表征焦炭顆粒的燃燒情況，將氣體反應(yīng)速率與固體反應(yīng)速率分離，在考慮O2濃度的同時也考慮了環(huán)境中CO2濃度對煤焦顆粒燃燒所產(chǎn)生的影響，因此在高O2濃度和低O2濃度的條件下都能和實驗數(shù)據(jù)較好吻合，實驗驗證證明了氣固反應(yīng)隨機孔模型較傳統(tǒng)隨機孔模型有著更好的適應(yīng)性。

通過以上對美國幼兒身體活動環(huán)境相關(guān)政策的綜合分析，結(jié)合我國幼兒體育政策的現(xiàn)狀，為今后我國幼兒體育的發(fā)展提出以下建議：

4 結(jié)果分析

4.1 焦炭顆粒內(nèi)部氣體傳輸特性

圖3(a)、(b)、(c)分別為O2、CO2和CO 3種氣體濃度隨燃燒時間的變化曲線。從圖中可以看出，焦炭顆粒燃燒初始，內(nèi)部氣體的濃度極不穩(wěn)定，存在一個短暫的振幅為10%～40%的波動。這是因為焦炭顆粒處在動力擴散控制區(qū)，化學反應(yīng)的生成與消耗，物理傳輸?shù)倪M入與導出都影響了孔隙內(nèi)部氣體濃度的變化，呈現(xiàn)出競爭的效應(yīng)。CO2和CO的波動持續(xù)時間更長，持續(xù)了20 ms左右，是因為二者在化學反應(yīng)中既是反應(yīng)氣體又是生成氣體，所產(chǎn)生競爭效應(yīng)也更為復(fù)雜。從顆粒核心、中間區(qū)域及外表面的濃度曲線可以看出，顆粒內(nèi)部沿半徑方向存在濃度差，這是由于顆粒多孔結(jié)構(gòu)導致的。

圖3 焦炭顆粒中氣體濃度隨燃燒時間的變化

4.2 焦炭顆粒的轉(zhuǎn)化特性

圖4(a)為焦炭顆粒整體的轉(zhuǎn)化率隨燃燒時間的變化曲線，圖4(b)為焦炭顆粒不同徑向長度的轉(zhuǎn)化速率隨焦炭整體轉(zhuǎn)化率的變化曲線，圖4(c)為碳基CHARH，CHARC和CHARG3種不同形式碳基轉(zhuǎn)化率隨燃燒時間的變化曲線，圖4(d)為碳基CHARH，CHARC和CHARG3種不同形式碳基轉(zhuǎn)化速率隨焦炭整體轉(zhuǎn)化率的變化曲線。

圖4 焦炭的轉(zhuǎn)化率及轉(zhuǎn)化速率的變化

結(jié)合圖4(a)、(b)可以看出，焦炭顆粒燃燒過程中有一個明顯的由外向內(nèi)的轉(zhuǎn)化趨勢，最外層碳基最先轉(zhuǎn)化，但是最內(nèi)層碳基轉(zhuǎn)化速率峰值最高，也就意味著越靠近顆粒核心，反應(yīng)性越強，這和顆粒的結(jié)構(gòu)參數(shù)有關(guān)。圖4(c)、(d)表明，燃燒初期CHARH的轉(zhuǎn)化速率變化最快。隨著燃燒進行，顆粒內(nèi)部氧氣濃度上升，CHARC的轉(zhuǎn)化速率也開始增加，并在焦炭轉(zhuǎn)化率為0.4時達到了3種碳基中最大峰值。CHARG的反應(yīng)性是3種碳基中最差的。Maffei等[17]和Gibbins等[29]的實驗證明了以上結(jié)論。燃燒初期3種碳基轉(zhuǎn)化速率的波動是因為初期顆粒內(nèi)部不同氣體濃度發(fā)生變化所產(chǎn)生的，具體見圖3。

圖5(a)、(b)分別為不同粒徑和不同環(huán)境O2濃度的條件下，焦炭顆粒轉(zhuǎn)化率隨時間的變化情況。圖5(a)表明，在環(huán)境O2濃度相同的條件下，粒徑大的焦炭顆粒燃燒情況較差，有短暫的滯燃現(xiàn)象，原因為氣體產(chǎn)生的競爭效應(yīng)阻礙燃燒。半徑為50 μm的顆粒轉(zhuǎn)化迅速，幾乎沒有滯燃。半徑為75 μm的顆粒燃燒初期存在20 ms左右的緩慢燃燒階段，且轉(zhuǎn)化完全需要接近200 ms。圖5(b)表明，在粒徑相同的條件下，環(huán)境O2濃度低的焦炭顆粒燃燒情況較差。環(huán)境O2濃度為20%時，顆粒的滯燃階段超過25 ms，轉(zhuǎn)化較慢，且在200 ms內(nèi)沒有達到完全轉(zhuǎn)化。環(huán)境O2濃度為40%時，焦炭顆粒有10 ms左右的滯燃階段，轉(zhuǎn)化速率最大。結(jié)合圖5(a)、(b)可以看出，環(huán)境O2濃度對于焦炭顆粒轉(zhuǎn)化率及轉(zhuǎn)化速率的影響要大于粒徑的影響，然而顆粒粒徑的變化對于焦炭顆粒燃燒初期的滯燃現(xiàn)象影響更大。

圖5 不同工況下焦炭顆粒轉(zhuǎn)化率的變化

5 結(jié) 論

本工作結(jié)合多孔介質(zhì)傳熱傳質(zhì)模型、焦炭燃燒化學反應(yīng)模型和氣固反應(yīng)隨機孔轉(zhuǎn)化模型在內(nèi)的焦炭顆粒燃燒的物理模型，通過FORTRAN語言的編譯，對焦炭顆粒燃燒的過程，內(nèi)部氣體生成及傳輸?shù)奶匦赃M行了數(shù)值模擬研究。得出以下結(jié)論。

（2）顆粒內(nèi)部的氣體濃度受到化學反應(yīng)與物理擴散的共同影響，所產(chǎn)生的競爭效應(yīng)會使焦炭顆粒初期的燃燒極不穩(wěn)定，形成滯燃現(xiàn)象，抑制焦炭顆粒的燃燒。減小焦炭顆粒的粒徑能夠明顯改變滯燃現(xiàn)象，原因在于粒徑縮小能夠加強物理擴散的控制力，削弱競爭效應(yīng)。提升O2濃度對于滯燃現(xiàn)象影響不大，但是能夠增大化學反應(yīng)速率，從而促進焦炭顆粒的燃燒。

（3）焦炭顆粒轉(zhuǎn)化過程本質(zhì)上是3種碳基轉(zhuǎn)化的綜合體現(xiàn)，CHARH含量高的焦炭在燃燒初期會體現(xiàn)出較好的反應(yīng)性，適用于前期助燃的煤焦所需。CHARC含量高對于焦炭整體反應(yīng)性的提高有所幫助，適用于作為燃燒主要燃料。CHARG含量會對焦炭燃燒產(chǎn)生負面影響。在選煤時應(yīng)注意3種碳基的相對含量。

符 號 說 明

Ai——第i步反應(yīng)的指前因子，mol·(m2·s)-1 a，b，c——3種存在形式在焦炭顆粒中各自所占的百分比，具體數(shù)值與煤種有關(guān)，本工作中假設(shè)3種存在形式比例相同，即abc0.33，其余為灰分，不參與反應(yīng) CCO——顆粒內(nèi)部CO體積分數(shù) ——顆粒內(nèi)部CO2體積分數(shù) Cj——顆粒內(nèi)部j氣體體積分數(shù) Cjb——環(huán)境中j氣體體積分數(shù) Cj0——顆粒內(nèi)部j氣體初始體積分數(shù) ——顆粒內(nèi)部O2體積分數(shù) Cp——顆粒有效比熱容，kJ·(kg·K)-1 CHARC——碳基在焦炭顆粒中3種存在形式之一，為無序結(jié)構(gòu)碳基 CHARG——碳基在焦炭顆粒中3種存在形式之一，為有序結(jié)構(gòu)碳基 CHARH——碳基在焦炭顆粒中3種存在形式之一，為碳氫比為2:1的化合物 Dbj——顆粒內(nèi)部j氣體的容積擴散系數(shù)，m2·s-1 Dej——顆粒內(nèi)部j氣體的有效擴散系數(shù)，m2·s-1 Dkj——顆粒內(nèi)部j氣體的努森擴散系數(shù)，m2·s-1 Ei——第i步反應(yīng)所需的活化能，kJ·mol-1 ΔHi——第i步反應(yīng)所需的反應(yīng)熱，kJ·mol-1 h——孔隙內(nèi)混合氣體的對流傳熱系數(shù)，W·(m2·K)-1 Kj——顆粒內(nèi)部j氣體的對流傳質(zhì)系數(shù)，m·s-1 L0——顆粒孔隙長度，m Mj，Mj'——分別為氣體j、j'相對分子質(zhì)量 m——碳基反應(yīng)速率重疊系數(shù)，本工作取1 n——反應(yīng)級數(shù) R——顆粒半徑，m RCHARH，RCHARC，RCHARG——焦炭顆粒的3種存在形式CHARH、CHARC、CHARG的化學反應(yīng)消耗速率，mol·(m2·s)-1 Ri——單種化學反應(yīng)的反應(yīng)速率，i1,2,,9 Rj——單種氣體的化學反應(yīng)速率，j為O2、CO2、CO RO2，RCO2，RCO——顆粒內(nèi)部反應(yīng)氣體O2、CO2、CO的化學反應(yīng)生成速率，mol·(m2·s)-1 Rs——單種碳基的化學反應(yīng)速率，mol·(m2·s)-1 RV——摩爾氣體常量，kJ·(mol·K)-1 r——沿顆粒徑向長度，m rp——顆?？讖?，m S——顆粒比表面積，m2·kg-1 S0——顆粒初始比表面積，m2·kg-1 T——顆粒溫度，K Tb——環(huán)境溫度，K T0——顆粒初始溫度，K t——反應(yīng)時間，s V——顆粒體積，m3 Vj，Vj'——分別為氣體j、j'在正常沸點時液態(tài)克摩爾容積，m3·(kg·mol)-1 X——焦炭顆粒的轉(zhuǎn)化率 Xs——單種碳基的轉(zhuǎn)化率 ε——顆?？紫堵?ε0——顆粒初始孔隙率 λ——顆粒熱導率，kJ·(m·K)-1 ρ——顆粒固相密度，kg·m-3 σ——孔結(jié)構(gòu)參數(shù)，σRS0/(1ε0) τ——顆?？紫队厍?τs——單種碳基轉(zhuǎn)化的量綱1時間，τsRsCnS0t/(1ε0) τ0——顆粒初始孔隙迂曲度 ——結(jié)構(gòu)參數(shù)， 下角標 i——化學反應(yīng)序號， j——O2、CO2、CO 3種氣體中任意一種 j'——O2、CO2、CO 3種氣體中除j以外另外兩種氣體的混合氣體 s——CHARH、CHARC、CHARG 3種碳基中的任意一種

References

[1] Hadjipaschalis I, Kourtis G, Poullikkas A. Assessment of oxyfuel power generation technologies [J]. Renewable and Sustainable Energy Reviews, 2009, 13 (9): 2637-2644.

[2] Toftegaard M B, Brix J, Jensen P A,. Oxy-fuel combustion of solid fuels [J]. Progress in Energy and Combustion Science, 2010, 36 (5): 581-625.

[3] 劉慧敏, 王春波, 黃星智, 等. 富氧燃燒方式下煤中砷的揮發(fā)行為 [J]. 化工學報, 2015, 66 (12): 5079-5087. LIU H M, WANG C B, HUANG X Z,. Volatilization of arsenic in coal during oxy-fuel combustion [J]. CIESC Journal, 2015, 66 (12): 5079-5087.

[4] 劉國偉, 董芃, 別如山. 富氧條件下煤粉氣流的著火性能變化 [J]. 化工學報, 2013, 64 (7): 2596-2603. LIU G W, DONG P, BIE R S. Ignition characteristics of coal powder airflow under various oxygen-enriched conditions [J]. CIESC Journal, 2013, 64 (7): 2596-2603.

[5] 段翠九, 譚力, 趙科, 等. 0.15MW循環(huán)流化床富氧燃燒試驗研究 [J]. 中國電機工程學報, 2012, 32 (S1): 138-142. DUAN C J, TAN L, ZHAO K,. Experimental study on combustion at oxygen-enriched atmosphere in 0.15 MW circulating fluidized bed [J]. Proceedings of the Chinese Society for Electrical Engineering, 2012, 32 (S1): 138-142.

[6] Struis R P W J, Scala C V, Stucki S,. Gasification reactivity of charcoal with CO2(Ⅰ): Conversion and structural phenomena [J]. Chemical Engineering Science, 2002, 57 (17): 3581-3592.

[7] Bejarano P A, Levendis Y A. Single-coal-particle combustion in O2/N2and O2/N2environment [J]. Combustion and Flame, 2008, 153: 270-287.

[8] ErnestoA L, Amundson N R. An analytical study of single particle char gasification [J]. AIChE J., 1978, 24 (1): 72-87.

[9] Morell J I, Amundson N R, Park S K. Dynamics of a single particle during char gasification [J]. Chemical Engineering Science, 1990, 45 (2): 387-401.

[10] Yamashita T, Fujii Y, Morozumi Y,. Modeling of gasification and fragmentation behavior of char particles having complicated structures [J]. Combustion and Flame, 2006, 146 (1/2): 85-94.

[11] Bhatia S K, Perlmutter D D. A random pore model for fluid-solid reactions (Ⅰ): Isothermal, kinetic control [J]. AIChE J., 1981, 26 (3): 379-385.

[12] Bhatia S K, Perlmutter D D. A random pore model for fluid-solid reactions (Ⅱ): Diffusion and transport effects [J]. AIChE J., 1981, 27 (2): 247-254.

[13] 陳明磊, 李保衛(wèi), 武文斐. 焦炭顆粒在不同控制區(qū)域中的燃燒特性 [J]. 過程工程學報, 2014, 14 (2): 291-295. CHEN M L, LI B W, WU W F. Combustion characteristics of char particles in different controlling regions [J]. The Chinese Journal of Process Engineering, 2014, 14 (2): 291-295.

[14] 代濤, 龔志軍, 李保衛(wèi), 等. 褐鐵礦顆粒低溫CO磁化還原焙燒的實驗研究與數(shù)值模擬 [J] . 過程工程學報, 2014, 14 (4): 624-630. DAI T, GONG Z J, LI B W,. Experiment and mathematical simulation on magnetization reduction roasting of limonite particles with CO at low temperature [J]. The Chinese Journal of Process Engineering, 2014, 14 (4): 624-630.

[15] Lu G Q, Do D D. A kinetic study of coal reject-derived char activation with CO2, H2O, and air [J]. Carbon, 1992, 30 (1): 21-29.

[16] Lu G Q, Do D D. Comparison of structural models for high-ash char gasification [J]. Carbon, 1994, 32 (2): 247-263.

[17] Maffei T, Khatami R, Pierucci S,. Experimental and modeling study of single coal particle combustion in O2/N2and oxy-fuel(O2/CO2) atmospheres [J]. Combustion and Flame, 2013, 160 (11): 2559-2572.

[18] Krishna R, Wesselingh J A. The maxwell-stefan approach to mass transfer [J]. Chemical Engineering Science, 1997, 52 (6): 861-911.

[19] Yu J, Lucas J A, Wall T F. Formation of the structure of chars during devolatilization of pulverized coal and its thermoproperties: a review [J]. Progress in Energy and Combustion Science, 2007, 33 (2): 135-170.

[20] Gomez-Barea A, Leckner B, Perales A L V,. Analytical solutions of sharp interface models withth order kinetics. Application to char conversion [J]. Chemical Engineering Journal, 2012, 183 (4):408-421.

[21] Wu S, Gu J, Zhang X,. Variation of carbon crystalline structures and CO2gasification reactivity of Shenfu coal chars at elevated temperatures [J]. Energy and Fuels, 2008, 22 (1):199-206.

[22] Oberlin A. 2. Carbonization and graphitization as studied by transmission electron microscopy (Charles E. Pettinos awrad) [J]. Carbon, 1984, 22 (2): 203.

[23] HURT R H, Calo J M. Semi-global intrinsic kinetics for char combustion modeling [J]. Combustion and Flame, 2001, 125 (3): 1138-1149.

[24] Niksa S, Liu G S, Hurt R H. Coal conversion submodels for design applications at elevated pressures (Ⅰ): Devolatilization and char oxidation [J]. Progress in Energy and Combustion Science, 2003, 29 (5): 425-477.

[25] Senneca O, Cortese L. Kinetics of coal oxy-combustion by means of different experimental techniques [J]. Fuel, 2012, 102 (6): 751-759.

[26] Hecht E S, Shaddix C R, Geier M,. Effect of CO2and steam gasification reactions on the oxy-combustion of pulverized coal char [J]. Combustion and Flame, 2012, 159 (11):3437-3447.

[27] Maffei T. Kinetic of coal combustion [D]. Milano: Politecnico di Milano, 2013.

[28] Everson R C, Neomagus H W J P, Kaitano R. The random pore model with intraparticle diffusion for the description of combustion of char particles derived from mineral- and inertinite rich coal [J]. Fuel, 2011, 90 (7): 2347-2352.

[29] Russell N V, Gibbins J R, Man C K,. Coal char thermal deactivation under pulverized fuel combustion conditions [J]. Energy & Fuels, 2002, 43 (4): 57.

Establishment of various char-RPM under O2/CO2atmosphere

QU Jian1,2, LI Baowei1, ZHENG Kuncan1,2, WU Wenfei1,2

(1Key Laboratory of Integrated Exploitation of Bayan Obo Multi-Metal Resources,Inner Mongolia University of Science & Technology, Baotou 014010,Inner Mongolia, China;2School of Environment & Energy,Inner Mongolia University of Science & Technology, Baotou 014010,Inner Mongolia, China)

The traditional random pore model is based on a simple one step reaction, which is not suitable to describe the complex gas-solid reaction of coke particles in O2/CO2atmosphere. Based on the various carbons in coke and the characteristics of single char combustion in O2/CO2atmosphere, the various char-RPM and the model of pore structure under the complex gas-solid reaction were established. The combustion process of char particles with a diameter of 100 μm in O2/CO2atmosphere was simulated, and the results were calculated and analyzed by FORTRAN. The research showed that the particles in the early stage of combustion had a competitive effect and the gas concentration in the pore was fluctuated. The reason for fluctuation was the competition between chemical reaction and physical diffusion, which can be improved by increasing the O2concentration or decreasing the particle size. The various char-RPM proposed in this work had a good adaptability to characterizing the combustion characteristics of char particles in O2/CO2atmosphere.

O2/CO2atmosphere; porous media; random pore model; single particle char combustion; mass transfer; mesoscale; numerical simulation

2016-04-22.

10.11949/j.issn.0438-1157.20160532

TQ 028.8

A

0438—1157（2016）10—4468—09

內(nèi)蒙古應(yīng)用技術(shù)研發(fā)與開發(fā)基金資助項目（20130310）；內(nèi)蒙古高校創(chuàng)新團隊研究計劃資助項目（NMGIRT1406）。

2016-04-22收到初稿，2016-06-06收到修改稿。

聯(lián)系人：武文斐。第一作者：曲踐（1990—），男，碩士研究生。

WU Wenfei, wwf@imust.com