想问下大家，我想要计算气液表面张力，模型为气液气盒子。将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 返回小木虫查看更多