[Lammps] 用lammps算热传导系数如何算virial部分的贡献？

[复制链接 Copy URL]

xunxun87

1 帖子	0 威望	14 eV

积分: 15

Level 1 能力者

发消息 Send PM

电梯直达 Go to

楼主

发表于 Post on 2016-3-27 00:03:45 | 只看该作者 Only view this author 回帖奖励

|倒序浏览 Reverse view |阅读模式 Reading model

rt. 本人想把kappa matrix分解成KE,PE 和virial部分，已搞定了标量的KE和PE部分，但由于stress 是6个element的vector，不太清楚如何算virial的heat flux vector S*v.

compute myStress all stress/atom virial
variable mySvx atom c_myStress*vx
variable mySvy atom c_myStress*vy
variable mySvz atom c_myStress*vz
compute mySvx all reduce sum v_mySvx
compute mySvy all reduce sum v_mySvy
compute mySvz all reduce sum v_mySvz
fix JJS all ave/correlate $s $p $d c_mySvx c_mySvy c_mySvz type auto file J0Jt.dat ave running
这样应该有问题吧？那么stress的维度该如何修改呢？不胜感激！

收藏 Add to favorites1

回复 Reply

举报 Report

返回列表 Return

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register