24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2771  |  回复: 20
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 luyaobao 的 64 个金币 ,回帖就立即获得 2 个金币,每人有 1 次机会

luyaobao

木虫 (著名写手)


[交流] 多体耗散粒子动力学(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
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名



luyaobao(金币+2): 谢谢参与
送红花一朵
本帖仅楼主可见
20楼2021-06-08 10:39:09
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页
简单回复
tzynew2楼
2020-12-31 10:38   回复  
luyaobao(金币+2): 谢谢参与
t
2020-12-31 12:31   回复  
luyaobao(金币+2): 谢谢参与
halomilu4楼
2020-12-31 12:33   回复  
luyaobao(金币+2): 谢谢参与
donghelei5楼
2020-12-31 12:37   回复  
luyaobao(金币+2): 谢谢参与
2020-12-31 13:18   回复  
luyaobao(金币+2): 谢谢参与
2020-12-31 16:07   回复  
luyaobao(金币+2): 谢谢参与
8
syhorchid8楼
2020-12-31 17:25   回复  
luyaobao(金币+2): 谢谢参与
2021-01-02 16:37   回复  
luyaobao(金币+2): 谢谢参与
2021-01-02 16:48   回复  
luyaobao(金币+2): 谢谢参与
suihg11楼
2021-01-03 06:13   回复  
luyaobao(金币+2): 谢谢参与
谢谢分享 发自小木虫Android客户端
2021-01-03 17:12   回复  
psylhh13楼
2021-01-04 00:21   回复  
luyaobao(金币+2): 谢谢参与
2021-01-04 13:07   回复  
luyaobao(金币+2): 谢谢参与
送红花一朵
感谢分享
假大空15楼
2021-01-04 14:00   回复  
luyaobao(金币+2): 谢谢参与
jiaoxg16楼
2021-01-05 09:59   回复  
luyaobao(金币+2): 谢谢参与
2021-01-06 09:39   回复  
luyaobao(金币+2): 谢谢参与
1
gter_wang18楼
2021-01-06 13:41   回复  
luyaobao(金币+2): 谢谢参与
guanlianwu19楼
2021-01-06 19:15   回复  
luyaobao(金币+2): 谢谢参与
1
Xu_dong21楼
2021-11-20 14:18   回复  
相关版块跳转 我要订阅楼主 luyaobao 的主题更新
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复(可上传附件)
信息提示
请填处理意见