Syntax:
pair_style sdpd/taitwater/isothermal temperature viscosity seed
Examples:
pair_style sdpd/taitwater/isothermal 300. 1. 28681 pair_coeff * * 1000.0 1430.0 2.4
Description:
The sdpd/taitwater/isothermal style computes forces between mesoscopic particles according to the Smoothed Dissipative Particle Dynamics model described in this paper by (Espanol and Revenga) under the following assumptions:
The third assumption is true for water in nearly incompressible flows. The fourth holds true for water for any reasonable size one can imagine for a mesoscopic particle.
The pressure forces between particles will be computed according to Tait's equation of state:
where gamma = 7 and B = c_0^2 rho_0 / gamma, with rho_0 being the reference density and c_0 the reference speed of sound.
The laminar viscosity and the random forces will be computed according to formulas described in (Espanol and Revenga).
IMPORTANT NOTE: Similar to brownian and dpd styles, the newton setting for pairwise interactions needs to be on when running LAMMPS in parallel if you want to ensure linear momentum conservation. Otherwise random forces generated for pairs straddling processor boundary will not be equal and opposite.
NOTE: The actual random seed used will be a mix of what you specify and other parameters like the MPI ranks. This is to ensure that different MPI tasks have distinct seeds.
The following coefficients must be defined for each pair of atoms types via the pair_coeff command as in the examples above.
Mixing, shift, table, tail correction, restart, rRESPA info:
This style does not support mixing. Thus, coefficients for all I,J pairs must be specified explicitly.
This style does not support the pair_modify shift, table, and tail options.
This style does not write information to binary restart files. Thus, you need to re-specify the pair_style and pair_coeff commands in an input script that reads a restart file.
This style can only be used via the pair keyword of the run_style respa command. It does not support the inner, middle, outer keywords.
Restrictions:
This pair style is part of the USER-SDPD package. It is only enabled if LAMMPS was built with that package. See the Build package doc page for more info.
Related commands:
Default:
The default seed is 0 (before mixing).
(Espanol and Revenga) Espanol, Revenga, Physical Review E, 67, 026705 (2003).