计算化学公社

标题: lammps中Muller-Plathe法计算粘度相关问题 [打印本页]

作者
Author: chemdumpling 时间: 2023-10-31 14:56
标题: lammps中Muller-Plathe法计算粘度相关问题
本人最近刚刚接触MD模拟以及lammps，在网上看模拟到计算粘度相关的教程之后存在如下问题：一般viscocity计算法通常有Muller-Plathe法以及Green-kubo法，但目前来看对于初期学习，Muller-Plathe法较为简便，但是现有教程感觉都并无完整计算流程，本人也参照论文原文Reversing the perturbation in nonequilibrium molecular dynamics: An easy way to calculate
the shear viscosity of fluids进行相关学习，计算粘度时一般参照下图公式：
(, 下载次数 Times of downloads: 31)
但教程中粘度计算部分如下：
compute layers all chunk/atom bin/1d z center 0.05 units reducedfix 4 all ave/chunk 20 50 1000 layers vx file velocity_profile.datfix 5 all viscosity 100 x z 20fix 6 nvevariable dVx equal f_4[11][3]-f_4[1][3]thermo 1000thermo_style custom step temp press etotal v_dVx f_5 run 50000

按照本人理解，此过程仅涉及δvx/δz的计算，对其他参量并无描述，请问各位是我的理解有误还是他确实没有将粘度算完？如果部分缺失的话，公式比值上方部分参量应如何计算？

作者
Author: lue611 时间: 2023-12-28 16:15
请问楼主解决这个问题了吗

作者
Author: chemdumpling 时间: 2024-1-1 21:39

lue611 发表于 2023-12-28 16:15
请问楼主解决这个问题了吗

基本解决了，不过也是搬运别人的代码

作者
Author: mumuwang 时间: 2024-1-25 16:32
本帖最后由 mumuwang 于 2024-1-25 16:36 编辑

你好，想加一下vx，请教一下Muller-Plathe方法粘度的计算，v-x是zls1137132886，非常感谢

作者
Author: SMJ999 时间: 2024-4-16 20:20

chemdumpling 发表于 2024-1-1 21:39
基本解决了，不过也是搬运别人的代码

你好！可以交流一下MP计算粘度的方法么？可以添加一下VX么？VX：18392102204

作者
Author: 759884279 时间: 2024-5-20 08:45

chemdumpling 发表于 2024-1-1 21:39
基本解决了，不过也是搬运别人的代码

你好，可以请教一下Muller-Plathe方法粘度的计算，VX：zxy759884279，望交流

作者
Author: Rainy_123 时间: 2024-7-1 20:35

chemdumpling 发表于 2024-1-1 21:39
基本解决了，不过也是搬运别人的代码

您好，请问能请教一下Muller-Plathe方法计算粘度吗，v-x是Rainy-YH，拜托了

作者
Author: chemdumpling 时间: 2024-7-2 16:38

Rainy_123 发表于 2024-7-1 20:35
您好，请问能请教一下Muller-Plathe方法计算粘度吗，v-x是Rainy-YH，拜托了

本人实践目前看来Muller-Plathe法存在每次在不同初始温度随机数测定下，粘度差距较大，差距大到随温度变化规律性非常差，就算多次测定取平均也只是稍微好一些，你可以试试看大佬的帖子，采用自扩散系数计算粘度，具体应该是爱因斯坦-斯托克斯公式，站里有相关的帖子。

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