From b3dacd6244dd77e179a04787cd428aa5943ad23f Mon Sep 17 00:00:00 2001
From: sjplimp <sjplimp@f3b2605a-c512-4ea7-a41b-209d697bcdaa>
Date: Thu, 16 Apr 2009 14:29:04 +0000
Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@2764
 f3b2605a-c512-4ea7-a41b-209d697bcdaa

---
 doc/pair_lubricate.html | 17 ++++++++++-------
 doc/pair_lubricate.txt  | 17 ++++++++++-------
 doc/units.html          | 10 +++++-----
 doc/units.txt           | 10 +++++-----
 4 files changed, 30 insertions(+), 24 deletions(-)

diff --git a/doc/pair_lubricate.html b/doc/pair_lubricate.html
index c7d54cf3ff..0561173318 100644
--- a/doc/pair_lubricate.html
+++ b/doc/pair_lubricate.html
@@ -15,7 +15,7 @@
 </P>
 <PRE>pair_style lubricate mu squeeze shear pump twist cutinner cutoff T_target seed 
 </PRE>
-<UL><LI>mu = viscosity (mass/distance/time units)
+<UL><LI>mu = dynamic viscosity (dynamic viscosity units)
 <LI>squeeze = 0/1 for squeeze force off/on
 <LI>shear = 0/1 for shear force off/on
 <LI>pump = 0/1 for pump force off/on
@@ -39,18 +39,21 @@ spherical particles via this formula from <A HREF = "#Ball">(Ball and Melrose)</
 <CENTER><IMG SRC = "Eqs/pair_lubricate.jpg">
 </CENTER>
 <P>which represents the dissipation W between two nearby particles due to
-their relative velocities in the presence of a background solvent.  Rc
-is the outer cutoff specified in the pair_style command, the
+their relative velocities in the presence of a background solvent with
+viscosity mu.  Note that this is dynamic viscosity which has units of
+mass/distance/time, not kinematic viscosity.
+</P>
+<P>Rc is the outer cutoff specified in the pair_style command, the
 translational velocities of the 2 particles are v1 and v2, the angular
 velocities are w1 and w2, and n is the unit vector in the direction
 from particle 1 to 2.  The 4 terms represent four modes of pairwise
 interaction: squeezing, shearing, pumping, and twisting.  The 4 flags
 in the pair_style command turn on or off each of these modes by
 including or excluding each term.  The 4 coefficients on each term are
-functions of the separation distance of the particles.  Details are
-given in <A HREF = "#Ball">(Ball and Melrose)</A>, including the forces and torques
-that result from taking derivatives of this equation with respect to
-velocity (see Appendix A).
+functions of the separation distance of the particles and the
+viscosity.  Details are given in <A HREF = "#Ball">(Ball and Melrose)</A>, including
+the forces and torques that result from taking derivatives of this
+equation with respect to velocity (see Appendix A).
 </P>
 <P>Unlike most pair potentials, the two specified cutoffs (cutinner and
 cutoff) refer to the surface-to-surface separation between two
diff --git a/doc/pair_lubricate.txt b/doc/pair_lubricate.txt
index 7b8438f6eb..c7b64d7adc 100644
--- a/doc/pair_lubricate.txt
+++ b/doc/pair_lubricate.txt
@@ -12,7 +12,7 @@ pair_style lubricate command :h3
 
 pair_style lubricate mu squeeze shear pump twist cutinner cutoff T_target seed :pre
 
-mu = viscosity (mass/distance/time units)
+mu = dynamic viscosity (dynamic viscosity units)
 squeeze = 0/1 for squeeze force off/on
 shear = 0/1 for shear force off/on
 pump = 0/1 for pump force off/on
@@ -36,18 +36,21 @@ spherical particles via this formula from "(Ball and Melrose)"_#Ball
 :c,image(Eqs/pair_lubricate.jpg)
 
 which represents the dissipation W between two nearby particles due to
-their relative velocities in the presence of a background solvent.  Rc
-is the outer cutoff specified in the pair_style command, the
+their relative velocities in the presence of a background solvent with
+viscosity mu.  Note that this is dynamic viscosity which has units of
+mass/distance/time, not kinematic viscosity.
+
+Rc is the outer cutoff specified in the pair_style command, the
 translational velocities of the 2 particles are v1 and v2, the angular
 velocities are w1 and w2, and n is the unit vector in the direction
 from particle 1 to 2.  The 4 terms represent four modes of pairwise
 interaction: squeezing, shearing, pumping, and twisting.  The 4 flags
 in the pair_style command turn on or off each of these modes by
 including or excluding each term.  The 4 coefficients on each term are
-functions of the separation distance of the particles.  Details are
-given in "(Ball and Melrose)"_#Ball, including the forces and torques
-that result from taking derivatives of this equation with respect to
-velocity (see Appendix A).
+functions of the separation distance of the particles and the
+viscosity.  Details are given in "(Ball and Melrose)"_#Ball, including
+the forces and torques that result from taking derivatives of this
+equation with respect to velocity (see Appendix A).
 
 Unlike most pair potentials, the two specified cutoffs (cutinner and
 cutoff) refer to the surface-to-surface separation between two
diff --git a/doc/units.html b/doc/units.html
index 2d9c086aa6..950d235e99 100644
--- a/doc/units.html
+++ b/doc/units.html
@@ -52,7 +52,7 @@ results from a unitless LJ simulation into physical quantities.
 <LI>torque = epsilon, where t* = t / epsilon
 <LI>temperature = reduced LJ temperature, where T* = T Kb / epsilon
 <LI>pressure = reduced LJ pressure, where P* = P sigma^3 / epsilon
-<LI>viscosity = reduced LJ viscosity, where eta* = eta sigma^3 / epsilon / tau
+<LI>dynamic viscosity = reduced LJ viscosity, where eta* = eta sigma^3 / epsilon / tau
 <LI>charge = reduced LJ charge, where q* = q / (4 pi perm0 sigma epsilon)^1/2
 <LI>dipole = reduced LJ dipole, moment where *mu = mu / (4 pi perm0 sigma^3 epsilon)^1/2
 <LI>electric field = force/charge, where E* = E (4 pi perm0 sigma epsilon)^1/2 sigma / epsilon 
@@ -68,7 +68,7 @@ results from a unitless LJ simulation into physical quantities.
 <LI>torque = Kcal/mole
 <LI>temperature = degrees K
 <LI>pressure = atmospheres
-<LI>viscosity = Poise
+<LI>dynamic viscosity = Poise
 <LI>charge = multiple of electron charge (+1.0 is a proton)
 <LI>dipole = charge*Angstroms
 <LI>electric field = volts/Angstrom 
@@ -84,7 +84,7 @@ results from a unitless LJ simulation into physical quantities.
 <LI>torque = eV
 <LI>temperature = degrees K
 <LI>pressure = bars
-<LI>viscosity = Poise
+<LI>dynamic viscosity = Poise
 <LI>charge = multiple of electron charge (+1.0 is a proton)
 <LI>dipole = charge*Angstroms
 <LI>electric field = volts/Angstrom 
@@ -100,7 +100,7 @@ results from a unitless LJ simulation into physical quantities.
 <LI>torque = Newton-meters
 <LI>temperature = degrees K
 <LI>pressure = Pascals
-<LI>viscosity = Pascal*second
+<LI>dynamic viscosity = Pascal*second
 <LI>charge = Coulombs
 <LI>dipole = Coulombs*meters
 <LI>electric field = volts/meter 
@@ -116,7 +116,7 @@ results from a unitless LJ simulation into physical quantities.
 <LI>torque = dyne-centimeters
 <LI>temperature = degrees K
 <LI>pressure = dyne/cm^2 or barye = 1.0e-6 bars
-<LI>viscosity = Poise
+<LI>dynamic viscosity = Poise
 <LI>charge = statcoulombs or esu
 <LI>dipole = statcoul-cm = 10^18 debye
 <LI>electric field = statvolt/cm or dyne/esu 
diff --git a/doc/units.txt b/doc/units.txt
index 0aebddeba8..f1f7888940 100644
--- a/doc/units.txt
+++ b/doc/units.txt
@@ -49,7 +49,7 @@ force = epsilon/sigma, where f* = f sigma / epsilon
 torque = epsilon, where t* = t / epsilon
 temperature = reduced LJ temperature, where T* = T Kb / epsilon
 pressure = reduced LJ pressure, where P* = P sigma^3 / epsilon
-viscosity = reduced LJ viscosity, where eta* = eta sigma^3 / epsilon / tau
+dynamic viscosity = reduced LJ viscosity, where eta* = eta sigma^3 / epsilon / tau
 charge = reduced LJ charge, where q* = q / (4 pi perm0 sigma epsilon)^1/2
 dipole = reduced LJ dipole, moment where *mu = mu / (4 pi perm0 sigma^3 epsilon)^1/2
 electric field = force/charge, where E* = E (4 pi perm0 sigma epsilon)^1/2 sigma / epsilon :ul
@@ -65,7 +65,7 @@ force = Kcal/mole-Angstrom
 torque = Kcal/mole
 temperature = degrees K
 pressure = atmospheres
-viscosity = Poise
+dynamic viscosity = Poise
 charge = multiple of electron charge (+1.0 is a proton)
 dipole = charge*Angstroms
 electric field = volts/Angstrom :ul
@@ -81,7 +81,7 @@ force = eV/Angstrom
 torque = eV
 temperature = degrees K
 pressure = bars
-viscosity = Poise
+dynamic viscosity = Poise
 charge = multiple of electron charge (+1.0 is a proton)
 dipole = charge*Angstroms
 electric field = volts/Angstrom :ul
@@ -97,7 +97,7 @@ force = Newtons
 torque = Newton-meters
 temperature = degrees K
 pressure = Pascals
-viscosity = Pascal*second
+dynamic viscosity = Pascal*second
 charge = Coulombs
 dipole = Coulombs*meters
 electric field = volts/meter :ul
@@ -113,7 +113,7 @@ force = dynes
 torque = dyne-centimeters
 temperature = degrees K
 pressure = dyne/cm^2 or barye = 1.0e-6 bars
-viscosity = Poise
+dynamic viscosity = Poise
 charge = statcoulombs or esu
 dipole = statcoul-cm = 10^18 debye
 electric field = statvolt/cm or dyne/esu :ul