Syntax:
pair_style style Kn gamma_n xmu dampflag
(mg/d units where m is mass, g is the gravitational constant, d is diameter of a particle)
Examples:
pair_style gran/history 200000.0 0.5 1.0 1
Description:
The gran styles use the following formula (Silbert) for frictional force between two granular particles that are a distance r apart when r is less than the contact distance d.
The 1st term is a normal force and the 2nd term is a tangential force. The other quantites are as follows:
The Kn and gamma_n coefficients are set as parameters to the pair_style command. Xmu is also specified which is the upper limit of the tangential force through the Coulomb criterion Ft = xmu*Fn. The tangential force between 2 particles grows according to a tangential spring and dash-pot model until Ft/Fn = xmu and then is held at Ft = Fn*xmu until the particles lose contact.
For granular styles there are no individual atom type coefficients that can be set via the pair_coeff command. All global settings are made via the pair_style command.
See the citation below for more discussion of granular potentials.
Mixing, shift, table, tail correction, per-atom energy/stress, restart, rRESPA info:
The pair_modify mix, shift, table, and tail options are not relevant for granular pair styles.
These pair styles do not calculate per-atom energy and stress, as used by the compute epair/atom, compute stress/atom, and dump custom commands.
These pair styles write their information to binary restart files, so a pair_style command does not need to be specified in an input script that reads a restart file.
These pair styles can only be used via the pair keyword of the run_style respa command. They do not support the inner, middle, outer keywords of the run_style command.
Restrictions: none
All the granular pair styles are part of the "granular" package. It is only enabled if LAMMPS was built with that package. See the Making LAMMPS section for more info.
Related commands:
Default: none
(Silbert) Silbert, Ertas, Grest, Halsey, Levine, Plimpton, Phys Rev E, 64, p 051302 (2001).