24小时热门版块排行榜     先丰纳米 买纳米材料就点这里

查看: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 别讨论后缀啦,告诉大家要特别注意的事 +33 氧还原反应 2018-06-14 36/1800 2018-06-19 19:56 by mrequ
[博后之家] 为什么好男人都是别人的? +104 沧海桑田小战士 2018-06-17 133/6650 2018-06-19 19:51 by willy111999
[基金申请] 不要说 +12 天天向上3035 2018-06-19 12/600 2018-06-19 19:15 by 霍元甲
[硕博家园] ZnO XPS能谱的O1s和C1s峰 +3 jamesding 2018-06-14 4/200 2018-06-19 16:09 by 借我一束光
[硕博家园] 我又想退学了 +68 leonsharon 2018-06-14 97/4850 2018-06-19 14:04 by CapHope
[教师之家] 了不起的冰岛队,逼平阿根廷 +14 leiming_0407 2018-06-17 20/1000 2018-06-19 12:42 by leiming_0407
[硕博家园] 色弱研究生求助 +3 有何不可5 2018-06-19 3/150 2018-06-19 10:46 by 胡一恒
[硕博家园] 上海药明康德求合租 +6 噶内特 2018-06-18 13/650 2018-06-19 10:41 by karas_xyf
[土建类考试] 最新的,陈印,陈大师的,葵花宝典!有需要的不? +7 今天是天晴 2018-06-15 9/450 2018-06-19 10:35 by 今天是天晴
[基金申请] 今天教育部人文社科会出结果吗? +46 木小虫.EDU 2018-06-15 105/5250 2018-06-17 09:50 by 改行卖白菜
[硕博家园] 模范导师 +31 科农741 2018-06-16 31/1550 2018-06-17 05:33 by SSTengfei
[找工作] 拟聘公示有影响吗? +4 jenniewang 2018-06-15 5/250 2018-06-16 21:01 by 好林100
[电化学] 纽扣半电池首次充放电 +5 锂电王子 2018-06-13 7/350 2018-06-16 15:54 by kaifeng838
[有机交流] 结晶状叔丁醇怎么变为液体的?(不加水) +5 235zhang 2018-06-15 6/300 2018-06-16 07:04 by Eric Li
[硕博家园] 十年饮冰,难凉热血 +11 杜鲁门.卡波特 2018-06-14 21/1050 2018-06-15 16:16 by SevenˉD
[硕博家园] 读博 +13 Elaine0618 2018-06-14 14/700 2018-06-15 12:42 by ZloveA
[论文投稿] 大家加油! +8 宇宝贝120 2018-06-15 8/400 2018-06-15 01:46 by WY0630
[有机交流] 知道有机物的英文名,如何知道它的CAS号和中文名称?(有机物较复杂,字典查不到) 50+4 李冬冬0415 2018-06-13 12/600 2018-06-14 20:22 by czyzsu
[基金申请] 能登陆项目中后期管理了,里面也有项目申报书,是否代表中了? +27 loveadz 2018-06-13 50/2500 2018-06-14 08:08 by nenugdi
[硕博家园] 联合培养对读博利弊 +5 18039292882 2018-06-13 10/500 2018-06-13 15:57 by YQ_coco
信息提示
请填处理意见