 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (3151)
- >基金申請(qǐng) (216)
- >導(dǎo)師招生 (209)
- >考研 (175)
- >文獻(xiàn)求助 (113)
- >蟲(chóng)友互識(shí) (106)
- >考博 (97)
- >休閑灌水 (44)
- >論文道賀祈福 (25)
- >論文投稿 (24)
- >招聘信息布告欄 (22)
- >博后之家 (15)
- >碩博家園 (13)
- >找工作 (13)
- >SciFinder/Reaxys (9)
- >教師之家 (7)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4042 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
Mayerss新蟲(chóng) (小有名氣)
|
[求助]
關(guān)于Green-Kubo方法模擬水的粘度 已有3人參與
|
||
|
求問(wèn)有沒(méi)有大神用lammps做過(guò)水的模擬？ 我現(xiàn)在改編了一下manual里面計(jì)算液氬粘度的那個(gè)算例，來(lái)模擬水的粘度，300k下得到的結(jié)果是0.2mPa.s，查得的標(biāo)準(zhǔn)值是0.8mPa.s，請(qǐng)路過(guò)的大神幫忙看一下是哪里出了問(wèn)題，多謝！ in文件貼在下面，data文件和log文件在附件中。 # LAMMPS input script for viscosity of purewater units real variable T equal 300 variable V equal vol variable dt equal 0.4 variable p equal 400 # correlation length variable s equal 5 # sample interval variable d equal $p*$s # dump interval # convert from LAMMPS real units to SI variable kB equal 1.3806504e-23 # [J/K/ Boltzmann variable atm2Pa equal 101325.0 variable A2m equal 1.0e-10 variable fs2s equal 1.0e-15 variable convert equal ${atm2Pa}*${atm2Pa}*${fs2s}*${A2m}*${A2m}*${A2m} # setup problem dimension 3 boundary p p p atom_style full bond_style harmonic angle_style harmonic pair_style lj/cut/tip4p/long 1 2 1 1 0.125 12 15 kspace_style pppm/tip4p 1e-6 read_data sketch.data pair_coeff 1 2 0.0 0.0 # O-H pair_coeff 2 2 0.0 0.0 # H-H pair_coeff 1 1 0.16275 3.16435 # O-O timestep ${dt} thermo $d # equilibration and thermalization velocity all create $T 102486 mom yes rot yes dist gaussian fix NVT all nvt temp $T $T 40 drag 0.2 run 100000 # viscosity calculation, switch to NVE if desired #unfix NVT #fix NVE all nve reset_timestep 0 variable pxy equal pxy variable pxz equal pxz variable pyz equal pyz fix SS all ave/correlate $s $p $d & v_pxy v_pxz v_pyz type auto file S0St.dat ave running variable scale equal ${convert}/(${kB}*$T)*$V*$s*${dt} variable v11 equal trap(f_SS[3])*${scale} variable v22 equal trap(f_SS[4])*${scale} variable v33 equal trap(f_SS[5])*${scale} thermo_style custom step temp press v_pxy v_pxz v_pyz v_v11 v_v22 v_v33 run 200000 variable v equal (v_v11+v_v22+v_v33)/3.0 variable ndens equal count(all)/vol print "average viscosity: $v [Pa.s/ @ $T K, ${ndens} /A^3" |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : log.lammps
- 附件 2 : sketch.data
2015-05-23 19:35:01, 21.42 K
2015-05-23 19:35:02, 361.88 K
» 猜你喜歡
285化工學(xué)碩求調(diào)劑（081700）
已經(jīng)有10人回復(fù)
-
求材料調(diào)劑
已經(jīng)有10人回復(fù)
-
本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑
已經(jīng)有4人回復(fù)
-
國(guó)自科面上基金字體
已經(jīng)有8人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有11人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 25°水蒸氣粘度
已經(jīng)有4人回復(fù)
- 液體粘度的計(jì)算
已經(jīng)有14人回復(fù)
- Green一Kubo公式的應(yīng)用
已經(jīng)有5人回復(fù)
- 1stopt能進(jìn)行隱函數(shù)非線性擬合嗎？誰(shuí)能幫擬合一下新型液體的粘度混合計(jì)算公式
已經(jīng)有5人回復(fù)
- 各位大俠請(qǐng)問(wèn)用fluent仿真液固兩相流時(shí)固相的粘度需要設(shè)置嗎?
已經(jīng)有12人回復(fù)
- 湍流粘度的意義
已經(jīng)有10人回復(fù)
- 計(jì)算精度不高，且出現(xiàn)湍流粘度限制怎么辦
已經(jīng)有8人回復(fù)
- 海水的動(dòng)力粘度u是多少？哪位大俠知道，，，十分感謝
已經(jīng)有10人回復(fù)
- 請(qǐng)問(wèn)：剪切粘度（shear viscosity）和容積粘度（bulk viscosity）的區(qū)別
已經(jīng)有5人回復(fù)
- 【資源】關(guān)于圓管內(nèi)氣體流場(chǎng)中雷諾數(shù)與氣流速度及動(dòng)力粘度的問(wèn)題
已經(jīng)有9人回復(fù)
- 【討論】Green-Kubo公式是基于第一性原理的嗎？ ？
已經(jīng)有4人回復(fù)
brucefan
專家顧問(wèn) (著名寫(xiě)手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲(chóng)號(hào): 3048965
- 注冊(cè): 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
3樓2016-10-05 01:10:05
kakaxicheng
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 26.3
- 帖子: 20
- 在線: 45.8小時(shí)
- 蟲(chóng)號(hào): 3523181
- 注冊(cè): 2014-11-06
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2016-10-04 15:58:53
kakaxicheng
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 26.3
- 帖子: 20
- 在線: 45.8小時(shí)
- 蟲(chóng)號(hào): 3523181
- 注冊(cè): 2014-11-06
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
4樓2016-10-25 09:56:35
xu_jiangxin
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 8
- 帖子: 2
- 在線: 7小時(shí)
- 蟲(chóng)號(hào): 5583966
- 注冊(cè): 2017-02-15
- 專業(yè): 熱物性與熱物理測(cè)試技術(shù)
【答案】應(yīng)助回帖
|
一樓的同學(xué)說(shuō)不用轉(zhuǎn)化為SI ，這個(gè)我保留意見(jiàn)。我在自己的模擬中使用了單位轉(zhuǎn)化。關(guān)于你的data文件和in文件我說(shuō)一下我的看法：你使用的是tip4p模型，卻使用了很多的spce模型的參數(shù)而不是tip4p模型的參數(shù)。 這是我自己的模擬結(jié)果：（300K,1.29E+5mol/m3） spce tip4p tip4p2005 0.567 0.420 0.721 使用的參數(shù)為lammps手冊(cè)“6.8 TIP4P water mode”中給出的參數(shù)，l |
5樓2017-08-28 09:13:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|