LIGGGHTS® - User Forum

Hi,

While setting up a case for a CFDEM simulation, we encountered a problem with the perceived packing density.

We first inserted the particles mid-air and let them settle. The bed height after settling was much higher than in reality, which is why we tried a higher gravity as well as shrinking the particles and cycling between growing and relaxing, to no avail.

Has anyone had a similar issue before? We would appreciate any response.

Hello guys,

I want to use DLVO forces in my simulation model. Therefore I need to use the pair_style colloid and yukawa, I guess.
I have read the manual of LAMMPS (http://lammps.sandia.gov/doc/Section_packages.html), but I still don't know, how to install the package colloid.

Hi,

I have a question to the combination of Hertz contact model and rolling model. Consider a particle with high overlap rotating in elastic region -> that results in a stiffer k_n and k_t -> results in a stiffer k_r in rolling behaviour. Now the overlap is reduced and the particles roll back -> kn kt is lower, also kr is lower -> results in a larger back rotation because of the lower k_r now. Attached is a sketch of this rolling case. So do I have to use hook if I want to use the epsd rolling model?

kind regards
Joshua

I'm trying to run a compilation using a hybrid granular molecular atom_style, but when i run LIGGGHTS, this error appears:

ERROR: Invalid atom_style command (../atom_vec.cpp:72)

I'm trying to solve this problem but even reading the documentation, i cannot find a solution.

PD: I'm using the 3.2.0 LIGGGHTS version

Cheers, Nicolás

Is there a way to simulate multispheres in LIGGGHTS without the particletemplate/multisphere command?. I'm currently working in a SAG mill simulation with a partner, and our boss ask to do the same simulation with multispheres... please, a hint would be more than useful.

hi
i am a new user of liggghts 3.x and want to pour my particles throgh:
fix ins nve_group pour/dev nparticles 16200 1 distributiontemplate pdd1 vol 0.3 10 massflowrate 20. vel uniform 0. 0. 0. 0. -0.5 region bc
or
fix ins nve_group pour/dev nparticles 16200 1 distributiontemplate pdd1 vol 0.3 10 vel uniform 0. 0. 0. 0. -0.5 region bc

Hello all,
any one came accros this issue? after compilation, then i launched a simulation and found an error : /lib64/libc.so.6: version `GLIBC_2.14' not found
This happend when i updated ligggghts to 3.X version.
thanks,
Zyan

Hello community,

I am interested in the new "modify_timing" command, but I can't find any clue in the manual of version 3.2.1. Is it a standalone fix or belongs "modify_timing" to another command??

Best regards
Niels

Hi,

I am trying to extract the shear stress on every mesh element of a surface/mesh/stress using "dump mesh/vtk" as instructed in the manual. The problem here is that I can't understand the structure of the output file. Do you have any tip or do you know any better way in this regard.

Thank you.
Mohammad

To simulate a scoop I would need the scoop mesh to move over a fixed distance and to rotate for a fixed angle sequentially. I tried to use a timestep calculated like this

variable deltatime equal sqrt(${rho}/($E/(2*(1+$v))))*c_minradi*PI/(0.1631*$v+0.876605)/5

and the set velocity of the mesh used in the

fix movecad1 all move/mesh mesh cad1 linear 0.0 0.0 ${scm1d}

to get the required number of timesteps