第二章1自由度分析及系统分解课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,主要内容,2.1,单元模型的自由度分析,2.2,化工单元模型和模块,2.3,过程系统的结构分析,主要内容2.1 单元模型的自由度分析,1,2.1,自由度分析,化工单元的数学模型是单元的输入变量、输出变量和过程参数用某种机理的、半机理的或统计的函数关系关联而成的线性的或非线性的方程组。,F(X)=0,其中,X=,（,x,1,x,2,x,n,)T,为完整描述对象所需的变量向量，,F=,（,f,1,f,2,，,f,m,）为函数关系向量。,2.1 自由度分析 化工单元的数学模型是单元的输入变,2,机理模型,由过程机理出发，经推导得到，并得到实验验证。一般只有对过程内在规律明确且相对比较简单的研究对象才能建立机理模型。,例如：流体在圆管中作层流流动时其压力,p,变化的关系式：,p,2,-p,1,=32,l u,/ d,2,机理模型是对实际过程直接的数学描述，是过程本质的反映，因此其结果可以外推。,机理模型 由过程机理出发，经推导得到，并得到实验验证。,3,统计模型,数学模型也可以根据实验装置、中型或大型工业装置的实测数据，通过数据回归分析得到纯经验的数学关系式，这就是统计模型。,Nu=0.023Re,0.8,Pr,0.3,统计模型与过程机理无关，是根据实验从输入和输出变量之间的关系，经分析整理得到。它只是在实验范围内有效，不具有外推性。,统计模型,4,混合模型,对实际过程进行抽象概括和合理简化，然后对简化的物理模型加以数学概述，这样得到数学关系式称为混合型模型，它是半经验半理论化性质的。,混合模型,5,自由度分析的主要目的：在系统求解之前，确定需要给定多少个变量，可以使系统有唯一确定的解。,F(X)=0,单元操作过程的数学模型含有,m,个独立方程，其中含有,n,个变量，且,nm,，则该模型具有的自由度为,d=n-m,即需要在,n,个变量中，给定,d,个变量的值，其余的,m,个变量可通过,m,个方程求解。,自由度分析的主要目的：在系统求解之前，确定需要给定多少个变量,6,为了求解单元数学模型而确定的,d,维独立变量称为决策变量。,在,m,个变量中究竟选取哪,d,个，具有一定的“自由度”。,d,个变量选择原则：,1,）选那些受限制较多的变量,.,如冷却水的温度、流量等，它们受当地气候和水资源条件的限制。又如高温状态下物料的温度将受设备材料耐温性能的限制。,2,）选出的变量，赋值后，可是系统模型方程的求解，最为方便、容易。,为了求解单元数学模型而确定的d维独立变量称为决策变量。,7,本节所讲的自由度概念与“物理化学”课程中提到的自由度不同。,在“物理化学”中，相律用下式表达体系的自由度,d=C-P +2,式中,C,组分数,P,想数目,“相律”中的自由度只涉及强度性质（,T,，,P,等）而不涉及系统的大小数量。,在化工单元操作模型中，必须考虑系统的大小量、如流股的质量流率、热负荷及压力变化等。,注意,本节所讲的自由度概念与“物理化学”课程中提到的自由度不同。注,8,在“化工热力学”中，杜赫姆（,Duhem,）定理指出，对于一个已知每个组分初始质量的封闭体系，其平衡状态取决于两个独立变量，而不论体系中有多少个相、多少个组分或多少个化学反应。,根据该定理，可推知一个独立流股具有（,C+2,）个自由度。,如规定流股的中,C,个组分的摩尔流量以及流股的温度,T,和压力,P,，则该流股就确定了。,注意：杜赫姆定理定义自由度时含有归一方程,xi=1.,所以不含归一方程，则独立流股的自由度为,C+3,在“化工热力学”中，杜赫姆（Duhem）定理指出，对于一个已,9,2.2,单元模型的自由度分析,单元模型方程的分类,物料衡算方程,指直接表达物料守恒关系式的方程,热量衡算方程,指直接表达能量守恒关系即热力学第一定律,设备约束方程,物料通过特定设备时受具体设备的约束，按照特定的关系发生变化。,其他,如混合物中各组分摩尔分率之和必须等于,1,的关系（摩尔分率约束方程）等。,2.2单元模型的自由度分析单元模型方程的分类,10,（,1,）混合器,图中有,2,个流股合并成,1,个流股，每个流股有,C+2,个独立变量，因此总的变量数有,3(C+2,）个。,（1）混合器图中有2个流股合并成1个流股，每个流股有C+2个,11,对该过程可以建立以下独立方程：,压力平衡方程,P,3,= min(P,1,P,2,),物料衡算方程,F,3,= F,1,+ F,2,热量衡算方程,F,1,H,1,+F,2,H,2,=F,3,H,3,式中,H,流股的比摩尔焓,F,流股的摩尔流量,x,流股中组分的摩尔分率,p,压力,独立方程数,C+2,对该过程可以建立以下独立方程：压力平衡方程 P3,12,混合器的独立方程数,m = C+ 2,混合器的自由度,d = n m =3(C+2)-(C+2)= 2(C+2),如果有,S,股输入物流，则自由度为,S(C+2),，即相当于,S,个输入流股变量后，混合器出口流股的变量也就确定。,混合器的独立方程数,13,由直观分析可知，当指定一股输入物流的变量,(C+2),以及一个分割率（其值在,0,1,之间），则该分割器的两股输出物流的变量就完全确定了，即该简单分割器的自由度为（,C+2,）,+1,。,（,2,）分割器,F1 C+2,F3,F2,r,由直观分析可知，当指定一股输入物流的变量(C+2)以及一个分,14,独立变量数,3(C+2)+1,过程参数,(,分流比,U) 1,个,2,个温度等式,T1=T2, T1=T3,2,个压力等式,P1=P2, P2=P3,2,个物料平衡方程,F2= U*F1, F3=(1-U)*F1,组份等式,x,1i,=x,2i,(i=1,2,，,C-1),x,1i,=x,3i,(i=1,2,，,C-1),独立方程数,2C+4,独立变量数3(C+2)+1, 过程参数(分流比U) 1个独立,15,当一个流股分成,S,个流股，指定输入流股变量,C+2,个以及,S-1,个分割分率值，则可由,S(C+2),个独立方程式解出,S,个分支流股包含的变量。这样该分割器的自由度为,d=(S+1)(C+2)+(S-1)-S(C+2)= (C+2)+(S-1),当一个流股分成S个流股，指定输入流股变量C+2个以及S-1个,16,（,3,）闪蒸器,F1,，,T1,，,P1,F,2,，,T,2,，,P,2,F3,，,T3,，,P3,Q,（3）闪蒸器F1，T1，P1F2，T2，P2F3，T3，P3,17,闪蒸器的加热量,Q,必须作为设备参数。故总变量数为,3,（,C+2,）,+1,，表示闪蒸器变量间关系的方程为：,（,3,）闪蒸器,物料衡算方程,x,1i,*,F,1,= x,2i,*,F,2,+ x,3i,*,F,3,(i=1,2,.,C),热量衡算方程,F,1,*,H,1,+F,2,*,H,2,=F,3,*,H,3,温度平衡方程,T,2,=T,3,压力平衡方程,P,2,=P,3,相平衡方程,x,2i,= k,i,*,x,3i,(i=1,2,.,C),共有,2C+3,个独立方程,故自由度为,d= 3(C+2)+1 (2C+3) =,（,C +2,）,+2,规定闪蒸压力，闪蒸温度，气化量，气化率等,闪蒸器的加热量Q必须作为设备参数。故总变量数为3（C+2）+,18,（,4,）换热器,F1,，,T1,F3,，,T3,Q,F4,，,T4,F2,，,T2,（4）换热器F1，T1F3，T3QF4，T4F2，T2,19,换热器的热负荷可作为设备参数由设计规定。对换热器的热流一侧和冷流一侧可各写,C+2,个关系方程式。故自由度为,d=4(C+2)+1 2(C+2)= 2(C+2) + 1,即当给定进口冷、热物流的变量,2(C+2),个以及换热器的热负荷,Q,（,1,个变量）后，出口流股的变量就完全可由,2(C+2),个独立方程式求出。,（,4,）换热器,换热器的热负荷可作为设备参数由设计规定。对换热器的热流一侧和,20,物料衡算,F1=F3,，,X,1i,= X,3i,（,i=1,2,，,，,C-1,）,F2=F4,，,X,2i,= X,4i,（,i=1,2,，,，,C-1,）,焓方程,H1+H3=H2+H4,传热方程,Q=K*A*T,压力方程,P3=f,1,(P1), P4=f,2,(P2),物料衡算,21,（,5,）反应器,F1,，,T1,，,P1,F2,，,T2,，,P2,Q,1,2,，,P,常用的反应器模型是规定出口反应程度的宏观模型。不假定反应达到平衡，而是规定了,r,个独立反应的反应程度（,1,，,2,，,，,r,),。向反应器提供的热量,Q,（移出时,Q,为负值）和反应器中的压降,P,是两个设备单元参数，所以共有,r+2,个设备单元参数，独立方程数为,C,个组分物料平衡方程，,1,个焓平衡方程，,1,个压力平衡方程，即独立方程数为,C+2,个。其自由度,d,d=2(C+2)+(r+2) - (C+2 ),=(C +2)+ r + 2,（5）反应器F1，T1，P1F2，T2，P2Q1 2,22,压力变化单元包括阀门、泵、压缩机等。压力单元中除了压降,P,作为设计参数规定外，对于泵、压缩机而言，与物料无关的能量流（轴功）也作为设计参数予以规定；独立方程数为,C,个组份物料平衡方程，,1,个温度相等（忽略温度变化）方程，,1,个压力方程，即独立方程总数为,C+2,个。,（,6,）压力变化单元,压力变化单元包括阀门、泵、压缩机等。压力单元中除了压降P作,23,阀门自由度为,d=2(C+2)+1-(C+2)=(C+2)+1,泵、压缩机自由度,d=2(C+2)+2-(C+2)= (C+2)+2,P,C+2,C+2,C+2,1,C+2,C+2,C+2,2,阀门自由度为 d=2(C+2)+1-(C+2)=(C+2),24,通过对上述过程单元的自由分析，可以归纳出过程单元的自由度计算通式：,通过对上述过程单元的自由分析，可以归纳出过程单元的自由度计算,25,过程系统的自由度分析,过程系统的自由度分析,26,混合器,5,压缩机,5,分流器,10,反应器,5,换热器,8,阀,5,闪蒸,9,放空,产物,3,5,5,5,1,3,2,1,4,1,5,5,5,1,3,5,5,1,）有箭头数字代表,Ci+2,2,）无箭头数字代表设备数,3,）方框内数字代表独立方程数,5,混合器5压缩机5分流器10反应器5换热器8阀 5闪蒸9放空产,27,第二章1自由度分析及系统分解课件,28,以上例子中系统的自由度,d=Nv-NE,（其中,Nv,，,Ne,分别为系统总变量数和总方程数）确定，这样对于复杂的系统，其自由度的确定相当繁琐。,过程系统的自由度可在过程单元自由度分析基础上，由下式确定：,以上例子中系统的自由度d=Nv-NE（其中Nv，Ne分别为系,29,这个结论是基于每增加一个联结流股，就相应地增加,C,j,+2,个联结方程这一事实得出的。,联结流股变量数可以流股组份,Cj,表示,k,j,(L),= C,j,+ 2,上述例子中联结流股有,7,条，每个流股的变量数是,5,，所以,k,j,(L,),=,5 7 =35,这个结论是基于每增加一个联结流股，就相应地增加Cj+2个联结,30,第二章1自由度分析及系统分解课件,31,该系统中自由度之和为,d,i,(u),=51,该系统中自由度之和为di(u)=51,32,实际建模的自由度分析：,实际建模的自由度分析：,33,平衡闪蒸的系统模型与分析,平衡闪蒸的系统模型与分析,34,第二章1自由度分析及系统分解课件,35,2.3,过程系统的结构分析,现代化的大型化工企业是一个规模庞大、结构复杂、循环嵌套、影响因素众多的大型过程系统。描述这样的系统要用成千上万个方程式，其中常常会出现某些必须同时求解的非线性，代数、微分方程混杂的方程组，当方程组的维数很高时，即使使用计算机求解也存在一定的困难。因此有必要采用结构分析的方法将系统分解，把大系统分解成若干个相互独立的子系统，然后按一定的次序计算、迭代求解。,2.3 过程系统的结构分析现代化的大型化工企业是一个规模庞大,36,谢谢你的阅读,知识就是财富,丰富你的人生,谢谢你的阅读知识就是财富,37,第二章1自由度分析及系统分解课件,38,