日韩av老熟女-久久播亚洲视频-国产精品久久久系列-日韩激情96久久-欧美亚洲熟女另类-久久人妻诱惑我-久久高清黄又色视频-日韩av中文字幕网站-日韩人妻精品中文字幕,日韩成人在线不卡av,国产精品一香蕉av,婷婷久久综合久久久

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

熱導(dǎo)率又稱(chēng)導(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ì)開(kāi)發(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ì)算過(guò)程中的模擬盒子示意圖如圖1所示，當(dāng)模擬盒子在某一溫度下平衡后，將模擬盒子中間作為熱端，施加熱流?ε；模擬盒子兩端作為冷端對(duì)熱流進(jìn)行抽取，在抽取的過(guò)程當(dāng)中體系中原子的動(dòng)量是從熱端到冷端不斷交換的，此時(shí)體系處于一個(gè)非平衡的狀態(tài)。因此該過(guò)程也叫非平衡分子動(dòng)力學(xué)(NEMD)過(guò)程。

在施加熱流后，將模擬盒子沿?zé)崃鞣较蚍殖稍S多層，每個(gè)特定的timesteps (如10000步)將每一層的溫度進(jìn)行統(tǒng)計(jì)，當(dāng)傳熱過(guò)程穩(wěn)定后 (這里是指統(tǒng)計(jì)的溫度保持基本穩(wěn)定，一般在加熱流數(shù)萬(wàn)步到十萬(wàn)步后)，對(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ì)算過(guò)程當(dāng)中，LAMMPS的in文件中的傳熱過(guò)程設(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命令取消的是前面作特定溫度下平衡的系綜。其余三段為熱流為傳熱過(guò)程的設(shè)置。

計(jì)算結(jié)束后，最關(guān)鍵的是對(duì)模擬盒子每一層原子溫度的統(tǒng)計(jì)，為此我們寫(xiě)了一個(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)盤(pán)下載：

鏈接：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í)前 四川
文字是人類(lèi)用符號(hào)記錄表達(dá)信息以傳之久遠(yuǎn)的方式和工具?，F(xiàn)代文字大多是記錄語(yǔ)言的工具。人類(lèi)往往先有口頭的語(yǔ)言后產(chǎn)生書(shū)面文字，很多小語(yǔ)種，有語(yǔ)言但沒(méi)有文字。文字的不同體現(xiàn)了國(guó)家和民族的書(shū)面表達(dá)的方式和思維不同。文字使人類(lèi)進(jìn)入有歷史記錄的文明社會(huì)。
點(diǎn)贊12
回復(fù)
全部
查看更多評(píng)論
相關(guān)文章

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

2020-05-03

手把手教你用ChemDraw 畫(huà)化學(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

熱門(mén)文章/popular

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

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

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

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

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

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

微信掃碼分享文章
欧美啊啊啊-91在线精品秘密秘 一区二区-porn日本人-日本一本道大香蕉 | 美女叉叉视频免费在线观看网站-高清欧洲熟妇 720在线手机播放-风满熟妇的名器-5252a我爱haose01亚洲 国产人成视频 | 成人黄色网站在线播放视-BD完整版观看 99热国产在线观看-成年女人毛片免费观看不卡-4k在线观看视频 热久久这里只有精品 | 人妻诱惑蜜臀-国产精品免费久久久久影院 -狠狠干2018最新版本-多P混交群体交乱视频网站导航 | 羞羞动漫成人a片在线观看-高清不卡完整在线观看,美国猛男做爱视频在线观看,91.n女免费在线破处,无水印-69A片-性摔跤xXx69性欧 | 台湾十八成人中文网-一本到视频免费-被C到喷水嗯H厨房交换视频-欧美色图片区 | 一本加勒比苍井空-在线播放“zzz”-中文字幕欧美亚洲-久久无码综合 日本一本免费一区二区三区 | 亚洲av岛国动作片在线观看-夫の上司人妻北条麻妃-激情床戏国产精品久久茄子影视-夜夜噜2019 | Chinese丰满东北熟女人妻-深喉影院导航-美女穿黑丝自慰-欧美一页9 | 亚洲午夜国产精品无码老牛影视-美女又爽 又黄 免费孕妇-夜本色狼人AV-一区二区三区AV夏目彩春 | EEUSS鲁丝片Av无码安全直达-国产真实偷乱视频在线观看-18禁美女裸体自慰网站-偶偶福利电影 | 秋霞EEUSS-国产午夜精品自在自线之la-约粉嫩女神美女在线-老王影院AV | 国产剧免费观看 国产午夜福利精品一区二区三区-欧美色图-1-271重口女m-搡老熟女-91Porn | 真人逼环AV-嫩白女特级黄片-国产在线精品二区韩国演艺界-日韩动态欧美 | 78摸成人视频-jjzzcom日韩-国产a级理论片免费网站-在线观看精品一区 | 一本热久-c逼好爽-国产丰满乱子伦无码专区-中国美女新婚之夜高潮 | 黄频网站大全-不卡的影视大全-干日本少妇-欧美淫妇 | 国产91人妻绿帽黑人-美女自慰网址-欧美丰满熟妇BBBBBB禁忌-扶她射成人 | 北条麻妃VS黑人巨大-久久偷偷做嫩草影院免费看-一久本道-DVD影视在线观看 中文字幕人妻无码乱精品 | 亚洲sss视频在线-XNXXChina国产实拍-国产同性人妖TS口直男-二人扑克剧烈运动软件直播下载 | 超碰1997-一本之道高清视频在线观看-亚洲精品手机在线观看电影-蜜美杏AV无码破解在线 | 国产在线观看免费播放电视剧美女图片-一区二区三区中文字幕脱狱者-91色久-gogogo高清免费观看完整版 | 国产精品制服一区二区-日本性久久-四季AV一区二区夜夜嗨-720免费全集观看 亚洲综合无码久久精品综合 | 黑人亚洲成人淫奸播放-涩av网站-Xxxxxxx色色色色电影Xxxxx-gogogo手机高清在线观看软件 | 亚洲爆乳无码一区二区三区-Chinese国产av-国产农村少妇性事-曰本毛茸茸性生话 | 一个人嗯啊视频-国产乱伦电影艳照门-公交车挺进朋友人妻的身体里-baomaav 国精品无码一区二区三区在线蜜臀 | ChinaAV在线-china熟女熟妇乱老女人-多人交换做爰g.gav国产-国产日韩欧美一区二区张家界 | 电影免费在线观看-干丰满老熟女视频-欧洲人玩亚洲人XXX-国产成人免费无码?v在线播放 | 国产又粗又黄又爽又硬P站-日韩在线 np-欧美日韩淫荡女人-3B肉蒲团之四虎影库之嫩嫩草 | www.少妇视频-久久国色Av无码免费看-开心情涩网-国产美女h在线 | 中文字幕精品一区二区绿巨人-裸体美女爱爱-av小说电影-色情A∨色人阁 | 亚洲精品国偷自产-牛夜精品久久久久久久99黑人-高清全集免费观看 在线观看一级毛片免费-国产永久免费看黄网站 | 一本加勒比苍井空-在线播放“zzz”-中文字幕欧美亚洲-久久无码综合 日本一本免费一区二区三区 | 熟女人妻精品猛烈进入-三级电影网站多少-亚洲第一色日韩激情在线-四虎884a | 亚洲你XX我XX网站-怡红院成永久免费人视频新的-肏屄毛片啊啊啊-亚洲精品久久久久毛片A级绿茶 | 被npc们啪哭高H(总受)-51看nba高清直播视频-亚洲国产日韩精品一区二区-台湾啪啪网 | wwwtubecom69日本高潮-超碰护士-妖精视频app在线观看ios-水野优香在线一区二区 | 香港,日韩,AV。-好骚好紧15p-91超级碰-日韩视频一区二区 | JIZZJIZZJIZZ亚洲熟妇-国产一二三区在线播放-伊人精品影院-www.se美女不卡 | 又粗又长又大又粗又硬毛片大香蕉-欧美黄视频-gogogo高清在线观看免费播放电视剧-一道本视频在线观看 | 午夜理伦精品一区二区三区在线成人-┅┅用力┅┅老熟女-国产精品三上悠亚在线观看网站-东北电影免费老师无码黄色片视屏网站大全 亚洲美女自慰-七十路熟女视频-日韩小穴-色女人综合 |
+

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

發(fā)送