24小时热门版块排行榜    

查看: 3777  |  回复: 6
【悬赏金币】回答本帖问题,作者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
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 7 个回答

dbsydxdm

铜虫 (小有名气)

有没有做表面张力的 可以交流下
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的回帖

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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 267化工调剂求助 +5 聪少OZ 2026-03-04 5/250 2026-03-04 20:33 by kakakapanpan
[考研] 一志愿中石油,化学工程,301求调剂 +3 豚馒枳枳 2026-03-04 3/150 2026-03-04 18:32 by houyaoxu
[考研] 材料专硕346求调剂 +3 旺一下 2026-03-04 3/150 2026-03-04 16:26 by sslc1985
[考研] 一志愿中南大学理学化学 +6 15779376950 2026-03-01 8/400 2026-03-04 16:19 by zhukairuo
[考研] 0857求调剂 +4 CDDD1013 2026-03-04 4/200 2026-03-04 16:08 by zhukairuo
[考研] 281求调剂 +3 Y?l?h 2026-03-04 3/150 2026-03-04 10:51 by 汪!?!
[考研] 材料工程269求调剂 +7 白刺玫 2026-03-02 7/350 2026-03-04 08:47 by zhyzzh
[考研] 0703化学学硕 288分 +4 旅行中的紫葡萄 2026-03-03 4/200 2026-03-03 23:41 by ch1979
[考研] 化工专硕调剂 +4 利好利好. 2026-03-03 7/350 2026-03-03 21:30 by L135790
[考研] 290求调剂 +9 ErMiao1020 2026-03-02 9/450 2026-03-03 18:03 by linlonghao
[考研] 266材料化工求调剂 +3 哇塞王帅 2026-03-03 3/150 2026-03-03 17:58 by houyaoxu
[考研] 化工专硕348,一志愿985求调剂 +8 弗格个 2026-02-28 11/550 2026-03-03 14:55 by relum
[考研] 267求调剂 +6 钓鱼佬as 2026-03-02 6/300 2026-03-03 13:59 by 13589
[考研] 299求调剂 +5 kkcoco25 2026-03-02 9/450 2026-03-03 12:55 by 公瑾逍遥
[考研] 材料工程专硕283求调剂 5+8 ,!? 2026-03-02 10/500 2026-03-03 12:01 by EBSD
[考研] 290分材料工程085601求调剂 数二英一 +8 llx0610 2026-03-02 9/450 2026-03-02 22:09 by 无际的草原
[考研] 材料工程274求调剂 +5 Lilithan 2026-03-01 5/250 2026-03-02 19:39 by caszguilin
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
信息提示
请填处理意见