最近又学习了 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
再补充一下,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,气相密度,误差,液相密度,误差项就可以了。