该in文件采用LAMMPS计算MDPD液体的表面张力,计算公式为
γ=1/2 ∫_(y_1)^(y_2)(p_yy-1/2(p_xx+p_zz)dy

Sample Text

y_1和y_2分别为模拟盒子在y方向的上下边界位置，如下图所示。公式中第一个1⁄2是因为积分跨越了两个界面，第二个1⁄2是pxx和pzz理论上是相等的所以相加除以2以提高精度。
通过修改in文件中的参数可直接进行模拟得到用以发表论文的结果，如不同MDPD参数下的表面张力值

使用方法将文件中代码拷贝粘贴入文本文件中，并将文本文件命名为in.tension（包括扩展名）
模拟完成将输出press.profile文件，该文件包含pxx，pyy和pzz沿y轴的分布，将文件中的数据粘贴入Excel中，pxx=-1*第四列*第三列/0.01/64/64，pyy=-1*第五列*第三列/0.01/64/64，pzz=-1*第六列*第三列/0.01/64/64。将pxx，pyy和pzz带入上述公式即可得到表面张力。

文件地址：https://www.jianguoyun.com/p/DaiMivYQ38aMCRjOhNYD 返回小木虫查看更多