24小时热门版块排行榜

返回列表

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

zjnu_yc

银虫 (小有名气)

应助: 3 (幼儿园)
金币: 423.9
散金: 16
红花: 2
帖子: 206
在线: 148.5小时
虫号: 1379993
注册: 2011-08-25
专业: 零件成形制造

[求助] 求助LAMMPS两次运算结果不一致问题！已有1人参与

本人打算采用LAMMPS计算高聚物的导热系数，由于是新手，首先按照LAMMPS例子里面的IN文件进行了练习（文件路径：LAMMPSExampleKAPPAin.heat），运算结果与所提供的参照结果比较吻合。而考虑到将来我的研究对象是复杂的高聚物体系，应该要用到MS来建模，为了先探探路，于是先将上述in.heat文件内描述的模型应用MS进行建模，建好模型后在Discover模块进行设置（energy里面选择compass），随后typing，并顺利导出、转换成LAMMPS所能识别的data文件（只有选择了compass才能顺利导出，后来在data文件中修改了pair_coeff）。
我碰到的问题是：从MS-->LAMMPS的方法计算结果仅为参照结果的一半（即按照LAMMPS例子中提供的in文件的计算结果）。我仔细对比过MS产生data文件与LAMMPS提供的in文件的设置，关于模型的参数（如原子数目、pair_style、pair_coeff、Mass等参数均相同），以及后续对物理过程的处理代码都是相同的，可是就是结果差别太大。（附上两次运算的LOG文件，其中Neighs数字有较大区别）
特请这方面经验丰富的朋友帮忙出出主意，或者提点宝贵意见，先感谢啦！@oxox6085@smutao

回复此楼

» 猜你喜欢

1楼 2014-03-22 11:07:22

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

ysshj

金虫 (小有名气)

应助: 3 (幼儿园)
金币: 669.6
红花: 4
帖子: 111
在线: 43.3小时
虫号: 2237709
注册: 2013-01-12
性别: GG
专业: 机械摩擦学与表面技术

★ ★ ★ ★ ★
月只蓝: 金币+5, 感谢指导 2014-04-13 09:28:51

引用回帖:

11楼: Originally posted by zjnu_yc at 2014-04-12 11:58:28
您好！您能将您发给我的 in.hydrogen 和 hydrogen.data文件完整运行一下吗？（不只是弛豫部分，包括后面的run），我想比较下最终的结果与例子是否接近。我的机子上一直运行不了，我现在又出现了新的问题（见附件LO ...

调试的话还是把in文件和data文件也发一下吧。
你用的什么版本的lammps?

logMS2与logMS的in文件相比，改动：

lattice fcc ${rho}
region hot block INF INF INF INF 0 1
region cold block INF INF INF INF 9.5 11

1 units lj下尺寸换算有点闹心，我不习惯。但是涉及尺寸的命令很可能影响结果。
2 可以把data文件里的原子z坐标都加0.470517500，再改in文件试试。