| 查看: 3826 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
热导率计算问题 已有1人参与
|
||
|
我是计算一个石墨烯纳米带300k下的热导率 下面是我的in文件 不知道为什么体系就是达不到平衡 整体温度一直在285到315之间震荡 局部温度震荡范围更大 请问我该怎么做 是什么问题 怎么样才能让体系弛豫达到平衡状态呢 求指教！ 附件是计算的log文件 boundary s s s units metal dimension 3 atom_style atomic atom_modify map array neighbor 2.0 nsq read_data data.graphene1550 pair_style airebo 3.0 1 0 pair_coeff * * CH.airebo C region box block 0.0 112.311 0.0 36.900 0.0 0.5 units box region begin block 7.21 8.63 INF INF INF INF units box region hot block 8.63 17.152 INF INF INF INF units box region a block 17.152 25.674 INF INF INF INF units box region b block 25.674 34.195 INF INF INF INF units box region c block 34.195 42.717 INF INF INF INF units box region d block 42.717 51.239 INF INF INF INF units box region e block 51.239 59.76 INF INF INF INF units box region e1 block 59.76 68.282 INF INF INF INF units box region d1 block 68.282 76.804 INF INF INF INF units box region c1 block 76.804 85.325 INF INF INF INF units box region b1 block 85.325 93.847 INF INF INF INF units box region a1 block 93.847 102.369 INF INF INF INF units box region cold block 102.369 110.891 INF INF INF INF units box region end block 110.891 112.311 INF INF INF INF units box region boundary union 2 begin end units box group begin region begin group hot region hot group a region a group b region b group c region c group d region d group e region e group e1 region e1 group d1 region d1 group c1 region c1 group b1 region b1 group a1 region a1 group cold region cold group end region end group boundary region boundary group 1 subtract all boundary velocity 1 create 300 7986875 rot yes dist gaussian units box velocity begin set 0.0 0.0 0.0 units box velocity end set 0.0 0.0 0.0 units box fix 1 begin setforce 0.0 0.0 0.0 fix 2 end setforce 0.0 0.0 0.0 fix a 1 nvt temp 300 300 0.025 compute ke 1 ke/atom variable temp atom c_ke*1.6e+4/(1.5*1.38) fix temp_profile 1 ave/spatial 1 100000 100000 x (0,0,0) 2.10202 v_temp file temp.profile units box compute c_temp c temp compute c1_temp c1 temp compute 1_temp 1 temp variable delta_temp equal c_c_temp-c_c1_temp fix delta_out 1 ave/time 1 100000 100000 v_delta_temp file delta_temp.dat thermo_style custom step c_1_temp etotal pe ke c_c_temp c_c1_temp #thermo_modify lost warn timestep 0.0005 thermo 1000 dump 1 all xyz 100 g.xyz unfix a fix 3 1 nve fix hot all heat 1 1 region hot fix cold all heat 1 -1 region cold variable thermal_conductivity equal 1*1.6e-19*40.448*1e-10/(36.9*3.4*1e-20*f_delta_out*1e-12) fix thermal_conductivity_out all ave/time 100000 1 100000 v_thermal_conductivity file thermal_conductivity.dat #run 5000000 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : log.lammps
2012-06-28 20:29:51, 30.64 K
» 猜你喜欢
303求调剂
已经有6人回复
-
306求调剂
已经有4人回复
-
求调剂
已经有6人回复
-
275求调剂
已经有5人回复
-
求调剂
已经有7人回复
-
考研调剂
已经有4人回复
-
285求调剂
已经有6人回复
-
一志愿北京化工大学070300 学硕336求调剂
已经有7人回复
-
318求调剂
已经有6人回复
-
305分求调剂（食品工程）
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
已经有7人回复
- 钢包用耐火材料壁厚的设计计算
已经有9人回复
- 请教关于多元（三种不同成分以上）复合材料热导率和热膨胀系数的问题
已经有5人回复
- SiC的热导率问题
已经有9人回复
- 复合材料比热容的确定
已经有16人回复
- 【求助】求解关于测热导率试样尺寸问题
已经有3人回复
- 【求助】第一性原理计算材料导电、导热、塞贝克系数问题
已经有18人回复
- 【资源】热物性的瞬态测量
已经有5人回复
tang_wei
金虫 (正式写手)
- 应助: 24 (小学生)
- 金币: 1353.5
- 红花: 28
- 帖子: 428
- 在线: 81.7小时
- 虫号: 4067192
- 注册: 2015-09-12
- 专业: 高分子物理与高分子物理化
5楼2015-11-20 09:24:56
维特的日记
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 3814.9
- 红花: 2
- 帖子: 169
- 在线: 55.5小时
- 虫号: 3159711
- 注册: 2014-04-24
- 专业: 工程热物理相关交叉领域
2楼2014-08-20 22:23:49
牛犊小子
铁杆木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 8786.5
- 散金: 275
- 红花: 1
- 帖子: 1165
- 在线: 169.8小时
- 虫号: 2659995
- 注册: 2013-09-16
- 性别: GG
- 专业: 传热传质学
4楼2014-08-27 09:35:53
6楼2015-12-02 19:57:34
