hallo everyone,
I would like to pack a mixture of two kinds of powders, 5000 for atom type 1 and 2000 for atom type 2. But the text file showed that there were only 5000 in total. Can anyone tell me how I can deal with it?
Best regards,
Sang
Text file:
Memory usage per processor = 12.5642 Mbytes
Step Atoms KinEng 1 ts[1] ts[2] Volume
1 5000 0.00084287269 1.7499841e-10 0 0 1.040504e-05
1000 5000 0.001850523 2.7148191e-07 0.043612943 0.025235621 1.1567814e-05
2000 5000 0.0018406285 2.6940846e-07 0.043612943 0.025235587 1.7389539e-05
3000 5000 0.0018394992 2.7286987e-07 0.043612943 0.025235553 2.9896268e-05
4000 5000 0.00184109 2.7289698e-07 0.043612943 0.02523552 5.0342628e-05
5000 5000 0.001842285 2.7289254e-07 0.043612943 0.025235486 7.8247758e-05
6000 5000 0.0018439144 2.7290271e-07 0.043612943 0.025235452 0.00011384118
7000 5000 0.0018455704 2.7290271e-07 0.043612943 0.025235418 0.0001570789
Input file:
atom_style granular
boundary m m m
newton off
atom_modify map array
communicate single vel yes
units si
processors 2 1 2
region domain block -0.0051 0.0051 0 0.1 -0.0051 0.0051 units box
create_box 2 domain
neighbor 0.0002 bin
neigh_modify delay 0
#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 1.e9 3.e8
fix m2 all property/global poissonsRatio peratomtype 0.3 0.23
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.64 0.68 0.68 0.64
fix m5 all property/global characteristicVelocity scalar 1.
fix m6 all property/global cohesionEnergyDensity peratomtypepair 2 1.e4 2.e4 2.e4 1.e4
fix m4 all property/global coefficientFriction peratomtypepair 2 0.01 0.05 0.05 0.01
fix m7 all property/global coefficientRollingFriction peratomtypepair 2 0.6 0.8 0.8 0.6
#New pair style
pair_style gran model hertz tangential history cohesion sjkr rolling_friction epsd2
pair_coeff * *
timestep 0.0000001
fix 1 all nve/sphere
fix 2 all gravity 9.81 vector 0.0 -1.0 0.0
#import triangular mesh
fix cad1 all mesh/surface file mesh4/pdcube.stl type 1 scale 0.001
#use the imported mesh as granular wall
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1
#group
group SS type 1
group SiN type 2
#distributions for insertion
fix pts1 SS particletemplate/sphere 10007 atom_type 1 density constant 7800 radius constant 0.00015
fix pts2 SS particletemplate/sphere 10009 atom_type 1 density constant 7800 radius constant 0.00021
fix pts3 SS particletemplate/sphere 10037 atom_type 1 density constant 7800 radius constant 0.000285
fix pts4 SS particletemplate/sphere 10039 atom_type 1 density constant 7800 radius constant 0.000345
fix pts5 SiN particletemplate/sphere 10061 atom_type 2 density constant 3200 radius constant 0.00010
fix pts6 SiN particletemplate/sphere 10067 atom_type 2 density constant 3200 radius constant 0.00014
fix pts7 SiN particletemplate/sphere 10069 atom_type 2 density constant 3200 radius constant 0.000165
fix pts8 SiN particletemplate/sphere 10079 atom_type 2 density constant 3200 radius constant 0.000265
fix pdd1 SS particledistribution/discrete/numberbased 10091 4 pts1 0.45 pts2 0.35 pts3 0.15 pts4 0.05
fix pdd2 SiN particledistribution/discrete/numberbased 10093 4 pts5 0.45 pts6 0.35 pts7 0.15 pts8 0.05
#region and insertion
region bc block -0.0048 0.0048 0.0005 0.09 -0.0048 0.0048
#particle insertion
fix ins1 SS insert/pack seed 32452843 distributiontemplate pdd1 vel constant 0. -1.0 0. &
insert_every once overlapcheck yes region bc ntry_mc 10000 particles_in_region 5000
fix ins2 SiN insert/pack seed 32452867 distributiontemplate pdd2 vel constant 0. -1.0 0. &
insert_every once overlapcheck yes region bc ntry_mc 10000 particles_in_region 2000
fix ts all check/timestep/gran 1000 0.1 0.1
compute 1 all erotate/sphere
compute VoroRadius all property/atom radius
variable r2 atom c_VoroRadius
compute voro all voronoi/atom radius v_r2
thermo_style custom step atoms ke c_1 f_ts[1] f_ts[2] vol
thermo 1000
thermo_modify lost ignore norm no
#dump commands
run 1
dump dmp all custom 100000 post/die_*.dump id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
#insert particles
run 1000000 upto
mschramm | Thu, 07/14/2022 - 00:24
insert/pack
Hello, your issue is with the insert/pack command. The pack version will only allow a maximum number of particles. It doesn't care what type of particle.
I was able to get your desired output by setting particles_in_region to 7000 for ins2.