催化反应工程华东理工大学多段换热式催化反应器

1、1会计学催化反应工程华东理工大学多段换热式催化反应工程华东理工大学多段换热式催化反应器催化反应器 平衡线，平衡线，Tm线不变线不变 冷却过程中，冷却过程中，x不变不变 各段绝热线斜率不变各段绝热线斜率不变 无返混，无返混，VR最小最小 付线调节温度付线调节温度 平衡线，平衡线，Tm线不变线不变 冷却过程中，冷却过程中，x变化变化 各段绝热线斜率不变各段绝热线斜率不变 有返混，有返混，VR大大 冷激线调节温度冷激线调节温度 逐段逐段Te， Tm线上移线上移 冷却过程中，冷却过程中，x不变不变 各段绝热线斜率上移各段绝热线斜率上移 一般放在最后两段之间一般放在最后两段之间 阀门调节温度阀门调节温度

2、 121进出NNNC 22211进进出出xNxN 32211进进出出TNTNTNCC 得式代入将出进,312NNNC 42211CCTTNTTN进进出出联立式联立式(2)、(4)求解，可得求解，可得CCCCTTTTxxTTxTTx进出进出进进出出21212211niiSiniiSniiRRVVVV10001, 01,1.解析法求解解析法求解n段，已知段，已知V0, x1进进，xn出出求求VR=最小时的各段最小时的各段x，T,00SRVVX1，进，进 ，XN，出出 已知已知 2个个X i+1，进，进=X i，出，出 n-1个方程个方程T i，出，出-T i，进，进=(X i,出出-X i，进，进

3、) n个方程个方程个变量出进出进nTTxxfViiiiR4,0 0标准接触时间，标准接触时间，hVS 0标准状态下初始态体积流量，标准状态下初始态体积流量，m3(stp)/h催化反应工程催化反应工程进进出出，选, 1, 1, 1nnTTxx进出,iiRTxfV 0,ijxTiRjixV出进进0,ijTxiRjiTV进出进独立变量数独立变量数2n-1反应器体积最小的必要条件反应器体积最小的必要条件为独立变量为独立变量催化反应工程催化反应工程第一类条件式确定段间间接换热过程第一类条件式确定段间间接换热过程第二类条件式确定各段出口转化率第二类条件式确定各段出口转化率催化反应工程催化反应工程出出,1,

4、 0, 0iiiixx第二类条件第二类条件式式0,1, 0, 0出出即iiiixx第一类条件式第一类条件式0,0进iT个进出进出nxxTTTiiiii,max,0,出iRxV0,1,出进进出出出iiiiRiiRiiRiRxTTVxVxVxV2. 出口温度受耐热温度限制出口温度受耐热温度限制独立变量数独立变量数n-1，取各段出口转化率取各段出口转化率xi，出出为独立变量为独立变量iiixT出进,0,1,进出出出iiRiiiRiiRiRTVxVxVxV0, 0,1, 0, 0, 0进出出出即iiiiiiiiiTxxx催化反应工程催化反应工程3. 直接搜索法直接搜索法最优化问题：最优化问题：1)反应

5、器催化床体积最小反应器催化床体积最小 2)产量最大产量最大 3)能耗最小能耗最小 4)经济效益最大经济效益最大解析法存在的问题：解析法存在的问题：1)VR=min，已推导；但产量已推导；但产量=max?经济效经济效 益益=max?2)有温度限制时，解析法求解很复杂。有温度限制时，解析法求解很复杂。3)直接换热式反应器，如何优化？直接换热式反应器，如何优化？4)不能判断多维极值的不能判断多维极值的“尖尖”与与“平平”催化反应工程催化反应工程直接搜索法直接搜索法1）分析独立变量数，确定独立变量）分析独立变量数，确定独立变量设计型：设计型：VR=f(x i,出出，T i,进进)=Min操作型：操作型

6、：Prod=f(x i,出出，T i,进进)=Max 必然存在一组必然存在一组(x i,出出，T i,进进)使使VR=Min，或，或Prod=Max3)给定初值给定初值-得得VR(或或Prod)-修正修正-直至直至VR=Min，或，或Prod=Max4)如何修正，最优化方法，如如何修正，最优化方法，如powell法法2)建立数学模型建立数学模型出出,1,0,0iiiixx5)若在改变初值的过程中发现下述现象，若在改变初值的过程中发现下述现象，改变搜索方向改变搜索方向T i,进进催化剂耐热温度催化剂耐热温度注意：直接搜索法求得的结果，判别注意：直接搜索法求得的结果，判别iiiiiRTTxxfV,

7、出进出进个，进出进出nxxTTiiiii,3 多段直接换热式反应器各段始末温多段直接换热式反应器各段始末温 度及转化率的度及转化率的 优化设计优化设计注意：原料气冷激，注意：原料气冷激，i i相同相同 非原料气冷激，非原料气冷激，i i不相同不相同催化反应工程催化反应工程共有共有5n-1个变量个变量进进出出, 1, 1,ipicpciipiTCNTCNTCN独立变量数独立变量数2n-1,当第一段进口温度给定，则独立变量数维当第一段进口温度给定，则独立变量数维2n-2进进出出, 1, 1,iiiixNxN冷激过程热量衡算冷激过程热量衡算 n-1 非原料气冷激非原料气冷激冷激过程中物料衡算冷激过程

8、中物料衡算 n-1n-1x i,进进，x n,出出已知已知 2 催化反应工程催化反应工程4 多段球形催化床甲醇合成多段球形催化床甲醇合成 反应器非均相数学模型及反应器非均相数学模型及 优化设计优化设计1 物料衡算物料衡算2 催化床非均相数学模型催化床非均相数学模型3 甲醇合成反应动力学甲醇合成反应动力学4 气固相间传质系数和给热系数气固相间传质系数和给热系数5 催化剂内扩散效率因子催化剂内扩散效率因子6 优化设计优化设计7 计算框图计算框图8 计算结果计算结果OHCHHCO322OHOHCHHCO23223独立反应数独立反应数2，关键组分数，关键组分数2，yCO，yCO2TinTinCOinC

9、OCOCOTinTNNCyyyyNN,22122122令22223232,COCOinCOinCOinHHyyCyyyy22,COCOinCOinCOinmmyyCyyyy2222,COinCOinOHOHyCyyyCyyinNN,22CyyinCHCH,44催化反应工程催化反应工程)(4,2,SCOgCOeCOggCOccSdRkRdN)(4,2,2222SCOgCOeCOggCOccSdRkRdN)(,SCOgCOeCOgsCOsCObccSkrCOR)(,22222SCOgCOeCOgsCOsCObccSkrCOR)(42gsesgpTTTSdRRdTCN)()()(222,gsesCO

10、sCOsCOCOsCOsCObTTSHrHrCOR物料衡算 气相固相热量衡算 气相固相COCOCOinCOinCOTCOTCOyyyyyNyNN221221221,222122221221,COCOCOinCOinCOTCOTCOyyyyyNyNN22221 COCOCOCOCOdyydyyDDdNCOCOCOCOCOdyydyyDDdN222221 2,221221221COCOinCOinCOTyyyyNDD式中式中催化反应工程催化反应工程gCOgCOsCOsCOgCOsCOsCOgCObyyyryryCORRdRdgCO,2222,2212)21 (4gCOgCOsCOsCOgCOsCO

11、sCOgCObyyyryryCORRdRdgCO,222222,22212)21 (4)()(4222,2COsCOsCOCOsCOsCOpTbgHrHrCNCORRdRdTssCOggCOgeCOgsCOsCObTyTyZRPSkrCOR,ssCOggCOgeCOgsCOsCObTyTyZRPSkrCOR,22222)()()(222,gsesCOsCOsCOCOsCOsCObTTSHrHrCORZ为混合为混合气体的压气体的压缩因子。缩因子。)()(4222,2COsCOsCOCOsCOsCOpTbgHrHrCNCORRdRdTssCOggCOgeCOgsCOsCObTyTyZRPSkrCO

12、R,ssCOggCOgeCOgsCOsCObTyTyZRPSkrCOR,22222)()()(222,gsesCOsCOsCOCOsCOsCObTTSHrHrCOR3121222221)1 (HHCOCOCOCOHCOCOCOfKfKfKffkdWdNr4221222222221)1 (HHCOCOCOCOHCOCOCOfKfKfKffkdWdNr386. 082. 032365. 0765. 0pepeBfgfgiDRRDGkJ35. 0323023. 0876. 2pepemfpmgpmmsHRRMCGCMJ上式中的Rep适用范围为0.015000。ERRVFEee)(3432311式中，

13、R1及R2分别是催化床外半径及内半径。、gCOy,、gCOy,2gTyes反应速率反应速率反应热反应热恒压热容恒压热容微分方程组微分方程组输入输入给定初值给定初值龙格龙格库塔法求库塔法求解方程解方程 求求产量产量=给定值给定值?VR最小最小输出输出调整调整R2调整调整Tb2,inTb3,inyesnono内扩散内扩散效率因子效率因子传质系数传质系数给热系数给热系数R2外半径外半径、sCOy,、sCOy,2sT催化反应工程催化反应工程Vin 1.8105 m3 (STP)/h ,yin CO 0.1053 CO2 0.0316 H2 0.7640 CH3OH 0.0055 H2O 0.0002

14、N2 0.0499 CH4 0.0435P 5 Mpa 生产能力为生产能力为 100 kt/year生产条件：生产条件：反应器反应器R/myi,gtg/yi,sts/ COCO2CH3OHCOCO2CH3OH11.3500.10530.03160.0055210.00.10520.03160.0056210.41.2620.10320.03120.0089217.10.10320.03120.0090217.51.1600.10090.03080.0128225.30.10080.03070.0130225.71.0370.09820.03020.0174234.70.09810.03020.

15、0175235.10.8730.09520.02960.0225245.00.09510.02950.0226245.40.6000.09200.02890.0278255.90.09190.02890.0278256.1反应器t/yiCOCO2CH3OH1进口出口215.0256.20.10530.09340.03160.02910.00550.02552进口出口218.0257.20.09340.08150.02910.02730.02550.04463进口出口223.0254.20.08150.07130.02730.02660.04660.0598个进出进出nxxTTTiiiii,ma

16、x,0,出iRxV0,1,出进进出出出iiiiRiiRiiRiRxTTVxVxVxV2. 出口温度受耐热温度限制出口温度受耐热温度限制独立变量数独立变量数n-1，取各段出口转化率取各段出口转化率xi，出出为独立变量为独立变量3. 直接搜索法直接搜索法最优化问题：最优化问题：1)反应器催化床体积最小反应器催化床体积最小 2)产量最大产量最大 3)能耗最小能耗最小 4)经济效益最大经济效益最大解析法存在的问题：解析法存在的问题：1)VR=min，已推导；但产量已推导；但产量=max?经济效经济效 益益=max?2)有温度限制时，解析法求解很复杂。有温度限制时，解析法求解很复杂。3)直接换热式反应器，如何优化？直接换热式反应器，如何优化？4)不能判断多维极值的不能判断多维极值的“尖尖”与与“平平”催化反应工程催化反应工程iiiiiRTTxxfV,出进出进个，进出进出nxxTTiiiii,3 多段直接换热式反应器各段始末温多段直接换热式反应器各段始末温 度及转化率的度及转化率的 优化设计优化设计注意：原料气冷激，注意：原料气冷激，i i相同相同 非原料气冷激，非原料气冷激，i i不相同不相同催化反应工程催化反应工程共有共有5n-1个变量个变量TinTinCOinCOCOCOTinTNNCyyyyNN,22122122令22223232,COCOinCOinCOinHHyyC