| 查看: 1263 | 回复: 3 | ||||
[交流]
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朵）
gx730» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问用MS中CASTEP计算弹性常数，可以加压力和改变温度吗？
已经有6人回复
- 拟合水的密度与温度的函数关系式
已经有4人回复
- 用什么软件可以计算植物在高温或低温胁迫中的半致死温度？具体怎么求？
已经有3人回复
- fluent 中二维模型中，材料属性随着温度变化的udf源代码？
已经有14人回复
- 一定量乙醇加入到密闭容器中，某一温度下压强怎么计算？
已经有7人回复
- pyrosim模拟中封闭的非火源房间的温度怎么测
已经有5人回复
- 把试样放到液氮中，怎么计算多长时间试样能达到液氮温度呢？
已经有4人回复
- 如何在Uspex计算中 设定压力和温度
已经有21人回复
- 高压反应釜中，由于反应温度较高，反应介质会汽化，如何计算反应物初始浓度？
已经有3人回复
- 用MS中的Forcite计算冰，结果发现温度波动异常的大，感觉像做周期性的震荡衰减
已经有6人回复
- 【求助】如何将使用水平集方法计算液滴铺展过程中的温度变化也绘制出来
已经有7人回复
★ ★
月只蓝: 金币+2, 感谢分享研究成果！ 2014-05-30 19:13:37
月只蓝: 金币+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，气相密度，误差，液相密度，误差项就可以了。 |
2楼2014-05-19 20:32:24
3楼2014-05-20 14:41:05
4楼2014-05-20 22:06:11
10