LIGGGHTS® - User Forum

liggghts tutorial "heatTransfer_1" set the temperature using the following codes defined to assign 800 degree to the particles at the right half packing.
see the official web: https://www.cfdem.com/heat-transfer

#set particle temperature for half the bed
region halfbed block 0 INF INF INF INF INF units box
set region halfbed property/atom Temp 800.

Dear LIGGGHTS users,

I tried to insert particles into a cubic simulation box with periodic boundary conditions. The particle sizes vary considerably (Dmax/Dmin=80, D: particle diameter). It is OK to construct the initial packing structure with a core, but when the number of particles become very large and need to use multiple cores, some domains are only filled with large particle with no small particles. I have used "communicate multi" and "neighbor multi" command.

Could you please tell me how to solve this problem?
Any help will really be appreciated.

Hello,
(Edit : Moved to liggghts users forum)
My simulations have been hanging when running them on several nodes. They work fine in a single node with several cpus. I was able to get a Minimal Working Example based on the packing example. Just by adding a single line
write_dump all custom out/dump_initial.gz id

Dear LIGGGHTS users,

I want to apply a constant pressure on the top wall and let the bottom wall being frozen.
Could you please tell me the appropriate command to achieve this?

Thank you in advance.

Best wishes,
June

Dear LIGGGHTS uses,

I am now simulating a system with millions of particles . Can GPU acceleration be used for LIGGGHTS?

Thank you very much.

Best wishes,
June

Dear LIGGGHTS users and developers,

I am trying to create an intial packing strucutre of three different size particles: 1, 10, 50. The numbers of particles are 100 (D=1), 10 (D=10), 3(D=50) and the initial volume fraction is 0.15. I have tried two insert command:

1. fix ins all insert/pack seed 32452867 distributiontemplate pdd1 vel contant 0.0 0.0 0.0 insert_every once overlapcheck yes all_in yes
particles_in_region 113 region bc

Hi
I am currently using LIGGGHTS-flexible fiber version.
My problem is somewhat similar to the shear box tutorial. I first create fiber templete by bonding the spheres; bonds and atoms are of type 1.
After fiber -fill in a space and depositing them onto a plate, I want to bond those fibers whose separation is less than a critical value. I want to create a bond of type 2 for inter-fiber bonds.
Since I have atoms of type 1 through out, how should I ensure that atoms are bonded with bond-1 within a fiber, however, with bond-2 inter-fibers (cross-link).

Hello
So far, i am able to generate fibers in random orientation using LIGGGHTS-flexible fibers.
However I wish to generate fibers with a preferred orientation, i.e. fiber crientation should follow a probability distribution of the form: f(theta)= 1/pi -A*cos(2theta). Whare A is a user defined parameter (=0 for random distribution).
Is there a way to implement above in liggghts

Thanks
Shubham

Hello
I am currently using LIGGGHTS-flexible fiber version
For my research needs, I have to make changes in the granular bond model (may be bond_gran2) such that it also incorporates elasto-plastic effects.
Unfortunately, I am not very fluent in C/C++ , however, I can code in python environment.
Is there a way that I can write my custom code code in python and run it with LIGGGHTS?

Hi,

I am performing a biaxial compression simulation where I have a granular array surrounded by 4 extremely large granules to mimic a wall. In this simulation, I start off by fixing the left and bottom "walls" (big granules) and compress the array by adding a force (fix addforce) to move the top and right wall to compress the system. The force I use is to get a desired pressure in the confined state. Up to here all is good. I get a stress state which has the desired pressure. I add viscosity to dampen out the kinetic energy.