LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Polydispsersity and random packing
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Polydispsersity and random packing


From: Steve Plimpton <sjplimp@...24...>
Date: Tue, 13 Jun 2017 08:39:00 -0600

Then you can setup any form of polydispersity you wish.
Each particle has its own diameter.

Steve

On Tue, Jun 13, 2017 at 8:36 AM, Roberts, Andrew <A.P.Roberts@...5884...> wrote:

Hi Steve,


My pair style is granular:


pair_style         gran/hooke/history 10000 5000 0.01672 0 0 0        


Kind regards,


Andrew


From: Steve Plimpton <sjplimp@...24...>
Sent: 13 June 2017 15:28:53
To: Axel Kohlmeyer
Cc: Roberts, Andrew; lammps-users@...396...sourceforge.net
Subject: Re: [lammps-users] Polydispsersity and random packing
 
What pair style are you planning to use?
Granular or something else?  Only a few
pair styles will work with randomized polydispersity.

Steve

On Mon, Jun 12, 2017 at 11:39 AM, Axel Kohlmeyer <akohlmey@...29....> wrote:


On Mon, Jun 12, 2017 at 12:24 PM, Roberts, Andrew <A.P.Roberts@...5884...> wrote:

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.​

 


The script I am working on is attached.


Kind regards,


Andrew


------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@...655....net
https://lists.sourceforge.net/lists/listinfo/lammps-users




--
Dr. Axel Kohlmeyer  akohlmey@...24...  http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.

------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@...655....net
https://lists.sourceforge.net/lists/listinfo/lammps-users