 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1717 | 回復: 3 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
/name-?=p/銅蟲 (初入文壇)
|
[交流]
lammps在npt平衡后換成nve溫度一直升高是怎么回事？ 已有2人參與
|
||
|
各位大神請教一個問題 用lammps模擬計算3C-SiC熱導率,時間步長1fs 在1000K下，先在NVT下模擬10萬步，使系統(tǒng)達到平衡然后換上NVE系綜后在運行125萬步。但是當換上NVE系綜不久后，溫度就開始升高，到結束時溫度都達到幾千甚至上萬K，這么高溫度碳化硅固體早就沒了。125萬步不是設定步數(shù)，125萬步時原子丟失了，在低溫時 溫度也會升高 想問這是哪里出問題了 這是in文件 晶格常數(shù)設置的讀坐標文件 盒子邊長大小是十個晶格常數(shù)也就是一共8000個原子 variable T equal 1000 # 溫度為 400 K #variable A equal 4.348 # 晶格常數(shù)大概為 6.35A variable DT equal 0.001 # 積分步長為 0.3 fs # 初始化坐標、勢函 數(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} # 設置積分步長 thermo_style custom step temp press # 觀察溫度和壓強的變化 thermo 1000 # 100 步輸出一次 # 根據(jù)以上設置跑若干步，目的是使系統(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 # 計算單原子熱流 # 在產出階段計算熱流自關聯(lián)函數(shù) <J(0)J(t)> variable Ns equal 100 # 取樣間隔 variable Nc equal 100 # 關聯(lián)數(shù)據(jù)量（最大關聯(lián)時間為 Ns*Nc*DT） variable Np equal ${Ns}*${Nc}*10 # 產出步數(shù)（這是一個好的選擇） 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 # 熱流自關聯(lián)函數(shù)的個數(shù) variable N equal ${M}*${Np} # 產出步數(shù) # 跑產出階段 |
回復此樓» 猜你喜歡
求調劑，數(shù)一英一274分
已經有4人回復
-
材料調劑，307分
已經有4人回復
-
320材料與化工，求調劑
已經有11人回復
-
296求調劑
已經有4人回復
-
308求調劑
已經有6人回復
-
085600材料與化工調劑 280分
已經有15人回復
-
求材料讀博院校
已經有8人回復
-
申博
已經有8人回復
-
復試調劑
已經有3人回復
-
中科大材料299求調劑
已經有16人回復
4樓2020-05-18 08:52:22
/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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 求材料讀博院校 +8 | yanglei131 2026-03-08 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 290 材料與化工求調劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
[考研] 086000生物與醫(yī)藥319分求調劑 +3 | Tolkien 2026-03-07 | 3/150 |
|
|
|
[考研] 求調劑，一志愿江南大學，食品科學與工程，總分，320 +3 | yyyyyukino 2026-03-07 | 3/150 |
|
|
|
[考研] 調劑 +3 | 13853210211 2026-03-08 | 3/150 |
|
|
|
[考研] 288求調劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
|
[考研] 303求調劑 +8 | forgman95 2026-03-05 | 10/500 |
|
|
|
[考研] 材料化工求調劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
|
[考研] 求調劑 一志愿蘇州大學，0856化工323分 | 本科應化 | 有專利/競賽/科研助手經歷 | +5 | 橙子cyx 2026-03-06 | 6/300 |
|
|
|
[考研] 266材料化工求調劑 +5 | 哇塞王帥 2026-03-03 | 5/250 |
|
|
|
[考研] 08工科求調劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
|
[考研] 不限學校專業(yè)的調劑同學看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
|
[考博] 2026申博自薦 六級440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
|
[考研] 274求調劑 +9 | 一個學習者 2026-03-04 | 9/450 |
|
|
|
[考研] 材料調劑 +5 | ws 上岸鴨 2026-03-05 | 5/250 |
|
|
|
[考研] 求材料調劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
|
[考研] 材料328求調劑 +10 | 一個蘿卜02 2026-03-03 | 10/500 |
|
|
|
[考研] 286 +6 | ksncj 2026-03-04 | 6/300 |
|
|
|
[考研] 0857求調劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
|
[考研] 11408，學碩276求調劑 +3 | 崔wj 2026-03-02 | 5/250 |
|