NEMD方法计算热导率之langevin控温法

阿磊的MD和CFD记录簿

已于 2022-04-30 12:15:01 修改

阅读量6.5k

点赞数 7

分类专栏： lammps学习总结文章标签： lammps

于 2021-11-02 09:00:43 首次发布

本文链接：https://blog.csdn.net/weixin_45896923/article/details/121093114

版权

lammps学习总结专栏收录该内容

17 篇文章 88 订阅

订阅专栏

NEMD方法计算热导率之langevin控温法

热导率反映物质的热传导能力，是影响材料热传输性质的关键参数，在纳米流体、相变材料等领域都需要对新型材料的热导率进行深入的研究。近代的观点把这种能量传输归因于原子运动导致的晶格波造成的。在非导体中，能量传输只依靠晶格波进行；在导体中（比如银、铁），除了晶格波还有自由电子的平移运动。Lammps一般计算的是晶格振动产生的声子热导率贡献，计算方法有EMD、NEMD以及RNEMD方法，NEMD又可分为加减能量法以及固定温度法，由于固定温度法操作简单，易于理解，因此这里介绍一下使用Langevin控制固定温度统计热流并构造温度梯度的NEMD方法。

NEMD方法基于傅里叶定律（公式1），其反应出热流密度（单位W·m-2）、温度梯度以及热导率之间的关系。

在这里插入图片描述
在NEMD中，Langevin恒温器与Nosé-Hoover恒温器相比具有随机性和局部局域性的优势。此外，对于不对称的碳基结构，Nosé-Hoover恒温器可能会产生导致不切实际的热整流的假象。因此，Langevin控温器是一种更好的选择。
Langevin控温计算热导率的思路（以周期性边界、热流流向Z方向为例）：
（1）模拟系统使用NVE系宗，体系沿Z方向分为20层。
（2）使用langvine控温方法将模拟系统的底部（第1层）温度固定在指定的较高值Thi，中间（第11层）固定在指定的较低温度Tlo。
（3）使用fix langevin中tally yes命令统计原子上加入/减去的能量（f_hot和f_cold）
（4）可根据公式2来计算热导率。
在这里插入图片描述

在这里插入图片描述

最后以lammps自带的例子KAPPA文件中的in.langevin为例，代码不再附上。
详细的计算过程如下：

原子数N=8000
总转移动能=8000*(f_hot+f_cold）/2=8000*（0.907+0.968）/2=7500
计算时间t为0.005*20000=100
A=Lx*Ly=18.82*18.82
温差平均=f_ave=0.57
温度梯度=温差/(Lz/2)=f_ave/lz=0.57/18.82
热导率=7500/2/100/18.82/0.57=3.49

阿磊的MD和CFD记录簿

关注

7
点赞
踩
37

收藏

觉得还不错? 一键收藏
打赏
9
评论
NEMD方法计算热导率之langevin控温法

NEMD方法计算热导率之langevin控温法热导率反映物质的热传导能力，是影响材料热传输性质的关键参数，在纳米流体、相变材料等领域都需要对新型材料的热导率进行深入的研究。近代的观点把这种能量传输归因于原子运动导致的晶格波造成的。在非导体中，能量传输只依靠晶格波进行；在导体中（比如银、铁），除了晶格波还有自由电子的平移运动。Lammps一般计算的是晶格振动产生的声子热导率贡献，计算方法有EMD、NEMD以及RNEMD方法，NEMD又可分为加减能量法以及固定温度法，由于固定温度法操作简单，易于理解，因此这
复制链接

扫一扫