反应器中的混合对反应的影响

1、反应器中的混合对反应的影响第1页，共96页，2022年，5月20日，3点55分，星期二 概述返混定义理想置换反应器（PFR）理想混合反应器（CSTR）两种流体流动极限模型反应器物料在反应器中的混合在空间上的混合在时间上的混合这种混合过程返混简单混合返混第2页，共96页，2022年，5月20日，3点55分，星期二理想反应器的流动模式 - 平推流 和 全混流。 间歇釜全混釜u = const平推流 理想的平推流和间歇釜停留时间均一，无返混。 全混釜反应器的返混无穷大，出口物料停留时间分布与釜内物料的停留时间相同。 第3页，共96页，2022年，5月20日，3点55分，星期二第4页，共96页，202

2、2年，5月20日，3点55分，星期二实际反应器流动形式的复杂性u沟流回流存在速度分布 存在死区和短路现象 存在沟流和回流 偏离理想流动模式，反应结果与理想反应器的计算值具有较大的差异。 DeadzoneShortcircuiting第5页，共96页，2022年，5月20日，3点55分，星期二 实际流动与PF的偏离(轴向返混) 1. 并行流存在速度分布2. 存在分子热运动DistributionH2空气TT0热反馈QEddy涡流DispersionVLaminar第6页，共96页，2022年，5月20日，3点55分，星期二 实际流动与CSTR的偏离短路活塞流滞流离子反应 非CSTR胃曝气池 生化

3、反应CSTR第7页，共96页，2022年，5月20日，3点55分，星期二 影响反应结果的三大要素： 停留时间分布（residence time distribution, RTD）凝集态（state of aggregation）早混或迟混（earliness and lateness of mixing）RTD对反应的影响 注入检测第8页，共96页，2022年，5月20日，3点55分，星期二平均停留时间为 实际停留时间ti不尽相同，转化率x1, x2, , x5亦不相同。出口转化率应为5个质点转化率的平均值，即 InjectionDetection第9页，共96页，2022年，5月20日，3

4、点55分，星期二聚集态的影响理想反应器假定混合为分子尺度，实际工程难以达到，如 结团弥散两种体系的反应程度显然应该是不同的。 工程中，尽量改善体系的分散尺度，以达到最有效的混合，从而改善反应效果。 鼓泡气体液体喷雾第10页，共96页，2022年，5月20日，3点55分，星期二混合迟早度的影响早混晚混即使两反应体系的空时相同，由于反应混合的迟早不同，反应结果也不相同。 第11页，共96页，2022年，5月20日，3点55分，星期二第一节 连续反应器中物料混合状态分析一、混合现象的分类按混合对象的年龄可以把混合分成两种：（1）同龄混合：相同年龄物料之间的混合（2）返混（back mixing） ：

5、不同年龄物料之间的混合混合尺度的大小混合也可分为两种类型：宏观(macro-)混合：设备尺度上的混合微观(micro-)混合：物料微团尺度上的混合第12页，共96页，2022年，5月20日，3点55分，星期二混合的机理总体流动：搅拌器旋转时使釜内液体产生一定途径的循环流动设备尺度上的宏观均匀高速旋转的旋涡与液体微团产生相对运动和剪切力更小尺度上的均匀分子扩散微团最终消失微观均匀均相第13页，共96页，2022年，5月20日，3点55分，星期二提高混合效果的措施：消除打旋现象：螺旋浆式涡轮式第14页，共96页，2022年，5月20日，3点55分，星期二提高混合效果的措施：加设导流筒螺旋式 涡轮式

6、第15页，共96页，2022年，5月20日，3点55分，星期二搅拌器的型式高转速搅拌器1、螺旋桨式搅拌器 螺旋桨式搅拌器的总体循环流动 推进式 三叶片式第16页，共96页，2022年，5月20日，3点55分，星期二2、涡轮式搅拌器涡轮式搅拌器的总体循环流动 a-直叶圆盘涡轮 b-弯叶圆盘涡轮 c-直叶涡轮 d-折叶涡轮 e-弯叶涡轮第17页，共96页，2022年，5月20日，3点55分，星期二大叶片低转速搅拌器1、桨式搅拌器2、框式和锚式搅拌器3、螺带式搅拌器 a-锚式 bc-框式 d-螺带式第18页，共96页，2022年，5月20日，3点55分，星期二二、连续反应过程的考察方法 在连续搅拌反

7、应釜或管式反应器中进行反应，如果反应物料的微观混合程度不同，则考察方法即研究方法就不同。微观混合有两种极限状态，完全混合和完全不混合，研究方法完全不同。第19页，共96页，2022年，5月20日，3点55分，星期二1、以反应容积（或微元）为对象当物料微观混合为完全混合时，物料呈分子状均匀分散，物料不存在微团。对于搅拌反应器，物料以反应器为边界，对于管式反应器，物料以dVR为边界，所以研究的对象分别为反应器容积VR和反应器微元容积dVR。第20页，共96页，2022年，5月20日，3点55分，星期二2、以反应物料为对象当物料微观混合为完全不混合时，物料呈微团独立运动，物料的边界为微团的边界，所以

8、以微团为研究对象，结合物料的停留时间分布函数和动力学方程方程可以有定量结果。第21页，共96页，2022年，5月20日，3点55分，星期二 如果微观混合处于中间状态，则几个微团可以组成微元，研究对象为微元，目前只有定性的认识，没有定量结果。综上所述，考察对象都是物料，不同的是按照微观混合的程度划分考察的范围：完全混合反应容积VR或dVR；中间状态微元（由微团组成）； 完全不混合微团第22页，共96页，2022年，5月20日，3点55分，星期二第二节 停留时间分布停留时间分布的数学描述停留时间分布的实验测定几种流型的停留时间分布函数与分布密度停留时间分布的应用第23页，共96页，2022年，5月

9、20日，3点55分，星期二 一、流体在反应器内的停留时间分布返混程度停留时间分布规律定量确定停留时间概率函数概率密度函数数学期望方差两个特征值两个函数描述第24页，共96页，2022年，5月20日，3点55分，星期二1、停留时间分布的数学描述全混流反应器：机械混合最大 逆向混合最大 返混程度无穷大平推流反应器：机械混合为零 逆向混合为零 返混程度等于零间歇反应器：机械混全最大 逆向混合为零 返混程度等于零 反应器内的返混程度不同停留时间不同浓度分布不同反应速率不同反应结果不同生产能力不同非理想流动反应器：介于两种理想情况之间停留时间是随机变量，因此停留时间分布是一种概率分布。第25页，共96页

10、，2022年，5月20日，3点55分，星期二 1）、停留时间分布的定量描述 借用人口统计学（Population）两个统计参数 a) 社会人口的年龄分布和 b) 死亡年龄分布，在反应工程中假设： 在反应器出口流体微元的停留时间分布：E(t), F(t)反应器注入检测原料出料a)b) 各微元保持 独立身份(identification), 即微元间不能混合c) 不研究微元在反应器内的历程, 只研究它在反应器内的停留时间, 即寿命。则定义：第26页，共96页，2022年，5月20日，3点55分，星期二停留时间（寿命）的概念？ 例：在连续操作的反应器内，如果在某一瞬间（t=0）极快地向入口物流中加入

11、100个红色粒子，同时在系统的出口处记下不同时间间隔流出的红色粒子数，结果如下表。如果假定红色粒子和主流体之间除了颜色的差别以外，其余所有性质都完全相同，那么就可以认为这100个粒子的停留时间分布就是主流体的停留时间分布。停留时间范围tt+t0-22-33-44-55-66-77-88-99-1010-1111-1212-14出口流中的红色粒子数02612182217126410分率N/N00.020.060.120.180.220.170.120.060.040.010第27页，共96页，2022年，5月20日，3点55分，星期二以时间t为横坐标，出口流中红色粒子数为纵坐标，将上表作图：若以

12、停留时间t为横坐标， 为纵坐标作图，则每一个长方形的面积为即表示停留时间为tt+t的物料占总进料的分率。第28页，共96页，2022年，5月20日，3点55分，星期二假如示踪剂改用红色流体，连续检测出口中红色流体的浓度，如果将观测的时间间隔缩到非常小，得到的将是一条连续的停留时间分布曲线。E(t)t t+dt t图中曲线下微小面积E(t)dt表示停留时间在t和t+dt之间的物料占t=0时进料的分率。第29页，共96页，2022年，5月20日，3点55分，星期二(1) 分布密度函数 E(t) E(t)定义为在t=0时刻进入反应器的流体微元，在t时离开反应器的概率根据定义，E(t)应具有归一性，即

13、因为当时间无限长时，t = 0时刻加入的流体质点都会流出反应器第30页，共96页，2022年，5月20日，3点55分，星期二(2) 分布函数F(t) F(t)函数定义为在t时刻离开反应器流体质点数占总示踪流体质点数的分率 E(t)与F(t)的关系： 第31页，共96页，2022年，5月20日，3点55分，星期二 二、停留时间分布规律的实验测定入口处加“激励” 出口处研究“响应”示踪法停留时间分布规律确定研究阶跃示踪法脉冲示踪法从原物料到示踪剂的阶跃变化用示踪剂代替原物料的瞬间变化第32页，共96页，2022年，5月20日，3点55分，星期二停留时间分布的测定一般采用示踪技术，示踪剂选用易检测其

14、浓度的物质，根据其光学、电学、化学及放射等特性，采用比色、电导、放射检测等测定浓度。选择示踪剂要求： 1) 与主流体物性相近；2) 高低浓度均易检测，以减少示踪剂的用量；3) 不产生相变或相转移；4) 示踪剂浓度易转变为光信号或电信号，以便于计算机数据采集和处理。停留时间的测定方法根据示踪剂的加入方式分为脉冲法、阶跃法和周期输入法，前两者应用较广。第33页，共96页，2022年，5月20日，3点55分，星期二1、阶跃示踪法从某一时刻起，将原物料全部切换为示踪物，示踪物浓度阶跃突变。01.0(a)“激励”曲线01.0(b)“响应”曲线图4-3 阶跃法测定停留时间分布函数第34页，共96页，202

15、2年，5月20日，3点55分，星期二由图可知，在t=0时，C=0；t ， C C0 时间为t时，出口物料中示踪剂浓度为C(t)，物料流量为V，所以示踪剂流出量为V C(t)，又因为在时间为t时流出的示踪剂，也就是反应器中停留时间小于t的示踪剂，按定义，物料中停留时间小于t的粒子所占的分率为F(t)，因此，当示踪剂入口流量为VC0时，出口流量VC0 F(t)，所以有：因此，用此法可直接方便地测定实际反应器的留时间分布函数。第35页，共96页，2022年，5月20日，3点55分，星期二2、脉冲示踪法在尽可能短的时间内，用示踪物瞬间代替不含示踪物的物料，然后又立刻恢复为原来的物料。即给物料一个脉冲讯

16、号，与之同时测定出口流的响应曲线图4-4 脉冲法测定停留时间分布密度函数01.0(b)“响应”曲线0(a)“激励”曲线第36页，共96页，2022年，5月20日，3点55分，星期二设t0时间内注入示踪剂的总量为M（mol），出口处浓度随时间变化为C(t)，在示踪剂注入后t t+dt时间间隔内，出口处流出的示踪剂量占总示踪剂量的分率：若在注入示踪剂的同时，流入反应器的物料量为N，在注入示踪剂后的t t+dt时间间隔内，流出物料量为dN，则在此时间间隔内，流出的物料占进料的分率为：第37页，共96页，2022年，5月20日，3点55分，星期二示踪剂的停留时间分布就是物料质点的停留时间分布，即：因此

17、：有：只要测得V，M和C（t），即可得物料质点的分布密度。第38页，共96页，2022年，5月20日，3点55分，星期二由于M=VC0 t0， C0 及t0难以准确测量，故示踪剂的总量可用出口所有物料的加和表示：因此，利用脉冲法可以很方便的测出停留时间分布密度。第39页，共96页，2022年，5月20日，3点55分，星期二脉冲法和阶跃法的比较脉冲法阶跃法示踪剂注入方法在原有的流股中加入示踪剂，不改变原流股流量将原有流股换成流量与其相同的示踪剂流股E(t)可直接测得F(t)可直接测得第40页，共96页，2022年，5月20日，3点55分，星期二研究不同流型的停留时间分布，通常是比较它们的统计特征

18、值。常用的特征值有两个：数学期望平均值方差离散程度平均停留时间它是指整个物料在设备内的停留时间，而不是个别质点的停留时间。不管设备型式和个别质点的停留时间，只要反应体积与物料体积流量比值相同，平均停留时间就相同。三、停留时间分布的数字特征第41页，共96页，2022年，5月20日，3点55分，星期二E(t)dt=dF(t)F(t)：所有停留时间为0-t的质点所占的分率F(t+dt)：所有停留时间为0-t+dt的质点所占的分率dF(t)= F(t+dt)- F(t)dF(t)：所有停留时间为t-t+dt的质点所占的分率数学期望：所有质点停留时间的“加权平均值”（4-17）,（4-18）第42页，

19、共96页，2022年，5月20日，3点55分，星期二在等时间间隔取样时：tt1t2t3.E(t)E(t1)E(t2)E(t3)对于离散型测定值，可以用加和代替积分值（4-19）第43页，共96页，2022年，5月20日，3点55分，星期二方差是停留时间分布离散程度的量度方差越小，越接近平推流对平推流，各物料质点的停留时间相等，故 方差为零。方差：各个物料质点停留时间t与平均停时间 差的平方的加权平均值。（4-20）第44页，共96页，2022年，5月20日，3点55分，星期二如果是离散型数据，将积分改为加和：取样为等时间间隔时：（4-20）第45页，共96页，2022年，5月20日，3点55分

20、，星期二平均对比时间：停留时间为t时， ，因此，和t一一对应，且有：F（ ）=F（t），此时：归一性：对比时间（无因次时间）：（4-22）第46页，共96页，2022年，5月20日，3点55分，星期二用表示的方差：（4-24）第47页，共96页，2022年，5月20日，3点55分，星期二描述停留时间分布的两个函数：（4-24）第48页，共96页，2022年，5月20日，3点55分，星期二四、理想流型反应器停留时间分布平推流模型全混流模型第49页，共96页，2022年，5月20日，3点55分，星期二1、平推流：“激励”曲线“响应”曲线图4-5 理想置换反应器“激励”与“响应”曲线两条曲线完全一样

21、，只是“响应”曲线比“激励”曲线平移了一段 时间 。第50页，共96页，2022年，5月20日，3点55分，星期二所有物料质点的停留时间都相同，且等于整个物料的平均停留时间tm，停留时间分布函数与分布密度为：由方差定义，第51页，共96页，2022年，5月20日，3点55分，星期二设进行阶跃注入实验，反应器的容积为VR，物料的体积流量为V，达到稳态后，从t=0开始，将进料切换为含示踪剂浓度为C0的物料，在切换后某dt时间内，对全釜作物料衡算：进入的示踪剂量=流出的示踪剂量+示踪剂的积累量全混流第52页，共96页，2022年，5月20日，3点55分，星期二给全混流反应器一个阶越“激励”“激励”曲

22、线“响应”曲线 全混流反应器“激励”与“响应”曲线第53页，共96页，2022年，5月20日，3点55分，星期二 对反应器做示踪物的物料衡算：第54页，共96页，2022年，5月20日，3点55分，星期二例：现有一个停留时间相同的全混釜和平推流反应器当二者以不同的方式进行串联时，其各自的最终转化率为多少？已知：停留时间为1min，反应速率常数为1.0m3/kmol.min,液相反应物浓度为1kmol/m解：全混釜在前，平推流在后全混釜停留1min后反应物的浓度：第55页，共96页，2022年，5月20日，3点55分，星期二反应物浓度为cA1后进入平推流在平推流停留1min后反应物的浓度：平推流

23、在前停留1min后反应物的浓度：第56页，共96页，2022年，5月20日，3点55分，星期二在全混釜中再停留1min后反应物的浓度：第57页，共96页，2022年，5月20日，3点55分，星期二 五、利用RTD诊断反应器内流动状况第58页，共96页，2022年，5月20日，3点55分，星期二第59页，共96页，2022年，5月20日，3点55分，星期二实际反应器中的流动状况总是偏离理想流动很难建立其真实方程可以先建立一种非理想流动模型，用它来描述实际反应器中的流动情况再通过对模型参数估值来确定偏离理想流动的具体程度 常用模型有：层流模型；轴向扩散模型； 多级串联全混流模型；组合模型非理想流动

24、是指物料在反应器内有部分返混发生，目前对非理想流动一般借助于模型加以描述。第三节 非理想流动模型第60页，共96页，2022年，5月20日，3点55分，星期二 CSTR(0维) PFR(一维) PD模型(一维) 径向扩散PFR(准二维)CoCC0CC0CC0CC1径向射流(圆柱形三维)四角射流(异形三维)涡流V(x,y,z)C(x,y,z)T(x,y,z)圆柱座标直角座标1、反应器流动模式 第61页，共96页，2022年，5月20日，3点55分，星期二2、停留时间分布的应用确定模型参数m或Pe用多级串联全混流模型或轴向扩散模型模拟实际反应器中的流动状况，关键是确定串联的级数m或Pe， m或Pe

25、又与方差有关，因此，可以通过实验确定停留时间分布，进而计算方差， m或Pe，然后求得转化率。定性分析流动状况活塞流全混流定量分析流动状况实际反应器中可能存在短路与死角，使实际的平均停留时间不等于VR/V，因此可以得用停留时间分布来定量估算死角与短路的程度。第62页，共96页，2022年，5月20日，3点55分，星期二例：某全混流反应器体积为100L，物料流率为1L/s，试求在反应器中停留时间为（1）90110s，（2）0100s，（3）100s的物料占总进料的比率。 解：出口物料的份额用F(t)表示， 所求比率90110s ：F(110) F(90) = 0.074 = 7.4% 小于平均停留

26、时间的物料占63.2% 大于 的物料占总物料的36.8% 第63页，共96页，2022年，5月20日，3点55分，星期二 实际反应器内的流动是很复杂的，仅用理想化的平推流或全混流进行计算是不够的，非常必要对实际的流型进行逼近模拟。常用的方法有如下几种。 第64页，共96页，2022年，5月20日，3点55分，星期二填料塔层流模型可对平推流模型进行修正，但比较粗糙。主要模型假设为：1) 反应混合物在管内作层流流动。2) 反应混合物在轴径向均不存在返混。 u第65页，共96页，2022年，5月20日，3点55分，星期二轴向混合模型在平推流模型的基础上再迭加一个轴向混合轴向扩散模型的基本假设：第66

27、页，共96页，2022年，5月20日，3点55分，星期二 流体沿轴向有参数变化，径向参数均一； 流体流动主流为平推流但叠加一逆向涡流扩散； 逆向涡流扩散遵循菲克定理，但整个反应器内扩散系数为一常数。第67页，共96页，2022年，5月20日，3点55分，星期二轴向扩散模型仿照分子扩散中用扩散系数分子运动表征有效扩散系数表征一维返混用平推流反应器+涡流扩散运动非理想流动的管式反应器主要用于湍流流动的管式反应器、固定床反应器、塔式反应器。第68页，共96页，2022年，5月20日，3点55分，星期二轴向混合模型对实际反应器，处理时在平推流的基础上迭加一个轴向混合来进行校正。适合于不存在死角、短路和

28、循环流、返混程度较小的非理想流动模型。 模型参数是轴向混合弥散系数EZ，停留时间分布可表示为EZ的函数。Peclet准数：当返混程度相当小时；数学期望=1 方差对一级不可逆反应，转化率可表示为：第69页，共96页，2022年，5月20日，3点55分，星期二第70页，共96页，2022年，5月20日，3点55分，星期二0 0.5 1.0 1.5 2.00.51.01.52.00.050.10.5110第71页，共96页，2022年，5月20日，3点55分，星期二0 0.5 1.0 1.5 2.00.20.40.60.81.00.0010.0120.0250. 25第72页，共96页，2022年，

29、5月20日，3点55分，星期二基本假设： 由m个体积相等的CSTR串联组成；从一个CSTR到下一个CSTR之间的管道内物料不发生反应。多级全混流串联模型多釜串联模型:用几个等体积的全混流反应器串联来模拟实际反应器中的流动状况。第73页，共96页，2022年，5月20日，3点55分，星期二0 0.5 1.0 1.5 0.51.0m=1m=2m=5m=10m0 0.5 1.0 1.5 0.51.0m=1m=2m=5m=10m图4-11 不同m值时的E(), F()曲线第74页，共96页，2022年，5月20日，3点55分，星期二每一级的停留时间ti=tm/m。模型参数为串联级数m。方差m=1时，

30、即为全混流模型m=时， 即为平推流模型对一级不可逆反应，转化率可表示为：假设实际反应器中的返混程度与m个等体积的全混流反应器串联时相同，m是虚拟釜数，不一定是整数。第75页，共96页，2022年，5月20日，3点55分，星期二 当m=1时，多级混合釜模型的停留时间分布函数与CSTR相同。当m时，多级混合釜模型的停留时间分布函数与PFR相同。当1m 时，多级混合釜模型的停留时间分布函数属于非理想流动反应器。第76页，共96页，2022年，5月20日，3点55分，星期二多级混合釜模型的应用1、多个CSTR串联操作时，当釜数由一增加到无限多个，其停留时间分布规律将按CSTR 到非理想流动反应器到PF

31、R的分布规律变化。2、如果釜数选择恰当，该模型的停留时间分布规律可以与任何一个非理想流动反应器的停留时间分布规律相同（即两者的返混程度相同）。3、可以用多级全混流模型计算任何一种非理想流动反应器的转化率。反应釜个数为：第77页，共96页，2022年，5月20日，3点55分，星期二组合模型适用于上述两种模型不能很好表达的情况将实际反应器的流动情况设想为平推流、全混流、死区、短路、循环流等部分组成组合模型的几种典型例子：第78页，共96页，2022年，5月20日，3点55分，星期二第五节 非理想流动反应器的计算准确地计算实际反应器的生产能力和转化率，需已知两个前提条件： 1) 本征动力学方程 2)

32、 反应器停留时间分布函数如实际反应器中诸微元具有独立身份，则可以想象为实际反应器由不同长度管式反应器并联组成： 入口出口为各并联反应器转化率的积分平均。第79页，共96页，2022年，5月20日，3点55分，星期二出口转化率 停留时间为ti的转化率 ti的流量分量 即或无论何种反应器，只要已知停留时间分布函数，即可接上式计算 出口第80页，共96页，2022年，5月20日，3点55分，星期二对全混流对2级动力学第81页，共96页，2022年，5月20日，3点55分，星期二对平推流对2级动力学第82页，共96页，2022年，5月20日，3点55分，星期二例：应用脉冲示踪法测定一容积为12L的反应装置，进入反应器的流体速度0.8L/min，在定常态下脉冲地输入80g示踪剂A，并同时在反应器出口处记录cA随时间的变化，其实测数据列于下表: t(min)05101520253035 cA(g/L)03554210若用其进行某一级液相不可逆反应(-rA)=k cA,k=0.307min-1,当平均停留时间为15min时，试求：（1）采用多级全混流模型时，其最终转化率为多