CSTR反应器开工过程动态模拟 指导教师：小组成员： 完成时间：2013年4月1日 目录引言 . 1问题的提出 . 2模型的建立 .2换热条件的确定.4连续反应器的开车与计算.5开车失败的原因分析.6解决方案 . 8CSTR 反应器稳定性分析 .10结论. .12附录 . .13引言 模型化 (Modeling)是现代化学工程方法论的重要组成部分，尤其是过程动态学的核心。 根据对过程系统中状态变量分布特征的不同描述方式，一般可以把数学模型分为集中参数模型、分布参数模型和多级集中参数模型；根据建立模型的不同方法，一般可以将数学摸型分为统计模型、确定性模型和介于二者之间的半经验模型 。 CSTR（Constant Str Tank Reactor）,即全混流反应器，又称全混釜或连续流动充分搅拌槽式反应器，是一类在工业生产中广泛应用的反应器。CSTR的特性有：（1）物料在反应器内充分返混；（2）反应器内各处物料参数均一；（3）反应器的出口组成与器内组成相同；（4）连续、稳定流动，是一定态操作。 针对CSTR反应器的特性，可以通过对所研究的反应器系统与系统内的关键组分，列出其质量、能量和动量守恒关系式，系统及关键组分内外质量、能量和动量交换速率系数计算式，相关的相平衡关系，以及化学反应速率表达式和化学反应平衡常数计算式（如果反应可逆的话），从而建立CSTR反应器内反应的确定性集中参数模型，以此对在CSTR内进行的反应进行描述。 要利用确定性模型来预测、揭示化工过程系统的内在规律，解决与动态学特征有关的工程实际问题，就必须对所涉及到的数学模型进行必要的数学处理。以此对CSTR反应器建立的确定性集中参数模型的数学表达式所构成的数学模型的正问题，涉及到代数方程组、常微分方程组和偏微分方程组，以及他们的混合方程组的求解，由于化工过程通常具有的非线性特征，求模型方程组的解析解往往是不可能的，不得不借助计算机求数值解。对于典型的常微分方程组的初值问题，通常可以利用龙格库塔（R-K）、基尔(Gear)法等通用程序求数值解。本篇阐述了我们通过文献调研，查取物性数据及相关资料，对一个在CSTR反应器内进行的二级不可逆放热反应进行分析、建立起反应的确定性集中参数模型，应用相关计算法则，编制计算机程序，并利用计算机辅助设计绘制了状态空间的相平面图，实施了对这一反应的状态空间分析和热平衡分析，并对所得结果进行讨论和综合的过程。在整个过程中，我们始终坚持用系统工程的思想、方法来解决化工过程系统的设计、开发、操作、控制等问题，取得了可喜的成绩，使我们加深了对化工过程分析与合成的认识以及对化工过程系统动态模拟与分析的理解，为我们在今后的学习和工作实践奠定了思想基础。一 、 问题的提出 工业上生产1,2-丙二醇主要采用环氧丙烷（PO）催化水合的方法，为一级不可逆放热反应，该方法以硫酸为催化剂，其化学反应方程式为： 该反应为一级反应： ，此反应速率常数为： 在反应温度范围内，可以假定反应热为常数，消耗每摩尔PO放出83.74 kJ 热量。根据化学反应工程的相关原理可知，在连续搅拌釜式反应（CSTR）中，该反应体系可能存在多定态现象。借助反应体系的多定态特性进行催化反应器的开车和停车，可为催化反应器的安全操作提供重要的指导，因此对其进行研究，对于加深对化学反应体系多定态体系的认识是十分必要的。 二 、 模型的建立（一） 基本假设1. 反应釜内处于分子级理想混合，且为液相均相反应，因此反应釜内混合物的温度和组成在反应区里是均匀的。2. 反应区的容积不随时间变化，则加料与排料的流量也认为是相等的，即 3. 反应釜内压力恒定。4. 反应物在进料管中不反应。5. 本反应不存在副反应，即反应只生成1,2-丙二醇，在反应釜中最多只有4种物质混合。6. 忽略混合热，混合物的物性为各物质物性的加和平均，也即是 （二） 衡算方程只着眼于A 组分：（1） .质量守恒方程为： (1)（2） 能量守恒为： (2)（3） 冷却水出口温度的计算 ： (3) 式中 (4) （三） 符号说明： 反应器示意图如右图所示： A 表观频率因子， Q 换热量，J Cp 恒压热容，J 平均停留时间，h Ea 活化能，kJ 其它物质与A 的进料物质的量比 F 摩尔流率，mol A 环氧丙烷 化学反应热，kJ B 水 k 化学反应速率常数， C 丙二醇 冷却水流量，kg eb 能量衡算 N 进料量，mol i i 种物质r 化学反应速率，mol mb 质量衡算R 气体常数，J M 甲醇t 时间，h s 稳态T 温度，K 1 进口冷却水UA 总传热系数与换热面积之积，kW 2 出口冷却水v 体积流量，m3 0 初始状态V 反应器体积，m3 X 转化率三 、换热条件的确定 （一）说明：1 由于反应放出热量太少，另外，所选取的夹套换热方式的弊端，衡算时采用一些经验公式；而且采用试差法进行计算；2 由于模拟是动态过程，所有相关的参数采用平均植；3 反应器与夹套换热装置为钢制；4 忽略污垢热阻； （二）热量衡算式： (5) （三）换热计算：牛顿冷却定律： (6)总对流传热系数： (7)反应釜内传热速率方程： (8)夹套侧传热系数（由于Re4400） (9)四 、连续反应器的开车与计算 目前国际上丙二醇单套装置的生产能力相差很大，如美国的陶氏化学和Lyondell 公司的丙二醇单套装置生产能力在1105 t/y 以上，而国内丙二醇单套装置年生产能力则在几千吨到几万吨不等。以国内投产的中海油壳牌6104 t/y 1,2-丙二醇项目为例，按年平均开工8000 h，平均转化率80.00%计算。 在PO 进料之前，反应器内应加入一定量质量分数为0.10%的H2SO4 水溶液。为保证PO 高的转化率，水溶液进料量应大于PO 进料量。经计算可知进料需由7.9 t/h 的PO，16.2 t/h 质量分数为0.10%的水溶液以及0.3 t/h 的甲醇组成。进料温度为297.15 K，冷却水入口温度为302.15 K，反应器体积为1.60 m3。据相关文献可知，环氧丙烷（A）、水（B）、丙二醇（C）和甲醇（M）的物质的量浓度分别为14.93，55.26，13.68 和24.67 kmol/m3。在297.00350.00 K 内，各物质的恒压热容可视为定值：CpA=146.53 J/(molK)，CpB=75.36 J/(molK)，CpC=192.58 J/(molK)，CpM=81.64 J/(molK)。 由于该反应为强放热反应，因此CSTR 应配备冷却器。取冷却器的换热能力UA 为2.90 kW/K。对于大的冷却介质流量值，换热量Q 的计算公式可简化为：