24小时热门版块排行榜    

查看: 1112  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 一志愿武理材料工程348求调剂 +5  ̄^ ̄゜汗 2026-03-19 7/350 2026-03-22 19:44 by 公瑾逍遥
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
信息提示
请填处理意见