| 查看: 1496 | 回复: 1 | |||
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 小弟在此期待大神的帮助! 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
277求调剂 数一104分
已经有7人回复
-
288环境专硕,求调材料方向
已经有16人回复
-
277工科求调剂
已经有8人回复
-
266求调剂
已经有19人回复
-
环境专硕调剂
已经有3人回复
-
081200-11408-276学硕求调剂
已经有5人回复
-
材料调剂
已经有13人回复
-
工科 22408 267求推荐
已经有4人回复
-
305求调剂
已经有3人回复
-
270分求调剂
已经有5人回复
awuli张
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 253
- 散金: 30
- 帖子: 97
- 在线: 22.5小时
- 虫号: 15793324
- 注册: 2019-07-31
- 专业: 化学工程及工业化学
2楼2019-10-18 18:42:12
