| 查看: 1255 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
ashioulou铜虫 (初入文坛)
|
[求助]
NEMD熱傳導性質 NVE能量不守衡問題
|
|
|
小弟目前在做複合材料的熱傳導性質，需要用到NVE系綜做系統平衡 但理論上的NVE為系統總能量需守?a，但目前小弟一直無法做到這點 在系統平衡時溫度會一直上升，能量也會越?碓酱螅恢焙茈y維持在300K 爬了一些文章，發現此種問題可能是跟lammps本身在做NVE積分的時候， 有些誤差累積所造成能量疊加上升。 但小弟的case上升的非常誇?垼瑢n文件貼出，希望有先進可以指點一二 ================================================= units real dimension 3 boundary p p p atom_style full kspace_style pppm 1e-6 pair_style lj/class2/coul/long 10.0 10.0 pair_modify mix sixthpower bond_style class2 angle_style class2 dihedral_style class2 improper_style class2 neighbor 3.0 nsq read_restart composite.restart.1200000 region 1 block INF INF INF INF 0 2.4740 units box region 2 block INF INF INF INF 2.4740 4.9480 units box region 3 block INF INF INF INF 4.9480 7.4220 units box region 4 block INF INF INF INF 7.4220 9.8960 units box region 5 block INF INF INF INF 9.8960 12.370 units box region 6 block INF INF INF INF 12.370 14.844 units box region 7 block INF INF INF INF 14.844 17.318 units box region 8 block INF INF INF INF 17.318 19.792 units box region 9 block INF INF INF INF 19.792 22.266 units box region 10 block INF INF INF INF 22.266 24.740 units box region 11 block INF INF INF INF 24.740 27.214 units box region 12 block INF INF INF INF 27.214 29.688 units box region 13 block INF INF INF INF 29.688 32.162 units box region 14 block INF INF INF INF 32.162 34.636 units box region 15 block INF INF INF INF 34.636 37.110 units box region 16 block INF INF INF INF 37.110 39.584 units box region 17 block INF INF INF INF 39.584 42.058 units box region 18 block INF INF INF INF 42.058 44.532 units box region 19 block INF INF INF INF 44.532 47.006 units box region 20 block INF INF INF INF 47.006 49.480 units box region 21 block INF INF INF INF 49.480 51.954 units box region 22 block INF INF INF INF 51.954 54.428 units box region 23 block INF INF INF INF 54.428 56.902 units box region 24 block INF INF INF INF 56.902 59.376 units box region 25 block INF INF INF INF 59.376 64.324 units box region 26 block INF INF INF INF 121.226 123.7 units box region cold union 2 1 26 units box group source region 25 group sink region cold group 1 region 2 group 2 region 3 group 3 region 4 group 4 region 5 group 5 region 6 group 6 region 7 group 7 region 8 group 8 region 9 group 9 region 10 group 10 region 11 group 11 region 12 group 12 region 13 group 13 region 14 group 14 region 15 group 15 region 16 group 16 region 17 group 17 region 18 group 18 region 19 group 19 region 20 group 20 region 21 group 21 region 22 group 22 region 23 group 23 region 24 fix 1 all nve compute temps sink temp compute temp1 1 temp compute temp2 2 temp compute temp3 3 temp compute temp4 4 temp compute temp5 5 temp compute temp6 6 temp compute temp7 7 temp compute temp8 8 temp compute temp9 9 temp compute temp10 10 temp compute temp11 11 temp compute temp12 12 temp compute temp13 13 temp compute temp14 14 temp compute temp15 15 temp compute temp16 16 temp compute temp17 17 temp compute temp18 18 temp compute temp19 19 temp compute temp20 20 temp compute temp21 21 temp compute temp22 22 temp compute temp23 23 temp compute temph source temp fix 2 source heat 4 0.01 fix 3 sink heat 4 -0.01 thermo 1 variable density equal mass(all)*10/(6.0221417930*vol) thermo_style custom step temp etotal c_temph dump out1 all custom 200 dump.lammpstrj id type x y z timestep 0.25 restart 100000 epon.restart run 1200000 ================================================== 我所選用的勢能函數為PCFF 先做NVT系統平衡到300K 再做NVE 有試著修改過ksapce pppm---> ksapce ewald 但結果還是一樣 也有更改過pair_style lj/class/coul/cut 10.0 8.0 結果一樣 最後做出?淼南到y溫度會上升到700多K ，如果在平衡更久溫度會在更高 不知道各位還有沒有什麼其他的想法，可以改進問題 能讓系統維持在300K左右(熱傳導與系統本身的溫度有關係) 我怕雖然做得出?淼膖emp gradient，但結果不是在300K時材料的熱傳導值 希望大家提出建議交流交流一下!! 謝謝!!!!!!!! |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
Tanya_snail
木虫 (小有名气)
- 应助: 25 (小学生)
- 金币: 3386.2
- 红花: 3
- 帖子: 298
- 在线: 333.3小时
- 虫号: 1656666
- 注册: 2012-03-01
- 性别: MM
- 专业: 金属材料的微观结构
3楼2013-09-02 22:31:53
2楼2013-09-02 21:36:28
hustli2014
新虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 13.2
- 散金: 102
- 红花: 2
- 帖子: 86
- 在线: 134.5小时
- 虫号: 3093527
- 注册: 2014-03-27
- 性别: GG
- 专业: 高分子材料的加工与成型
4楼2016-07-18 16:03:17
5