24小时热门版块排行榜    

查看: 827  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 心辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +4 于海海海海 2026-03-19 4/200 2026-03-19 11:11 by kakakapanpan
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 321求调剂 +6 何润采123 2026-03-18 6/300 2026-03-19 10:55 by daliangogo
[考研] 271材料工程求调剂 +5 .6lL 2026-03-18 5/250 2026-03-19 03:07 by 无懈可击111
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 0703化学调剂 ,六级已过,有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
信息提示
请填处理意见