应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps MP法求热导率,为何体系的温度下降了。
21/1返回列表
查看: 1364  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

雪雨星风17

木虫 (小有名气)

无

[求助] lammps MP法求热导率,为何体系的温度下降了。

采用lammps的MP方法求热导，关键部分代码如下：
units             metal
variable            T equal 300
velocity          all  create  $T 458127641 mom yes  rot yes dist gaussian units box
fix               1 all  temp/berendsen $T $T 0.001
compute           ke  all  ke/atom
variable          temp atom  c_ke/(1.5*1.3806504e-4)*1.602763
fix               temp_profile    all    ave/spatial  1  100000  100000  z  lower  2    v_temp  file  temp.profile    units  lattice
thermo_style      custom step temp etotal vol
thermo_modify     lost warn
thermo            100
# Run
dump                         id all atom 10000 dump.lammpstrj
timestep          0.001
run               100001
unfix             1
fix               nve  all  nve
fix               heat_swap   all  thermal/conductivity  10  z  400
fix               e_exchange  all  ave/time  10  10000  100000  f_heat_swap  file  e_exchange.dat   
run               10000000
开始原子的温度设定为300K,热浴之后得到的结果   99600          300   -132337.22     12788141
   99700          300   -132337.22     12788141
   99800          300   -132337.22     12788141
   99900          300   -132337.22     12788141
  100000          300   -132337.22     12788141
  100001          300   -132337.22     12788141
紧接着进行热导计算的结果
  100001          300   -132337.22     12788141
  100100    214.80808   -132337.05     12788141
  100200    198.52124   -132336.87     12788141
  100300     221.0091   -132336.93     12788141
  100400    264.85103   -132336.99     12788141
  100500    270.57875   -132336.95     12788141
  100600    304.28282   -132337.09     12788141
在连续的波动之后，温度在  246K左右稳定了，请问大神这是什么原因？
445400    246.03241   -132336.74     12788141
  445500    245.74197   -132336.74     12788141
  445600    246.37433   -132336.74     12788141
  445700    244.57574   -132336.74     12788141
  445800    246.93508   -132336.74     12788141
  445900    245.90909   -132336.74     12788141
  446000    246.30204   -132336.74     12788141
  446100    246.08624   -132336.74     12788141
  446200    244.93367   -132336.74     12788141
  446300    245.27753   -132336.74     12788141
  446400    245.15068   -132336.74     12788141
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
可以可以啊
1楼 2015-04-10 14:40:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

电容2016

新虫 (初入文坛)
楼主问题解决了吗？我出现的问题和你的类似，而且fix  thermal/conductivity  10  z  400的结果为零，不知道什么原因。如果楼主问题解决了，希望分享一下经验！万分感谢！
赞一下(1人) 回复此楼
2楼2016-01-19 19:23:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 雪雨星风17 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我，科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 3/150 2026-02-08 03:05 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 3/150 2026-02-08 02:52 by vs90ilomwc
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 4/200 2026-02-08 02:45 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 3/150 2026-02-08 02:32 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:27 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:25 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 02:07 by vs90ilomwc
[教师之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 02:05 by vs90ilomwc
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 01:46 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 01:32 by vs90ilomwc
[教师之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 01:26 by vs90ilomwc
[教师之家] 有院领导为了换新车，用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人） +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
关闭
请填处理意见
关闭