第2讲：精馏严格计算.pptVIP

精馏计算 精馏计算 2.3.1精馏的定态数学模型 塔板模型 热量衡算H方程的讨论－ Heat Balance 2.3.2 精馏定态模拟计算算法概述 非线性方程组求解要点 2.3.3 三对角线矩阵算法 W-H法的计算框图 W-H的求解方法: 2. 流率加和（SR）法（弗莱迪，史密斯） 流率加和(SR)法计算框图 BP法与SR法的算法特点及适用场合 2.3.4 全变量迭代法（Naphtali-Sandhom） MEH方程的剩余函数形式 迭代变量与函数的排列顺序： 2．联立解一般方法 全变量迭代算法的特点： 总结： 课 堂 练 习 * 精馏 精馏严格计算 多元精馏与二元精馏的区别： xDA，xDB xFA，xFB xBA，xBB 苯－甲苯 xDA，xDB，xDC xFA，xFB，xFC xBA，xBB，xBC 苯－甲苯－二甲苯 分离要求： xDA=0.99，xBB=0.985 多元精馏计算主要采用的方法是: 设计型计算—给定分离任务，选择工艺条件； 2．操作型计算—设备给定，进行工艺条件核算。 定态数学模型： 多级平衡级模型 多级非平衡级模型 M方程－Material Balance E方程－Phase Equilibrium S方程－Mole Fraction Summation H方程－Heat Balance Subscript denotes tray number Tj Pj Qj Uj Wj Vj yj _ Lj-1 xj-1 _ Lj xj _ Vj+1 yj+1 _ Fj XF _ Overbar denotes component vectors: e.g., x = (x1, x2, ...xNC) _ _ Lj , Vj Fj Qj xj , yj XF hj , Hj Tj , Pj Uj Wj Liquid, vapor flowrate Feed flowrate Heater/cooler duty Liquid, vapor mole frac Feed mole fractions Liquid, vapor enthalpies Temperature, pressure Liquid Draw rate Vapor Draw rate _ _ j＝1(塔顶冷凝器) 以物料衡算M方程为例讨论－Component material balance 对 i 组分物料衡算： j＝2，…，N-1塔板 小写字母－i 组分的流率； 大写字母－总流率 j＝N(塔釜再沸器) Lj-1 xi,j-1 Lj xij Uj Qj Vj+1 yij+1 Wj Fj zi Vj yi,j Lj-1 xi,j-1 hj-1 Lj hj xij Uj Qj Vj+1 Hj+1 yij+1 Wj Fj zi HFj Vj Hj, yi,j MESH方程 M方程 E方程 S方程 H方程 变量分析 自由度数＝（6C+11)N-2C-1-(5C+6)N-2C=(C＋5)N-1 (6C+11)N-2C-1 (5C+6)N-2C 总计 1 板数N N Qj (N-1)C (2-110)－u N Hj N hFj (N-1)C (2-109)－w N hj N Fj N (2-108)－H N Tj NC fij N (2-107)－h N Pj NC Kij NC (2-106)－K N-1 Wj (N-1)C wij N (2-105)－H-B N-1 Uj (N-1)C uij 3N (2-104)－S N Vj NC vij NC (2-103)－E N Lj NC lij NC (2-102)－M 变量数 变量 变量数 变量 方程数 方程编号 若N=50，C=3，自由度＝399 操作型计算一般命题，通常 精馏模拟计算的算法 1）迭代变量的选择； 2）迭代算法的组织；（分块解还是联立解） 3）变量圆整或归一及迭代加速方法的选择。 设有如下非线性方程组： 1、分块求解 将方程用泰勒级数在初值点 展开，并取一阶导数的近似值，得： 2、联立求解（即牛顿－拉普生法） 写成矩阵形式： （5) 求解步骤： f t 结束 由于精馏系统的复杂性，MESH方程也有两种解法。 1.W-H法的解题思路(BP法) 将MESH方程分成三块，M+E方程，S方程，H方程; M+E方程，S方程，内层迭代求解，与 lij，Tj 匹配; H方程在外层迭代求解，与 Lj(Vj) 匹配。 收敛判据： 1) ME方程 （对塔内所有板上的某 i 组分作物料衡算） M+E方程