既然你需要的只是数据(液相0.37mol分率的甲醛),那我干脆计算后列出来给你:
Temperature [C] Pressure [bar]
10.0000000000000 1.30563074383307
20.0000000000000 1.82419464454581
30.0000000000000 2.48804937239319
40.0000000000000 3.32154129280952
50.0000000000000 4.35048920012608
60.0000000000000 5.60197455987817
70.0000000000000 7.10454940576685
80.0000000000000 8.88834823071947
90.0000000000000 10.9851133979955
100.000000000000 13.4289104256675
本来以为是一道很简单的题目,用aspen plus和我的程序算没有什么疑问。
选用nrtl和wilson方程,发现没有活度系数参数(我的数据来自于a+ 10),刚才用 a+ 2006,都没有发现此参数,莫非此体系可以认为是理想溶液?因此,活度系数常为1,我11楼的计算则可以简化为
y_i*P=x_i*Ps_i,
此方程的解过程则很简单了,只要计算出相应温度下的水和甲醛的Ps_i值即可。
问题是我在用hysys时,发现选取nrtl,其值相当反常,50度时的压力只有0.13bar左右,与4.3bar的值相差也太远了一点,换了好几种物性包,就extended nrtl 和general nrtl方程的解与a+的结果接近。
怀疑hysys的nrtl数据库有问题,由于我手头查不到甲醛-水体系的相平衡参数,如果有哪位网友有此数据,麻烦对比验证一下,看看究竟hysys和a+的数据究竟谁正确。
由于我的程序来自于a+,结果与a+相同。
我更倾向于a+的计算结果,因此以上所列数据为hysys 用extended nrtl方程时的结果。 |