多组分物系的泡点和露点计算.ppt

2.2多组分物系的泡点和露点计算,泡点计算：泡点温度计算：已知xi，p，求解T，yi泡点压力计算：已知xi，T，求解P，yi,露点计算：露点温度计算：已知yi，P，求解T，xi露点压力计算：已知yi，T，求解P，xi,2.2.1泡点温度和泡点压力计算,已知xi(i=1,2,c),P或T,计算方程有：,yi,T或P,相平衡关系,浓度和式,汽液平衡常数关系式,yi=Kixi(i=1,2,c),C个,1个,1个,Ki=f(P,T,xi,yi),C个,2C+2个,2.2.1泡点温度和泡点压力计算,总变量：C-xi，C-yi，C-Ki，T，P。共3C+2个已知C个变量：（C-1）个组分和T或P未知变量为2C+2个由上述2C+2个方程求解（唯一解）,需试差计算,2.2.1泡点温度和泡点压力计算,1、泡点温度计算,泡点方程：,已知xi,P,或,1、泡点温度计算,计算方法,1.简化计算（手算）（Ki与组成无关)(1)试差法(2)相对挥发度法2.严格计算法（计算机计算）（Ki与组成有关)Ki=f(T,P,xi,yi),1、泡点温度计算,(1)简化计算,Ki=f(T,P),P已知，Ki=f(T),气相理想气体或理想溶液常压或中压下烃类混合物液相理想溶液各组分分子结构相似,Ki与组成无关,已知xi,p,1、泡点温度计算,1）试差法,已知xi,P,设T,P-T-K图Ki=f(T，P，xi，yi)iKi,Ki,设yi,Y,T，yi,N,1、泡点温度计算,2）相对挥发度法方法依据：相对挥发度随温度变化比Ki小得多，在一定温度范围内可视为常数.,相对挥发度定义：,ki组分中任选的一个基准组分,计算式推导：,1、泡点温度计算,相对挥发度法,公式推导：,1、泡点温度计算,相对挥发度法,由,可计算,T,通常不需要试差！,1、泡点温度计算,相对挥发度法计算步骤,已知xi,P,设T,Pis,Pks,.i=1,2,kc,N,Y,Tb,1、泡点温度计算,（2）严格计算法（计算机计算）,对于非理想性较强的物系：,Ki=f(T,P,xi,yi),1、泡点温度计算,严格计算法计算过程,(一般受P的影响较小),计算或yi时需迭代（因T未知）。当P不太高时，P的影响不太大，收敛较快。,一般受压力的影响较小,T,P,T,P，yi,T,P，xi,1、泡点温度计算,严格计算法计算方法,牛顿迭代法：,一阶导数,二阶导数,T初值可由相对挥发度求得,也可取某组分i的沸点为初值。,收敛快！,泡点温度计算框图（P50图2-5）,开始,调整T,输入P、及有关参数,计算Ki和yi,是否第一次迭代,yi有无变化,输出T、,计算yi,lnyi,结束,设定T并令1作第一次迭代,圆整计算,计算Pis、ViL、is、,纯组分i维里方程，wilson(T,P,x),G（1T）非线性,高压下的烃类（非理想）T的调整采用抛物线法（Muller）；G（1T）线性,P20，用牛顿拉普森（Neuton-Raphson）迭代法（线性）切线法。,先圆整后计算,例2-3,例2-4确定含正丁烷（1）1.15、正戊烷（2）0.4和正己烷（3）0.45（均为摩尔分率）之烃类混合物在0.2MPa压力下的泡点温度。,解：,因各组分都是烷烃，所以汽、液相均可看成理想溶液，平衡常数只取决于温度和压力。,如计算要求不甚高，可使用烃类的P-T-K图（见图2-1）。,例2-3,假设T=50,因P=0.2MPa,查图求平衡常数：,说明所设温度偏低。重设T=58.7,例2-3,假设T=58.7,查图求平衡常数：,故泡点温度：T=58.7,例2-4,例2-4某厂氯化法合成甘油车间，氯丙烯精馏二塔的釜液组成为：3-氯丙烯0.0145，1,2-二氯丙烷0.3090，1,3-二氯丙烯0.6765（摩尔分率）。塔釜压力为常压，试求塔釜温度。各组分的饱和蒸汽压数据为：（Ps:kPa；T:）：,3-氯丙烯,1,2-二氯丙烷,1,3-二氯丙烯,例2-4,已知：3-氯丙烯x1=0.0145,1,2-二氯丙烷x2=0.3090,1,3-二氯丙烯x3=0.6765，P常压。求：Tb,yi（i=1,2,3）,解：,气相为理想气体,常压,液相为理想溶液,结构相似,参考解法：（1）试差法（手算）（2）相对挥发度法（3）计算机求解,例2-4,（1）试差法-计算（见下表）,塔釜温度应为100。,例2-4（2）相对挥发度法-计算,以1,2-二氯丙烷为相对挥发度之基准组分k，则根据98时各组分之蒸汽压值可得出相对挥发度：,t=99.7,antoin方程,例2-4,相对挥发度法,虽然采用了98的相对挥发度数据，但计算结果与前法十分接近，完全满足工程计算的要求，而不需要试差。相对挥发度也可用于修定平衡常数法泡点计算的迭代温度。,由式（2-53）导出：,例2-4,（3）计算机计算按单参数牛顿迭代求解,Antoine常数,第k+1次和第k次迭代温度,例2-4,可选组分1或组分3的沸点为初值开始计算：,（3）计算机计算按单参数牛顿迭代求解,例2-4,牛顿迭代计算结果为：,达到迭代精度要求，故泡点温度为99.812。,迭代5次！,（3）计算机计算按单参数牛顿迭代求解,例2-4,若用Richmond算法，还需求二阶导数,（3）计算机计算按单参数牛顿迭代求解,例2-4,计算结果,已达到牛顿法t5的精度，故t3为所求。,迭代3次！,（3）计算机计算按单参数牛顿迭代求解,例2-5,例2-5丙酮（1）-丁酮（2）-乙酸乙酯（3）三元混合物所处压力为2026.5kPa,液相组成为x1=x2=0.3,x3=0.4（摩尔分率）。试用活度系数法计算泡点温度（逸度系数用维里方程计算；活度系数用Wilson方程计算）。,Wilson参数为（J/mol）,各组分的液相摩尔体积cm3/mol,例2-5,安托尼方程常数,Ps：kPa，t：K,组分的临界参数和偏心因子,例2-5,解,用Abbott公式计算第二维里系数,交叉临界性质用Lorentz-Berthelot规则求出,逸度系数计算用式（2-30）,由于系统压力接近于各组分的饱和蒸汽压，普瓦延因子近似等于1，故平衡常数公式简化为,计算步骤见图2-2,例2-5,设泡点温度初值为500K，计算机计算中间结果如下：,例2-5,最终迭代泡点温度468.7K下各变量数值如下表所示：,2、泡点压力计算,泡点压力方程：,已知xi,T，计算yi,P,2、泡点压力计算,（1）简化计算,完全理想体系,气相理想气体，液相非理想溶液,2、泡点压力计算,（2）严格计算,压力不太高,T已知,T已知，P，yi试差！当压力不太高时,P对的影响不太大,故收敛较快.,T,xi已知（P影响不大）,外循环迭代Pb,图2-3泡点压力计算框图,内循环迭代yi,2、泡点压力计算,（2）严格计算,压力较高时,状态方程求解,图2-4泡点压力计算框图,已知氯仿（1）-乙醇（2）溶液浓度为x1=0.3445(摩尔分率)，温度为55。试求泡点及气相组成。该系统的Margules方程式常数为：A12=0.59，A21=1.42。55时，纯组分的饱和蒸汽压P1S=82.37kPa,P2S=37.31kPa,第二维里系数：B11=-963cm3/mol,B22=-1523,B12=-1217。指数校正项可以忽略。,例2-6（P52）,例2-6,解：将式（2-35）代入（2-50）中，并忽略指数项，得：,同理,因为1及2是P及y的函数，而P及y又未知，故需用数值方法求解。,在x1=0.3445时，由Margules方程式求得：,例2-6,同理,例2-6,先假设1=2=1,由式（B）求出P及Py1和Py2,Pyi/P=yi,故得y1及y2,P=Py1+Py2,以P及y1、y2为初值,(C)和(D)算出1及2,1及2算出新的P及y1、y2,反复进行，直至算得P与y和假设值相等（或差值小于规定值）时为止,以1=2=1和P1S、P2S代入式（B），得：P1.64750.344582.371.04020.655537.3146.7525.4472.19y1=46.75/72.19=0.6476;y2=0.3524将P、y1和y2值代入式（C）和（D）求出：11.0038；20.9813,例2-6,以上述1和2的值代入式（B），得：y1P=46.57故y1=0.6424y2=0.3576由于此次计算结果与第一次试算结果相差甚小，故不再继续算下去。因此P75.50kPa,y1=0.6424,y2=0.3576。,例2-6,2.2.2露点温度和压力计算,露点方程：,yi，TorP,xi,PorT,2.2.2露点温度和压力计算,计算方法：,（1）简化计算,（2）严格计算,具体计算方法同泡点计算！,例2-7,例2-7乙酸甲酯（1）-丙酮（2）-甲醇（3）三组分蒸汽混合物的组成为y1=0.33,y2=0.34，y3=0.33(摩尔分率)。试求50时该蒸汽混合物之露点压力。,解：汽相假定为理想气体，液相活度系数用Wilson方程表示。由有关文献查得或回收的所需数据为：,50时各纯组分的饱和蒸汽压，kPa,例2-7,50时各组分液体摩尔体积，cm3/mol,由50时各两组分溶液的无限稀释活度系数回归得到的Wilson常数：,取xy初值,内层：迭代xi外层：迭代P,例2-7计算框图,由x和求,例2-7-解,从多组分Wilson方程,同理,假定x值，取x1=0.33，x2=0.34，x3=0.33。按理想溶液确定p初值,例2-7-解,求,同理,例2-7-解,求,圆整得,在P=71.916kPa内层迭代汇总如下：,例2-7-解,调整P,在新的P下重复上述计算，迭代至P达到所需精度。最终结果：露点压力85.101kPa平衡液相组成：x1=0.28958x2=0.33889x3=0.37153,上述计算一般不能依靠手算，而必须利用计算机。若省略普瓦延因子，则可节省机时。Ki对xi的变化敏感，对压力的变化不敏感，因此，内层迭代xi,外层迭代P的计算是合理的。,例2-7-解,2.2小结,1.泡点计算：已知xi，p或T，求解T或p，yi,2.泡点方程：,或,2.2小结,或,4.露点方程,3.露点计算：已知yi，P或T，求解T或P，xi,2.2小结,5.计算方法,（1）简化计算（Ki与组成无关)试差法相对挥发度法（2）严格计算法（计算机计算）（Ki与组成有关)Ki=f(T,P,xi,yi),2.2小结,已知xi,P,设T,P-T-K图Ki=f(T，P，xi，yi)iKi,Ki,设yi,Y,T，yi,N,试差法,2.2小结,已知xi,P,设T,Pis,Pks,.i=1,2,kc,N,Y,Tb,相对挥发度法,作业,P86(7)，(11)_每周五交作业！,