二元汽液平衡数据测定实验数据处理.docVIP

二元汽液平衡数据测定实验数据处理： 原理： 以循环法测定气液平衡数据的平衡釜基本原理相同，如图1所示，体系达到平衡时，两个容器的组成不随时间变化，这时从A和B 两容器中取样分析，即可得到一组平衡数据。 当达到平衡时，除了两相的压力和温度分别相等外，每一组分的化学位也相等，即逸度相等，其热力学基本关系为： ……………………………………….(1) 常压下，气相可视为理想气体，；再忽略压力对液体逸度的影响，从而得出低压下气液平衡关系为： ……………………………………(2) 式中 p——体系压力（总压）； pi0——纯组分I在平衡温度下饱和蒸汽压，可用安托尼（Antoine）公式计算； xi、yi——分别为组分I在液相和气相中的摩尔分率； γi——组分i的活度系数 由实验测得等压下气液平衡数据，则可用 ……………………………….. ….. (3) 计算出不同组成下的活度系数。 本实验中活度系数和组成关系采用Wilson方程关联。Wilson方程为： lnγ1=-ln………(4) lnγ2=-ln………(5) Wilson方程二元配偶参数和采用非线形最小二乘法,由二元气液平衡数据回归而得. 目标函数选为气相组成误差的平方和,即 F=…………………………………(6) 数据处理（以表1的实验数据为例子） Ⅰ.实测温度及苯的气液相组成 大气压：101.33Kpa 表1 苯—正庚烷混合液平衡温度及气、液相组成的测定数据 组号 主温度 t主/（℃） 辅助温度 ts/（℃） 气相(苯) 液相（苯） 折射率 yi/mol% 折射率 xi/mol% 1 82.0 22.0 1.4610 81.9 1.4510 74.6 2 83.4 23.0 1.4475 72.9 1.4361 63.6 3 84.8 24.0 1.4398 67.0 1.4235 53.7 4 86.2 25.0 1.4310 59.4 1.4155 46.8 5 87.6 28.1 1.4238 52.6 1.4038 38.4 表2 水银温度计检定结果 温度计示值（℃） 80 85 90 95 100 修正值（℃） +0.08 +0.07 +0.07 +0.09 +0.08 Ⅱ. 平衡温度及纯物质的饱和蒸气压的计算： ① 平衡温度的计算：（使用热电偶测温时略去此步骤） t实际值= t主+t校正 t校正=kn(t主-ts) matlab运算结果如下： t0=[82.0 83.4 84.8 86.2 87.6]; t2=[22.0 23.0 24.0 25.0 28.1]; t1=0.00016*60*(t0-t2); t3=t0+t1 结果： t3 = 82.5760 83.9798 85.3837 86.7875 88.1712 其中t3 代表平衡温度,结果填入表3中 ② 饱和蒸气压的计算 lgPi0=Ai – Bi/(Ci+t) 式中： t—温度，℃ P0—饱和蒸汽压，mmHg Pi0(苯)=10^(6.87987-(1196.76/(219.161+t))) Pi0(正庚烷)=10^(6.89386-(1264.37/(216.64+t))) 用matlab运算结果如下 t=[82.5760 83.9798 85.3837 86.7875 88.1712]; for i=1:5 p1(i)=10^(6.87987-(1196.76/(219.161+t(i)))) p2(i)=10^(6.89386-(1264.37/(216.640+t(i)))) end 结果： p1 =819.6614 855.0698 891.6626 929.4612 967.9270 p2 =465.8546 487.5090 509.9562 533.2131 556.9509 Ⅲ．用Wilson方程计算如下: 表3 平衡温度及饱和蒸气压 组号 平衡温度t /(℃) Pi0(苯) /(mmHg) Pi0(正庚烷) /(mmHg) 1 82.58 819.6614 465.8546 2 83.98 855.0698 487.5090 3 85.38 891.6626 509.9562 4 86.79 929.4612 533.2131 5 88.17 967.9270 556.9509 ① Wilson方程参数的确定 方程格式： lnPy1/P10x1= -ln(x1+λ12x2)+x2(λ12/(x1+λ12x2)- λ21/ (x2+λ21x1)) lnPy2/P20x2= -ln(x2+λ21x1)+x1(λ21/(x