| 24小時(shí)熱門(mén)版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 1106 | 回復(fù): 3 | ||||
[交流]
towhee中計(jì)算臨界溫度和密度 已有1人參與
|
|
最近又學(xué)習(xí)了 gx730的文章《學(xué)習(xí)TOWHEE的經(jīng)驗(yàn)和感悟》,收獲頗豐,對(duì)于其中使用towhee自帶軟件計(jì)算 純物質(zhì) 臨界溫度密度進(jìn)行了簡(jiǎn)單的學(xué)習(xí),并且在這里分享一下。 http://www.gaoyang168.com/bbs/viewthread.php?tid=6669171&fpage=1 ================以下內(nèi)容說(shuō)明書(shū)里也有============================ 1、首先進(jìn)入到安裝目錄下,進(jìn)入/Utils/文件夾,打開(kāi)終端,輸入以下命令: make fitcoex 該命令為安裝fitcoex軟件,速度很快,并且生成了名為fitcoex的文件。 2、之后進(jìn)入到/Examples/VLCC_Fit/里面找到后綴為.sim和.exp的兩個(gè)文件夾。這兩個(gè)文件夾就是原始的輸入文件,且.sim文件為必須的文件,.exp文件為可選擇文件。Examples里面給出的是ethane的例子,因此將這兩個(gè)文件復(fù)制到任意其他位置,并對(duì)文件名和內(nèi)容進(jìn)行修改。 3、假如我現(xiàn)在要計(jì)算butenal,將.sim和.exp文件命名為butenal。 打開(kāi)butenal.sim文件,該文件包含了命令行,臨界指數(shù)“beta”的值,以及至少三組不同溫度下的towhee_vlcc數(shù)據(jù),以及一個(gè)特殊的文件結(jié)束零線。 例如以下格式(括號(hào)為注解): 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 (分別對(duì)應(yīng)上面內(nèi)容,最后多處兩個(gè)數(shù),不清楚是什么,感覺(jué)像是壓力) 197.0 .56274E-02 .51137E-03 .52686E+00 .15355E-02 144.0 11.0 (本人計(jì)算時(shí)直接將后面兩個(gè)數(shù)字都刪掉了) 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 (零線,即代表結(jié)束) ---------------------------------------------------------------------------------------------------------------------------------------------------------- .exp文件應(yīng)該是實(shí)驗(yàn)數(shù)據(jù)文件,因?yàn)槭强蛇x文件,直接就沒(méi)有用,其他人使用后可以補(bǔ)充。 4、將之前安裝好的fitcoex文件copy到.sim文件所在位置,之后打開(kāi)終端,輸入./fitcoex開(kāi)始計(jì)算。因?yàn)槊看斡?jì)算都需要將fitcoex文件copy過(guò)來(lái)比較麻煩,這里嘗試了將fitcoex文件放入到bin文件下作為整個(gè)系統(tǒng)的命令,但是并不能被識(shí)別,本人對(duì)linux不太了解,還望高手賜教。 開(kāi)始計(jì)算后,會(huì)出現(xiàn)幾個(gè)問(wèn)題: 如讓輸入文件名,我是直接輸入的.sim文件的名字,即butenal; 處理氣體壓力數(shù)據(jù)? 因?yàn)橹鞍衙菜茷閴毫Φ臄?shù)據(jù)都刪掉了,這里我選的F; 輸入最小溫度點(diǎn):178 是否希望轉(zhuǎn)換格式?F 這里所有的同意都用T,不同意用F。而不可以用Y/N 然后會(huì)出現(xiàn)計(jì)算公式,不同的計(jì)算公式計(jì)算結(jié)果不同,自己可以進(jìn)行選擇。 其他小軟件還有待學(xué)習(xí),學(xué)習(xí)后會(huì)繼續(xù)進(jìn)行補(bǔ)充,也希望各位大神能夠分享自己的經(jīng)驗(yàn)。 再次感謝gx730起到的帶頭作用。 其文章鏈接:http://www.gaoyang168.com/bbs/viewthread.php?tid=6669171&fpage=1 |
towhee的運(yùn)用討論!! |

|
再補(bǔ)充一下,gx730 同學(xué)建議使用原始乙烷文件對(duì)該小程序進(jìn)行測(cè)試,測(cè)試內(nèi)容為:1、包含ethane.sim ethane.exp 文件,且ethane.sim文件中數(shù)據(jù)部分是完整的,即包括了上面所說(shuō)的不清楚其意義的倆數(shù);2、包含ethane.sim ethane.exp 文件,但去掉了上面所說(shuō)的不清楚其意義的倆數(shù);3、只包含ethane.sim,且ethane.sim文件中數(shù)據(jù)部分是完整的,即包括了上面所說(shuō)的不清楚其意義的倆數(shù);4、只包含ethane.sim,但去掉了上面所說(shuō)的不清楚其意義的倆數(shù); 最后的計(jì)算結(jié)果四組完全相同,除了出現(xiàn)過(guò)一點(diǎn)不同,即第四組第一次計(jì)算的時(shí)候出的數(shù)略有偏差,又算了第二次就跟其他的一樣了,沒(méi)找到原因。 另外該程序提供了三種計(jì)算臨界溫度和臨界密度的公式,第一種比較簡(jiǎn)單,后兩種里面需要輸入計(jì)算的點(diǎn)數(shù),這個(gè)往大了填吧,比如400、500,因?yàn)楹芸炀湍苡?jì)算完,如果數(shù)字較小可能尚未達(dá)到最后的結(jié)果。 總之,目前來(lái)看只需要.sim文件,且該文件只需要 溫度K,氣相密度,誤差,液相密度,誤差項(xiàng)就可以了。 |


| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|