大家好,我是小马老师。
本文分享使用lammps NEMD算法模拟水分子热导率的案例。
关于热导率的模拟,公众号已经有好几个案例分析。
在之前的石墨烯热导模拟案例中,通过郎之万方式输入热量,使用这种方法,需要根据输入的热量计算热流。
lammps还有另外一种输入热量的方式:fix ehex命令。
这个命令可以直接设置热流大小,不需要再画热量曲线求斜率,因此,相对简单。
本文介绍的案例来自lammps官方案例,对应的参考文献为:http://arxiv.org/pdf/1507.07081。
模拟的原理和石墨烯热导一样,在传热方向设置一个热源和一个冷源,在热源持续输入热量,在冷源持续抽走热量,由此产生温度梯度,模拟完成将热流和温度梯度等参数带入公式即可算出热导率。
所用的模型如下图所示,沿z方向进行热传导:
in文件关键语句已经注释,供大家参考。
# 定义热流
variable F equal 0.