| 5 | 1/1 | 返回列表 |
| 查看: 3854 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
熱導(dǎo)率計算問題 已有1人參與
|
||
|
我是計算一個石墨烯納米帶300k下的熱導(dǎo)率 下面是我的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ā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +6 | 勇敢的小吳 2026-04-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 372求調(diào)劑 +3 | jj涌77 2026-04-02 | 3/150 |
|
|
[考研] 085801 總分275 本科新能源 求調(diào)劑 +10 | bradoner 2026-04-01 | 12/600 |
|
|
[考研] 0710生物學(xué),325求調(diào)劑 +3 | mkkkkkl 2026-04-01 | 3/150 |
|
|
[考研] 377求調(diào)劑 +3 | RASKIN 2026-04-02 | 3/150 |
|
|
[考研] 一志愿北京科技,085601總分305求調(diào)劑 +9 | 半生瓜! 2026-04-01 | 11/550 |
|
|
[考研] 求調(diào)劑 +4 | DADA怪 2026-03-31 | 4/200 |
|
|
[考研] 348求調(diào)劑 +9 | zzzzyk123 2026-04-01 | 9/450 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +6 | 篤行致遠. 2026-03-28 | 6/300 |
|
|
[考研] 333求調(diào)劑 +4 | 阿科逸 2026-03-31 | 4/200 |
|
|
[考研] 初試301,代碼085701環(huán)境工程,本碩一致,四六級已過,有二區(qū)一作,共發(fā)表5篇論文 +3 | axibli 2026-04-01 | 3/150 |
|
|
[考研] 土木304求調(diào)劑 +5 | 頂級擦擦 2026-03-31 | 5/250 |
|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +4 | 1段時間 2026-03-31 | 5/250 |
|
|
[考研] 已決定調(diào)劑院校 +8 | JKSOIID 2026-03-26 | 8/400 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
[考研]
|
Gymno 2026-03-30 | 6/300 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進合金團隊招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 一志愿雙一流機械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|