| 查看: 1494 | 回复: 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客户端 |
回复此楼» 猜你喜欢
304求调剂
已经有4人回复
-
326求调剂
已经有4人回复
-
一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项
已经有4人回复
-
环境工程调剂
已经有5人回复
-
本人考085602 化学工程 专硕
已经有14人回复
-
一志愿211 0703方向310分求调剂
已经有3人回复
-
0703化学336分求调剂
已经有3人回复
-
318求调剂
已经有3人回复
-
化学工程321分求调剂
已经有9人回复
-
0703化学调剂
已经有9人回复
awuli张
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 253
- 散金: 30
- 帖子: 97
- 在线: 22.5小时
- 虫号: 15793324
- 注册: 2019-07-31
- 专业: 化学工程及工业化学
2楼2019-10-18 18:42:12
20