24小时热门版块排行榜    

查看: 8125  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:24 by 9i7PCSZyMzot
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:12 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:12 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:06 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:01 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 07:01 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 4/200 2026-09-20 06:55 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
信息提示
请填处理意见