Hi everyone,
I am going to simulate granite using LIGGGHTS-with-bonds with particles of different radii. I need a dense packing of particles. When I used fix particlestemplate/sphere in LIGGGHTS to fill a desired volume, even with a very high volume fraction, there was still a lot of void space between particles, and they were not in contact at all.
To address this, I tried using a MATLAB code to fill the desired irregular volume with spherical particles and then passed the coordinates and radii to LIGGGHTS using read_data command. However, to achieve dense packing, I would need to insert a large number of particles, which is not computationally efficient.
Now, my plan is to first insert particles into the desired volume and then gradually increase their radii to a target value, allowing the particles to settle by running some cycles. Once they reach equilibrium, the bonds between them will be defined.
For this, I need to define fixed walls with an irregular shape around the particles that match the volume shape. Is it possible to define an irregularly shaped wall around particles in LIGGGHTS?
Thanks,
Tannaz