70 lines
2.4 KiB
ReStructuredText
70 lines
2.4 KiB
ReStructuredText
Granular models
|
|
===============
|
|
|
|
Granular system are composed of spherical particles with a diameter,
|
|
as opposed to point particles. This means they have an angular
|
|
velocity and torque can be imparted to them to cause them to rotate.
|
|
|
|
To run a simulation of a granular model, you will want to use
|
|
the following commands:
|
|
|
|
* :doc:`atom_style sphere <atom_style>`
|
|
* :doc:`fix nve/sphere <fix_nve_sphere>`
|
|
* :doc:`fix gravity <fix_gravity>`
|
|
|
|
This compute
|
|
|
|
* :doc:`compute erotate/sphere <compute_erotate_sphere>`
|
|
|
|
calculates rotational kinetic energy which can be :doc:`output with thermodynamic info <Howto_output>`.
|
|
The compute
|
|
|
|
* :doc:`compute fabric <compute_fabric>`
|
|
|
|
calculates various versions of the fabric tensor for granular and non-granular pair styles.
|
|
|
|
Use one of these 4 pair potentials, which compute forces and torques
|
|
between interacting pairs of particles:
|
|
|
|
* :doc:`pair_style gran/history <pair_gran>`
|
|
* :doc:`pair_style gran/no_history <pair_gran>`
|
|
* :doc:`pair_style gran/hertzian <pair_gran>`
|
|
* :doc:`pair_style granular <pair_granular>`
|
|
|
|
These commands implement fix options specific to granular systems:
|
|
|
|
* :doc:`fix freeze <fix_freeze>`
|
|
* :doc:`fix pour <fix_pour>`
|
|
* :doc:`fix viscous <fix_viscous>`
|
|
* :doc:`fix wall/gran <fix_wall_gran>`
|
|
* :doc:`fix wall/gran/region <fix_wall_gran_region>`
|
|
|
|
The fix style *freeze* zeroes both the force and torque of frozen
|
|
atoms, and should be used for granular system instead of the fix style
|
|
*setforce*\ .
|
|
|
|
To model heat conduction, one must add the temperature and heatflow
|
|
atom variables with:
|
|
* :doc:`fix property/atom <fix_property_atom>`
|
|
a temperature integration fix
|
|
* :doc:`fix heat/flow <fix_heat_flow>`
|
|
and a heat conduction option defined in both
|
|
* :doc:`pair_style granular <pair_granular>`
|
|
* :doc:`fix wall/gran <fix_wall_gran>`
|
|
|
|
For computational efficiency, you can eliminate needless pairwise
|
|
computations between frozen atoms by using this command:
|
|
|
|
* :doc:`neigh_modify <neigh_modify>` exclude
|
|
|
|
.. note::
|
|
|
|
By default, for 2d systems, granular particles are still modeled
|
|
as 3d spheres, not 2d discs (circles), meaning their moment of inertia
|
|
will be the same as in 3d. If you wish to model granular particles in
|
|
2d as 2d discs, see the note on this topic on the :doc:`Howto 2d <Howto_2d>`
|
|
doc page, where 2d simulations are discussed.
|
|
|
|
To add custom granular contact models, see the
|
|
:doc:`modifying granular sub-models page <Modify_gran_sub_mod>`.
|