 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 861 | 回復(fù): 3 | |||
wuqiqilwj新蟲 (小有名氣)
|
[求助]
文本中倒數(shù)第二行的熱導(dǎo)率公式怎么來的？
|
|
# MD simulation of Ar thermal conductivity # Initialization units lj dimension 3 newton on boundary p p p atom_style atomic neighbor 0.3 bin neigh_modify check yes lattice fcc 0.844 region box block -4 4 -4 4 -4 4 units lattice create_box 1 box create_atoms 1 box region up1 block INF INF INF INF -0.5 -0.25 units lattice region up2 block INF INF INF INF 0.5 0.75 units lattice region up union 2 up1 up2 region down1 block INF INF INF INF -3.5 -3.25 units lattice region down2 block INF INF INF INF 3.5 3.75 units lattice region down union 2 down1 down2 region hot block INF INF INF INF 0.0 0.25 units lattice group hot region hot region cold block INF INF INF INF -4.0 -3.75 units lattice group cold region cold mass 1 1.0 #mass0 6.633e-26 #epsilon0 1.67e-21 #sigma0 3.405e-10 velocity all create 0.71 458127641 mom yes rot yes dist gaussian units box # Tersoff potential ********************************************************* pair_style lj/cut 2.8 pair_coeff 1 1 1.0 1.0 # LJ parameters for Ar-Ar fix temp all temp/berendsen 0.71 0.71 0.0466 fix nve all nve compute ke all ke/atom variable temp atom c_ke/(1.5*1.0) fix temp_profile all ave/spatial 1 100000 100000 z lower 0.25 v_temp file temp.profile units lattice compute up_temp all temp/region up compute down_temp all temp/region down variable delta_temp equal c_up_temp-c_down_temp fix delta_out all ave/time 1 100000 100000 v_delta_temp file delta_temp.dat thermo_style custom step temp etotal vol thermo_modify lost warn thermo 100 # Run timestep 0.000466 run 100001 unfix temp fix hot all heat 1 50 region hot fix cold all heat 1 -50 region cold variable thermal_conductivity equal 50.0*0.5*1.67e-21/3.405e-10/sqrt(6.633e-26/1.67e-21)/((4.0*8.0*8.0*8.0/0.844)^(1.0/3.0)*3.405e-10*2.0*f_delta_out*1.67e-21/1.3806504e-23)*6.0/8.0 fix thermal_conductivity_out all ave/time 100000 1 100000 v_thermal_conductivity file thermal_conductivity.dat # Run run 10000000 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 分子動(dòng)力學(xué)模擬 |
» 猜你喜歡
22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù)
已經(jīng)有4人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有7人回復(fù)
-
南京大學(xué)化學(xué)376求調(diào)劑
已經(jīng)有6人回復(fù)
-
南昌大學(xué)材料專碩311分求調(diào)劑
已經(jīng)有5人回復(fù)
-
295求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿南昌大學(xué)，327分，材料與化工085600
已經(jīng)有9人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有15人回復(fù)
-
308求調(diào)劑
已經(jīng)有3人回復(fù)
-
321求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有8人回復(fù)
wuqiqilwj
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 51.6
- 帖子: 53
- 在線: 16.3小時(shí)
- 蟲號(hào): 4249120
- 注冊(cè): 2015-11-26
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2016-09-18 11:12:15
3樓2016-09-21 21:03:24
wuqiqilwj
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 51.6
- 帖子: 53
- 在線: 16.3小時(shí)
- 蟲號(hào): 4249120
- 注冊(cè): 2015-11-26
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
4樓2016-09-22 19:44:14
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wuqiqilwj 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
|
[考研] 316求調(diào)劑 +5 | 梁茜雯 2026-03-19 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +19 | 不知道叫什么！ 2026-03-15 | 21/1050 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|