 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3584 | 回復(fù): 14 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
清風(fēng)染翠金蟲 (小有名氣)
|
[求助]
液體粘度的計(jì)算 已有2人參與
|
|||
|
最近在用MD模擬計(jì)算一種液體的粘度，使用fortran編的程序，模擬區(qū)域是一個(gè)長方體，總共有4096個(gè)粒子，三個(gè)方向都是周期邊界，通過熱浴保持體系溫度恒定。用Green-Kubo公式計(jì)算粘度。模擬的參數(shù)都按照文獻(xiàn)中給定的，所以我是在復(fù)現(xiàn)文獻(xiàn)中的結(jié)果。但是遇到一個(gè)問題。 粘度的積分步數(shù)在一萬步時(shí)，得到的結(jié)果與文獻(xiàn)中符合很好，當(dāng)積分步數(shù)小于十萬步時(shí)，結(jié)果看起來還可以，但是積分步數(shù)在100萬步時(shí)，結(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分子動(dòng)力學(xué) | 動(dòng)力學(xué)或熱力學(xué)計(jì)算 |
» 猜你喜歡
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有9人回復(fù)
-
085601一志愿西北工業(yè)大學(xué)初試346
已經(jīng)有4人回復(fù)
-
材料與化工（0856）304求B區(qū)調(diào)劑
已經(jīng)有10人回復(fù)
-
本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
070300化學(xué)279求調(diào)劑
已經(jīng)有7人回復(fù)
-
085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑
已經(jīng)有5人回復(fù)
-
293分求調(diào)劑，外語為俄語
已經(jīng)有4人回復(fù)
-
0856調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 可以讓液體粘度變化的材料
已經(jīng)有5人回復(fù)
- 求助：軟膠囊膠液水分和粘度的操作方法
已經(jīng)有7人回復(fù)
- 1stopt能進(jìn)行隱函數(shù)非線性擬合嗎？誰能幫擬合一下新型液體的粘度混合計(jì)算公式
已經(jīng)有5人回復(fù)
- 復(fù)數(shù)黏度和表觀粘度一樣嗎？
已經(jīng)有3人回復(fù)
- 石油運(yùn)動(dòng)粘度的測量
已經(jīng)有7人回復(fù)
- [求助]液體或者水的揮發(fā)速度的推導(dǎo)計(jì)算公式
已經(jīng)有7人回復(fù)
- 測低粘度溶液的粘度用哪種粘度計(jì)比較好，或者說那些
已經(jīng)有6人回復(fù)
- 固體溶解在液體中溶液的密度和粘度
已經(jīng)有10人回復(fù)
- 搖床上的培養(yǎng)瓶中的液體能計(jì)算剪切力嗎？粗略一點(diǎn)的就行？
已經(jīng)有5人回復(fù)
- 剪切粘度和拉伸粘度有何區(qū)別？
已經(jīng)有7人回復(fù)
- 支個(gè)招！高粘度高溫膠狀液體 如何測量粘度呢？水分？
已經(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ù)
清風(fēng)染翠
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1030
- 帖子: 51
- 在線: 20.1小時(shí)
- 蟲號(hào): 2048991
- 注冊: 2012-10-08
- 性別: GG
- 專業(yè): 多相流熱物理學(xué)
3樓2013-12-17 11:17:11
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲號(hào): 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(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
|
時(shí)間長的話噪音會(huì)有些大，可能只需積分到某個(gè)時(shí)間就可以了。 這篇文章里談到了這個(gè)問題：http://pubs.acs.org/doi/abs/10.1021/jp211952y 如果想進(jìn)一步研究的話這篇文章的參考文獻(xiàn)可能講的更詳細(xì)一些。 |
2樓2013-12-17 10:03:56
zlf9452006
金蟲 (小有名氣)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1365.7
- 散金: 539
- 紅花: 3
- 帖子: 59
- 在線: 111.1小時(shí)
- 蟲號(hào): 587860
- 注冊: 2008-08-22
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
清風(fēng)染翠: 金幣+30, ★★★很有幫助 2013-12-17 17:02:07
月只藍(lán): 金幣+1, 鼓勵(lì)交流！ 2013-12-18 08:57:07
感謝參與，應(yīng)助指數(shù) +1
清風(fēng)染翠: 金幣+30, ★★★很有幫助 2013-12-17 17:02:07
月只藍(lán): 金幣+1, 鼓勵(lì)交流！ 2013-12-18 08:57:07
|
個(gè)人認(rèn)為是由周期性設(shè)置引起的。GK方法本來就是由宏觀過程推導(dǎo)出來的，加之其本身信噪比極低，其漲落大是必然的。模擬中的周期性在宏觀上導(dǎo)致了液體性質(zhì)在個(gè)別方面的顯著變化。 建議樓主考慮下非平衡方法，要穩(wěn)定很多，也要準(zhǔn)確得多。 |
4樓2013-12-17 12:38:59
清風(fēng)染翠
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1030
- 帖子: 51
- 在線: 20.1小時(shí)
- 蟲號(hào): 2048991
- 注冊: 2012-10-08
- 性別: GG
- 專業(yè): 多相流熱物理學(xué)
5樓2013-12-17 17:09:09
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293分求調(diào)劑，外語為俄語 +4 | 加一一九 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085601 329分調(diào)劑 +3 | yzsa12 2026-03-31 | 3/150 |
|
|
|
[考研] 085600 295分求調(diào)劑 +7 | W55j 2026-03-30 | 9/450 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
|
[考研] 309求調(diào)劑 +15 | 誰不是少年 2026-03-29 | 15/750 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級(jí)已過，有科研經(jīng)歷，求調(diào)劑！ +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +14 | 北道巷 2026-03-24 | 14/700 |
|
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
|
[考研] 291求調(diào)劑 +8 | HanBeiNingZC 2026-03-24 | 8/400 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +10 | 1124361 2026-03-24 | 10/500 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 085600，材料與化工321分，求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|