《化工过程分析与合成2》由会员分享,可在线阅读,更多相关《化工过程分析与合成2(63页珍藏版)》请在金锄头文库上搜索。
1、模拟是对过程系统模型的求解,过程系统的模拟分析 对某个给定的过程系统模型进行模拟求解,可得出该系统的全部状态变量,从而可以对该过程系统进行工况分析,过程系统设计当对某个或某些系统变量提出设计规定要求时,通过调整某些决策变量使模拟结果满足设计规定要求,过程系统参数优化过程系统模型与最优化模型联解得到一组使工况目标函数最佳的决策变量(优化变量)。从而实施最佳工况,思考题,结合上节课发酵过程的例子,说明过程系统的模拟、设计和优化。 改变上节课例子中的参数,分别进行过程系统模拟、设计和优化。,2.1 过程系统模拟的三种基本方法,序贯模块法(Sequential Modular Method) 面向方程
2、法(Equation Oriented Method) 联立模块法(Simultaneously Modular Method),过程系统模拟的序贯模块法,SMM的基本部分是模块(子程序),用以描述物性、单元操作以及系统其它功能。 SMM对过程系统的模拟以单元模块的模拟计算为基础。 SMM按照由各种单元模块组成的过程系统的结构,序贯的对各单元模块进行计算,从而完成该过程系统模拟计算。,优点: 与实际过程的直观联系强 模拟系统软件的建立、维护和扩充都很方便,易于通用化 计算出错时易于诊断出错位置,缺点: 计算效率较低,尤其是解决设计和优化问题时计算效率更低,过程系统模拟的面向方程法,EOM又称联
3、立方程法,将描述整个过程系统的数学方程式联立求解,从而得出模拟计算结果,面向方程法可以根据问题的要求灵活地确定输入、输出变量,而不受实际物流和流程结构的影响 EOM解算过程系统模型快速有效,对设计、优化问题灵活方便。效率较高 EOM的形成通用软件比较困难;不能利用现有大量丰富的单元模块;缺乏实际流程的直观联系;计算失败之后难于诊断错误所在;对初值的要求比较苛刻;计算技术难度较大,过程系统模拟的联立模块法,联立模块法又称双层法,将过程系统的近似模型方程与单元模块交替求解 联立模块法兼有序贯模块法和面向方程法的优点。既能使用序贯模块法积累的大量模块,又能将最费计算时间的流程收敛和设计约束收敛等迭代
4、循环合并处理,通过联立求解达到同时收敛,过程系统稳态模拟三种方法的比较,思考题,结合上节课的例子,说明什么是序贯模块法,什么是面向方程法,什么是联立模块法。 结合上例,比较三种过程系统模拟方法的优缺点,2.2 序贯模块法,2.2.1 序贯模块法的基本原理,序贯模块法的基础是单元模块,通常单元模块与过程单元是一一对应的 单元模块是依据相应过程单元的数学模型和求解算法编制而成的子程序 单元模块具有单向性特点。,过程单元与单元模块,序贯模块法的基本思想,从系统入口物流开始,经过接受该物流变量的单元模块的计算得到输出物流变量,作为下一个相邻单元的输入物流变量。 依此逐个的计算过程系统中的各个单元,最终
5、计算出系统的输出物流。计算得出过程系统中所有的物流变量值,也即状态变量值,序贯模块法的求解与过程系统的结构是有关的。 具有反馈联结的系统(不可分割子系统),需要用到断裂(Tearing)和收敛(Convergence)技术,具有反馈的系统与收敛单元,通过断裂技术可以打开回路,以便采用序贯模块法进行求解。在断裂物流处设置一个收敛单元。 对于复杂系统,收敛单元设置的位置不同,其效果也将不同。 最优设置要通过断裂技术去解决。 如何得到新的变量值,如何保证计算收敛,如何加快收敛,取决于收敛算法,与断裂物流变量的特性有关。,2.2.2 再循环物流的断裂,(1)断裂的基本概念,求解方法,1. 联立方程组
6、2. 通过断裂进行降维 X2 X3 X4 X1 X2 把一个四维求解问题降阶成为了四个一维问题,从而减化了计算难度。 这种通过迭代把高维方程组降阶为低维方程组的办法称为“断裂”。,有向图,选择不同的断裂物流,相应的迭代序列也不一样:,由于系统中各物流及其变量特性的不同,在收敛计算上常是有很大差异的。 如何选择断裂物流、确定迭代序列,是实施序贯模块法进行过程系统模拟计算中必须要解决的问题。,(2)断裂方法的研究,六十年代初,Rubin就提出了断裂的思想 判断最佳断裂的准则分为四类 断裂的物流数最少; 断裂物流的变量数最少; 断裂物流的权重因子之和最少; 断裂回路的总次数最少。,i=1,m,代表回
7、路;j =1,n,代表物流,约束方程的含义是每个回路至少要被断裂一次。 准则设定j =1 ; 准则令j为物流变量数; 准则中j为可根据物流性质而取的选择值,如物流变量对计算过程灵敏度大小的估计值; 准则的j等于每个断裂物流所切断的回路总数。,(3)回路矩阵,要断裂再循环物流,必须先识别再循环回路,并借助一定的方法描述它们。 一个不可分隔子系统包含若干个再循环回路。包含两个以上再循环物流,且其中的任何单元只被通过一次,称作简单回路(Simple Cycle)。,S1S2S4S2S5构成的回路不是一个简单回路,因为其中的单元和单元被通过了两次。,过程系统中的简单回路可以用回路矩阵(Loop/Str
8、eam Matrix)表示。 表示方法: 矩阵中的行代表回路,列代表物流。 若某回路I中包括有物流j则相应的矩阵元素aji=1,否则为空白或零。,(4)UpadyheGrens断裂法,美国加州大学的Upadhye等提出的,一种类似动态规划法的寻求最佳断裂物流的算法。 为了对该不可分隔子系统的高维求解进行降维运算,须将该子系统中的某些回路进行断裂。达到断裂的方案并不是唯一的。 需要解决的两个问题:一是要有一种能把所有的有效断裂物流组都能搜索出来的办法;二是要能把最优断裂组从中选择出来。,替代规则:,令D1为一有效断裂组,Ai为全部输入流均属于D1的单元(至少有一个这样的单元存在,否则D1为无效断
9、裂组)。将Ai的所有输入流用Ai的全部输出流替代,构成新的断裂组。令得到的新的断裂组为D2,则,D2也是有效断裂组; 对于直接迭代,D2与D1具有相同的收敛性质。,由替代规则联系起来的所有断裂组的集合,定义为断裂族,(5)寻求最优断裂组的算法,从任一有效断裂开始,运用替代规则: 如果在任何一步中出现有两次被断裂的物流(二次断裂组),则消去其中的重复物流。消去重复后断裂组则作为进行下一步的新起点。 重复步骤、,直到不再有二次断裂组出现,且每个“树枝”上有重复的断裂组出现时为止。从最后一个新的起点开始,其后出现的所有不重复的断裂组成为非多余断裂族。 非多余断裂族中总权最小的断裂组为最优断裂组。,非
10、多余断裂族: 断裂组S1,S4,S7为最优断裂组。,通过断裂可以把不可分隔子系统中的回路物流打开,从而可以利用序贯模块法对该过程系统进行模拟计算。这种模拟计算的开始是首先要设定起始物流变量的猜值,计算的终点则在于该猜值与计算值的收敛。,2.2.3 断裂物流变量的收敛,执行断裂物流变量收敛功能的模块称收敛单元模块 不可分隔子系统的断裂物流 断裂物流变量的收敛问题实际上是个迭代求解非线性方程组的问题x = y = G(x),当断裂物流变量猜值x与计算值y之差小于收敛容差时 y x = G(x) x 则x为断裂物流变量的收敛解 收敛单元的功能总计有如下作用: ()获取猜值的初值x0 ()根据计算值y
11、,以一定的方法确定新的猜值x ()比较猜值x和计算值y,若其结果满足给定精度要求则结束迭代计算,否则继续迭代计算过程,可见,收敛单元实质上就是一个数值迭代求解非线性方程组的子程序。求解非线性方程组的数值计算方法很多,适合于收敛单元的数值计算方法一般应尽可能满足下列要求: 1. 对初值的要求不高。 1).初值易得,不易引起迭代计算的发散; 2)初值的组数少。,例: 用直接迭代法求解下列方程组,解:令猜值为X12;X210;X35,解:令猜值为X16;X23.5;X35,2. 数值稳定性好,3. 收敛速度快 对收敛速度的影响主要有三个因素:迭代次数;函数G(x)的计算次数;矩阵求逆的次数。 每计算
12、一次函数值就相当于做一次流程回路的模拟计算,每求一次导数就要做两次流程模拟计算 对于断裂物流的收敛,好的非线性方程组的数值迭代次数少,而且应该尽量避免导数计算和矩阵求逆 4. 占用计算机存储空间少,直接迭代法,直接迭代法是将计算值yk作为下一轮迭代的猜值xk+1而实施迭代计算 xk+1G(xk) 非线性方程组x = y = G(x) 的另外一种形式为 F(x)xG(x)= 0 xk+1xkF(xk) 与牛顿公式相比较:,直接迭代法的雅可比矩阵为单位矩阵 优点:方法简单, 只需要一组初值,不需计算导数和逆矩阵 弱点:迭代次数多、收敛速度慢, 对初值要求较高 为改善直接迭代法的收敛行为,提出了阻尼
13、直接迭代法,或称加权直接迭代法:,为改善直接迭代法的收敛行为,提出了阻尼直接迭代法,或称加权直接迭代法: q为阻尼因子,可以人为给定 q=0 直接迭代 0q1 加权直接迭代,可改善收敛的稳定性 q0 外推直接迭代,加速收敛,但稳定性下降 q1 无意义,当闪蒸温度分别为以下值时: 分别用直接迭代法和阻尼直接迭代法计算汽相和液相产品的流量和组成,阻尼因子分别取值为0.5,0.3,-0.2,-0.3,-0.7,-0.9 解:依据闪蒸条件,设该闪蒸过程为理想体系,三个闪蒸器均为等温闪蒸过程,建成相应的单元模块。并将其改绘为如下三级闪蒸过程模拟模块流程,平衡闪蒸单元模型,从上例可见,阻尼因子q值的选取具
14、有较大的任意性和经验性。1958年Wegstien提出了一种简便的方法,可以弥补这种阻尼因子取值困难的弱点,Wegstein法,一维Wegstein法 求解一维代数方程 (1) Wegstein迭代公式为: (2) 其中: (3),对于隐式一维代数方程: (4) 相应的迭代公式称作割线法,其迭代公式可从Wegstein迭代公式导出 从(4)式可得出: (5) 将上式代入(2)式,(6),从(2)和(3)式得到: (7) 上式代入(5)式,则有: (8) 式(8)就是割线法的迭代公式。由此可见,Wegstein法与割线法是相通的,隐式方程具有更大的普遍性,所以割线法常为人们所熟知。 在流程模拟领域中,物流回路多用显式方程描述的。因而多用Wegstein法。,由(2)式可见,一维Wegstein法需要有两个初值,其中第一个初值是设置的猜值,第二个初值可根据第一个初值按直接迭代法得到。,