LIGGGHTS® - Developer Forum

Topics related to developing with LIGGGHTS® can be discussed here: discussion about implementation details, C++, MPI and debugging tools

Incompatible restart

Submitted by SHUBHAM AGARWAL on Thu, 05/11/2023 - 09:47

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,

Error in codes (wr1 in rolling_model_cdt.h)

Submitted by min1km on Sat, 04/08/2023 - 05:26

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.

How to delete particle at certain temperature (Tm)

Submitted by emahvelati on Tue, 02/07/2023 - 22:30

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

Aspherix not available for India

Submitted by talk2tarunde@gm... on Mon, 12/19/2022 - 12:13

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

adding new force model | math libraries

Submitted by dusty_sikka on Wed, 07/27/2022 - 22:45

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?

Differences in the stress tensor when using primitive or mesh walls

Submitted by ivanmg94 on Mon, 04/04/2022 - 10:03

Hi, I ran some simulations with dense granular material where the particles are stored in different containers ( cube, cylinder and annular cylinder) and the only external force on them is gravity. I model the particle-particle and the particle-wall contacts with the Hertz contact model with the history model enabled (without adding cohesion nor rolling friction). Every case has been run once using .stl meshes and replicated using primitive walls.

Pages

Subscribe to RSS - LIGGGHTS® - Developer Forum