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

» 猜你喜欢

LWT投已经有3人回复
真的好想读博！已经有22人回复
博士白读了已经有49人回复
基金开始函评了吗？已经有16人回复
没学上已经有12人回复
研二光催化6月底四篇二区什么水平已经有15人回复
如何分离原料与产品已经有12人回复
环肽的合成已经有5人回复
学硕专硕已经有5人回复
0854－0855调剂已经有12人回复

1楼2020-11-11 16:09:02

luyaobao

木虫 (著名写手)

应助: 1 (幼儿园)
金币: 2443.5
散金: 889
红花: 22
沙发: 8
帖子: 1524
在线: 149.9小时
虫号: 7847640
注册: 2018-01-30
性别: GG
专业: 凝聚态物性I:结构、力学和

用pair_style dpd/tstat

发自小木虫Android客户端

赞一下(1人)

2楼2020-11-11 18:47:35

TW九九

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 76.5
帖子: 14
在线: 11.9小时
虫号: 23899513
注册: 2020-10-16
性别: GG
专业: 传热传质学

好的，这几天来我这看到了这个命令，但还是一直不确定，因为之前控温都是用的fix命令，而且复现的文献中说采用的LJ势，是不是可以用 pair_style hybrid lj/cut …… dpd/tstat ……命令

3楼2020-11-15 10:51:34

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

TW九九

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 76.5
帖子: 14
在线: 11.9小时
虫号: 23899513
注册: 2020-10-16
性别: GG
专业: 传热传质学

引用回帖:

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

谢谢，那对应于pair_style dpd/tstat的pair_coeff 参数是什么呢，我看手册说分别是A (force units)、γ(force/velocity units)、cutoff (distance units)，其中的γ是阻尼系数吗，这个数应该怎么得到呢，我做的是氩原子的。

5楼2020-11-16 19:10:28

luyaobao

木虫 (著名写手)

应助: 1 (幼儿园)
金币: 2443.5
散金: 889
红花: 22
沙发: 8
帖子: 1524
在线: 149.9小时
虫号: 7847640
注册: 2018-01-30
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

5楼: Originally posted by TW九九 at 2020-11-16 19:10:28
谢谢，那对应于pair_style dpd/tstat的pair_coeff 参数是什么呢，我看手册说分别是A (force units)、γ(force/velocity units)、cutoff (distance units)，其中的γ是阻尼系数吗，这个数应该怎么得到呢，我做的是氩 ...

温控没有系数A，γ不同液体粘度不同，看你的体系是什么了

发自小木虫Android客户端

6楼2020-11-16 19:37:35

TW九九

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 76.5
帖子: 14
在线: 11.9小时
虫号: 23899513
注册: 2020-10-16
性别: GG
专业: 传热传质学

引用回帖:

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

TW九九

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 76.5
帖子: 14
在线: 11.9小时
虫号: 23899513
注册: 2020-10-16
性别: GG
专业: 传热传质学

引用回帖:

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
...

写错了，我的in文件没有fix nvt那一条命令，然后原子始终不动，后来我试着加上nvt，又动的太快了，几千步就原子丢失了

8楼2020-11-22 10:40:33