LIGGGHTS® - User Forum

Hi

Does anyone have any experience with the "lubrication force" between particles? Is the "lubricate" pair_style in LAMMPS/LIGGGHTS the same as what many papers use as the lubrication force? It looks like a totally different formula in the code.

e.g.: Zhang, W., Noda, R., and Horio, M. “Evaluation of Lubrication force on colliding particles for DEM simulation of fluidised beds”. Powder Technology, 158, pp 92-101, 2005.

Thanks
Evan

There was a previous post where I discussed a method of making fairly dense random packings of rigid molecules. Here is the python script I wrote. It is fairly well commented for the procedure.

I have not tested it on LIGGGHTS 3.10 yet, but the output file would work as a restart file for LIGGGHTS 2.3.8.

Hi,

I try to dump out the particles information with

dump id all custom 2000 post/dump.fill id type x y z ....

but it stops writing out dump files over 2.1 GB. The dump files of the meshes are still created but not the particle ones.
what could be the problem?

thanks in advance
Joshua

Hello all,

I'm simulating a particle system with multiple particle sizes (D [mm]: 0.5, 1.0, 2.0, 3.0, 4.0, 5.0). I'm wondering what kind of settings can reduce computational time? Has anyone done some research on this?

Hi all,
I looked in manual and found "partition" command for dividing the total area into multiple partitions, each partition is assigned a different number processors. But I could not exactly understand how? If anyone knows, please tell me how!

Regards
Elham

I have to write a script that generates a cylindrical container (radius 12.5 mm and height 150 mm) filled with a random packing of 500 glass spheres having the same radius (for example 1 mm).
I am interested only in obtaining the spatial coordinates (x,y,z) of the spheres after running the simulation.
How could I write this script?

Hi, All

I am trying to use the servo controlled wall (shear cell example in http://www.cfdem.com/forums/liggghts-tutorial) in LIGGGHTS 3.2.0 (the latest), there is always an error for the below lines (line 66-69), the file is attached:

fix cad3 all mesh/surface/stress/servo file top_plate.stl type 1 scale 0.001 &
stress on com 0 0 0.06 dim z ctrlPV force target_val -10.0 vel_max 1.0
fix geometry2 all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad3

Dear all,
Could someone tell me how to show the width of the lines corresponding to the magnitude of the force chains in paraview? Thanks a lot!

Jian

I am starting out with LIGGGHTS and would appreciate if anyone could take time to review the code that I have written for simulation. There is no error with it. The simulation has been running for last 14 hours and is on 61000 step as of now. I was just wondering if this is efficient way of doing it or it could be made better?

Thanks
--------------------------------------------------------------------------------------------
#################################################

Hi all,

I am going to insert particles in a box as follows, but this take too much time ( about 10 days with 120 cores)!
I would be pleased if anyone knows why!

--------------------------------------------------------------------------my input----------------------------------------------------------------
# insert granular particles into channel

echo both

atom_style granular
atom_modify map array
communicate single vel yes

boundary f f f
newton off
units si