日本久久久久久一区-国产麻豆乱视频av-日韩黄色一级免费电影-精品久久久久久久二区-久久亚洲精品尤物-久久久久偷看国产亚洲87-日韩精品资源在线观看,-久久草视频免费观看-亚洲少妇久久嫩13,91精品福利小视频,日韩欧美熟妇在线播放,亚洲成熟女人毛茸茸

預(yù)存
Document
當(dāng)前位置：文庫(kù)百科 ? 文章詳情
分子動(dòng)力學(xué)(LAMMPS)——實(shí)例講解：NEMD方法計(jì)算晶格熱導(dǎo)率
來源：科學(xué)10分鐘 時(shí)間：2021-05-18 17:58:03 瀏覽：14307次

熱導(dǎo)率又稱導(dǎo)熱系數(shù)，反映物質(zhì)的熱傳導(dǎo)能力，其定義為單位溫度梯度（在1m長(zhǎng)度內(nèi)溫度降低1K）在單位時(shí)間內(nèi)經(jīng)單位導(dǎo)熱面所傳遞的熱量。材料的熱導(dǎo)率是由電子熱導(dǎo)率和晶格熱導(dǎo)率共同貢獻(xiàn)的。在溫度較高的情況下，材料的熱導(dǎo)率主要由電子熱導(dǎo)率主導(dǎo)；在溫度較低的情況，主要由晶格熱導(dǎo)率主導(dǎo)。

實(shí)驗(yàn)中測(cè)材料熱導(dǎo)率是在材料兩端施加熱流，此時(shí)可以測(cè)得沿?zé)崃鞣较虻臋M截面上的熱流通量J以及溫度梯度?T ，而J和?T 具有公式（1）的線性關(guān)系，比例系數(shù)K 就是材料的熱導(dǎo)率。

J = -k ?T （1）

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

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


圖1. NEMD計(jì)算熱導(dǎo)率的三維模型示意圖

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

 在計(jì)算過程當(dāng)中，LAMMPS的in文件中的傳熱過程設(shè)置腳本段如下：

unfix            1

fix              2 all nve

fix              hot all heat 10 ${power} region hot

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

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

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

最后，介紹幾個(gè)可以使結(jié)果更優(yōu)的方法：

（1）將模擬盒子三個(gè)方向邊長(zhǎng)，建成W單包三個(gè)方向邊長(zhǎng)的倍數(shù)，這樣模擬盒子具有完美晶體的周期性。

（2）將采集溫度的模擬時(shí)間盡量增大，這樣能夠采集更多的樣本，統(tǒng)計(jì)的溫度就會(huì)更加準(zhǔn)確。

該算例以及溫度梯度處理腳本可以在網(wǎng)盤下載：

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

提取碼：CSGO

參考文獻(xiàn)：

[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.
文章精選 


你的點(diǎn)贊和在看，我都認(rèn)真當(dāng)成喜歡
評(píng)論 / 文明上網(wǎng)理性發(fā)言
12條評(píng)論
全部評(píng)論 / 我的評(píng)論
最熱 /  最新
全部 3小時(shí)前 四川
文字是人類用符號(hào)記錄表達(dá)信息以傳之久遠(yuǎn)的方式和工具。現(xiàn)代文字大多是記錄語(yǔ)言的工具。人類往往先有口頭的語(yǔ)言后產(chǎn)生書面文字，很多小語(yǔ)種，有語(yǔ)言但沒有文字。文字的不同體現(xiàn)了國(guó)家和民族的書面表達(dá)的方式和思維不同。文字使人類進(jìn)入有歷史記錄的文明社會(huì)。
點(diǎn)贊12
回復(fù)
全部
查看更多評(píng)論
相關(guān)文章

基礎(chǔ)理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結(jié)構(gòu)）

2020-05-03

手把手教你用ChemDraw 畫化學(xué)結(jié)構(gòu)式：基礎(chǔ)篇

2021-06-19

晶體結(jié)構(gòu)可視化軟件 VESTA使用教程（下篇）

2021-01-22

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（上）

2019-10-25

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（下）

2019-10-25

XRD的基本原理與應(yīng)用

2020-11-03

熱門文章/popular

基礎(chǔ)理論丨一文了解XPS（概念、定性定量分析、分析方法、譜線結(jié)構(gòu)）

手把手教你用ChemDraw 畫化學(xué)結(jié)構(gòu)式：基礎(chǔ)篇

晶體結(jié)構(gòu)可視化軟件 VESTA使用教程（下篇）

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（上）

電化學(xué)實(shí)驗(yàn)基礎(chǔ)之電化學(xué)工作站篇 （二）三電極和兩電極體系的搭建 和測(cè)試

【科研干貨】電化學(xué)表征：循環(huán)伏安法詳解（下）

微信掃碼分享文章
亚洲精品,欧美日韩-蜜臀av性久久久久蜜臀aⅴ-日韩殴美在线视频-91大神天堂福利在线 | 99的爱精品免费视频-国产av一区二区3区-超91福利在线观看-日韩中文字幕人妻免费 | 久久综合综合久久22色-免费精品国自产拍在线不卡-日韩专区 在线观看-熟女av一区二区三区 | 精品少妇一区二区三区-国产美女av一区二区三区-中文字幕亚洲天堂中文字幕-中文字幕久久j日 | 久久91精品久久久久久水蜜桃-日本一本东京热久久-国产精品久久久久999-蜜臀精品国产91久久 | 日韩熟女在线播放电影-亚洲精品国自产在线-91精品 粉嫩 久久-一本一道久久综合狠狠老_ | 日韩一区二区三区欧美-欧美日韩激情在线观看-亚洲精品99久91在线水蜜桃-40路熟女在线观看 | 91久久精品国产原创91-久久久久熟妇熟女-久一久久久久久久久-人妻系列中文字幕大乳丰满人妻 | 999国产精品久久久久-国产超碰人人做人人爱-日韩av一级大片在线观看-久久草草视频在线 | 国产欧美日韩综合区-性色av一区二区三区蜜臀-成人好爽网站免费毛片av-日韩欧美美女一区二区三区 久久男人的天堂成人-国产一区二区在线红桃-久久精品国产99久久99久久免费-超碰97人人做人人爱亚洲尤物 | 91九色蝌蚪蜜桃臀-亚洲熟妇av在线免费观看-中文字幕在线高清视频-热久久综合在线观看 | 国产成人精品99久久久-婷婷久久av综合一区二区三区-狠狠做婷婷激情四射-国产成人精品视频a | 国模一区二区三区四区视频-999国产精品亚洲-日本久久久三级-玩弄超骚少妇一区二区三区 | 99精品视频在线观看全部-狠狠久久之综合久久-大胆丰满熟妇一区二区三区-99热99re在线播放 | av岛国片免费在线观看-五月六月婷婷综合激情-久久国产综合精品五月天-国产麻豆精品人妻av | 午夜激情视频资源网-亚洲视频在线观看一二三四区-精品日韩视频在线看-色综合久久久久久久久8噜啦噜 | 99热成人精品热久久6网站-亚洲人成电影在线网站-欧美中文字幕成人在线观看-日韩免费精品视频在线 | 亚洲少妇久久嫩13-国产在线成人小视频-91亚洲精品99久91在线传媒-蜜臀av噜噜一区二区三区 | 成人精品在线免费视频-日韩av中文字幕在线播放免费观看-久久人妻少妇中文字幕少妇-免费日韩资源福利一区二区视频 | 久久男人的天堂成人-国产一区二区在线红桃-久久精品国产99久久99久久免费-超碰97人人做人人爱亚洲尤物 | 亚洲欧美日韩永久在线-中文字幕在线看日本-精品视频在线你懂得-色婷婷在线视频在线观看 | 亚洲欧美另类色婷婷-亚洲欧美熟女一区-婷婷激情五月天影院-四季av人妻精品一区二区三区 | 久久精品国产亚洲av久野外-欧美精品亚洲精品日韩精品-日韩人妻中文字幕一区二区久久-日韩人妻精品字幕免费 | 亚洲精品成av人在线观看夫-亚洲va在线va天堂va第1集-欧美日韩中国一区久久-欧美,日韩,亚洲,丝袜 日韩欧美在线视频观看免费-亚洲av乱码毛片在线播放-欧美日韩国产自拍av-91久久国产综合久久久久 | 天天干天天摸天天操天天操-2012中文字幕在线资源-日韩一区二区欧美一区-97久久久精品综合88久久 国产91色综合久久麻豆-999久久久精品免费-日韩美女a级片免费视频-91久久久久久久久精 | 久久久久久亚洲精品熟女-91成人午夜在线视频-亚洲欧洲avtt在线-国产精品久久国产精麻豆99网站 | 婷婷激情久久综合-午夜精品久久婷婷-精品一区二区三区四区在线视频-麻豆精品传媒一二三区蜜桃 | 99精品,国产精品-中文字幕一区二区三区四区精品-一区二区美女现场直播-日韩蜜乳一区二区三区 | 国产人妻福利一区二区三区-国产日韩亚洲中文字幕-美女视频一区二区三区四区-日韩熟女www | 91精品久久久aaaaa-久久精品人人爽a∨-欧美高清…一区二区三区-91精品国产综合久久久久婷婷 | 日本五十路熟女xx-久久99人妻免费精品一区-91精品久久久久久久亚洲国产-中文字幕av在线人妻 | 久久视频成人国产91-六月丁香婷久久-亚洲永久精品日韩成人av-久久最近最新高清中文字幕 | 无套内射高清视频-91久久99久久91熟女-久久久99婷婷久久久久久-午夜剧场成人自拍亚洲品质 | 婷婷开心网狠狠爱-国产精品 日韩精品 中文字幕-国产一区二区三区四区五区加勒比-97超碰在线观 | 久久久久精品中文字幕-亚洲熟女一区二区三区在线观看-国产精品中文字幕在线不卡-精品少妇一区二区三区在线播放 | 中文字幕一区二区蜜桃-ww久久综合久中文字幕-精品久久久久久久密桃-激情综合网五月天五月激情 | 日韩翘臀诱惑美女主播-99精品视频这里免费看-蜜臀久久99精品久久久老-久久精品三级视频 | 成人激情自拍视频在线观看-精品毛片av一区二区三区-久久精品国产亚洲av亚-国产麻豆黄色一级视频 国产99视频精品免视看9-激情综合网激情五月网-成人av在线不卡一区-蜜臀久久99精品,fv | 日本久久高清免费福利-久久精品国产99久久6动漫-精品人妻一区二区专区-都市激情 中文字幕 校园春色 | 久久久人妻在线观看-av中文字幕乱码-黑人 hdhdhd-久久亚洲国产一区二区三区 | beeg老熟妇喷水-av黄色一级片在线观看第二区-亚洲av日韩综合一区在线观看-婷婷www久久 |
+

你好，很高興為您服務(wù)！

發(fā)送