第25卷第3期南京工業(yè)大學(xué)學(xué)報2003年5月JOURNAL OF NANJING UNIVERSITY OF TECHNOLOGYMay 2003變壓吸附空分制氧非等溫過程模擬周漢濤,馬正飛,姚虎卿南京工業(yè)大學(xué)化學(xué)化工學(xué)院江蘇南京21000摘要∶就變壓吸附空氣分離制氧過程對接近真實情況的非線性、非等溫模型構(gòu)成的偏微分方程組采用正交配置進行空間離散化和三階半隱式龍格庫塔法的數(shù)值計箅方法研究了變壓吸附過程中床層內(nèi)溫度和濃度的動態(tài)行為考察了清洗比、吸附壓力、進氣流量、吸附時間等操作參數(shù)對過程性能的影響為過程優(yōu)化設(shè)計建立基礎(chǔ)。關(guān)鍵詞:變壓吸附;空分制氧;非等溫;模擬中圖分類號:TO028文獻標識碼:A文章編號:1671-764x2003)3-0039-06自20世紀70年代基于吸附選擇性的空分體系驟及時間分配如表1所示。商業(yè)化以來,SA制氧技術(shù)在煉鋼、焊接、醫(yī)療、制表1過程循環(huán)步驟及時間藥、發(fā)酵、造紙、環(huán)保、水產(chǎn)等領(lǐng)域都得到了廣泛的應(yīng)Table 1 Cycle sequence and step time用。變壓吸附氣體分離過程是一個動態(tài)過程描述層1ADDP其動態(tài)過程的行為是偏微分方程組構(gòu)成的模型操時間/作性能和參數(shù)之間為復(fù)雜的非線性關(guān)系實驗研究床層2PGPP時間/與數(shù)學(xué)模擬均不簡單工作量大所以變壓吸附氣體AD吸附DE均壓降DP卸壓PG清洗PPE均壓升FP原料氣沖分離技術(shù)存在著商業(yè)開發(fā)領(lǐng)先于理論研究的問題。壓。對它的理論研究遠不如變溫吸附過程深入,尤其對多組分、多吸附劑、非線性吸附等溫線、非等溫和非吸附過程流程示意圖如圖1所示。平衡吸附的變壓吸附過程的理論研究較為缺乏1-41不能適應(yīng)工業(yè)應(yīng)用的要求。采用計算機數(shù)值模擬計算的方法研究變壓吸附過程的動態(tài)行為揭示操作參數(shù)和過程性能之間的關(guān)系對于指導(dǎo)試驗工作進行過程優(yōu)化設(shè)計有著重要的意義。本文將采用非等溫模型的模擬計算來研究影響變壓吸附制氧過程性能的諸多因素考慮所用模型及數(shù)值計算方法的通用性和高效性。模擬內(nèi)容反映吸附床層內(nèi)溫度和濃度的動態(tài)行為以揭示變壓吸圖1二塔制氧變壓吸附流程示意圖附過程操作參數(shù)清洗比P/F、吸附壓力、進氣流量、Fig. I Schematic diagram of the apparatus for a two-bed O, PSa process吸附時間與氧氣純度、回收率之間的關(guān)系為變壓吸附制氧過程的優(yōu)化設(shè)計提供了基礎(chǔ)。建立用于描述變壓吸附過程的數(shù)學(xué)模型采用如下假設(shè)氣體為理想氣體流體流動模型為軸向分散數(shù)學(xué)模型活塞流模型油向流速在床層內(nèi)可變化涊忽略吸附時床層軸向壓降;傳質(zhì)速率方程可用線性推動力變壓吸附過程選取經(jīng)典的二床六步驟循環(huán),(DF模型描述等溫線模型為負載比關(guān)聯(lián)LRC舶的南京工業(yè)大學(xué)學(xué)報第25卷Langmuir吸附模型煆假設(shè)氣相和固相熱傳遞瞬時平B=k3ex(k4/To)n=ks+k6/To衡變壓吸附過程中變壓步驟的壓力呈指數(shù)型變LDF傳質(zhì)速率模型化5見式1P= Pend +((1)變壓吸附循環(huán)過程中壓力的變化情況如圖2。床內(nèi)壁的傳熱系數(shù)6和外壁的自然對流傳熱系數(shù)71分別通過下面二式關(guān)聯(lián)h;=3.404ag ega e)e-6d/D8(9軸向分散系數(shù)D1通過 Edward8l關(guān)聯(lián)式計算得到:D1=0.73Dx+1+9.7D/udn(10)00如4如熱分散系數(shù)由以下經(jīng)驗關(guān)聯(lián)式得到91:KL/Ag= Kio/Ag+oPrRe(11)本文模型采用 Danckwerts0邊界條件圖2一個循環(huán)過程床層壓力變化(1)壓和吸附階段Fig. 2 Pressure swing history during one cycled yi對床層進行物料及熱量衡算可得以下模型方程p Cpgl流動相模型(2清洗和均壓升階段DL0dtDIaz l:sl總傳質(zhì)平衡方程az.=pg Gpg TIaCl(uc)+aC(3均壓降和卸壓階段(3)氣固相能量平衡模型aTd70TEKl az2 pg Cx(ep, Cg +Pb Cs)T床層的初始狀態(tài)為潔凈床層?！?)∑T-T)=0(4)物性數(shù)據(jù)、吸附劑、床層特性參數(shù)及動力學(xué)參數(shù)床壁能量衡算方程取自文獻1112由于本文采用的模型為非等溫情形考慮了吸pC、A=2xR1h(T-Tw)-附等溫線的非線性、吸附引起的流速變化、傳質(zhì)阻2πRnh(Tw-T)(5)力傳質(zhì)及傳熱的軸向擴散、物性參數(shù)是溫度和壓力其中A=R2-R)的函數(shù)因此模型接近真實情況。負載比關(guān)聯(lián)LBC舶 Langmuir吸附平衡模型B P2計算方法(6)第2期周漢濤等變壓吸附空分制氧非等溫過程模擬41然后再應(yīng)用常微分方程組的數(shù)值求解方法進行數(shù)值再生越完全氧氣純度也越高而反過來清洗耗氧量求解。上述模型中的偏微分方程組采用正交配置的增大使得回收率降低。法13在軸向上進行離散化使之成為一系列常微分方程組。常微分方程組的數(shù)值積分采用三階半隱式Runge-Kut法去。采用正交配置的優(yōu)點在于在獲得同樣的準確度下只需要較少的變量離散點洏半隱式sRunge-Kutta法既適合求解顯式 Runge Kutta不能有效求解或不能求解的陡峭常微分方程組體系,又可3105同時求解微分代數(shù)方程組體系,計算時間較少。計00000算程序用 FORTRAN語言編寫任何復(fù)雜的PA分離過程都可以分解為一系0550600.5650.700750.800.85090095列相對簡單的步驟而且每個床層都經(jīng)歷同樣的過程步驟因此解一個吸附床層的模型方程并循環(huán)起來就可以得到整個過程的解。這種方法不受床數(shù)圖39.0( STP )min時不同壓力及P/F下的氧氣純度與回收率的限制可以用于任意床數(shù)的過程模擬。計算過程Fig.3O2 purity and recovery at2×103,3×103,4×105 Pa and vari中每一步結(jié)束時保存床層軸向的組成和溫度分布作ous P/F ratios with the feed flow rate of 9.0 I( STP )min為下一個步驟的初始值。吸附和均降步驟中不同時間離開床層的氣體組成、溫度和流速以數(shù)組的形式3.2吸附壓力保存經(jīng)過線性插值作為清洗和均升階段的進氣條圖3為不同清洗比值下壓力對氧氣的純度與回件收率的影響關(guān)系。從圖3中可看出壓力在低清洗比值下的影響比在高清洗比值的大。隨著壓力的增3結(jié)果與討論大氧氣的純度升高而回收率降低,這與Famo1的試驗和模擬結(jié)果一致?；厥章实慕档褪怯捎谖侥M計算為跟蹤床層濃度及溫度的動態(tài)行為。過程吸附氧氣的絕對量比較大但卸壓過程的損失模擬計算結(jié)果表明層中溫度從波動到達穩(wěn)態(tài)比也較多的緣故。濃度從波動到達穩(wěn)態(tài)慢溫度一般要在15個循環(huán)后圖4為不同吸附壓力下兩種清洗比時吸附階段逐漸達到穩(wěn)定而濃度在幾個循環(huán)就能達到穩(wěn)態(tài)。床層內(nèi)氧氣的軸向摩爾分數(shù)分布。從圖4上可看床層內(nèi)濃度的變化是很難測定的因此模擬床層內(nèi)的濃度變化是有意義的工作。為考察清洗比(P/F)吸附壓力、進氣流速、吸附時間等操作參數(shù)對過程性能的影響即對回收率與純度的影響氧氣純度為吸附過程出氣的體積平均值)定義回收率與清洗比如下O3回收率=(吸附過程產(chǎn)品氣氧氣總量-清洗過程耗氧總量)(原料氣沖壓和吸附階段進氣氧氣02PF=060總量)×100%=KO2)P/F=清洗過程耗氧量/吸附過程進氣氧總量Z/L3.1清洗比清洗比P/F對氧氣純度與回收率的影響見圖圖49.0sm)m,0.60P/F和090P/F下吸附壓3、5、x氧氣純度用R(O2表示在圖35、7中隨力對氧氣在吸附階段軸向摩爾分數(shù)分布的影響著P/F的增加氧氣的純度上升,回收率降低。但Fig, 4 Effects of adsorption pressure on O2 concentration profiles alo42南京工業(yè)大學(xué)學(xué)報第25卷出在較高吸附壓力時傳質(zhì)區(qū)變得陡峭床層在相同在相同清洗量下純度自然升高。回收率隨吸附時間時間內(nèi)利用率提高純度也越高但在壓力大于3×的延長而增加因為在相同的清洗比下短吸附時間I03Pa時分布的差別很小,同文獻15]一致。一般所得到的產(chǎn)品氣量小當(dāng)沖壓階段的時間不變時根工業(yè)變壓吸附制氧過程的吸附壓力選在3×105Pa據(jù)回收率的定義回收率必然下降。左右的原因也在此。3.3進氣流量圖5為不同流速和清洗比下氧氣的純度與回收率的變化情況。在圖5中我們可以看出高流速的PF=090純度比低流速的純度低在13.0I( STP Ymin時純度顯著降低。另一方面氰氣的回收率卻隨流速的降低而降低因為低流速下吸附階段吸附的氮氣相11.0 L STP) hnin13.0L( STP)Amin對多那么在清洗階段需要相對多的產(chǎn)品氣清洗P/F=060流速高回收率也相對較高0.00.10203040.506圖64×105Pa,0.60P/F和0.90P/F下進氣流量對85氣在吸附階段軸向摩爾分數(shù)分布的影響Fig 6 Eifects of feed flow rate on O concentration profiles along9.0L( STP) Aminthe adsorption bed11.0 L(STP)Aminpressure of4×1030.60,0.90P/Fati0.550600.650.700750.80085090095圖54×105Pa時不同進氣流量及P/F下氧氣純度與回收率AD step 30Fig 5 O2 purity and recovery between 9.0 I( STP)min11.0 I( STP )min, 13. 0 I( STP )min at various♂xP/F ratios under 4x 10 Pa85506006507075080085309095圖6為相同壓力下不同流速時吸附階段床層內(nèi)氧氣的軸向摩爾分數(shù)分布?？梢钥闯鲅鯕獾膫髻|(zhì)區(qū)隨著流速的升高而變寬,即吸附階段產(chǎn)品氣隨著圖74×105Pa90( STP )min不同P/F條件下流速的增加床層尾端污染加重純度下降。在低流兩種吸附時間的氧氣純度與回收率速下吸附階段床層尾端相對潔凈但是由于氧氣摩Fig 7 O2 purity and recovery between two different adsorptio爾分數(shù)高在隨后步驟中損失更多回收率下降。在step times at 4 10 Pa 9.0 I( STP ymin, and various高清洗比下床層尾端氧氣摩爾分數(shù)更高回收率下P/F ratios降3.4吸附時間上述模擬計算考察了清洗比(P/F入吸附壓如圖7為不同吸附時間對氧氣純度與回收率的力進氣流速吸附時間等主要參數(shù)對過程性能的影影響?？梢娂兌入S吸附時間的增加而降低回收率響可見影響過程的因素較多是一復(fù)雜的非線性關(guān)隨著吸附時間的延長而升高。原因在吸附時間縮短系。通過計算機模擬計算的結(jié)果來了解各變量之第2期周漢濤等變壓吸附空分制氧非等溫過程模擬T氣相溫度T.床壁溫度結(jié)論T環(huán)境溫度氣體流速x氧氣摩爾分率y氣相組分摩爾分率本文應(yīng)用非等溫模型模擬計算空分制氧變壓吸床層軸向坐標密度床層空隙率E1床層總空隙率附過程研究過程的特征與性能得到如下結(jié)論λ。氣體導(dǎo)熱系數(shù)λ。吸附劑導(dǎo)熱系數(shù)4.1模擬了變壓吸附過程床層中氧氣濃度與溫度H氣體粘度Re雷諾準數(shù)的動態(tài)變化。濃度達到穩(wěn)態(tài)比溫度快濃度達到穩(wěn)P普朗特準數(shù)態(tài)大約需要10個循環(huán)而溫度在15個循環(huán)后變化才減小參考文獻4.2隨著清洗比(P/F)的增加氧氣純度上升回[1李化治李寶善曾鳳鳴變壓吸附富氧過程的數(shù)學(xué)模拗J北收率下降在高清洗比(P/F)時純度變化很小回京科技大學(xué)學(xué)報1995172):110-114收率持續(xù)下降;[2]馮晉晢張玉文陳流芳等,FA法分離空氣基本原理和模型[J].低溫工程997,427-324.3隨著壓力的增大氧氣的純度升高而回收率降[3]楊春育馬寧寧焦玉海簿變壓吸附空氣分離過程的動態(tài)模低。壓力在低清洗比(P/F)下的影響比在高清洗攬J石油煉制與化工1998291)9.53比P/F)在壓力大于4×105Pa后純度基本不【4]楊春育馬寧寧焦玉海箐變壓吸附空氣分離過程中變壓步驟的數(shù)學(xué)模擬J]石油化工A9982X4)335-344.4氧氣的純度隨進氣流量增大而減小回收率卻5] Farooq S, Rathor M N, Hidajat K. A predictive model for a kineticallpressure swing adsorption separation proces[ J ]. Chem Er隨流速的增加而升高484129-41414.5純度隨吸附時間的增加而降低回收率隨著吸[6]袁渭康化學(xué)工程手冊—一反應(yīng)動力學(xué)及反應(yīng)M]北京北附時間的延長而升高學(xué)工業(yè)出版社1996本文的工作為變壓吸附制氧過程的優(yōu)化設(shè)計提[7]陳學(xué)俊化學(xué)工程手冊一傳熱及傳熱設(shè)備M北京北學(xué)工供了基礎(chǔ)。業(yè)出版社1996[8] Edwards M F, Richardson J F. Gas dispersion in packed bed[ J]Chem Eng Sci 1968 23( 1): 109-123符號說明[9] Suzuki M. Adsorptionind M ]. Elsevier, Amsterdam 1990a變壓時間常數(shù)床壁截面積[10 Danckwerts P V, Continuous flow systems distribution of residenceB等溫線模型參數(shù)c氣相濃度time J]. Chem Eng Sci,, 1953,X 1): 1-13C氣相總濃度C氣相熱容[11]時鈞汪家鼎余國琮籌化學(xué)工程手冊化工基礎(chǔ)數(shù)摭M]C.固相熱容C。塔壁熱容北京北學(xué)工業(yè)出版社,996d顆粒直徑D1床層內(nèi)直徑12 Jee JG, Lee JS, Lee C H. Air separation by a small-scale two-bedD3分子擴散系數(shù)D軸向擴散系數(shù)medical O2 pressure swing adsorptior[ J ]. Ind Eng Chem Res, 20h;對流傳熱系數(shù)h自然對流傳熱系數(shù)403647-3658k傳質(zhì)系I 13 1 Sristava R K, Joseph B. Simulation of packed bed separation processes等溫線模型參K1熱軸向擴散系數(shù)using orthogonal collocation[ J ] Computers Chemical EngineerinL床層長度n等溫線模型參數(shù)P床層總壓組分平衡吸附量14 I Farooq S, Ruthven D M, Boniface H A. Numerical simuq:組分固相平均濃度qm組分飽和吸附量orption oxygen unit J 1. Chem Eng Sci2809-2816R床層內(nèi)徑R床層外徑R氣體常數(shù)t吸附時間15 Rege S U, Yang R T, Limits for air separation by adsorption with LiXzeolit J]. Ind Eng Chem Res 1997 36 5 358-5 365南京工業(yè)大學(xué)學(xué)報第25卷Simulation of non-isothermal pressure swing adsorptionfor separating oxygen from airZHOU Han-tao MA Zheng-fei, YAO Hu-qingCollege of Chemistry and Chemical Engineering Nanjing University of Technology, Nanjing 210009, ChinaAbstract:A more realistic mathematical model, in which axial dispersion, heat effects, mass transferresistancetion of velocity and nonlinear of isothermal were taken into account was developed to the pressure swing adsorption of aiseparation for oxygen. The orthogonal collocation method and third order semi-implicitKutta integration routinewere used to solve the resulting model of partial differential equations. The dynamic behavior of temperature and concen-tration in the bed and effects of operating parameters on performance of process were discussed including purge/feed ratio, adsorption pressure, feed flow rate and adsorption time. The simulation results are helpful in the process design andKey words: pressure swing adsorption oxygen air non-isothermal model simulation簡訊《生物加工過程》第一屆編委會召開《生物加工過程》是由國家新聞出版總署批準出版南京工業(yè)大學(xué)主辦、中國化工學(xué)會生物化工專業(yè)委員會協(xié)辦的技術(shù)類全國性期刊第一屆編委會于2003年4月4日在江浦龍華大酒店由主編歐陽平凱院土主持召開。參加這次會議的編委為來自清華大學(xué)、浙江大學(xué)、北京化工大學(xué)、華東理工大學(xué)、深圳國家生化中心、中科院有關(guān)研究所等國內(nèi)生物化工知名學(xué)者、專家。專家們就刊物的定位、辦刊宗旨、辦刊欄目、刊物發(fā)展等問題展開了熱烈討論?！渡锛庸み^程》將于2003年5月份創(chuàng)刊發(fā)刊詞由中國科學(xué)院外籍院土、美國國家工程院院士黎念之博士撰寫。(學(xué)報編輯部)