| 查看: 609 | 回复: 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左右，差的有点多。 |
回复此楼» 猜你喜欢
换导师
已经有8人回复
-
来用亲身体会，一起说说年轻老师的辛苦
已经有21人回复
-
求助在硝基苯上进行F-C反应甲基化的方案
已经有8人回复
-
博三一直没文章怎么办
已经有22人回复
-
推荐转投（ transfer pending）是否有用？
已经有3人回复
-
白天不懂夜的黑，90后青椒与60后老板
已经有4人回复
-
求期刊推荐
已经有5人回复
-
基金委也是打工的
已经有21人回复
-
博士去高校就是为了有寒暑假吗？
已经有22人回复
-
这个模块怎么成了烧香拜佛的地方了
已经有5人回复
2楼2022-08-28 17:57:59
50