LIGGGHTS® - User Forum

Hello all,

Why in LIGGGHTS the cohesion density energy is in J/m3. When Johnson-Kendall-Roberts theory say that is the surface energy, and that is in J/m2?. What is the relation betwen both of that definitions?. Because in the literature the definition of the cohesion energy density is other than the adhesion, and that is related with chemestry, it seems here is a mistake in definitions. So, the process to calibrate this values with experimental data is confuse.

Please say to me what paper could I can read for understand the approach did for LIGGGHTS with this terms

Thanks

I want to count two sizes of particles separately with fix massflow/mesh command. I thought it should be very straightforward but I didn't figure out how to do that.

Thanks and best regards,
Min

Dear all,
I want to impose an isotropic compression on a granular sample. I am trying to use the servo wall for it. But in spite of many test (playing with KP, KD, KI,mode auto, etc) the walls are not getting the fixed force and also not moving. I am attaching an example script and the related cad file, in case any of you can help me.

Hello,
I installed LIGGGHTS WITH BONDS and I have problem with the visualization of the example.
I found program bondsfalloncone and I ran in.template.
https://github.com/Polyun/LIGGGHTS-PUBLIC/tree/master/bondspackage/bonds...

In this program, parameters are write to output file by command:
dump dmp all custom 250 post/dump*.bondcone id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

Hello,
Can anyone tell me how to get the elastic potential energy due to the particle-wall contact. pe/atom or pe gives only gravitational potential energy. (RESTRICTION: This compute does not include the potential energy due to the overlap of granular particles.).
There is an option 'computeElasticPotential' as a keyword for the hertz contact which saves the elastic potential energy, but I don't know how to use the saved value.
My wall is not a primitive wall it is a mesh.
Can anyone guide me through the procedure?
Thanks in advance.

Hello Liggghters,

I am trying to understand the physics behind the tutorial of Hydrogel_multicontact but while lunch it give me the following error attached in the log file below how do I make work ??

Salutations,

Dear all,
I want to simulate a case in which a gate move up. I used move mesh command.
However, the mesh did not move.

#Multisphere

atom_style sphere
atom_modify map array sort 0 0
boundary f f f
newton off

communicate single vel yes

#PUBLIC version only supports serial execution
processors 1 1 1

units si

region reg block -0.5 0.3 -0.2 0.2 0. 0.5 units box
create_box 1 reg

neighbor 0.004 bin
neigh_modify delay 0

#Material properties required for new pair styles

Hi all,

I am simulating a 3D geometry, but in the end I would like to visualize results in only a section on the geometry.
e.g. I have a cylinder with particles in in, I want to visualize a plane in the centre of the cylinder. Is there is a way this can be done in LIGGGHTS + paraview?

Thanks a lot.

Regards,
Satyajeet

Hi,

Can anyone tell whether the glyph to visualize superquadric particles from LIGGGHTS has been released or available? Any response will be greatly appreciated.

Thanks,
Sounik