| 查看: 1393 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者CarbonO2将赠送您 20 个金币 | ||
[求助]
lammps nemd计算热导率问题
|
||
|
大家好，我在使用nemd计算带有凸起的石墨烯热导率时（冷热源均取在凸起上），出现凸起部分有较大温度变化，但是主体部分温度几乎不变的情况，请问该如何解决该问题？以下为我的in文件 units metal dimension 3 newton on boundary p s s atom_style full neighbor 2.0 nsq neigh_modify delay 2 check yes #processors 4 4 1 #echo log # read data of graphite read_data Result_With_Pillar_y.data #read_restart restart.nve #reset_timestep 0 # define groups # total=5112 group downb id 1:72 group upb id 5041:5112 group source id 73:216 group sink id 4897:5040 group caculate subtract all upb downb source sink group tot subtract all upb downb compute avelength caculate chunk/atom bin/1d y lower 3 # tersoff potential pair_style tersoff pair_coeff * * SiC.tersoff C C #2 atoms types velocity all create 50.0 65877893 mom yes rot yes dist gaussian units box min_style sd minimize 1.0e-10 1.0e-10 10000 10000 dump 0 all custom 500000 dump.*.lammpstrj id type x y z dump 1 all xyz 100000 dump.xyz thermo 50000 thermo_style custom step temp pe ke etotal vol lx ly lz fix 2 upb setforce 0.0 0.0 0.0 fix 3 downb setforce 0.0 0.0 0.0 timestep 0.0005 fix momentum all momentum 1 linear 1 1 1 angular #fix npt all npt temp 50 50 1 iso 1.0 1.0 1 #run 1000000 #write_restart restart.npt #unfix npt fix nvt tot nvt temp 50 50 0.05 drag 0.2 run 2000000 write_restart restart.nvt unfix nvt fix nve tot nve fix 9 all temp/berendsen 50 50 1.0 #fix nvt all langevin 50 50 0.1 36102036 run 2000000 unfix nve write_restart restart.nve fix nve1 tot nve fix 5 source heat 1 1.5 fix 6 sink heat 1 -1.5 fix 7 caculate ave/chunk 50 100 5000 avelength temp file Tprofile.txt run 2000000@smutao@oxox6085 |
回复此楼» 猜你喜欢
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
2026国自然函评费到账
已经有16人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2021-03-19 09:52:35
40