<P>keyword = <I>*DENSITY_SUMMATION</I> or <I>*DENSITY_CONTINUITY</I> and <I>*VELOCITY_GRADIENT</I> or <I>*NO_VELOCITY_GRADIENT</I> and <I>*GRADIENT_CORRECTION</I> or <I>*NO_GRADIENT_CORRECTION</I>
<P>The <I>smd/ulsph</I> style computes particle interactions according to continuum mechanics constitutive laws and an updated Lagrangian Smooth-Particle Hydrodynamics algorithm.
</P>
<P>This pair style is invoked similar to the following command:
<P>Here, <I>i</I> and <I>j</I> denote the <I>LAMMPS</I> particle types for which this pair style is
defined. Note that <I>i</I> and <I>j</I> can be different, i.e., <I>ulsph</I> cross interactions
between different particle types are allowed. However, <I>i</I>--<I>i</I> respectively <I>j</I>--<I>j</I> pair_coeff lines have to preceed a cross interaction.
In contrast to the usual <I>LAMMPS</I> <I>pair coeff</I> definitions, which are given solely a
number of floats and integers, the <I>ulsph</I> <I>pair coeff</I> definition is organised using
keywords. These keywords mark the beginning of different sets of parameters for particle properties,
material constitutive models, and damage models. The <I>pair coeff</I> line must be terminated with
the <I>*END</I> keyword. The use the line continuation operator <I>&</I> is recommended. A typical
invocation of the <I>ulsph</I> for a solid body would consist of an equation of state for computing
the pressure (the diagonal components of the stress tensor), and a material model to compute shear
stresses (the off-diagonal components of the stress tensor).
</P>
<P>Note that the use of *GRADIENT_CORRECTION can lead to severe numerical instabilities. For a general fluid simulation, *NO_GRADIENT_CORRECTION is recommended.
</P>
<P>Please see the <A HREF = "USER/smd/SMD_LAMMPS_userguide.pdf">SMD user guide</A> for a complete listing of the possible keywords and material models.