24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2255  |  回复: 7
【悬赏金币】回答本帖问题,作者TW九九将赠送您 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

TW九九

铜虫 (初入文坛)

[求助] lammps如何实现耗散粒子动力学控温

想问下lammps如何采用耗散粒子动力学控温方式,nvt、npt这种好像都用的Nose-Hoover控温,我看了几个关于dpd的命令感觉都不是。

lammps如何实现耗散粒子动力学控温
微信图片_20201111160502.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TW九九

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by luyaobao at 2020-11-16 12:37:35
温控没有系数A,γ不同液体粘度不同,看你的体系是什么了
...

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

luyaobao

木虫 (著名写手)

2楼2020-11-11 18:47:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TW九九

铜虫 (初入文坛)

好的,这几天来我这看到了这个命令,但还是一直不确定,因为之前控温都是用的fix命令,而且复现的文献中说采用的LJ势,是不是可以用 pair_style hybrid  lj/cut  ……  dpd/tstat ……命令
3楼2020-11-15 10:51:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luyaobao

木虫 (著名写手)

引用回帖:
3楼: Originally posted by TW九九 at 2020-11-15 10:51:34
好的,这几天来我这看到了这个命令,但还是一直不确定,因为之前控温都是用的fix命令,而且复现的文献中说采用的LJ势,是不是可以用 pair_style hybrid  lj/cut  ……  dpd/tstat ……命令



发自小木虫Android客户端
4楼2020-11-15 14:37:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见