第二章化工系统的定常态模拟与分析2.ppt

1、2.2 序贯模块法,2.2.1 序贯模块法的基本原理,序贯模块法的基础是单元模块，通常单元模块与过程单元是一一对应的单元模块是依据相应过程单元的数学模型和求解算法编制而成的子程序单元模块具有单向性特点。,F,序贯模块法的基本思想,从系统入口物流开始，经过接受该物流变量的单元模块的计算得到输出物流变量，作为下一个相邻单元的输入物流变量。依此逐个的计算过程系统中的各个单元，最终计算出系统的输出物流。计算得出过程系统中所有的物流变量值，也即状态变量值。,通常把输入物流变量及单元模块参数（回流比、几何尺寸等)作为决策变量。序贯模块法的求解与过程系统的结构是有关的。具有反馈联结的系统（不可分割子系统），

2、需要用到断裂（Tearing）和收敛（Convergence）技术。,图2-9 具有反馈的系统与收敛单元,2,3,5,收敛单元位置的选取选取新的S4数值，保证计算收敛，加速收敛,2.2.2 再循环物流的断裂,（1）断裂的基本概念,2.2.2 再循环物流的断裂,（1）断裂的基本概念,求解方法,1.联立方程组2.通过断裂进行降维 X2 X3 X4 X1 X2 把一个四维求解问题降阶成为了四个一维问题，从而减化了计算难度。这种通过迭代把高维方程组降阶为低维方程组的办法称为“断裂”。,求解方法,1.联立方程组2.通过断裂进行降维 X2 X3 X4 X1 X2 把一个四维求解问题降阶成为了四个一维问题，

3、从而减化了计算难度。这种通过迭代把高维方程组降阶为低维方程组的办法称为“断裂”。,有向图,断裂点,图2-10 有向图,选择不同的断裂物流，相应的迭代序列也不一样：,图2-11 不可分隔子系统,由于系统中各物流及其变量特性的不同，在收敛计算上常是有很大差异的。如何选择断裂物流、确定迭代序列，是实施序贯模块法进行过程系统模拟计算中必须要解决的问题。,（2）断裂方法的研究,六十年代初，Rubin就提出了断裂的思想判断最佳断裂的准则分为四类 断裂的物流数最少；断裂物流的变量数最少；断裂物流的权重因子之和最少；断裂回路的总次数最少。,（2）断裂方法的研究,六十年代初，Rubin就提出了断裂的思想判断最佳

4、断裂的准则分为四类 断裂的物流数最少；断裂物流的变量数最少；断裂物流的权重因子之和最少；断裂回路的总次数最少。,j,权重因子；i=1，m，代表回路；j=1，n,代表物流,约束方程的含义是每个回路至少要被断裂一次。准则设定j=1；准则令j为物流变量数；准则中j为可根据物流性质而取的选择值，如物流变量对计算过程灵敏度大小的估计值；准则的j等于每个断裂物流所切断的回路总数。,（3）回路矩阵,要断裂再循环物流，必须先识别再循环回路，并借助一定的方法描述它们。一个不可分隔子系统包含若干个再循环回路。包含两个及以上再循环物流，且其中的任何单元只被通过一次，称作简单回路（Simple Cycle）。,S1S

5、2S4S2S5构成的回路不是一个简单回路，因为其中的单元和单元被通过了两次。,图2-12 含有四个简单回路的不可分隔子系统,过程系统中的简单回路可以用回路矩阵（Loop/Stream Matrix）表示。表示方法：矩阵中的行代表回路，列代表物流。若某回路I中包括有物流j则相应的矩阵元素aji=1，否则为空白或零。图2-12中简单回路可用回路矩阵表示如下：,美国加州大学的Upadhye等提出的，一种类似动态规划法的寻求最佳断裂物流的算法。为了对该不可分隔子系统的高维求解进行降维运算，须将该子系统中的某些回路进行断裂。达到断裂的方案并不是唯一的。需要解决的两个问题：一是要有一种能把所有的有效断裂物

6、流组都能搜索出来的办法；二是要能把最优断裂组从中选择出来。,（4）Upadyhe-Grens,替代规则：,令D1为一有效断裂组，Ai为全部输入流均属于D1的单元（至少有一个这样的单元存在，否则D1为无效断裂组）。将Ai的所有输入流用Ai的全部输出流替代，构成新的断裂组。令得到的新的断裂组为D2，则,D2也是有效断裂组；对于直接迭代，D2与D1具有相同的收敛性质。,由替代规则联系起来的所有断裂组的集合，定义为断裂族,（5）寻求最优断裂组的算法,从任一有效断裂开始，运用替代规则：如果在任何一步中出现有两次被断裂的物流（二次断裂组），则消去其中的重复物流。消去重复后断裂组则作为进行下一步的新起点。重

7、复步骤、，直到不再有二次断裂组出现，且每个“树枝”上有重复的断裂组出现时为止。从最后一个新的起点开始，其后出现的所有不重复的断裂组成为非多余断裂族。非多余断裂族中总权值最小的断裂组为最优断裂组。,物流 权重因子,图2-13 不可分隔子系统,6,4,2,3,7,1,5,单元1,单元2,单元3,单元4,(2),(3),(4),(9),(3),(2),(2),权重因子,回路,流股,用Upadhye-Grens断裂法求图2-13中的最优断裂组,从断裂组（S1,S2,S3）开始，反复利用替代规则,S2,S1 S2 S6 S7,S3,S3,S1 S4 S5 S6 S7,S1 S4 S7,S2 S5 S6,

8、S5 S6,S1 S4 S7,S2,S3 S4 S5,S2,S3,S4 S5 S6 S7,S5 S6,S1 S4 S7,非多余断裂族:断裂组S1，S4，S7为最优断裂组。,通过断裂可以把不可分隔子系统中的回路物流打开，从而可以利用序贯模块法对该过程系统进行模拟计算。这种模拟计算的开始是首先要设定起始物流变量的猜值，计算的终点则在于该猜值与计算值的收敛。,物流 权重因子,图2-13 不可分隔子系统,6,4,2,3,7,1,5,单元1,单元2,单元3,单元4,(2),(3),(4),(9),(3),(2),(2),权重因子,回路,流股,用Upadhye-Grens断裂法求图2-13中的最优断裂组,

9、从断裂组（S1,S2,S3）开始，反复利用替代规则,S2,S1 S2 S6 S7,S3,S3,S1 S4 S5 S6 S7,S1 S4 S7,S2 S5 S6,S5 S6,S1 S4 S7,S2,S3 S4 S5,S2,S3,S4 S5 S6 S7,S5 S6,S1 S4 S7,非多余断裂族:断裂组S1，S4，S7为最优断裂组。,通过断裂可以把不可分隔子系统中的回路物流打开，从而可以利用序贯模块法对该过程系统进行模拟计算。这种模拟计算的开始是首先要设定起始物流变量的猜值，计算的终点则在于该猜值与计算值的收敛。,物流 权重因子,图2-13 不可分隔子系统,6,4,2,3,7,1,5,单元1,单元

10、2,单元3,单元4,(2),(3),(4),(9),(3),(2),(2),权重因子,回路,流股,用Upadhye-Grens断裂法求图2-13中的最优断裂组,从断裂组（S1,S2,S3）开始，反复利用替代规则,S2,S1 S2 S6 S7,S3,S3,S1 S4 S5 S6 S7,S1 S4 S7,S2 S5 S6,S5 S6,S1 S4 S7,S2,S3 S4 S5,S2,S3,S4 S5 S6 S7,S5 S6,S1 S4 S7,非多余断裂族:断裂组S1，S4，S7为最优断裂组。,通过断裂可以把不可分隔子系统中的回路物流打开，从而可以利用序贯模块法对该过程系统进行模拟计算。这种模拟计算的

11、开始是首先要设定起始物流变量的猜值，计算的终点则在于该猜值与计算值的收敛。,2.2.3 断裂物流变量的收敛,执行断裂物流变量收敛功能的模块称收敛单元模块 X 是断裂物流变量的猜值，y是经过程系统模型计算得出的断裂物流变量的计算值。断裂物流变量的收敛问题实际上是个迭代求解非线性方程组的问题x=y=G(x),当断裂物流变量猜值x与计算值y之差小于收敛容差时y x=G（x）x 则x为断裂物流变量的收敛解收敛单元的功能总计有如下作用：()获取猜值的初值x0()根据计算值y，以一定的方法确定新的猜值x()比较猜值x和计算值y，若其结果满足给定精度要求则结束迭代计算，否则继续迭代计算过程,可见，收敛单元实

12、质上就是一个数值迭代求解非线性方程组的子程序。求解非线性方程组的数值计算方法很多，适合于收敛单元的数值计算方法一般应尽可能满足下列要求：1.对初值的要求不高。1)初值易得，不易引起迭代计算的发散；2)初值的组数少。,2.数值稳定性好,2.数值稳定性好,2.数值稳定性好,3.收敛速度快对收敛速度的影响主要有三个因素：迭代次数；函数G（x）的计算次数；矩阵求逆的次数。每计算一次函数值就相当于做一次流程回路的模拟计算,每求一次导数就要做两次流程模拟计算对于断裂物流的收敛，好的非线性方程组的数值迭代次数少，而且应该尽量避免导数计算和矩阵求逆4.占用计算机存储空间少,直接迭代法,直接迭代法是将计算值yk

13、作为下一轮迭代的猜值xk+1而实施迭代计算xK+1=yK,即xk+1G（xk）非线性方程组x=y=G(x)的另外一种形式为F(x)xG(x)=0 xk+1xkF(xk)与牛顿公式相比较：,直接迭代法的雅可比矩阵为单位矩阵优点：方法简单，只需要一组初值，不需计算导数和逆矩阵弱点：迭代次数多、收敛速度慢，对初值要求较高,例:用直接迭代法求解下列方程组,解：令猜值为X12；X210；X35,解：令猜值为X16；X23.5；X35,为改善直接迭代法的收敛行为，提出了阻尼直接迭代法，或称加权直接迭代法：q为阻尼因子，可以人为给定q=0 直接迭代0q1 加权直接迭代，可改善收敛的稳定性q0 外推直接迭代，

14、加速收敛，但稳定性下降q1 无意义,当闪蒸温度分别为以下值时:分别用直接迭代法和阻尼直接迭代法计算汽相和液相产品的流量和组成，阻尼因子分别取值为0.5，0.3，-0.2，-0.3，-0.7，-0.9解：依据闪蒸条件,设该闪蒸过程为理想体系，三个闪蒸器均为等温闪蒸过程，建成相应的单元模块。并将其改绘为如下三级闪蒸过程模拟模块流程,平衡闪蒸单元模型,从上例可见，阻尼因子q值的选取具有较大的任意性和经验性。1958年Wegstien提出了一种简便的方法，可以弥补这种阻尼因子取值困难的弱点,Wegstein法,一维Wegstein法求解一维代数方程（1）Wegstein迭代公式为：其中：（3）,（2）

15、,对于隐式一维代数方程：（4）相应的迭代公式称作割线法，其迭代公式可从Wegstein迭代公式导出从（4）式可得出：（5）将上式代入（2）式,（6）,从（2）和（3）式得到：（7）上式代入（5）式，则有：（8）式（8）就是割线法的迭代公式。由此可见，Wegstein法与割线法是相通的,隐式方程具有更大的普遍性，所以割线法常为人们所熟知。在流程模拟领域中，物流回路多用显式方程描述的。因而多用Wegstein法。,由（2）式可见，一维Wegstein法需要有两个初值，其中第一个初值是设置的猜值，第二个初值可根据第一个初值按直接迭代法得到。,阻尼因子q的取值不当可使迭代计算收敛缓慢甚至发散Wegstein法虽然无须人为选定q值，但是也会因为q值不当导致坏的收敛行为有界Weigstein法就是凭借经验人为地把q值限定在一定的范围内，以改善收敛行为,有界Wegstein法,FLOWTRAN流程模拟系统中取为-5、0CHESS系统中当q0或q-10时令q=0,