找回密码
 注册
搜索
热搜: 超星 读书 找书
楼主: 醉乡常客

[【理工类】] (无奖求助)二元非线性方程求根问题

[复制链接]
发表于 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,我帖的图的上半截就是书上的解题过程。)
回复

使用道具 举报

发表于 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方程是分段函数,且无甲乙酮数据。
回复

使用道具 举报

发表于 2007-7-24 10:39:25 | 显示全部楼层
醉兄,我从移动硬盘里找到软件了,需要哪个,我给你传^_^

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

发表于 2007-7-25 00:51:55 | 显示全部楼层
醉兄,我试算了你的题目,几乎可以肯定,题目给出的数据不是他原来计算时用的数据。

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

我取 gamma_1 = 1.031,gamma_2 = 1.667,得到 Lambda12 = 1.0818,Lambda21 = 0.3778,用此数据和你所给的安东尼常数计算得到:

x1  y1  y1*  t (°C)  t* (°C)
0.076  0.193  0.188  75.3  75.5
0.197  0.377  0.381  70.7  71.1
0.356  0.528  0.541  67.5  67.6
0.498  0.622  0.641  65.9  65.7
0.622  0.695  0.714  65.1  64.6
0.747  0.777  0.784  64.4  64.0
0.844  0.832  0.846  64.3  63.8
0.936  0.926  0.922  64.4  64.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.161737  1.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


回复

使用道具 举报

发表于 2007-7-25 10:43:07 | 显示全部楼层
不要用matlab那个巨无霸,用Scilab吧免费的。
回复

使用道具 举报

发表于 2007-9-22 23:31:37 | 显示全部楼层
难得一见的理工讨论贴,可为论坛理工讨论的模板,经典!

不知醉版后事如何?
回复

使用道具 举报

 楼主| 发表于 2007-9-23 00:07:57 | 显示全部楼层
引用第49楼brook于2007-09-22 23:31发表的 :
难得一见的理工讨论贴,可为论坛理工讨论的模板,经典!

不知醉版后事如何?

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

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

使用道具 举报

发表于 2007-9-24 13:57:56 | 显示全部楼层
额 做了个验算的工具
怎么解我真的忘了 - -+
这里高手真多

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|网上读书园地

GMT+8, 2024-11-14 00:08 , Processed in 0.256848 second(s), 4 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表