| 查看: 4033 | 回复: 14 | ||||
[求助]
液体粘度的计算 已有2人参与
|
||||
|
最近在用MD模拟计算一种液体的粘度，使用fortran编的程序，模拟区域是一个长方体，总共有4096个粒子，三个方向都是周期边界，通过热浴保持体系温度恒定。用Green-Kubo公式计算粘度。模拟的参数都按照文献中给定的，所以我是在复现文献中的结果。但是遇到一个问题。 粘度的积分步数在一万步时，得到的结果与文献中符合很好，当积分步数小于十万步时，结果看起来还可以，但是积分步数在100万步时，结果看起来是要发散的样子，而且不同次模拟得到的100万步图像是不宜氧的，为什么积分步数一长就无法保持结果的收敛和稳定性，有没有人做过或者有什么看法的可以交流一下。  1万.png  10万.png  100万.png [ Last edited by 清风染翠 on 2013-12-16 at 14:16 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 MD分子动力学 | 动力学或热力学计算 |
» 猜你喜欢
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 可以让液体粘度变化的材料
已经有5人回复
- 求助：软胶囊胶液水分和粘度的操作方法
已经有7人回复
- 1stopt能进行隐函数非线性拟合吗？谁能帮拟合一下新型液体的粘度混合计算公式
已经有5人回复
- 复数黏度和表观粘度一样吗？
已经有3人回复
- 石油运动粘度的测量
已经有7人回复
- [求助]液体或者水的挥发速度的推导计算公式
已经有7人回复
- 测低粘度溶液的粘度用哪种粘度计比较好，或者说那些
已经有6人回复
- 固体溶解在液体中溶液的密度和粘度
已经有10人回复
- 摇床上的培养瓶中的液体能计算剪切力吗？粗略一点的就行？
已经有5人回复
- 剪切粘度和拉伸粘度有何区别？
已经有7人回复
- 支个招！高粘度高温胶状液体 如何测量粘度呢？水分？
已经有7人回复
- 高炉熔渣粘度问题
已经有6人回复
- 如何降低流体的粘度？
已经有14人回复
- 粘度测量问题(标准物的选择)
已经有9人回复
- 磷酸粘度的测定
已经有5人回复
- 【求助】粘度的定义
已经有5人回复
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
清风染翠: 金币+50, ★★★很有帮助 2013-12-17 17:01:58
月只蓝: 金币+2, 感谢指导！ 2013-12-18 08:56:53
感谢参与，应助指数 +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
6楼2013-12-17 17:09:24
3楼2013-12-17 11:17:11
zlf9452006
金虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 1365.7
- 散金: 539
- 红花: 3
- 帖子: 59
- 在线: 111.1小时
- 虫号: 587860
- 注册: 2008-08-22
- 专业: 理论和计算化学
4楼2013-12-17 12:38:59
5楼2013-12-17 17:09:09
zlf9452006
金虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 1365.7
- 散金: 539
- 红花: 3
- 帖子: 59
- 在线: 111.1小时
- 虫号: 587860
- 注册: 2008-08-22
- 专业: 理论和计算化学
7楼2013-12-18 04:36:33
sciencejoy
新虫 (著名写手)
- 应助: 436 (硕士)
- 金币: 11208.1
- 红花: 91
- 帖子: 2974
- 在线: 498小时
- 虫号: 802149
- 注册: 2009-07-02
- 性别: GG
- 专业: 高分子物理与高分子物理化
8楼2013-12-18 11:15:09
9楼2013-12-18 13:17:31
qdykswang
金虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 6421.1
- 散金: 43
- 红花: 5
- 帖子: 1199
- 在线: 310.1小时
- 虫号: 1382019
- 注册: 2011-08-28
- 性别: GG
- 专业: 数值预报与数值模拟
10楼2013-12-19 07:23:02
