|
|
是的,就是用了这个公式。
现在的问题是,用了perry 手册里的实验数据和用兰式化学手册里的abc数据用公式计算出来结果吻合度好一些。但是用crc物理化学手册里的abc参数计算则相差极大。
问题是,我只是随便选取了Mg进行了计算,我觉得即使用crc数据也不应差别这么大啊。不知道我哪里弄错了。下面是matlab程序。
%计算物质的饱和蒸气压,参数可以从CRC物理化学手册及兰式化学手册及fundimentals of vacuum technology
%上查。另外X轴要取对数坐标。兰式化学手册+perry化学手册数据相对较吻合。
clear;
t = 1:50:2400;
A =8.489;%前面是crc参数,12.993兰式化学参数;
B = -7813;%前面是crc参数,13579.8兰式化学参数;
C = -0.8253;%crc有c和d两个参数,兰式化学没有。
D = 0;
y = exp(2.881+A+B./(t+273)+C*log(t+273));%crc手册给出的公式,但是是torr
%y = exp((A-B./(t+273)))*51.715;% 兰式化学给出的公式,自己换算成了torr。
x=[1 5 10 20 40 60 100 200 400 760; 621 702 743 789 838 868 909 967 1034 1107];% perry手册里面的实验数据。
plot(y, t,'-',x(1, , x(2, ,'o');
%plot()
set(gca,'XScale','log');
那位大侠帮我研究研究? |
|