24小时热门版块排行榜    

查看: 3558  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

清风染翠

金虫 (小有名气)

[求助] 液体粘度的计算 已有2人参与

最近在用MD模拟计算一种液体的粘度,使用fortran编的程序,模拟区域是一个长方体,总共有4096个粒子,三个方向都是周期边界,通过热浴保持体系温度恒定。用Green-Kubo公式计算粘度。模拟的参数都按照文献中给定的,所以我是在复现文献中的结果。但是遇到一个问题。

粘度的积分步数在一万步时,得到的结果与文献中符合很好,当积分步数小于十万步时,结果看起来还可以,但是积分步数在100万步时,结果看起来是要发散的样子,而且不同次模拟得到的100万步图像是不宜氧的,为什么积分步数一长就无法保持结果的收敛和稳定性,有没有人做过或者有什么看法的可以交流一下。
液体粘度的计算
1万.png


液体粘度的计算-1
10万.png


液体粘度的计算-2
100万.png



[ Last edited by 清风染翠 on 2013-12-16 at 14:16 ]
回复此楼

» 收录本帖的淘帖专辑推荐

MD分子动力学 动力学或热力学计算

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清风染翠

金虫 (小有名气)

直观上的理解应该是积分步数越多,波动越小,越趋于恒定值。我的三个方向边界都有差不多10个纳米,周期性边界会影星很大吗?非平衡方法已经算过了,就是用couette流算剪切力,现在就是想用Green-Kubo或者Einstein公式计算一遍。
5楼2013-12-17 17:09:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
清风染翠: 金币+50, ★★★很有帮助 2013-12-17 17:01:58
月只蓝: 金币+2, 感谢指导! 2013-12-18 08:56:53
时间长的话噪音会有些大,可能只需积分到某个时间就可以了。

这篇文章里谈到了这个问题:http://pubs.acs.org/doi/abs/10.1021/jp211952y
如果想进一步研究的话这篇文章的参考文献可能讲的更详细一些。
2楼2013-12-17 10:03:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清风染翠

金虫 (小有名气)

引用回帖:
2楼: Originally posted by lsloneil at 2013-12-17 10:03:56
时间长的话噪音会有些大,可能只需积分到某个时间就可以了。

这篇文章里谈到了这个问题:http://pubs.acs.org/doi/abs/10.1021/jp211952y
如果想进一步研究的话这篇文章的参考文献可能讲的更详细一些。...

我找到相关文献了,和你说的差不多,谢了.
3楼2013-12-17 11:17:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zlf9452006

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
清风染翠: 金币+30, ★★★很有帮助 2013-12-17 17:02:07
月只蓝: 金币+1, 鼓励交流! 2013-12-18 08:57:07
个人认为是由周期性设置引起的。GK方法本来就是由宏观过程推导出来的,加之其本身信噪比极低,其涨落大是必然的。模拟中的周期性在宏观上导致了液体性质在个别方面的显著变化。
建议楼主考虑下非平衡方法,要稳定很多,也要准确得多。
4楼2013-12-17 12:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +11 妮妮ninicgb 2026-03-15 15/750 2026-03-21 19:15 by ColorlessPI
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见