24小时热门版块排行榜    

Znn3bq.jpeg
查看: 7870  |  回复: 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一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:55 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:30 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:15 by ue3ir18jc3
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:10 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 xx7gd5zq4e 2026-05-15 5/250 2026-05-16 23:57 by ue3ir18jc3
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 23:10 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +6 水和泥不是水泥 2026-05-13 8/400 2026-05-16 19:50 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 3/150 2026-05-16 17:11 by zznnnj
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[硕博家园] 考博自荐 +3 科研狗111 2026-05-13 4/200 2026-05-16 11:45 by 科研狗111
[文学芳草园] 裁员滚滚,退居二线 +4 J_wei 2026-05-10 4/200 2026-05-16 10:52 by zh10246
[考博] 西南大学考核制博士 +4 lijunjie84 2026-05-11 7/350 2026-05-15 23:20 by 同仁堂教主
[基金申请] 青C资助名额大幅增加! +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见