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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps中气液表面张力计算，得不到稳定的结果
【调剂】北京石油化工学院2024年16个专业接受调剂
61/1返回列表
查看: 969  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者dbsydxdm将赠送您 15 个金币

dbsydxdm

铜虫 (小有名气)

[求助] lammps中气液表面张力计算，得不到稳定的结果

想问下大家，我想要计算气液表面张力，模型为气液气盒子。将z方向分块，使用ave/chunk命令输出每块法向力和切向力的差值rk最后加和得到表面张力，但是在气相和液相区域的切向力和法向力波动比较大，特别是气相。该怎么得到稳定的表面张力值？还有就是得到每块的stress/atom需要乘以每块原子数得到应力值吗，还是直接输出的每块stress/atom就是应力值。

compute p all stress/atom tmid
variable pn atom -c_p[3]/(20*20)
variable pt atom -0.5*(c_p[1]+c_p[2])/(20*20)
variable rk atom (v_pn-v_pt)

compute cc1 zall chunk/atom bin/1d z center 0.25 units box
fix st zall ave/chunk 10 1000 10000 cc1 v_rk v_pn v_pt file st.dat
回复此楼

» 猜你喜欢
1楼2020-12-01 15:39:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbsydxdm

铜虫 (小有名气)
有没有做表面张力的 可以交流下
赞一下(1人) 回复此楼
2楼2020-12-11 15:42:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

santiago℃

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dbsydxdm at 2020-12-11 15:42:11
有没有做表面张力的 可以交流下

我也在做表面张力，不过我是刚安装上软件的那种程度
赞一下 回复此楼
3楼2021-01-04 13:09:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbsydxdm

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by santiago℃ at 2021-01-04 13:09:04
我也在做表面张力，不过我是刚安装上软件的那种程度...

大概就是这么算的用上面的程序
赞一下 回复此楼
4楼2021-01-05 20:53:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mpiyq

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by dbsydxdm at 2020-12-11 15:42:11
有没有做表面张力的 可以交流下

有问题可以交流下  2289968992
赞一下 回复此楼
5楼2021-03-15 17:01:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

本帖仅楼主可见
赞一下
6楼2021-06-07 15:59:39
 已阅   申请模拟EPI   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 dbsydxdm 的主题更新
61/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博求助 +4 dskabdh 2024-04-24 10/500 2024-04-25 20:31 by 药嘿1233
[论文投稿] Nature一直在编辑手里，考虑好几天了，是悬了吗 +12 彩虹初见 2024-04-24 12/600 2024-04-25 19:21 by 雪径踏青
[教师之家] 刚博士毕业，入职普通一本高校 +18 保持肌肤 2024-04-20 19/950 2024-04-25 18:30 by lyfbangong
[考博] 25年博士申请 +5 Changzixuan 2024-04-25 5/250 2024-04-25 17:37 by 那年的旧城
[基金申请] 基金开始函评了吗？ +11 wych1103 2024-04-25 11/550 2024-04-25 16:54 by 未岸亡灵
[教师之家] 期末给学生划重点都是什么话术啊 +16 luokereng 2024-04-20 18/900 2024-04-25 15:46 by BusyGer
[考博] 求博导 +6 好okjh 2024-04-21 10/500 2024-04-25 14:04 by 好okjh
[论文投稿] 一直找不到审稿人 +5 lizhengke06 2024-04-21 6/300 2024-04-25 14:01 by chongdong
[硕博家园] 聊天 +11 暮色恋伊人 2024-04-22 12/600 2024-04-25 13:53 by UCTS
[基金申请] 国社科项目，你们学校都限额申报吗？ +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[博后之家] 南京大学-广州大学联合招聘博士后 欢迎广大优秀人才！！！ +4 黑魔变身啾 2024-04-20 12/600 2024-04-25 11:18 by dodonaomi
[论文投稿] Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？ 10+4 yifeng11 2024-04-20 10/500 2024-04-25 02:14 by yifeng11
[访问学者] CSC的访问学者申请，没有个评审意见，也不知道怎么改，还有必要申请吗 +4 flyingship 2024-04-20 4/200 2024-04-24 21:09 by 59038mute
[基金申请] 化学结清有情发出来了 +3 starboy7286 2024-04-20 3/150 2024-04-24 15:52 by 嘿，黑贝
[有机交流] 紧急求助，有谁用过三甲基硅醇钾TMSOK吗？ 1000+3 genius2008 2024-04-19 10/500 2024-04-24 12:37 by genius2008
[教师之家] 大家访学都是怎么找的啊？ +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[论文投稿] 无under review，直接DIP，咋么办。 +10 lizhengke06 2024-04-19 16/800 2024-04-23 19:15 by 化学程序员
[教师之家] 河海大学常州校区新进的讲师待遇如何？考核要求如何？ +3 德国小球迷 2024-04-19 18/900 2024-04-22 19:23 by ZHONGWU_U
[论文投稿] 研究光催化的，好中的三四区 20+3 sl.0117 2024-04-20 3/150 2024-04-22 09:53 by bobvan
[论文投稿] 编辑是选国外的好还是国内的好。 +8 lizhengke06 2024-04-20 8/400 2024-04-22 08:58 by cuiyunjian
信息提示
关闭
请填处理意见
关闭