LIGGGHTS® - User Forum

I am trying to study superquadric particles in a rotating drum.

I am uploading my input file for reference. I want to write dump information for the two kind of particle templates (cube and sphere) in different files. I want to apply different filters on them during visualisation. I am trying to use group command but it isn't working. Please let me know if there is a tutorial on this or if this is possible.

Thanks,
Apoorva

Dear community!

I use the Fix deform variable command to perform controlled compression process. But when I put in equation current value of pressure I get ERROR: Per-atom virial was not tallied on needed timestep.

Does anybody meet the same problem? Any idea to solve this issues?

Below my script...

Im doing tutorial and got an error and not able to solve. here is my code can anyone help me with it.
units si
atom_style sphere
atom_modify map array
boundary f f f
newton off
communicate single vel yes
#processors 2 2 3

# Declare domain
region domain block -0.138 0.138 -0.138 0.138 -0.0045 0.43 units box
create_box 2 domain

# Neighbor listing
neighbor 0.003 bin
neigh_modify every 1 check no

Hi all,
Im currently using LIGGGHTS to simulate the multisphere particle system. However, there are some difficulties in sample generation because of high Young modulus.
My question is: are there any functions similar to nve/limit to restrict the motion of particles within an expected value?
Thank you!
Cheers,
Neiht

Dear all,

I just implemented the Van Der Waals force model on LIGGGHTS. In order to validate my numerical model, I need to know the non-contacting van der Waals force between two interacting particles in LIGGGHTS.

Actually, I tried to compute the aforementioned force using 'compute 1 all pair/gran/local id force' in LIGGGHTS, however, the results only contain some contacting forces.

I would like to know whether LIGGGHTS can output or compute this type of non-contacting force now? If not, I need to write a 'fix' function to write this kind of forces?

Hi, all,
I'm running a simulation with some problems on the mesh. I use a STL file as the wall of my model. I got this STL from the CFD ICEM software. However, the STL exported from this software seems to be bad. The LIGGGHTS gives neither warning or error on this STL file.

The warning is like the following line. It just warns, but the script can run without a hitch. Does it matter? Does this warning have impact on the result?

WARNING: Fix mesh: Mesh contains highly skewed element, moving mesh (if used) will not parallelize well (../surface_mesh_I.h:657)

Hello, all,
I'm trying to use the fix massflow/mesh command to count the particle flowing trough the outlet of my model. I use the following code in my script.
# Mass flow measurement
fix massflow all mesh/surface file Outlet.stl type 3
fix mass all massflow/mesh mesh massflow vec_side 1.0 0.0 0.0 count once file ../DEM/post/Particlemassflow*.txt

Hello
I am trying to print mass of the particle in a particular region but I am unable to get. Please guide me

Hi all,

I already read the documentation about the contact model but I still can't get it.

(https://www.cfdem.com/media/DEM/docu/gran_model_hertz.html)

I understand the equation about the contact force but how can we calculate the normal overlap, tangential overlap,

normal relative vel, tangential relative vel in the LIGGGHTS?

I searched the code 'contact_model_hertz.h' for detail but can't understand how those are calculated.

Can anyone explain this for me?

post deleted