該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文件中的參數(shù)可直接進(jìn)行模擬得到用以發(fā)表論文的結(jié)果，如不同MDPD參數(shù)下的表面張力值

使用方法將文件中代碼拷貝粘貼入文本文件中，并將文本文件命名為in.tension（包括擴展名）
模擬完成將輸出press.profile文件，該文件包含pxx，pyy和pzz沿y軸的分布，將文件中的數(shù)據(jù)粘貼入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 返回小木蟲查看更多