lammps案例：水分子热导率模拟

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

in文件关键语句已经注释，供大家参考。

# 定义热流
 variable F         	equal 0.