|
LAMMPS讲解38-SPC/E模型下水的表面张力计算本例发布平直液面的表面张力计算过程。使用的是水的SPCE模型。方法是计算垂直液面和平行液面的压强分布，即pxx，pyy和pzz在垂直液面方向上的分布曲线。如需原始的in和data文件以及输出的p.profile和计算用的Excel文件，请加微信lmp_zhushou。 计算命令是 compute myp all stress/atom NULL compute cc1 all chunk/atom bin/1d x lower 1.0 units box fix 1 all ave/chunk 1 1000000 1000000 cc1 c_myp[1] c_myp[2] c_myp[3] norm sample file p.profile 计算完成后打开输出的p.profile文件，文件中第一列是bin的编号，第二列是bin的位置，第三列是bin中的平均原子数，第四列到第六列是c_myp[1]，c_myp[2] ，c_myp[3]。用第四列乘以第三列再除以bin的体积取相反数就是pxx，用第五列乘以第三列除以bin的体积取相反数就是pyy，用第六列乘以第三列除以bin的体积取相反数就是pzz。将得到的pyy和pzz取平均就得到了平行于液面的压强pt，而pxx就是垂直于液面的压强pn。将pn列减去pt列，并将得到的结果乘以bin的厚度，此处是1埃。将最终得到的结果列求和除以2就是表面张力值。本例计算出的表面张力值是54mN/m，文献值是52mN/m，实验值是72mN/m。由此可见本例的计算结果是可靠的。 pxx，pyy和pzz沿垂直液面方向的分布曲线如下图
感谢鲍路瑶老师的分享，内容来自于鲍老师分享出来的资料 如有需要添加微信：lmp_zhushou 进入微信群，帮助他人，共建社区 获取完整版lammps讲义可以加微信lmp_zhushou或加入QQ群994359511 |
