Notice: Undefined index: portal_style in /data/wwwroot/emuch.net/t.php on line 5
多体耗散粒子动力学(MDPD)液体表面张力的LAMMPS模拟计算 - 分子模拟 - 小木虫 - 学术 科研 互动社区
当前位置: 首页 > 分子模拟 >多体耗散粒子动力学(MDPD)液体表面张力的LAMMPS模拟计算

多体耗散粒子动力学(MDPD)液体表面张力的LAMMPS模拟计算

作者 luyaobao
来源: 小木虫 900 18 举报帖子
+关注

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

今日热帖
  • 精华评论
猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓