24小时热门版块排行榜    

查看: 1029  |  回复: 3

豆豆sk8

银虫 (小有名气)

[交流] towhee中计算临界温度和密度 已有1人参与

最近又学习了 gx730的文章《学习TOWHEE的经验和感悟》,收获颇丰,对于其中使用towhee自带软件计算  纯物质  临界温度密度进行了简单的学习,并且在这里分享一下。
http://muchong.com/bbs/viewthread.php?tid=6669171&fpage=1
================以下内容说明书里也有============================
1、首先进入到安装目录下,进入/Utils/文件夹,打开终端,输入以下命令:
make fitcoex
该命令为安装fitcoex软件,速度很快,并且生成了名为fitcoex的文件。
2、之后进入到/Examples/VLCC_Fit/里面找到后缀为.sim和.exp的两个文件夹。这两个文件夹就是原始的输入文件,且.sim文件为必须的文件,.exp文件为可选择文件。Examples里面给出的是ethane的例子,因此将这两个文件复制到任意其他位置,并对文件名和内容进行修改。
3、假如我现在要计算butenal,将.sim和.exp文件命名为butenal。
打开butenal.sim文件,该文件包含了命令行,临界指数“beta”的值,以及至少三组不同温度下的towhee_vlcc数据,以及一个特殊的文件结束零线。
例如以下格式(括号为注解):
TraPPE-UA butenal      (名称)         
beta
0.325    (beta值,一般用0.325)
   T       RHOG      ERROR      RHOL       ERROR           (分别为温度K,气相密度,误差,液相密度,误差)
   178.0 .21766E-02 .11795E-03 .55108E+00 .87617E-03 49.0 0.2      (分别对应上面内容,最后多处两个数,不清楚是什么,感觉像是压力)
   197.0 .56274E-02 .51137E-03 .52686E+00 .15355E-02 144.0 11.0   (本人计算时直接将后面两个数字都删掉了)
   217.0 .99620E-02 .31008E-03 .49902E+00 .11414E-02 362.0 12.0
   236.0 .19879E-01 .16763E-02 .46898E+00 .17548E-02 774.0 144.0
   256.0 .30961E-01 .44995E-02 .43228E+00 .52414E-02 1118.0 146.0
   275.0 .55561E-01 .79753E-02 .39554E+00 .52225E-02 1916.0 157.0
   0.00  0.00       0.00       0.00       0.00 0.00 0.00                                 (零线,即代表结束)
----------------------------------------------------------------------------------------------------------------------------------------------------------
.exp文件应该是实验数据文件,因为是可选文件,直接就没有用,其他人使用后可以补充。
4、将之前安装好的fitcoex文件copy到.sim文件所在位置,之后打开终端,输入./fitcoex开始计算。因为每次计算都需要将fitcoex文件copy过来比较麻烦,这里尝试了将fitcoex文件放入到bin文件下作为整个系统的命令,但是并不能被识别,本人对linux不太了解,还望高手赐教。
开始计算后,会出现几个问题:
如让输入文件名,我是直接输入的.sim文件的名字,即butenal;
处理气体压力数据?  因为之前把貌似为压力的数据都删掉了,这里我选的F;
输入最小温度点:178
是否希望转换格式?F     这里所有的同意都用T,不同意用F。而不可以用Y/N
然后会出现计算公式,不同的计算公式计算结果不同,自己可以进行选择。

其他小软件还有待学习,学习后会继续进行补充,也希望各位大神能够分享自己的经验。   再次感谢gx730起到的带头作用。
其文章链接:http://muchong.com/bbs/viewthread.php?tid=6669171&fpage=1
回复此楼

» 收录本帖的淘帖专辑推荐

towhee的运用讨论!!

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

Justskate
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

豆豆sk8

银虫 (小有名气)

★ ★
月只蓝: 金币+2, 感谢分享研究成果! 2014-05-30 19:13:37
再补充一下,gx730 同学建议使用原始乙烷文件对该小程序进行测试,测试内容为:1、包含ethane.sim  ethane.exp 文件,且ethane.sim文件中数据部分是完整的,即包括了上面所说的不清楚其意义的俩数;2、包含ethane.sim  ethane.exp 文件,但去掉了上面所说的不清楚其意义的俩数;3、只包含ethane.sim,且ethane.sim文件中数据部分是完整的,即包括了上面所说的不清楚其意义的俩数;4、只包含ethane.sim,但去掉了上面所说的不清楚其意义的俩数;    最后的计算结果四组完全相同,除了出现过一点不同,即第四组第一次计算的时候出的数略有偏差,又算了第二次就跟其他的一样了,没找到原因。  
另外该程序提供了三种计算临界温度和临界密度的公式,第一种比较简单,后两种里面需要输入计算的点数,这个往大了填吧,比如400、500,因为很快就能计算完,如果数字较小可能尚未达到最后的结果。
总之,目前来看只需要.sim文件,且该文件只需要    温度K,气相密度,误差,液相密度,误差项就可以了。
Justskate
2楼2014-05-19 20:32:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gx730

金虫 (正式写手)

送红花一朵
送红花  
3楼2014-05-20 14:41:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

豆豆sk8

银虫 (小有名气)

引用回帖:
3楼: Originally posted by gx730 at 2014-05-20 14:41:05
送红花  

谢谢啦~
Justskate
4楼2014-05-20 22:06:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 豆豆sk8 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见