化工系统工程稳态模拟一3.pptVIP

* * * 4） 牛顿法 迭代公式： 特点：收敛速度很快，具有二次收敛性。 对维数很大的非线性方程组，直接利用比较困难。 基本思想：将非线性方程组逐次线性化。 和 离散型牛顿法 Broyden法 近似牛顿法： 2.3 过程系统模拟的联立方程法 基本思想： 将描述过程系统的所有变量（含单元内部的模型方程）全部联立起来，在数学上看作一个庞大的非线性方程组； ◆方程组的变量包括: 内部变量：物性变量、单元模型变量； 外部变量：描述系统结构的流股变量和过程参数； ◆建立方程组后，稳态模拟就可以视为一个与过程系统结构无关，甚至与化工背景无关的纯数学问题。 过程系统数学模型的特点 （1）高维大型非线性方程组 系统的模型方程组： X —— 状态变量； W —— 设计变量； F—— 系统模型方程组，其中包括： 物性 估算方程：计算混合物在各种状态下的物性。 单元模型方程：包括物料平衡方程、能量平衡方程、反应动力学方程、传递方程等。 过程单元间的联结方程：描述过程系统结构的拓扑关系。 设计规定方程。 （2）方程组的复杂性 过程系统数学模型所包含的方程的种类非常复杂，是一个极其复杂的混合型方程组： 过程系统由若干单元组成，整个系统的总变量数很大 现代化工物料、热量充分利用，热集成度高，组成了高度交联的流程拓扑结构 线性代数方程，如物料平衡方程 非线性代数方程，如化学平衡方程、热力学方程、物性估算方程； 微分方程，如化学反应动力学方程、描述反应器内部传递特征的偏微分方程 （3）方程组的稀疏性 无论描述全系统的方程组的维数有多高，函数关系有多复杂 每个过程单元所涉及的仅局限于几个方程，且每个方程往往仅涉及少数几个变量。 ? 化工系统的模型方程组是 超大型 稀疏矩阵。 描述方程组的稀疏程度－稀疏比：α＝N0/N2 N－描述系统的方程组的阶数 N0－为方程组线性化后的系数矩阵含非零元素的数目。 大型化工系统的数学模型，稀疏比一般小于1％， 约为0.2～0.5％。 （4）变量初值设定困难 目前求解非线性方程组的方法都是在开始迭代之前给出高维变量的初值（初始假设值），且这些变量的化工背景决定了变量在物理意义上有所限制，加上数学函数定义域的限制以及某些算法对初值的要求较高，使得初值合理给定比较困难。 利用非线性方程组的常用解法直接求解含有部分微分方程的联立方程组是有困难的，解决的方法有两种： （1）利用方程组分块技术，将含有微分方程的方程块，包括边界条件（代数形式的或微分形式的）分隔独立出来，用数值积分的方法求解； （2）将微分方程、包括导数形式的边界条件改写成差分方程，或直接从严格微分模型开发近似代数模型。（近似代数模型可以是线性或非线性的，在迭代过程中随时更新方程参数，以保持各迭代点附件近似模型与严格模型的等效性） 联立方程组中微分方程的处理 初值的选取 困难：－变量维数大 －迭代的收敛稳定性对初值要求苛刻 产生初值的方法： － 利用严格模型在基点附近产生一组简化模型，利用序贯模块法产生一组初值 － 产生一组近似线性模型，联立求解线性方程组得到初值 － 直接用严格模型按序贯模块法直接迭代数次后，得到初值 变量物理意义的限制 迭代计算中通过增补必须的等式约束或不等式约束对变量的取值加以适当限制，以避免变量在迭代过程中失去物理意义或超出计算的定义域引起的计算失败。 联立方程法的核心问题：求解超大型稀疏非线性方程组。 求解稀疏方程组的特殊的解算方法： 只对非零元素进行运算，以提高运算速度； 压缩贮存系数矩阵的非零元素，减少存储空间； …… …… 求解方法大致可分为两类： 降维求解法和线性联立求解法。 2.3.1 大型稀疏非线性方程组的降维解法 降维的必要性：所有方程联立求解困难。 降维的可能性：每一个方程并不含所有变量。 矩阵的稀疏性 降维的方法： （1）建立独立的物性数据估算模块； （2）取消单元间的联结方程； （3）进行方程组的分解。 2.3.1.1 建立独立的物性估算模块 原因： 化工系统数学模型中一半以上的方程为物性估算方程，约70-80%的模拟机时消耗在物性计算上。 办法： （1）建立单独物性估算模块，从系统模型方程组中独立出来。 （2）将焓值和平衡常数的计算建立简化线性模型，加入方程组的联立求解中，其余的物性估算仍然保留在物性计算模块中。 X1 X2 Y1 1 2 Y2 联结方程： X1 X2 1 2 取消联结方程后，模型方程组的阶数减少。