| 4 | 1/1 | 返回列表 |
| 查看: 1473 | 回復: 3 | ||||
wang06x木蟲 (小有名氣)
|
[求助]
熱傳導邊界問題
|
|
最近看論壇上一些熱傳導例子,幫助很大。首先感謝前輩們的無私奉獻精神。 同時我有些問題很困惑,希望大家給予解答。十分感激! 1.我發(fā)現(xiàn)例子中3個方向都是周期性邊界條件,那么請問G-K公式、MP法 以及fix heat法這三種方法可以使用非周期邊界條件么?如 f p p 2.我自己嘗試了把條件改為f p p可是 出現(xiàn)了丟失原子(lost atoms)的情況,請問這是為什么?看相關帖子說可能是有覆蓋的原子,應該用minimize 命令,可是我用了之后,雖然不像原來的那樣一開始就出現(xiàn)lost atoms ,但是跑了幾步之后仍然出現(xiàn)了。 3.還有的人說可以用fix wall/reflect 命令,我分別嘗試了命令:fix 1 all wall/refiect xlo -150 xhi 150 units box(盒子x方向尺寸是-100-100.單位是metal) 或者fix 2 all wall/reflect xlo EDGE xhi EDGE ?墒沁是會出現(xiàn)lost atoms.是我的命令用錯了么? 4.論壇上例子中NEMD法貌似都是先給定熱流密度的,那么請問有沒有使用分子動力學模擬給定溫度梯度的熱傳導方面相關理論文獻? 問題有點多。水平有限,真的想不出來。希望大家給予解答,十分感謝! |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +218 |
|
第一,用自由邊界條件是可能丟失原子的,這就類似于初中學過的升華現(xiàn)象。所以,用nemd時最好用固定邊界條件,即固定兩端一小部分原子。用nemd方法時也可以用周期邊界條件,此時熱流會從熱源左右兩端流出,從熱匯左右兩端流入。即使用了這兩種邊界條件,系統(tǒng)在高溫下還是有可能不穩(wěn)定,因為這是非平衡模擬,就是不那么穩(wěn)定。第二,可以先用速度標度法或動量交換法產(chǎn)生熱流,再測量溫度梯度,也可以用熱浴,如郎之萬熱浴或者nose-hoover chain熱浴控制熱源和熱匯的溫度。它們都大同小異。第三,如果用emd方法,即格林-久保公式,在輸運方向必須用周期邊界條件,在其它方向的邊界條件視模擬體系而定。以上所說都是熱導率模擬的最基本知識,文獻有幾百篇,其中重要的有幾十篇,樓主可自行查閱,不可看看lammps的例子就拉到。我看過lammps作者關于熱導率模擬的一個ppt,感覺他并不是很懂。 [ 發(fā)自手機版 http://www.gaoyang168.com/3g ] |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +218 |
|
第一,用自由邊界條件是可能丟失原子的,這就類似于初中學過的升華現(xiàn)象。所以,用nemd時最好用固定邊界條件,即固定兩端一小部分原子。用nemd方法時也可以用周期邊界條件,此時熱流會從熱源左右兩端流出,從熱匯左右兩端流入。即使用了這兩種邊界條件,系統(tǒng)在高溫下還是有可能不穩(wěn)定,因為這是非平衡模擬,就是不那么穩(wěn)定。第二,可以先用速度標度法或動量交換法產(chǎn)生熱流,再測量溫度梯度,也可以用熱浴,如郎之萬熱浴或者nose-hoover chain熱浴控制熱源和熱匯的溫度。它們都大同小異。第三,如果用emd方法,即格林-久保公式,在輸運方向必須用周期邊界條件,在其它方向的邊界條件視模擬體系而定。以上所說都是熱導率模擬的最基本知識,文獻有幾百篇,其中重要的有幾十篇,樓主可自行查閱,不可看看lammps的例子就拉到。我看過lammps作者關于熱導率模擬的一個ppt,感覺他并不是很懂。 [ 發(fā)自手機版 http://www.gaoyang168.com/3g ] |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 北科281學碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 求調(diào)劑 +5 | Mqqqqqq 2026-03-19 | 5/250 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 材料學求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 08工學調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-17 | 8/400 |
|
|
[考研] 0856調(diào)劑,是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +4 | 申子申申 2026-03-19 | 8/400 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|