关于Green-Kubo方法模拟水的粘度

最新推荐文章于 2024-07-25 17:07:39 发布

ditigao0965

最新推荐文章于 2024-07-25 17:07:39 发布

阅读量1.5k

点赞数

文章标签： sketch

原文链接：http://www.cnblogs.com/Simulation-Campus/p/8743104.html

版权

博主尝试使用LAMMPS和Green-Kubo方法模拟300K下水的粘度，得到结果0.2mPa.s，远低于标准值0.8mPa.s。讨论涉及模拟设置、单位转换、模型选择（如TIP4P）及误差分析，指出模拟需要考虑关联时间和多次运行以减小误差，并分享了不同水分子模型的粘度比较。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

关于Green-Kubo方法模拟水的粘度

求问有没有大神用lammps做过水的模拟？
我现在改编了一下manual里面计算液氩粘度的那个算例，来模拟水的粘度，300k下得到的结果是0.2mPa.s，查得的标准值是0.8mPa.s，请路过的大神帮忙看一下是哪里出了问题，多谢！

in文件贴在下面，data文件和log文件在附件中。

# LAMMPS input script for viscosity of purewater

units real
variable T equal 300
variable V equal vol
variable dt equal 0.4
variable p equal 400 # correlation length
variable s equal 5 # sample interval
variable d equal $p*$s # dump interval

# convert from LAMMPS real units to SI

variable kB equal 1.3806504e-23 # [J/K/ Boltzmann
variable atm2Pa equal 101325.0
variable A2m equal 1.0e-10
variable fs2s equal 1.0e-15
variable convert equal ${atm2Pa}*${atm2Pa}*${fs2s}*${A2m}*${A2m}*${A2m}

# setup problem

dimension 3
boundary p p p
atom_style full
bond_style harmonic
angle_st