| 5 | 1/1 | 返回列表 |
| 查看: 3834 | 回復: 8 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
熱導率計算問題 已有1人參與
|
||
|
我是計算一個石墨烯納米帶300k下的熱導率 下面是我的in文件 不知道為什么體系就是達不到平衡 整體溫度一直在285到315之間震蕩 局部溫度震蕩范圍更大 請問我該怎么做 是什么問題 怎么樣才能讓體系弛豫達到平衡狀態(tài)呢 求指教! 附件是計算的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 |
木蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
金蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調劑 +3 | α____ 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 328求調劑,英語六級551,有科研經歷 +6 | 生物工程調劑 2026-03-17 | 10/500 |
|
|
[考研] 298求調劑 +6 | 上岸6666@ 2026-03-20 | 6/300 |
|
|
[考研] 環(huán)境學碩288求調劑 +6 | 皮皮皮123456 2026-03-22 | 6/300 |
|
|
[考研] 廣西大學材料導師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 求調劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 0805 316求調劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 297求調劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 一志愿深大,0703化學,總分302,求調劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 求調劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 求調劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 350求調劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 材料與化工專碩調劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 085601材料工程專碩求調劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 材料專碩306英一數二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 085601求調劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|