compute fabric command

Syntax

compute ID group-ID fabric cutoff attribute ... keyword values ...

ID, group-ID are documented in compute command
fabric = style name of this compute command

cutoff = type or radius

type = cutoffs determined based on atom types
radius = cutoffs determined based on atom diameters (atom style sphere)

one or more attributes may be appended

attribute = contact or branch or force/normal or force/tangential

contact = contact tensor
branch = branch tensor
force/normal = normal force tensor
force/tangential = tangential force tensor

zero or more keyword/value pairs may be appended

keyword = type/include

type/include value = arg1 arg2
  arg = separate lists of types (see below)

Examples

compute 1 all fabric type contact force/normal type/include 1,2 3*4
compute 1 all fabric radius force/normal force/tangential

Description

Define a compute that calculates various fabric tensors for pairwise interaction (Ouadfel). Fabric tensors are commonly used to quantify the anisotropy or orientation of granular contacts but can also be used to characterize the direction of pairwise interactions in general systems. The type and radius settings are used to select whether interactions cutoffs are determined by atom types or by the sum of atomic radii (atom style sphere), respectively. Calling this compute is roughly the cost of a pair style invocation as it involves a loop over the neighbor list. If the normal or tangential force tensors are requested, it will be more expensive than a pair style invocation as it will also recalculate all pair forces.

Four fabric tensors are available: the contact, branch, normal force, or tangential force tensor. The contact tensor is calculated as

C_{a b} = \frac{15}{2} (ϕ_{a b} - \frac{1}{3} Tr (ϕ) δ_{a b})

where $a$ and $b$ are the $x$ , $y$ , $z$ directions, $δ_{a b}$ is the Kronecker delta function, and the tensor $ϕ$ is defined as

ϕ_{a b} = \sum_{n = 1}^{N_{p}} \frac{r_{a} r_{b}}{r^{2}}

where $n$ loops over the $N_{p}$ pair interactions in the simulation, $r_{a}$ is the $a$ component of the radial vector between the two pairwise interacting particles, and $r$ is the magnitude of the radial vector.

The branch tensor is calculated as

B_{a b} = \frac{15}{2 Tr (D)} (D_{a b} - \frac{1}{3} Tr (D) δ_{a b})

where the tensor $D$ is defined as

D_{a b} = \sum_{n = 1}^{N_{p}} \frac{1}{N_{c} (r^{2} + C_{c d} r_{c} r_{d})} \frac{r_{a} r_{b}}{r}

where $N_{c}$ is the total number of contacts in the system and the subscripts $c$ and $d$ indices are summed according to Einstein notation.

The normal force fabric tensor is calculated as

F_{a b}^{n} = \frac{15}{2 Tr (N)} (N_{a b} - \frac{1}{3} Tr (N) δ_{a b})

where the tensor $N$ is defined as

N_{a b} = \sum_{n = 1}^{N_{p}} \frac{1}{N_{c} (r^{2} + C_{c d} r_{c} r_{d})} \frac{r_{a} r_{b}}{r^{2}} f_{n}

and $f_{n}$ is the magnitude of the normal, central-body force between the two atoms.

Finally, the tangential force fabric tensor is only defined for pair styles that apply tangential forces to particles, namely granular pair styles. It is calculated as

F_{a b}^{t} = \frac{5}{Tr (N)} (T_{a b} - \frac{1}{3} Tr (T) δ_{a b})

where the tensor $T$ is defined as

T_{a b} = \sum_{n = 1}^{N_{p}} \frac{1}{N_{c} (r^{2} + C_{c d} r_{c} r_{d})} \frac{r_{a} r_{b}}{r^{2}} f_{t}

and $f_{t}$ is the magnitude of the tangential force between the two atoms.

The type/include keyword filters interactions based on the types of the two atoms. Interactions between two atoms are only included in calculations if the atom types are in the two lists. Each list consists of a series of type ranges separated by commas. The range can be specified as a single numeric value, or a wildcard asterisk can be used to specify a range of values. This takes the form “*” or “*n” or “m*” or “m*n”. For example, if $M$ is the number of atom types, then an asterisk with no numeric values means all types from 1 to $M$ . A leading asterisk means all types from 1 to n (inclusive). A trailing asterisk means all types from m to $M$ (inclusive). A middle asterisk means all types from m to n (inclusive). Multiple type/include keywords may be added.

Output info

This compute calculates a global vector of doubles and a global scalar. The vector stores the unique components of the first requested tensor in the order $x x$ , $y y$ , $z z$ , $x y$ , $x z$ , $y z$ followed by the same components for all subsequent tensors. The length of the vector is therefore six times the number of requested tensors. The scalar output is the number of pairwise interactions included in the calculation of the fabric tensor.

Restrictions

This fix is part of the GRANULAR package. It is only enabled if LAMMPS was built with that package. See the Build package doc page for more info.

Currently, compute fabric does not support pair styles with many-body interactions. It also does not support models with long-range Coulombic or dispersion forces, i.e. the kspace_style command in LAMMPS. It also does not support the following fixes which add rigid-body constraints: fix shake, fix rattle, fix rigid, fix rigid/small. It does not support granular pair styles that extend beyond the contact of atomic radii (e.g., JKR and DMT).

Default

none

(Ouadfel) Ouadfel and Rothenburg “Stress-force-fabric relationship for assemblies of ellipsoids”, Mechanics of Materials (2001). (link to paper)