| 查看: 1497 | 回复: 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客户端 |
回复此楼» 猜你喜欢
生物学308分求调剂（一志愿华东师大）
已经有12人回复
-
调剂求助（生物与医药）
已经有3人回复
-
一志愿吉大化学327求调剂
已经有7人回复
-
一志愿生物与医药，296分，求调剂
已经有13人回复
-
327求调剂
已经有5人回复
-
初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文
已经有8人回复
-
327考研调剂推荐
已经有6人回复
-
328求调剂
已经有3人回复
-
304求调剂
已经有9人回复
-
化工求调剂！
已经有16人回复
awuli张
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 253
- 散金: 30
- 帖子: 97
- 在线: 22.5小时
- 虫号: 15793324
- 注册: 2019-07-31
- 专业: 化学工程及工业化学
2楼2019-10-18 18:42:12
