LIGGGHTS® - User Forum

Hi, All:

I was trying to use the dynamic load balancing of LIGGHTS to minimize the computational time. I found the following document: http://cfdem.dcs-computing.com/media/DEM/benchmarks/LIGGGHTS_Benchmarks.pdf, the LIGGGHTS version I am using is: Version LIGGGHTS-PUBLIC 3.0.6

However if I try to use the load balancing command (e.g. in page 43 of the document) :

fix lb all balance 0 xyz 10 1.1 out balance.txt

The program prompts the error:

Hello, can you give me any info about how to use multiprocessing in Liggghts 3.0.6? i mean, usually in 3.0.2 i was using mpirun for 4 my procs and it worked well! and after the installation of the new version, it does not work anymore... is something changed? Thanks for your answers!

Dear all,

I am simulating two spheres (sphere1 and sphere2) which are dropped in rotating cylinder.

I'm looking for a command that give me the total number of interaction that sphere1 participates(total number of interact with wall and with sphere 2) .

I tried the commands (compute wall/gran/local id …..and pair/gran/local id and then dump local), but I got the IDs of particles in contact at each timesteps.However I am looking for the total number of interactions.

Thank you in advance.

Best,

Abdul

Hello everyone,

I'm trying to pack rashig rings to create a packed bed reactor.
I don't know from where I have to start. I have an STL of the ring, I think I have to create a clump of that shape and then use the multisphere particle template. Is it right?
I'm not able to create the dump file from the STL, can anyone help me?

Any hint or suggestion is very helpful

Thanks

Hi all,

As we all know, neighbor cutoff = force cutoff+skin. However, I don't know how LIGGGHTS determines the force cutoff value.

Are there any command in LIGGGHTS that we could use to modify force cutoff (haven't found anyone yet)?

Best regards,

Su

Hi,guys

I want to calculate the average contactarea between particles,below is my command:

compute area particle pair/gran/local contactArea
variable avearea equal ave(c_area)
thermo_style custom step atoms ke v_avearea

but, it pointed out that mismatched compute in variable formula.........

I don't know why , who can tell me ?

thanks

hi,

if I use the command

variable zero equal 0
variable an equal -0.1
fix movement all move/mesh mesh mole linear/variable v_zero v_zero v_an
run 10000

and try to change now the vertical velocity to

variable an equal -10

the velocity to the mesh is set positiv doesnt mind if I set it to 10 or -10.
To test it just create a mesh and do

variable zero equal 0
variable an equal -0.1
fix movement all move/mesh mesh mole linear/variable v_zero v_zero v_an
run 10000
variable an equal -10

Hi all,

When unfixing a move wiggle fixed for a substrate made of particles, the kinetic energy of these particles wouldn't go to zero unless zeroing velocity by hand, i.e., by using "velocity group zero linear". Any idea?

Selim,

Hi all,

I am trying to understand the 6 components of stress tensor computed by stress/atom.
What are the physical meaning of kinetic energy contribution, angle, dihedral and improper component of the stress tensor?

What should be the keyword listed if I want to compute the stress tensor of a particle packing (sand foundation) under the gravitational force field?

Many thanks

Hi,

I'm trying to understand the values that are calculated with:

dump custom fx fy fz
compute pair/gran/local force
compute stress/atom

therefor I insert 2 particles, but just one particle is integrated. One particle lies on top of the other and gravity loading is on both of them acting.