granular wall

raguelmoon's picture
Submitted by raguelmoon on Thu, 04/28/2011 - 18:04

Hi LIGGGHTS users,
I want to make walls of my simulation box with granular media but couldn't do that. Can anybody help me doing so?
Thanks a lot,
R

raguelmoon's picture

raguelmoon | Fri, 04/29/2011 - 18:42

Hi,
I mean wall of spherical particles of radius more than the radius of granular particles in box. I poured particles in simulation box and then deleted atoms by region which is a little smaller ( bye one particle diameter) than actual box size. This might leave a thin layer (wall) of one particle thickness around box boundary. But i see that there is gap between particles and no regular arrangement was observed which may be used as wall. Then i can apply fix rigid command to make it non-dynamical.
Thanks,
Ram

Ram

ckloss's picture

ckloss | Mon, 05/02/2011 - 08:44

to do what you want, typically the lattice command is used. alternatively, you can use read_data to read a list of particles at pre-defined positions

Christoph

raguelmoon's picture

raguelmoon | Mon, 05/09/2011 - 19:02

Dear Christoph,
It is very difficult to use lattice command in LIGGGHTS. Do you know how can we define radius and density of the particles? SET command is written after create_atoms. So create_atoms command has no effect of changing its number of atoms if we change radius in SET command. It doesnt care if we put "set" command or not. And, how to use lattice with SI? Also reading data by read_data has some difficulties. I want to make restart_file first and then can change some positions of particles if i get data file from restart_file. But tool came with LIGGGHTS which is used for converting restart2data is not working...
Thanks,
Ram

Ram

raguelmoon's picture

raguelmoon | Tue, 05/10/2011 - 07:52

Hi,
let me try to code restart2data.cpp for LIGGGHTS... Can you tell me structure of restart file? I mean what is written out in restart file?
Thanks,
Ram

Ram

ckloss's picture

ckloss | Tue, 05/10/2011 - 08:47

>>It is very difficult to use lattice command in LIGGGHTS
I would not say that. you can control both the lattice spacing and the radius via variables, so I do not see a problem here

>>So create_atoms command...It doesnt care if we put "set"
Why should create_atoms care about the set command? There is simply no reason for that.
You have to choose a lattice spacing that fits the radius - it's as simple as that

Cheers,
Christoph