04序贯模块法课程.pptVIP

第四章 序贯模块法（Sequential modular method）;4.1 系统自由度的确定;混合器;用环路切断的方法来解算不可再分块，以切断流股上 的变量作为切断变量，通过迭代计算得到不可再分块 的解。;迭代收敛框的作用：修正迭代变量、判断是否达到收敛，要求满足收敛准则。;序贯模块法步骤 1、过程系统分块和计算顺序确定； 2、不可再分块的切断； 3、序贯求解各单元（或不可再分块）;4.3 不可再分块计算示例;模拟步骤;function [S2] = mix(S1,S5,S6) %混合器模型 % S1, S5, S6 – 进口物流 % S2 – 出口物流 S2.P = min([S1.P, S5.P, S6.P]); % 出口压力为最小压力 S2.C = (S1.C*S1.F + S5.C*S5.F + S6.C*S6.F) / (S1.F + S5.F + S6.F); % 出口物料组成 S2.T = (S1.T*S1.F + S5.T*S5.F + S6.T*S6.F) / (S1.F + S5.F + S6.F); % 出口物料温度 S2.F = S1.F + S5.F + S6.F; % 出口物料流量;闪蒸模型（等温平衡）;function [Fl, Cl, Tl, Fg, Cg, Tg] = flash(F,zs,ks,Tf,Pf,b0) % F - 进口物流摩尔流率 (mol/h) % zs - 进口物流组成 % ks - 进口物流各组分气液平衡常数 % Tf - 闪蒸器温度 (C) % Pf - 闪蒸器压力 (kPa) % b0 - 气相分率初始值 for i = 1:100 [f,df] = flash(b10, S20.C, flashKs(:,1)); b1 = b10 - f/df; if abs(b1-b10)1e-3 break; end b0 = b; end % 液相出料信息 Fl = F*(1-b); Tl = Tf; Pl = Pf; Cl = zs ./ ( 1+b*(ks-1) ); Cl = Cl/sum(Cl); % 气相出料信息 Fg = F*b; Tg = Tf; Pg = Pf; Cg = ks.*zs; Cg = Cg/sum(Cg); function [f,df] = flashEq(b, zs, ks) %闪蒸器方程 tmp = zs.*(1-ks) ./ ( 1+b*(ks-1) ); f = sum(tmp); tmp = ( zs.*(1-ks).^2 )./( ( 1+(ks-1)*b ).^2 ); df = sum(tmp);;2、根据模拟要求，列出所需物性数据种类和计算方法；;根据Antoine方程，可计算各组分在不同温度下饱和蒸汽压，如下表所示：;3、不可再分块的切割（无多余切断）;4、切断物流的初始值确定;5、列出系统模型的显式非线性方程组 x =φ (x);混合器;6、选定算法，求解上述方程组。;ASPEN模拟步骤;1、添加模块与物流;;未命名物流;对物流命名;全局设置;选择被切断物流;设置收敛标准;模拟计算方法的选择;添加过程涉及的纯组分;查询涉及到的组分;组分添加完毕;选择物性方法;添加输入物流信息;系统输入物流添加完毕;添加各模块操作条件;运行该流程;最终模拟结果;收敛信息;4.4 简化的加界Wegstein法;二、显式方程的求解-Wegstein法;x;Wegstein法即自适应确定松弛因子的部分迭代法； 几何上，Wegstein法的新点是迭代函数的割线与直线 y = x 的交点； 显然割线斜率sk不能接近1； 为保持迭代过程的稳定性，有人开发了加界Wegstein法和间歇Wegstein法（Delayed Wegstein method）;;三、显式方程组求解的Wegstein法;% 加界Wegstein法 X0 = [0.3; 0.3]; % 初始点 X1 = f_42(0,X0,0); for k = 1:1000 X2 = [0;0]; % 确立X2为列向量 for i = 1:2 xi = X1(i); x = f_42(xi,X1,i); ki = (x - X1(i)) / (X1(i) - X0(i)); oi = 1/(1-ki); if oi10 % 松弛因子加界 oi = 10; elseif oi-10 oi = -10; end X2(i) = oi*