查看: 250  |  回复: 5
【悬赏金币】回答本帖问题,作者Mayerss将赠送您 20 个金币

Mayerss

新虫 (初入文坛)

[求助] 关于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"
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:emuch2018@163.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

» 猜你喜欢

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

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

kakaxicheng

新虫 (初入文坛)

【答案】应助回帖

哥们,应该不需要把real转换成SI 这一步,即不要单位转换的因子convert
2楼2016-10-04 15:58:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

不要光看最后一个结果,要仔细分析跑动粘度,即粘度随关联时间的变化情况。通常,跑动粘度需要一定的关联时间才能收敛。另外,一次模拟也不够,因为算粘度误差会很大,每次得到的结果也会差别较大。

发自手机版WAP
3楼2016-10-05 01:10:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kakaxicheng

新虫 (初入文坛)

【答案】应助回帖

水分子选的模型会影响剪切粘度:
TIP3P TIP4P TIP5P SPC/E TIP4P/2005 Expt.a
0.321 0.494  0.699 0.729  0.855           0.896
4楼2016-10-25 09:56:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xu_jiangxin

新虫 (初入文坛)

【答案】应助回帖

一楼的同学说不用转化为SI ,这个我保留意见。我在自己的模拟中使用了单位转化。关于你的data文件和in文件我说一下我的看法:你使用的是tip4p模型,却使用了很多的spce模型的参数而不是tip4p模型的参数。
这是我自己的模拟结果:(300K,1.29E+5mol/m3)
spce            tip4p       tip4p2005
0.567          0.420          0.721
使用的参数为lammps手册“6.8 TIP4P water mode”中给出的参数,l
5楼2017-08-28 09:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

SDXM

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by xu_jiangxin at 2017-08-28 09:13:22
一楼的同学说不用转化为SI ,这个我保留意见。我在自己的模拟中使用了单位转化。关于你的data文件和in文件我说一下我的看法:你使用的是tip4p模型,却使用了很多的spce模型的参数而不是tip4p模型的参数。
这是我自 ...

你好,方便把您的data文件和input文件给我传一份吗?我也在做这个模拟,数据都来在6.8 TIP4P water mode,但是最终的结果却对不上,谢谢您啦 我的邮箱是  1178408727@qq.com  非常感谢
6楼2018-03-14 04:04:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Mayerss 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[留学DIY] 美国密苏里大学 材料力学及先进制造 (能源及医疗应用) 课题组招2-3名全奖博士生 +3 wadlut1123 2018-09-20 4/200 2018-09-21 06:49 by wadlut1123
[文学芳草园] ex结婚,难过你说我矫情,不难过你说我没用真心 +3 斯往 2018-09-19 3/150 2018-09-20 23:39 by txq272
[硕博家园] 假若爱有天意,有点想要追求你 +26 樱花翩翩 2018-09-19 91/4550 2018-09-20 23:28 by zsl668013
[职场人生] 高校行政不太理想,何去何从。 +36 夏蒟蒻just 2018-09-16 86/4300 2018-09-20 23:01 by 夏蒟蒻just
[考博] 32岁高龄,Mba毕业,考博有竞争力吗,出来以后有前途吗 +12 武汉小博士 2018-09-14 14/700 2018-09-20 21:44 by tmgczb
[有机交流] 手都在颤抖! +30 t-Bu型心态 2018-09-18 40/2000 2018-09-20 20:22 by 啄木鸟先森
[论文投稿] IEEE ACCESS重投祈福 +4 tendertime 2018-09-20 4/200 2018-09-20 16:20 by wang31aa
[考博] 清华复试询问导师结果导师回复 +3 wmdzcsxcdh 2018-09-20 6/300 2018-09-20 16:13 by pourchopin
[论文投稿] osa审稿人状态看不到,这是为什么? +5 Hanlaser 2018-09-15 9/450 2018-09-20 13:02 by qh2008723
[基金申请] 教育部人文社会科学课题 +6 fengchunshu 2018-09-14 7/350 2018-09-20 11:51 by Peter_Gu
[硕博家园] 利用Materials Studio怎么将两个物质放到一起 +3 有梦为马333 2018-09-16 10/500 2018-09-20 08:23 by 有梦为马333
[硕博家园] 各位有没有遇到这种让人窝火的事。 +7 Sweet_135 2018-09-19 7/350 2018-09-20 00:26 by whitebox
[有机交流] 苯胺与酸酐的反应 10+3 (☆_☆) 2018-09-18 4/200 2018-09-19 10:54 by binzy
[论文投稿] 修改后,编辑上午直接又送审啦 +8 卜卜卜卜卜 2018-09-16 11/550 2018-09-18 12:43 by 若冲93
[找工作] 青岛石油大学教师待遇如何 +3 萤火虫27 2018-09-16 5/250 2018-09-18 07:57 by 萤火虫27
[论文投稿] 文章被接收,开心 (金币+5) +76 化学小星星 2018-09-15 81/4050 2018-09-18 06:41 by QHJ100100
[考研] 考研狗 +4 hmily1 2018-09-16 7/350 2018-09-17 23:36 by 易爆是官配
[电化学] 电解液无法润湿隔膜 3+4 xuxudeshuai 2018-09-14 13/650 2018-09-17 16:22 by fuxun1979
[硕博家园] 选择困难症 +15 路在脚下z 2018-09-16 17/850 2018-09-16 22:31 by winner3
[硕博家园] 博士是一般只有一两个课题吗 +5 第三个 2018-09-15 7/350 2018-09-15 12:09 by 硬币啊硬币
信息提示
请填处理意见