24小时热门版块排行榜

返回列表

雪雨星风17

木虫 (小有名气)

无

应助: 1 (幼儿园)
金币: 5327.9
散金: 76
红花: 2
帖子: 236
在线: 371小时
虫号: 1368355
注册: 2011-08-14
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 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

回复此楼