气固相催化反应固定床反应器

1、会计学1气固相催化反应固定床反应器气固相催化反应固定床反应器23颗粒密度床层堆积密度，床层体积颗粒体积床层体积空隙体积PBPBBPB111VV45V31S3S66dVdV球形体积：a21S2SdSdS球形外表面积：6SSVS32SSV6666SVSddddVSS球形比表面积：iidxd1m7：床层空隙率气体密度颗粒当量直径：床层高度空塔气速：平均流速修正的雷诺数，式中：BgsmBggmsmms2mg3BBm:1Re:Re175. 1Re150dddluuddulP 89粒径 ds/mm3.404.606.90质量分率 w0.600.250.1510m1096. 3mm96. 390. 615.

2、 060. 425. 040. 360. 0131Siidxd190644. 01103 . 22 . 61096. 31Re53BgSmGd11Pa10898. 1444. 044. 0146. 21096. 32 . 675. 11903150175. 1Re150175. 1Re15053323BBgS2m3BBSg2mmLdGLdup12131415 操作方式： 绝热、换热两种；操作方式的不同，反应器的结构就不同。 操作方式由反应的热效应和操作范围的宽窄及反应的经济效益等决定。 从反应器的设计、制造及操作考虑，绝热型比较简单。 从设计上讲，基本方程是一样的。1617模型化 对于一个过程

3、，进行合理的简化，利用数学公式进行描述，在一定的输入条件下，预测体系输出的变化。 对同一个体系，根据不同的简化和假定，可以构造不同的模型。 不同的简化和假定，也决定了模型必然含有一些参数，以修正模型与实际体系的差异。 根据不同的简化和假定，分为几种不同层次的模型。18 对于固定床反应器，一般有以下模型： 一维拟均相平推流模型 一维拟均相带有轴向返混的模型 二维拟均相模型 二维非均相模型 二维非均相带有颗粒内梯度的模型 1920一维拟均相平推流模型 质量衡算 在管式反应器中垂直于流动方向取一个微元，以这个微元对A组份做物料衡算：dv输入 输出= 反应 积累FA FA+dFA (-RA)(1-B)

4、Aidl 021 整理得： 对照平推流反应器模型 二者相同出：催化剂堆密度或者以催化剂质量计：以催化剂体积计A0A0A00A00BAAB0A0BAAAA00BAAdddddd1ddxBALRxculculRxcuRlxRcuRlx出A0AA0A0RdxrxcVV22 热量衡算：（仍然是那块体积）输入热量输出热量+反应热效应=与外界的热交换+积累输入：G cp T G质量流量, cp恒压热容输出：G cp(T+dT)反应热效应：(-RA)(1-B)(-H)Aidl热交换：U(T-Tr)didl di反应器直径积累：0U:气流与冷却介质之间的换热系数Tr：环境温度23 将各式代入，得 动量衡算：仍

5、然是Ergun方程gprBA41dducTTdUHRlTi s2mg3BBm175. 1Re150dddulp24 将三个方程联立： 边界条件：L=0, p=p0, xA=xA0, T=T0gprBA41dducTTdUHRlTi0A0BAA1ddcuRlx s2mg3BBm175. 1Re150dddulp25需要注意的问题 1 从解题的角度看，一般壁温恒定，实际情况并非如此。 2 对于低压系统，压降十分重要。 3 U不是物性参数，需实验确定。 4 注意u0, u, um 的关系。 5 如果多根管子并联，体系将自动调节各管的流量，使压降相同，此时各管的处理量不同，转化率不同，造成生产能力和产

6、品质量下降。26典型模拟结果27。为常数，计算大大简化常数中，反应速率的同时：kRcuRlxucTTdUHRlTiA0A0BAAgprBA1dd041dd28 2绝热：若绝热，则T=Tr，或者认为U=0。 此时，将物料衡算式与热量衡算式合并，可得： ：绝热温升，如果在一定范围内视物性为常数， 将不随x及T变化。则： TT0=(xx0)温度与转化率形成一一对应关系， 中，温度可以由T= T0 +(xx0)代替。p0Agp0A0ddmcHFAAuccuHxTiiA00BAA1ddcuRlxAA0AAA0APA0AP0A0TBTAA0A0B0A0BAAdddddd1dd11ddRxFWRxFVRFV

7、cuAlARxlcucuRlx质量为基准，则：如果动力学方程以颗粒准，则：学方程以颗粒体积为基力中的基准对应，如果动要求与反应动力学方程积：左端上下同乘床层截面变形：2930可逆放热反应绝热反应器的最优化（以SO21/2O2=SO3为例）xT平衡线等速率线0二氧化硫氧化反应Tx图示意31 二氧化硫氧化反应气固相催化反应，用于硫酸生产，可逆，强放热，绝大多数生产过程采用多段绝热操作。 最优化目的：在完成一定生产任务的条件下，使用的催化剂最少。 已知条件：第一段入口和最后一段出口转化率；第一段入口反应物浓度，各物性参数；段与段间采用间接冷却。 可以改变的参数：各段的入口温度；段与段之间的转化率。3

8、2以四段为例： 催化剂用量为：（基于拟均相平推流模型） 基于某一动力学方程，适当选取各段的入口温度；段与段之间的转化率共7个（N段为2N1个）参数，使W最小。4outin4out3in3out2in21outin1,d,d,d,d43210AxxxxxxxxTxrxTxrxTxrxTxrxFW33x1in,T1inx1out,T2inx2outT3inx3outT4inx4out第一段第二段第三段第四段34斜线为段内操作线，斜率为1/。水平线表示段间为间接冷却，只是温度降低，转化率不变。xT0二氧化硫氧化反应Tx图示意在Tx图上看：35 调用最优化程序，就可以求得W最小值？ 可以，但很困难。

9、进一步数学处理： 在任意一段内，当xin及xout确定之后，应选取适当的进口温度Tin，使催化剂量最小。0d,1,din2in0AoutinoutinouinxTTxrrTxrdxTTxrxFWxxxxxxt36 在任意相邻两段间：速率。等于前一段的出口反应后一段的入口反应速率少，无论中间转化率是多即：变上限定积分的偏微分。使两段催化剂之和最小选取中间转化率TxrTxrTxrTxrTxrxTxrxxxxTxrxTxrxFWFWiiiiiiiixxixxiiiixxixxiiiiiiiiiii,0,1,10,d,d,d,din11outin11out1outin1out1A01A0out1in1

10、outinout1in1outin37 汇总：01,01,01,01442444433333233332222222222111112144332211dxTrrTxrTxrdxTrrTxrTxrdxTrrTxrTxrdxTrrinxxininoutoutinxxininoutoutinxxininoutoutinxxoutinoutinoutinoutin第四段：第三、四段之间：第三段：第二、三段之间：第二段：第一、二段之间：第一段：38 七个方程，七个未知数，可能是唯一解。 讨论：从Tx图上看：xT0二氧化硫氧化反应Tx图示意39404142 SPBAO3HOCHC3OCHHC232462

11、2346TkpkpLR13636837.19exphkgkmol11BARA4344112tWSWKhmkJ88dhhUS45催化剂有效因子 总括给热系数 U 催化剂比活性 LR 管间冷载体熔盐温度 TS 进料组成（原料气中 邻二甲苯的初摩尔分数）yA0 操作压力 p 管内径 dt 管数 n 混合气入口温度 T0 原料气体混合物单管入 口质量流速 G 460A0mBAA1ddcuRlxWt0BAPBm411ddTTdhRHculT4748 一、带有轴向返混的一维模型 非理想模型，当平推流模型描述不够满意时采用。 修正轴向热量、质量返混带来的与平推流模型的偏离。 物理模型： 在拟均相平推流模型上

12、迭加一个轴向返混，与 非理想流动中介绍的返混模型相同，但增加热扩散的考虑。49 稳态，在dVR体积中对A组份做物料衡算： 输入 输出 反应 输入输出反应LdlcA0,FA0,xA0=0,V0cA,FA,xA,VFA, xAFA+dFA,xA+dxAdVRlARcuAFAllcclEFFAlcEFiiiid1dd ddddddddBAAAAAZAAAZA50 将以上三式合并，得： 式中，EZ为轴向有效扩散系数。 相应，在同样条件下，对dVR做热量衡算：01ddddAA2A2ZBRlculcElTATTcAullTTlATcAuiiiiddddddddZpgZpg输出：输入：51 反应： 散热：

13、输入放热输出散热 整理得： Z为轴向有效导热系数lTThdlAHRiidd1r0BA041ddddrBAp22ZTTdhHRlTculTiog52 边值条件： 二阶常微分方程组，两点边值问题。 可调用程序求解0dd0dddddd0AZ0pg0AZAA00lTlcLllTTTculcEccul53讨论： 1 轴向扩散的引入，可以导致温度、浓度分布趋于平缓。 2 许多不确定因素可以归结到轴向扩散中。 3 轴向扩散可能会造成多重态。 4 轴向扩散系数与轴向导热系数有一定的函数关系。 5 经验证明，当床层厚度大于50倍颗粒直径时，轴向热质扩散（轴向返混）对出口转化率所造成的影响可以忽略不计。54 6

14、轴向扩散系数和轴向导热系数都不是物性参数。其中都包含了流体和固体颗粒双重的贡献。 7 轴向扩散系数和轴向导热系数需通过实验求取或参考文献值及通过经验公式求取。55二、二维拟均相模型 二维：轴向和径向 对于径向存在较大的温度差、浓度差的反应器，一维模型有时不能满足要求，需要考虑径向的温度浓度分布。 与一维模型相比，考虑的因素更多，得到的结果更复杂，各有优缺点。56模型假定： 1 反应在圆管式反应器中进行。 2 流体在催化剂管内为非理想流动，存在着轴、径向的质量和热量扩散。 3 流固相之间没有温度、浓度差。 4 扩散遵循Fick扩散定律。57 在管式反应器中取一微元：drdlRr58 定常态条件下

15、就环形微元对A做物料衡算：0dd21d2dd2dd2ddd2d2d2AArAAlAAAArAlA积累：反应：输出：输入：lrrRrcElrllcclErrllccurrrrccrElrrlcErrucrrB59 输入输出=反应 整理得：BAAA2A2rBAA2A2lA2A2r1111RlcurcrrcERlculcErcrrcE散项被忽略：在多数情况下，轴向扩60 热量衡算：lrrHRrTlrlTrrllTTcurrrrTTrlrrllTTlrrTcurrdd21d2d2dd2ddd2dd2d2BArlpgrlpg微元内反应热：输出微元的热量：输入微元的热量：61 输入输出反应 与质量衡算类似

16、，轴向热扩散项可以忽略： 动量衡算方程与一维模型相同。HRlTculTrTrrTBApg22l22r11HRlTcurTrrTBApg22r11 s2mg3BBm175. 1Re150dddulp62 边界条件：wwrAA00A0A, 000, 000,00TThrTrcLlRrrTrcLlrppTTccRrll=0l=L63 在任意截面上流体的平均温度浓度RRRRrrcRRrcrcrrTRcuRrTcruT0A22A0A02pg2pg0d2d2d2d264关于模型参数 模型参数是模型的一个重要组成部分，与模型紧密结合。 模型参数包含轴径向有效导热系数与扩散系数及流体与管壁之间的给热系数。 模

17、型参数的取得，与实验条件有关，在具体应用时，要选择尽可能接近应用条件的文献值。656667热量传递热量传递拟均相一维平推流模型热量传递热量传递带有轴向返混的拟均相一维模型热质传递热量传递热量传递拟均相二维模型热质传递68固体相热量传递热量传递二维非均相模型热质传递抽象成为热量传递二维非均相模型热质传递热量传递流体相69热量传递考虑颗粒内梯度的二维非均相模型热质传递热量传递流体相固体相707172模型化 对于一个过程，进行合理的简化，利用数学公式进行描述，在一定的输入条件下，预测体系输出的变化。 对同一个体系，根据不同的简化和假定，可以构造不同的模型。 不同的简化和假定，也决定了模型必然含有一些参数，以修正模型与实际体系的差异。 根据不同的简化和假定，分为几种不同层次的模型。73 整理得： 对照平推流反应器模型 二者相同出：催化剂堆密度或者以催化剂质量计：以催化剂体积计A0A0A00A00BAAB0A0BAA