| 查看: 4237 | 回复: 6 | |||
fengguanfeng木虫 (初入文坛)
|
[求助]
lammps自带的例子中用 langevin热浴计算 liquid LJ热导率的问题，有些看不大懂。 已有2人参与
|
|
# sample LAMMPS input script for thermal conductivity of liquid LJ # thermostatting 2 regions via fix langevin # settings variable x equal 10 variable y equal 10 variable z equal 20 variable rho equal 0.6 variable t equal 1.35 variable rc equal 2.5 variable tlo equal 1.0 variable thi equal 1.70 #variable rho equal 0.85 #variable t equal 0.7 #variable rc equal 3.0 #variable tlo equal 0.3 #variable thi equal 1.0 # setup problem units lj atom_style atomic lattice fcc ${rho} region box block 0 $x 0 $y 0 $z create_box 1 box create_atoms 1 box mass 1 1.0 velocity all create $t 87287 pair_style lj/cut ${rc} pair_coeff 1 1 1.0 1.0 neighbor 0.3 bin neigh_modify delay 0 every 1 # heat layers region hot block INF INF INF INF 0 1 region cold block INF INF INF INF 10 11 compute Thot all temp/region hot compute Tcold all temp/region cold # 1st equilibration run fix 1 all nvt temp $t $t 0.5 thermo 100 run 1000 velocity all scale $t unfix 1 # 2nd equilibration run fix 1 all nve fix hot all langevin ${thi} ${thi} 1.0 59804 tally yes # 问题:为什么这里group_id 是all（all不应该是系统所有的原子吗）？是否是因为这里的id是hot的原因，这样跟前面region部分的id是一样的，所以这里的all其实指的就是边界加热浴的部分，并不是系统所有的原子。 fix cold all langevin ${tlo} ${tlo} 1.0 287859 tally yes fix_modify hot temp Thot #问题:修改之前的fix命令的参数，意思是什么？？？ fix_modify cold temp Tcold variable tdiff equal c_Thot-c_Tcold thermo_style custom step temp c_Thot c_Tcold f_hot f_cold v_tdiff thermo 1000 run 10000 # thermal conductivity calculation # reset langevin thermostats to zero energy accumulation # 问题:为什么重设？？ compute ke all ke/atom variable temp atom c_ke/1.5 fix hot all langevin ${thi} ${thi} 1.0 59804 tally yes # 问题:上面这条语句已经出现了，表示加郎之万热浴，这里为什么要再出现一次？？ fix cold all langevin ${tlo} ${tlo} 1.0 287859 tally yes fix_modify hot temp Thot fix_modify cold temp Tcold fix ave all ave/time 10 100 1000 v_tdiff ave running thermo_style custom step temp c_Thot c_Tcold f_hot f_cold v_tdiff f_ave compute layers all chunk/atom bin/1d z lower 0.05 units reduced fix 2 all ave/chunk 10 100 1000 layers v_temp file profile.langevin run 20000 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps运行报错
已经有2人回复
- 一个语句的问题
已经有0人回复
- 为什么在创建原子时有时候会出现多余的键？
已经有3人回复
- 关于lammps自带算例in.crack的几点疑问，望前辈们指教
已经有2人回复
- 关于lammps模拟结果可视化的问题
已经有0人回复
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
2楼2016-08-01 17:59:29
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
3楼2016-08-01 18:01:14
fengguanfeng
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2263
- 帖子: 29
- 在线: 187.4小时
- 虫号: 2060465
- 注册: 2012-10-14
- 性别: MM
- 专业: 基础物理学
4楼2016-08-02 16:36:09
5楼2021-01-04 17:40:15
luxujie3
金虫 (著名写手)
- 应助: 54 (初中生)
- 金币: 4974.1
- 散金: 27
- 红花: 4
- 帖子: 1221
- 在线: 266.5小时
- 虫号: 756922
- 注册: 2009-04-24
- 性别: GG
- 专业: 高分子合成化学
6楼2021-01-16 11:34:47
7楼2022-08-28 19:12:51
40