 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 860 | 回復(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ù)此樓» 收錄本帖的淘帖專輯推薦
 分子動力學(xué)模擬 |
» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有23人回復(fù)
-
265求調(diào)劑
已經(jīng)有9人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有10人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有12人回復(fù)
-
296求調(diào)劑
已經(jīng)有8人回復(fù)
-
招收博士1-2人
已經(jīng)有3人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有11人回復(fù)
-
274求調(diào)劑
已經(jīng)有8人回復(fù)
-
290求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑
已經(jīng)有4人回復(fù)
wuqiqilwj
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 51.6
- 帖子: 53
- 在線: 16.3小時
- 蟲號: 4249120
- 注冊: 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小時
- 蟲號: 4249120
- 注冊: 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ā)表 | |
|---|---|---|---|---|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +6 | 吃吃吃才有意義 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-20 | 4/200 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +7 | sllhht 2026-03-19 | 8/400 |
|
|
[考研] 304求調(diào)劑 +5 | 曼殊2266 2026-03-18 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +6 | 梁梁校校 2026-03-17 | 6/300 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +17 | enenenhui 2026-03-13 | 18/900 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|