| 查看: 2476 | 回复: 7 | |||
| 【悬赏金币】回答本帖问题，作者TW九九将赠送您 5 个金币 | |||
[求助]
lammps如何实现耗散粒子动力学控温
|
|||
想问下lammps如何采用耗散粒子动力学控温方式，nvt、npt这种好像都用的Nose-Hoover控温，我看了几个关于dpd的命令感觉都不是。 微信图片_20201111160502.png |
回复此楼» 猜你喜欢
303求调剂
已经有6人回复
-
07化学303求调剂
已经有5人回复
-
332求调剂
已经有6人回复
-
考研调剂
已经有7人回复
-
0856求调剂
已经有5人回复
-
总分293求调剂
已经有6人回复
-
275求调剂
已经有5人回复
-
085600 材料与化工 329分求调剂
已经有7人回复
-
材料与化工328分调剂
已经有6人回复
-
291 求调剂
已经有6人回复
luyaobao
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2457.9
- 散金: 1091
- 红花: 23
- 沙发: 8
- 帖子: 1533
- 在线: 150.7小时
- 虫号: 7847640
- 注册: 2018-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2020-11-11 18:47:35
3楼2020-11-15 10:51:34
4楼2020-11-15 14:37:59
5楼2020-11-16 19:10:28
luyaobao
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2457.9
- 散金: 1091
- 红花: 23
- 沙发: 8
- 帖子: 1533
- 在线: 150.7小时
- 虫号: 7847640
- 注册: 2018-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
6楼2020-11-16 19:37:35
|
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 fix 1 all nvt temp 86 86 0.02 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 2 wall setforce 0.0 0.0 0.0 fix 3 flow ave/chunk 100 5 500 1 density/mass file c_0.5.profile timestep 0.005 thermo 2000 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文件输出的结果从始至终原子都不动，我文件中也设置了velocity命令了，您知道吗 |
7楼2020-11-20 20:34:05
8楼2020-11-22 10:40:33
