 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 864 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
wuqiqilwj新蟲(chóng) (小有名氣)
|
[求助]
文本中倒數(shù)第二行的熱導(dǎo)率公式怎么來(lái)的？
|
||
|
# 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é)模擬 |
» 猜你喜歡
268求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑
已經(jīng)有4人回復(fù)
-
302求調(diào)劑
已經(jīng)有11人回復(fù)
-
306求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué)
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
313求調(diào)劑
已經(jīng)有3人回復(fù)
-
南昌大學(xué)材料專碩311分求調(diào)劑
已經(jīng)有6人回復(fù)
-
316求調(diào)劑
已經(jīng)有6人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有7人回復(fù)
bojue327
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 3
- 在線: 1.2小時(shí)
- 蟲(chóng)號(hào): 4993394
- 注冊(cè): 2016-09-08
3樓2016-09-21 21:03:24
wuqiqilwj
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 51.6
- 帖子: 53
- 在線: 16.3小時(shí)
- 蟲(chóng)號(hào): 4249120
- 注冊(cè): 2015-11-26
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2016-09-18 11:12:15
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 280求調(diào)劑 +7 | 咕嚕曉曉 2026-03-18 | 8/400 |
|
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +8 | 吃吃吃才有意義 2026-03-19 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無(wú)極 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|