版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者TW九九將贈(zèng)送您 5 個(gè)金幣

TW九九

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76.5
帖子: 14
在線: 11.9小時(shí)
蟲號(hào): 23899513
注冊(cè): 2020-10-16
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

[求助] lammps如何實(shí)現(xiàn)耗散粒子動(dòng)力學(xué)控溫

想問下lammps如何采用耗散粒子動(dòng)力學(xué)控溫方式，nvt、npt這種好像都用的Nose-Hoover控溫，我看了幾個(gè)關(guān)于dpd的命令感覺都不是。

微信圖片_20201111160502.png

回復(fù)此樓

» 猜你喜歡

0856材料與化工，270求調(diào)劑已經(jīng)有15人回復(fù)
085600 材料與化工 298 已經(jīng)有12人回復(fù)
312求調(diào)劑已經(jīng)有7人回復(fù)
304分材料專碩求調(diào)劑已經(jīng)有9人回復(fù)
320材料與化工，求調(diào)劑已經(jīng)有5人回復(fù)
成果系統(tǒng)訪問量大，請(qǐng)15分鐘后再嘗試。由此給您造成的不便，敬請(qǐng)諒解。已經(jīng)有12人回復(fù)
281求調(diào)劑已經(jīng)有5人回復(fù)
347求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
一志愿314求調(diào)劑已經(jīng)有6人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

luyaobao

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 2457.9
散金: 1091
紅花: 23
沙發(fā): 8
帖子: 1533
在線: 150.7小時(shí)
蟲號(hào): 7847640
注冊(cè): 2018-01-30
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

用pair_style dpd/tstat

發(fā)自小木蟲Android客戶端

贊一下(1人)

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

TW九九

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76.5
帖子: 14
在線: 11.9小時(shí)
蟲號(hào): 23899513
注冊(cè): 2020-10-16
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

好的，這幾天來我這看到了這個(gè)命令，但還是一直不確定，因?yàn)橹翱販囟际怯玫膄ix命令，而且復(fù)現(xiàn)的文獻(xiàn)中說采用的LJ勢，是不是可以用 pair_style hybrid lj/cut …… dpd/tstat ……命令

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

luyaobao

木蟲 (著名寫手)

引用回帖:

3樓: Originally posted by TW九九 at 2020-11-15 10:51:34
好的，這幾天來我這看到了這個(gè)命令，但還是一直不確定，因?yàn)橹翱販囟际怯玫膄ix命令，而且復(fù)現(xiàn)的文獻(xiàn)中說采用的LJ勢，是不是可以用 pair_style hybrid lj/cut …… dpd/tstat ……命令

是

發(fā)自小木蟲Android客戶端

4樓2020-11-15 14:37:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

TW九九

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76.5
帖子: 14
在線: 11.9小時(shí)
蟲號(hào): 23899513
注冊(cè): 2020-10-16
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

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)該怎么得到呢，我做的是氬原子的。

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

luyaobao

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 2457.9
散金: 1091
紅花: 23
沙發(fā): 8
帖子: 1533
在線: 150.7小時(shí)
蟲號(hào): 7847640
注冊(cè): 2018-01-30
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

5樓: Originally posted by TW九九 at 2020-11-16 19:10:28
謝謝，那對(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)該怎么得到呢，我做的是氬 ...

溫控沒有系數(shù)A，γ不同液體粘度不同，看你的體系是什么了

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

TW九九

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76.5
帖子: 14
在線: 11.9小時(shí)
蟲號(hào): 23899513
注冊(cè): 2020-10-16
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

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命令了，您知道嗎

贊一下

回復(fù)此樓

7樓2020-11-20 20:34:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

TW九九

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 76.5
帖子: 14
在線: 11.9小時(shí)
蟲號(hào): 23899513
注冊(cè): 2020-10-16
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

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)的太快了，幾千步就原子丟失了

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 TW九九的主題更新

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎(jiǎng)勵(lì)，但不允許灌水，必須填寫15個(gè)字符以上)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 304分材料專碩求調(diào)劑 +7	qiuzhigril 2026-03-03	9/450	2026-03-04 14:25 by kakakapanpan
[考研] 322,求調(diào)劑 +3	菜菜愛玩 2026-03-04	3/150	2026-03-04 12:15 by xiongyaxuan
[考研] 0857求調(diào)劑 +3	CDDD1013 2026-03-04	3/150	2026-03-04 11:26 by lature00
[考研] 085601 材料工程 320 +3	和樂瑤 2026-03-03	3/150	2026-03-04 09:09 by gaoxiaoniuma
[考研] 材料工程269求調(diào)劑 +7	白刺玫 2026-03-02	7/350	2026-03-04 08:47 by zhyzzh
[考研] 291求調(diào)劑 +4	Afy123456 2026-03-03	7/350	2026-03-03 22:29 by barlinike
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5	Qwertyuop 2026-03-01	5/250	2026-03-03 17:50 by 梅子鴨鴨
[考研] 283求調(diào)劑 +9	鹿沫笙 2026-03-02	9/450	2026-03-03 15:58 by tgxtgxtgx9
[考研] 求調(diào)劑 +4	Guo_yuxuan 2026-03-02	5/250	2026-03-03 14:39 by xiaomc_gzh
[考研] 計(jì)算機(jī)學(xué)碩分?jǐn)?shù)285求調(diào)劑 +4	glwshine 2026-03-02	5/250	2026-03-03 14:27 by king呀
[考研] 284求調(diào)劑 +6	天下熯 2026-03-02	6/300	2026-03-03 13:46 by Iveryant
[考研] 材料工程專碩283求調(diào)劑 5+8	，��？ 2026-03-02	10/500	2026-03-03 12:01 by EBSD
[考研] 材料工程求調(diào)劑 +3	1431251 2026-03-03	3/150	2026-03-03 11:58 by EBSD
[考研] 295求調(diào)劑。一志愿報(bào)考鄭州大學(xué)化學(xué)工藝學(xué)碩，總分295分 +8	yl1 2026-03-02	9/450	2026-03-02 21:21 by sunny81
[考研] 一志愿中石油（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-02	3/150	2026-03-02 18:54 by caszguilin
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5	半截的詩0927 2026-03-02	5/250	2026-03-02 18:37 by 明亮9527
[考研] 284求調(diào)劑 +10	天下熯 2026-02-28	11/550	2026-03-02 11:03 by 無際的草原
[基金申請(qǐng)] 本子寫完了，給DS兄弟看了，得了92分 +3	Doma 2026-03-01	7/350	2026-03-02 00:00 by jnzsy
[基金申請(qǐng)] 剛錄用，沒有期刊號(hào)，但是在線可看的論文可以放為代表作嗎 10+3	arang1 2026-03-01	3/150	2026-03-01 16:43 by babero
[論文投稿] Optics letters投稿被拒求助 30+3	luckyry 2026-02-26	4/200	2026-03-01 09:06 by babero