 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 4045 | 回復: 5 | |||
[求助]
關于Green-Kubo方法模擬水的粘度 已有3人參與
|
|
求問有沒有大神用lammps做過水的模擬？ 我現在改編了一下manual里面計算液氬粘度的那個算例，來模擬水的粘度，300k下得到的結果是0.2mPa.s，查得的標準值是0.8mPa.s，請路過的大神幫忙看一下是哪里出了問題，多謝！ 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" |
回復此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯系郵箱：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
» 猜你喜歡
324分 085600材料化工求調劑
已經有3人回復
-
317求調劑
已經有6人回復
-
0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章
已經有11人回復
-
307求調劑
已經有5人回復
-
求調劑
已經有3人回復
-
求調劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分
已經有3人回復
-
化學求調劑
已經有4人回復
-
招收調劑碩士
已經有9人回復
-
0703化學調劑
已經有8人回復
-
0703化學調劑
已經有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 25°水蒸氣粘度
已經有4人回復
- 液體粘度的計算
已經有14人回復
- Green一Kubo公式的應用
已經有5人回復
- 1stopt能進行隱函數非線性擬合嗎？誰能幫擬合一下新型液體的粘度混合計算公式
已經有5人回復
- 各位大俠請問用fluent仿真液固兩相流時固相的粘度需要設置嗎?
已經有12人回復
- 湍流粘度的意義
已經有10人回復
- 計算精度不高，且出現湍流粘度限制怎么辦
已經有8人回復
- 海水的動力粘度u是多少？哪位大俠知道，，，十分感謝
已經有10人回復
- 請問：剪切粘度（shear viscosity）和容積粘度（bulk viscosity）的區(qū)別
已經有5人回復
- 【資源】關于圓管內氣體流場中雷諾數與氣流速度及動力粘度的問題
已經有9人回復
- 【討論】Green-Kubo公式是基于第一性原理的嗎？ ？
已經有4人回復
kakaxicheng
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 26.3
- 帖子: 20
- 在線: 45.8小時
- 蟲號: 3523181
- 注冊: 2014-11-06
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
2樓2016-10-04 15:58:53
brucefan
專家顧問 (著名寫手)
-
專家經驗: +218 - 模擬EPI: 4
- 應助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 分子模擬
3樓2016-10-05 01:10:05
kakaxicheng
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 26.3
- 帖子: 20
- 在線: 45.8小時
- 蟲號: 3523181
- 注冊: 2014-11-06
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
4樓2016-10-25 09:56:35
5樓2017-08-28 09:13:22
SDXM
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 120.5
- 帖子: 5
- 在線: 34.9小時
- 蟲號: 2453408
- 注冊: 2013-05-07
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動化
|
你好，方便把您的data文件和input文件給我傳一份嗎？我也在做這個模擬，數據都來在6.8 TIP4P water mode，但是最終的結果卻對不上，謝謝您啦 我的郵箱是 1178408727@qq.com 非常感謝 |
6樓2018-03-14 04:04:02
相關版塊跳轉
我要訂閱樓主 Mayerss 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調劑 +5 | 冷笙123 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料考研調劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 材料080500調劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
|
[考研] 一志愿吉林大學材料學碩321求調劑 +6 | Ymlll 2026-03-18 | 9/450 |
|
|
|
[考研] 330求調劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 材料專業(yè)求調劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 295求調劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[基金申請] 被我言中：新模板不強調格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 材料工程專碩274一志愿211求調劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 考研調劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研]
 11408 一志愿西電，277分求調劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 085600材料與化工 求調劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 288求調劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[基金申請] 現在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 085600材料與化工 309分請求調劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
|
[考研] 070303一志愿西北大學學碩310找調劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|