| 查看: 1010 | 回复: 2 | ||
| 【悬赏金币】回答本帖问题，作者我是东成西就将赠送您 50 个金币 | ||
[求助]
lammps计算热导率自带算例求助已有1人参与
|
||
|
请教各位，如下lammps自带算例中的convert是用于计算转换因子scale的，请问convert和scale是如何确定的？非常感谢各位的帮助！ 以下是自带算例的输入脚本： # Sample LAMMPS input script for thermal conductivity of solid Ar units real variable T equal 70 variable V equal vol variable dt equal 4.0 variable p equal 200 # correlation length variable s equal 10 # sample interval variable d equal $p*$s # dump interval # convert from LAMMPS real units to SI variable kB equal 1.3806504e-23 # [J/K] Boltzmann variable kCal2J equal 4186.0/6.02214e23 variable A2m equal 1.0e-10 variable fs2s equal 1.0e-15 variable convert equal ${kCal2J}*${kCal2J}/${fs2s}/${A2m} # setup problem dimension 3 boundary p p p lattice fcc 5.376 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region box block 0 4 0 4 0 4 create_box 1 box create_atoms 1 box mass 1 39.948 pair_style lj/cut 13.0 pair_coeff * * 0.2381 3.405 timestep ${dt} thermo $d # equilibration and thermalization velocity all create $T 102486 mom yes rot yes dist gaussian fix NVT all nvt temp $T $T 10 drag 0.2 run 8000 # thermal conductivity calculation, switch to NVE if desired #unfix NVT #fix NVE all nve reset_timestep 0 compute myKE all ke/atom compute myPE all pe/atom compute myStress all stress/atom NULL virial compute flux all 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 all ave/correlate $s $p $d & c_flux[1] c_flux[2] c_flux[3] type auto file J0Jt.dat ave running variable scale equal ${convert}/${kB}/$T/$T/$V*$s*${dt} variable k11 equal trap(f_JJ[3])*${scale} variable k22 equal trap(f_JJ[4])*${scale} variable k33 equal trap(f_JJ[5])*${scale} thermo_style custom step temp v_Jx v_Jy v_Jz v_k11 v_k22 v_k33 run 100000 variable k equal (v_k11+v_k22+v_k33)/3.0 variable ndens equal count(all)/vol print "average conductivity: $k[W/mK] @ $T K, ${ndens} /A^3"@smutao@oxox6085 |
回复此楼» 猜你喜欢
Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？
已经有10人回复
-
0854－0855调剂
已经有10人回复
-
南京大学-广州大学联合招聘博士后 欢迎广大优秀人才！！！
已经有11人回复
-
召集青年编委
已经有6人回复
-
申博求助
已经有8人回复
-
Nature一直在编辑手里，考虑好几天了，是悬了吗
已经有7人回复
-
博士生招生！！
已经有5人回复
-
国社科项目，你们学校都限额申报吗？
已经有8人回复
-
企业博后是否能申请CSC博士后项目？谢谢
已经有6人回复
-
博士白读了
已经有40人回复
2楼2019-02-17 12:13:38
S-Loser
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1126.4
- 散金: 218
- 红花: 3
- 沙发: 1
- 帖子: 613
- 在线: 231.6小时
- 虫号: 7545220
- 注册: 2017-11-29
- 专业: 特种医学(航空、航天、航海
3楼2019-02-22 15:45:30
66