www.毛片在线观看_怡红院亚洲红怡院天堂麻豆_欧美jizz18欧美_国产精品免费久久久久影院

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

    熱導率又稱導熱系數，反映物質的熱傳導能力，其定義為單位溫度梯度（在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使用教程（下篇）

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

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

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

    微信掃碼分享文章
    www.毛片在线观看_怡红院亚洲红怡院天堂麻豆_欧美jizz18欧美_国产精品免费久久久久影院 
    

        

                

欧美亚洲国产一区二区三区|
国产精品丝袜久久久久久app|
国产精品久久久久久久久快鸭
|
国产日韩精品视频一区|
一区二区久久久久|
国产成人免费xxxxxxxx|
中文字幕一区二区三区色视频
|
99久久伊人久久99|
亚洲一区二区三区四区在线观看
|
精品一区二区三区免费|
久久久久久一二三区|
色综合久久天天|
国产亚洲成年网址在线观看|
日韩精品一卡二卡三卡四卡无卡|
92国产精品观看|
婷婷开心久久网|
中文字幕巨乱亚洲|
久久国产成人午夜av影院|
国产女主播视频一区二区|
欧美视频一区二区三区|
国产精品免费视频观看|
国产老妇另类xxxxx|
亚洲美女屁股眼交|
精品国产免费人成电影在线观看四季|
亚洲一区二区影院|
91丝袜呻吟高潮美腿白嫩在线观看|
亚洲国产婷婷综合在线精品|
国产亚洲精品bt天堂精选|
精品一区二区三区蜜桃|
亚洲乱码国产乱码精品精可以看|
日韩精品一区二区在线|
日韩精品欧美精品|
国产精品麻豆久久久|
日韩一区二区三区av|
天天操天天干天天综合网|
久久久www成人免费毛片麻豆
|
成人av动漫网站|
色综合久久久网|
自拍偷拍国产精品|
caoporen国产精品视频|
欧美午夜精品理论片a级按摩|
亚洲精选视频在线|
91麻豆精品视频|
欧美丰满少妇xxxbbb|
亚洲v中文字幕|
中文一区二区在线观看|
欧美不卡一区二区三区|
久久99蜜桃精品|
亚洲综合色在线|
中文字幕中文乱码欧美一区二区|
成人在线视频首页|
欧美性生活影院|
亚洲午夜久久久久久久久电影网|
久久精品夜色噜噜亚洲aⅴ|
7777精品伊人久久久大香线蕉完整版
|
国产精品一区二区在线观看不卡|
午夜精品久久久久久久久|
亚洲天堂中文字幕|
久久九九久久九九|
欧美一级精品在线|
久久精品久久99精品久久|
亚洲一区二区中文在线|
亚洲男人都懂的|
国产欧美一区二区三区在线看蜜臀|
日韩欧美在线123|
国产自产2019最新不卡|
美美哒免费高清在线观看视频一区二区
|
久久精品这里都是精品|
成人免费视频视频|
欧美电影在线免费观看|
免费欧美在线视频|
亚洲成a人片综合在线|
亚洲免费av观看|
欧美国产97人人爽人人喊|
久久九九99视频|
91美女片黄在线|
精品99999|
成人av在线播放网站|
91精品久久久久久久91蜜桃|
狠狠色丁香九九婷婷综合五月|
一本一道波多野结衣一区二区|
亚洲成人自拍偷拍|
亚洲综合免费观看高清在线观看|
又紧又大又爽精品一区二区|
一区在线播放视频|
17c精品麻豆一区二区免费|
国产欧美一区二区三区在线看蜜臀
|
中文字幕在线一区免费|
亚洲国产高清aⅴ视频|
欧美国产1区2区|
国产日韩高清在线|
日本一区二区三区视频视频|
久久色视频免费观看|
久久精品视频一区|
国产亚洲福利社区一区|
国产三区在线成人av|
久久日一线二线三线suv|
久久久精品2019中文字幕之3|
99精品一区二区|
久久久亚洲综合|
久久久国产综合精品女国产盗摄|
久久精品亚洲一区二区三区浴池
|
国产一区不卡在线|
欧美一区二区在线视频|
成人精品鲁一区一区二区|
精品蜜桃在线看|
久久这里只有精品视频网|
欧美韩国日本一区|
国产精品美女久久福利网站|
亚洲精品一二三四区|
亚洲精选视频在线|
香港成人在线视频|
91久久香蕉国产日韩欧美9色|
精品在线视频一区|
欧美一区二区三区四区久久|
波多野结衣中文一区|
国产亚洲一区二区三区|
国产精品无码永久免费888|
亚洲精品免费看|
亚洲大片精品永久免费|
久久国产精品露脸对白|
欧美一级高清片|
久久一留热品黄|
亚洲视频在线观看三级|
亚洲一区二区欧美日韩
|
在线观看免费视频综合|
国产精品自产自拍|
久久综合五月天婷婷伊人|
亚洲国产精品精华液ab|
亚洲五码中文字幕|
欧美伊人久久大香线蕉综合69|
久久久久久久性|
中文字幕欧美一|
日韩精品视频网站|
欧美一区三区二区|
久久久精品欧美丰满|
怡红院av一区二区三区|
91福利精品第一导航|
国产suv精品一区二区三区|
国产三级一区二区|
亚洲女厕所小便bbb|
麻豆91在线播放免费|
日韩你懂的电影在线观看|
欧美国产综合一区二区|
亚洲超丰满肉感bbw|
91精品国产手机|
欧美国产丝袜视频|
首页综合国产亚洲丝袜|
日韩西西人体444www|
国产清纯白嫩初高生在线观看91
|
精品国一区二区三区|
国产精品色哟哟网站|
日本大胆欧美人术艺术动态|
日韩一级二级三级精品视频|
国产蜜臀97一区二区三区|
天天做天天摸天天爽国产一区|
91精品国产丝袜白色高跟鞋|
中文字幕高清不卡|
奇米一区二区三区|
久久久久成人黄色影片|
亚洲综合一区二区|
国产精品一区二区三区乱码|
国产精品毛片高清在线完整版|
午夜精品久久久久久久久久久|
成人黄色免费短视频|
亚洲美女免费在线|
欧美精品一二三四|
国产精品少妇自拍|
激情综合网天天干|
亚洲欧洲成人精品av97|
欧美综合久久久|
久久久精品国产99久久精品芒果|
丝袜美腿亚洲一区|
久久久夜色精品亚洲|
亚洲h动漫在线|
99精品在线免费|
日韩精品亚洲一区|
久久蜜臀中文字幕|
亚洲h精品动漫在线观看|
91视频一区二区|
日韩电影在线观看一区|
国产亚洲福利社区一区|
在线视频一区二区三区|
国产色一区二区|
国精产品一区一区三区mba桃花|
国产精品毛片久久久久久久|
欧美日韩一区三区四区|
国产精品久久久久影院亚瑟|
国产福利精品导航|
一区二区理论电影在线观看|
精品99久久久久久|
91久久国产综合久久|
中文字幕电影一区|
国产福利91精品一区二区三区|
夜夜嗨av一区二区三区网页|
精品国产一区二区在线观看|
一本一道久久a久久精品综合蜜臀|
日本一区二区视频在线观看|
精品无人区卡一卡二卡三乱码免费卡|
日韩理论在线观看|











                
    
    
                
        
                
        
                +
                
    
                

    
    
                
        
                
            
                
                
                
            
                
            
                
                
                你好，很高興為您服務！
                
                
            
                
        
                
    
                

    
    
                
        
        
                
            
                
            
                發送