Notice: Undefined index: portal_style in /data/wwwroot/emuch.net/t.php on line 5
lammps中气液表面张力计算,得不到稳定的结果 - 分子模拟 - 小木虫 - 学术 科研 互动社区
当前位置: 首页 > 分子模拟 >lammps中气液表面张力计算,得不到稳定的结果

lammps中气液表面张力计算,得不到稳定的结果

作者 dbsydxdm
来源: 小木虫 200 4 举报帖子
+关注

想问下大家,我想要计算气液表面张力,模型为气液气盒子。将z方向分块,使用ave/chunk命令输出每块法向力和切向力的差值rk最后加和得到表面张力,但是在气相和液相区域的切向力和法向力波动比较大,特别是气相。该怎么得到稳定的表面张力值?还有就是得到每块的stress/atom需要乘以每块原子数得到应力值吗,还是直接输出的每块stress/atom就是应力值。

compute p all stress/atom tmid
variable pn atom -c_p[3]/(20*20)
variable pt atom -0.5*(c_p[1]+c_p[2])/(20*20)
variable rk atom (v_pn-v_pt)

compute cc1 zall chunk/atom bin/1d z center 0.25 units box
fix st zall ave/chunk 10 1000 10000 cc1 v_rk v_pn v_pt file st.dat 返回小木虫查看更多

今日热帖
  • 精华评论
  • dbsydxdm

    有没有做表面张力的 可以交流下

  • santiago℃

    引用回帖:
    2楼: Originally posted by dbsydxdm at 2020-12-11 15:42:11
    有没有做表面张力的 可以交流下

    我也在做表面张力,不过我是刚安装上软件的那种程度

  • dbsydxdm

    引用回帖:
    3楼: Originally posted by santiago℃ at 2021-01-04 13:09:04
    我也在做表面张力,不过我是刚安装上软件的那种程度...

    大概就是这么算的用上面的程序

  • mpiyq

    引用回帖:
    2楼: Originally posted by dbsydxdm at 2020-12-11 15:42:11
    有没有做表面张力的 可以交流下

    有问题可以交流下  2289968992

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓