| 查看: 1463 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用Green-Kubo自相关函数计算热导率，lammps如何引用Green-Kubo公式，in文件应如何写
已经有4人回复
- lammps求热导率的问题
已经有8人回复
- 用lammps模拟填充水的碳管的热导率，能设置碳管为固定边界，水为周期性边界条件吗？
已经有6人回复
- 用LAMMPS模拟金属热导率时要考虑电子的贡献吗？
已经有6人回复
- lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
已经有7人回复
2楼2016-01-19 19:23:38
投票:
10