Mixing two types of powders

Submitted by Sang on Thu, 06/30/2022 - 15:38

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

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.