化工过程分析与合成3.ppt
- 文档编号:2644980
- 上传时间:2022-11-05
- 格式:PPT
- 页数:64
- 大小:447KB
化工过程分析与合成3.ppt
《化工过程分析与合成3.ppt》由会员分享,可在线阅读,更多相关《化工过程分析与合成3.ppt(64页珍藏版)》请在冰豆网上搜索。
有界有界WegsteinWegstein法法阻尼因子q的取值不当可使迭代计算收敛缓慢甚至发散Wegstein法虽然无须人为选定q值,但是也会因为q值不当导致坏的收敛行为有界Weigstein法就是凭借经验人为地把q值限定在一定的范围内,以改善收敛行为FLOWTRAN流程模拟系统中取为-5、0CHESS系统中当q0或q-10时令q=0多维多维WegsteinWegstein法法分别用于每一个分量。
令初始猜值为x0,则第二个初值可由直接迭代得到严格多维严格多维WegsteinWegstein法法用向量代替变量,通过矩阵运算进行迭代求解对于n维方程,这一方法需要n+1组初始猜值(x0,x1,xn)(4)4)收敛判据收敛判据断断断断裂裂裂裂物物物物流流流流迭迭迭迭代代代代计计计计算算算算的的的的收收收收敛敛敛敛判判判判据据据据通通通通常常常常是是是是猜猜猜猜值值值值与与与与计计计计算算算算值值值值的的的的绝绝绝绝对对对对误误误误差差差差或相对误差或相对误差或相对误差或相对误差2.2.4.2.2.4.序贯模块法解设计问题序贯模块法解设计问题序序贯贯模模块块法法具具有有计计算算方方向向不不可可逆逆的的特特点点,单单元元模模块块的的计算只能按从输入到输出的方向进行计算只能按从输入到输出的方向进行只能通过调整某些决策变量或系统参数使计算结果满足设计要求D设计规定向量H过程系统方程组p决策变量与系统参数向量估计反应单元的温度为估计反应单元的温度为TT估计再循环物流估计再循环物流S4S4依次计算混合单元、反应单元、分离单元,得到新的依次计算混合单元、反应单元、分离单元,得到新的S4S4的的比较比较S4S4与与S4S4,若两者相等则进行下一步,若不相等则返回若两者相等则进行下一步,若不相等则返回在在收收敛敛单单元元内内比比较较S5S5和和设设计计值值,若若两两者者不不相相等等则则返返回回,若若相等则计算结束相等则计算结束控制模块的设置增加了迭代循环圈,导致计算量的增加为了提高收敛速度可以联立求解再循环物流方程和设计方程,这就是同时收敛。
使断裂物流变量x和系统参数p同时逼近收敛解,从而大大的提高了收敛速度2.32.3过程系统模拟的面向方程法过程系统模拟的面向方程法序贯模块法由于具有收敛计算的循环圈以致大大的增加了计算量。
对于过程系统的设计计算问题和参数优化问题,情况将更为严重。
因此,人们把注意力投向了面向方程法2.3.12.3.1面向方程法的原理面向方程法的原理把描述过程系统的所有数学模型汇集到一起,形成一个非线性方程组进行求解x状态变量向量w决策变量向量F系统模型方程组,其中包括物性方程物料、能量、化学平衡方程过程单元间的联结方程设计规定方程等等比之序贯模块法,在决策变量的确定上要随意的多,决策变量和状态变量的地位是等同的通常可以把设计规定的变量(如系统出口浓度)直接指定为决策变量。
面向方程法在求解一般模拟问题和设计问题上是没有差异的通常过程系统模型方程组总是稀疏方程组过程系统模型的方程数和变量数往往都很大,但每个方程涉及的变量数一般只有几个面向方程法的核心问题是求解超大型稀疏非线性方程组,求解方法大致分为两类降维求解法联立求解法2.3.22.3.2大型稀疏非线性方程组的降维解法大型稀疏非线性方程组的降维解法把大型稀疏方程组分解成若干个小的非稀疏方程组,然后依次分别求解,从而达到降维和增大稀疏比的目的(11)方程组的分解概念方程组的分解概念对于nnnn阶阶阶阶稀疏方程组,常常可以找到一个包含有kkkk1111个变量的kkkk1111阶阶阶阶子方程组。
这个kkkk1111阶阶阶阶子方程组可以单独求解。
其余的n-kn-kn-kn-k1111个方程中还可以再找出包含有kkkk2222个变量的kkkk2222阶阶阶阶子方程组,这个子方程组也可以单独求解。
重复这一过程,最终将把原方程组分解成一系列可顺序求解的子方程组
(2)回路搜索法分解方程组在描述方程组的有向图上进行回路搜索。
为了用有向图表示方程组的结构,首先必须对每个方程指定一个变量作为其输出变量输出变量是可通过其所存在的方程中其它变量求解的变量,且每个变量只能被指定一次作为输出变量步骤:
步骤:
步骤:
步骤:
选事件矩阵中元素最少的行和元素最少的列的交点处元素对应的变量作为优先指定的输出变量,然后从事件矩阵中删去该输出变量对应的行和列重复上述过程直至矩阵中所有的行和列都被删掉有向图图中每个节点代表一个方程。
如果方程fi的输出变量存在于fj中,则从节点fi向fj作一有向边这个图代表了方程间的信息流动方向4回路搜索不可分解稀疏方程组的断裂降维解法断裂与收敛是相辅相成的,断裂后的系统必须通过收敛得以求解。
为了易于收敛,因而总是希望断裂的变量数最少。
所以,总是要选择包含变量数最少的方程中的变量作为断裂变量,断裂变量数等于该方程中的变量数减1。
然后给断裂变量赋初值,再进行迭代计算直至收敛f3,f4,f5行的变量数最少,都只有两个。
选择f3中的x5为断裂变量。
从而解出x6把f3行和x5,x6列删去,得到左式该式为五行四列,有一个多余方程(它是由删除断裂变量x5产生的)。
对其余的四行,四列进行重排,可得到右式2.3.3联立拟线性方程组法解联立拟线性方程组法解大型稀疏非线性方程组大型稀疏非线性方程组大型稀疏非线性方程组的另一种求解方法是把非线性方程组线性化。
然后联立求解线性方程组。
由于线性化引入了误差,所以要借助迭代使线性化方程组的解逐渐逼近非线性方程组的解线性化方法对于n维非线性方程组用n维线性方程组逼近该拟线性方程组的解(用下标QL表示)为:
作台劳展开可得到牛顿迭代解(下标NR)把(235)式代入(237)式,得到令J=A牛顿迭代具有二阶收敛特性。
下面方程也具有二阶收敛。
系数A和B均是向量x的函数。
从x的第k次近似解xk可以计算得到Jk、F(xk),从而得到Ak和Bk。
将Ak和Bk代入,得到线性方程组。
过程系统的模型方程组一般由线性方程和非线性方程组成,因而线性化的对象应该是非线性方程例例例例组分A的稀溶液在常温下离解:
质量平衡热力学平衡求当k=2,A的初始浓度1时平衡态的组分浓度解解解解:
质量平衡式是线性方程,热力学平衡式是非线性方程,首先利用对热力学平衡式线性化此外,还可以得到原方程的另一种线性化方程(即直接迭代式)两种方法都可以收敛到解。
第一种方法的收敛速度明显比第二种方法快。
这是由于牛顿迭代法具有二次收敛的特点,而直接迭代法只是线性收敛稀疏线性方程组的解法稀疏非线性方程组经线性化后得到的线性方程组仍然是稀疏的,从而把求解稀疏非线性方程组的问题转化成求解稀疏线性方程组的问题常规的消去法是不经济的,且计算效率低。
为了减少计算时间和存储空间,常用下列两方面的技术只对非零元素进行计算只对非零元素进行计算只存储非零元素(如压缩存储技术)只存储非零元素(如压缩存储技术)填充量填充量填充量填充量用高斯消去法进行消元过程的同时,会在原来零元素处引入非零元素新新出出现现的的非非零零元元素素称称作作填填充充量量,填填充充时时与与消消元元成成零零的的非非零零元元素素之之差差称称作作填填充充增增量量。
填填充充量量与与主主元元选取的次序有关选取的次序有关在求解大型稀疏线性方程组时,应尽可能减少填充,否则会使计算效率下降。
减少填充与提高数值稳定性和计算精度是矛盾的。
如,为减少填充,需把55作为主元素,但如果它的绝对值很小,会引入较大的误差,使计算精度、数值稳定性变差主元容限主元容限主元容限主元容限通常把绝对值最大的元素作为主元,进行消元。
目的是提高计算精度。
但如果这样选取的主元导致较大的填充,将引起计算效率的下降往往选择一个绝对值不是最大,且不会引起填充量过大的元素作为主元人为规定一个界限e0。
当矩阵元素的绝对值大于e,该元素就具备了作为主元的资格,若它引入的填充量也不是很大,就可定为主元。
这这这这个个个个界界界界限限限限称称称称为为为为主主主主元元元元容限容限容限容限经验给定,但应满足提高计算精度和减少填充量的统一要求Bending-HutchisonBending-Hutchison算法算法算法算法该算法是在全元消去法的基础上派生出来的一种求解稀疏线性方程组的算法。
其其其其核核核核心心心心是是是是避避避避免免免免填填填填充充充充,同同同同时时时时保保保保证计算的精度证计算的精度证计算的精度证计算的精度用用用用过过过过的的的的:
凡与被选作主元的元素有关的方程和变量都称作“用过的”,反之为“未用过的”;橫列(橫列(橫列(橫列(rankrank):
):
):
):
未用过的方程中包含的未用过的变量数;纵纵纵纵列列列列(filefile):
未用过的变量在未用过的方程中出现的次数11选择纵列最小的变量,如不止一个,任选其一;选择纵列最小的变量,如不止一个,任选其一;22在与此变量有关的方程中,选择橫列最小的方程在与此变量有关的方程中,选择橫列最小的方程所对应的元素作为主元所对应的元素作为主元33如果橫列最小的方程不止一个,则选择绝对值最如果橫列最小的方程不止一个,则选择绝对值最大的元素作为主元大的元素作为主元44检验选出主元的绝对值是否大于用户给出的主元检验选出主元的绝对值是否大于用户给出的主元容限。
不大于,则返回容限。
不大于,则返回.否则进行下一步;否则进行下一步;55用这样选择出的主元进行常规的高斯消元,然后用这样选择出的主元进行常规的高斯消元,然后返回返回。
上述过程中,步骤上述过程中,步骤上述过程中,步骤上述过程中,步骤和和和和都是为了避免填充。
而步都是为了避免填充。
而步都是为了避免填充。
而步都是为了避免填充。
而步骤骤骤骤和和和和是为了保证计算精度和系数矩阵非奇异。
是为了保证计算精度和系数矩阵非奇异。
是为了保证计算精度和系数矩阵非奇异。
是为了保证计算精度和系数矩阵非奇异。
例例一个物流分割器及混合器构成的简化流程。
列2和列8只含一个元素,即纵列=1。
这两个元素分别为方程1和8的主元。
这两列中无其它元素,不用执行消元过程。
第3,5,7,9列均含两个非零元素,即纵列=2。
选列3,非零元素存在于方程2和9中,方程2橫列=2,方程9橫列=3,选方程2中的该元素为主元。
消去方程9中第3列的元素,这将导致方程9中的第一列产生一个非零元素。
反复进行上述过程,然后进行回代过程2.4过程系统模拟的联立模块法过程系统模拟的联立模块法两种系统模拟方法的比较两种系统模拟方法的比较联立模块法与序贯模块法的共同之处在于面向模块;联立模块法与序贯模块法的共同之处在于面向模块;与面向方程法共同在于联立求解过程过程系统模型与面向方程法共同在于联立求解过程过程系统模型方程方程联立模块法利用严格模块产生相应的简化模型方程的系数,然后把所有的简化模型方程汇集到一起进行联解,得到系统的一组状态变量。
由于简化模型是严格模块的近似,所以计算结果往往不是问题的解,必须用严格模块对这组解进行计算,修正简化模型的系数。
重复这一过程,直到收敛到原问题的解把序贯模块法中最费时、收敛最慢的回路迭代计算,用由简化模型组成的方程组的联解而代之,使计算加速,尤其是处理有多重再循环流或有设计规定要求的问题时具有较好的收敛行为。
因此,联联立模块法计算效率较高立模块法计算效率较高由于单元模块数比之过程方程数要少得多,所以简化模型方程组的维数比面向方程法也小得多,求解起来也容易得多。
能利用大量原有的丰富的序贯模块软件。
可在原有序贯模块模拟器上修改得到联立模块模拟器。
特点:
特点:
计算效率较高;对初值要求较低;迭代循环圈较少;计算出错时诊断较容易;能利用大量原有的软件。
优点:
优点:
12n严格模块简化模型方程联立解状态变量图2-24联立模块法S赋初值k=1开始扰动各模块入口变量,求简化模型系数Ak按流程结构组合AK建立系统简化模型求解简化模型得到Sik结束TF以过程单元为基本单位建立简化模型;以回路为基本单位建立简化模型这两种
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 化工 过程 分析 合成