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

合并得： 用矩阵表示为： 过程系统稳态模拟三种方法的比较 线性化方法 对于n维非线性方程组 用n维线性方程组逼近 过程系统的模型由线性方程组和非线性方程组组成， 线性化的对象是非线性方程组。 例 组分A的稀溶液在常温下离解： 质量平衡 热力学平衡 求当k=2，A的初始浓度＝1时平衡态的组分浓度 解：质量平衡式是线性方程，热力学平衡式是非线性方程，首先利用对热力学平衡式线性化 此外，还可以得到原方程的另一种线性化方程（即直接迭代式） 两种方法都可以收敛到解。第一种方法的收敛速度明显比第二种方法快。这是由于牛顿迭代法具有二次收敛的特点，而直接迭代法只是线性收敛 稀疏线性方程组的解法 稀疏非线性方程组经线性化后得到的线性方程组仍然是稀疏的，从而把求解稀疏非线性方程组的问题转化成求解稀疏线性方程组的问题 常规的消去法是不经济的，且计算效率低。为了减少计算时间和存储空间，常用下列两方面的技术 只对非零元素进行计算 只存储非零元素（如压缩存储技术） 填充量 用高斯消去法进行消元过程的同时，会在原来零元素处引入非零元素 新出现的非零元素称作填充量，填充时与消元成零的非零元素之差称作填充增量。填充量与主元选取的次序有关 在求解大型稀疏线性方程组时，应尽可能减少填充，否则会使计算效率下降。 减少填充与提高数值稳定性和计算精度是矛盾的。如，为减少填充，需把55作为主元素，但如果它的绝对值很小，会引入较大的误差，使计算精度、数值稳定性变差 主元容限 通常把绝对值最大的元素作为主元，进行消元。目的是提高计算精度。但如果这样选取的主元导致较大的填充，将引起计算效率的下降 往往选择一个绝对值不是最大，且不会引起填充量过大的元素作为主元 人为规定一个界限 e＞0。当矩阵元素的绝对值大于e，该元素就具备了作为主元的资格，若它引入的填充量也不是很大，就可定为主元。这个界限称为主元容限 经验给定，但应满足提高计算精度和减少填充量的统一要求 Bending-Hutchison算法 该算法是在全元消去法的基础上派生出来的一种求解稀疏线性方程组的算法。其核心是避免填充，同时保证计算的精度 用过的：凡与被选作主元的元素有关的方程和变量都称作“用过的”，反之为“未用过的”； 橫列（rank）：未用过的方程中包含的未用过的变量数； 纵列（file）：未用过的变量在未用过的方程中出现的次数 1 选择纵列最小的变量，如不止一个，任选其一； 2 在与此变量有关的方程中，选择橫列最小的方程所对应的元素作为主元 3 如果橫列最小的方程不止一个，则选择绝对值最大的元素作为主元 4 检验选出主元的绝对值是否大于用户给出的主元容限。不大于，则返回①.否则进行下一步； 5 用这样选择出的主元进行常规的高斯消元，然后返回①。 上述过程中，步骤①和②都是为了避免填充。而步骤③和④是为了保证计算精度和系数矩阵非奇异。 例 一个物流分割器及混合器构成的简化流程。 列2和列8只含一个元素，即纵列=1。这两个元素分别为方程1和8的主元。这两列中无其它元素，不用执行消元过程。 第3，5，7，9列均含两个非零元素，即纵列=2。选列3，非零元素存在于方程2和9中，方程2橫列=2，方程9橫列=3，选方程2中的该元素为主元。 消去方程9中第3列的元素，这将导致方程9中的第一列产生一个非零元素。 反复进行上述过程，然后进行回代过程 2.4 过程系统模拟的联立模块法 两种系统模拟方法的比较 联立模块法与序贯模块法的共同之处在于面向模块； 与面向方程法共同在于联立求解过程过程系统模型方程 联立模块法利用严格模块产生相应的简化模型方程的系数，然后把所有的简化模型方程汇集到一起进行联解,得到系统的一组状态变量。 由于简化模型是严格模块的近似，所以计算结果往往不是问题的解，必须用严格模块对这组解进行计算，修正简化模型的系数。 重复这一过程，直到收敛到原问题的解 把序贯模块法中最费时、收敛最慢的回路迭代计算，用由简化模型组成的方程组的联解而代之，使计算加速，尤其是处理有多重再循环流或有设计规定要求的问题时具有较好的收敛行为。因此，联立模块法计算效率较高 由于单元模块数比之过程方程数要少得多，所以简化模型方程组的维数比面向方程法也小得多，求解起来也容易得多。 能利用大量原有的丰富的序贯模块软件。可在原有序贯模块模拟器上修改得到联立模块模拟器。 特点： ① 计算效率较高； ② 对初值要求较低； ③ 迭代循环圈较少； ④ 计算出错时诊断较容易； ⑤ 能利用大量原有的软件。 优点： 1 2 n 严格模块 简 化 模 型 方 程 联 立 解 { 状 态 变 量 } 图2-24 联 立 模 块 法 S’ 赋初值 k =1 开 始 扰