LIGGGHTS® - User Forum

Hi,

I now have output the vtk file for my liggghts code's results. In my code, I have two types of particles with different size, i.e. d=1mm and 2mm. Because I want to see the movement of the particle with d=1mm. Does anyone know how to visulize only particles with d=1mm in PARAVIEW while the other one is invisible? Thanks ahead.

Best,
Rachel

I have to simulate particles with size less than or equal to 1 micron. When I give radius in input script file as 0.000001, I get run time error i.e. volume_limit is exceeded , change the limit. I change the volume_limit in fix_template_sphere .cpp file and build the liggghts again. But it doesn't work out properly while running the script with new volume limit as expected.

Is there any specific for this much particle size?

Govind

Hi there, LIGGGHTS forum. I've been learning to use LIGGGHTS, and I'd like to run a simulation where moving parts stay still and move at specific intervals. Maybe somewhat like the 'pulsed' behavior that's possible by using the duration keyword with the insert_every keyword of the insert/stream fix.

I imagine I could accomplish this by using a repeating pattern of commands like:

run ... upto
fix all spinner move/mesh ...
run ... upto
unfix spinner

Hello,

I have seen similar questions posted from quite a while back (https://www.cfdem.com/forums/restart-heat-transfer-simulation) but still I can't see how to do this successfully.

In the first .in file I use heat/gran to run a heat conduction simulation, with this:

fix heattransfer temp_group heat/gran initial_temperature 200

Hi all,

Is there any command or method to output the cohesion force only?

Kind regards,
Linhan

Hi all,

I want to obtain a vector or array of particle radius.
I use the following commends:

compute rad all property/atom radius
variable ra atom c_rad
variable rb equal ra[1]
if "${rb} >= 0." then "run 1"

But it does not work.

Who can tell me how to use per-atom vectors?

Thank you very much.

Hao XIONG

I am trying to simulate ellipsoidal particles using superquadrics atom_style. Therefore, I started testing with the input file that was released by LIGGGHTS in the example folder (It is attached). For visualization in Paraview, I have compiled it with the Tensorglyph plugins.
I have also attached the video of visualization. I am not able to understand, why after the collision the semi-axes values of both the ellipsoids are changing. For getting a closer look at the problem, I calculated the eigenvalues of both the ellipsoids and found that they are

Hello,

I have a question regarding my build of the 'superquadric' module. I went to the 'LIGGGHTS-PUBLIC/src/MAKE' folder in the 'Makefile.auto' file and modified the USE_SUPERQUADRICS? = "OFF" to "ON". When starting the compilation was warned in the terminal (Ubuntu 16.04 LTS):

Could I simulate the breakage of a distribution of the particles using the LIGGGHTS. Kindly, help me!!

Dear liggghts users,
I want to break multisphere through liggghts. I want to break multispheres when the z coordinate of the multsphere exceeds certain height(here .4 m).
i considered gravity in the +ve z coordinate.
here is my input code.
#Multisphere

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

communicate single vel yes

#PUBLIC version only supports serial execution
processors * * *

units si

region reg block -0.3 0.3 -0.3 0.3 0.0 0.5 units box
create_box 1 reg