| 查看: 3869 | 回复: 8 | ||
| 【悬赏金币】回答本帖问题，作者邹卓民将赠送您 5 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
邹卓民新虫 (小有名气)
|
[求助]
lammps采用Green-Kubo计算沥青分子的粘度已有1人参与
|
|
|
lammps采用Green-Kubo计算沥青分子的粘度 出现以下问题：ERROR on proc 0: Bond atoms 3078 3115 missing on proc 0 at step 4275849 (../ntopo_bond_all.cpp:63) #采用Green-Kubo方法计算沥青分子粘度 #采用real单位系统（记住：长度：A；时间：fs） #定义几个参数 variable T equal 533.15 #温度为533.15K variable DT equal 10 #积分步长为10fs variable V equal vol #convert from LAMMPS real units to SI variable kB equal 1.3806504e-23 #J/K Boltzman 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} #初始化 dimension 3 units real boundary p p p atom_style full bond_style harmonic angle_style harmonic dihedral_style harmonic improper_style cvff pair_style lj/cut/coul/long 10.0 kspace_style ewald 1.0e-4 neighbor 2.0 bin read_data data.vigrin region box block -50 100 -50 100 -50 100 velocity all create ${T} 12345 mom yes rot yes dist gaussian #在NVT系综下平衡体系 fix NVT all nvt temp ${T} ${T} 10 #Nose-Hoover热浴 timestep 0.001 thermo 1000 run 100000 #在产出阶段将用NVE系综 unfix NVT fix NVE all nve #在产出阶段计算自关联函数 variable Ns equal 100 #取样间隔 variable Nc equal 100000 #关联数据量 variable Np equal ${Ns}*${Nc} #产出步数 variable pxy equal pxy variable pxz equal pxz variable pyz equal pyz fix SS all ave/correlate ${Ns} ${Nc} ${Np} & v_pxy v_pxz v_pyz type auto file SS.txt ave running variable M equal 1 #粘度自关联函数的个数 variable N equal ${M}*${Np} #跑产出阶段 run ${N} variable scale equal ${convert}/(${kB}*$T)*$V*${Ns}*${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 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" |
回复此楼» 猜你喜欢
博士白读了
已经有45人回复
-
25年博士申请
已经有7人回复
-
381求调剂
已经有5人回复
-
Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？
已经有13人回复
-
考研，求职还是考编?
已经有23人回复
-
"颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人
已经有5人回复
-
博士论文被抄袭
已经有41人回复
-
刚博士毕业，入职普通一本高校
已经有20人回复
-
没学上
已经有8人回复
-
取博导收留
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
邹卓民
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 602.2
- 散金: 18
- 帖子: 209
- 在线: 124.9小时
- 虫号: 4881684
- 注册: 2016-07-30
- 性别: GG
- 专业: 交通工程
6楼2018-06-14 20:58:30
邹卓民
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 602.2
- 散金: 18
- 帖子: 209
- 在线: 124.9小时
- 虫号: 4881684
- 注册: 2016-07-30
- 性别: GG
- 专业: 交通工程
7楼2018-06-14 21:02:28
邹卓民
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 602.2
- 散金: 18
- 帖子: 209
- 在线: 124.9小时
- 虫号: 4881684
- 注册: 2016-07-30
- 性别: GG
- 专业: 交通工程
8楼2018-06-14 22:10:58
5