944f4fe875fb39d 发表于 2007-7-23 20:16:41

那么,你原来给出的 1.0082 和 1.6719 是什么?

gamma_1 = 1.031

gamma_2 = 1.667

把这两个数据代进去试试,看看对不对?

醉乡常客 发表于 2007-7-23 20:33:01

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,我帖的图的上半截就是书上的解题过程。)

944f4fe875fb39d 发表于 2007-7-23 21:14:57

按1.013和1.667,得到 Lambda12=1.08182135,Lambda21=0.377766147

是不是 1.0818、0.3778?用它们计算 t 和 y1 结果是否更好些?

另外,你再查一下安东尼常数给得对不对。

有时候,人们在处理实验数据时用一套数据得到了较好的结果,写书时就忘了。

醉乡常客 发表于 2007-7-23 22:22:29

浙大陈钟秀教材上的甲醇和甲乙酮数据和例题一致。

天大马沛生教材上Antoine方程是分段函数,且无甲乙酮数据。

louselice^_^ 发表于 2007-7-24 10:39:25

醉兄,我从移动硬盘里找到软件了,需要哪个,我给你传^_^

944f4fe875fb39d 发表于 2007-7-25 00:51:55

醉兄,我试算了你的题目,几乎可以肯定,题目给出的数据不是他原来计算时用的数据。

他原来的计算结果是合理的。

我取 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 无正确的解。

按最可能犯的错误判断,他们当初计算该例题时用的是另外的安东尼常数,或者是单位搞错了,或者是单位转换时什么地方搞混淆了……

醉乡常客 发表于 2007-7-25 00:56:52

我的答案

x1
T
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


醉乡常客 发表于 2007-7-25 01:16:53


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


ereree 发表于 2007-7-25 10:43:07

不要用matlab那个巨无霸,用Scilab吧免费的。

brook 发表于 2007-9-22 23:31:37

难得一见的理工讨论贴,可为论坛理工讨论的模板,经典!

不知醉版后事如何?

醉乡常客 发表于 2007-9-23 00:07:57

引用第49楼brook于2007-09-22 23:31发表的 :
难得一见的理工讨论贴,可为论坛理工讨论的模板,经典!

不知醉版后事如何?

我的结论和bookish老师说的一样。该书此题未验算!

当然,我不敢肯定。毕竟这本书算得上是一本非常优秀的热力学教材。

ly188 发表于 2007-9-24 13:57:56

额 做了个验算的工具
怎么解我真的忘了 - -+
这里高手真多
页: 1 2 [3]
查看完整版本: (无奖求助)二元非线性方程求根问题