LIGGGHTS® - User Forum

Dear All:

I am trying to run a parallel packing as a simulation starting point just using the examples provided by LIGGGHTS 3.3.1 which lies in:

LIGGGHTS-PUBLIC/examples/LIGGGHTS/Tutorials_public/packing

As I am running on a cluster, I want to utilize all the 16 cores on a node (each node has 16 cores), thus LIGGGHTS 3.3.1 is compiled using Intel Compiler 14.0.2 using OpenMPI 1.6.3 (other versions of MPIs, e.g. mvapich2 are tested as well, with similar results), the only change I made to the in.packing is that I added one line at line 10:

processors 2 2 4

Hi There

It appears that all the links to to the LPP files are broken on git hub - or is this just me ?

Have installed an older version with python 35 and are getting this problem

File "C:/cygwin64/home/Admin/liggghts/liggghts/LPP/src/lpp.py", line 45
raise ValueError, "Invalid or no argument given for chunksize"

Tried manually setting chunksize with no result

Any ideas ??

Best Regards

John

Hi All

New installation on 8 core HP workstation works fine in serial mode with Windows / Cygwin install

In parallel the analysis completes but the program hangs with all 8 processors working.

Only way I can exit is to close down Cygwin and end the processes

The output post files appear fine

Help please !

John

Hello all,

I'm trying to simulate heat transfer in a packed bed. When I reduce the simulation box dimension and particle size, I get the following error:

ERROR on proc 0: Too many neighbor bins

Does anyone know why?

Thanks!
Joel

Hi all,

Having faced the issure of re-assignment of atom IDs for the LIGGGHTS post-processing. As I want to record the position or velocity of some certain particles using
the following commands:
---------------------------------------------------
...
variable x1 equal x[512]
variable y1 equal y[512]
variable z1 equal z[512]
fix pos all print 1000 "${time} ${x1} ${y1} ${z1}" file post/pos_512.txt
....
_____________________________

here 512 is the atom ID (in Paraview as PointID) at time = 0.

Good Day,

I have two different particles, m1 (oranges) and m2 (Apples) in transportation with a conveyor belt (m3). I want to be able to identify m1 and m2 by their type (group) in the simulation. Also, m1 and m2 have different mesh/surface/planar fixes as they are getting into the domain at different positions and times as shown in the attached screenshot. Some sections of the script have:

create_box 3 reg (for m1, m2 and m3)

Hello all,

I am calculating the pair wise contact forces using the compute pair/gran/local command.

If my output of the simulation is : id1 id2 fx fy fz

then, what is the direction of the force between particles of id1 and id2? Is it towards id1's (particles one) centre or id2's (particle two) centre?

I understand fx, fy and fz are the x, y and z components of the force between the particles.

Can someone please help me with this?

Best,

Naren

Hi everyone,

I tried to run my simulations in liggghts and open it in Paraview. I tried to modify the dump file using custom/vtk and I tried to install the liggghts plugin in paraview as well but all my visualization seems to appear in squarish appearance as shown in the image i provided in this forum post.

I heard some used pizza.py and lpp for the visualization but I am unsure on how to use them. I installed pizza.py by using the steps:

Hi all

I am using ParaView to view LIGGGHTS dump files after converting to vtk format. However, the rendering is not so good, as I see only the centres particles. How can I change this, so that I can view complete particles? (I have attached an image showing what I mean)

Thanks,
Joel

Hi all

I was wondering if LIGGGHTS had the feature of a conducting wall, where I could set the wall conductivity. I could only find the option to set wall temperature, but not a conducting wall.

I don't want to measure wall properties. Instead, I just want to change boundary conditions on the particles near the wall, so that, instead of having to interact with a wall at fixed temperature, the particles are now interacting with a wall whose heat flux is given by q=-k(dT/dx)