Hello,
I am new to LAMMPS, but I wondered if someone could tell me how to specify this for a 2D Couette flow:
10% polydispersity on sphere diameter so that all the diameters (flow and walls) are randomly distributed between 0.9 and 1.1?
Random packing of spheres in the flow region, but keep the walls packed in a straight line?
The script I am working on is attached.
Kind regards,
Andrew
in.flow2.couette (1.93 KB)
Hello,
I am new to LAMMPS, but I wondered if someone could tell me how to specify
this for a 2D Couette flow:- 10% polydispersity on sphere diameter so that all the diameters
(flow and walls) are randomly distributed between 0.9 and 1.1?
this is easily done with an atom style variable and the set diameter
command. after changing diameters, you will need to run some kind of
(restricted) relaxation to remove overlaps.
- Random packing of spheres in the flow region, but keep the walls
packed in a straight line?
just use two separate create_atoms calls on the different regions and
request positions on lattice points in the first (the default) and random
positions in the second. for random positions running some unoverlap
relaxation is essentially. possibly you may need to use a different pair
style, e.g. soft, and fix adapt to handle large overlaps smoothly. see the
micelle example for how to do the latter.
axel.
What pair style are you planning to use?
Granular or something else? Only a few
pair styles will work with randomized polydispersity.
Steve
Hi Steve,
My pair style is granular:
pair_style gran/hooke/history 10000 5000 0.01672 0 0 0
Kind regards,
Andrew
Then you can setup any form of polydispersity you wish.
Each particle has its own diameter.
Steve