| 6 | 1/1 | 返回列表 |
| 查看: 1645 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 081700 調(diào)劑 267分 +3 | 迷人的哈哈 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|