| 查看: 5063 | 回复: 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 |
回复此楼» 猜你喜欢
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有9人回复
-
版面费该交吗
已经有9人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
面上可以超过30页吧？
已经有4人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有5人回复
-
什么是人一生最重要的？
已经有4人回复
【答案】应助回帖
|
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