查看: 1283  |  回复: 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的回帖

gx730

金虫 (正式写手)

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

豆豆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的回帖

豆豆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的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +6 Isak0320 2026-10-12 9/450 2026-10-12 21:13 by Isak0320
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见