东东（完）

diyifan001

按照上述思路在matlab和C#中进行了实际操作，下面就用计算结果检验一下（由于c#中的未完成  故全部用matlab中的程序运行）。
1、泡点检验：
   甲烷、乙烷、丙烷组成的三元系统检验
  随便捏一个组成和压力  对比程序求出的泡点和hysys求出的泡点（压力   5000Kpa  组成：甲烷 ：0.7，乙烷：0.2，丙烷： 0.1）
  wz=[ 4599000 190.564 0.0115478 0.7;
     4872000  305.32  0.099493 0.2;
     4248000  369.83 0.15229 0.1];
Kij=[0 -0.0026 0.014;
     -0.0026 0 0.0011;
     0.014 0.0011  0];
p=5000000;
[Bt,yi,Ki]=findBT(p,wz,Kij)

其中  findBT函数是根据上述思路编的。wz是物质临界压力、临界温度、偏心因子、在系统中的摩尔组成形成的二维矩阵Kij是二元交互系数形成的矩阵，p是压力：5000Kpa。函数返回泡点温度（Bt）、与泡点温度对应的气相组成（yi），以及相平衡常数。计算结果如下：
>> paodian
Bt =
  213.5819

yi =
    0.9294
    0.0601
    0.0105

Ki =
    1.3278
    0.3004
    0.1047
>>
hysys计算截图：
组成：

计算结果：

hysys算出的泡点温度为：213.3  程序算出的为213.58，结果可以接受。
在甲烷、乙烷、丙烷、氮气组成的体系中，在很宽的范围内验证，结果和hysys的差别不超过2℃。

沙发进行PT闪蒸的验证。

diyifan001

这次加入氮气 程序如下：
wz=[ 4599000 190.564 0.0115478 0.5;
     4872000  305.32  0.099493 0.2;
     4248000  369.83 0.15229 0.1;
     3394000  126.2  0.04  0.2];
Kij=[0 -0.0026 0.014  0.0311;
     -0.0026 0 0.0011  0.0515;
     0.014 0.0011  0  0.0852;
     0.0311  0.0515  0.0852 0];
p=4000000;t=200;
[e,xi,yi,ki]=PTFlash(p,t,wz,Kij)

和上楼一样，只是PTFlash是根据前述思路编的函数，函数返回 气化率(e)、平衡时气象（yi）和液相（xi）组成及平衡时各组分的平衡常数。计算结果如下：

>> test
e =
    0.4905

xi =
    0.4059
    0.3462
    0.1915
    0.0564

yi =
    0.5977
    0.0481
    0.0050
    0.3492

ki =
    1.4732
    0.1391
    0.0259
    6.1922
>>


同样条件下hysys计算截图：









从气化率、相平衡常数，及平衡时气液相组成看，结果那是相当的符合。
看来教课书上的东西还有实际作用。

下一楼通过编的计算泡点的函数，利用matlab画图的功能  做出甲烷的饱和蒸汽压曲线。

diyifan001

单质的泡点和露点及为温度及为其在某个压力下的饱和温度，根据这个思路，由于编findBT函数时考虑的通用性，所以一个单质也可以计算，给出一串压力，求出一串泡点，以泡点温度为横坐标，压力为纵坐标，就可作出甲烷的饱和蒸汽压线
程序如下：
wz=[4599000 190.564 0.0115478 1];
Kij=[0];
p=100000:200000:4100000
for i=1:21
    t(i)=findBT(p(i),wz,Kij);
end
t
plot(t,p);

结果如下
>> jiawantu

p =

  Columns 1 through 8

      100000      300000      500000      700000      900000     1100000     1300000     1500000

  Columns 9 through 16

     1700000     1900000     2100000     2300000     2500000     2700000     2900000     3100000

  Columns 17 through 21

     3300000     3500000     3700000     3900000     4100000


t =

  Columns 1 through 10

  111.4158  126.6213  135.2480  141.5965  146.7320  151.1098  154.9105  158.3410  161.4625  164.3049

  Columns 11 through 20

  166.9501  169.4155  171.7609  173.9843  176.0362  178.0204  179.9426  181.7324  183.4698  185.1472

  Column 21

  186.7626

>>

matlab画出的图如下：


横坐标为温度   K
纵坐标为压力  pa


连个向量画图，在matlab里太简单了。

diyifan001

小小的收获：
1、虽然这些东西很幼稚，但通过亲身实践，终于知道了数值计算的复杂。看来没有数学的发展，今天的化工计算就不可想象有多复杂。
2、终于浅显知道了hysys物流框里下面的条条从黄到绿的变化过程中，开发人员用了多少心血。
3、文献里介绍的方法有局限，比如按照书上的方法  就算自己根据理解加一些制约手段，但很多点上很难收敛于正确数值。这一点就很棘手，加上很多设备放在一起模拟，贯续法、EO法等很陌生很可怕的数学方面的东西，模拟软件卖那么多钱还真是有道理地。
4、个人觉得hysys做的比Aspen plus好。
5、写文献的人太不负责任，一个状态方程，抄都要抄错，真是误导人。
6、学校里的日子，白混了。

Horse

我回复过的，但是没看懂过

diyifan001

饱和蒸汽压总看得懂呀
晕你

Horse

饱和蒸汽压倒是懂啊，，编程的这一块都是不懂得

lidaxue1987

建议楼主把完整的程序给出一下，呵呵  仅仅是建议

weihuawind

完整的matlab 程序有吗 可以传上来嘛 谢谢啦

chenhaijunsjy

感谢版主的资料

浅露

楼主威武，受教

wklbuct

强大，不懂编程的路过