| 查看: 5112 | 回复: 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 |
回复此楼» 猜你喜欢
291分调剂
已经有9人回复
-
调剂求收留
已经有34人回复
-
291 求调剂
已经有38人回复
-
22408 312求调剂
已经有17人回复
-
一志愿华中农业071010，320求调剂
已经有6人回复
-
290调剂生物0860
已经有41人回复
-
291求调剂
已经有3人回复
-
211本科材料化工求调剂
已经有23人回复
-
山东省基金2026
已经有9人回复
-
药学求调剂
已经有13人回复
【答案】应助回帖
|
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
