《化工过程模拟与分析(第七章联立模块法)(精品).ppt》由会员分享，可在线阅读，更多相关《化工过程模拟与分析(第七章联立模块法)(精品).ppt（27页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第七章第七章 联立模块法（联立模块法（Simultaneous modular approach）主要内容主要内容1.联立模块法的基本思想2.Jacobi矩阵的近似计算3.联立模块法示例4.联立模块法与联立方程法的区别5.联立模块法系统分析的特点6.联立模块三层法7.1 联立模块法的基本思想联立模块法的基本思想一、联立方程法与序贯模块法的比较一、联立方程法与序贯模块法的比较 序贯模块法：更容易被理解和接受、迭代层次多导致收敛慢 联立方程法：迭代层次少、难以被理解和接受、占内存多二、如何互补联立方程法与序贯模块法？二、如何互补联立方程法与序贯模块法？继承已开发的模块 仅对不可再分块联立求解 降低
2、方程维数和求解难度三、联立模块法基本思路三、联立模块法基本思路混合器反应器分离器X1X2X3X4X5设体系中共有c个组分，Xi为描述流股的c+2维向量（或称外部变量外部变量外部变量外部变量，包括温度、压力和各组份摩尔流率）1、求解不可再分块的简化线性模型简化线性模型方程组如何获得上述线性方程组？如何获得上述线性方程组？1.分流分率模型（Rosen）2.非线性方程组的线性展开2、交替进行序贯模块迭代和联立线性方程组的求解。可改写为可改写为AX=b形式形式混合器反应器分离器X1X2X3X4X5联立模块法步骤（仅对于不可再分块）联立模块法步骤（仅对于不可再分块）1.赋切割变量初值；2.利用严格模块运
3、行序贯模块法一次，获得内部变量；3.对各模块方程和设计规定方程进行一阶展开，获得只含外部变量只含外部变量的联立线性方程组；4.求解联立线性方程组，获得切割变量的计算值，若收敛则退出，否则返回2继续。优点是：在切割变量值更新时减少了计算量优点是：在切割变量值更新时减少了计算量优点是：在切割变量值更新时减少了计算量优点是：在切割变量值更新时减少了计算量联立模块法联立模块法联立模块法联立模块法I I I I 全切断去连接形式全切断去连接形式全切断去连接形式全切断去连接形式外部变量外部变量 描述各单元输入输出流股的变量（X1X4）和设计规定方程所约束的设备参数（u）,模拟过程所涉及的其余变量为内部变量
4、内部变量。mixreactorseparatorsplitter对所有联立方程在初始值附近做一阶展开可得：以g1为例：将模块的雅可比矩阵分为2块整理上式，则可得线性方程组的一般形式：X2 X3 X4 X1 u 求解上述线性方程组之后，则可更新X1X4，u稀疏矩稀疏矩阵阵联立模块法联立模块法联立模块法联立模块法II II II II 环路切断形式环路切断形式环路切断形式环路切断形式mixreactorseparatorsplitter虚拟模块虚拟模块输入输出方程数和变量数最少，占用内存最少。方程数和变量数最少，占用内存最少。联立模块法联立模块法联立模块法联立模块法III III III III
5、非线性联立模块法非线性联立模块法非线性联立模块法非线性联立模块法 对各单元建立非线性简单模型，不可再分块可通过联立求解这些单元的联立非线性简单方程组求解。对简单非线性模型的要求对简单非线性模型的要求1.方程数比严格模型少得多；2.必须包含输入、输出变量，设备参数、热物性及其他；3.模型自由度必须与严格模型相同；4.应是机理模型；5.每种单元模型必须唯一确定，且可解析；6.简化模型参数由严格模型给出；优点优点1.强稀疏性；2.方程数目少；3.更容易求解。合适的求解方法：牛顿-拉弗逊法联立模块法的实质联立模块法的实质1.将包含非线性描述方程和相互影响的循环流的困难问题转变成一系列易于求解的线性描述
6、方程的简单子问题；2.单元内部变量和方程从线性系统消去，线性方程组只包含输入、输出信息流变量，因而方程求解难度和占用内存远小于联立方程法；3.相对于序贯模块法，求解设计型问题更有效。目前联立模块法的研究重点目前联立模块法的研究重点1.如何产生近似Jacobi矩阵；2.联立线性方程组涉及的变量；3.不可再分块切断准则；4.如何提出更多层的联立模块法；5.近似Jacobi矩阵的更新。五、五、I、II两种方式的比较两种方式的比较切断方式；方程组Jacobi矩阵稀疏性；与序贯模块法和联立方程法的比较；7.2 单元模块单元模块Jacobi矩阵的计算矩阵的计算一、解析法（一、解析法（Analytical
7、differentiation）最准确、计算效率最高，但难以实现。最准确、计算效率最高，但难以实现。unitu该单元数学模型为 H(Z,Y)=0，其中 Z=X；u，则其Jacobi矩阵为：二、全块扰动法（二、全块扰动法（Full-block perturbation）即用差商代替偏导数，如下所示：自变量为n维的非线性方程组，计算该矩阵须运算n+1次方程组单元模块的差分近似Jacobi矩阵：ei 单位矩阵的第i列；hi 第i个自变量的步长；fu 单元模型方程组；n 输入流股自变量数+设计方程约束的设备参数数目。该计算方法要求每计算一次Jacobi矩阵需要运行n+1次严格模型的计算。注注1.效果最
8、好；2.再循环流、变量多和单元需要迭代收敛导致计算量大；3.相对序贯模块法，该方法用于模拟型问题并不具有优势；4.该法用于设计型问题更有优势。三、对角块摄动法（三、对角块摄动法（Diagonal-block perturbation）计算效率高，误差大。计算效率高，误差大。假设第k个输出流股向量的取值只受第j个输入流股向量的影响，则：模块模拟次数：输入流股数+设计方程约束的设备参数数目四、直接差分近似法（四、直接差分近似法（Direct difference approximation）将组合单元看作一个虚拟模块，随后用全块扰动法求取该虚拟模块和设计规定方程的Jacobi矩阵。mixerrea
9、ctorseparatorsplitter虚拟模块虚拟模块输入输出随着切割变量数目增加，直接差分近似法计算随着切割变量数目增加，直接差分近似法计算Jacobi矩阵的工作量增加。矩阵的工作量增加。ABCD123456为提高虚拟模块Jacobi矩阵的计算效率，可根据欧拉链关系和虚拟模块中各模块的子Jacobi矩阵来计算虚拟模块的Jacobi矩阵。mixreactorseparatorsplitter虚拟模块虚拟模块输入输出Euler链关系（链关系（Euler chain relation）对于函数g=g(x,y)，其全微分为：若由f(x,y)=0可以推出y=(x)，并将之代入函数g全微分的表达式，
10、则：而因此根据欧拉链关系,(q)的Jacobi矩阵如下所示：如此，则虚拟模块的Jacobi矩阵计算量与切割变量数与切割变量数目无关目无关，且各模块内部变量可以方便地保存。7.3 联立模块法的结构分析联立模块法的结构分析对于联立模块法I，结构分析仅包含过程系统分隔以及分隔块之间计算次序；联立模块法II，不可再分块的切断以变量最少为基准；注注1.全模块扰动（模块摄动）法线性化效率最高，收敛快；2.对角块扰动法线性化效率高，收敛慢；3.仅对设计型问题，联立模块法优于序贯模块法；4.“组合单元”的全模块扰动法占用内存最少。闪蒸器1闪蒸器2闪蒸器3闪蒸器4混合器1混合器2P1S1R1Z1R2S2R3S3Z2P2F联联立立模模块块法法实实例例闪蒸器1闪蒸器2闪蒸器3闪蒸器4混合器1P1S1R1Z1R2S2R3S3Z2P2F13个组分个组分进料温度、压力进料温度、压力和组成已知和组成已知闪蒸器闪蒸器2，3，4温度压力已知温度压力已知闪蒸器闪蒸器1温度未知温度未知气相分率须满足设计值气相分率须满足设计值