LIGGGHTS® - User Forum

LIGGGHTS® related topics can be discussed here: discussion about models, installation, feature requests and general discussion

Avoid moding mesh using force fields

Submitted by antoniosereno on Wed, 10/16/2019 - 11:24

Hello everyone!
I'm trying to simulate a ball mill and I'm using at the moment a very simple geometry to obtain acceptable computation times since I'm learning now LIGGGHTS. While looking at the documentation page, I read about the possibility to assign variables instead of exact values to both the magnitude and the direction of hte gravity vector and I was thus wondering if there's a way to avoid heavy computation by using this trick instead of moving he mesh.

As a first attempt, I first inserted the particle into my geometry using a classical gravity vector, using:

Please suggest any way to compile the in. file without copy executable 'lmp_auto' into the same folder

Submitted by zen on Tue, 10/15/2019 - 16:47

Hi all,

I have got question about how to compile the 'in.filename' in parallel without copy executable 'lmp_auto' into the same folder. As suggested by the Liggghts documentation that the final step to make it compiled is

"cp lmp_auto ../examples/LIGGGHTS/Tutorials_public/chute_wear
cd ../examples/LIGGGHTS/Tutorials_public/chute_wear
mpirun -np 4 lmp_auto < in.chute_wear"

Large-Large overlapping in pair_colloid

Submitted by Lowered on Fri, 10/11/2019 - 17:43

Hello,

I'm facing a problem in pair_colloid.cpp
I designed a small testcase containing just 2 particles, radius 0.01µm and varying the initial distance between them. They are underlaying DLVO potential by hybrid overlay colloid and yukawa/colloid
For larger initial distances the simulation is working. At a distance of 2µm I get an unspecific error (seems like division by zero). For distances smaller than 2µm I get the error: "Overlapping large/large in pair colloid" (line 187 of pair_colloid.cpp).

Missing quad1, qaud2, qaud3, quad4 from dump file

Submitted by johnr on Thu, 10/10/2019 - 15:18

Hi, after running :

examples/LIGGGHTS/Tutorials_public/superquadric/in.particle_particle

When I open the dump files in Paraview the variables quad1, qaud2, qaud3, quad4 are missing, while shape x,y,z and blockiness 1,2 are present. I would like to check the initial settings for quad. Does someone know why these values are missing?

(The Liggghts version is 3.8.0 and I used 'make auto' in the src folder then copied the lmp_auto executable to run the code.)

Many thanks, John

How to assign number of multisphere particles into the insert region ?

Submitted by Michelle Phan on Wed, 10/09/2019 - 08:46

This might be a idiot question but I have been stuck for a while. I would like to ask how to assign number of multisphere particles into the insert region. I used the "multisphere_stone_restitution" in LIGGGHTS' example. There are 50 coordinations in the data file, and only 2 clumped particles were created. I wonder how to increase the number of clumped particles ? Should I increase the number of coordination? I read the LIGGGHTS' manual but it does not (or I missed it) mention about this problem.
Thank you for your time.

Clump particle in LIGGGHTS 3.2

Submitted by Michelle Phan on Fri, 10/04/2019 - 05:51

Hello everyone;
By some reasons, I need to use LIGGGHTS 3.2 to simulate my model. I need to clump the particles to make the irregular shape. I know the function fix particletemplate/multisphere can do that job, but it is only available in LIGGGHTS 3.3 onward. I tried to copy the fix_multisphere in src 3.3 to 3.2, but I could not compile. Can anyone suggest me a solution?
Thank you

how to make particles position fixed?

Submitted by zhouzunxie@gmail.com on Mon, 09/23/2019 - 06:21

Hello everyone,

I am working on liggghts, and I generated a particle packing in the geometry by command insert/stream, after particles become stationary in the bottom, i would like to make the packing position fixed and not be moved anymore. In other words, after packing formed, we treated it as a wall.

I tried the command fix freeze and fix setforce, but the particles not acted force tend to overlap with each other, and it is unexpected.

Do you know how to realize this ?

thanks.

Number of processors

Submitted by jvvelez on Sun, 09/22/2019 - 05:47

I am simulating shaking grains in the vertical direction and I would like to speed up my simulation. I can use up to 24 cores. Since, the motion is predominantly in the z direction, should I maximize number of cores in this direction? Also, I noticed if I reduced the number of cores and let liggghts decide the optimum distribution, 20 cores make it faster instead of 24. I thought it was counter-intuitive but then again it could just be how liggghts distributes it. What should I do?

How to relax configuration with overlapcheck no

Submitted by bern on Thu, 09/19/2019 - 15:34

Hello everyone,

I want to simulate with a volume fraction of 0.4288 but this would lead to unreasonable high numbers for maxattempt to work.
So deactivated overlapcheck, but then the manual says the overlapping configuration has to be relaxed.
My problem: there is no word on how to do this.

Can anyone tell me how to relax properly? ;)

Thank you in advance.

Regarding contact with particles

Submitted by Rohith tvn on Thu, 09/19/2019 - 13:55

Hello everyone,
I am a student doing a project work, as part of my task I need to calculate traction when my sample come in contact with soil. Here the problem is particles are inserted and now there is no contact of my sample with my particles as sample is just at the surface but due to clearance between sample and particles there is no contact. Can anyone help me with this to make contact between sample and particles.
Thanks & regards.
Rohith.

Pages

Subscribe to RSS - LIGGGHTS® - User Forum