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
SHUBHAM AGARWAL | Fri, 07/23/2021 - 23:09
Try group
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,
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.
mschramm | Sat, 07/24/2021 - 22:16
first in first out
Hello,
fixes are completed in the order that they are set. So ins1 will run before ins2.
If there is no more space or 500 particles are inserted, ins2 will not run.