活性炭固定床吸附苯蒸汽的數(shù)值模擬

宋穎韜，黨明巖，王鼎雯

(沈陽理工大學 環(huán)境與化學工程學院，沈陽 110159)

活性炭固定床吸附苯蒸汽的數(shù)值模擬

宋穎韜，黨明巖，王鼎雯

(沈陽理工大學 環(huán)境與化學工程學院，沈陽 110159)

以守恒定律為基礎，結合吸附平衡關系和傳質(zhì)速率方程，建立了活性炭吸附苯蒸汽固定床的數(shù)學模型，采用有限差分法對模型方程進行離散，并借助于Matlab軟件編程求解，考察了傳質(zhì)系數(shù)、進料濃度、空塔氣速和溫度的變化對吸附過程的影響。計算結果表明，采用本文建立的數(shù)學模型可以很好地預測吸附床的負荷變化和穿透曲線，可為吸附床的設計和最優(yōu)化操作提供理論指導。

活性炭；固定床；吸附；數(shù)值模擬

揮發(fā)性有機物對環(huán)境及人體健康危害極大，目前控制揮發(fā)性有機化合物的排放已成為環(huán)境保護的熱點問題。吸附作為一種揮發(fā)性有機物回收技術具有凈化效率高、處理量大、能耗低、設備簡單等優(yōu)點，在工業(yè)中應用廣泛[1-3]。工業(yè)上活性炭動態(tài)吸附有機化合物的設備主要是固定床，吸附操作時流體流過固定床吸附劑顆粒層，床層中吸附劑的吸附量隨著操作過程的進行而逐漸增加，同時床層內(nèi)各處濃度和溫度分布也隨時間而變化，相比于一般的固定床反應器，吸附固定床操作的影響因素更加復雜。鑒于動態(tài)吸附過程的復雜性和各種參數(shù)的耦合性，如果只通過實驗方法對其進行研究，既耗時又不經(jīng)濟，而數(shù)值模擬是對實驗研究的有力補充。國內(nèi)外對于吸附固定床的研究主要集中在高效吸附劑的開發(fā)、吸附性能的實驗研究等[4-8]，由于模型的建立及求解的困難，對于固定床吸附器的模擬方面則略顯不足。揮發(fā)性有機物中比較重要的有苯、甲苯、二甲苯、丙酮等，本文以空氣—苯體系為例，采用數(shù)值模擬的方法研究活性炭對苯的動態(tài)吸附過程，以此為實際吸附過程的操作和設計提供一定的理論依據(jù)。

1 吸附固定床數(shù)學模型的建立和求解

固定床吸附過程比較復雜，影響因素較多，建立數(shù)學模型時首先假設吸附床的軸向氣體流速不變，徑向速度、壓力和濃度分布忽略不計；吸附過程為等溫過程；軸向擴散系數(shù)為常數(shù)。在吸附床軸向上取dz長度的微元體，以其為衡算對象，建立物料衡算方程如下：

(1)

式中：C為氣相中苯的濃度，mol/m3；q為吸附量，mol/kg；u為空塔氣速，m/s；Dz為軸向擴散系數(shù)，m2/s；ρb為床層的堆積密度，kg/m3；ε為床層的空隙率。

目前吸附模擬中使用最為廣泛、較為可靠的吸附速率模型是線性動力模型[9]，采用該模型描述吸附速率的方程為

(2)

式中：k為傳質(zhì)系數(shù)，s-1；q?為平衡吸附量，mol/kg。

活性炭吸附苯的平衡關系符合Langmuir吸附等溫式，即

(3)

式中：qmax為最大吸附量，mol/kg；B為Langmuir吸附平衡常數(shù)，m3/mol。

根據(jù)吸附過程的特征，采用初始條件和邊界條件分別為

t=0,C=C0,q=0

(4)

(5)

(6)

由式(1)～(6)組成的微分方程組即為吸附床的數(shù)學模型。

由于待求解的模型為偏微分方程組，本文采用有限差分法進行離散。借助于Matlab軟件對離散方程進行編程求解。模型參數(shù)主要參考文獻中提供的實驗數(shù)據(jù)[10]。

2 計算結果和分析

2.1 吸附床層內(nèi)吸附質(zhì)濃度變化

在溫度T=293.15K、苯進料濃度Cin=0.5mol/m3、床層空隙率ε=0.35、吸附時間t=300min、空塔氣速u=0.3m/s、傳質(zhì)系數(shù)k=1.5×10-3s-1條件下，模擬計算得到吸附床內(nèi)苯濃度和吸附量的動態(tài)分布規(guī)律見圖1和圖2，圖中軸向長度和濃度均采用無量綱數(shù)值。由圖1和圖2可見，隨著吸附的進行，在同樣的床層位置處苯濃度逐漸增大，相應地吸附量也逐漸增大，而在吸附時間相同時，沿床層長度方向苯濃度和吸附量逐漸減小。從圖中可明顯看到傳質(zhì)區(qū)的形成過程，吸附初始階段固定床出口附近苯濃度和吸附量均為0，隨著傳質(zhì)區(qū)向出口方向移動，一段時間后吸附床被穿透，床層出口處苯濃度逐漸增大，直至傳質(zhì)區(qū)移出固定床，此時固定床出口處苯濃度接近進料濃度，整個固定床達到吸附飽和，失去吸附能力。以二維圖形表示出吸附負荷曲線和透出曲線如圖3和圖4所示，圖3表示出了吸附時間分別為30min、50min和80min時吸附負荷曲線的變化，可明顯看到吸附曲線隨著時間的變化，當曲線下端到達床層末端時，開始有吸附質(zhì)的流出。從圖4可見，當吸附時間為100min左右時，床層被穿透，此時計算得到床層末端穿透吸附量為0.18mol/kg，此后流出物中苯濃度逐漸增加，到220min時床層出口濃度接近于入口濃度，此時飽和吸附量達到2.86mol/kg。

2.2 傳質(zhì)系數(shù)對床層內(nèi)吸附質(zhì)濃度分布的影響

為簡化計，本文計算中將傳質(zhì)系數(shù)k視為常數(shù)，而實際生產(chǎn)情況下影響k的因素多而復雜，其數(shù)值需要大量的實驗數(shù)據(jù)擬合確定，本文以定性分析為主，僅考察單因素變化的影響。分別取傳質(zhì)系數(shù)k為0.8×10-3s-1、1.5×10-3s-1、2.5×10-3s-1，其它參數(shù)同2.1，計算得到吸附負荷曲線和透出曲線見圖5和圖6。從圖中可以看出，傳質(zhì)系數(shù)越大，傳質(zhì)阻力越小，傳質(zhì)速率越大，吸附越快，故傳質(zhì)區(qū)越短，濃度分布曲線越陡。當傳質(zhì)系數(shù)較小時，傳質(zhì)區(qū)較長，此時若采用的床層較短，則吸附床就會很快被穿透，這樣會造成大量的吸附空間不能利用。因此吸附床的操作條件若發(fā)生變化，床層高度可能會成為影響吸附效果的因素，這在設計和生產(chǎn)中尤為重要，采用數(shù)值模擬的方法可以方便地預測這些影響的結果。模擬計算結果顯示，當傳質(zhì)系數(shù)從0.8×10-3s-1增加到2.5×10-3s-1時，飽和吸附量由2.88mol/kg變化到2.86mol/kg，說明傳質(zhì)系數(shù)的變化對飽和吸附量幾乎沒有影響。

圖1 吸附床內(nèi)氣相中苯的濃度分布

圖2 吸附床內(nèi)固相中苯的吸附量分布

圖3 不同時間下的吸附負荷曲線

圖4 吸附床透出曲線

圖5 不同傳質(zhì)系數(shù)下的吸附負荷曲線(t=60min)

圖6 不同傳質(zhì)系數(shù)下的透出曲線

2.3 進料濃度對床層內(nèi)吸附質(zhì)濃度分布的影響

其它條件同2.1，分別取苯進料濃度Cin=0.35mol/m3、0.5mol/m3、0.65mol/m3，計算得到進料濃度對固定床吸附過程的影響見圖7和圖8，為更好地表達吸附結果，圖8采用以吸附量表示的透出曲線。從圖中曲線分布可以看出，進料濃度越高，吸附床相同位置處苯的濃度越大，負荷曲線越靠近右側。從透出曲線可見，進料濃度越高，穿透時間越短，飽和吸附量越大。這是因為，進料濃度增大，相同時間內(nèi)進入吸附床的苯量也增大，傳遞推動力也增加，因此越快達到吸附平衡，使得更多的苯進入下一段床層，從而穿透時間變短。模擬計算結果顯示，當進料濃度從0.35mol/m3增加到0.65mol/m3時，飽和吸附量由2.56mol/kg增加到3.07mol/kg。

圖7 不同進料濃度下的吸附負荷曲線(t=60min)

圖8 不同進料濃度下的透出曲線

2.4 空塔氣速對床層內(nèi)吸附質(zhì)濃度分布的影響

其它條件同2.1，分別取空塔氣速u為0.2m/s、0.3m/s、0.4m/s，模型計算結果見圖9和圖10。

圖9 不同空塔氣速下的吸附負荷曲線(t=60min)

從圖中可以看出，空塔氣速對吸附過程影響很大，文中數(shù)學模型的左邊界采用第三類邊界條件，故在入口處苯濃度隨空塔氣速的不同而不同，空塔氣速增加，入口處苯濃度也會有所增加。此外，空塔氣速越大，苯濃度分布曲線越平緩，同時由于空塔氣速增加，即負荷曲線沿吸附床軸向移動速度加快，同樣床層長度下，吸附床越快被穿透。同時由于流速的增大導致氣體與吸附劑接觸時間減少，苯來不及被吸附就進入下一段床層。從圖中還可明顯看到，當空塔氣速較大時，傳質(zhì)區(qū)較長，此時若吸附床較短，也會造成吸附劑利用率低的結果。實際生產(chǎn)中流速加快會使得傳質(zhì)系數(shù)也增加，當增加到一定程度時，會造成吸附過程的控制步驟發(fā)生變化，使得吸附過程更為復雜。

圖10 不同空塔氣速下的透出曲線

2.5 溫度對床層內(nèi)吸附質(zhì)濃度分布的影響

其它條件同2.1，在溫度T分別取293.15K、303.15K、313.15K時，固定床吸附過程模擬計算結果如圖11和圖12所示。可以看出，在相同條件下，溫度越高時，傳質(zhì)區(qū)越長，傳質(zhì)區(qū)末端越靠近吸附床出口，床層越快被穿透。從靜態(tài)吸附分析結果可知活性炭吸附苯是放熱過程[10]，隨著溫度的升高，一方面吸附速率有所增大，另一方面物理吸附過程減弱，兩種影響共存下導致傳質(zhì)區(qū)位置的不同。因為文中建立數(shù)學模型時假設床層內(nèi)溫度不變，不考慮吸附熱引起的溫度效應，按該數(shù)學模型計算時，在一定溫度差異范圍內(nèi)，溫度對吸附過程影響不是很大，所以濃度曲線總體差異不大。但從透出曲線來看，溫度較高時對吸附過程的影響較大，計算結果顯示，溫度由293K升溫至313K時飽和吸附量由2.87mol/kg降低到2.56mol/kg，這與活性炭對苯的吸附為放熱過程相符合。

圖11 不同溫度下的吸附負荷曲線(t=60min)

圖12 不同溫度下的透出曲線

3 結論

(1)通常條件下，隨著吸附的進行，同樣床層位置處的苯濃度和吸附量都逐漸增大，吸附時間相同時，沿床層軸向苯濃度和吸附量逐漸減小，模擬計算結果顯示出傳質(zhì)區(qū)的形成和移動過程。

(2)隨著傳質(zhì)系數(shù)的增大，吸附速率加快，濃度分布曲線變陡，傳質(zhì)區(qū)變短，穿透時間增加；隨著進料濃度的提高，穿透時間減少，飽和吸附量增加；隨著空塔氣速的增大，濃度分布曲線變緩，穿透時間減少；隨著溫度的升高，傳質(zhì)區(qū)變長，穿透時間變短，飽和吸附量降低。

(3)所建立的數(shù)學模型能夠較好地對固定床內(nèi)的吸附過程進行描述和預測，模擬計算結果比較合理可靠。

[1]陶宇鵬,陳中明,代茂節(jié),等.揮發(fā)性有機物(VOCs)污染治理中吸附技術的相關應用[J].低溫與特氣,2015,33(6):45-50.

[2]陸豪,吳祖良,高翔.吸附法凈化揮發(fā)性有機物的研究進展[J].環(huán)境工程,2013,31(3):93-97.

[3]席勁瑛,武俊良,胡洪營,等.工業(yè)VOCs 氣體處理技術應用狀況調(diào)查分析[J].中國環(huán)境科學,2012,32(11):1955-1960.

[4]范忠雷,查會平,王玲慧.乙二胺硅膠材料固定床吸附銅離子的研究[J].鄭州大學學報,2014,35(1):68-72.

[5]劉娟,王小藝,劉載文.殼聚糖衍生物固定床中Cu(Ⅱ)的吸附性能研究[J].北京工商大學學報:自然科學版,2009,27(2):4-7.

[6]于偉華,劉鵬,龍超,等.疏水性超高交聯(lián)吸附樹脂對氯代烴蒸汽的固定床吸附特性研究[J].環(huán)境科學,2011,32(9):2805-2809.

[7]張會平,劉劍,鄢瑛.苯在結構化固定床的吸附動力學[J].化工進展,2011,30(Z):41-45.

[8]蘇青青,楊嘉謨.活性炭固定床吸附SO2的機理及動力學研究[J].三峽大學學報:自然科學版,2010,32(5):97-101.

[9]王丹莉,劉向軍,劉應書.氣體吸附動力學模型的研究現(xiàn)狀[J].化工進展,2010,29(Z):5-11.

[10]胡祖美,張艷,王金渠,等.苯在活性炭纖維上吸附等溫線的測定及分析[J].石油學報:石油加工,2008,24(4):484-487.

(責任編輯：馬金發(fā))

Numerical Simulation of the Benzene Vapor Adsorption onto the Activated Carbon Fixed-bed

SONG Yingtao,DANG Mingyan,WANG Dingwen

(Shenyang Ligong University,Shenyang 110159,China)

The mathematical model for the benzene vapor adsorption onto the activated carbon fixed-bed has been set up on the basis of conservation law along with the adsorption equilibrium relationship and mass transfer velocity equation.The finite differential method was adopted to discrete the equation,while MATLAB programming was also used so as to investigate the effects of mass transfer coefficient,inlet material concentration,superficial gas velocity and the change of temperature on the adsorption process.The calculations show that the mathematical model developed could forecast the load change and the breakthrough curve of the adsorption bed,then provide theoretical guidance to the design and optimization of the adsorption operation.

activated carbon；fixed-bed；adsorption；numerical simulation

2016-03-07

宋穎韜(1972—)，女,副教授,研究方向：化工數(shù)值模擬。

1003-1251(2017)02-0088-05

TQ028.1

A