7固定床反应器69.pptVIP

5 压力沿床层变化的等温或变温固定床反应器计算过程 设计问题 V=f(X) 计算问题 X=f ’(V) 对A作摩尔衡算 PFR CSTR (代数式) BATCH 求A的反应速率 求A的化学计量关系 代入上二式才能求得V与X关系 * 对液相反应可认为: V=V0或 即得 对气相反应要考虑: 1.等温等压下因反应引起的膨胀率 2.因温度,压力变化引起的气体密度变化 或是 计算中温度变化要求求出T(X),(一般是非等温非绝热关系) 压力变化要求求出P(X) , 前者要用热衡算,后者用Ergun等公式二者均需用计算机算 * 管式反应器计算时,为便于列物料衡算和热量衡算,常用摩尔流量Fi作变量.而反应速率方程必需用i组分浓度Ci耒表示. 为解联立方程组,需用Fi[mol/s]置换掉速率方程中的Ci[mol/m3].以便求解Fi沿床层的变化. 即 由气体方程得 即 式中P,T沿层床长L是变化的,即它们与L的关系需由方程式给出 P/Po 称相对压力 T/To称相对温度 * 上机题4 平均分子量为100g/mol的烷烃，在一反应器内进行脱氢反应，以生产辛烷值较高的烯烃，反应式为：烷——烯＋H2是一级气相催化反应。床层空隙率为ε=0.4，催化剂密度ρc=2.6 [克/cm3]。反应在等温下进行，反应速率 (-rA)W=kCA，k=0.02 [cm3/s/克催化剂]。原料气粘度μ=1.5·10-6 [Kg/s/dm]，并假设反应过程中粘度不变化，反应器进口压力P0=2000 [kPa]，此时气体密度ρ0=0.032 [克/cm3]，进料摩尔流量FA0=440 [mol/s]。反应器为直径2.4米，长度为25米的管式反应器时，内装催化剂173870 [Kg]，床层压降符合Ergun公式： 式中G为单位截面上的质量流率 [Kg/s/m2]，Dp为催化剂颗粒直径0.002米，z为管长，ρ为气体密度 [克/cm3]。 * 试求反应器内每一点处的无因子压力值y=P/P0和反应转化率x，并以管内装催化剂量W做自变量，画出y(W)和x(W)图。 为减少压力损失，并提高转化率，将上述管式反应器换成如下的球形反应器，内装催化剂与管式反应器相等。同样画出y(W)和x(W)图并说明与管式反应器的差别。球半径R=3米，半轴长L=2.7米。 R L L’ Z0=0 Z0=L Z0=L+L’ (y为等温下压力损耗系数) (x为转化率, 由气态方程当n变化时有 ) 解： * 对管式反应器而言： 对球形反应器而言： 由上述各式，代入数值并用matlab求解如下： 对球形反应器： 对管式反应器： * 近似计算时，在T0－－T范围内，可取Cpt，平均热容。 积分显热为 若流量单位用G[Kg/hr]，则总显热为 应是T温度下的反应热 为生成热，为负值。 等压反应热 反应前 反应后 反应热为负值时，表示 放热反应。 沿轴向G不变(G不受反应影响) 但C’p因组分变化而变化，并因温度变化而变化。 * 反应物 产物 反应物 产物 反应物 产物 反应物 产物 热衡方程 面积 总扩给热系数[千焦/米2.0C.时] 等温反应 绝热反应 非绝热，非等温反应 (U=const) 再依据Hess定律 原料在降温过程中的焓变。 原料在升温过程中的焓变。 T * 2. 等温反应 物衡 热衡 反应速率式 图解法： * 3.绝热反应 绝热温升： 可以数值解（T，C分布）。 xA 0 0.1 0.2 ………… T T0 －rA L * x A T 1 T 0 ( - r A ) = 0 ( - r A ) + d - d 0 , 0 t t F F = = d 例题：p196，p213－6 求一段，二段催化剂装填量。 * 数值解法 xA 0 0.1 0.2 0.3 0.4 0.5 T ℃ 370 -rA * 1/(rA) * W kg 0 * 4.非绝热非等温反应 数值解法： （递推公式） （1） （2） 由2－1－2－1迭代至稳定值。 或者： 给出T0，X0。由龙格库塔法求解。 U=const * 0，2，3，4，5……………….N 热斑hot spot Tm=const Tm=const * 上机题3 固定床反应器的设计计算： 异丁烯催化叠合成辛烯的气相反应，在非等温非绝热有冷却夹套的管式固定床内进行。反应器直径D＝0.05ｍ，长L＝2ｍ，冷却介质温度Tｍ＝100℃，总传热系数U＝3600[cal/hr?m2