24小时热门版块排行榜

返回列表

zhangyanzhi

木虫 (小有名气)

应助: 2 (幼儿园)
金币: 1632.8
散金: 988
红花: 2
帖子: 216
在线: 294.4小时
虫号: 1787309
注册: 2012-04-30
性别: GG
专业: 多相流热物理学

[交流] 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

小弟在此期待大神的帮助!

发自小木虫Android客户端

回复此楼