化工过程模拟与分析(联立方程法).pptVIP

* 第六章 联立方程法 主要内容 过程系统数学模型及其特点 过程系统非线性方程组及其求解的一些处理 大型稀疏非线性方程组的分解 联立方程法的潜在优势 6.1 过程系统数学模型及其特点 一、组成 物性估算方程 单元模型方程 连接流方程 设计规定方程 二、特点 高维大型非线性方程 方程组的复杂性 （线性、非线性、常微分、偏微分） 方程组的高稀疏性 变量初值给定的困难性 闪 蒸 器 1 闪 蒸 器 2 闪 蒸 器 3 闪 蒸 器 4 混合器1 混合器2 P1 S1 R1 Z1 R2 S2 R3 S3 Z2 P2 F 联立方程法示例 闪 蒸 器 1 闪 蒸 器 2 闪 蒸 器 3 闪 蒸 器 4 混合器1 P1 S1 R1 Z1 R2 S2 R3 S3 Z2 P2 F 13个组分 进料温度、压力 和组成已知 闪蒸器2，3，4 温度压力已知 闪蒸器1温度未知 气相分率须满足设计值 模拟步骤 1、根据需要，列出系统所包含单元种类，根据已知条件建立各单元模型的方程组； 混合器单元模型方程组 假设： 只需计算各流股组成 闪蒸器单元模型的基本方程组 2、建立联立方程中各单元形成的子方程组 各流股由一个13维向量描述，例如对于进料流股F： 流股摩尔流量 第1个组分摩尔组成 混合器1 混合器1 R1 Z1 R2 F 14个方程 闪蒸器1 P1 S1 闪 蒸 器 1 R1 29个方程 设计规定方程 闪蒸器2 28个方程 S1 闪 蒸 器 2 S2 Z1 闪蒸器3 R2 闪 蒸 器 3 S3 Z2 28个方程 闪蒸器4 R3 闪 蒸 器 4 P2 S3 28个方程 混合器2 混合器2 R3 Z2 S2 已知22个变量： 未知141个变量： 14个方程 方程数：141个 3、联立各单元形成的子方程组，分析特点并确定初值和选择解法 注 建立方程组后，必须检查方程数和变量数是否相等； 平衡常数的计算可独立在方程组之外； 各流股组成变量可能在求解过程中失去物理意义； 显然流股连结方程并不包含在联立方程组中； 方程组的雅可比矩阵显然稀疏； 显然141个初值确定比较困难。 6.2 过程系统数学模型的求解 一、微分方程的处理 将微分方程分隔出来数值积分 将微分方程改写成差分方程 直接从严格微分模型开发近似代数模型 二、初值的选取 利用严格模型生成简化模型，并用序贯模块法计算 通过求解线性化产生近似线性方程组并求解得到 按照序贯模块法迭代数次后得到初始解 三、求解中变量取值的限制 四、方程组的降阶 建立独立的物性估算模块 取消连接流方程 方程组的分隔和切断 五、稀疏矩阵的处理 压缩存储 稀疏线性方程组的消去法 6.3 大型稀疏非线性方程组的分解 概述 对象：数学模型的结构，而不是过程系统的拓扑结构。 分隔：确定必须联立求解的不可再分子方程块。 切断：切断的对象是不可再分子方程块中的变量。 方程组与其关联矩阵和邻接矩阵 关联矩阵 f x 邻接矩阵 从 f 到 f 邻接矩阵 在选出输出变量 集之后才可得到 联立方程的分解步骤 选取输出变量集； 构造方程组邻接矩阵； 利用回路搜索法等方法对方程组分块。 联立方程的切割原理 写出方程组的关联矩阵，调换该矩阵行和列的顺序，使出现在矩阵对角线上的元素为非0元素，然后选取矩阵上三角部分中含有非0元素的列所对应的变量，即切割变量。 有时不分块，直接切割再求解更为方便。 *