 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 7500 | 回復(fù): 13 | ||||
zxzj05榮譽版主 (著名寫手)
|
[交流]
【轉(zhuǎn)帖】熱導率計算的in文件 已有12人參與
|
|||
|
看到有不少人在找熱導率計算方面的in文件，我就貢獻三個in文件吧，僅供參考。 同時，附件里貼出了我的計算結(jié)果。EMD的輸出結(jié)果（compute heat/flux command+compute tc command的計算結(jié)果）中， “ac.dat”（見附件中的"ac.wmf"）是熱流自相關(guān)函數(shù)（我已經(jīng)修改了compute_tc.cpp，目前輸出的是normalized HCACF，但結(jié)果中給出的還是沒有歸一化的熱流自相關(guān)函數(shù)，但形狀和歸一化的是一樣的，請注意�。╇Sm的變化，"tc.dat"(見附件中的"tc.wmf"  是熱導率隨m的變化（m的涵義請參看熱導率計算的Green Kubo離散化公式,見附件"Comparison of atomic-level simulation methods for computing thermal conductivity”中的(9)式），"tc_time.dat"(見附件中的"tc_time.wmf"是熱導率隨時間的變化；NEMD的結(jié)果中，"temp.profile"(見附件中的"temp_distribution.wmf"是z方向上的溫度分布，"thermal_conductivity.dat“(見附件中的"NEMD.wmf"表示熱導率隨時間的變化，兩者都包含了fix heat command 和fix thermal/conductivity command的計算結(jié)果。1.用compute heat/flux command+compute tc command得到熱流自相關(guān)函數(shù)和熱導率（EMD方法） # 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 mass 1 1.0 velocity all create 0.71 458127641 mom yes rot yes dist gaussian units box # LJ 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.000466 fix nve all nve thermo_style custom step temp etotal vol thermo_modify lost warn thermo 100 # Run timestep 0.000466 run 200000 reset_timestep 0 # -------------- Flux calculation in nve --------------- compute myKE all ke/atom compute myPE all pe/atom compute myStress all stress/atom virial variable factor_ac equal 1.0 variable factor_tc equal 1.3806504e-23*sqrt(1.67e-21/6.633e-26)/3.405e-10^2 compute jflux all heat/flux myKE myPE myStress compute tc all tc c_thermo_temp c_jflux v_factor_ac v_factor_tc iso first 10000 900000 100000 fix tc_out all ave/time 1 1 1 c_tc file tc_time.dat thermo_style custom step temp restart 100000 restart.* run 1000000 2. 用fix thermal/conductivity command得到溫度梯度，進而得到熱導率(NEMD方法) # 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 mass 1 1.0 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 heat_swap all thermal/conductivity 10 z 32 fix e_exchange all ave/time 10 10000 100000 f_heat_swap file e_exchange.dat variable thermal_conductivity equal f_e_exchange/(0.000466*10.0*4.0*f_delta_out)*1.3806504e-23/3.405e-10/3.405e-10*sqrt(1.67e-21/6.633e-26)*6.0/8.0 # 以上variable命令需要特別注意，因為我所模擬的系統(tǒng)，盒子邊長Lx=Ly=Lz，熱導率計算公式經(jīng)過推導變成為e_exchange/(4.0*t*L*delta_T)， # 為了不在in文件里給L賦值，我修改了fix_thermal_conductivity.cpp文件(見附件)，將e_exchange修改成了 e_exchange += force->mvv2e * (all[0].value - all[1].value) / (domain->zprd); 同時在end_of_step() 里添加了一句 “ e_exchange = 0.0;“,詳見附件中的fix_thermal_conductivity.cpp，這樣所得的e_exchange曲線基本上是一條水平曲線，而不是用原來的fix thermal/conductivity command所得到的斜向上的曲線，請注意�。�！ # 所以才出現(xiàn)以上variable的表達式。 # 請看明白后再做計算，免得算出錯誤的結(jié)果�。�！ fix thermal_conductivity_out all ave/time 100000 1 100000 v_thermal_conductivity file thermal_conductivity.dat # Run run 10000000 3. 用fix heat command建立溫度梯度，進而得到熱導率（NEMD方法） # 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 【轉(zhuǎn)自】http://www.mdbbs.org/thread-26451-1-1.html |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 收藏 |
» 猜你喜歡
能源材料化學課題組招收碩士研究生8-10名
已經(jīng)有6人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
268求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料，紡織，生物（0856、0710），化學招生啦
已經(jīng)有7人回復(fù)
-
301求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有7人回復(fù)
-
材料與化工求調(diào)劑
已經(jīng)有6人回復(fù)
-
299求調(diào)劑
已經(jīng)有4人回復(fù)
-
277調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- DFT Scan之后使用MP2計算能量，如何建立批量輸入文件？
已經(jīng)有6人回復(fù)
- 如何通過聲子計算求解熱導率
已經(jīng)有5人回復(fù)
- 大家來幫我看看這兩個INCAR 文件有什么區(qū)別哦。。為啥計算結(jié)果不一樣呢
已經(jīng)有4人回復(fù)
- 【求助】關(guān)于VASP計算材料的晶格熱導率
已經(jīng)有8人回復(fù)
- 關(guān)于計算完后 輸出文件 _BandStr.castep文件里內(nèi)容的疑問
已經(jīng)有4人回復(fù)
- SiC的熱導率問題
已經(jīng)有9人回復(fù)
- 復(fù)合材料熱導率的理論計算方法及軟件
已經(jīng)有5人回復(fù)
- 【討論】lammps 用sw勢計算硅的熔點怎么這么高（2400k),附in文件
已經(jīng)有13人回復(fù)
- 【求助】關(guān)于QHA模塊的使用
已經(jīng)有17人回復(fù)
babaleo
木蟲 (正式寫手)
- 應(yīng)助: 18 (小學生)
- 金幣: 2391.1
- 散金: 40
- 紅花: 1
- 帖子: 375
- 在線: 477.1小時
- 蟲號: 981394
- 注冊: 2010-03-25
- 專業(yè): 人工晶體
5樓2011-01-17 10:12:19
c_friends
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1012.5
- 散金: 26
- 帖子: 251
- 在線: 134.7小時
- 蟲號: 1156964
- 注冊: 2010-11-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2010-12-12 18:35:59
zxzj05
榮譽版主 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.242
- 金幣: 5912.5
- 紅花: 2
- 帖子: 1528
- 在線: 144.5小時
- 蟲號: 406147
- 注冊: 2007-06-18
- 性別: GG
- 專業(yè): 金屬功能材料
- 管轄: 第一性原理
3樓2010-12-13 15:16:07
354895492
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5835.2
- 散金: 50
- 紅花: 4
- 帖子: 375
- 在線: 221.5小時
- 蟲號: 1165513
- 注冊: 2010-12-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
，我的分子模擬論壇也沒金幣了 4樓2010-12-20 22:19:54
WJ-MS
新蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1559.2
- 散金: 100
- 帖子: 194
- 在線: 92.5小時
- 蟲號: 1697144
- 注冊: 2012-03-17
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
6樓2013-05-27 09:25:48
某某人
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 38
- 帖子: 9
- 在線: 5.8小時
- 蟲號: 2396451
- 注冊: 2013-04-02
- 專業(yè): 熱物性與熱物理測試技術(shù)
7樓2013-05-27 14:09:47
WJ-MS
新蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1559.2
- 散金: 100
- 帖子: 194
- 在線: 92.5小時
- 蟲號: 1697144
- 注冊: 2012-03-17
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
8樓2013-05-28 08:32:49
燦爛的幸福
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 906.3
- 帖子: 93
- 在線: 101.5小時
- 蟲號: 1335565
- 注冊: 2011-07-01
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
9樓2013-08-08 10:41:20
xiaohongxing
金蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 3948.3
- 帖子: 994
- 在線: 350.4小時
- 蟲號: 1484719
- 注冊: 2011-11-09
- 性別: GG
- 專業(yè): 核技術(shù)及其應(yīng)用
10樓2014-12-10 22:02:48
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zxzj05 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進步 2026-03-11 | 5/250 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 341求調(diào)劑 +3 | 番茄頭--- 2026-03-10 | 3/150 |
|
|
|
[考研] 一志愿中科院，化學方向，295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
|
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 一志愿211化學學碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123！ 2026-03-11 | 4/200 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|