| 查看: 623 | 回复: 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左右，差的有点多。 |
回复此楼» 猜你喜欢
迷茫的毕业季
已经有3人回复
-
这样的课题组值得去吗
已经有5人回复
-
一般很空的审稿意见大家一般是怎么回复的
已经有6人回复
-
蹲一个6月7号江苏省青基公示
已经有25人回复
-
间接经费用于绩效支持比例人员等问题
已经有4人回复
-
席夫碱反应ph敏感 控制ph在多少呢？
已经有5人回复
-
难上的课安排给我，杂事安排给我。评优时我就没份儿了，即便我做得不比绝大多数人差
已经有11人回复
-
看到大家一直在说代表作NC，不明觉厉
已经有22人回复
-
前些天开会有个人见到人就搞关系，一查此人全是MDPI/Hindawi论文，鄙视!
已经有49人回复
-
副教授跳槽
已经有16人回复
2楼2022-08-28 17:57:59
5