应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 多体耗散粒子动力学（MDPD）液体表面张力的LAMMPS模拟计算
211/1返回列表
查看: 3227  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 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朵）

santiago℃ ymsw_zzw

» 猜你喜欢
1楼 2020-12-31 09:28:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名


★
luyaobao(金币+2): 谢谢参与
送红花一朵
本帖仅楼主可见
赞一下
20楼2021-06-08 10:39:09
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页
简单回复
tzynew2楼
2020-12-31 10:38   回复  
luyaobao(金币+2): 谢谢参与
t
xiaoqiuqiu203楼
2020-12-31 12:31   回复  
luyaobao(金币+2): 谢谢参与
halomilu4楼
2020-12-31 12:33   回复  
luyaobao(金币+2): 谢谢参与
。
donghelei5楼
2020-12-31 12:37   回复  
luyaobao(金币+2): 谢谢参与
xishandian6楼
2020-12-31 13:18   回复  
luyaobao(金币+2): 谢谢参与
miaojiabing7楼
2020-12-31 16:07   回复  
luyaobao(金币+2): 谢谢参与
8
syhorchid8楼
2020-12-31 17:25   回复  
luyaobao(金币+2): 谢谢参与
mikeshen1019楼
2021-01-02 16:37   回复  
luyaobao(金币+2): 谢谢参与
Cynthia198610楼
2021-01-02 16:48   回复  
luyaobao(金币+2): 谢谢参与
顶
suihg11楼
2021-01-03 06:13   回复  
luyaobao(金币+2): 谢谢参与
谢谢分享 发自小木虫Android客户端
Cynthia198612楼
2021-01-03 17:12   回复  
顶
psylhh13楼
2021-01-04 00:21   回复  
luyaobao(金币+2): 谢谢参与
santiago℃14楼
2021-01-04 13:07   回复  
luyaobao(金币+2): 谢谢参与
送红花一朵
感谢分享
假大空15楼
2021-01-04 14:00   回复  
luyaobao(金币+2): 谢谢参与
jiaoxg16楼
2021-01-05 09:59   回复  
luyaobao(金币+2): 谢谢参与
苏州华博17楼
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 的主题更新
211/1返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-10-01 21:59 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 CfXuS1rDhLYN 2026-09-28 4/200 2026-10-01 09:23 by Equinoxhua
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
信息提示
关闭
请填处理意见
关闭