Dear Lammps users,
Good day. I am facing a problem during creating a cylinder inside the rectangular simulation box. The cylinder is not showing anywhere in the box. I am looking for some suggestions regarding this. Thank you.
Here is the code:
atom_style sphere #Spheres has per-particle diameter and mass
boundary p p fm
newton off #Two interacting atoms are on different processors and the resulting force information is not communicated.
comm_modify vel yes #Conserving properties to two ghost atoms
region reg block -20 20 -10 10 -0.5 16 units box
create_box 1 reg
region cyl cylinder y 0 0 5 0 5 units box
Regards,
Yeasir Mohammad Akib
MS Student
Manufacturing Engineering
University of Texas Rio Grande Valley
Dear Lammps users,
Good day. I am facing a problem during creating a cylinder inside the rectangular simulation box. The cylinder is not showing anywhere in the box. I am looking for some suggestions regarding this. Thank you.
what do you mean by “not showing”?
axel.
Dear Dr. Axel,
There is no cylinder present in the simulation box.
Regards,
Yeasir Mohammad Akib
MS Student
Manufacturing Engineering
University of Texas Rio Grande Valley
Dear Dr. Axel,
There is no cylinder present in the simulation box.
how do you determine this?
Dear Dr. Axel,
There are more parts of the code, after running the whole code and opening the file in the OVITO, I found no cylinder in the simulation box.
Here is the full code:
atom_style sphere #Spheres has per-particle diameter and mass
boundary p p fm
newton off #Two interacting atoms are on different processors and the resulting force information is not communicated.
comm_modify vel yes #Conserving properties to two ghost atoms
region reg block -20 20 -10 10 -0.5 16 units box
create_box 1 reg
region cyl cylinder y 0 0 5 0 5 units box
neighbor 0.2 bin
neigh_modify delay 0
IMPORTANT NOTE: these values of stiffness (2000) and timestep (0.001)
are used in this example file to produce a quick simulation and movie.
More appropriate values for realistic simulations would be
k = 2.0e5 and dt = 0.0001, as in bench/in.chute.
pair_style gran/hooke/history 2000.0 NULL 50.0 NULL 0.5 0
pair_coeff * *
timestep 0.001
fix 1 all nve/sphere
fix 2 all gravity 1.0 spherical 0.0 -180.0
fix zlower all wall/gran hooke/history 2000.0 NULL 50.0 NULL 0.5 0 &
zplane 0.0 2000.0
region slab block -19.0 -9.5 -9.0 9.0 10.0 15.0 units box
fix ins all pour 10000 1 300719 vol 0.13 50 region slab
compute 1 all erotate/sphere #erotate/sphere comutes the rotational kinetic energy
thermo_style custom step atoms ke c_1 vol #c_1 is coumpted with compute id
thermo 100
thermo_modify lost ignore norm no #ignore value does not check for lost atoms, norm value normalized the thermodyanmic paprameters
compute_modify thermo_temp dynamic/dof yes
dump id all atom 100 dump.pour
dump 1 all custom 100 positions.lammpstrj id x y z vx vy vz
axes yes 0.8 0.02 view 80 -30
#dump_modify 2 pad 5
dump 2 all movie 100 movie.mpg type type &
axes yes 0.8 0.02 view 80 -30
#dump_modify 3 pad 5
run 100000
unfix ins
fix 2 all gravity 1.0 chute 26.0
run 5000
Regards,
Yeasir Mohammad Akib
MS Student
Manufacturing Engineering
University of Texas Rio Grande Valley
Dear Dr. Axel,
There are more parts of the code, after running the whole code and opening the file in the OVITO, I found no cylinder in the simulation box.
where in the documentation for the dump command(s) does it say that the dump files will contain information about regions?
a region is just an abstract construct for use with other LAMMPS commands, so you can only “see” that region indirectly.
e.g. when you add a command like:
create_atoms 1 region cyl
it will fill the intersection of your box with that region with atoms according to the currently defined lattice. so you can see the presence of the region from the added atoms.
you can also list all defined regions with their boundary information.
info region
command.
however, in your complete input I don’t see anything making use of the region and thus it will not do anything.
axel.