24小时热门版块排行榜

>论坛更新日志 (3768)
>休闲灌水 (372)
>虫友互识 (344)
>导师招生 (210)
>博后之家 (157)
>招聘信息布告栏 (108)
>文献求助 (89)
>论文投稿 (77)
>硕博家园 (49)
>育儿交流 (32)
>考博 (32)
>基金申请 (31)
>教师之家 (26)
>学术会议 (25)
>考研 (17)
>有机交流 (16)

返回列表

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 luyaobao 的 64 个金币，回帖就立即获得 2 个金币，每人有 1 次机会

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

luyaobao

木虫 (著名写手)

应助: 1 (幼儿园)
金币: 2457.9
帖子: 1533
在线: 150.7小时
虫号: 7847640

[交流] 多体耗散粒子动力学（MDPD）液体表面张力的LAMMPS模拟计算

该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

回复此楼