24小时热门版块排行榜    

查看: 828  |  回复: 1

心辉

铜虫 (小有名气)

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

greenkubo公式计算导热率时体系没有涨落
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

心辉

铜虫 (小有名气)

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的回帖
相关版块跳转 我要订阅楼主 心辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856调剂,是学校就去 +3 sllhht 2026-03-19 4/200 2026-03-19 16:55 by allen-yin
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 274求调剂 +6 S.H1 2026-03-18 6/300 2026-03-19 09:34 by 花店相见
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
信息提示
请填处理意见