漫谈分子动力学计算热导率的五种方法
本仓库提供了一个资源文件,详细介绍了使用分子动力学(Molecular Dynamics, MD)计算材料热导率的五种方法。这些方法包括平衡态的GK法和非平衡态的NEMD/RNEMD(速度重标法、局部热浴法、动量交换法)。资源文件中包含了动态group设置、脚本获取、后处理分析和结果对比等内容,适用于不同工况下的热导率研究。
内容概述
- 平衡态分子动力学(EMD)
- GK法(Green-Kubo方法)
- 方法说明
- 报道结果
- 非平衡态分子动力学(NEMD)
- 速度重标法
- 局部热浴法
- 动量交换法
- 方法说明
- 后处理分析
- 计算结果
- 非平衡态分子动力学(RNEMD)
- 速度重标法
- 局部热浴法
- 动量交换法
- 方法说明
- 后处理分析
- 计算结果
使用说明
- 动态group设置
- 针对流体体系,热源和热汇需要设置为动态group。
- 代码说明
- 以固态60K下的Argon为例进行测试,脚本分为两个:Run in以及几种不同方法。
- 代码获取
- 全部in文件已经开源在博客,可以直接复制。
- Python绘图脚本以及In文件需要打包文件请下载。
后处理分析
- 温度梯度拟合
- 能量输出
- 温度拟合
- 报道结果
结果对比
- 五种方法结果对比
- 注意动态group设置
适用范围
本资源文件适用于需要进行热导率计算的研究人员,特别是使用LAMMPS进行分子动力学模拟的科研工作者。通过本资源文件,用户可以了解并应用不同的热导率计算方法,以满足不同工况下的研究需求。
贡献与反馈
如果您在使用过程中有任何问题或建议,欢迎通过GitHub的Issue功能提出。我们非常欢迎您的反馈和贡献,以帮助我们不断改进和完善本资源文件。