|
LAMMPS讲解52-润湿in文件讲解上一节中的水分子模型采用TIP4P/2005模型。在弛豫完成后通过将水滴向下移动一段距离使其与下壁面接触从而自发润湿。在移动水滴时要调试移动距离以确保水滴下端距离下壁面足够近,但又不会引发原子重叠,一般距离3埃左右为宜。 in文件内容为 ########基本模拟规则########
dimension 3 boundary p p f units real atom_style full pair_style lj/cut/tip4p/long 1 2 1 1 0.1546 8.5 bond_style harmonic angle_style harmonic kspace_style pppm/tip4p 0.0001 kspace_modify slab 3.0
########定义几何模型########
read_data "system.data"
#########定义物理模型######### #type 1=upwall type 2=downwall type 3=O type 4=H
pair_coeff 1 1 0.1852 3.1589 pair_coeff 2 2 0.0 2.058 pair_coeff 3 3 0.1852 3.1589 pair_coeff 4 4 0.1852 3.1589 bond_coeff 1 600.0 0.9572 angle_coeff 1 75.0 104.52 ######定义分组信息#############
group water type 3 4 group upwall type 1 group downwall type 2 #########弛豫过程#########
fix 1 upwall setforce 0 0 0 fix 2 downwall setforce 0 0 0 minimize 1.0e-5 1.0e-7 1000 10000
unfix 1 unfix 2
compute myT water temp timestep 1.0 thermo 1000 thermo_style custom step temp thermo_modify temp myT flush yes fix fShakeSPCE water shake 0.0001 10 0 b 1 a 1 velocity water create 300 372748 temp myT rot yes dist gaussian dump mydump all custom 10000 dump.lammpstrj id type x y z fix fxnvt water nvt temp 300.0 300.0 100.0 run 1000000 ######润湿过程######## displace_atoms water move 0 0 -14 units box
run 1000000
感谢鲍路瑶老师的分享,内容来自于鲍老师分享出来的资料 如有需要添加微信:lmp_zhushou 进入微信群,帮助他人,共建社区 获取完整版lammps讲义可以加微信lmp_zhushou或加入QQ群994359511 |