在之前的专栏中,曾介绍过动态统计某一个固定区域内原子数量的方法:
也介绍过动态统计某一个原子周围原子数量的方法:
下面介绍第三种类型:动态统计某一个原子周围不同类型原子数量的方法。
以小球的随机碰撞为例,原子类型为1的红色小球(id也为1)周围分布原子类型为2、3的蓝色和黄色小球。
分别统计红色小球周围蓝色原子和红色原子的数量,需要用到count(type,region)命令,改命令可统计区域region内type类型原子的数量。
对应的代码为:
fix 1 all nve
group type2 type 2
group type3 type 3
variable x1 equal x[1]
variable y1 equal y[1]
variable z1 equal z[1]
region ss sphere v_x1 v_y1 v_z1 15
group mobile dynamic all region ss every 10
variable number_all equal count(mobile)
variable number_type2 equal count(type2,ss)
variable number_type3 equal count(type3,ss)
thermo_style custom step v_number_all v_number_type2 v_number_type3
thermo 100
run 10000
统计结果如下图所示:
报名一对一培训课程可领完整版《lammps中文教程》,配套学习更轻松!
扫码领取粉丝专属优惠券↓↓↓