| 6 | 1/1 | 返回列表 |
| 查看: 1641 | 回復(fù): 5 | |||
tsadaye金蟲 (初入文壇)
|
[求助]
fix viscosity命令求助
|
|
剛剛開始學(xué)習(xí)lammps,想先試著用dpd方法計(jì)算一個(gè)單一珠子體系的粘度,其中粘度由fix viscosity命令計(jì)算。寫的in文件如下: dimension 3 units lj atom_style atomic read_data data.initial neighbor 1.2 bin neigh_modify delay 3 pair_style dpd 1.0 1.0 1234567 pair_coeff 1 1 25.0 4.5 communicate single vel yes velocity all create 1.0 4928259 dist gaussian fix 1 all nve fix 2 all viscosity 50 x z 10 timestep 0.05 thermo 10 dump 1 all local 1000 initial.dump f_2 run 30000 write_restart initial.restart 問題是fix viscosity命令將計(jì)算一個(gè)粘度,但這個(gè)粘度應(yīng)該由什么命令讀取出來?我試著用dump local,但是貌似這個(gè)是個(gè)group attribute,讀不出來。請(qǐng)高人解惑,謝謝。 第一次寫in文件,難免很多錯(cuò)誤,請(qǐng)輕噴 |
動(dòng)力學(xué)或熱力學(xué)計(jì)算 |
新蟲 (小有名氣)
| 樓主manual沒有仔細(xì)看呀,這個(gè)命令還算不出粘度呢,只是得到體系的總的動(dòng)量交換值,除以時(shí)間和模擬盒子的橫截面積得到動(dòng)量通量,再除以什么velocity profile的斜率就得到粘度了,這個(gè)命令的計(jì)算值可以用print f_2命令在屏幕上顯示出來,然后在log文件里面也能找到了。我目前也剛弄出這些,對(duì)于velocity profile一籌莫展啊,希望能互相探討探討~~~ |

新蟲 (初入文壇)

新蟲 (初入文壇)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 招08考數(shù)學(xué) +4 | laoshidan 2026-03-20 | 8/400 |
|
|
[考研] 306求調(diào)劑 +6 | chuanzhu川燭 2026-03-18 | 6/300 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +3 | jiajunX 2026-03-22 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗(yàn)) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|