| 查看: 624 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者yingqian50将赠送您 20 个金币 | ||
yingqian50金虫 (初入文坛)
|
[求助]
计算热导率，得到的结果不理想，大了几倍
|
|
|
我是要计算烷烃的热导率，照着lammps里面计算G-K的例子修改后计算的，试着修改了关联长度，虽然得到的结果会有变化，但是还是大了实际几倍，希望能够指点一下，感谢！ 附上in文件和输出文件 variable t equal 303 variable p equal 1000 # correlation length variable s equal 8 # sample interval variable d equal $p*$s # dump interval variable kB equal 1.3806504e-23 # [J/K] Boltzmann variable kCal2J equal 6.948e-21 variable A2m equal 1.0e-10 variable fs2s equal 1.0e-15 variable convert equal ${kCal2J}*${kCal2J}/${fs2s}/${A2m} # setup problem units real atom_style full boundary p p p read_data 1.data region 1 block INF INF INF INF INF INF units box group flow region 1 #----- initialization ----- minimize 1.0e-6 1.0e-8 1000 1000 velocity flow create 303 9567969 #-----relax process----- # 1st equilibration run fix 1 flow npt temp 303 303 100 iso 0 0 1000 thermo_style custom step temp pe etotal lz thermo 10000 run 100000 velocity flow scale $t # thermal conductivity calculation unfix 1 fix 1 flow npt temp 303 303 100 iso 1 1 1000 thermo_style custom step temp pe etotal lz thermo 10000 run 200000 velocity flow scale $t reset_timestep 0 compute myKE flow ke/atom compute myPE flow pe/atom compute myStress flow stress/atom NULL virial compute flux flow heat/flux myKE myPE myStress variable Jx equal c_flux[1]/vol variable Jy equal c_flux[2]/vol variable Jz equal c_flux[3]/vol fix JJ flow ave/correlate $s $p $d & c_flux[1] c_flux[2] c_flux[3] type auto & file profile.heatflux ave running variable scale equal ${convert}/${kB}*$s*dt/$t/$t/vol variable k11 equal trap(f_JJ[3])*${scale} variable k22 equal trap(f_JJ[4])*${scale} variable k33 equal trap(f_JJ[5])*${scale} variable k44 equal (v_k11+v_k22+v_k33)/3 thermo $d thermo_style custom step temp v_Jx v_Jy v_Jz v_k11 v_k22 v_k33 v_k44 run 1400000 variable kappa equal (v_k11+v_k22+v_k33)/3.0 print "running average conductivity: ${kappa}" 输出文件结果稳定在0.38附近，实际应该是0.13左右，差的有点多。 |
回复此楼» 猜你喜欢
请各位老师谈谈身边非升即走的例子
已经有12人回复
-
中者为王-观某校级项目为领导瓜分有感
已经有25人回复
-
美国对华为的制裁会影响英文论文发表吗
已经有7人回复
-
函评意见：既然你有中科院****人才项目了，就不要申请面上了。
已经有32人回复
-
难上的课安排给我，杂事安排给我。评优时我就没份儿了，即便我做得不比绝大多数人差
已经有12人回复
-
前些天开会有个人见到人就搞关系，一查此人全是MDPI/Hindawi论文，鄙视!
已经有51人回复
-
B01化学口的，近几年要什么样的评价才能上会呀？
已经有5人回复
-
Analytica Chimica Acta 投稿系统中 abstract 字符数限制600什么意思？
已经有5人回复
-
博士找工作真难
已经有17人回复
-
审稿交流
已经有5人回复
2楼2022-08-28 17:57:59
5