化工过程分析与合成

化工过程分析与合成第一页，共六十三页，编辑于2023年，星期日FLOWTRAN流程模拟系统中取为-5、0CHESS系统中当q>0或q<-10时令q=0第二页，共六十三页，编辑于2023年，星期日多维Wegstein法分别用于每一个分量。令初始猜值为x0,则第二个初值可由直接迭代得到第三页，共六十三页，编辑于2023年，星期日第四页，共六十三页，编辑于2023年，星期日严格多维Wegstein法用向量代替变量，通过矩阵运算进行迭代求解对于n维方程，这一方法需要n+1组初始猜值（x0,x1,…,xn）第五页，共六十三页，编辑于2023年，星期日第六页，共六十三页，编辑于2023年，星期日(4)收敛判据断裂物流迭代计算的收敛判据通常是猜值与计算值的绝对误差或相对误差第七页，共六十三页，编辑于2023年，星期日2.2.4.序贯模块法解设计问题序贯模块法具有计算方向不可逆的特点，单元模块的计算只能按从输入到输出的方向进行只能通过调整某些决策变量或系统参数使计算结果满足设计要求D－设计规定向量H－过程系统方程组p－决策变量与系统参数向量第八页，共六十三页，编辑于2023年，星期日①估计反应单元的温度为T②估计再循环物流S4③依次计算混合单元、反应单元、分离单元，得到新的S4’的④比较S4与S4’，若两者相等则进行下一步，若不相等则返回②在收敛单元内比较S5和设计值，若两者不相等则返回①，若相等则计算结束控制模块图2-19具有再循环物流的过程系统混合器S1S2反应器分离器S3S4S5

第九页，共六十三页，编辑于2023年，星期日控制模块的设置增加了迭代循环圈，导致计算量的增加为了提高收敛速度可以联立求解再循环物流方程和设计方程，这就是同时收敛。使断裂物流变量x和系统参数p同时逼近收敛解，从而大大的提高了收敛速度第十页，共六十三页，编辑于2023年，星期日2.3过程系统模拟的面向方程法序贯模块法由于具有收敛计算的循环圈以致大大的增加了计算量。对于过程系统的设计计算问题和参数优化问题，情况将更为严重。因此，人们把注意力投向了面向方程法第十一页，共六十三页，编辑于2023年，星期日2.3.1面向方程法的原理把描述过程系统的所有数学模型汇集到一起，形成一个非线性方程组进行求解

x－状态变量向量w－决策变量向量

F－系统模型方程组，其中包括① 物性方程② 物料、能量、化学平衡方程③ 过程单元间的联结方程④ 设计规定方程等等第十二页，共六十三页，编辑于2023年，星期日比之序贯模块法，在决策变量的确定上要随意的多，决策变量和状态变量的地位是等同的通常可以把设计规定的变量（如系统出口浓度）直接指定为决策变量。面向方程法在求解一般模拟问题和设计问题上是没有差异的第十三页，共六十三页，编辑于2023年，星期日通常过程系统模型方程组总是稀疏方程组过程系统模型的方程数和变量数往往都很大，但每个方程涉及的变量数一般只有几个第十四页，共六十三页，编辑于2023年，星期日面向方程法的核心问题是求解超大型稀疏非线性方程组，求解方法大致分为两类① 降维求解法② 联立求解法第十五页，共六十三页，编辑于2023年，星期日2.3.2大型稀疏非线性方程组的降维解法把大型稀疏方程组分解成若干个小的非稀疏方程组，然后依次分别求解，从而达到降维和增大稀疏比的目的（1）方程组的分解概念对于n阶稀疏方程组，常常可以找到一个包含有k1个变量的k1阶子方程组。这个k1阶子方程组可以单独求解。其余的n-k1个方程中还可以再找出包含有k2个变量的k2阶子方程组，这个子方程组也可以单独求解。重复这一过程，最终将把原方程组分解成一系列可顺序求解的子方程组第十六页，共六十三页，编辑于2023年，星期日第十七页，共六十三页，编辑于2023年，星期日第十八页，共六十三页，编辑于2023年，星期日（2）回路搜索法分解方程组在描述方程组的有向图上进行回路搜索。为了用有向图表示方程组的结构，首先必须对每个方程指定一个变量作为其输出变量输出变量是可通过其所存在的方程中其它变量求解的变量，且每个变量只能被指定一次作为输出变量步骤：选事件矩阵中元素最少的行和元素最少的列的交点处元素对应的变量作为优先指定的输出变量，然后从事件矩阵中删去该输出变量对应的行和列重复上述过程直至矩阵中所有的行和列都被删掉第十九页，共六十三页，编辑于2023年，星期日第二十页，共六十三页，编辑于2023年，星期日有向图图中每个节点代表一个方程。如果方程fi的输出变量存在于fj中，则从节点fi向fj作一有向边这个图代表了方程间的信息流动方向4第二十一页，共六十三页，编辑于2023年，星期日回路搜索第二十二页，共六十三页，编辑于2023年，星期日第二十三页，共六十三页，编辑于2023年，星期日不可分解稀疏方程组的断裂降维解法1第二十四页，共六十三页，编辑于2023年，星期日断裂与收敛是相辅相成的，断裂后的系统必须通过收敛得以求解。为了易于收敛，因而总是希望断裂的变量数最少。所以，总是要选择包含变量数最少的方程中的变量作为断裂变量，断裂变量数等于该方程中的变量数减1。然后给断裂变量赋初值，再进行迭代计算直至收敛f3，f4，f5行的变量数最少，都只有两个。选择f3中的x5为断裂变量。从而解出x6第二十五页，共六十三页，编辑于2023年，星期日把f3行和x5,x6列删去，得到左式该式为五行四列，有一个多余方程（它是由删除断裂变量x5产生的）。对其余的四行，四列进行重排，可得到右式第二十六页，共六十三页，编辑于2023年，星期日2.3.3

联立拟线性方程组法解

大型稀疏非线性方程组大型稀疏非线性方程组的另一种求解方法是把非线性方程组线性化。然后联立求解线性方程组。由于线性化引入了误差，所以要借助迭代使线性化方程组的解逐渐逼近非线性方程组的解第二十七页，共六十三页，编辑于2023年，星期日线性化方法对于n维非线性方程组用n维线性方程组逼近该拟线性方程组的解（用下标QL表示）为：作台劳展开可得到牛顿迭代解（下标NR）第二十八页，共六十三页，编辑于2023年，星期日把（2－35）式代入（2－37）式，得到令J=A第二十九页，共六十三页，编辑于2023年，星期日牛顿迭代具有二阶收敛特性。下面方程也具有二阶收敛。系数A和B均是向量x的函数。从x的第k次近似解xk可以计算得到Jk、F(xk)，从而得到Ak和Bk。将Ak和Bk代入，得到线性方程组。过程系统的模型方程组一般由线性方程和非线性方程组成，因而线性化的对象应该是非线性方程j第三十页，共六十三页，编辑于2023年，星期日例

组分A的稀溶液在常温下离解：质量平衡热力学平衡求当k=2，A的初始浓度＝1时平衡态的组分浓度解：质量平衡式是线性方程，热力学平衡式是非线性方程，首先利用对热力学平衡式线性化第三十一页，共六十三页，编辑于2023年，星期日此外，还可以得到原方程的另一种线性化方程（即直接迭代式）两种方法都可以收敛到解。第一种方法的收敛速度明显比第二种方法快。这是由于牛顿迭代法具有二次收敛的特点，而直接迭代法只是线性收敛第三十二页，共六十三页，编辑于2023年，星期日第三十三页，共六十三页，编辑于2023年，星期日稀疏线性方程组的解法稀疏非线性方程组经线性化后得到的线性方程组仍然是稀疏的，从而把求解稀疏非线性方程组的问题转化成求解稀疏线性方程组的问题常规的消去法是不经济的，且计算效率低。为了减少计算时间和存储空间，常用下列两方面的技术只对非零元素进行计算只存储非零元素（如压缩存储技术）第三十四页，共六十三页，编辑于2023年，星期日填充量

用高斯消去法进行消元过程的同时，会在原来零元素处引入非零元素新出现的非零元素称作填充量，填充时与消元成零的非零元素之差称作填充增量。填充量与主元选取的次序有关第三十五页，共六十三页，编辑于2023年，星期日在求解大型稀疏线性方程组时，应尽可能减少填充，否则会使计算效率下降。减少填充与提高数值稳定性和计算精度是矛盾的。如，为减少填充，需把55作为主元素，但如果它的绝对值很小，会引入较大的误差，使计算精度、数值稳定性变差第三十六页，共六十三页，编辑于2023年，星期日主元容限通常把绝对值最大的元素作为主元，进行消元。目的是提高计算精度。但如果这样选取的主元导致较大的填充，将引起计算效率的下降往往选择一个绝对值不是最大，且不会引起填充量过大的元素作为主元人为规定一个界限e＞0。当矩阵元素的绝对值大于e，该元素就具备了作为主元的资格，若它引入的填充量也不是很大，就可定为主元。这个界限称为主元容限经验给定，但应满足提高计算精度和减少填充量的统一要求第三十七页，共六十三页，编辑于2023年，星期日Bending-Hutchison算法该算法是在全元消去法的基础上派生出来的一种求解稀疏线性方程组的算法。其核心是避免填充，同时保证计算的精度用过的：凡与被选作主元的元素有关的方程和变量都称作“用过的”，反之为“未用过的”；橫列（rank）：未用过的方程中包含的未用过的变量数；纵列（file）：未用过的变量在未用过的方程中出现的次数第三十八页，共六十三页，编辑于2023年，星期日1选择纵列最小的变量，如不止一个，任选其一；2在与此变量有关的方程中，选择橫列最小的方程所对应的元素作为主元3如果橫列最小的方程不止一个，则选择绝对值最大的元素作为主元4检验选出主元的绝对值是否大于用户给出的主元容限。不大于，则返回①.否则进行下一步；5用这样选择出的主元进行常规的高斯消元，然后返回①。上述过程中，步骤①和②都是为了避免填充。而步骤③和④是为了保证计算精度和系数矩阵非奇异。第三十九页，共六十三页，编辑于2023年，星期日

例

一个物流分割器及混合器构成的简化流程。第四十页，共六十三页，编辑于2023年，星期日列2和列8只含一个元素，即纵列=1。这两个元素分别为方程1和8的主元。这两列中无其它元素，不用执行消元过程。第3，5，7，9列均含两个非零元素，即纵列=2。选列3，非零元素存在于方程2和9中，方程2橫列=2，方程9橫列=3，选方程2中的该元素为主元。消去方程9中第3列的元素，这将导致方程9中的第一列产生一个非零元素。反复进行上述过程，然后进行回代过程第四十一页，共六十三页，编辑于2023年，星期日第四十二页，共六十三页，编辑于2023年，星期日V第四十三页，共六十三页，编辑于2023年，星期日2.4过程系统模拟的联立模块法两种系统模拟方法的比较联立模块法与序贯模块法的共同之处在于面向模块；与面向方程法共同在于联立求解过程过程系统模型方程第四十四页，共六十三页，编辑于2023年，星期日联立模块法利用严格模块产生相应的简化模型方程的系数，然后把所有的简化模型方程汇集到一起进行联解,得到系统的一组状态变量。由于简化模型是严格模块的近似，所以计算结果往往不是问题的解，必须用严格模块对这组解进行计算，修正简化模型的系数。重复这一过程，直到收敛到原问题的解第四十五页，共六十三页，编辑于2023年，星期日把序贯模块法中最费时、收敛最慢的回路迭代计算，用由简化模型组成的方程组的联解而代之，使计算加速，尤其是处理有多重再循环流或有设计规定要求的问题时具有较好的收敛行为。因此，联立模块法计算效率较高由于单元模块数比之过程方程数要少得多，所以简化模型方程组的维数比面向方程法也小得多，求解起来也容易得多。能利用大量原有的丰富的序贯模块软件。可在原有序贯模块模拟器上修改得到联立模块模拟器。特点：第四十六页，共六十三页，编辑于2023年，星期日① 计算效率较高；② 对初值要求较低；③ 迭代循环圈较少；④ 计算出错时诊断较容易；⑤ 能利用大量原有的软件。优点：第四十七页，共六十三页，编辑于2023年，星期日12n严格模块简化模型方程联立解{状态变量}图2-24联立模块法第四十八页，共六十三页，编辑于2023年，星期日S’

赋初值k=1开始扰动各模块入口变量，求简化模型系数Ak按流程结构组合AK建立系统简化模型求解简化模型得到ΔSik结束TF第四十九页，共六十三页，编辑于2023年，星期日以过程单元为基本单位建立简化模型；以回路为基本单位建立简化模型这两种划分策略分别与两种切断方式相对应联结物流全切断方式；回路切断方式。2.4.2建立简化模型的切断方式第五十页，共六十三页，编辑于2023年，星期日这种方式相当于把所有过程单元之间的联结物流全部切断，形成一系列互相独立的过程单元联结物流全切断方式yx图2－26联结物流全切断方式4123yxxyyx

第五十一页，共六十三页，编辑于2023年，星期日例

用联立模块法对三级闪蒸过程进行稳态模拟图2-18三级闪蒸过程的模拟模块流程入料混合器1闪蒸器3液相产品闪蒸器1闪蒸器2汽相产品收敛单元混全器2第五十二页，共六十三页，编辑于2023年，星期日解：①建立简化模型严格单元模块的输入流股变量向量x与输出流股变量y之间有严格模型：上式的一阶台劳展开式为即第五十三页，共六十三页，编辑于2023年，星期日令便可得到严格模型的线性增量简化模型别对每个过程单元写出其简化模型：混合器：闪蒸器1：闪蒸器2：闪蒸器3：第五十四页，共六十三页，编辑于2023年，星期日由于混合器的严格模型为线性模型，且系统入料流股变量为给定值，所以有把上述线性简化模型写成矩阵形式的迭代格式，则有：第五十五页，共六十三页，编辑于2023年，星期日②从严格模块计算简化模型的系数对每个单元建立简化模型，然后把单元简化模型、联结方程、设计规定方程集合到一起组成过程系统的简化模型，由于切断了全部联结物流，描述整个过程系统的简化模型方程数为：ne系统简化模型方程数；nc联结物流数nd设计规定方程数；ci联结物流组分数第五十六页，共六十三页，编辑于2023年，星期日流股全切断方式很类似于面向方程法。主要区别在于后者是严格模型方程，变量数也要大得多（包括单元内部变量）在处理实际问题时。联结物流全切断方式的维数往往还是很大的。因此人们又提出了回路切断方式。第五十七页，共六十三页，编辑于2023年，星期日回路切断方式相当于把若干个