第四章序贯模块法.pptVIP

第四章 序贯模块法 第四章 序贯模块法 本章要点 1）序贯模块法——直接选代，加权直接送代，严格Wegstein法，简化Wegstein法，Newton法，一维搜索。 2）序贯模块法的选代层次；系统的自由度；迭代格式；收敛准则，算法框图；算法的几何意义，线性收敛，平方收敛。 3）不同的加速收敛方法；各种收敛方法之间的比较。 4）序贯模块法解设计问题；自由变量（过程参数）；收敛方法；序贯模块法结构分析的切断准则。 5）数学方法应用于解模拟问题、设计问题。 过程系统的模拟分析 过程系统设计 过程系统参数优化 过程系统模拟的三种基本方法 序贯模块法（Sequential Modular Method） 联立方程法（Equation Oriented Method） 联立模块法（Simultaneously Modular Method） 过程系统模拟的序贯模块法 SMM的基本部分是模块（子程序），用以描述物性、单元操作以及系统其它功能。 SMM对过程系统的模拟以单元模块的模拟计算为基础。 SMM按照由各种单元模块组成的过程系统的结构，序贯的对各单元模块进行计算，从而完成该过程系统模拟计算。 过程系统模拟的联立方程法 EOM又称面向方程法，将描述整个过程系统的数学方程式联立求解，从而得出模拟计算结果 EOM解算快速有效，对设计、优化问题灵活方便。效率较高 EOM的形成通用软件比较困难；不能利用现有大量丰富的单元模块；缺乏实际流程的直观联系；计算失败之后难于诊断错误所在；对初值的要求比较苛刻；计算技术难度较大 过程系统模拟的联立模块法 过程系统稳态模拟三种方法的比较 1 序贯模块法迭代层次与系统自由度 经典序贯模块法的数学含义 由于化工流程中，虽然变量多，但大部分变量之间没有直接的函数关系，故可借助分隔的手段将系统分解成若干个必须联立求解的子方程组，然后对各个方程组分别求解。 涉及方程有： 1）物性估算方程：估算系统中各物质物性和混合物性质的模型方程； 2）单元模型方程：过程单元的模型方程； 3）流股连接方程：描述系统中各单元之间联结情况的方程； 4）设计规定方程：对单元输出、系统输出的约束方程。 迭代法 迭代法是数值计算中一类典型方法，应用于方程求根，方程组求解，矩阵求特征值等方面。其基本思想是逐次逼近，先取一个粗糙的近似值，然后用同一个递推公式，反复校正此初值，直至达到预定精度要求为止。 1 序贯模块法迭代层次与系统自由度 迭代层次 1 序贯模块法迭代层次与系统自由度 系统的自由度 系统自由度Dsys的确定是为了知道在系统模拟时应设定哪些必要的决策变量。 Di代表系统中单元 i的单元自由度；Kj为系统中单元之间第j个联结的联结限制数，即独立的联结流方程数。（每股流量的变量数:C+2） 2 不可再分块迭代的收敛方法 收敛：经过反复迭代，使迭代值接近初值的过程。 在切断处设置迭代收敛框，其作用是： 1）修正迭代变量；2）判别是否达到收敛。 迭代收敛：即当满足一定的收敛准则时模拟问题得到近似解。 修正切断迭代变量的常用方法： 1、直接迭代法；2、Wegstein法；3、Newton法；4、拟Newton法；5、综合算法。 2 不可再分块迭代的收敛方法 收敛准则 相邻两次迭代的迭代变量绝对误差的平方和小于某一预定的误差限ε1 相邻两次迭代的迭代变量相对误差的平方和小于误差限ε2 相邻两次迭代的迭代变量的加权平方和小于误差限ε3 2 不可再分块迭代的收敛方法 直接迭代法 直接迭代法是以切断流股的初始猜算值为起点，按X(i+1)=F(X(i))的迭代格式构成迭代序列进行迭代运算的。当满足迭代收敛准则式时，得到不可再分块的模拟解。 变量X的维数为切断流股的总变量数Nst。显然Nst大大小于整个不可再分块的总变量数。 序贯模块法将待解方程组进行了降阶处理，使求解过程变得更加容易。 2 不可再分块迭代的收敛方法 直接迭代法：收敛性与具体描述化工过程系统模型的非线性特征有关。 稳定单调下降收敛； 为振荡衰减收敛； 为振荡发散； 为直接发散。 2 不可再分块迭代的收敛方法 直接迭代法 一维x=f(x)直接迭代收敛的充分条件是f(x)一阶导数绝对值小于1。 多维是一阶偏导数矩阵绝对值小于1。 2 不可再分块迭代的收敛方法 加权直接迭代法 由于直接选代可能会发生不收敛的情况，为了改善收敛性能，采用加权直接送代的迭代格式：X(i+1)=QX(i)+(I-Q)F(X(i)) I——为单位矩阵，Q——为对角权矩阵。 2 不可再分块迭代的收敛方法 严格Wegstein法 严格Wegstein法是一种加速迭代收敛的方法，与加权直接迭代法的迭代格式在形式上完全相同