| 2 | 1/1 | 返回列表 |
| 查看: 1448 | 回復(fù): 1 | ||
| 【懸賞金幣】回答本帖問題,作者cuteCarl將贈(zèng)送您 5 個(gè)金幣 | ||
cuteCarl新蟲 (初入文壇)
|
[求助]
lammps 周期邊界條件為啥會(huì)原子丟失啊
|
|
|
新人,金幣不多,請(qǐng)見諒 我什么我使用了周期邊界條件,還會(huì)出現(xiàn)院子丟失的現(xiàn)象啊 下面是我的in文件 echo screen dimension 3 boundary p p fm units metal atom_style atomic neighbor 2.0 bin neigh_modify every 1 delay 5 check yes region box block 0.0 100.0 0.0 25.0 0.0 200.0 create_box 2 box region wall_down block 0.0 100.0 0.0 25.0 0.0 10.0 region wall_up block 0.0 100.0 0.0 25.0 190.0 200.0 region wall_shape1 cone z 12.5 12.5 10.0 0.0 10.0 20.0 region wall_shape2 cone z 37.5 12.5 10.0 0.0 10.0 20.0 region wall_shape3 cone z 62.5 12.5 10.0 0.0 10.0 20.0 region wall_shape4 cone z 87.5 12.5 10.0 0.0 10.0 20.0 region flow block 0.0 100.0 0.0 25.0 10.0 30.0 lattice fcc 5.738 create_atoms 1 region flow lattice fcc 3.483 create_atoms 1 region wall_down create_atoms 1 region wall_up create_atoms 1 region wall_shape1 create_atoms 1 region wall_shape2 create_atoms 1 region wall_shape3 create_atoms 1 region wall_shape4 mass 1 39.948 mass 2 39.948 pair_style hybrid lj/cut 8.5 dpd/tstat 10.0 98.0 8.5 34387 pair_coeff 1 1 lj/cut 0.0103125 3.4 pair_coeff 1 2 lj/cut 0.0505 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_shape1 group 2 region wall_shape2 group 3 region wall_shape3 group 4 region wall_shape4 group 5 region wall_down group 5 region wall_up group wall union 1 2 3 4 5 6 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 10.0 98.0 0.02 compute 1 Ar chunk/atom bin/2d y lower 2 z 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 new.profile timestep 0.005 thermo 20 thermo_style custom step temp pe press density etotal dump 1 all custom 20 new.dump id x y z run 310000 各位大佬,幫忙看一下吧 |
新蟲 (初入文壇)
|
再請(qǐng)教一下,我加了fix 4 flow wall/reflect zhi 200.0還是不行,是不是我寫錯(cuò)了,我看手冊(cè)是這樣呀 新手太難了發(fā)自小木蟲Android客戶端 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩290求調(diào)劑 +5 | 杰尼龜aaa 2026-03-04 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +4 | Bahati 2026-03-05 | 4/200 |
|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +3 | 鋼鐵大炮 2026-03-04 | 3/150 |
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 3/150 |
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +9 | 土土小蟲 2026-03-03 | 11/550 |
|
|
[考研] 304求調(diào)劑 +7 | 曼殊2266 2026-02-28 | 8/400 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
[考研] 267化工調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
[考研] 成績276,專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
[考研] 0855機(jī)械工程324求調(diào)劑 +3 | 鄭SK 2026-03-02 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩0927 2026-03-02 | 5/250 |
|
|
[考研] 290求調(diào)劑 +9 | 材料專碩調(diào)劑; 2026-02-28 | 11/550 |
|