想問下大家，我想要計(jì)算氣液表面張力，模型為氣液氣盒子。將z方向分塊，使用ave/chunk命令輸出每塊法向力和切向力的差值rk最后加和得到表面張力，但是在氣相和液相區(qū)域的切向力和法向力波動(dòng)比較大，特別是氣相。該怎么得到穩(wěn)定的表面張力值？還有就是得到每塊的stress/atom需要乘以每塊原子數(shù)得到應(yīng)力值嗎，還是直接輸出的每塊stress/atom就是應(yīng)力值。

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 返回小木蟲查看更多