| 查看: 5292 | 回复: 3 | ||
zhangyanzhi木虫 (小有名气)
|
[求助]
Lammps如何计算气液界面的表面张力及其分布？ 已有1人参与
|
|
求助大神，如何计算气液界面的表面张力及其分布。下面是我in文件中的相应部分，但总是计算不出文献中表面张力的分布规律。 #--------------------calculate surface tension----------------- variable N equal 100 variable box equal (xhi-xlo)*(yhi-ylo)*(zhi-zlo) # volume of the box #--------------------surface tension(x-direction)-------------- compute p all stress/atom NULL compute k all reduce sum c_p[1] c_p[2] c_p[3] variable pn atom -c_p[1] variable pt atom -(c_p[2]+c_p[3])/2.0 variable rk atom ((v_pn-v_pt)/(v_box/v_N))*((xhi-xlo)/v_N)*1.0e-2 # unit mN/m #--------------------surface tension(global)------------------- variable tension equal 0.5*(xhi-xlo)*(c_k[1]/v_box-(c_k[2]+c_k[3])/(2.0*v_box))*1.0e-2 # unit mN/m fix mean all ave/time 100 5 1000 v_tension mode scalar ave one file tension.dat 小弟在此期待大神的帮助!@smutao@oxox6085 |
回复此楼» 猜你喜欢
朋友圈看到的
已经有4人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
【答案】应助回帖
|
compute p all stress/atom NULL 上面的all如果改成液体的group,体积也只使用液体的体积，是不是更准确一些 参考 http://blog.sina.com.cn/s/blog_c0468c8f0101tfea.html 以及 http://www.afs.enea.it/software/ ... e_voronoi_atom.html 另外我也在尝试计算表面张力，，但是compute voronoi/atom 这个命令还没有完全搞清楚，如果你会和我也说一下啊 |
2楼2017-10-28 20:38:38
3楼2020-04-18 21:37:06
4楼2020-11-04 14:50:22
