24小时热门版块排行榜

返回列表

心辉

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 837.2
散金: 5
帖子: 87
在线: 33.9小时
虫号: 3505939
注册: 2014-10-29
专业: 有机化学

[求助] greenkubo公式计算导热率时体系没有涨落

greenkubo公式计算导热率时体系没有涨落

回复此楼

» 猜你喜欢

参与限项已经有3人回复
假如你的研究生提出不合理要求已经有7人回复
实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

1楼 2018-08-09 09:33:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

心辉

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 837.2
散金: 5
帖子: 87
在线: 33.9小时
虫号: 3505939
注册: 2014-10-29
专业: 有机化学

unfix  1

compute KE all ke/atom
compute PE all pe/atom
compute V all stress/atom NULL virial
compute J all heat/flux KE PE V

variable Ns equal 10
variable Nc equal 200
variable Np equal ${Ns}*${Nc}*20
variable Np equal 10*${Nc}*20
variable Np equal 10*200*20
fix  3 all ave/correlate ${Ns} ${Nc} ${Np}       c_J[1] c_J[2] c_J[3] type auto file hac.txt
fix  3 all ave/correlate 10 ${Nc} ${Np}       c_J[1] c_J[2] c_J[3] type auto file hac.txt
fix  3 all ave/correlate 10 200 ${Np}       c_J[1] c_J[2] c_J[3] type auto file hac.txt
fix  3 all ave/correlate 10 200 40000       c_J[1] c_J[2] c_J[3] type auto file hac.txt
timestep 0.01
variable M equal 200
variable N equal ${M}*${Np}
variable N equal 200*${Np}
variable N equal 200*40000

run  ${N}
run  8000000
PPPM initialization ...
  using 12-bit tables for long-range coulomb (../kspace.cpp:321)
  G vector (1/distance) = 1.15486
  grid = 90 90 270
  stencil order = 5
  estimated absolute RMS force accuracy = 1.67538e-07
  estimated relative force accuracy = 1.16349e-08
  using double precision FFTs
  3d grid and FFT values/proc = 147500 113400
Setting up Verlet run ...
  Unit style : metal
  Current step  : 200046
  Time step    : 0.01
Per MPI rank memory allocation (min/avg/max) = 38.52 | 38.7 | 38.89 Mbytes
Step Temp Press
  200046 301.13981 7774.4738
  200500 301.13981 7774.4738
  201000 301.13981 7774.4738
  201500 301.13981 7774.4738
  202000 301.13981 7774.4738
  202500 301.13981 7774.4738
  203000 301.13981 7774.4738
  203500 301.13981 7774.4738
  204000 301.13981 7774.4738
  204500 301.13981 7774.4738
  205000 301.13981 7774.4738
  205500 301.13981 7774.4738
  206000 301.13981 7774.4738
  206500 301.13981 7774.4738
  207000 301.13981 7774.4738
  207500 301.13981 7774.4738

赞一下

回复此楼

2楼2018-08-09 09:34:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主心辉的主题更新

返回列表