关于多元混合液相饱和蒸汽压的问题

29楼，我的最终结果也是8.2m3, 你看到的8.23，是第一次B的结果，迭代第二次是，结果已经是4421kg，即8.206m3（不好意思，原扫描件未列体积只是列了质量）了，应该说我们的结果相同。
我的数学功底不好，直接解方程组，我经常出错，后来有了excel，发现用迭代法对于我这种数学功底差的人来说，是个偷懒的好办法，所以能用迭代我就不再用直接解方程的路子了。

另，我是用打印机扫描的，相机拍摄以前试过，效果不好，好象理光出了一款相机，有个扫描模式，排出来的效果据说与扫描相差无几，没有见过。

[ 本帖最后由 csldg 于 2009-1-19 19:46 编辑 ]

楼上的朋友你果然是“数学功底差”啊，我说怎么你的思路方法和我一样怎么结果差这么多，你计算N2分压的时候算错了，1.5*50/41.77你怎么算出来2.7？

晕哦。哈哈

郁闷啊，好不容易有个不用程序计算的地方，还犯了个低级错误。
本来想把帖子修改过来，不过一想还是不改了，放在这里当个反面教材吧 - 越是简单的地方越要注意不要犯低级错误，做技术的不能有丝毫马虎！

各位高人，没想到这个帖子能引起热议，挺感谢大家给予的帮助，让我对这个问题的理论计算有了深入的认识。
不过有些问题还是想澄清一下。这个问题的原形是3000立方的LPG球罐，组成我也列出来吧：丙烯（0.07％），丙烷（0.08％），异丁烯（17.14％），正丁烯（14.45％），反丁烯（18.13％），顺丁烯（12.73），异丁烷（26.66％），正丁烷（8.13％），异戊烯（1.23％），异戊烷（1.34％），正戊烷（0.03％）。
大家应该知道LPG球罐正常情况下存储是不需要氮封的，所以气相空间仅为LPG蒸汽，球罐的装填系数可以取到0.8，所以可以知道气相空间的体积。但是我认为在平衡状态下（就是前面说得液相汽化速率等于气相冷凝速率的时候），只要球罐内有液相存在，气相的压力应该是一定的，就是说，是3000立的液化气球罐也好，家用的液化气钢瓶也好，只要组成相同，环境温度相同，他们内部的气相压力应该相等（前提是都有液相存在），与气相空间的体积无关。不知道我的这一观点是否正确，请各位高人指正。

楼上的，我认为严格说起来，你上面的观点是错误的。汽液平衡时，Ps=f(T, Vf)，对于固定的罐容积，汽液相的体积和为常数，但汽液相的体积却是可以变化的，即，Vf不是固定值，随汽相体积变化而变划，那么相应的，Ps肯定也要变。

但实际中生产时，我猜应该可以观察到，排除温度的影响，压力接近常数，原因是：
实际操作时，由于汽相密度远小于液相密度，也就是说，汽相的体积变化对分率Vf的影响不大，所以Ps变化也不大，除非罐快要空，此时液相的量也很少，或者罐几乎满（当然实际操作时是不能这么操作），不大的汽相变化也可以引起较明显的Vf变化，此时压力变化会比较明显。

[ 本帖最后由 dafugege 于 2009-1-23 09:51 编辑 ]

看到这篇帖子，特别是zjmche 、csidg和dafugege的辩论受益匪浅。 这个问题我觉得用aspen静态模拟也可以计算。
具体步骤：
（1）、首先计算50m3罐内N2在25°，0.5BarG下的质量为85.543kg（可用aspen模拟，亦可手算）。
（2）、把N2和5t烃类都作为进料，用闪蒸或换热器模块都可以，设定温度为25°，通过改变压力设定值，使得当出口物料的体积为50m3时的压力值即为罐内的平衡压力。
（3）、计算的压力为6.825barG。罐内液体体积为8.22m3。
状态方程选用的是SRK。

这是一个有些意思的话题，对模拟软件没有一定深度的理解，并不是很容易将这个题作出来。
巧用模拟软件的控制器就可以轻松搞定了,同意35楼的观点，但是由此想到汽液平衡实验时加入的物料的量和平衡釜的体积严格上讲都会影响到实验结果。

[ 本帖最后由 lty 于 2009-7-13 08:49 编辑 ]

这是个汽液平衡的问题。按汽液项都是理想溶液来说明问题。假设有n个组分，那么罐内压力P=p1+p2+p3+...+pn;pi=pi0*xi;pi0=f(ti);可以用安托因公式计算纯组分的饱和蒸汽压，用拉乌尔定律计算各组分的汽相分压（因为前面假设是理想溶液），一定温度下总压就可以计算出来。因为pi0和ti的关系一般是非线性关系，所以温度 t一定，LPG组分不同，pi会不一样，从而P也不一样；反之，LPG组分一定，不管进料状态如何，一定的温度下，总有一定的总压，因为pi是唯一确定的。汽相组分yi=pi/P。实际计算时，系统可能偏离理想溶液很大，这时压力要用逸度或活度来代替。根据我的经验，汽相用状态方程法求逸度，液相用活度系数法求活度，往往能得到较好的结果，当然这要用计算机求解。

你可以通过ASPEN软件完成。