| 查看: 2152 | 回复: 4 | |||
[交流]
为什么我的模拟结果原子不动 已有3人参与
|
|
units metal dimension 3 boundary f f p atom_style atomic neighbor 2.0 bin neigh_modify every 1 delay 5 check yes region box block 0 169 0 168 0 400 create_box 2 box region wall_left block 0 4 0 168 0 400 region wall_right block 164 169 0 168 0 400 region wall_down block 4 164 0 4 0 400 region flow block 4 164 4 44 0 400 lattice fcc 5.738 create_atoms 1 region flow lattice fcc 3.483 create_atoms 1 region wall_left create_atoms 1 region wall_right create_atoms 1 region wall_down mass 1 39.948 mass 2 39.948 pair_style hybrid lj/cut 8.5 dpd/tstat 86 86 8.5 34387 pair_coeff 1 1 lj/cut 0.0103125 3.4 pair_coeff 1 2 lj/cut 0.00515625 3.4 pair_coeff 2 2 lj/cut 0.00515625 3.4 pair_coeff 1 1 dpd/tstat 0.0672 8.5 comm_modify vel yes group 1 region wall_left group 2 region wall_right group 3 region wall_down group wall union 1 2 3 group flow region flow group Ar subtract all wall set group wall type 2 velocity flow create 86 482748 mom yes rot yes compute 1 Ar chunk/atom bin/2d x lower 2 y lower 2 units box velocity wall set 0.0 0.0 0.0 fix 1 wall setforce 0.0 0.0 0.0 fix 5 flow ave/chunk 100 5 500 1 density/mass file c_0.5.profile timestep 0.005 thermo 1000 thermo_style custom step temp pe press density etotal dump 1 all custom 1000 c_0.5.dump id type x y z run 310000 想问下为什么我这个in文件跑出来的结果原子完全不动，我明明也对flow原子施加了速度， 后来我又试着加上fix nvt命令，结果原子动的太快也不是想要的结果，有人知道吗。 |
回复此楼» 猜你喜欢
欢迎211本科同学，过A区国家线，A区非偏远一本，交叉学科课题组
已经有12人回复
-
成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。
已经有14人回复
-
281电子信息求调剂
已经有8人回复
-
材料277分求调剂
已经有7人回复
-
求调剂院校
已经有12人回复
-
293求调剂
已经有3人回复
-
化学 0703求调剂 总分293 一志愿211
已经有10人回复
-
320材料一志愿西工大专硕英二数二 有两年光伏工作经验
已经有7人回复
-
材料328求调剂
已经有8人回复
-
304分材料专硕求调剂
已经有12人回复
luyaobao
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2457.9
- 散金: 1091
- 红花: 23
- 沙发: 8
- 帖子: 1533
- 在线: 150.7小时
- 虫号: 7847640
- 注册: 2018-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2020-11-23 18:35:05
jluchemlib
铁杆木虫 (正式写手)
- 应助: 119 (高中生)
- 金币: 11607.4
- 散金: 200
- 红花: 19
- 帖子: 521
- 在线: 356.1小时
- 虫号: 589656
- 注册: 2008-08-30
- 性别: GG
- 专业: 高分子物理与高分子物理化
3楼2020-11-24 10:03:29
4楼2020-11-24 15:04:56
5楼2020-12-24 09:14:05
