化工过程合成和分析.pptVIP

隐式方程具有更大的普遍性，所以割线法常为人们所熟知。 在流程模拟领域中，物流回路多用显式方程描述的。因而多用Wegstein法。 由（2）式可见，一维Wegstein法需要有两个初值，其中第一个初值是设置的猜值，第二个初值可根据第一个初值按直接迭代法得到。 （3）回路矩阵 要断裂再循环物流，必须先识别再循环回路，并借助一定的方法描述它们。 一个不可分隔子系统包含若干个再循环回路。包含两个以上再循环物流，且其中的任何单元只被通过一次，称作简单回路（Simple Cycle）。 Ⅰ－S1－Ⅱ－S2－Ⅲ－S4－Ⅱ－S2－Ⅲ－S5－Ⅰ构成的回路不是一个简单回路，因为其中的单元Ⅱ和单元Ⅲ被通过了两次。 图2-12 含有四个简单回路的不可分隔子系统 s7 s1 s5 Ⅰ Ⅱ Ⅲ Ⅳ s6 s4 s2 s3 过程系统中的简单回路可以用回路矩阵（Loop/Stream Matrix）表示。 表示方法： 矩阵中的行代表回路，列代表物流。 若某回路I中包括有物流j则相应的矩阵元素aji=1，否则为空白或零。 （4）Upadyhe－Grens断裂法 美国加州大学的Upadhye等提出的，一种类似动态规划法的寻求最佳断裂物流的算法。 为了对该不可分隔子系统的高维求解进行降维运算，须将该子系统中的某些回路进行断裂。达到断裂的方案并不是唯一的。 需要解决的两个问题：一是要有一种能把所有的有效断裂物流组都能搜索出来的办法；二是要能把最优断裂组从中选择出来。 图2-13 不可分隔子系统 6 4 2 3 7 1 5 单元1 单元2 单元3 单元4 (2) (3) (4) (9) (3) (2) (2) 流股序号 （权重因子 ） W f 替代规则： 令｛D1｝为一有效断裂组，Ai为全部输入流均属于｛D1｝的单元（至少有一个这样的单元存在，否则｛D1｝为无效断裂组）。将Ai的所有输入流用Ai的全部输出流替代，构成新的断裂组。令得到的新的断裂组为｛D2｝，则 ①｛D2｝也是有效断裂组； ② 对于直接迭代，，｛D2｝与｛D1｝具有相同的收敛性质。 由替代规则联系起来的所有断裂组的集合，定义为断裂族 （5）寻求最优断裂组的算法 ① 从任一有效断裂开始，运用替代规则： ② 如果在任何一步中出现有两次被断裂的物流（二次断裂组），则消去其中的重复物流。消去重复后断裂组则作为进行下一步的新起点。 ③ 重复步骤①、②，直到不再有二次断裂组出现，且每个“树枝”上有重复的断裂组出现时为止。从最后一个新的起点开始，其后出现的所有不重复的断裂组成为非多余断裂族。 ④ 非多余断裂族中总权最小的断裂组为最优断裂组。 s5,s6 s1,s4,s7 s3 s5,s6 s2 图2-14 替代过程 s1,s2,s3 s1,s4,s5,s3,s3 { s1,s4,s5,s6,s7 } { s1,s2,s6,s7 } s1,s4,s7,s1 { s2,s6,s7,s5,s6 } { s2 } { s3,s4,s5 } { s1,s4,s7 } { s4,s5,s6,s7 } s2 s3 s3 s1,s4,s7 非多余断裂族: 断裂组｛S1，S4，S7｝为最优断裂组。 通过断裂可以把不可分隔子系统中的回路物流打开，从而可以利用序贯模块法对该过程系统进行模拟计算。这种模拟计算的开始是首先要设定起始物流变量的猜值，计算的终点则在于该猜值与计算值的收敛。 2.2.3 断裂物流变量的收敛 执行断裂物流变量收敛功能的模块称收敛单元模块 不可分隔子系统的断裂物流 断裂物流变量的收敛问题实际上是个迭代求解非线性方程组的问题 x = y = G(x) 图2－15　 不可分隔子 系统的断裂 物流 过程系统 (a) 断裂流股 (b) 当断裂物流变量猜值x与计算值y之差小于收敛容差ε时 y － x = G（x）－ x ε 则x为断裂物流变量的收敛解 收敛单元的功能总计有如下作用： (Ⅰ)获取猜值的初值x0 (Ⅱ)根据计算值y，以一定的方法确定新的猜值x (Ⅲ)比较猜值x和计算值y，若其结果满足给定精度要求则结束迭代计算，否则继续迭代计算过程 可见，收敛单元实质上就是一个数值迭代求解非线性方程组的子程序。求解非线性方程组的数值计算方法很多，适合于收敛单元的数值计算方法一般应尽可能满足下列要求： 1. 对初值的要求不高。 1).初值易得，不易引起迭代计算的发散； 2)初值的组数少。 例: 用直接迭代法求解下列方程组 解：令猜值为X1＝2；X2＝10；X3＝5 解：令猜值为X1＝6；X2＝3.5；X3＝5 2. 数值稳定