24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 7783  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +4 关忆北. 2026-04-14 4/200 2026-04-15 09:20 by heqizheng
[考研] 求助调剂,跨调 +18 X十甫寸Y 2026-04-11 19/950 2026-04-14 19:26 by Art1977
[考研] 药学305求调剂 +10 玛卡巴卡boom 2026-04-10 10/500 2026-04-14 15:55 by zs92450
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 人工智能320调剂08工类还有机会吗 +18 振—TZ 2026-04-10 19/950 2026-04-14 10:34 by screening
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 求调剂,262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 296求调剂 +14 汪!?! 2026-04-08 15/750 2026-04-11 20:28 by dongdian1
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
信息提示
请填处理意见