| 3 | 1/1 | 返回列表 |
| 查看: 1877 | 回復(fù): 2 | ||
放逐孤島銅蟲 (小有名氣)
|
[求助]
lammps超聲波模擬設(shè)置 已有1人參與
|
|
我的題目是使用lammps模擬超聲空化,我把體系分為兩部份,一部分為邊界區(qū),一部分為主區(qū)。通過邊界分子的簡(jiǎn)諧運(yùn)動(dòng)來實(shí)現(xiàn)超聲波的傳遞進(jìn)而完成空化過程。但是我在計(jì)算了大約20多個(gè)簡(jiǎn)諧運(yùn)動(dòng)周期后,主區(qū)內(nèi)的分子依然不動(dòng),僅僅是邊界區(qū)的分子在進(jìn)行簡(jiǎn)諧運(yùn)動(dòng)。而且出現(xiàn)了 ERROR: Out of range atoms - cannot compute PPPM (src/KSPACE/pppm_tip4p.cpp:106) 這個(gè)錯(cuò)誤(我是用的是TIP4P水分子)。希望大家能幫助我,這是哪里出了問題呢,應(yīng)該如何修改,也歡迎各位做超聲空化的同僚前來交流,謝謝大家 ![]() 我的in文件: dimension 3 boundary p p p units real atom_style full read_data b125.data group 1 read_data b2000.data add append group 2 neighbor 2 bin neigh_modify delay 0 every 1 pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5 pair_coeff 1 1 0.1521 3.154000 8.5 pair_coeff * 2 0.000000 0.000000 8.5 bond_style harmonic bond_coeff 1 1000 0.9572 angle_style harmonic angle_coeff 1 3000 104.52 kspace_style pppm/tip4p 0.00001 minimize 1.0e-12 1.0e-12 1000 1000 variable vel equal 0.1*sin(2*PI*0.01*time) fix 1 1 move variable NULL NULL NULL v_vel NULL NULL timestep 0.001 dump 1 all atom 100 dump.lammpstrj run 500000 文件中1即為邊界區(qū),2為主區(qū)。 |
銅蟲 (小有名氣)
新蟲 (初入文壇)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工,270求調(diào)劑 +12 | YXCT 2026-03-01 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 085601 材料305分求助 +3 | 泡泡郵件 2026-03-03 | 5/250 |
|
|
[考研] 347求調(diào)劑 +3 | 啊歐歐歐 2026-03-03 | 3/150 |
|
|
[基金申請(qǐng)] 請(qǐng)問大家,研究風(fēng)險(xiǎn)與應(yīng)對(duì)措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 材料284求調(diào)劑,一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
[考研] 302材料工程求調(diào)劑 +5 | Doleres 2026-03-01 | 6/300 |
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
[考博] 26申博 +4 | 想申博! 2026-02-26 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|
|
[考研] 調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-02-28 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-02-28 | 5/250 |
|