variable x0 equal 0.0
variable x1 equal 10
variable y0 equal 0.0
variable y1 equal 10
variable z0 equal 0.0
variable z1 equal 10
variable regx0 equal -0.0001+${x0}
variable regx1 equal 0.0001+${x1}
variable regy0 equal -0.0001+${y0}
variable regy1 equal 0.0001+${y1}
variable regz0 equal -0.0001+${z0}
variable regz1 equal 0.0001+${z1}
lattice sc .6
region sim_reg block ${regx0} ${regx1} ${regy0} ${regy1} ${regz0} ${regz1} units box ##simulation box
region insParticles block ${x0} ${x1} ${y0} ${y1} ${z0} ${z1} units box ##region in which our particles are generated
create_box 2 sim_reg
variable binSize equal .1
neighbor ${binSize} bin
variable neighDelay equal 0
variable neighOne equal 500000
variable neighPage equal 10000000
neigh_modify delay ${neighDelay} one ${neighOne} page ${neighPage}
#######################################################################
variable dia_NMC_1 equal 0.1
variable dia_NMC_2 equal 4
#############################################################
#Create atoms in the sim_reg box
create_atoms 1 region sim_reg
#mass 1 1.0
set atom 1 diameter ${dia_NMC_1}
run 1
################################################################
#Delete atoms part one
region sphere_outside sphere 5.0 5.0 5.0 5 side out #try it out with 4 may be
delete_atoms region sphere_outside
run 1
################################################################
#Delete atoms part two
region sphere_inside sphere 5.0 5.0 5.0 4.5 side in
delete_atoms region sphere_inside
run 1
##########################################################
# Create a big particle
create_atoms 2 single 5.0 5.0 5.0
# Set diameter of the big particle
set atom 2 diameter ${dia_NMC_2}
run 1
As you can see in the picture, the small particles should cover the big particle but is not happening. A particle is generated in the centre. But when I am changing the diameter of the particle, it is not changing. Another particle from the corner which is already deleted, that diameter is changing. When I comment on the delete atom parts, I can create the big particle properly.
Can someone please help me how can I solve the problem?
Attachment | Size |
---|---|
![]() | 119.21 KB |