 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學(xué)術(shù)科研互動平臺 » 計算模擬區(qū) » 分子模擬 » 專家會診 » 液體粘度的計算
| 查看: 3545 | 回復(fù): 14 | ||||
清風(fēng)染翠金蟲 (小有名氣)
|
[求助]
液體粘度的計算 已有2人參與
|
|||
|
最近在用MD模擬計算一種液體的粘度，使用fortran編的程序，模擬區(qū)域是一個長方體，總共有4096個粒子，三個方向都是周期邊界，通過熱浴保持體系溫度恒定。用Green-Kubo公式計算粘度。模擬的參數(shù)都按照文獻(xiàn)中給定的，所以我是在復(fù)現(xiàn)文獻(xiàn)中的結(jié)果。但是遇到一個問題。 粘度的積分步數(shù)在一萬步時，得到的結(jié)果與文獻(xiàn)中符合很好，當(dāng)積分步數(shù)小于十萬步時，結(jié)果看起來還可以，但是積分步數(shù)在100萬步時，結(jié)果看起來是要發(fā)散的樣子，而且不同次模擬得到的100萬步圖像是不宜氧的，為什么積分步數(shù)一長就無法保持結(jié)果的收斂和穩(wěn)定性，有沒有人做過或者有什么看法的可以交流一下。 1萬.png 10萬.png 100萬.png [ Last edited by 清風(fēng)染翠 on 2013-12-16 at 14:16 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 MD分子動力學(xué) | 動力學(xué)或熱力學(xué)計算 |
» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 可以讓液體粘度變化的材料
已經(jīng)有5人回復(fù)
- 求助：軟膠囊膠液水分和粘度的操作方法
已經(jīng)有7人回復(fù)
- 1stopt能進(jìn)行隱函數(shù)非線性擬合嗎？誰能幫擬合一下新型液體的粘度混合計算公式
已經(jīng)有5人回復(fù)
- 復(fù)數(shù)黏度和表觀粘度一樣嗎？
已經(jīng)有3人回復(fù)
- 石油運(yùn)動粘度的測量
已經(jīng)有7人回復(fù)
- [求助]液體或者水的揮發(fā)速度的推導(dǎo)計算公式
已經(jīng)有7人回復(fù)
- 測低粘度溶液的粘度用哪種粘度計比較好，或者說那些
已經(jīng)有6人回復(fù)
- 固體溶解在液體中溶液的密度和粘度
已經(jīng)有10人回復(fù)
- 搖床上的培養(yǎng)瓶中的液體能計算剪切力嗎？粗略一點(diǎn)的就行？
已經(jīng)有5人回復(fù)
- 剪切粘度和拉伸粘度有何區(qū)別？
已經(jīng)有7人回復(fù)
- 支個招！高粘度高溫膠狀液體 如何測量粘度呢？水分？
已經(jīng)有7人回復(fù)
- 高爐熔渣粘度問題
已經(jīng)有6人回復(fù)
- 如何降低流體的粘度？
已經(jīng)有14人回復(fù)
- 粘度測量問題(標(biāo)準(zhǔn)物的選擇)
已經(jīng)有9人回復(fù)
- 磷酸粘度的測定
已經(jīng)有5人回復(fù)
- 【求助】粘度的定義
已經(jīng)有5人回復(fù)
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗: +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時
- 蟲號: 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
清風(fēng)染翠: 金幣+50, ★★★很有幫助 2013-12-17 17:01:58
月只藍(lán): 金幣+2, 感謝指導(dǎo)！ 2013-12-18 08:56:53
感謝參與，應(yīng)助指數(shù) +1
清風(fēng)染翠: 金幣+50, ★★★很有幫助 2013-12-17 17:01:58
月只藍(lán): 金幣+2, 感謝指導(dǎo)！ 2013-12-18 08:56:53
|
時間長的話噪音會有些大，可能只需積分到某個時間就可以了。 這篇文章里談到了這個問題：http://pubs.acs.org/doi/abs/10.1021/jp211952y 如果想進(jìn)一步研究的話這篇文章的參考文獻(xiàn)可能講的更詳細(xì)一些。 |
2樓2013-12-17 10:03:56
清風(fēng)染翠
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1030
- 帖子: 51
- 在線: 20.1小時
- 蟲號: 2048991
- 注冊: 2012-10-08
- 性別: GG
- 專業(yè): 多相流熱物理學(xué)
6樓2013-12-17 17:09:24
清風(fēng)染翠
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1030
- 帖子: 51
- 在線: 20.1小時
- 蟲號: 2048991
- 注冊: 2012-10-08
- 性別: GG
- 專業(yè): 多相流熱物理學(xué)
3樓2013-12-17 11:17:11
zlf9452006
金蟲 (小有名氣)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1365.7
- 散金: 539
- 紅花: 3
- 帖子: 59
- 在線: 111.1小時
- 蟲號: 587860
- 注冊: 2008-08-22
- 專業(yè): 理論和計算化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
清風(fēng)染翠: 金幣+30, ★★★很有幫助 2013-12-17 17:02:07
月只藍(lán): 金幣+1, 鼓勵交流！ 2013-12-18 08:57:07
感謝參與，應(yīng)助指數(shù) +1
清風(fēng)染翠: 金幣+30, ★★★很有幫助 2013-12-17 17:02:07
月只藍(lán): 金幣+1, 鼓勵交流！ 2013-12-18 08:57:07
|
個人認(rèn)為是由周期性設(shè)置引起的。GK方法本來就是由宏觀過程推導(dǎo)出來的，加之其本身信噪比極低，其漲落大是必然的。模擬中的周期性在宏觀上導(dǎo)致了液體性質(zhì)在個別方面的顯著變化。 建議樓主考慮下非平衡方法，要穩(wěn)定很多，也要準(zhǔn)確得多。 |
4樓2013-12-17 12:38:59
清風(fēng)染翠
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1030
- 帖子: 51
- 在線: 20.1小時
- 蟲號: 2048991
- 注冊: 2012-10-08
- 性別: GG
- 專業(yè): 多相流熱物理學(xué)
5樓2013-12-17 17:09:09
zlf9452006
金蟲 (小有名氣)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1365.7
- 散金: 539
- 紅花: 3
- 帖子: 59
- 在線: 111.1小時
- 蟲號: 587860
- 注冊: 2008-08-22
- 專業(yè): 理論和計算化學(xué)
7樓2013-12-18 04:36:33
sciencejoy
新蟲 (著名寫手)
- 應(yīng)助: 436 (碩士)
- 金幣: 11207.1
- 紅花: 90
- 帖子: 2974
- 在線: 498小時
- 蟲號: 802149
- 注冊: 2009-07-02
- 性別: GG
- 專業(yè): 高分子物理與高分子物理化
8樓2013-12-18 11:15:09
清風(fēng)染翠
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1030
- 帖子: 51
- 在線: 20.1小時
- 蟲號: 2048991
- 注冊: 2012-10-08
- 性別: GG
- 專業(yè): 多相流熱物理學(xué)
9樓2013-12-18 13:17:31
qdykswang
金蟲 (著名寫手)
- 應(yīng)助: 14 (小學(xué)生)
- 金幣: 6401.1
- 散金: 43
- 紅花: 5
- 帖子: 1195
- 在線: 309.8小時
- 蟲號: 1382019
- 注冊: 2011-08-28
- 性別: GG
- 專業(yè): 數(shù)值預(yù)報與數(shù)值模擬
10樓2013-12-19 07:23:02
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 清風(fēng)染翠 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工（085600）總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|