| 查看: 1385 | 回复: 1 | |||
雪雨星风17木虫 (小有名气)
无
|
[求助]
lammps MP法求热导率,为何体系的温度下降了。
|
|
采用lammps的MP方法求热导，关键部分代码如下： units metal variable T equal 300 velocity all create $T 458127641 mom yes rot yes dist gaussian units box fix 1 all temp/berendsen $T $T 0.001 compute ke all ke/atom variable temp atom c_ke/(1.5*1.3806504e-4)*1.602763 fix temp_profile all ave/spatial 1 100000 100000 z lower 2 v_temp file temp.profile units lattice thermo_style custom step temp etotal vol thermo_modify lost warn thermo 100 # Run dump id all atom 10000 dump.lammpstrj timestep 0.001 run 100001 unfix 1 fix nve all nve fix heat_swap all thermal/conductivity 10 z 400 fix e_exchange all ave/time 10 10000 100000 f_heat_swap file e_exchange.dat run 10000000 开始原子的温度设定为300K,热浴之后得到的结果 99600 300 -132337.22 12788141 99700 300 -132337.22 12788141 99800 300 -132337.22 12788141 99900 300 -132337.22 12788141 100000 300 -132337.22 12788141 100001 300 -132337.22 12788141 紧接着进行热导计算的结果 100001 300 -132337.22 12788141 100100 214.80808 -132337.05 12788141 100200 198.52124 -132336.87 12788141 100300 221.0091 -132336.93 12788141 100400 264.85103 -132336.99 12788141 100500 270.57875 -132336.95 12788141 100600 304.28282 -132337.09 12788141 在连续的波动之后，温度在 246K左右稳定了，请问大神这是什么原因？ 445400 246.03241 -132336.74 12788141 445500 245.74197 -132336.74 12788141 445600 246.37433 -132336.74 12788141 445700 244.57574 -132336.74 12788141 445800 246.93508 -132336.74 12788141 445900 245.90909 -132336.74 12788141 446000 246.30204 -132336.74 12788141 446100 246.08624 -132336.74 12788141 446200 244.93367 -132336.74 12788141 446300 245.27753 -132336.74 12788141 446400 245.15068 -132336.74 12788141 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
用Green-Kubo自相关函数计算热导率，lammps如何引用Green-Kubo公式，in文件应如何写
已经有4人回复
- lammps求热导率的问题
已经有8人回复
- 用lammps模拟填充水的碳管的热导率，能设置碳管为固定边界，水为周期性边界条件吗？
已经有6人回复
- 用LAMMPS模拟金属热导率时要考虑电子的贡献吗？
已经有6人回复
- lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
已经有7人回复
2楼2016-01-19 19:23:38
