qujinfeng

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 608.9
散金: 105
紅花: 1
帖子: 130
在線: 105.9小時(shí)
蟲號(hào): 1396671
注冊(cè): 2011-09-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 熱導(dǎo)率計(jì)算問(wèn)題已有1人參與

我是計(jì)算一個(gè)石墨烯納米帶300k下的熱導(dǎo)率下面是我的in文件不知道為什么體系就是達(dá)不到平衡  整體溫度一直在285到315之間震蕩局部溫度震蕩范圍更大請(qǐng)問(wèn)我該怎么做是什么問(wèn)題怎么樣才能讓體系弛豫達(dá)到平衡狀態(tài)呢  求指教！  附件是計(jì)算的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ù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : log.lammps

2012-06-28 20:29:51, 30.64 K

» 猜你喜歡

一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有13人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有6人回復(fù)
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷已經(jīng)有10人回復(fù)
298求調(diào)劑已經(jīng)有6人回復(fù)
310求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有6人回復(fù)
307求調(diào)劑已經(jīng)有11人回復(fù)
材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑已經(jīng)有5人回復(fù)
315分，誠(chéng)求調(diào)劑，材料與化工085600 已經(jīng)有3人回復(fù)
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) 已經(jīng)有10人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

lammps求助：計(jì)算含有不同原子體系的熱導(dǎo)率的in文件該怎么修改已經(jīng)有7人回復(fù)
鋼包用耐火材料壁厚的設(shè)計(jì)計(jì)算已經(jīng)有9人回復(fù)
請(qǐng)教關(guān)于多元（三種不同成分以上）復(fù)合材料熱導(dǎo)率和熱膨脹系數(shù)的問(wèn)題已經(jīng)有5人回復(fù)
SiC的熱導(dǎo)率問(wèn)題已經(jīng)有9人回復(fù)
復(fù)合材料比熱容的確定已經(jīng)有16人回復(fù)
【求助】求解關(guān)于測(cè)熱導(dǎo)率試樣尺寸問(wèn)題已經(jīng)有3人回復(fù)
【求助】第一性原理計(jì)算材料導(dǎo)電、導(dǎo)熱、塞貝克系數(shù)問(wèn)題已經(jīng)有18人回復(fù)
【資源】熱物性的瞬態(tài)測(cè)量已經(jīng)有5人回復(fù)

1樓 2012-06-28 20:30:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

維特的日記

木蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 3814.9
紅花: 2
帖子: 169
在線: 55.5小時(shí)
蟲號(hào): 3159711
注冊(cè): 2014-04-24
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

【答案】應(yīng)助回帖

我感覺(jué)是你nvt系綜的damp值設(shè)置有問(wèn)題
另外你計(jì)算熱導(dǎo)率公式里的36.9*3.4*10e-20是那個(gè)區(qū)域的面積有問(wèn)題把

贊一下

回復(fù)此樓

2樓2014-08-20 22:23:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

牛犢小子

鐵桿木蟲 (著名寫手)

應(yīng)助: 14 (小學(xué)生)
金幣: 8786.5
散金: 275
紅花: 1
帖子: 1165
在線: 169.8小時(shí)
蟲號(hào): 2659995
注冊(cè): 2013-09-16
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

幫頂…………

回復(fù)此樓

3樓2014-08-22 10:46:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

牛犢小子

鐵桿木蟲 (著名寫手)

應(yīng)助: 14 (小學(xué)生)
金幣: 8786.5
散金: 275
紅花: 1
帖子: 1165
在線: 169.8小時(shí)
蟲號(hào): 2659995
注冊(cè): 2013-09-16
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

那個(gè)40.448是什么呀

贊一下

回復(fù)此樓

4樓2014-08-27 09:35:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tang_wei

金蟲 (正式寫手)

應(yīng)助: 24 (小學(xué)生)
金幣: 1353.5
紅花: 28
帖子: 428
在線: 81.7小時(shí)
蟲號(hào): 4067192
注冊(cè): 2015-09-12
專業(yè): 高分子物理與高分子物理化

你好，請(qǐng)問(wèn)一下你的問(wèn)題解決了嗎？我也是用fix heat 方法，發(fā)現(xiàn)即使我所加的熱流為0的話，體系每一層的溫度也不是想的的，波動(dòng)較大

贊一下(1人)

回復(fù)此樓

5樓2015-11-20 09:24:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

青春逍遙啊

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 414.6
紅花: 1
帖子: 202
在線: 28.9小時(shí)
蟲號(hào): 2465810
注冊(cè): 2013-05-15
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

樓主您好，很抱歉我無(wú)法回答您的問(wèn)題，我是剛剛?cè)腴T做熱導(dǎo)率的，我想問(wèn)您variable temp atom c_ke*1.6e+4/(1.5*1.38)是什么意思，我看手冊(cè)沒(méi)弄明白，如果您有時(shí)間，麻煩您可以教我。謝謝。

贊一下

回復(fù)此樓

6樓2015-12-02 19:57:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

計(jì)算化學(xué)

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.8
帖子: 13
在線: 46.1小時(shí)
蟲號(hào): 8764310
注冊(cè): 2018-05-15
性別: MM

引用回帖:

5樓: Originally posted by tang_wei at 2015-11-20 09:24:56
你好，請(qǐng)問(wèn)一下你的問(wèn)題解決了嗎？我也是用fix heat 方法，發(fā)現(xiàn)即使我所加的熱流為0的話，體系每一層的溫度也不是想的的，波動(dòng)較大

我的也是，就是兩端設(shè)置了冷熱區(qū)能看到，但中間部分就像一條直線，且波動(dòng)特大，請(qǐng)教一下，你知道原因嗎

贊一下

回復(fù)此樓

7樓2018-05-24 09:38:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

青春逍遙啊

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 414.6
紅花: 1
帖子: 202
在線: 28.9小時(shí)
蟲號(hào): 2465810
注冊(cè): 2013-05-15
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

7樓: Originally posted by 計(jì)算化學(xué) at 2018-05-24 09:38:47
我的也是，就是兩端設(shè)置了冷熱區(qū)能看到，但中間部分就像一條直線，且波動(dòng)特大，請(qǐng)教一下，你知道原因嗎...

可以用其他的方法計(jì)算，比如mp之類的呢

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

8樓2018-06-12 18:46:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

機(jī)智康

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46.2
帖子: 24
在線: 9.3小時(shí)
蟲號(hào): 8691806
注冊(cè): 2018-05-07
專業(yè): 固體力學(xué)

問(wèn)題解決了嗎？

發(fā)自小木蟲Android客戶端

回復(fù)此樓

9樓2020-05-07 13:35:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +6	幸運(yùn)哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 306求調(diào)劑 +6	chuanzhu川燭 2026-03-18	6/300	2026-03-22 19:39 by 公瑾逍遙
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-20	3/150	2026-03-22 16:00 by ColorlessPI
[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 生物學(xué)調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	4/200	2026-03-22 11:04 by 搏擊518
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 268求調(diào)劑 +9	簡(jiǎn)單點(diǎn)0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂(lè)29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao