| 2 | 1/1 | 返回列表 |
| 查看: 1285 | 回復(fù): 1 | |||
| 【懸賞金幣】回答本帖問題,作者CarbonO2將贈(zèng)送您 20 個(gè)金幣 | |||
[求助]
lammps nemd計(jì)算熱導(dǎo)率問題
|
|||
|
大家好,我在使用nemd計(jì)算帶有凸起的石墨烯熱導(dǎo)率時(shí)(冷熱源均取在凸起上),出現(xiàn)凸起部分有較大溫度變化,但是主體部分溫度幾乎不變的情況,請問該如何解決該問題?以下為我的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 |
|
lammps系統(tǒng)學(xué)習(xí)內(nèi)容可以了解一下 https://docs.qq.com/pdf/DWG9YdEloVWxJeWhQ 發(fā)自小木蟲Android客戶端 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿江南大學(xué)085701環(huán)境工程專碩總分287求調(diào)劑 +5 | 18266118446 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿安徽大學(xué)材料工程專碩313分,求調(diào)劑的學(xué)校 +7 | Yu先生 2026-03-10 | 9/450 |
|
|
[考研] 材料371求調(diào)劑 +7 | 鱷魚? 2026-03-11 | 9/450 |
|
|
[考研] 材料工程專碩,一志愿中國礦業(yè)大學(xué),總分314,求調(diào)劑 +4 | 無懈可擊的巨人 2026-03-10 | 4/200 |
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 10/500 |
|
|
[考研] 一志愿天津大學(xué),英一數(shù)二305分求調(diào)劑,四六級(jí)已過 +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
[考研]
|
likeihood 2026-03-06 | 16/800 |
|
|
[考研]
|
簡木ChuFront 2026-03-09 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | 辛仁豆腐 2026-03-08 | 5/250 |
|
|
[考研] 313求調(diào)劑 +4 | Yyt楊1 2026-03-07 | 5/250 |
|
|
[考研] 求調(diào)劑,數(shù)一英一274分 +4 | 小菲會(huì)努力 2026-03-08 | 4/200 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 0817找調(diào)劑 +6 | kk扛 2026-03-08 | 6/300 |
|
|
[考研] 297求調(diào)劑 +3 | 胡達(dá)靈 2026-03-05 | 5/250 |
|
|
[考研] 理學(xué)287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
[考研] 085701環(huán)境工程專業(yè),初試305,均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級(jí)有競賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
[考研] 材料277分求調(diào)劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|