 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1720 | 回復(fù): 3 | |||
/name-?=p/銅蟲 (初入文壇)
|
[交流]
lammps在npt平衡后換成nve溫度一直升高是怎么回事？ 已有2人參與
|
|
各位大神請(qǐng)教一個(gè)問題 用lammps模擬計(jì)算3C-SiC熱導(dǎo)率,時(shí)間步長(zhǎng)1fs 在1000K下，先在NVT下模擬10萬步，使系統(tǒng)達(dá)到平衡然后換上NVE系綜后在運(yùn)行125萬步。但是當(dāng)換上NVE系綜不久后，溫度就開始升高，到結(jié)束時(shí)溫度都達(dá)到幾千甚至上萬K，這么高溫度碳化硅固體早就沒了。125萬步不是設(shè)定步數(shù)，125萬步時(shí)原子丟失了，在低溫時(shí) 溫度也會(huì)升高 想問這是哪里出問題了 這是in文件 晶格常數(shù)設(shè)置的讀坐標(biāo)文件 盒子邊長(zhǎng)大小是十個(gè)晶格常數(shù)也就是一共8000個(gè)原子 variable T equal 1000 # 溫度為 400 K #variable A equal 4.348 # 晶格常數(shù)大概為 6.35A variable DT equal 0.001 # 積分步長(zhǎng)為 0.3 fs # 初始化坐標(biāo)、勢(shì)函 數(shù)、速度等數(shù)據(jù) boundary p p p read_data SiC-beta.txt pair_style tersoff pair_coeff * * SiC.tersoff C Si min_style cg #minimize 1.0e-9 1.0e-9 10000 30000 velocity all create ${T} 123456 dist gaussian # 根據(jù)溫度初始化速度 # 在NVT系綜下平衡體系 fix NVT all nvt temp ${T} ${T} 0.1 # Nose-Hoover 熱浴 timestep ${DT} # 設(shè)置積分步長(zhǎng) thermo_style custom step temp press # 觀察溫度和壓強(qiáng)的變化 thermo 1000 # 100 步輸出一次 # 根據(jù)以上設(shè)置跑若干步，目的是使系統(tǒng)達(dá)到熱力學(xué)平衡態(tài) run 100000 # 跑 10000 步 # 在產(chǎn)出階段將用NVE系綜 unfix NVT # 撤掉之前的 NVT 系綜 fix NVE all nve # 換上 NVE 系綜 # 在產(chǎn)出階段計(jì)算熱流 compute KE all ke/atom # 計(jì)算單原子動(dòng)能 compute PE all pe/atom # 計(jì)算單原子勢(shì)能 compute V all stress/atom NULL virial # 計(jì)算單原子位力 compute J all heat/flux KE PE V # 計(jì)算單原子熱流 # 在產(chǎn)出階段計(jì)算熱流自關(guān)聯(lián)函數(shù) <J(0)J(t)> variable Ns equal 100 # 取樣間隔 variable Nc equal 100 # 關(guān)聯(lián)數(shù)據(jù)量（最大關(guān)聯(lián)時(shí)間為 Ns*Nc*DT） variable Np equal ${Ns}*${Nc}*10 # 產(chǎn)出步數(shù)（這是一個(gè)好的選擇） fix HAC all ave/correlate ${Ns} ${Nc} ${Np} c_J[1] c_J[2] c_J[3] & type auto file hac.txt #！�。。。。。。。。。。。。�！ variable M equal 14 # 熱流自關(guān)聯(lián)函數(shù)的個(gè)數(shù) variable N equal ${M}*${Np} # 產(chǎn)出步數(shù) # 跑產(chǎn)出階段 |
回復(fù)此樓» 猜你喜歡
材料371求調(diào)劑
已經(jīng)有10人回復(fù)
-
311求調(diào)劑
已經(jīng)有7人回復(fù)
-
341求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿中科大材料與化工。353分能調(diào)劑到哪？
已經(jīng)有8人回復(fù)
-
337一志愿華南理工0805材料求調(diào)劑
已經(jīng)有7人回復(fù)
-
26讀博
已經(jīng)有4人回復(fù)
-
復(fù)試調(diào)劑
已經(jīng)有3人回復(fù)
-
420求調(diào)劑
已經(jīng)有6人回復(fù)
-
308 085701 四六級(jí)已過求調(diào)劑
已經(jīng)有5人回復(fù)
-
研究生招生
已經(jīng)有7人回復(fù)
/name-?=p/
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 186.9
- 紅花: 1
- 帖子: 49
- 在線: 9.4小時(shí)
- 蟲號(hào): 7570339
- 注冊(cè): 2017-12-05
- 專業(yè): 物理學(xué)I
2樓2020-05-12 18:06:31
lqpzh 
新蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 111.6
- 散金: 5
- 帖子: 90
- 在線: 5.8小時(shí)
- 蟲號(hào): 9150633
- 注冊(cè): 2018-06-21
- 專業(yè): 金屬材料的合金相、相變及
3樓2020-05-13 07:15:26
4樓2020-05-18 08:52:22
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 /name-?=p/ 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +10 | dream…… 2026-03-10 | 12/600 |
|
|---|---|---|---|---|
|
[基金申請(qǐng)]
NSFC申報(bào)書里申請(qǐng)人簡(jiǎn)歷中代表性論著還需要在申報(bào)書最后的附件里面再上傳一遍嗎
 20+4
|
NSFC2026我來了 2026-03-10 | 11/550 |
|
|
|
[考研] 282分材料專業(yè)求調(diào)劑院校 +17 | 楓橋ZL 2026-03-09 | 24/1200 |
|
|
|
[考研] 277求調(diào)劑 +3 | anchor17 2026-03-12 | 3/150 |
|
|
|
[考研] 293求調(diào)劑，一志愿陜師大生物學(xué) +3 | ??????.?.??? 2026-03-09 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑 +11 | 與冬清寧 2026-03-07 | 14/700 |
|
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
|
[碩博家園] 木蟲好像不熱鬧了，是不是？ +4 | 偏振片 2026-03-10 | 4/200 |
|
|
|
[考研] 化工284求調(diào)劑 +12 | 小黑想變白 2026-03-07 | 12/600 |
|
|
|
[碩博家園] 2026級(jí)碩士研究生招生/調(diào)劑 +3 | 知足常樂的樂 2026-03-06 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
|
[考研] 0856材料與化工290求調(diào)劑 +7 | Nebulala 2026-03-08 | 8/400 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-08 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級(jí)有競(jìng)賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
|
[考研] 化學(xué)290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
|
[考研] 0307化學(xué)求調(diào)劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工，前三科206，總分283，求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|