分离工程II--07多级分离计算一三对角矩阵法选编.pptVIP

* 式(7-71a)则保持原来各项之值，可写为下列形式 ： (7-72a) (7-72c) * 由公式可知，松弛法是一种操作型算法，要预先规定平衡级数N及各进料位置，此外，还必须预先设定汽液相流率(Vj及Lj)，温度(Tj)及组成(xji)的初值。流率初值按恒摩尔流假设进行估算，温度初值在塔顶温度及塔釜温度间做线性内插，而各组分各级的组成初值，一律取进料平均组成，按下式算出： (7-73) * 松弛系数?j先按下式假定： 式中，分母(2~5)?Fj表示全部进料量的2~5倍间的适当值。 ?j值在运算过程中重新调优。 松弛法计算各级温度用泡点法。 汽液相流率由热量平衡法求之。第j级热平衡式为： (7-74) (7-75) * 根据第j级的物料平衡，有： 将式(7-76)代入(7-75)中，消去Lj并整理可以得到： 由第j级向上至第1级做物料平衡，得： (7-76) (7-76) (7-77) * 式(7-76)及式(7-77)就是计算汽液相流率的递推公式，将L1=RD，V2=D(R+1)代入式(7-76)算出V3，将V3代入式(7-77)算出L2，如此交替使用式(7-76)及式(7-77)，顺次得出V4，L3，V5，L4，…，VN，LN-1，而LN可用下式得到： (7-78) * 松弛法计算步骤如下： (1). 输入原始数据，设定Vj，Lj，Tj，xji初值，算出Kji及yji，初设?j值。 (2). 用式(7-72)，逐级计算出新的xji值，计算中，若某一级的xji产生负值，用减小?j的办法改正它。 (3). 用下列判断式(7-79)及式(7-80)，判断是否收敛，如果满足两式则输出计算结果，否则向下进行。 (7-79) * (4). 将xji(k+1)逐级圆整，得到圆整后的xji(k+1)值。 (5). 逐级求泡点算出新的Tj，同时求得新的Kji及yji，用式(7-76)及(7-77)算出新的Vj，Lj。 (6). 如果不满足收敛条件，返回第(2)步，重复计算。 (7-80) * 典型精馏计算法优缺点比较 算法 优点 缺点 泡(露)点法 a. 算法简单 b. 对初值要求不高 c. 占用内存较小 a. 当临近解时，出现收敛速度慢，迭代次数过多。 b. 算法包含泡(露)点计算，因它也是迭代过程，费机时 c. 对高度非理想系难收敛 SR法 a. 算法中不包含泡(露)点计算，收敛速度快。 b. 适用吸收塔及宽沸程精馏。 a. 不适用于窄沸程精馏系统。 b. 板数较多时稳定性差。 N-R法 a. 适用于宽、窄沸程的精馏系统及吸收塔 b. 不包含泡(露)点计算 c. 收敛速度快，稳定性好。 a. 对初值要求高。 b. 计算雅可比矩阵，耗时长 c. 占内存单元多。 * 续表： 算法 优点 缺点 松弛法 a. 适用范围宽 b. 对初值要求不高 c. 稳定性好 收敛速度慢，不宜用做常规算法。 * The end * * 当用于表示液液平衡关系时，以代表了高密度物流即比重大的液体，以代表低密度物流即比重小的液体。 * * 3. 用托马斯算法求解出液相摩尔分率xji，不经圆整，直接代入下式（7-27），计算出液相流率Lj。 相应的汽相流率可由总物料衡算式(7-28)求得： 然后直接迭代至收敛。 4. 由H方程校正Tj。具体方法是： 令j级热平衡的偏差函数为： (7-27) (7-28) (7-29) * 按牛顿近似法将Ej展开成泰勒级数且只保留一阶偏导数得： 对j=1： 对2? j ?N-1 (7-30) (7-30a) (7-30b) * 对j=N 对于2? j ?N-1，对照式(7-30)与式(7-30b)并略去后一式中的上标，则有： (7-30c) (7-31a) (7-31b) (7-31c) * 式中，偏导数 和 实际上表示汽、液相在温度T时的比热cV，cL，对于理想混合物： 式中，ci为纯组分的比热；角标V，L分别代表汽、液相。 所求的Tj应满足Ej(k+1) = 0。由式(7-31)和(7-32)可写出类似ME方程的形式： (7-32a) (7-32b) * 同三对角矩阵法的泡点法相类似，用追赶法可解 (7-33) * 出Tj；第k+1次同第k次Tj的关系为： 流率加和法液可以用类似前面所介绍的?加速迭代法，对第一级组分流率都用校正因子，其定义是： 三对角矩阵解出lij，vij后用?j校正，得出新的lij，vij再加和求Lj，Vj。 (7-34) (7-35) * 输入物性数据及N，进料位置 假设Tj(0), Vj(0), xjj(0), 计算组分