| 查看: 4290 | 回复: 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 |
回复此楼» 猜你喜欢
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
现代”学阀”该如何界定
已经有8人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
实名举报：华南理工大学物理与光电学院副院长李志远婚内两次出轨女学生（博士）
已经有16人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有12人回复
-
我的奶奶
已经有3人回复
-
国社科系统bug了，是不是要放榜了？
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps运行报错
已经有2人回复
- 一个语句的问题
已经有0人回复
- 为什么在创建原子时有时候会出现多余的键？
已经有3人回复
- 关于lammps自带算例in.crack的几点疑问，望前辈们指教
已经有2人回复
- 关于lammps模拟结果可视化的问题
已经有0人回复
7楼2022-08-28 19:12:51
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