LIGGGHTS® - Developer Forum

I have been trying run the command "liggghts -in input.liggghts". However, I'm getting this error:
(1)ERROR: Mesh Silo: The minimum angle of mesh element 0 (line 2) is lower than the specified curvature. Increase mesh quality or decrease curvature (currently 0256235)

(2)ERROR on proc 0: Fix mesh: Bad mesh, cannot continue. You can try setting 'curvature' to 1e-5 or lower or use 'curvature_tolerant yes'

Hi Everyone,

I have been working with LIGGGHTS for the past three months. I am trying to understand how exaclty is the heat flux calculated in the program.

Hey everyone,

I'm a beginner to LIGGGHTS and C++. I'm currently looking to modify LIGGGHTS to compute and store frictional dissipation at the contacts.

Hi !

The last question on the GPU support with LIGGGHTS was in 2020.
I just wanted to know if since there are some updates.
Is LIGGGHTS-Public compatible with GPU or not?
Does the non public version provide the GPU support ?

Thanks in advance.
Best regards,
Ainur

Hello
I am extending LIGGGHTS-flexible fibers version to include nonlinear bond model. To correctly calculate bond energies in bondgran.cpp, I have increase the size of n_granhistory from 14 to 18, which I update in increments at every step.
If I now use some of the old simulations' restarts (based on previous version, i.e, 14 n_granhistory), it throws the following error:
"Íncompatibel restart file: file was created using a bond model with a different number of history values" at line 767 in atom_vec_bond_gran.cpp,

Hello,

I plan to use rolling resistance models.
Studying those models and related codes, I feel that there is a mistake in this header file (ROLLING_MODEL_CDT.H).
I think that for rolling resistance models, wr1 should be sidata.omega_i[0]-sidata.omega_j[0] based on Ai et al. (2011)'s paper instead of using
wr1=(cri * sidata.omega_i[0] + crj * sidata.omega_j[0]) * rinv in the header file (surface_model_default.h).

Could anyone confirm this for me?

Thanks in advance.

Hello All,

I am trying to develop a simple shear test simulation, but I want to set the particles as a wall. How can I implement that using Ligggghts?

See the attached image for clarification.

Best

Saeed

Hi all LIGGGHTS experts,

I am trying to customize LIGGGHTS code such that DEM model can delete particles which has reached a temperature threshold. I am wondering if someone has already succeeded to implement such a temperature threshold on LIGGGHTS code?

Thanks

Hello all,
I am Tarun from India. I wanted to install Aspherix basic in my system. But, while purchasing, I noticed that there is no option for "India" in the "Region" drop-down. Thus, I could not buy it. Can anyone please help me with this???

I really appreciate any help you can provide.

Regards,
Tarun

I am trying to add a spherical harmonics solver model (Keaveny and Maxey) for magnetic force calculation between paramagnetic particles in a uniform magnetic field. Currently, I have the solver written in python using many scipy libraries like linear system solver, combination formula, and associated Legendre functions. Is there a simple way to write these in c++ for LIGGGHTS, or can I include some other libraries like Eigen which may help me?