24小时热门版块排行榜    

查看: 7756  |  回复: 2

cl1994861

金虫 (小有名气)

[求助] 关于fix ave/chunk command问题 已有1人参与

在运行lammps中examples 中的例子KAPPA中in.mp时出现  ”error :no value in fix ave/chunk command",我都没改它的代码,只是照着运行一下,这个例子的in文件的文本如下:
echo screen
variable        x equal 10
variable        y equal 10
variable        z equal 20

variable        rho equal 0.6
variable        t equal 1.35
variable        rc equal 2.5

units                lj
atom_style        atomic

lattice                fcc ${rho}
region                box block 0 $x 0 $y 0 $z
create_box        1 box
create_atoms        1 box
mass                1 1.0

velocity        all create $t 87287

pair_style        lj/cut ${rc}
pair_coeff        1 1 1.0 1.0

neighbor        0.3 bin
neigh_modify        delay 0 every 1

# 1st equilibration run

fix                1 all nvt temp $t $t 0.5
thermo                100
run                1000

velocity        all scale $t

unfix                1

# 2nd equilibration run

compute                ke all ke/atom
variable        temp atom c_ke/1.5

fix                1 all nve

compute layers all chunk/atom bin/1d z lower 1 units reduced
fix                2 all ave/chunk 10 100 1000 layers v_temp file profile.mp

fix                3 all thermal/conductivity 10 z 20

variable        tdiff equal f_2[11][3]-f_2[1][3]
thermo_style        custom step temp epair etotal f_3 v_tdiff

thermo                1000
run                20000

# thermal conductivity calculation
# reset fix thermal/conductivity to zero energy accumulation

fix                3 all thermal/conductivity 10 z 20

fix             ave all ave/time 1 1 1000 v_tdiff ave running
thermo_style        custom step temp epair etotal f_3 v_tdiff f_ave

run                20000

关于fix ave/chunk  command问题
M5B_V5XE5VZ9S0C`_7[8~U2.png
回复此楼

» 猜你喜欢

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

小小波bi

铁虫 (初入文坛)

楼主,compute chunk/atom 这个命令的问题解决了吗?
2楼2016-11-15 21:31:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2828512

铁虫 (初入文坛)

【答案】应助回帖

compute layers all chunk/atom bin/1d z lower 1 units reduced
这条命令中的1改成0.02,因为这条命令的意思是将all沿着z轴分成许多份,而你如果写的是1的话就是只有一份,应该输入小于1的值,例如输入0.02就是分成了1/0.02=50份
3楼2023-03-22 17:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cl1994861 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 生物学一志愿985,分数349求调剂 +3 zxts12 2026-03-21 3/150 2026-03-21 16:34 by 33来了真来了
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 0856材料专硕353求调剂 +3 NIFFFfff 2026-03-20 3/150 2026-03-21 10:23 by luoyongfeng
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
信息提示
请填处理意见