\(\renewcommand{\AA}{\text{Å}}\)

fix rheo/pressure command

Syntax

fix ID group-ID rheo/pressure type1 pstyle1 args1 ... typeN pstyleN argsN

  • ID, group-ID are documented in fix command

  • rheo/pressure = style name of this fix command

  • one or more types and pressure styles must be appended

  • types = lists of types (see below)

  • pstyle = linear or taitwater or cubic

    linear args = none
taitwater args = none
cubic args = cubic prefactor \(A_3\) (pressure/density^2)

Examples

fix 1 all rheo/pressure * linear
fix 1 all rheo/pressure 1 linear 2 cubic 10.0

Description

Added in version 29Aug2024.

This fix defines a pressure equation of state for RHEO particles. One can define different equations of state for different atom types. An equation must be specified for every atom type.

One first defines the atom types. A wild-card asterisk can be used in place of or in conjunction with the types argument to set the coefficients for multiple pairs of atom types. This takes the form “*” or “*n” or “m*” or “m*n”. If \(N\) is the number of atom types, then an asterisk with no numeric values means all types from 1 to \(N\). A leading asterisk means all types from 1 to n (inclusive). A trailing asterisk means all types from m to \(N\) (inclusive). A middle asterisk means all types from m to n (inclusive).

The types definition is followed by the pressure style, pstyle. Current options linear, taitwater, and cubic. Style linear is a linear equation of state with a particle pressure \(P\) calculated as

\[P = c (\rho - \rho_0)\]

where \(c\) is the speed of sound, \(\rho_0\) is the equilibrium density, and \(\rho\) is the current density of a particle. The numerical values of \(c\) and \(\rho_0\) are set in fix rheo. Style cubic is a cubic equation of state which has an extra argument \(A_3\),

\[P = c ((\rho - \rho_0) + A_3 (\rho - \rho_0)^3) .\]

Style taitwater is Tait’s equation of state:

\[P = \frac{c^2 \rho_0}{7} \biggl[\left(\frac{\rho}{\rho_0}\right)^{7} - 1\biggr].\]

Restart, fix_modify, output, run start/stop, minimize info

No information about this fix is written to binary restart files. None of the fix_modify options are relevant to this fix. No global or per-atom quantities are stored by this fix for access by various output commands. No parameter of this fix can be used with the start/stop keywords of the run command. This fix is not invoked during energy minimization.

Restrictions

This fix must be used with an atom style that includes density such as atom_style rheo or rheo/thermal. This fix must be used in conjunction with fix rheo. The fix group must be set to all. Only one instance of fix rheo/pressure can be defined.

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

Default

none