Problem with model two groups of particles with two atom types

Submitted by tuongnhu1990 on Fri, 07/23/2021 - 14:35

I try to simulate an assembly with two group of particles. However, there is only one group of particles (ins1) added to the simulation while the other (ins2) is ignored. My code is given below. Could anyone tell where I did wrong? Thank you!

atom_style granular
atom_modify map array
boundary m m m
newton off

communicate single vel yes

units si

region reg block 0.0 2.0 0.0 0.1 0.0 1.0 units box
create_box 2 reg
group nve_group region reg

neighbor 0.001 bin
neigh_modify delay 0

group pave type 1
group grain type 2

#distributions for insertion
fix pts1 pave particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.004
fix pdd1 pave particledistribution/discrete 15485867 1 pts1 1.0

fix pts2 grain particletemplate/sphere 15485863 atom_type 2 density constant 2500 radius constant 0.002
fix pdd2 grain particledistribution/discrete 15485867 1 pts2 1.0

#region for insertion
region bc block 0.0 0.2 0.0 0.1 0.0 0.7 units box
group bcg region bc

#particle insertion
fix ins1 pave insert/pack seed 32452867 distributiontemplate pdd1 &
maxattempt 10000 insert_every once overlapcheck no all_in yes vel constant 0.0 0.0 0.0 &
region bc particles_in_region 500 check_dist_from_subdomain_border no

fix ins2 grain insert/pack seed 32452843 distributiontemplate pdd2 &
maxattempt 10000 insert_every once overlapcheck no all_in yes vel constant 0.0 0.0 0.0 &
region bc particles_in_region 500 check_dist_from_subdomain_border no

#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 5.e6 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.7 0.7 0.7 0.7
fix m4 all property/global coefficientFriction peratomtypepair 2 0.5 0.5 0.5 0.5

Forums:

SHUBHAM AGARWAL | Fri, 07/23/2021 - 23:09

Hello
I don't know the exact answer, but you may try grouping the two fixes and try to insert that in single fix command.
Let me know if it works

Shubham

tuongnhu1990 | Sat, 07/24/2021 - 03:08

Thanks Shubham,

However, my final aim is to simulate a system with a group of particles running for a certain time, then the second group with different type and properties will be added. So, I guess I can't avoid having two separate inserts.