 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1718 | 回復: 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ù) # 跑產出階段 |
回復此樓» 猜你喜歡
新疆大學地質與礦業(yè)工程學院招生
已經有32人回復
-
0856材料與化工353分求調劑
已經有4人回復
-
282分材料專業(yè)求調劑院校
已經有13人回復
-
315食品工程求調劑
已經有5人回復
-
環(huán)境調劑
已經有3人回復
-
接受26屆調劑生
已經有20人回復
-
考研材料與化工，求調劑
已經有3人回復
-
歡迎211本科同學，過A區(qū)國家線，A區(qū)非偏遠一本，交叉學科課題組
已經有35人回復
-
304求調劑（085602一志愿985）
已經有4人回復
-
294 英二數(shù)二物化 求調劑
已經有6人回復
lqpzh 
新蟲 (小有名氣)
- 應助: 4 (幼兒園)
- 金幣: 111.6
- 散金: 5
- 帖子: 90
- 在線: 5.8小時
- 蟲號: 9150633
- 注冊: 2018-06-21
- 專業(yè): 金屬材料的合金相、相變及
3樓2020-05-13 07:15:26
/name-?=p/
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 186.9
- 紅花: 1
- 帖子: 49
- 在線: 9.4小時
- 蟲號: 7570339
- 注冊: 2017-12-05
- 專業(yè): 物理學I
2樓2020-05-12 18:06:31
4樓2020-05-18 08:52:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境調劑 +3 | 曉看天暮看云 2026-03-09 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 0817化學工程319求調劑 +5 | lv945 2026-03-08 | 6/300 |
|
|
|
[碩博家園] 2026級碩士研究生招生/調劑 +3 | 知足常樂的樂 2026-03-06 | 5/250 |
|
|
|
[考研] 一志愿天津大學，英一數(shù)二305分求調劑，四六級已過 +5 | 小小番的茄 2026-03-09 | 5/250 |
|
|
[考研] 求調劑，數(shù)一英一274分 +4 | 小菲會努力 2026-03-08 | 4/200 |
|
|
|
[考研] 337求調劑 +3 | 睡醒，。 2026-03-09 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程專碩求調劑 +5 | 幼稚園大班打手 2026-03-08 | 5/250 |
|
|
|
[考研] 調劑 +3 | 13853210211 2026-03-08 | 3/150 |
|
|
|
[考研] 求調劑，一志愿華中科大0702，數(shù)一英一，293 +4 | 小羅露一二 2026-03-07 | 4/200 |
|
|
|
[考研] 0307化學求調劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
|
[考研] 278求調劑 +5 | Gale1314 2026-03-06 | 5/250 |
|
|
|
[考研] 材料專碩290求調劑 +8 | 杰尼龜aaa 2026-03-04 | 8/400 |
|
|
|
[考研] 化學專業(yè)調劑 +5 | 好好好1233 2026-03-04 | 6/300 |
|
|
|
[考研] 0856材料與化工求調劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 材料考研339求調劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
|
[考研] 334求調劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|
|
[考研] 085602 293分求調劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
|
[考研] 求調劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|