久久超碰99,91一区二区三区四区,精品捆绑调教一区二区三区,日韩国产欧美三级

預存
Document
當前位置：文庫百科 ? 文章詳情
分子動力學(LAMMPS)——實例講解：NEMD方法計算晶格熱導率
來源：科學10分鐘 時間：2021-05-18 17:58:03 瀏覽：15391次

熱導率又稱導熱系數，反映物質的熱傳導能力，其定義為單位溫度梯度（在1m長度內溫度降低1K）在單位時間內經單位導熱面所傳遞的熱量。材料的熱導率是由電子熱導率和晶格熱導率共同貢獻的。在溫度較高的情況下，材料的熱導率主要由電子熱導率主導；在溫度較低的情況，主要由晶格熱導率主導。

實驗中測材料熱導率是在材料兩端施加熱流，此時可以測得沿熱流方向的橫截面上的熱流通量J以及溫度梯度?T ，而J和?T 具有公式（1）的線性關系，比例系數K 就是材料的熱導率。

J = -k ?T （1）

本文我們將使用由美國能源部桑迪亞國家實驗室設計開發的分子動學軟件包LAMMPS[1]復現文獻[2]中用分子動力學計算金屬鎢晶格熱導率的方法。文獻[2]采用的是非平衡分子動力學(NEMD) [3-6]方法計算晶格熱導率。計算過程中的模擬盒子示意圖如圖1所示，當模擬盒子在某一溫度下平衡后，將模擬盒子中間作為熱端，施加熱流?ε；模擬盒子兩端作為冷端對熱流進行抽取，在抽取的過程當中體系中原子的動量是從熱端到冷端不斷交換的，此時體系處于一個非平衡的狀態。因此該過程也叫非平衡分子動力學(NEMD)過程。

在施加熱流后，將模擬盒子沿熱流方向分成許多層，每個特定的timesteps (如10000步)將每一層的溫度進行統計，當傳熱過程穩定后 (這里是指統計的溫度保持基本穩定，一般在加熱流數萬步到十萬步后)，對特定timesteps間隔統計的分層溫度做平均并進行線性擬合得到如圖2所示的溫度梯度曲線（黑色）。得到溫度梯度后根據公式（2）得到晶格熱導率k。公式（2）中A為模擬盒子的橫截面積，M為熱交換的頻率，t為單位時間。最后，結合公式（1）和公式（2）得出金屬鎢的晶格熱導率。筆者取的?ε/t =1eV/ps，A=16?2（筆者這里設置的有些不合理但對計算結果影響不大），M =10。筆者計算得出W在300K時延<100>方向的熱導率是16.7 W/KM，和文獻的計算結果17.28034 W/KM基本一致[2]。


圖1. NEMD計算熱導率的三維模型示意圖

圖2. 模擬盒子分層溫度統計曲線(紅色),溫度梯度擬合曲線(黑色)

 在計算過程當中，LAMMPS的in文件中的傳熱過程設置腳本段如下：

unfix            1

fix              2 all nve

fix              hot all heat 10 ${power} region hot

fix              cold all heat 10 -${power} region cold

其中，unfix 1命令取消的是前面作特定溫度下平衡的系綜。其余三段為熱流為傳熱過程的設置。

計算結束后，最關鍵的是對模擬盒子每一層原子溫度的統計，為此我們寫了一個python小腳本，可以直接運行python3 thermal_cond.py profile.txt 200 25 45 > result.txt輕松處理，python的用法，及函數的定義，腳本thermal_cond.py中都有注釋。處理后得到文本文件result.txt，該文件中的信息是模擬盒子每一層溫度與層中心到模擬盒子左端距離的關系。直接將result.txt拖到作圖軟件中，或用gnuplot等程序作圖即可獲得火山峰狀的溫度統計曲線。最后截取其中的直線段擬合即可得溫度梯度?T，本文的?T=-0.6K/?。

最后，介紹幾個可以使結果更優的方法：

（1）將模擬盒子三個方向邊長，建成W單包三個方向邊長的倍數，這樣模擬盒子具有完美晶體的周期性。

（2）將采集溫度的模擬時間盡量增大，這樣能夠采集更多的樣本，統計的溫度就會更加準確。

該算例以及溫度梯度處理腳本可以在網盤下載：

鏈接：https://pan.baidu.com/s/1zkpULJ2T12NFkgnQOt5KOA

提取碼：CSGO

參考文獻：

[1] Plimpton S 1995 J. Comput. Phys. 117 1

[2] Baoqin Fu,Wensheng Lai,Yue Yuan, et al. 2012 J. Nucl. Mater. 427 1-3

[3] P. K. Schelling, S.R. Phillpot, P. Keblinski. 2002 Phys. Rev. B 65 144306.

[4] A. Maiti, G.D. Mahan, S.T. Pantelides. 1997 Solid State Commun. 102 517.

[5] A. Maeda. 1995 Phys. Rev. E 52 34.

[6] A.R. Abramson, C.L. Tien, A. Majumdar. 2002 J. Heat Transfer 124 963.
文章精選 


你的點贊和在看，我都認真當成喜歡
評論 / 文明上網理性發言
12條評論
全部評論 / 我的評論
最熱 /  最新
全部 3小時前 四川
文字是人類用符號記錄表達信息以傳之久遠的方式和工具?，F代文字大多是記錄語言的工具。人類往往先有口頭的語言后產生書面文字，很多小語種，有語言但沒有文字。文字的不同體現了國家和民族的書面表達的方式和思維不同。文字使人類進入有歷史記錄的文明社會。
點贊12
回復
全部
查看更多評論
相關文章

基礎理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結構）

2020-05-03

手把手教你用ChemDraw 畫化學結構式：基礎篇

2021-06-19

晶體結構可視化軟件 VESTA使用教程（下篇）

2021-01-22

【科研干貨】電化學表征：循環伏安法詳解（上）

2019-10-25

【科研干貨】電化學表征：循環伏安法詳解（下）

2019-10-25

XRD的基本原理與應用

2020-11-03

熱門文章/popular

基礎理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結構）

手把手教你用ChemDraw 畫化學結構式：基礎篇

晶體結構可視化軟件 VESTA使用教程（下篇）

電化學實驗基礎之電化學工作站篇 （二）三電極和兩電極體系的搭建 和測試

【科研干貨】電化學表征：循環伏安法詳解（上）

【科研干貨】電化學表征：循環伏安法詳解（下）

微信掃碼分享文章
久久超碰99,91一区二区三区四区,精品捆绑调教一区二区三区,日韩国产欧美三级
日韩成人综合| 欧美一级二级三级视频| 国产亚洲观看| 亚洲综合二区| 亚洲视频播放| 国产精品日韩欧美一区| 午夜在线播放视频欧美| 亚洲视频播放| 亚洲专区视频| 日韩精品视频中文字幕| 青青草精品视频| 国产欧美日韩精品一区二区免费| 欧美久久久网站| 国产精品亚洲综合久久| 麻豆精品在线视频| 四虎成人av| 亚洲成人精品| 中文一区二区| 中文字幕一区二区三区四区久久| 亚洲一二三区视频| 国产另类在线| 福利在线免费视频| 欧美+日本+国产+在线a∨观看| 精品一区毛片| 免费人成精品欧美精品| 日韩高清二区| 久久精品九色| 99精品国产一区二区三区| 9国产精品视频| 日韩精品午夜视频| 成人在线免费观看91| 久久国产直播| 模特精品在线| 91成人在线网站| 国产一区二区三区久久| 亚洲91视频| 亚洲精品黄色| 里番精品3d一二三区| 成人精品天堂一区二区三区| 亚洲欧美日韩国产一区二区| 日韩精品国产欧美| 精品五月天堂| 女同性一区二区三区人了人一 | 国产精品久久久免费| 国产999精品在线观看| 国产一区二区三区自拍| 亚洲自啪免费| 国产精品美女午夜爽爽| 电影亚洲精品噜噜在线观看| 久久先锋影音| 免费视频一区二区三区在线观看 | 国产精品99久久免费观看| 欧美精品资源| 日日夜夜免费精品视频| 精品国产乱码久久久久久1区2匹| 亚洲小说欧美另类婷婷| 日韩高清在线一区| 给我免费播放日韩视频| 国产视频一区在线观看一区免费| 国产精品亚洲产品| 亚洲国产成人精品女人| 国产伦理久久久久久妇女| 美女少妇全过程你懂的久久| 日本电影久久久| 欧美精品日日操| 日韩欧美美女在线观看| 日韩欧美视频专区| 日本aⅴ精品一区二区三区| 91亚洲一区| 日韩有码av| 国产理论在线| 日韩高清中文字幕一区| 啪啪国产精品| 欧美精品中文字幕亚洲专区| 久久激情婷婷| 国产精品久久久久毛片大屁完整版| 青青久久av| 欧美日韩亚洲一区二区三区在线| 亚洲午夜视频| 国产一区二区三区久久 | 精品视频久久| 亚洲精品欧洲| 99久久婷婷这里只有精品| 国产极品一区| 日韩视频不卡| 婷婷激情一区| 欧美中文高清| 国产午夜精品一区二区三区欧美 | 97精品久久| 国产一区二区高清| 亚洲精品国产嫩草在线观看| 国产精品一区二区av日韩在线| 日韩在线卡一卡二| 天堂资源在线亚洲| 国产精品久久久久久久久妇女| 日本大胆欧美人术艺术动态| 欧美精选视频一区二区| 精品欧美视频| 久久国产精品免费精品3p | 国产日产精品_国产精品毛片 | 黄色av一区| 日韩欧美国产精品综合嫩v| 久久99性xxx老妇胖精品| 蜜桃视频在线观看一区| 欧美国产91| 99成人在线视频| 日韩精品中文字幕第1页| 国产精品国码视频| 美女91精品| 亚洲电影在线一区二区三区| 久久久久久黄| 在线天堂中文资源最新版| 国产精品一区毛片| 日韩精品视频在线看| 爽好久久久欧美精品| 久久国产电影| 91亚洲无吗| 亚洲精品女人| 欧美一区成人| 日韩国产在线观看一区| 亚洲天堂免费| 欧美专区在线| 久久99伊人| 免费精品一区| 国产日本精品| 国产成人a视频高清在线观看| 国产精品欧美在线观看| 日韩精品欧美| 久久理论电影| 日韩欧美中文| 国内精品麻豆美女在线播放视频| 69精品国产久热在线观看| 欧美日本精品| 国产精品magnet| 九九99久久精品在免费线bt| 久久久久黄色| 91日韩免费| 亚洲天堂1区| 欧美粗暴jizz性欧美20| 不卡在线一区| 久久成人一区| 亚洲精品一级| 日韩激情综合| 国产精品一区毛片| 成人台湾亚洲精品一区二区| 红杏一区二区三区| 亚洲精品成人图区| 99久久www免费| 国产韩日影视精品| 午夜在线精品偷拍| 日本成人在线视频网站| 国产精品久久久网站| 九九久久国产| 亚洲h色精品| 免费在线欧美视频| 日本欧美一区| 国产精品丝袜在线播放| 国产成人免费视频网站视频社区| 亚洲精品.com| 99在线|亚洲一区二区| 午夜久久av| 国产精选一区| 亚洲黄色免费看| 99国产精品| 视频一区日韩精品| 精品国产一区二区三区2021| 综合日韩av| 精品中文一区| 日韩精彩视频在线观看| 岛国精品一区| 欧美精品一卡| 久久国产乱子精品免费女| 亚洲美女久久精品| 天堂va蜜桃一区二区三区| 国产激情久久| 亚洲91视频| 蜜臀av一区二区三区| 国产精品1区在线| 欧美日韩色图| 日本成人在线网站| 婷婷综合六月| 日韩精品导航| 亚洲va中文在线播放免费| 亚洲日产av中文字幕| 成人污污视频| 日韩影院在线观看| 国产一区二区精品福利地址| 亚洲精华国产欧美| 六月丁香综合在线视频| 红桃视频国产精品| 久久伊人久久| 久久国产中文字幕| 欧美亚洲三区| 婷婷亚洲五月| 国产精品久久久久av蜜臀| 在线亚洲精品| 日本久久综合| 亚洲精品少妇| 日韩高清中文字幕一区二区| 欧美一级网址|
+

你好，很高興為您服務！

發送