| 查看: 3900 | 回复: 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 |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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 (初中生)
- 金币: 4952.7
- 散金: 27
- 红花: 4
- 帖子: 1213
- 在线: 266.2小时
- 虫号: 756922
- 注册: 2009-04-24
- 性别: GG
- 专业: 高分子合成化学
6楼2021-01-16 11:34:47
7楼2022-08-28 19:12:51