文档介绍:气液平衡-实验报告
图1 平衡法测定气液平衡原理图
当达到平衡时,除了两相的压力和温度分别相等外,每一组分的化学位也相等,即逸度相等,其热力学基本关系为:
常压下,气相可视为理想气体x=,修正值:)
经过计算得到液体平衡时实际的温度如表3 所示:
表3 气液平衡时实际温度
t主(℃)
t辅(℃)
t修正(℃)
t校正(℃)
t实际(℃)
苯和正庚烷纯组分在本实验不同温度下的饱和蒸汽压和活度系数的计算
(1)苯和正庚烷纯组分在本实验不同温度下的饱和蒸汽压的计算
由Antoine(安托尼)公式:
lgPi0=Ai–Bi/(Ci+t)
式中:
t—温度,℃(即计算所得的实际温度) P0—饱和蒸汽压,mmHg
所以可得到:
Pi0=10(Ai-Bi/(Ci+t))
纯物料的Antoine常数见表4
表4:纯物料的安托尼常数
组分
A
B
C
t范围(℃)
苯
15~84
正庚烷
3~127
℃为例:
苯的饱和蒸气压p10=10(-÷(+))=
正庚烷的饱和蒸气压P20=10(-÷(+)=
计算五组纯组分饱和蒸汽计算结果见表5
表5:纯组分饱和蒸汽在实验温度下计算结果
组号
1
2
3
4
5
P10/ mmHg
P20/ mmHg
(2)苯和正庚烷纯组分在本实验不同温度下的活度系数的计算
由实验测得等压下气液平衡数据,则可用下式计算不同组成下的活度系数:
℃为例计算:
表6 苯和正庚烷活度系数计算表
t实际(℃)
苯的液相组成(%)
苯的气相组成(%)
苯的饱和蒸气压 (mmHg)
正庚烷饱和蒸气压 (mmHg)
苯的活度系数
正庚烷的活度系数
用非线形最小二乘法回归配偶参数Λ12 、Λ21 并求液相组成
活度系数和组成关系采用Wilson方程关联
lnγ1=-ln
lnγ2=-ln
目标函数选为气相组成误差的平方和,即
F=(y1实-y1计)2j+(y2实-y2计)2j
用非线性最小二乘法拟合:matlab拟合程序见下:
function F=li(bb)
x1=[ ]
x2=[ ]
y1=[ ]
y2=[ ]
p10=[ ]
p20=[ ]
fori=1:4 F(i)=y1(i)-p10(i)*x1(i)/760*exp(-log(x1(i)+bb(1)*x2(i))+x2(i)*(bb(1)/(x1(i) +bb(1)*