24小时热门版块排行榜    

查看: 3162  |  回复: 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分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 6/300 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +6 启萌科技 2026-08-12 21/1050 2026-08-14 19:48 by 启萌科技
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +6 Tide man 2026-08-14 7/350 2026-08-14 18:36 by 正能量1斤
[基金申请] filecode +6 cratir 2026-08-14 10/500 2026-08-14 18:29 by 笑叹辞穷
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见