 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3573 | 回復(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分子動力學(xué) | 動力學(xué)或熱力學(xué)計(jì)算 |
» 猜你喜歡
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有7人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有7人回復(fù)
-
289求調(diào)劑
已經(jīng)有13人回復(fù)
-
299求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有9人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
調(diào)劑求院校招收
已經(jīng)有5人回復(fù)
-
調(diào)劑310
已經(jīng)有8人回復(fù)
-
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有5人回復(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)動粘度的測量
已經(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í)
- 蟲號: 2048991
- 注冊: 2012-10-08
- 性別: GG
- 專業(yè): 多相流熱物理學(xué)
9樓2013-12-18 13:17:31
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲號: 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í)間長的話噪音會有些大，可能只需積分到某個(gè)時(shí)間就可以了。 這篇文章里談到了這個(gè)問題：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小時(shí)
- 蟲號: 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小時(shí)
- 蟲號: 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, 鼓勵交流！ 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
|
個(gè)人認(rèn)為是由周期性設(shè)置引起的。GK方法本來就是由宏觀過程推導(dǎo)出來的，加之其本身信噪比極低，其漲落大是必然的。模擬中的周期性在宏觀上導(dǎo)致了液體性質(zhì)在個(gè)別方面的顯著變化。 建議樓主考慮下非平衡方法，要穩(wěn)定很多，也要準(zhǔn)確得多。 |
4樓2013-12-17 12:38:59
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +8 | 張zz111 2026-03-27 | 9/450 |
|
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|