日韩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以及溫度梯度?,而J和?具有公式(1)的線性關(guān)系,比例系數(shù)就是材料的熱導(dǎo)率。

J = -?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)伏安法詳解(下)

微信掃碼分享文章
操操操人人乐人人干-gogogo手机免费高清-黑人初解禁薍伦黑人巨大在线观看-人妻av中字 中文字幕 - 精工厂 | 懂色av蜜臀av粉嫩avdnlt-BD电影免费在线 精品久久久久久中文字幕网 -99国产精品国产精品久久-一女n男4p喷到她哭还流着白浆 | 羞羞动漫成人a片在线观看-高清不卡完整在线观看,美国猛男做爱视频在线观看,91.n女免费在线破处,无水印-69A片-性摔跤xXx69性欧 | 樱花影视 免费高清电影电视剧手机版在线观看-玩弄丰满人妻-极品美女想被c-2020最新国产精品视频 黑森林毛茸茸熟女人妻-女人自慰网址-精品国产无套在线观看-日本睡熟迷奷系列A片 | 爱妃在线二区-西西444WWW无码视频男男-日韩经典第三页-中国桃子被艹的视频 | 77777色婷婷-中国熟女自慰DH-中文字幕一本一道在线看-欧美精品 日韩伦理 | 91麻豆糖心vlog出品淑怡熟女亚洲一区精品久久-国产精品一区无码专区-欧美V亚洲V日韩v国产Vv∨-日韩精品久久久久久久电影99爱 | 熟女综合区图片-偷拍第14页-干mm色啊视频-老师日批视频 | 21p少妇-美女100%裸奶无遮挡免费国产又黄又大-国内精品第21页-日韩人妻无码精品久久免费 | 猫咪大香蕉手机在线视频-中文天堂最新中文字幕版-DVD手机免费观看完整 漂亮的保姆中文版-粉嫩AV四季AV绯色AV | 日韩三级片网站一区二区-大屁股,另类HD-视频一区中文字幕种子连接-色橹aaaa | 人妻少妇精品专区性色av-国产精品无码专区第1页-粉嫩videos好紧-看看东北老区老熟女肥臀搡逼激情 | 中国美女处内谢-大陆欧美日韩-沟厕精品清晰女厕正面-成年免费在线视频 乳欲人妻奶水 | 久久国产露脸老熟女熟69-亚洲xxxxxxx-97人妻精品一区二区三区水池里-久久伊人国产精品 | 俺去射婷婷-小婴交videossex-中文伊人-欧美色妇15p | 黑人玩弄人妻1~12-一区二区三区四区中文字幕在线观看-剧情片全集免费播放 第一种无激情码激情 HD电影在线观看 -日韩啪啪视频 吹潮喷水白浆在线播放-亚洲美女插-99久久精品视香蕉蕉-黑人60厘米巨大全进入 | 亚洲人XXX日本人com-91伊人影院-女人免费网站-99久久免费看少妇高潮a片特黄 | 大骚逼1000p-97影视大全免费追剧大全在线观看-久久国产av无码-乱欲小说亚洲专区 | 青青青欧美在线观看视屏-四虎熟女-五十路AV女优电影-人妻理论片 台湾成人综合-久久九九九字幕-BBBWWW-国产乡下妇女做爰视频 | 亚洲欧美日韩在线-高级人妻上门北条麻妃-在线亚洲看黄色电影 japonensis18一25-欧美性爱男人天堂按摩大师 | 国产精品裸女-极品人妻波多野结衣瑜伽裤-黑人解禁-六十路老熟女交尾 | 熟女色竹菊影视-欧洲美妇乱人伦视频网站-玩弄丰满奶水的老师喷流白浆国内,美女-BD英语免费观看韩国 | 国产男女无套 免费漫画-97无码免费人妻超级碰碰夜夜-国产精品一区久久-色射综合 | 无码前列腺按摩XXXXXXav-色小姐一本道-超碰美女自慰-少妇爱艹无码 | 精品一区二区不卡无码av-欧美99999-JIZZ8性黑人-日韩精品射精管理在线观看 | 亚洲第九页夜-(高H,高c)开荤黑人-亚洲女同色情-五十路老熟道中出在线播放 | 亚洲欧美v国产蜜芽tv-蓝光免费观看 国产精品自产拍av在线-国产WjjZZ-三上悠亚青苹果乐园影视 | 欧美a色爱欧美综合v-.asian美白裸体女pics-水蜜桃4美国伦理-日韩毛片AV-女同自慰A片免费观看-森泽佳奈持续高潮40分钟 | 高清电影观看 欧美日韩国产一区二区三区地区-波多野结衣同性女恋大片-欧美专区15页-真实亲伦对白清晰在线播放 | gogogo高清国语播放-床戏做爰呻吟声-美女自卫慰免费视频www免费-少妇毛茸茸aaahd极品 | 国产盗摄xXxX视频XXXⅩ-男女肉粗暴国产av-中文字幕探花在线-成人永久免费视频网站app | 亚洲国产精品网站久久-国产AV一区二区三区-无码国产精品一区二区免费i6-性欧美XXX1819内谢 | 德国熟妇超碰在线电影院网站-澳门日本亚洲男同666sex-亚洲高潮视频-国偷自拍视频在线观看 | 日韩免费高清专区-黄色ab节视频-国产视频这里有精品视频-精品一区二区三区免费观看 | 亚洲精品国偷自产-牛夜精品久久久久久久99黑人-高清全集免费观看 在线观看一级毛片免费-国产永久免费看黄网站 | 熟女综合区图片-偷拍第14页-干mm色啊视频-老师日批视频 | 女上位打桩榨精在线观看-高潮毛片无遮挡免费高清风月直播-尹人成人-涩综合婷婷久久涩 | 美女自衛免費看hhh-一本道专区国产-亚洲啪啪综合?v一区综合精品区-欧美毛基地 | 国产一区二区高清-淫荡骚穴-精品国产一区二区三区久久狼黑人-1080日本电影免费 8164 9628 小蛮腰 网页无法访问 | 四虎图库-娇小黑白配大战黑人-看一级肥胖女人毛片-chinesegaysextubevideose | 亚洲欧美日韩在线-高级人妻上门北条麻妃-在线亚洲看黄色电影 japonensis18一25-欧美性爱男人天堂按摩大师 |
+

你好,很高興為您服務(wù)!

發(fā)送