用pair_style dpd/tstat

好的，這幾天來我這看到了這個(gè)命令，但還是一直不確定，因?yàn)橹�前控溫都是用的fix命令，而且復(fù)現(xiàn)的文獻(xiàn)中說采用的LJ勢(shì)，是不是可以用 pair_style hybrid  lj/cut  ……  dpd/tstat ……命令

引用回帖:

4樓: Originally posted by luyaobao at 2020-11-15 07:37:59
是
...

謝謝，那對(duì)應(yīng)于pair_style dpd/tstat的pair_coeff 參數(shù)是什么呢，我看手冊(cè)說分別是A (force units)、γ(force/velocity units)、cutoff (distance units)，其中的γ是阻尼系數(shù)嗎，這個(gè)數(shù)應(yīng)該怎么得到呢，我做的是氬原子的。

引用回帖:

6樓: Originally posted by luyaobao at 2020-11-16 12:37:35
溫控沒有系數(shù)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


想問您以下為什么我這個(gè)in文件輸出的結(jié)果從始至終原子都不動(dòng)，我文件中也設(shè)置了velocity命令了，您知道嗎
，

引用回帖:

7樓: Originally posted by TW九九 at 2020-11-20 13:34:05
units                  metal
dimension                     3
boundary                     f f p

atom_style                   atomic
neighbor                     2.0 bin
neigh_modify                 every 1  delay 5 check yes
...

寫錯(cuò)了，我的in文件沒有fix nvt那一條命令，然后原子始終不動(dòng)，后來我試著加上nvt，又動(dòng)的太快了，幾千步就原子丟失了