| 查看: 618 | 回复: 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左右，差的有点多。 |
回复此楼» 猜你喜欢
每到中夜，情难自抑
已经有58人回复
-
实验室太吵闹，无法安静学习，怎么办？
已经有12人回复
-
申请2024或2025年博士研究生
已经有12人回复
-
在大地上我们只过一生---看完我的阿勒泰上头了好几天，完结那天晚上几乎失眠
已经有13人回复
-
双非本科毕业论文，气人
已经有12人回复
-
化学B02口青基 代表作都是什么水平的？向大佬求助
已经有10人回复
-
化学口B0109(高分子合成)，拿青年基金一般需要怎样的文章水平？
已经有22人回复
-
每天学术时间不能保证，能保证的只有：
已经有10人回复
-
研0二导师分到新来的博士后靠谱吗
已经有7人回复
-
研究生在毕业答辩时挂了，遗憾
已经有13人回复
2楼2022-08-28 17:57:59
10