| 5 | 1/1 | 返回列表 |
| 查看: 1483 | 回復(fù): 3 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||
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)生熱流,再測(cè)量溫度梯度,也可以用熱浴,如郎之萬熱浴或者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)生熱流,再測(cè)量溫度梯度,也可以用熱浴,如郎之萬熱浴或者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 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工調(diào)劑,339 +5 | 10213207 2026-03-31 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 已決定調(diào)劑院校 +8 | JKSOIID 2026-03-26 | 8/400 |
|
|
[考研] 一志愿華東理工大學(xué),080500學(xué)碩,317分,求調(diào)劑 +9 | s1145 2026-03-31 | 9/450 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +11 | Nadiums 2026-03-25 | 12/600 |
|
|
[考研] 286分調(diào)劑 +11 | Faune 2026-03-30 | 13/650 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +20 | Mr. Z 2026-03-25 | 21/1050 |
|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
[考研] 求調(diào)劑 +8 | 11ggg 2026-03-30 | 8/400 |
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +10 | 一定上岸aaa 2026-03-27 | 14/700 |
|
|
[考研] 085600,專業(yè)課化工原理,320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|