gamma_1 = 1.031
gamma_2 = 1.667
把这两个数据代进去试试,看看对不对? 1.0082和1.6719是我自己用Antoine方程算出的活度系数。
教材上给出的1.013和1.667代回去不对。
Antoine方程(压力单位是bar)算得的64.3℃饱和蒸汽压
甲醇 1.004817bar
甲基乙基酮 0.609528bar
我觉得书上给出的0.982×10^5、0.6067×10^5两个数有问题,此后的活度系数也就有点问题。
(这是陈钟秀书上的例7-3,我帖的图的上半截就是书上的解题过程。) 按1.013和1.667,得到 Lambda12=1.08182135,Lambda21=0.377766147
是不是 1.0818、0.3778?用它们计算 t 和 y1 结果是否更好些?
另外,你再查一下安东尼常数给得对不对。
有时候,人们在处理实验数据时用一套数据得到了较好的结果,写书时就忘了。 浙大陈钟秀教材上的甲醇和甲乙酮数据和例题一致。
天大马沛生教材上Antoine方程是分段函数,且无甲乙酮数据。 醉兄,我从移动硬盘里找到软件了,需要哪个,我给你传^_^ 醉兄,我试算了你的题目,几乎可以肯定,题目给出的数据不是他原来计算时用的数据。
他原来的计算结果是合理的。
我取 gamma_1 = 1.031,gamma_2 = 1.667,得到 Lambda12 = 1.0818,Lambda21 = 0.3778,用此数据和你所给的安东尼常数计算得到:
x1y1y1*t (°C)t* (°C)
0.0760.1930.18875.375.5
0.1970.3770.38170.771.1
0.3560.5280.54167.567.6
0.4980.6220.64165.965.7
0.6220.6950.71465.164.6
0.7470.7770.78464.464.0
0.8440.8320.84664.363.8
0.9360.9260.92264.464.0
*号为计算值。
我觉得问题在安东尼常数和恒沸点的x1的不一致上,按该组常数和恒沸点的x1,使得Lambda12 和 Lambda21 无正确的解。
按最可能犯的错误判断,他们当初计算该例题时用的是另外的安东尼常数,或者是单位搞错了,或者是单位转换时什么地方搞混淆了……
我的答案
x1T
y1
y1
实验值
温度
实验值
0
352.75
0.00000
0.05
349.5
0.14303
0.076
348.12
0.20220
0.193
348.45
0.1
346.94
0.24963
0.197
343.37
0.39687
0.337
343.85
0.2
343.3
0.40086
0.3
340.83
0.50452
0.356
339.78
0.54999
0.528
340.65
0.4
339.1
0.58165
0.498
337.95
0.64329
0.622
339.05
0.5
337.93
0.64442
0.6
337.13
0.69859
0.622
337
0.71004
0.695
338.25
0.7
336.65
0.75052
0.747
336.49
0.77484
0.777
337.55
0.8
336.42
0.80538
0.829
336.41
0.82330
0.932
337.45
0.842
336.43
0.83230
0.9
336.59
0.87696
0.936
336.8
0.91129
0.926
337.55
0.95
336.95
0.92798
1
337.65
0.99978
0.1617371.26246
x1
T
y1
y1实验值
温度实验值
0
352.75
0
0.05
346.27
0.222321
0.076
344.73
0.273694
0.193
348.45
0.1
343.73
0.30786
0.197
341.39
0.395499
0.337
343.85
0.2
341.34
0.397617
0.3
340
0.46128
0.356
339.44
0.494634
0.528
340.65
0.4
339.07
0.520999
0.498
338.39
0.581303
0.622
339.05
0.5
338.38
0.582616
0.6
337.88
0.648874
0.622
337.8
0.664437
0.695
338.25
0.7
337.57
0.722145
0.747
337.49
0.75959
0.777
337.55
0.8
337.43
0.803842
0.829
337.42
0.829335
0.932
337.45
0.842
337.42
0.84105
0.9
337.45
0.895483
0.936
337.5
0.931404
0.926
337.55
0.95
337.55
0.946743
1
337.65
0.999775
不要用matlab那个巨无霸,用Scilab吧免费的。 难得一见的理工讨论贴,可为论坛理工讨论的模板,经典!
不知醉版后事如何? 引用第49楼brook于2007-09-22 23:31发表的 :
难得一见的理工讨论贴,可为论坛理工讨论的模板,经典!
不知醉版后事如何?
我的结论和bookish老师说的一样。该书此题未验算!
当然,我不敢肯定。毕竟这本书算得上是一本非常优秀的热力学教材。 额 做了个验算的工具
怎么解我真的忘了 - -+
这里高手真多