I have three questions about LAMMPS:
First of all, the units used in LAMMPS for granular flow: it is unit mass
or unit density?
In the code, it reads density, then calculate mass and gravity. The data
format in the data file is as follows: granular atom-ID atom-type diameter
density x y z.
For pour command: the default is "The option defaults are diam = 1.0 1.0,
dens = 1.0 1.0, vol = 0.25 50, rate = 0.0, vel = 0.0 0.0 0.0 0.0 0.0. "
Second, Why is there "*PI/6.0" in the follow line?
In the subroutine of fix_gran_diag, we have:
for (m = 0; m < nmax; m++)
fprintf(fpden,"%d %g %g\n",m+1,(m+1)*stepz+boxzlo,dendens[m]*PI/6.0);
Last but not the least, anybody simulate the granular system with liquid
using lattice boltzman method or some other method?
Thank you so much.