In hydrodynamics, the Perrin friction factors are multiplicative adjustments to the translational and rotational friction of a rigid spheroid, relative to the corresponding frictions in spheres of the same volume. These friction factors were first calculated by Jean-Baptiste Perrin.
These factors pertain to spheroids (i.e., to ellipsoids of revolution), which are characterized by the axial ratio p = (a/b), defined here as the axial semiaxis a(i.e., the semiaxis along the axis of revolution) divided by the equatorial semiaxis b. In prolate spheroids, the axial ratio p > 1 since the axial semiaxis is longer than the equatorial semiaxes. Conversely, in oblate spheroids, the axial ratio p < 1 since the axial semiaxis is shorter than the equatorial semiaxes. Finally, in spheres, the axial ratio p = 1, since all three semiaxes are equal in length.
The formulae presented below assume "stick" (not "slip") boundary conditions, i.e., it is assumed that the velocity of the fluid is zero at the surface of the spheroid.
For brevity in the equations below, we define the Perrin S factor. For prolate spheroids (i.e., cigar-shaped spheroids with two short axes and one long axis)
S \stackrel{def
where the parameter
\xi
\xi \stackrel{def
Similarly, for oblate spheroids (i.e., discus-shaped spheroids with two long axes and one short axis)
S \stackrel{def
For spheres,
S=2
p → 1
The frictional coefficient of an arbitrary spheroid of volume
V
ftot=fsphere fP
where
fsphere
fsphere=6\piηReff=6\piη\left(
3V | |
4\pi |
\right)(1/3)
and
fP
fP \stackrel{def
The frictional coefficient is related to the diffusion constant D by the Einstein relation
D=
kBT | |
ftot |
Hence,
ftot
There are two rotational friction factors for a general spheroid, one for a rotation about the axial semiaxis (denoted
Fax
Feq
Fax \stackrel{def
Feq \stackrel{def
for both prolate and oblate spheroids. For spheres,
Fax=Feq=1
p → 1
These formulae may be numerically unstable when
p ≈ 1
p → 1
1 | |
Fax |
=1.0+\left(
4 | |
5 |
\right)\left(
\xi2 | |
1+\xi2 |
\right)+\left(
4 ⋅ 6 | |
5 ⋅ 7 |
\right)\left(
\xi2 | |
1+\xi2 |
\right)2+\left(
4 ⋅ 6 ⋅ 8 | |
5 ⋅ 7 ⋅ 9 |
\right)\left(
\xi2 | |
1+\xi2 |
\right)3+\ldots
for oblate spheroids.
The rotational friction factors are rarely observed directly. Rather, one measures the exponential rotational relaxation(s) in response to an orienting force (such as flow, applied electric field, etc.). The time constant for relaxation of the axial direction vector is
\tauax=\left(
1 | |
kBT |
\right)
Feq | |
2 |
whereas that for the equatorial direction vectors is
\taueq=\left(
1 | |
kBT |
\right)
FaxFeq | |
Fax+Feq |
These time constants can differ significantly when the axial ratio
\rho
It may seem paradoxical that
\tauax
Feq
\taueq