反应器设计原理-第四章固定床反应器-.ppt

 第一节 概 述                    当 时， 温度增加， 值减小， 值上升，即： 而 （6） 若xA 沿绝热操作线变化， 以 对 作图，得讲义p114，ch4-11 若方程式（6）得到满足，则ch4-11中x轴上下两块阴影面积应相等。由于 xAi 已知，即面积 一定，通过试差，使面积 等于 ， 即可求出出口转化率： 至于 可由动力学方程式直接对T求导而得，也可用图解法求 l = 0时， l = L时， (4-68) (4-67) 4.5.3 拟均相二维模型 所谓拟均相二维模型就是假定催化剂颗粒与气流之间没有温度和浓度梯度，可作均相处理，但需同时考虑床层的轴向和径向的浓度分布和温度分布。 如果在床层中取一个环状的微元体积，如图4-5所示，对反应组分A的物料衡算，以及热量衡算，可建立二维模型微分方程。假定 等参数均不随位置而变，当作常数。 图4-5 固定床中环状微元体积示意图 对反应组分A的物料衡算： 轴向流动带入量： (4-69) 径向混和扩散带入量： (4-70） 轴向混和扩散带入量： (4-71) 轴向流动带出量： (4-72) 径向混和扩散带出量： (4-73) 轴向混和扩散带出量： (4-74) 单位时间内消耗反应物A的量： (4-75) 定常态条件下，带入总量－带出总量 = 反应消耗量，即： 根据微分方程中值定理，有如下关系式： (4-77) (4-78) (4-79) 上述关系代入式(4-76)，经化简，并略去二次微分项得到 (4-80) 用类似的方法，可导出微元体积内热量衡算式 (4-81) 对于实用的固定床反应器，轴向导热和扩散可以忽略，这样就得 (4-82) (4-83) 还可以利用转化率的关系 又有 (4-84) (4-85) 式中yAO为进料中组分A的摩尔分率，M为混合组分平均分子量。将式(4-84)和(4-85)代入式(4-82)得到 (4-86) 方程组的边界条件为： (4-87) (4-90) (4-89) (4-88) 微分方程组结合边界条件，可用计算机求解得到管式固定床反应器内径向和轴向温度分布和浓度分布。方程组中包含二类参数，其中表征反应特征的参数在反应动力学方程表达式中，传递特征参数 为Der、λer 及边界条件中的ɑw。 式(4-82)和式(4-83)可经坐标变换成无因次形式，令： (4-91) 代入式(4-82)和(4-83)后得到 (4-93) (4-92) 式(4-92)和式(4-93)是二阶非线性偏微分方程组，可以用显式差分法、隐式差分法以及正交配置法数值求解。 由上述分析表明：二维模型合理的考虑了径向传递因素，使模型计算与实际操作更为贴切。对于强放热反应，应该采用二维模型对反应器进行稳态操作分析。而一维模型则主要应用于在线计算和过程控制研究。 4.5.4 非均相模型 对于热效应很大而且反应又非常快的反应过程，需要考虑流体和催化剂颗粒之间的温度梯度和浓度梯度，甚至在颗粒内部也存在不同的温度和浓度。描述这类情况的模型称为非均相模型。非均相模型也分为一维和二维。 一维非均相模型：与平推流拟均相模型的基本假定相同，但考虑了催化剂颗粒外表面和流体主体温度和浓度的差异，相间存在传热和传质。一维非均相模型为： 1. 流体 (4-94) (4-95) 2. 固体 (4-96) (4-97) 边界条件： (4-98) 式中下标S表示固体颗粒，上标S表示颗粒外表面，所以 、 分别为催化剂颗粒外表面的浓度和温度。kG、h为相间传质和传热系数。 二维非均相模型：在考虑流体和催化剂颗粒之间温度和浓度的区别后，得到如下的二维非均相模型： 1. 流体 (4-99) (4-100) 2. 固体 (4-101) (4-102) 边界条件： (4-107) (4-106) (4-105) (4-104) (4-103) 因此只要把催化剂的有效系数和床层的有效导热系数解决好，固定床反应器的设计和