 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1714 | 回復: 3 | |||
/name-?=p/銅蟲 (初入文壇)
|
[交流]
lammps在npt平衡后換成nve溫度一直升高是怎么回事？ 已有2人參與
|
|
各位大神請教一個問題 用lammps模擬計算3C-SiC熱導率,時間步長1fs 在1000K下，先在NVT下模擬10萬步，使系統(tǒng)達到平衡然后換上NVE系綜后在運行125萬步。但是當換上NVE系綜不久后，溫度就開始升高，到結束時溫度都達到幾千甚至上萬K，這么高溫度碳化硅固體早就沒了。125萬步不是設定步數，125萬步時原子丟失了，在低溫時 溫度也會升高 想問這是哪里出問題了 這是in文件 晶格常數設置的讀坐標文件 盒子邊長大小是十個晶格常數也就是一共8000個原子 variable T equal 1000 # 溫度為 400 K #variable A equal 4.348 # 晶格常數大概為 6.35A variable DT equal 0.001 # 積分步長為 0.3 fs # 初始化坐標、勢函 數、速度等數據 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 # 根據溫度初始化速度 # 在NVT系綜下平衡體系 fix NVT all nvt temp ${T} ${T} 0.1 # Nose-Hoover 熱浴 timestep ${DT} # 設置積分步長 thermo_style custom step temp press # 觀察溫度和壓強的變化 thermo 1000 # 100 步輸出一次 # 根據以上設置跑若干步，目的是使系統(tǒng)達到熱力學平衡態(tài) run 100000 # 跑 10000 步 # 在產出階段將用NVE系綜 unfix NVT # 撤掉之前的 NVT 系綜 fix NVE all nve # 換上 NVE 系綜 # 在產出階段計算熱流 compute KE all ke/atom # 計算單原子動能 compute PE all pe/atom # 計算單原子勢能 compute V all stress/atom NULL virial # 計算單原子位力 compute J all heat/flux KE PE V # 計算單原子熱流 # 在產出階段計算熱流自關聯函數 <J(0)J(t)> variable Ns equal 100 # 取樣間隔 variable Nc equal 100 # 關聯數據量（最大關聯時間為 Ns*Nc*DT） variable Np equal ${Ns}*${Nc}*10 # 產出步數（這是一個好的選擇） 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 # 熱流自關聯函數的個數 variable N equal ${M}*${Np} # 產出步數 # 跑產出階段 |
回復此樓» 猜你喜歡
085600材料與化工 292分 求調劑
已經有3人回復
-
接受26屆調劑生
已經有8人回復
-
273求調劑
已經有7人回復
-
A區(qū)一本交叉課題組，低分調劑，招收機械電子信息通信等交叉方向
已經有45人回復
-
334求調劑
已經有3人回復
-
求調劑推薦
已經有3人回復
-
求調劑
已經有4人回復
-
一志愿211，0860總分286食品工程專業(yè)求調劑
已經有6人回復
-
招調劑
已經有14人回復
-
誠求調劑，323分 有專利/科研/實習經歷
已經有6人回復
/name-?=p/
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 186.9
- 紅花: 1
- 帖子: 49
- 在線: 9.4小時
- 蟲號: 7570339
- 注冊: 2017-12-05
- 專業(yè): 物理學I
2樓2020-05-12 18:06:31
lqpzh 
新蟲 (小有名氣)
- 應助: 4 (幼兒園)
- 金幣: 111.6
- 散金: 5
- 帖子: 90
- 在線: 5.8小時
- 蟲號: 9150633
- 注冊: 2018-06-21
- 專業(yè): 金屬材料的合金相、相變及
3樓2020-05-13 07:15:26
4樓2020-05-18 08:52:22
相關版塊跳轉
我要訂閱樓主 /name-?=p/ 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調劑 +12 | 帆船哥 2026-03-04 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工 298 調劑 +10 | 小西笑嘻嘻 2026-03-03 | 10/500 |
|
|
|
[考研] 材料277分求調劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
|
[考研] 材料專碩290求調劑 +8 | 杰尼龜aaa 2026-03-04 | 8/400 |
|
|
|
[考研] 0703 東華大學 理學化學均過 a 線，總分 281 求調劑 +6 | 香香 Lu 2026-03-03 | 11/550 |
|
|
|
[考研] 274求調劑 +6 | cgyzqwn 2026-03-01 | 13/650 |
|
|
|
[考研] 復試調劑 +5 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 304分材料專碩求調劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
|
[考研] 化學 0703求調劑 總分293 一志愿211 +9 | 土土小蟲 2026-03-03 | 11/550 |
|
|
|
[考研] 材料與化工，291分，求調劑 +9 | 咕嚕咕嚕123123 2026-03-04 | 11/550 |
|
|
|
[基金申請] 沒有青基直接申請面上，感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
|
[考研] 322分 085600求調劑，有互聯網+國金及主持省級大創(chuàng)經歷 +6 | 熊境喆 2026-03-04 | 6/300 |
|
|
|
[考研] 295求調劑 +6 | 等春來， 2026-03-04 | 6/300 |
|
|
|
[考研] 能動297求調劑，本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
|
[考研] 281求調劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 環(huán)境調劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 化工專碩調劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
|
[考研] 一志愿中科大能動297求調劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 0857調劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
|
[考研] 272求調劑 +6 | 田智友 2026-02-28 | 6/300 |
|