24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

relonfbq

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1294.3
散金: 1140
紅花: 2
帖子: 301
在線: 99.7小時
蟲號: 744333
注冊: 2009-04-09
專業(yè): 金屬材料的微觀結(jié)構(gòu)

[求助] 我這樣的程序到底對不對

我的程序是做分子動力學(xué)的，我編程的思想是這樣的，用的p-r恒壓模型，積分算法是gear算法。
我現(xiàn)在通過一個算例來驗(yàn)證我的程序的正確性（首先我查不出我程序的錯誤了）
我引入了5*5*5的249個原子的體系（完美晶體應(yīng)該是250個原子），設(shè)置初始溫度0K，外壓保持0pa，然后運(yùn)行，運(yùn)行結(jié)果是，前面28000步（步長是1fs）體系的溫度在50K以下震蕩（或說漲落），因而動能也在很小的數(shù)值下漲落，同樣，在前面的28000步長里勢能，和總能也在約百分之一的范圍里漲落。（注意我的數(shù)據(jù)是200步一平均）。
但是在約30000步左右，體系的溫度開始逐步向上升，因而引起動能和總能量逐漸上升（勢能和總能的絕對值下降），我運(yùn)行了40000步，沒有看見停止上升的趨勢。
當(dāng)然體系的體積卻一直在百分之10的范圍內(nèi)漲落！
我這樣的結(jié)果有什么問題呢?而最關(guān)鍵的我的問題（我的程序）出在哪里呢？？？？

[ Last edited by relonfbq on 2011-4-26 at 11:13 ]

回復(fù)此樓

» 猜你喜歡

調(diào)劑310 已經(jīng)有11人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有9人回復(fù)
296求調(diào)劑已經(jīng)有3人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有3人回復(fù)
299求調(diào)劑已經(jīng)有9人回復(fù)
22408 359分調(diào)劑已經(jīng)有6人回復(fù)
085600 材料與化工 329分求調(diào)劑已經(jīng)有12人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有14人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有11人回復(fù)
總分293求調(diào)劑已經(jīng)有9人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2011-04-22 22:25:49

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 relonfbq 的主題更新