成人3d动漫一区二区三区_五月天婷婷在线观看视频_999精品视频在线_久久免费看毛片_日本女优爱爱视频_裸体大乳女做爰69_成人午夜视频免费在线观看_777久久精品一区二区三区无码_欧美 日韩 激情_蜜臀av性久久久久蜜臀av_国产九九在线视频_日本丰满少妇xxxx

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

熱導率又稱導熱系數，反映物質的熱傳導能力，其定義為單位溫度梯度（在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小時前 四川
文字是人類用符號記錄表達信息以傳之久遠的方式和工具。現代文字大多是記錄語言的工具。人類往往先有口頭的語言后產生書面文字，很多小語種，有語言但沒有文字。文字的不同體現了國家和民族的書面表達的方式和思維不同。文字使人類進入有歷史記錄的文明社會。
點贊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使用教程（下篇）

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

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

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

微信掃碼分享文章
成人3d动漫一区二区三区_五月天婷婷在线观看视频_999精品视频在线_久久免费看毛片_日本女优爱爱视频_裸体大乳女做爰69_成人午夜视频免费在线观看_777久久精品一区二区三区无码_欧美 日韩 激情_蜜臀av性久久久久蜜臀av_国产九九在线视频_日本丰满少妇xxxx
中文字幕日韩亚洲| 久久婷婷国产| 视频精品一区| 国产色播av在线| 亚洲欧美在线综合| 亚洲精品一级二级| 久久中文欧美| 国产精品最新自拍| 日本v片在线高清不卡在线观看| 深夜视频一区二区| 97人人精品| 免费看久久久| 国产精品日本一区二区三区在线 | 在线亚洲欧美| 日韩另类视频| 欧美国产偷国产精品三区| 国产日韩精品视频一区二区三区| 免费中文字幕日韩欧美| 亚洲欧美日韩国产| 丝袜诱惑制服诱惑色一区在线观看| 久久精品国产亚洲夜色av网站| 精品国产鲁一鲁****| 国产欧美另类| 久久99国产精品视频| 精品视频在线观看网站| 高清日韩欧美| 国产一在线精品一区在线观看| 精品国产第一福利网站| 九九99久久精品在免费线bt| 久久香蕉精品香蕉| 亚洲一区资源| 99精品电影| 中文字幕av一区二区三区四区| 在线日韩成人| 国产精品chinese| 精品午夜视频| 天堂资源在线亚洲| 爽好久久久欧美精品| 日本视频一区二区| 美女视频黄免费的久久| 成人羞羞在线观看网站| 欧美专区一区二区三区| 欧美亚洲免费| 精品欧美一区二区三区在线观看| 欧美日韩三区| 久久99久久久精品欧美| 狠狠操综合网| 欧美日本精品| 婷婷丁香综合| 91精品啪在线观看国产爱臀| 美女av在线免费看| 蜜桃久久久久久| 成人午夜在线| 自拍自偷一区二区三区| 精品视频91| 爽好多水快深点欧美视频| 欧美精品91| 影音国产精品| 卡一卡二国产精品| 免费中文字幕日韩欧美| 成人国产精品一区二区免费麻豆| 夜久久久久久| 91亚洲一区| 欧美日韩一视频区二区| 激情综合网站| 97精品一区二区| 日韩av在线免费观看不卡| 亚洲午夜在线| 精品久久国产一区| 一区二区三区国产盗摄| 色爱综合av| 国内自拍视频一区二区三区| 日韩视频1区| 黄色日韩在线| 成人日韩在线| 国产成人77亚洲精品www| 久久精品 人人爱| 六月天综合网| 国产二区精品| 91精品久久久久久久久久不卡| 久久精品三级| 精品日韩在线| 国产精品一区二区精品视频观看| 六月婷婷一区| 视频一区在线视频| 国产视频一区三区| 久久精品99久久无色码中文字幕| 国产一区二区视频在线看| 在线精品国产亚洲| 激情婷婷久久| 天堂网在线观看国产精品| 在线精品亚洲欧美日韩国产| 精品美女在线视频| 国产情侣一区| 国产精品激情电影| 美女视频一区在线观看| 久久的色偷偷| 精品精品99| 黄色在线观看www| 欧美亚洲国产激情| 亚洲成人日韩| 丝袜美腿一区二区三区| 日本少妇精品亚洲第一区| 青青草国产精品亚洲专区无| 日韩va亚洲va欧美va久久| 欧美一区二区三区久久精品| 免费一级欧美片在线观看网站| 国产美女久久| 国产福利电影在线播放| 久久激情中文| 亚洲欧美日韩综合国产aⅴ| 亚洲一区二区三区高清不卡| 最新国产精品| 美女性感视频久久| 亚洲精品福利电影| 欧美va亚洲va日韩∨a综合色| 国产一区日韩一区| 久久亚洲影院| 国产精品18| 欧美日韩一区二区三区视频播放| 丝袜美腿成人在线| 国产日韩视频| 日韩高清欧美| 五月激激激综合网色播| 水蜜桃精品av一区二区| 免费成人网www| 日韩精品社区| 亚洲国产福利| 日韩中文av| 久久国产日韩| 国产精品嫩模av在线| 欧美香蕉视频| 日本视频一区二区| 最新中文字幕在线播放| 首页国产欧美久久| 日韩综合在线| 色8久久久久| 欧美片第1页| 国产日韩欧美| 日韩午夜电影| 久久精品国内一区二区三区| 日韩一级欧洲| 色婷婷亚洲mv天堂mv在影片| 日本 国产 欧美色综合| 美女av在线免费看| 91福利精品在线观看| 另类小说一区二区三区| 日韩影院在线观看| 久久亚洲精品中文字幕蜜潮电影| 欧美一级网站| 亚洲伦乱视频| 欧美成人精品一级| 日本免费新一区视频| 亚洲欧美综合| 午夜精品成人av| 精品国产亚洲日本| 久久精品99国产精品| 久久午夜精品| 欧美午夜不卡| 成人日韩在线| 国产成人在线中文字幕| 国产欧美一区二区精品久久久| 红桃视频国产精品| 99久久久久| 欧美少妇精品| 一区二区三区四区日本视频| 久久在线91| 久久一区国产| 欧美国产另类| 国产乱码精品| 日本色综合中文字幕| 伊人久久亚洲| 免费欧美在线视频| 午夜宅男久久久| 中文无码日韩欧| 日韩精品一区第一页| 久久亚洲美女| 婷婷五月色综合香五月| 视频一区视频二区中文| 蜜臀91精品一区二区三区| 99香蕉国产精品偷在线观看| 欧美午夜不卡| 久久国产精品99国产| 日韩中文一区二区| 久久精品xxxxx| 欧美国产中文高清| 精品国产不卡一区二区| 欧美韩日一区| 日本欧美不卡| 午夜电影亚洲| 视频一区视频二区中文字幕| 蜜臀久久久99精品久久久久久| 国产精品毛片| 一区二区国产在线观看| 日本午夜精品视频在线观看| 欧美日韩一区二区国产| jizzjizz中国精品麻豆| 久久久久91| 综合激情五月婷婷| 久久在线91|
+

你好，很高興為您服務！

發送