| 4 | 1/1 | 返回列表 |
| 查看: 1474 | 回復(fù): 3 | ||||
wang06x木蟲 (小有名氣)
|
[求助]
熱傳導(dǎo)邊界問題
|
|
最近看論壇上一些熱傳導(dǎo)例子,幫助很大。首先感謝前輩們的無私奉獻(xiàn)精神。 同時(shí)我有些問題很困惑,希望大家給予解答。十分感激! 1.我發(fā)現(xiàn)例子中3個(gè)方向都是周期性邊界條件,那么請(qǐng)問G-K公式、MP法 以及fix heat法這三種方法可以使用非周期邊界條件么?如 f p p 2.我自己嘗試了把條件改為f p p可是 出現(xiàn)了丟失原子(lost atoms)的情況,請(qǐng)問這是為什么?看相關(guān)帖子說可能是有覆蓋的原子,應(yīng)該用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 ?墒沁是會(huì)出現(xiàn)lost atoms.是我的命令用錯(cuò)了么? 4.論壇上例子中NEMD法貌似都是先給定熱流密度的,那么請(qǐng)問有沒有使用分子動(dòng)力學(xué)模擬給定溫度梯度的熱傳導(dǎo)方面相關(guān)理論文獻(xiàn)啊? 問題有點(diǎn)多。水平有限,真的想不出來。希望大家給予解答,十分感謝! |

新蟲 (初入文壇)
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +218 |
|
第一,用自由邊界條件是可能丟失原子的,這就類似于初中學(xué)過的升華現(xiàn)象。所以,用nemd時(shí)最好用固定邊界條件,即固定兩端一小部分原子。用nemd方法時(shí)也可以用周期邊界條件,此時(shí)熱流會(huì)從熱源左右兩端流出,從熱匯左右兩端流入。即使用了這兩種邊界條件,系統(tǒng)在高溫下還是有可能不穩(wěn)定,因?yàn)檫@是非平衡模擬,就是不那么穩(wěn)定。第二,可以先用速度標(biāo)度法或動(dòng)量交換法產(chǎn)生熱流,再測量溫度梯度,也可以用熱浴,如郎之萬熱浴或者nose-hoover chain熱浴控制熱源和熱匯的溫度。它們都大同小異。第三,如果用emd方法,即格林-久保公式,在輸運(yùn)方向必須用周期邊界條件,在其它方向的邊界條件視模擬體系而定。以上所說都是熱導(dǎo)率模擬的最基本知識(shí),文獻(xiàn)有幾百篇,其中重要的有幾十篇,樓主可自行查閱,不可看看lammps的例子就拉到。我看過lammps作者關(guān)于熱導(dǎo)率模擬的一個(gè)ppt,感覺他并不是很懂。 [ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +218 |
|
第一,用自由邊界條件是可能丟失原子的,這就類似于初中學(xué)過的升華現(xiàn)象。所以,用nemd時(shí)最好用固定邊界條件,即固定兩端一小部分原子。用nemd方法時(shí)也可以用周期邊界條件,此時(shí)熱流會(huì)從熱源左右兩端流出,從熱匯左右兩端流入。即使用了這兩種邊界條件,系統(tǒng)在高溫下還是有可能不穩(wěn)定,因?yàn)檫@是非平衡模擬,就是不那么穩(wěn)定。第二,可以先用速度標(biāo)度法或動(dòng)量交換法產(chǎn)生熱流,再測量溫度梯度,也可以用熱浴,如郎之萬熱浴或者nose-hoover chain熱浴控制熱源和熱匯的溫度。它們都大同小異。第三,如果用emd方法,即格林-久保公式,在輸運(yùn)方向必須用周期邊界條件,在其它方向的邊界條件視模擬體系而定。以上所說都是熱導(dǎo)率模擬的最基本知識(shí),文獻(xiàn)有幾百篇,其中重要的有幾十篇,樓主可自行查閱,不可看看lammps的例子就拉到。我看過lammps作者關(guān)于熱導(dǎo)率模擬的一個(gè)ppt,感覺他并不是很懂。 [ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 298求調(diào)劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|