版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂

返回列表

【悬赏金币】回答本帖问题，作者Mayerss将赠送您 20 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Mayerss

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 115
帖子: 63
在线: 22小时
虫号: 3218673
注册: 2014-05-19
专业: 工程热物理相关交叉领域

[求助] 关于Green-Kubo方法模拟水的粘度已有3人参与

求问有没有大神用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_style harmonic

pair_style    lj/cut/tip4p/long 1 2 1 1 0.125 12 15
kspace_style pppm/tip4p 1e-6
read_data    sketch.data

pair_coeff  1 2  0.0 0.0 # O-H
pair_coeff  2 2  0.0 0.0 # H-H
pair_coeff  1 1  0.16275 3.16435 # O-O

timestep    ${dt}
thermo       $d

# equilibration and thermalization

velocity    all create $T 102486 mom yes rot yes dist gaussian
fix          NVT all nvt temp $T $T 40 drag 0.2
run          100000

# viscosity calculation, switch to NVE if desired

#unfix       NVT
#fix          NVE all nve

reset_timestep 0
variable    pxy equal pxy
variable    pxz equal pxz
variable    pyz equal pyz
fix          SS all ave/correlate $s $p $d &
            v_pxy v_pxz v_pyz type auto file S0St.dat ave running
variable    scale equal ${convert}/(${kB}*$T)*$V*$s*${dt}
variable    v11 equal trap(f_SS[3])*${scale}
variable    v22 equal trap(f_SS[4])*${scale}
variable    v33 equal trap(f_SS[5])*${scale}
thermo_style custom step temp press v_pxy v_pxz v_pyz v_v11 v_v22 v_v33
run          200000
variable    v equal (v_v11+v_v22+v_v33)/3.0
variable    ndens equal count(all)/vol
print       "average viscosity: $v [Pa.s/ @ $T K, ${ndens} /A^3"

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
附件 1 : log.lammps

2015-05-23 19:35:01, 21.42 K

附件 2 : sketch.data

2015-05-23 19:35:02, 361.88 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

请教平衡分子动力学模拟和非平衡分子动力学模拟的区别已经有11人回复
25°水蒸气粘度已经有4人回复
Green-Kubo公式计算热导率，自相关函数问题已经有18人回复
液体粘度的计算已经有14人回复
Green一Kubo公式的应用已经有5人回复
Green-Kubo线性响应理论应该是在平衡态用还是非平衡态？已经有5人回复
1stopt能进行隐函数非线性拟合吗？谁能帮拟合一下新型液体的粘度混合计算公式已经有5人回复
各位大侠请问用fluent仿真液固两相流时固相的粘度需要设置吗? 已经有12人回复
湍流粘度的意义已经有10人回复
MD计算粘度，结果分析不稳定已经有17人回复
计算精度不高，且出现湍流粘度限制怎么办已经有8人回复
海水的动力粘度u是多少？哪位大侠知道，，，十分感谢已经有10人回复
请问：剪切粘度（shear viscosity）和容积粘度（bulk viscosity）的区别已经有5人回复
【求助】在计算的时候出现湍流粘度比受到限制已经有12人回复
【资源】关于圆管内气体流场中雷诺数与气流速度及动力粘度的问题已经有9人回复
【讨论】Green-Kubo公式是基于第一性原理的吗？？已经有4人回复

1楼2015-05-23 19:35:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 Mayerss 的主题更新