河北科技大學(xué)反應(yīng)工程課設(shè)

1、化學(xué)反應(yīng)工程課程設(shè)計(jì)任務(wù)書某日產(chǎn)1000噸合成氨的裝置中的中溫變換反應(yīng)器在30.5atm下操作，進(jìn)口氣體流量NTO=1O600kmol/h，采用C9*9mm催化劑，堆密度1500kg/m3，催化床直徑為4m，高度 為1.5m，催化劑有效因子為0.6。組份%0,0*20,0%02,0，H2,0，N2M摩爾分?jǐn)?shù)0.09900.35550.05880.34350.14328300已知本征動(dòng)力學(xué)方程:mol/(h.g)竺CO = k p p 一0.5(1 PCO 七 ) M%。=3 x 106 x exp(24900)R Tmol/(h.g.atm0.5)要求:1、采用龍格庫塔方法求微分方程的數(shù)值解;

2、2、進(jìn)口溫度的初值可設(shè)定為346.5，在初始溫度與400之間取10個(gè)溫度點(diǎn)進(jìn)行計(jì)算， 繪制（1）在不同進(jìn)口溫度下床層溫度隨床層高度的變化曲線；（2）床層出口溫度隨入 口溫度變化曲線。并對(duì)結(jié)果進(jìn)行討論。3、以給定的進(jìn)口流量為中心，選擇10個(gè)不同進(jìn)口流量，間隔為100kmol/h，繪制（1） 在不同進(jìn)口流量下床層溫度隨床層高度的變化曲線；（2）床層出口溫度隨進(jìn)口流量變化 曲線。并對(duì)結(jié)果進(jìn)行討論。4、可用VB程序或EXCEL進(jìn)行計(jì)算5、比熱、熱效應(yīng)、平衡常數(shù)可以查閱化學(xué)工程手冊(cè)1.1工藝流程工藝流程設(shè)計(jì)依據(jù)工藝流程設(shè)計(jì)的依據(jù)，首先是原料氣中CO的含量。CO含量高則先采用中（高） 溫變換，因?yàn)橹校ǜ撸?/p>

3、變催化劑操作溫度范圍較寬，而且價(jià)廉易得，壽命長，大多數(shù)合 成氨原料氣中CO均高于10%，故都先通過中（高）變除去大部分CO。根據(jù)系統(tǒng)反應(yīng) 溫度的升高，為使催化劑在允許活性溫度范圍操作，對(duì)CO含量高于15%者，一般應(yīng)考 慮將反應(yīng)器分為二段或三段。其次，是根據(jù)進(jìn)入系統(tǒng)的原料氣溫度及濕含量，當(dāng)溫度及 水蒸氣含量低，則應(yīng)考慮氣體的預(yù)熱和增濕，合理利用余熱。第三是將CO變換與脫除 殘余CO的方法結(jié)合考慮。如脫除方法允許殘余CO含量較高，則僅采用中（高）變即 可，否則，可將中變與低變串聯(lián)，以降低變換氣中CO含量。工藝流程概述由于原料氣中CO含量較低，中變催化劑只需配置一段。含9.50%CO的原料氣經(jīng) 廢熱

4、鍋爐降溫，在壓力30.5atm下、溫度340C下進(jìn)入高變爐，因原料氣中水蒸氣含量 較高，一般不需添加蒸汽。經(jīng)反應(yīng)后氣體中CO含量降低，溫度升高。氣體通過高變廢 熱鍋爐，溫度下降，產(chǎn)生高壓的飽和蒸汽，由于氣體溫度尚高，一般用來加熱其它工藝 氣體而變換氣繼續(xù)被冷卻。1.2工藝條件1.2.1 壓力壓力對(duì)變換反應(yīng)的平衡幾乎沒有什么影響。但提高壓力將使析炭和生成甲烷等副反 應(yīng)易于進(jìn)行。單就平衡而言，加壓并無好處。但從動(dòng)力學(xué)角度，加壓可提高反應(yīng)速率， 如設(shè)常壓下反應(yīng)速率為r，加壓下為或可用校正系數(shù)斗r來表示壓力對(duì)反應(yīng)速率的 影響。對(duì)于一級(jí)反應(yīng)，在一定溫度下，對(duì)同一尺寸的催化劑，隨著壓力的升高，校正系 數(shù)寸

5、值增大，即反應(yīng)速率提高。當(dāng)壓力小于2MPa時(shí)，催化劑校正系數(shù)增加比較明顯， 在更高壓力下，由于內(nèi)擴(kuò)散的影響，校正系數(shù)增加比較緩慢。如果在相同的溫度和壓力下，則小顆粒的催化劑具有較大的W值。而在一定的壓 力下，對(duì)同一尺寸的催化劑，溫度升高，則寸值越低。上述情況均與內(nèi)擴(kuò)散的影響有 關(guān)。實(shí)際生產(chǎn)中，以煤為原料的合成氨廠，常壓下中變催化劑的十氣空速僅為300500h-i， 壓力為12MPa時(shí)可達(dá)8001500 h-1。以烴類為原料的大型合成氨廠，由于原料氣中CO 含量較低，壓力3MPa時(shí)空速可達(dá)25002800 h-i。從能量消耗上看，加壓是有利的。由于干原料氣的摩爾數(shù)小于干變換氣的摩爾數(shù)， 所以，先

6、壓縮原料氣后再進(jìn)行變換的能耗，比常壓變換再壓縮變換氣的能耗低。根據(jù)原 料氣中CO含量的差異，其能耗約可降低15%30%。當(dāng)然，加壓變換需用壓力較高的 蒸汽，對(duì)設(shè)備材質(zhì)的要求也較高，但綜合起來，優(yōu)點(diǎn)還是主要的。具體操作壓力的數(shù)值， 應(yīng)根據(jù)大、中、小型氨廠的不同特點(diǎn)，特別是工藝蒸汽的壓力及壓縮機(jī)各段壓力的合理 配置而定。一般小型氨廠操作壓力為0.71.2MPa，中型氨廠為1.21.8MPa，以天然氣 為原料的大型氨廠變換壓力由蒸汽轉(zhuǎn)化的壓力決定。1.2.2 溫度溫度是CO變換最重要的工藝條件。由于CO變換為放熱反應(yīng)，隨著CO變換反應(yīng) 的進(jìn)行，溫度不斷升高，速率常數(shù)k值增加的影響大于的影響，故對(duì)反應(yīng)

7、速率有利。 繼續(xù)增加溫度，二者的影響互相抵消，反應(yīng)速率隨溫度的增值為零。在提高溫度時(shí)，K P 的不利影響大于k值的增益影響，此時(shí)反應(yīng)速度會(huì)隨溫度升高而下降。對(duì)一定類型的催 化劑和一定的氣體組成而言，必將出現(xiàn)最大的反應(yīng)速率值，與其對(duì)應(yīng)的溫度，稱為最佳 溫度或最適宜溫度。由于平衡溫度隨系統(tǒng)組成而改變，不同催化劑活化能也不相同，最佳溫度Tm隨系 統(tǒng)組成與催化劑的不同而變化。對(duì)一定初始組成的反應(yīng)系統(tǒng)，隨著CO變化率x的增加， 平衡溫度Te及最佳溫度Tm均降低。對(duì)同一變換率，最佳溫度一般比相應(yīng)的平衡溫度低 幾十度。如果工業(yè)反應(yīng)器中按最佳溫度進(jìn)行反應(yīng)，則反映速率最大，即在相同的生產(chǎn)能 力下所需催化劑用量最

8、少。對(duì)中溫變換來說：（1）應(yīng)在催化劑活性溫度范圍內(nèi)操作，反應(yīng)開始溫度應(yīng)高于催化劑起始活性溫度20 C左右。（2）隨著催化劑使用年限的增長，由于中毒、老化等原因，催化劑活性降低，操 作溫度應(yīng)適當(dāng)提高。（3） 為了盡可能接近最佳溫度線進(jìn)行反應(yīng)，可采用分段冷卻。段數(shù)越多，則越接 近最佳反應(yīng)溫度線，但流程也復(fù)雜。根據(jù)原料氣中的CO含量，一般多將催化劑床層分 為一段、二段或多段，段間進(jìn)行冷卻。1.2.3 汽氣比汽氣比一般指H2O/CO比值或水蒸氣/干原料氣（摩爾比）。改變水蒸氣比例是工業(yè) 變換反應(yīng)中最主要的調(diào)節(jié)手段。增加水蒸氣用量，提高了CO的平衡變換率，從而有利 于降低CO殘余含量，加速變換反應(yīng)的進(jìn)行

9、。由于過量水蒸氣的存在，保證催化劑中活 性組分Fe3O4的穩(wěn)定而不被還原。過量的水蒸氣還起到熱載體的作用。提高水蒸氣比例，合濕原料氣中CO含量下降， 催化劑床層的溫升將減少，即改變水蒸氣的用量是調(diào)節(jié)床層溫度的有效手段。但是，水 蒸氣用量是變換過程中最主要消耗指標(biāo)，盡量減少其用量對(duì)過程的經(jīng)濟(jì)性具有重要意 義。水蒸氣比例過高，還將造成催化劑床層阻力增加，CO停留時(shí)間縮短，余熱回收設(shè) 備負(fù)荷加重等。中溫變換操作時(shí)適宜的水蒸氣比例一般為：H2O/CO=35。經(jīng)反應(yīng)后， 中變氣中H2O/CO可達(dá)15以上（水蒸氣/十氣Q0.51.0），不必再添加蒸氣即可滿足低溫 變換的要求。2動(dòng)力學(xué)模型的選取和反應(yīng)器數(shù)學(xué)

10、模型的建立2.1動(dòng)力學(xué)模型的選取根據(jù)不同原料、不同原料氣制造方法，合成氨生產(chǎn)中使用了中(高)變催化劑、低 變催化劑和耐硫變換催化劑。由于催化劑型號(hào)眾多、性能的差異以及實(shí)驗(yàn)條件的不同， 經(jīng)整理歸納后動(dòng)力學(xué)方程有數(shù)十個(gè)。在此設(shè)計(jì)計(jì)算中，選用變換反應(yīng)動(dòng)力學(xué)方程為幕函 數(shù)型動(dòng)力學(xué)方程式形式，使用比較方便，本設(shè)計(jì)中采用國產(chǎn)B系列的中變催化劑的本征 動(dòng)力學(xué)方程：dNc。= k p p-0.5 (1 co2七2 )mol / (h.g) dt t co co 2K P P式中, c f , 24900、 k = 3 x 106 x exp( )mol / (h.g.atm.0.5)TR T其中Rg = 1.

11、987 kcal /(kmol - K)2.2反應(yīng)器數(shù)學(xué)模型的建立合成氨裝置中，氣固相催化反應(yīng)本是非均相反應(yīng)，宏觀反應(yīng)過程包括反應(yīng)物分子 的外擴(kuò)散、內(nèi)擴(kuò)散、反應(yīng)物分子在催化劑活性表面上的化學(xué)吸附與反應(yīng)，生成物分子的 脫附，以及生成物分子的粒內(nèi)、粒外擴(kuò)散等步驟。對(duì)反應(yīng)過程及反應(yīng)器的特點(diǎn)進(jìn)行分析，可以合理假定而選用簡化模型，即采用一維、 擬均相、平推流模型。氣體以平推流通過催化床，不存在徑向流速分布和軸向流體的返混不存在徑向的溫度差和濃度差，只考慮軸向的溫度差和濃度差按氣流主體中反應(yīng)組分濃度和溫度來計(jì)算反應(yīng)速率而將氣體與顆粒外表面之間的界面濃度差和溫度差、顆粒內(nèi)部的濃度差和溫度差等因素都計(jì)入校正系

12、數(shù)。軸向傳質(zhì)傳熱只是由平推流的總體流動(dòng)引起的。2.2.1物料衡算NA0 xANAi r11F1Xa + dXAF1Na-f卜dNA物料衡算選微元dL高度床層作為物質(zhì)衡算基礎(chǔ)，在穩(wěn)態(tài)時(shí)，根據(jù)質(zhì)量守恒定律在某一時(shí)刻 反應(yīng)物CO在單位時(shí)間內(nèi)進(jìn)入微元體的量等于單位時(shí)間內(nèi)CO流出微元體的量加上在微 元體中消耗的量，即：CO進(jìn)入微元體的量=CO流出微元體的量+ CO在微元體中消耗的物質(zhì)的量用符號(hào)表示為：NA = NA + dN a + rdw公式（1）其中：dw =pBdVB公式（2）丸整理得：-dN =p r （d2）dl公式（3）又連續(xù)流動(dòng)條件下轉(zhuǎn)化率定義為：Nao -NA NaoNa = NAo（1

13、-七）公式（4）對(duì)式（4）微分后代入式（3）得：一兀 一N dx = p r （ d2）dl公式（5）ao a b a 4 t經(jīng)整理后，可以得到如下的微分方程式：dx（- r ） MA = p A dlB G能量衡算同樣取微元體d匕為熱量衡算基礎(chǔ)，設(shè)在微元d匕中的溫度變化為dT。進(jìn)入物料的 總摩爾流量為F0,總體積流量為，換熱介質(zhì)的溫度為TS，則在等壓下對(duì)于放熱反應(yīng)熱量 衡算關(guān)系為：單位時(shí)間d匕中反應(yīng)放出的熱量=單位時(shí)間體系升溫吸收的熱量+單位時(shí)間體系向 外傳出的熱量用公式表示為：r x (-AH ) x dw = VC dT + K (T - T ) Sdl公式(6)AA0 pmqS式中，S

14、為單位長度反應(yīng)管的管壁傳熱面積，S = 3.14ddw=pb胃 )dl帶入反應(yīng)器的物料衡算方程式(5)整理得：wC 空=N (-AH)竺-K (T -T)S公式(7)pm dlA0r dl q s由于選取的反應(yīng)器為單段絕熱式，所以傳向環(huán)境的熱量為零，即：K (T - T) Sdl = 0qs即：NMCdT = NA0dxA(-AHA)公式(8)經(jīng)整理后，可以得到如下的微分方程式：dT = p B( - r) (-AHA )dGcp ,m公式推導(dǎo)x - A0_Z = 1 -公式(9)n A0A0(-r ) = k P P-0.5(1-)公式(10)A公式（11）公式（12）d(1-業(yè)bJ一p p

15、乙M P b k P P -0.5 ( 1 - 一CO H )公式(13)dl Gy t co co2K P PCO, 0P CO H 2Odycodl-4- p Q k P P-0.5(1PCOPH2N b T CO CO 2K P P公式（14）公式（15）公式（16）公式（17）公式（18）公式（19）Nt 0 Ft = GMDdT P b( -rj(-AHA) 頑GcP ,m7冗 7 “N dx p r ( d2)dlA。 a b a 4 tr - NaodyAA 冗d2p y dldT 工(-AH.)(- 5)皂 6件3 dlGCpmA: d 2 P ByA0 dlCpmM* cd4

16、 t3模型參數(shù)選取3.1摩爾定壓熱容由參考文獻(xiàn)7附錄八某些氣體的摩爾定壓熱容與溫度的關(guān)系：Cp m = a + bT + cT2 7 cal - mol-i - K-1公式（20）各種氣體與溫度的關(guān)系如下：c = 26.537 + 7.6831 x 10-3T -1.172 x 10-6 T2c h2O = 29.16 +14.49 x 10-3T - 2.022 x 10-6T2c 2 = 26.75 + 42.258 x 10-3 T -14.25 x 10-6 T2c h2 = 26.88 + 4.347 x10-3T - 0.3265 x 10-6T2c n2 = 27.32 + 6.

17、226 x 10-3 T - 0.9502 x 10-6 T2c pm = E Ypi3.2反應(yīng)熱已知在25C下摩爾反應(yīng)熱為AHr25 = -9838cal/mol,由參考文獻(xiàn)9查知 反應(yīng)熱與溫度的關(guān)系符合下式關(guān)系：T 2T 3AH =-41868 -1.2184 x T + 2.3822 x 10-3 一 -12.1875 x 10-6 一 cal / mol 凹R233.3平衡常素Kp由參考文獻(xiàn)查知平衡常數(shù)隨溫度的變化可根據(jù)下式:1914T-1.782其中 T 溫度，K4計(jì)算方法的確定通過分析，我們需要對(duì)常微分方程組求數(shù)值解，常用方法有歐拉法、改良?xì)W拉法和-9-龍格-庫塔(Runge-Ku

18、tta )法。通過對(duì)以上方法比較可知：若計(jì)算一次f(x,y)的值(歐拉法)，截?cái)嗾`差o(h2)； 若計(jì)算一次f(x,y)的值(改良?xì)W拉法)，截?cái)嗾`差o(h3)；增加計(jì)算f(x,y)的次數(shù)， 可以提高截?cái)嗾`差的階；龍格-庫塔(Runge-Kutta )方法是一種在工程上應(yīng)用廣泛的 高精度單步算法，屬于單步驟法，此算法精度高，截?cái)嗾`差比較小。所謂四階 Runge-Kutta公式，就是采用四次計(jì)算f(x,y)的值的方法，適當(dāng)選擇參數(shù)，使得誤 差盡可能小，相應(yīng)的四階 Runge-Kutta公式的階段誤差為o(h5)。由于本設(shè)計(jì)過程 中，需要解微分方程組，且參數(shù)的推導(dǎo)比較復(fù)雜，我們常采用四階Runge-

19、Kutta法解微 分方程組，即：dyCOdl兀八彳 P。kP P-0.5 ( 1 勺 P2)B T CO CO2K P P1 D 2dT=5=(-豈p。kP pi-%性N B T CO CO2K P PcT 0P CO H OPdlc dlyco,n+1=y + (k + 2k + 2k + k ) co,n 61234_ 1、T = T + (q + 2q + 2q + q )n+1 n 6 123 4k = f (l , y , T ) Al1nco,nnqi = F (l , y , T ) Al = kAl+ 丁，y2 co , nk qA H+ -21-, T + f) A l =

20、k rP=f (l + , y + 才，T + n 2 co, n 2 nV + 匕 T + 婦 A l - k AHr,y i, i)Pk4 = f (l+ A l, y+ k , T + q ) A lq = F (l +A l, y + k , T + q ) A l = k4nco, n3 n 35計(jì)算過程5.1不同進(jìn)口溫度下，床層溫度與床層高度關(guān)系以進(jìn)口溫度341.5C為例:公式（21）也=_ 4 p。k P P -0.5(1-PCOPH 2 )dl N b T CO CO 2K P P式中D = 4mN = 10600kmol / h p = 1500000 g / m3 & =0

21、.63x106 exp(-24900RT)mol / (h.g .atm0.5)gP = 30.5 atmR = 1.987kcal / kmol.h1914T-1.782各組分初始濃度為:組分yCOyHOy COyHy N摩爾濃度0.099020.355520.058820.343520.1432代入公式（10）得:dyco = -3.2 x106 x exp( dl12531.45-)(0.1498 + yco )-。.5 yc。(0.1498 - y )(0.4365 - y )L914c(孑 T78)yco (yco + 0.2685)當(dāng) y = 0.099, Al = 0.075 時(shí)

22、，COk = f (l , y , T )Al = -0.00018 1n co ,n nq = F (l , y , T ) Al = 0.220732 K1n co), n nA lk qk 2 = f (l + , y + 寸,T + 號(hào))A l = -0.00014136由于溫度在微單元內(nèi)變化甚微，所以在這里近似認(rèn)為其為固定值與上面的Z y,c,和AHr均不發(fā)生變化。PAlk qq = F(l + ，y+t,T + 寫)Al = 0.174827 K2 n 2co ,n2 n 2Al k qk3 = f (l + , y + 寸,T + 專)Al = -0.00014138A lk q

23、q3 = F (l + , y + 寸，T + 號(hào))Al = 0.174846 Kk = f (l + A l, y+ k , T + q ) A l = 0.00014132 TOC o 1-5 h z 4nco, n3 n 3q = F (l + A l, y+ k , T + q ) A l = 0.174772 K4nco, n3 n 3即在0.075高度處CO含量和溫度分別為：1y = y + (k + 2k + 2k + k ) = 0.0985734co ,1 co ,061234T = T0 + 6 (q 1 + 2q2 + 2q3 + q4) = 615.20按照此方法依次類

24、推可分別求得在每段的CO出口含量及溫度。5.2不同進(jìn)口流量下，床層溫度與床層高度關(guān)系以進(jìn)口流量Nt0=10600 kmol/h為例，具體計(jì)算過程如下：將公式(21)化簡為：L = -3.7518x106 xexp(-12581.78)(0.1498 + y )-0.5y -E.1498 吃)(0.4365 一吃)dlTco co (914 -1.78)y (y + 0.2685)Tco co當(dāng)yco = 0.099，Al =0.075時(shí)，用用樣的方法可分別求得：k = f (l , y , T )Al = -0.00018 1n co,n nq = F (l , y , T ) Al = 0.

25、220732 K1n co, n nA lk ak2 = f (l + 亍 y + 苛,T + 號(hào))Al = -0.00014136A lk aa = F (l + ，y +T, T + 紅)Al = 0.1748272n 2 co ,n2 n 2Al k ak3 = f (l + , y + 寸,T + 號(hào))Al = -0.00014138A lk aa3 = F (l + , y + 寸，T + 號(hào))Al = 0.174846 Kk = f (l + A l, y+ k , T + a ) A l = 0.00014132 TOC o 1-5 h z 4nco ,n3n3a = F (l

26、+ A l, y+ k , T + a ) a l = 0.174772 K4nco, n3n31y = y + (k + 2k + 2k + k ) = 0.0985734co ,1 co ,061234t = t + 1(a + 2a + 2a + a ) = 615.20 1061234按照此方法依次類推可分別求得在每段的CO出口含量及溫度。6結(jié)果與討論6.1不同進(jìn)口溫度下床層溫度隨床層高度的變化在不同進(jìn)口溫度下，根據(jù)龍格-庫塔(Runge-Kutta )法計(jì)算床層溫度與床層高度關(guān) 系見以下圖表:表一不同進(jìn)口溫度下床層溫度隨床層高度的變化進(jìn)口溫度床層溫(K)度(K)、614.65621.

27、15627.65634.15640.65床層高度(m)0.000615.20621.81628.45635.10641.780.075615.74622.47629.24636.05642.900.150616.29623.13630.02636.98644.000.225616.83623.78630.80637.90645.080.300617.37624.43631.58638.81646.140.375617.91625.07632.34639.71647.180.450618.44625.71633.10640.59648.190.525618.98626.35633.84641.4

28、5649.170.600619.51626.98634.58642.30650.130.675620.04627.60635.31643.13651.050.750620.56628.22636.02643.94651.950.825621.08628.83636.73644.73652.810.900621.60629.44637.42645.51653.650.975622.12630.04638.10646.26654.441.050622.63630.63638.76646.99655.211.125623.13631.21639.42647.69655.941.200623.6363

29、1.78640.06648.38656.641.275624.13632.35640.68649.04657.301.350624.62632.91641.29649.67657.931.425625.11633.45641.88650.29658.531.500625.59633.99642.46650.88659.09床層溫度、(K)、進(jìn)口溫 度(K) 647.15653.65660.15666.65673.15床層高度(m) A0.000648.49655.22661.99668.78675.610.075649.80656.76663.76670.81677.910.150651.09

30、658.24665.46672.73680.040.225652.34659.67667.07674.52681.990.300653.56661.05668.60676.18683.750.375654.73662.37670.04677.71685.330.450655.87663.62671.38679.11686.720.525656.97664.81672.63680.37687.930.600658.02665.93673.78681.49688.980.675659.02666.98674.84682.50689.880.750659.98667.96675.80683.3869

31、0.650.825660.89668.88676.67684.16691.290.900661.75669.73677.45684.84691.830.975662.57670.51678.16685.42692.271.050663.33671.23678.78685.93692.641.125664.05671.88679.34686.36692.951.200664.72672.48679.83686.73693.201.275665.34673.03680.27687.05693.411.350665.92673.52680.65687.32693.571.425666.46673.9

32、7680.99687.54693.711.500666.96674.37681.28687.73693.82溫層床700.00690.00680.00670.00660.00650.00640.00630.00620.00610.00600.00590.00580.00570.00不同進(jìn)口溫度下床層溫度隨床層高度的變化曲線床層高度(m)0000.1500.3000.4500.6000.7500.9001.0501.2001.3501.614.65621.15627.65634.151640.05-647.15653.65660.15666.65673.15圖1不同進(jìn)口溫度下床層溫度隨床層高度的

33、變化由上圖可以看出：同一進(jìn)口溫度下，隨著床層高度的增加，出口溫度逐漸升高；在同一床層高度，隨著進(jìn)口溫度的升高，出口溫度也逐漸升高。6.2床層出口溫度隨入口溫度的變化運(yùn)用龍格-庫塔(Runge-Kutta )法計(jì)算得CO中溫變換過程中，床層出口溫度與 進(jìn)口溫度的關(guān)系如下圖表所示：表2床層出口溫度隨入口溫度的變化入口溫度(K)床層出口溫度614.65625.59621.15633.99627.65642.46634.15650.88640.65659.09647.15666.96653.65674.37660.15681.28666.65687.73673.15693.82床層出口溫度隨入口溫度的

34、變化曲線617.15623.15629.15635.15641.15647.15653.15659.15665.15673.15 入口溫度(K)圖2床層出口溫度隨入口溫度的變化曲線由上圖可以看出：在同一進(jìn)口流量下，床層出口溫度隨著進(jìn)口溫度的升高而 增加，在進(jìn)口溫度不太高時(shí)，出口溫度隨進(jìn)口溫度變化較明顯；當(dāng)進(jìn)口溫度達(dá)到 某一范圍后，出口溫度隨進(jìn)口溫度變化緩慢。6.3不同進(jìn)口流量下床層溫度隨床層高度的變化在不同進(jìn)口流量下，根據(jù)龍格-庫塔(Runge-Kutta )法計(jì)算床層溫度與床層高度關(guān) 系如下圖表所示：表3不同進(jìn)口流量下床層溫度隨床層高度的變化進(jìn)口床層溫度流量(欄(kmolh)98501000

35、0101501030010450床層高度(m) -0.000A615.24615.23615.22615.21615.200.075615.83615.81615.79615.77615.760.150616.41616.38616.36616.33616.310.225616.99616.96616.92616.89616.860.300617.57617.53617.49617.45617.410.375618.15618.10618.05618.00617.950.450618.73618.67618.61618.55618.500.525619.30619.23619.17619.10

36、619.040.600619.87619.79619.72619.65619.580.675620.43620.35620.27620.19620.110.750621.00620.90620.81620.73620.640.825621.55621.45621.36621.26621.170.900622.11622.00621.90621.79621.700.975622.65622.54622.43622.32622.221.050623.20623.08622.96622.84622.731.125623.74623.61623.48623.36623.251.200624.27624

37、.13624.00623.88623.751.275624.80624.66624.52624.39624.261.350625.32625.17625.03624.89624.751.425625.83625.68625.53625.39625.251.500626.34626.18626.03625.88625.73床層溫度(K進(jìn)口流量(kmol/h)1060010750109001105011200床層高度(m)-0.000615.20615.19615.18615.18615.170.075615.74615.73615.71615.70615.680.150616.29616.266

38、16.24616.22616.200.225616.83616.80616.77616.74616.710.300617.37617.33617.30617.26617.230.375617.91617.86617.82617.78617.740.450618.44618.39618.34618.29618.250.525618.98618.92618.86618.81618.750.600619.51619.44619.38619.32619.260.675620.04619.96619.89619.82619.760.750620.56620.48620.40620.33620.260.825621.08621.00620.91620.83620.750.900621.60621.51621.42621.33621.240.975622.12622.02621.92621.82621.731.050622.63622.52622.42622.32622.221.125623.13623.02622.91622.80622.701.200623.63623.51623.40623.29623.181.275624.13624.01623.89623.77623.651.350624.62624.49624.376