| 查看: 3607 | 回复: 7 | ||
[求助]
lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
|
|
我是lammps新手，在学校做一个关于纳米管填充的计算，我只有计算同一原子体系热导率的in文件，代码如下： units metal atom_style atomic #neighbor 3.0 nsq boundary s s p read_data data.55 pair_style tersoff pair_coeff * * SiC.tersoff C dump ribbon all custom 10000 all.xyz type x y z velocity all create 300 17287 loop geom minimize 0.0 0.0 5000 50000 compute r1 all temp compute r2 all pe compute r3 all ke variable r4 equal c_r2+c_r3 fix a1 all ave/time 2 100 1000 c_r1 c_r2 c_r3 v_r4 file relax.dat thermo 1000 thermo_style custom step temp ke pe etotal fix s1 all nvt temp 300.0 300.0 10 drag 0.2 timestep 0.0002 run 340000 unfix s1 fix s2 all nve timestep 0.0002 run 60000 write_restart fnve.%.* unfix a1 unfix s2 fix s3 all nve fix s4 all thermal/conductivity 80 z 50 swap 1 compute ke all ke/atom variable temp atom c_ke/0.0001292355 fix a2 all ave/spatial 5 200 1000 z 0.0 0.02 v_temp file tmp.profile units reduced fix a3 all ave/atom 5 200 1000 x y z dump a3 all custom 1000 dump.position id type x y z fix a4 all ave/atom 5 200 1000 vx vy vz dump a4 all custom 1000 dump.velocity id type vx vy vz thermo_style custom step temp press vol ke f_s4 thermo 10000 timestep 0.0002 run 200000 unfix s4 unfix s3 不过这个只适合单原子体系，如果要计算多原子体系热导率，要用到fix heat指令，请问各位大神该怎么修改？求指点 先谢谢了！ |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 学习第一性原理计算的准备
已经有13人回复
- MS的CASTEP计算一定要联网吗？
已经有4人回复
- 原子吸收调取以前测的标线的问题
已经有13人回复
- 不同压力下的焓值计算
已经有8人回复
- 在Lammps中安装xmovie工具时，Makefile文件怎么修改X库路径
已经有3人回复
- 【求助】LAMMPS-运行报错
已经有9人回复
- 【转帖】热导率计算的in文件
已经有13人回复
- 【求助】将MS数据文件转化成lammps data
已经有25人回复
- 【求助】关于LAMMPS 运行的问题
已经有6人回复
- 【求助】MS建模导入lammps问题求救
已经有14人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
sdzjc1988
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 1695.4
- 红花: 4
- 帖子: 208
- 在线: 237.2小时
- 虫号: 1164293
- 注册: 2010-12-07
- 专业: 工程热力学
2楼2012-04-01 14:25:11
|
这个in文件的计算方法是通过对热端最低温原子和冷端最高温原子进行互换地方法产生热流。具体方法为：在热流方向将模拟结构分为N层，选取第1层作为冷端，第N/2+1层作为热端，只在径向方向上采用周期性边界条件，然后从冷端选n个最高温原子（n一般取1、2），从热端选n个最低温原子，交换原子的速度（如图一），这样会使热端更热、冷端更冷，但是由于结构自身的导热能量从热端流向冷端，当两者达到平衡时管子将会呈现一定的温度梯度▽T，而且此时热端和冷端每次交换的能量恰好等于两次交换之间的热流J。然后利用公式k=-J/▽T便可求得热导率。这种方法只适合计算单原子体系这种方法只适用于单个碳纳米管或者单一原子体系，由于Muller-Plathe方法仅仅是交换原子的速度，如果交换原子的原子质量不同，那么会导致交换的时候能量不守恒，从而使得系统能量越来越低。另一种方法是，先将模拟结构分成n层，然后在其中一层加能量（该层做热端），另一层抽相同的能量（该层做冷端），这种方法会用到fix heat指令吧，不知道该怎么修改 |
3楼2012-04-09 15:46:31
4楼2012-04-19 11:27:29
5楼2012-04-19 11:28:55
sdzjc1988
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 1695.4
- 红花: 4
- 帖子: 208
- 在线: 237.2小时
- 虫号: 1164293
- 注册: 2010-12-07
- 专业: 工程热力学
6楼2012-04-25 12:21:18
tanxjezh
铁杆木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 6847.5
- 红花: 2
- 帖子: 77
- 在线: 271.2小时
- 虫号: 605576
- 注册: 2008-09-18
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2012-05-17 16:54:03
8楼2013-05-07 09:20:18