串联CSTR动态模型仿真.docVIP

串联CSTR的动态仿真 建模 图1 CSTR）工艺流程图 图1所示[1~4]。在两个反应器里都发生着一个单独的一阶不可逆放热反应A→B[]，中间混合槽用来引进第二种进料在建模时作如下简化假定：混合物的密度为常数，与其温度和浓度无关。热容也为常数，与其温度和成分无关两个CSTR的液体体积为恒值。: CSTR中的物料平衡方程： CSTR中的能量平衡方程： CSTR的体积恒定： 每个反应器被加在其周围的冷却罩冷却，单位时间内反应器与冷却罩之间的传热为 忽略混合槽动态行为，物料和能量平衡下 联立上述方程，解 (1.1) 式中 符号说明 C —— 浓度，mol m-3 Cool —— 冷却量，W cP —— 比热，J m-3 K-1 Dh —— 反应热，J mol-1 E —— 活化能，J mol-1 K —— 化学反应速率 k0 —— 化学反应速率常数，s-1 Q —— 流量，m3 s-1 R —— 理想气体常数，8.314 J mol-1 K-1 T —— 温度，K Ua —— 传热系数，W K-1 V —— 反应器内反应物的体积，m3 上下角标 C —— 冷却 ci —— 进入的冷却液 cw —— 进入的冷却液 f —— CSTR进料 0 —— 外界进料 1,2 —— CSTR编号 常数项参数值 符号 数值 说明 R 8.314 理想气体常数，J mol-1 K-1 E 6000*R 活化能，J mol-1 k0 2.7*10^8 化学反应速率常数，s-1 Ua 0.35 传热系数，W K-1 Dh 5 反应热，J mol-1 cP 1 反应物比热，J m-3 K-1 cpc 1 冷却水比热，J m-3 K-1 变量稳态值及其波动上下限 符号 稳态值 波动下限 波动上限 说明 物料1 Q01 0.339311 外界进物料1的流量，m3 s-1 T01 300 295 320 外界进物料1的温度，K C01 20 19.5 22 外界进物料1中A的浓度，mol m-3 物料2 Q02 0.216402 外界进物料2的流量，m3 s-1A T02 300 295 320 外界进物料2的温度，K C02 20 19.5 22 外界进物料2中A的浓度，mol m-3 冷却液1 Qcw1 0.35 进入CSTR1冷却液的流量，m3 s-1 Tci1 250 进入CSTR1冷却液的温度K 冷却液2 Qcw2 0.8 进入CSTR2冷却液的流量，m3 s-1 Tci2 250 进入CSTR2冷却液的温度K 中间产品 Q1 反应器1出的中间产品流量，m3 s-1 T1 反应器1出的中间产品温度，K C1 反应器1出的中间产品中A的浓度，mol m-3 产品 Q2 反应器2出的产品流量，m3 s-1 T2 反应器2出的产品温度，K C2 反应器2出的产品中A的浓度，mol m-3 状态变量 V1 5 4.5 5.5 反应器1内反应物的体积，m3 V2 5 4.5 5.5 反应器2内反应物的体积，m3 2利润与生产约束 过程的目标函数为整个流程的净利润： J= 产值- 冷却水成本- 原料成本 maxJ= 150. 00 (Q f1Cf1+ Q f2Cf2- Q 2C2)- 0. 15 (V c1+ V c2) - 1. 50 (Q f1+ Q f2)式中150. 00 （￥/mol）、0.15(￥/ m3) 和1.5(￥/ m3)分别为产品、冷却水和原料价格。式中计算的是每秒的利润。利用M 软件的Optimization工具箱得静态最优值J 0= 1563. 08元秒 过程的目标函数为整个流程的净利润： J= 产值- 冷却水成本- 原料成本 (3.2) 优化变量为第一个和第二个进料流量(Q01,Q02)，两股进料的温度和组分可认为是干扰(T01, C01, T0, C02)，其下限、标称值、上限分别为T01(295,300,320), C01(19.5,20,22), T02(295,300,320), C02(19.5,20,22)，模型中的不确定参数为传热系数Ua，其范围为Ua(0.30,0.35)，其中，上限值为标称值。为了避免失控和生成不希望的副产品，这个过程受到的约束如下两个CSTR的温度（即两个CSTR出口产品的温度）限制T1≤350K、T2≤350K；两个CSTR出口产品中原料A的含量限制C1≤1mol m-3、C2≤0.3 mol m-3；进料限制Q01≥0.05m3 s-1、Q02≥0.05m3 s-1、Q01+Q02≤0.8m3 s-1；