计算化学公社

标题: 求助，如何根据应力自相关函数求粘度？ [打印本页]

作者
Author: truemantwo 时间: 2015-5-21 18:08
标题: 求助，如何根据应力自相关函数求粘度？
各位大侠：

在分子动力学模拟里，可以根据应力自相关函数的积分求出粘度。

问题是：这个积分如何用程序实现？

附件中是原始文献，请各位大侠帮忙，谢谢！

作者
Author: sobereva 时间: 2015-5-21 20:34
应该就是对图2的曲线进行积分。用梯形法积分就可以，所有数值算法的书里讲数值积分的地方都会讲到。简单来说对每个t值都要算一次自相关函数，然后乘上步长dt，然后加和。

作者
Author: EricCQU 时间: 2022-10-30 19:36
楼主，这个问题解决了吗？我模拟出来和实际差很多

欢迎光临计算化学公社 (http://bbs.keinsci.com/)