全頻蒙地卡羅材料熱傳模擬工具之研發

Development of a full-spectrum Monte-Carlo simulation tool for phonon flow

簡培盛(Pei-Sheng Chien)

指導教授 ： 黃美嬌

國立臺灣大學/工學院/機械工程學研究所/碩士(2015年)

https://doi.org/10.6342/NTU.2015.02332

摘要

本論文主要研究目的為建立一全頻蒙地卡羅法(MC)數值工具來分析半導體材料內的熱傳導行為，藉由此數值工具可追蹤並記錄聲子的運動及分布情形，進而求得其不同頻率及不同極化之聲子對材料熱傳導係數的貢獻。
蒙地卡羅法的理論基礎為聲子波茲曼傳輸方程式(PBTE)，在求解聲子波茲曼傳輸方程式的過程中，聲子性質部分是由實驗量測之材料色散關係以及溫度來決定，每顆聲子皆擁有各自獨立的速度、頻率及極化，並因而有各自各種散射機制的散射頻率，其中散射機制使用的是Holland所提出之經驗方程式；同時亦建立及設計了適當的數值邊界條件，並利用存在於系統中的統計對稱性條件來減少計算量。
本論文選擇鍺(Ge)塊材內的聲子傳輸現象為測試例。首先比較有無灰介質假設下之聲子性質。研究結果說明在灰介質假設下之聲子群速比多數頻率之聲子群速低，但卻高於聲子量最多的高頻橫波聲子群速；灰介質假設下之聲子平均散射頻率也普遍較低，但高於熱傳貢獻較大的低頻縱波聲子的散射頻率；綜合以上兩者，灰介質假設下之聲子平均自由路徑較長，但低於熱傳貢獻較大的低頻縱波聲子的平均自由路徑。最後在熱傳導係數的貢獻上，雖然縱波低頻聲子對比熱的貢獻不及高頻聲子，但其速度較快且散射頻率較低，故在材料內其熱傳的貢獻最大。
在模擬工具的驗證方面，研究結果顯示模擬過程中有三種情況需重設聲子的性質:當聲子發生本質散射時、當聲子撞擊材料邊界發生邊界散射時、以及為了滿足能量守恆而需增加聲子時。重設聲子性質的方法為影響模擬準確性的關鍵因素，決定聲子的性質時，應分別依聲子發生事件之機率為比例進行聲子性質的重新設定，因此前述三種情況應分別以本質散射頻率為權重、聲子速率為權重、及無權重下之平衡聲子分佈函數為比例，配合電腦亂數隨機挑選之。最後研究設計一熱流控制邊界條件來產生非平衡系統，進而量測得鍺塊材之熱傳導係數，結果顯示理論及模擬之值相吻合，本研究成功開發出一全頻蒙地卡羅法數值模擬工具。

關鍵字

全頻 ； 頻率相依 ； 蒙地卡羅法 ； 權重頻率分布 ； 熱傳導係數

並列摘要

This study aims at developing a full-spectrum Monte-Carlo (MC) simulation tool for phonon flow in semiconductors. The phonon Boltzmann transport equation (PBTE) was solved based on the bulk phonon dispersion relation and Holland’s empirical relations for various scattering mechanisms. Most of all, this research reveals proper ways to reset the phonon properties for different scattering mechanisms, namely, it must be done based on the actual probabilities of scattering events.
A 1D system with either periodic boundary conditions, adiabatic boundary conditions, or constant boundary temperatures was employed to verify the accuracy of the proposed simulation tool. The study demonstrates that the phonon properties must be randomly selected from the un-weighted local equilibrium phonon population distribution, the weighted one by the three-phonon scattering rate, or the weighted one by the phonon group velocity respectively when a local energy imbalance, a three-phonon scattering event, or a boundary scattering event occurs. The accuracy of the proposed simulation tool was finally confirmed via a transient problem and a steady problem.

並列關鍵字

Full-spectrum ； Monte-Carlo simulation ； phonon flow ； thermal conductivity

參考文獻

[20] 蔡東峻, "奈米複合材料聲子傳輸現象蒙地卡羅模擬法之研發" , 碩士, 臺灣大學機械工程學研究所學位論文, 台灣大學, 2009.

[22] 莊璧躍, "不規則奈米顆粒複合物熱傳導性質研究" , 碩士, 臺灣大學工學院機械工程學研究所學位論文, 台灣大學, 2012.

[2] H. J. Goldsmid and R.W. Douglas, "The use of semiconductors in thermoelectric refrigeration" , British Journal of Applied Physics, vol.5, pp. 386-390, 1954.

[3] J. Yang and T. Caillat, "Thermoelectric materials for space and automotive power generation" , MRS Bulletin, vol.31, pp. 224-229, 2006.

[4] R. P. Chasmar and R. Stratton, "The Thermoelectric Figure of Merit and its Relation to Thermoelectric Generators" , Journal of Electronic and Control, vol.7, pp. 52-72, 1959.