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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 多体耗散粒子动力学（MDPD）液体表面张力的LAMMPS模拟计算
211/1返回列表
查看: 2842  |  回复: 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分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭