| 8 | 1/1 | 返回列表 |
| 查看: 2412 | 回復: 7 | ||
| 【懸賞金幣】回答本帖問題,作者TW九九將贈送您 5 個金幣 | ||
TW九九銅蟲 (初入文壇)
|
[求助]
lammps如何實現(xiàn)耗散粒子動力學控溫
|
|
|
想問下lammps如何采用耗散粒子動力學控溫方式,nvt、npt這種好像都用的Nose-Hoover控溫,我看了幾個關(guān)于dpd的命令感覺都不是。 微信圖片_20201111160502.png |
木蟲 (著名寫手)
銅蟲 (初入文壇)
銅蟲 (初入文壇)
木蟲 (著名寫手)
銅蟲 (初入文壇)
|
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文件輸出的結(jié)果從始至終原子都不動,我文件中也設(shè)置了velocity命令了,您知道嗎 |
銅蟲 (初入文壇)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 博斯特525 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 298 +11 | 小西笑嘻嘻 2026-03-03 | 11/550 |
|
|
[考研] 0703化學 一志愿東北大學 279分求調(diào)劑 +3 | 聶樂平 2026-03-03 | 3/150 |
|
|
[考研] 080500材料科學與工程 +5 | 202114020319 2026-03-03 | 5/250 |
|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
[考研] 085601 材料305分求助 +3 | 泡泡郵件 2026-03-03 | 5/250 |
|
|
[考研] 接收調(diào)劑 +15 | 津萌津萌 2026-03-02 | 23/1150 |
|
|
[考研] 成績276,專業(yè)代碼0856求調(diào)劑 +7 | 小陳朵 2026-03-03 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 282求調(diào)劑 +4 | 2103240126 2026-03-02 | 7/350 |
|
|
[考研] 材料化工調(diào)劑 +13 | 今夏不夏 2026-03-01 | 16/800 |
|
|
[考研] 0856材料求調(diào)劑 +12 | hyf hyf hyf 2026-02-28 | 13/650 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 化工專碩342,一志愿大連理工大學,求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-02-28 | 7/350 |
|
|
[考研] 調(diào)劑 +3 | 簡木ChuFront 2026-02-28 | 3/150 |
|