| 查看: 5071 | 回复: 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 |
回复此楼» 猜你喜欢
材料与化工304求B区调剂
已经有7人回复
-
0703化学336分求调剂
已经有4人回复
-
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
-
化学工程321分求调剂
已经有12人回复
-
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有4人回复
-
材料专硕326求调剂
已经有5人回复
-
东南大学364求调剂
已经有5人回复
-
国自科面上基金字体
已经有7人回复
-
药学383 求调剂
已经有4人回复
-
286求调剂
已经有5人回复
【答案】应助回帖
|
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
20