| 查看: 4588 | 回复: 6 | ||
| 【悬赏金币】回答本帖问题,作者dbsydxdm将赠送您 15 个金币 | ||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||
[求助]
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 |
» 猜你喜欢
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
现代”学阀”该如何界定
已经有14人回复
求合成方法
已经有8人回复
师弟论文见刊大半年才想起来申请专利,还能抢救一下吗?
已经有4人回复
有机合成以后会不会被AI改变?做科研的虫友怎么看
已经有8人回复
我的奶奶
已经有3人回复
6楼2021-06-07 15:59:39
2楼2020-12-11 15:42:11
santiago℃
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 375.4
- 散金: 92
- 沙发: 1
- 帖子: 204
- 在线: 13.7小时
- 虫号: 12896718
- 注册: 2018-11-18
- 性别: GG
- 专业: 多相流热物理学
3楼2021-01-04 13:09:04
4楼2021-01-05 20:53:28










回复此楼