化工系统工程_01_2系统模拟方法.ppt

系统模拟方法概述 模拟是对过程系统模型的求解 过程系统的模拟分析 对某个给定的过程系统模型进行模拟求解，可得出该系统的全部状态变量，从而可以对该过程系统进行工况分析 过程系统设计当对某个或某些系统变量提出设计规定要求时，通过调整某些决策变量使模拟结果满足设计规定要求 过程系统参数优化过程系统模型与最优化模型联解得到一组使工况目标函数最佳的决策变量（优化变量）。从而实施最佳工况 过程系统模拟的三种基本方法 序贯模块法（SMM, Sequential Modular Method） 面向方程法（EOM, Equation Oriented Method） 联立模块法（Simultaneously Modular Method） 过程系统模拟的序贯模块法 SMM的基本部分是模块（子程序），用以描述物性、单元操作以及系统其它功能。 SMM对过程系统的模拟以单元模块的模拟计算为基础。 SMM按照由各种单元模块组成的过程系统的结构，序贯的对各单元模块进行计算，从而完成该过程系统模拟计算。 过程系统模拟的面向方程法 EOM又称联立方程法，将描述整个过程系统的数学方程式联立求解，从而得出模拟计算结果 EOM解算快速有效，对设计、优化问题灵活方便。效率较高 EOM的形成通用软件比较困难；不能利用现有大量丰富的单元模块；缺乏实际流程的直观联系；计算失败之后难于诊断错误所在；对初值的要求比较苛刻；计算技术难度较大 过程系统稳态模拟三种方法的比较 通过断裂技术可以打开回路，以便采用序贯模块法进行求解。在断裂物流处设置一个收敛单元。 对于复杂系统，收敛单元设置的位置不同，其效果也将不同。 最优设置要通过断裂技术去解决。 如何得到新的变量值，如何保证计算收敛，如何加快收敛，取决于收敛算法，与断裂物流变量的特性有关。 由于系统中各物流及其变量特性的不同，在收敛计算上常是有很大差异的。 如何选择断裂物流、确定迭代序列，是实施序贯模块法进行过程系统模拟计算中必须要解决的问题。 例: 用直接迭代法求解下列方程组 Wegstein法的几何意义 考虑到连续函数都可进行Taylor展开，也即在足够小区域内用有限项多项式（特别是一次多项式）近似原函数可获得足够精确性。 Wegstein法的几何意义 如采用光滑曲线上某点的切线近似则导出牛顿（Newton）法。 如采用光滑曲线上两点连线近似则导出割线法。 现考虑一次方程一般情形： Wegstein法的几何意义 方程的根计算的是函数与y = 0直线的交点，即下面方程组的解： Wegstein法的几何意义 x = kx + b属于不动点迭代标准形式。此情况下方程的根为 FLOWTRAN流程模拟系统中取为-5、0 CHESS系统中当q＞0或q＜-10时令q=0 断裂与收敛是相辅相成的，断裂后的系统必须通过收敛得以求解。 为了易于收敛，因而总是希望断裂的变量数最少。 所以，总是要选择包含变量数最少的方程中的变量作为断裂变量，断裂变量数等于该方程中的变量数减1。然后给断裂变量赋初值，再进行迭代计算直至收敛 f3，f4，f5行的变量数最少，都只有两个。选择f3中的x5为断裂变量。从而解出x6 令J = A 牛顿迭代具有二阶收敛特性。下面方程也具有二阶收敛。 系数A和B均是向量x的函数。 从x的第k次近似解xk可以计算得到Jk、F(xk)，从而得到Ak和Bk。将Ak和Bk代入，得到线性方程组。 过程系统的模型方程组一般由线性方程和非线性方程组成，因而线性化的对象应该是非线性方程 例 组分A的稀溶液在常温下离解： 质量平衡 热力学平衡 求当k=2，A的初始浓度＝1时平衡态的组分浓度 解：质量平衡式是线性方程，热力学平衡式是非线性方程，首先利用对热力学平衡式线性化 对于隐式一维代数方程： （4） 相应的迭代公式称作割线法，其迭代公式可从Wegstein迭代公式导出 从（4）式可得出： （5） 将上式代入（2）式 （6） 从（2）和（3）式得到： （7） 上式代入（5）式，则有： （8） 式（8）就是割线法的迭代公式。由此可见，Wegstein法与割线法是相通的 隐式方程具有更大的普遍性，所以割线法常为人们所熟知。 在流程模拟领域中，物流回路多用显式方程描述的。因而多用Wegstein法。 由（2）式可见，一维Wegstein法需要有两个初值，其中第一个初值是设置的猜值，第二个初值可根据第一个初值按直接迭代法得到。 有界Wegstein法 阻尼因子q的取值不当可使迭代计算收敛缓慢甚至发散； Wegstein法虽然无须人为选定q值，但是也会因为q值不当导致坏的收敛行为； 有界Weigstein