24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3408  |  回复: 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的回帖

sodendrimer

新虫 (初入文坛)

correlation length 的大小对粘滞系数的结果影响这么大,不是很明白,
15楼2015-12-23 20:57:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 清风染翠 的主题更新
信息提示
请填处理意见