拟三元体系汽液平衡测定Word文件下载.docx
《拟三元体系汽液平衡测定Word文件下载.docx》由会员分享，可在线阅读，更多相关《拟三元体系汽液平衡测定Word文件下载.docx（17页珍藏版）》请在冰点文库上搜索。
对
作图，如图1所示，然后计算偏差
，
图1气—液平衡数据校验图
将
与另一数量
进行比较，
，对于乙醇-水体系来说，
为最低沸点温度与水的沸点差的绝对值，
为
之间的最低沸点。
采用的判别标准为：
若
一般认为该数据是符合热力学同一性的。
这种判别标准虽然是近似的，仍然是很有用的经验方法而被广泛采用，但是该法对数据进行的是整体检验而非逐点检验，不仅需要整个浓度范围内的数据，而且有可能因为实验误差的相互抵消造成面积检验虽然符合要求，但数据并非热力学一致的情况。
对此，可以采用微分检验法（逐点检验法）可以克服这一缺点，但本实验中不要求。
对于乙醇-水-乙二醇（KAc）体系，涉及4种组分，为简化后续的实验数据处理，实验安排和设计中存在以下假设：
（1）KAc为不挥发物质，因此，本实验中四元体系可看作拟三元体系；
（2）乙二醇的饱和蒸汽压很高，并与乙醇，水的饱和蒸汽压差别比较大，气相组成中一般不含有乙二醇；
（3）气液平衡釜的体积比较大，液相量大，气相取样量很少，对釜液相组成影响很小，因此，简化实验数据测试，只测定气相组成，液相组成不变。
本实验测定常压（1atm）条件下，当温度稳定时，通过气相色谱测定平衡时的汽相组成。
图2实验装置图
3.实验药品和仪器
实验药品：
无水乙二醇（KAc）（分析纯）、无水乙醇（分析纯）、蒸馏水。
实验仪器：
Othmer气液平衡仪（见图2）、气相色谱仪、可调压电热套、250mL量筒、100mL量筒、5mL移液管、1mL移液管、25mL移液管、洗耳球。
Othmer气液平衡仪使用介绍
Othmer气液平衡仪是目前在科研单位广泛使用的气相单循环式气液平衡测定装置，具有平衡快（约40分钟），取样方便等特点。
取样时，用1mL的移液管插入冷凝器上口，迅速地从气相冷凝液积存处吸取少量液体后拔出，将液体挤入样品瓶中即可。
4.实验步骤
（1）在气液平衡仪内加入100ml无水乙二醇（KAc）；
（2）准确用量筒量取250ml分析纯的无水乙醇，加入到气液平衡仪中；
（3）准确用移液管量取4ml蒸馏水，加入到气液平衡仪中；
（4）加入适量干燥过的沸石，并加热至沸腾；
（由于时间关系，以上操作步骤由本实验的指导教师完成）；
（4）读取温度计温度后，取①样，然后加入9.7ml蒸馏水，记录时间并分析样品；
（5）40分钟后，第二次读取温度，取②样，加入19.2ml蒸馏水，记录时间分析样品；
（6）40分钟后，第三次读取温度，取③样，加入18.6ml蒸馏水，记录时间分析样品；
（7）40分钟后，第四次读取温度，取④样，加入25ml蒸馏水，记录时间分析样品；
（8）40分钟后，第五次读取温度，取⑤样，停止加热，分析样品；
（9）关闭冷凝水，完成实验。
3、实验记录
原始数据记录：
表1原始数据
实验数据点
x1
y1（质量分数）
y1（摩尔分数）
T/oC
1
78.34
2
0.95
0.9836
0.9590
78.85
3
0.85
0.9572
0.8974
78.75
4
0.70
0.9256
0.8295
78.60
5
0.60
0.9056
0.7896
78.70
6
0.50
0.8689
0.7216
78.92
注：
计算时取值为M（乙醇）=46.07，M（水）=18.015，并近似认为气相为二元体系。
书上给出了乙醇和水的Antione常数：
表2乙醇和水的Antoine常数
Antoine常数
A
B
C
乙醇
18.9119
3803.98
-41.68
水
18.3036
3816.44
-46.13
其中，计算公式为
。
4、数据处理
1.数据处理
由Antione方程
可以计算出乙醇和水在实验温度下的饱和蒸汽压，列于下表。
表3乙醇和水的饱和蒸汽压
温度/oC
饱和蒸汽压/
775.61
339.09
772.55
337.71
767.98
335.65
771.03
337.02
777.76
340.06
772.98
337.90
下面计算组分的活度系数。
假设气相中不含乙二醇，那么就只含有乙醇和水两种组分，活度的计算公式为
其中p=760mmHg
于是可得到各数据点时两组分的活度系数值。
表3乙醇和水的活度系数
——
0.05
0.9892
1.8368
0.15
1.0386
1.5390
0.7
0.3
1.1727
1.2865
0.6
0.4
1.2971
1.1863
0.5
1.4103
1.2444
画出活度系数关于
的变化关系对原始数据进行初步的观察：
图3γ1与γ2的实验点
通过观察明显地发现了一个异常点，即第五组数据存在异常，根据热力学的普遍规律，在本实验中，γ应是关于
的单调函数，而第五组数据不满足这一基本规律，与其他组数据相比较也有较大的差距，故该组数据不可靠，所以在之后的计算中将略去该组数据，关于数据异常的详细分析见文后。
二元体系的Wilson方程：
利用W1式拟合Wilson参数用于乙醇的活度系数计算；
利用W2式拟合Wilson参数用于水的活度系数计算；
利用MATLAB编写程序进行参数拟合得到的结果如下：
（所有程序均附于文后）
表4乙醇和水的Wilson参数拟合结果
0.2775
0.3632
1.1156
0.4354
利用拟合得到的参数可以绘制全浓度范围内的乙醇和水的x-y相图，
y的计算方法为：
温度取实验中的平均温度T=78.76℃，此温度下乙醇和水的饱和蒸汽压分别为：
由此做出相图如下：
图4乙醇-水Wilson方程参数拟合x-y相图
从图中看出，实验数据点与拟合曲线非常接近，这说明了Wilson方程的精确性，但实验数据点较少，拟合出的参数没有太大的参考价值。
关于这张相图之后会再做更详细的讨论。
为了对比乙醇-乙二醇-水拟三元体系与乙醇-水二元体系的差别，利用aspen得到上述温度下的乙醇-水体系的x-y相图，计算方法仍然设定为Wilson方程，与实验数据拟合得到的相图相比较如下：
图5乙醇-乙二醇-水体系与乙醇-水体系比较
从图中可以明显看出，添加了乙二醇的拟三元体系的曲线比单纯的乙醇-水体系的x-y曲线明显的向上偏移了，这说明添加了乙二醇的体系有更大的相对挥发度。
另外，在添加了乙二醇的体系中原有的共沸点消失了，这可能是由于乙二醇对体系的影响，增加了乙醇的相对挥发度，改变了体系内的相互作用，这个体系使得精馏方法得到高纯度乙醇成为了可能，在工业上这一方法也已经得到了应用。
乙二醇在工业上也常用作防冻剂，可以降低某些物质的熔点，这点也似乎能够说明乙二醇有比较活泼的理化性质，能有效的改变一些体系的性质。
当然，在这个拟三元体系中，盐的作用也是不可忽视的，加入某种电解质（在此实验中为KAC）也可以显著提高体系中易挥发组分的相对挥发度。
最后再对实验数据进行热力学一致性检验，采用Herrington面积法进行检验：
计算方法如下：
作图：
图6
图解积分图像
经过分析可知，
利用此式计算出：
由于乙醇-水体系形成沸点为78.60℃的最低恒沸物，因此
℃，而
℃。
得到
所以
，所以根据Herrington检验方法可以判断出这几组数据对热力学一致性的偏差较大。
5、实验结果与讨论总结
实验结论如下：
实验通过参数拟合得到了Wilson方程的两个参数分别为：
并用该参数作出了乙醇-水体系的相图，关于相图已经做过了比较详细的分析。
最后通过热力学一致性检验分析实验数据的质量，但未能达到检验要求。
误差分析：
1、实验室的压力不是一哥标准大气压，但处理数据的过程中将气压视为一个大气压进行计算，可能会有一定的误差；
2、温度计未进行校准，不能保证温度的准确性；
3、实验过程中可能会有部分液体蒸发散失，使计算数据与真实情况有所差距；
4、每次添加水的过程可能存在人为因素引起误差；
5、做气相色谱时可能出现了一些问题，导致最后一组数据异常；
6、实验中只采集了0.5~1的浓度范围内的数据，数据样本不够大；
7、这些原因都不足以说明为何实验会大幅偏离热力学一致性，所以猜测系统误差存在的可能性。
实验小结：
这次实验虽然操作非常简单，没有繁琐的操作过程，但是最终得到的实验数据质量还是不太满意，甚至有一组数据不能使用，当然这样舍弃一组数据是否合理也值得思考，但既然判断出是一个错误的数据就不应该采用，否则结论也是没有意义的。
但这样的错误应该是可以避免的，因为在实验过程中就已经计算出了各组数据的γ值，但并没有留意其变化趋势，导致忽视了这一个明显的异常数据，如果当时及时发现或许可以再补做一组，但现在已经没有机会了。
这是以后实验必须注意的一个问题。
在处理数据的过程中，对比了实验数据所拟合的相图与乙醇-水理论数据计算得出的相图进行了比较，发现了很多不同，为此查阅了相关文献，有幸查阅到了两篇我系数十年前的文献，了解了这一工艺的更多的知识，也是实验的一大收获。
6、参考文献
[1]王保国.化工过程综合实验.清华大学出版社.北京,2004.
[2]朱文涛.物理化学（上册）.清华大学出版社.北京,1999.
[3]高光华.化工热力学.清华大学出版社
[4]段占庭.加盐萃取精馏的研究_用乙二醇加醋酸钾制取无水乙醇[J].石油化工,1980,9（6）:
350-353.
7、附录
1、全部MATLAB计算程序：
函数文件：
%利用LN（γ1）关联式拟合
functionf=fwilson1（beta,x1）
x2=1-x1;
f=-log（x1+beta
（1）.*x2）+x2.*（beta
（1）./（x1+beta
（1）.*x2）-beta
（2）./（x2+beta
（2）.*x1））;
end
------------------------------------------------------------------------------------------------------------------------------------
functionf=fwilson（beta,x1）
f=-log（x2+beta
（2）.*x1）+x1.*（beta
（2）./（x2+beta
（2）*x1）-beta
（1）./（x1+beta
（1）.*x2））;
主脚本程序：
%化工实验5-拟三元体系MATLAB程序%
%第一部分，参数拟合
clearall;
x1=[0.950.850.70.6];
gama1=[0.989181351.0386375591.1727437421.297149698];
%数据由EXCEL计算得到
gama2=[1.8367574951.5389645341.2865347871.18625591];
lngama1=log（gama1）;
lngama2=log（gama2）;
f1=log（gama1）;
f2=log（gama2）;
wbeta1=nlinfit（x1,f1,@fwilson1,[0.6,0.6]）;
wbeta2=nlinfit（x1,f2,@fwilson2,[0,0]）;
%第二部分，利用拟合参数作x-y相图，ASPEN相图数据较多，未在此处列出
p1sat=772.98;
p2sat=334.28;
x11=0:
0.01:
1;
x22=1-x11;
gama11=exp（feval（@fwilson1,wbeta1,x11））;
gama22=exp（feval（@fwilson2,wbeta2,x11））;
p=x11.*gama11*p1sat+x22.*gama22*p2sat;
y11=gama11.*x11*p1sat./p;
y22=gama22.*x22*p2sat./p;
p=x1.*gama1*p1sat+x2.*gama2*p2sat;
y1=gama1.*x1*p1sat./p;
y2=gama2.*x2*p2sat./p;
plot（x11,y11）;
holdon;
plot（x11,y22,'
r'
）;
line（[01],[01]）;
line（[10],[01],'
color'
'
axis（[0101]）;
plot（x1,y1,'
*'
plot（x1,y2,'
xlabel（'
x1'
ylabel（'
y'
holdoff;
pause;
%第三部分，热力学一致性检验及作图
loggama11gama22=log10（gama22./gama11）;
loggama1gama2=log10（gama2./gama1）;
plot（x11,loggama11gama22）;
line（[01],[00]）;
plot（x1,loggama1gama2,'
log（γ1/γ2）’）;
D1=trapz（loggama11gama22）;
D2=trapz（abs（loggama11gama22））;
D=D1/D2*100;
tm=78.60;
seita=tm-100;
J=150*seita/tm;
%结束
升级会员 