网站地图 关于 合作 招聘 帮助 会员中心

    热线:400-152-6858

    测试狗科研服务

    预存 定制
    Document
    当前位置:文库百科文章详情
    分子动力学(LAMMPS)——实例讲解:NEMD方法计算晶格热导率
    来源:科学10分钟 时间:2021-05-18 17:58:03 浏览:6642次

    热导率又称导热系数,反映物质的热传导能力,其定义为单位温度梯度(在1m长度内温度降低1K)在单位时间内经单位导热面所传递的热量。材料的热导率是由电子热导率晶格热导率共同贡献的。在温度较高的情况下,材料的热导率主要由电子热导率主导;在温度较低的情况,主要由晶格热导率主导。

    实验中测材料热导率是在材料两端施加热流,此时可以测得沿热流方向的横截面上的热流通量J以及温度梯度,而J和具有公式(1)的线性关系,比例系数就是材料的热导率。

    J = -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
    回复
    全部
    查看更多评论
    相关文章

    晶体结构可视化软件 VESTA使用教程(下篇)

    2021-01-22

    【科研干货】电化学表征:循环伏安法详解(上)

    2019-10-25

    手把手教你用ChemDraw 画化学结构式:基础篇

    2021-06-19

    【科研干货】电化学表征:循环伏安法详解(下)

    2019-10-25

    晶体结构可视化软件 VESTA使用教程(上篇)

    2021-01-22

    电化学实验基础之电极基础篇 (四)参比电极的使用和养护

    2020-08-21

    热门文章/popular

    晶体结构可视化软件 VESTA使用教程(下篇)

    【科研干货】电化学表征:循环伏安法详解(上)

    电化学实验基础之电化学工作站篇 (二)三电极和两电极体系的搭建 和测试

    手把手教你用ChemDraw 画化学结构式:基础篇

    【科研干货】电化学表征:循环伏安法详解(下)

    晶体结构可视化软件 VESTA使用教程(上篇)

    微信扫码分享文章

    意见反馈

    有奖举报

    商务合作

    公众号

    关注我们 了解更多

    小程序

    随时预约 掌握进度

    举报有奖

    TEL: 191-3608-6524

    如:在网络上恶意使用“测试狗”等相关关键词误导用户点击、恶意盗用测试狗商标、冒称官方工作人员等情形,请您向我们举报,经查实后,我们将给予您奖励。

    举报内容:

    200

    上传附件:
    文件格式不正确,请重新上传文件格式不正确,请重新上传文件格式不正确,请重新上传
    文件格式:jpg、jpeg、png、gif、tif、doc、docx、ppt、pptx、xls、xlsx、pdf、zip、rar
    联系方式
    姓名
    电话
    提交意见

    意见反馈

    Suggestions

    您可以在此留下您宝贵的意见,您的意见或问题反馈将会成为我们不断改进的动力。

    意见类型
    测试服务
    网站功能
    财务报账
    其他类型
    意见内容

    200

    联系方式
    姓名
    电话
    提交意见

    收起

    Document
    关于我们 新手帮助 测试干货 商务合作 基金查询 相关资质 模拟计算 现场测试 服务项目 科研绘图 同步辐射 电池行业

    联系方式/contact

    400-152-6858

    工作时间/work time

    09:00-18:00

    测试狗公众号

    关注我们 了解更多

    测试狗小程序

    随时预约 掌握进度

    蜀公网安备51010602000648号

    蜀ICP备17005822号-1

    成都世纪美扬科技有限公司

    Copyright@测试狗·科研服务