Einstein viscosity formula for a dilute suspension Created 2026-10-07 Updated 2026-10-14
A dilute suspension of identical rigid spheres raises the dynamic viscosity by a fraction , where is the particle volume fraction. At imposed trace-free rate-of-strain tensor , a sphere of radius adds viscous dissipation . Multiplication by the number density gives , which added to proves the formula. This is the leading noninteracting-sphere term; hydrodynamic interactions enter at higher concentration.
In the free-draining approximation, each separated component of a body experiences its isolated viscous drag, and hydrodynamic interactions between components are neglected. For widely separated spheres of radii in an unbounded Newtonian fluid of dynamic viscosity , Stokes law gives , with . This is a leading approximation in sphere size divided by separation. Any connecting rod must also have negligible drag; a small diameter alone does not guarantee this for an arbitrarily long rod.
Hookean dumbbell model 2026-10-24
Represent a dilute polymer by two beads connected by a spring obeying Hooke's law, with stiffness and bead drag coefficient . Neglect bead inertia and hydrodynamic interactions. For component-first velocity gradient , the configuration Fokker-Planck probability current is
The factor two comes from relative motion of both beads. At rest its normalized probability density is . The conformation tensor has dimensions of length squared; is its dimensionless version. The Hookean dumbbell in planar extension is a particular imposed flow in this model.
Hydrodynamic mobility matrix Created 2026-10-06 Updated 2026-10-25
With a consistent force convention, the mobility matrix maps applied generalized forces to rigid-body velocities; it is the inverse of the hydrodynamic resistance matrix. For a sphere with force on the fluid , its isolated mobility is . Coupling between bodies represents hydrodynamic interactions.
Oseen tensor Created 2026-10-20 Updated 2026-10-25
The Oseen tensor is the unbounded three-dimensional incompressible Stokes flow Green tensor relating a point force to the velocity it produces: . At separation , its longitudinal mobility is twice its transverse mobility. It generates the leading far-field hydrodynamic interaction between externally forced particles. The singular value at is not a particle's self mobility; a no-slip sphere instead has self mobility . Walls and near-contact lubrication require different or corrected tensors.
Use the deterministic Hookean dumbbell model, with a zero-rest-length spring and negligible connector drag. For , its elastic forces on the two beads are
Here is the spring constant. No independent thermal forcing from Brownian motion or bead inertia is included. The unperturbed linear velocity field is , with by incompressibility.
The fluid velocity used in each bead's drag law is the imposed flow plus the disturbance from the other bead, with its own singular field omitted. At leading order in , the evenness of the Oseen tensor, , gives
These are the pair hydrodynamic mobility matrix relations. Their signs follow from balancing each spring force with , the fluid drag on that bead.
Adding the equations cancels both the direct spring slips and the cross-induced velocities, since . For the midpoint ,
The midpoint therefore follows the undisturbed linear velocity field. Subtracting the equations instead gives the hydrodynamic interaction of a Hookean dumbbell:
The dimensionless interaction matrix is
Its transverse and longitudinal eigenvalues are and , respectively. Since the spring force is longitudinal, the separation equation can also be written
The relative correction to the free-draining spring relaxation rate is therefore , explicitly proportional to . Entrainment by the other bead reduces their relative approach for given opposing spring forces.
This is a far-field hydrodynamic interaction, valid for ; it is not a near-contact law and omits higher reflections and lubrication. The printed target also fixes the spring's rest length implicitly. For an ordinary spring obeying Hooke's law with nonzero rest length , the elastic force would be , and every spring term above would acquire the factor . The displayed result uses the zero-rest-length spring convention appropriate to the ideal Hookean dumbbell model.
Take positive upward, let sphere A's centre be , and write its downward velocity as , with . At B's centre the displacement from A is . The vorticity of translating-sphere Stokes flow therefore gives
B is held in translation but has zero applied torque. By Faxén's rotational law, it rotates with half the incident vorticity:
Integrating the rotation induced by a passing Stokes sphere along the complete path cancels :
Its magnitude is radians, with direction about the negative axis by the right-hand rule. All instantaneous rotations have this same axis, so the integrated angular velocity gives the net angle without a noncommuting-rotation correction. The result is independent of B's radius in the prescribed incident-flow calculation.
Figure 1.
Angular velocity and cumulative rotation from an isolated sphere passing a torque-free target
. Details.
The inequality ensures no overlap. It does not by itself make mutual hydrodynamic interactions negligible: the stated flow is the isolated-A solution, and a held finite B also disturbs A's flow. Thus the displayed angle is the result for that supplied incident flow, or the leading one-way approximation to separated bodies. A fully coupled two-sphere calculation uses the method of reflections for Stokes flow or an exact pair resistance solution.
Use the free-draining approximation: neglect the hydrodynamic interaction between the widely separated spheres and neglect the connector's viscous drag. Stokes law gives for each sphere. Their resultant is
This is the leading two-sphere approximation, not an exact finite-separation hydrodynamic resistance matrix. Writing and , the first Stokeslet reflection gives, if the connector remains negligible,
The negative correction expresses reduced drag from fluid entrainment by the other sphere. A very thin but arbitrarily long rod need not have negligible total drag: its typical transverse drag scales as for rod radius . Thus the simple boxed result uses the intended negligibly dragging connector limit as well as ; the two printed size inequalities alone do not specify a uniform limit in which rod drag vanishes.
Hold the two-sphere object fixed in a uniform stream . In the free-draining approximation, its drag forces are , with . The hydrodynamic torque about is
To make it zero for every incoming-stream direction, choose the hydrodynamic centre of a two-sphere dumbbell:
It lies closer to the larger sphere, whose larger drag requires a shorter lever arm. Uniform translation does not spin an isolated sphere, so no additional isolated-sphere spin torque enters this leading approximation.
Figure 1.
Drag-weighted hydrodynamic centre of two widely separated spheres
. Details.
Finite-separation hydrodynamic interactions and connector drag can change the displayed location; the calculation has the same approximation as part c.
Put and . In the Unscaled Papkovich–Neuber representation, take the harmonic potentials and . Differentiating gives the Stokeslet velocity and its pressure:
Indeed, writing ,
The antisymmetric first two terms cancel in the rate-of-strain tensor, whereas their curl supplies the vorticity. Thus
The rate-of-strain tensor has zero trace, as required by incompressibility. The stress is ; integrating its traction with normal over a sphere gives . This fixes the force normalization of the Stokeslet.
For a sphere in a uniform straining Stokes flow, write the ambient velocity as , where is symmetric and trace-free, and set . Symmetry gives zero translation and rotation of the centred force-free, torque-free sphere. Seek a decaying disturbance using the Unscaled Papkovich–Neuber representation with
The components of are derivatives of , and is a trace-free harmonic quadrupole, so both potentials satisfy Laplace's equation outside the sphere. Direct differentiation yields
The no-slip boundary condition at requires the coefficient of to be and that of to vanish: , . Hence , , and
This disturbance cancels the ambient velocity at , decays at infinity, and has a leading stresslet proportional to .
For the two-sphere hydrodynamic interaction, let point from the forced sphere to the second sphere. The first sphere's leading far field is the Stokeslet with . Its incident rate-of-strain tensor at the second sphere is
The force-free sphere translates with the local uniform incident flow, and the torque-free sphere rotates with half the incident vorticity. Those parts cause no leading Stokeslet or rotlet disturbance. The remaining incident strain produces the stresslet just calculated. At displacement from the second sphere its velocity is
The first sphere translates in this returned incident flow. Finite-size sampling and higher incident multipoles enter only beyond the retained reflection order. Thus the stresslet reflection between two force-free and forced spheres gives
The factor makes explicit the velocity scale implicit in the printed remainder.
The resulting input power, equal to the viscous dissipation, is
It decreases at fixed force. The Minimum-dissipation theorem for Stokes flow instead compares admissible fields at fixed first-sphere velocity: extending a two-sphere field rigidly into the second sphere gives an admissible one-sphere field, so the extra rigid inclusion can only increase the minimum resistance. Allowing the second sphere to choose its translation and rotation minimizes over its rigid motions, producing zero force and couple there. An increase of this effective hydrodynamic resistance matrix decreases its inverse hydrodynamic mobility matrix, and hence decreases fixed-force power. The negative mobility correction is therefore consistent with minimum dissipation.
If the second sphere is prevented from rotating, its angular velocity relative to the incident fluid becomes
Use the rotating sphere in Stokes flow at . The additional reflected rotlet gives
Adding this rotation-clamped reflection between two spheres to the strain reflection,
The holding couple supplies no power because the held sphere's angular velocity is zero.
For a Newtonian fluid, write , where is the rate-of-strain tensor. The Stokes equation without body force gives . Since the stress tensor is symmetric, its contraction with the antisymmetric part of the velocity gradient is zero. Thus
The divergence theorem proves the boundary-work identity
For an incompressible flow, , so is precisely the viscous dissipation.
Apply this identity to the particle-free volume and then to the fluid outside the added particle. On the rigid particle, the velocity has the form . Its contribution to the boundary work is , which vanishes because the particle is force-free and torque-free. The outer velocity is unchanged. Subtracting the two boundary-work identities therefore gives
To move this expression to the particle surface, apply the Lorentz reciprocal theorem to and in the punctured fluid domain. For completeness, its local proof is
On the outer boundary . Let be the particle-outward normal, opposite to the fluid-outward normal there. Integrating the preceding divergence identity gives
This expression includes the removal of the original fluid inside the particle; that removal must not be counted separately again.
For a sphere in a uniform straining Stokes flow, is symmetric and trace free, and . Symmetry gives zero rigid translation and rotation. Put and . Use the Unscaled Papkovich–Neuber representation printed on the paper, with the decaying harmonic potentials
Their harmonicity follows from and from being a trace-free second derivative of . Direct differentiation gives
At , the no-slip boundary condition requires . Matching the independent and terms fixes and . The disturbance velocity and pressure are consequently
The pressure follows from . The disturbance decays, the total velocity vanishes at the sphere, and the harmonic representation verifies the Stokes equation. The boundary traction has zero total force by odd symmetry, and zero total torque because is symmetric, confirming the assumed force-free and torque-free motion.
At the sphere , so the extra dissipation due to a rigid inclusion becomes . Contracting the supplied total stress tensor with cancels its terms proportional to , leaving . Thus
Using proves the single-sphere dissipation increase
One may perform the comparison in a large bounded vessel with prescribed straining velocity and then take its size to infinity; the given near-sphere approximation makes this limit compatible with the fixed-outer-velocity definition of .
The sphere number density is . Neglecting hydrodynamic interactions, the additional viscous dissipation per unit volume is . Adding the particle-free value gives
This is the Einstein viscosity formula for a dilute suspension, to first order in the particle volume fraction.
For an incompressible flow, write the Newtonian fluid stress tensor as , where is the rate-of-strain tensor. The Stokes equation gives . Symmetry of the Cauchy stress tensor therefore gives
Integrating and applying the divergence theorem expresses the viscous dissipation as boundary power:
At a moving rigid body, the surface velocity is . Its boundary-power contribution is , with the force and torque evaluated using the fluid's outward normal vector. Both resultants vanish for a force-free, torque-free inclusion. Thus the new viscous dissipation comes entirely from the unchanged outer boundary velocity. Subtracting the particle-free boundary power gives
Apply the Lorentz reciprocal theorem to and in the fluid outside the inclusion. On the outer boundary , so
Writing for the particle's outward normal vector gives the required extra dissipation due to a rigid inclusion:
This subtraction already accounts for the fluid volume displaced by the inclusion; it is not just the integral of the disturbance's local viscous dissipation.
For the sphere, the ambient rate-of-strain tensor is symmetric and trace free. The sphere in a uniform straining Stokes flow has zero translational and rotational velocity: inversion symmetry eliminates its force, and symmetry of eliminates its torque. To derive the disturbance, put , and seek
The incompressible flow condition and Stokes equation reduce to
A decaying family satisfying these equations is , , . The no-slip boundary condition requires and , so and . Consequently
The total velocity is zero at , and the disturbance decays as . The pressure constant has been set to zero. These boundary and far-field conditions, together with Uniqueness of Stokes flow, establish the solution. There is no rotational background in a pure strain flow.
On the sphere . Contracting the supplied Newtonian fluid stress tensor with makes its two terms proportional to cancel, leaving . Hence the extra dissipation due to a rigid inclusion is
Here . Taking the large outer boundary limit only after using the fixed-boundary power identity avoids replacing that prescribed boundary by an uncontrolled boundary at infinity.
If denotes the number of spheres per unit volume, . Thus
Add this to the particle-free viscous dissipation density . Comparing with gives the Einstein viscosity formula for a dilute suspension:
to first order in . The hydrodynamic interactions neglected here contribute beyond that dilute order.
In resistive-force theory (RFT), the hydrodynamic force per unit arclength exerted by the fluid on a slender filament is
Here is the unit tangent and the local velocity relative to the background fluid. The parallel and perpendicular drag coefficients of a slender filament are positive and generally satisfy ; their leading logarithmic ratio is about two. This local approximation assumes a small filament radius compared with length and curvature scales, negligible fluid inertia, and a Newtonian fluid. It represents drag by the local tangent direction and neglects nonlocal hydrodynamic interactions between separated filament segments. Boundaries, close approaches and end corrections may require slender-body theory or a more complete flow calculation. In this problem the background fluid is at rest, and the drag coefficients are taken uniform along the filament.
The rigid-body velocity in a deforming swimmer frame is the sum of translation, rotation and material deformation. With ,
so, in instantaneous swimmer-frame components,
The reference-frame conditions attach the origin and orientation to the filament's base and base tangent. They prevent arbitrary shape translations or tilts from being absorbed into the definition of .
The exact tangent and arclength element are
Since , the local velocity is . Changing the drag tensor by its tangent correction therefore changes only at . The arclength correction is smaller still at this order. Also . Thus is sufficient at order . This assumes the small-slope expansion is uniform, .
Put and , with . The first-order local force is
Its total hydrodynamic forces are
Taking the moment about the swimmer-frame origin, , the second term is beyond first order. Hence
For force-free and torque-free motion, the three leading coefficients vanish. The longitudinal force immediately gives . The two transverse equations are
Their determinant is . Solving yields the first-order free swimming of a planar filament:
The resistive-force theory coefficient cancels because both equations use the same transverse drag. Another interpretation is the least-squares projection of filament deformation velocity: is the best affine approximation to on . Zero transverse force and torque mean that the residual is orthogonal to and .
If is sufficiently differentiable and periodic with period , each has zero temporal mean. More explicitly,
Both bracketed quantities return to their initial values over a period. Thus the first-order periodic transverse swimming velocity satisfies
These results concern first order and swimmer-frame components. Changes of orientation can affect laboratory displacements at second order; higher-order net swimming is not ruled out.
Use the same resistive-force theory convention as before: the force density is the fluid's force on the filament. Let be the usual cylindrical azimuth. Choose the orientation of a left-handed helix so that, as increases, decreases. With , ,
Reversing the orientation of the tangent leaves the local drag tensor unchanged. Define . From resistive-force theory, the relevant force-density components are
These scalar axial and azimuthal components are uniform along the helix. Integration over its contour length , and use of the axial torque density , give the axial resistance matrix of a slender helix:
Thus , as also required by the Lorentz reciprocal theorem. The signs here belong to , with forces and moments on the helix. The handedness reversal of helical hydrodynamic resistance reverses but leaves unchanged. If the positive rotation direction or handedness convention is reversed, the coupling sign reverses with it.
As a check on physical admissibility, the determinant of helical resistance in resistive-force theory is
Together with , this makes the hydrodynamic resistance matrix positive definite and the viscous power loss positive. Isotropic local drag would have , so rotating a helix would not propel it in this approximation.
For the pair, use a local additive resistive-force theory model: both helices have the same drag coefficients, and inter-helix hydrodynamic interactions and unresolved end effects are omitted. The right-handed helix has resistance entries . Define the signed motor motion by
The printed relative rotation specifies a magnitude rather than which helix rotates relative to which. Choosing instead reverses all signed speeds below.
The common translation speed and zero total force and torque obey
Internal motor forces and torques cancel in these totals. Introduce
Solving the three linear equations yields the opposite-handed counterrotating helical swimmer:
For , and , , so , and . The opposite-handed helices counterrotate but contribute thrust in the same axial direction. The formula satisfies the specified relative rotation without identifying either laboratory rotation rate with the motor rate.
At , the second helix supplies no hydrodynamic resistance. The remaining helix must have , and invertibility of its hydrodynamic resistance matrix forces . Formally is the rotation of a zero-resistance motor shaft or vanishing second rotor; there is no physical finite second helix to propel or to provide a reaction torque. This is the vanishing reaction rotor in a helical swimmer.
At , equal-length opposite-handed helices have
The translation-generated axial torques cancel between the two helices, and equal counterrotation balances the rotational torques. Their propulsive forces add, rather than cancel, because handedness and rotation both reverse. These are equal-length opposite-handed helices.
For , the large reaction helix limit is
The increasingly long second helix nearly anchors the whole assembly: its very small translation and rotation suffice to balance the finite force and torque generated by the first helix. The first helix rotates at almost the full motor rate, but moving the large resistive second helix makes the common translation tend to zero. These are idealized limits of the additive local model, rather than a resolution of short-helix end effects or infinitely extended interacting filaments.
Figure 1.
Counterrotation rates and common translation versus the right-to-left helix length ratio, showing no propulsion at zero or infinite ratio and symmetric counterrotation for equal lengths
. Details.
In the imposed free-draining polymer friction model, the drag coefficients of the monomers add, giving and . Taking the same three-dimensional centre-of-mass RMS convention as part (ii),
The polymer is a one-dimensional connected contour, but that alone does not restrict its centre of mass to one spatial axis. If the intended coordinate is instead literal one-dimensional curvilinear motion along an entanglement tube, its RMS contour displacement is . That is a different coordinate and should not be compared with a three-dimensional displacement without stating the conversion/model.
Additive local drag is the free-draining Rouse model approximation. It neglects long-range solvent-mediated hydrodynamic interactions and is not a complete theory of an entangled chain. In a tube/reptation description, curvilinear diffusion can scale as while long-time spatial centre-of-mass diffusion scales as after the chain renews its tube. Conversely a dilute unentangled coil can exhibit Zimm model hydrodynamic coupling. Monomer motion, centre-of-mass motion and contour sliding therefore cannot all be identified with the same displacement law.
In a dilute solution, solvent-mediated hydrodynamic interactions couple distant chain segments. The Zimm model describes the resulting cooperative dynamics: an intermediate-time segment mean-square displacement scales as , while whole-chain translation has . Once internal modes relax, each tagged segment follows the chain's ordinary long-time diffusion, with mean-square displacement asymptotic to .
In a concentrated solution or polymer melt, hydrodynamic coupling is screened. The unentangled Rouse model has an intermediate segment regime and whole-chain diffusivity proportional to . Long chains with polymer entanglement are additionally confined by neighbouring chains: reptation gives a slower intermediate regime and an asymptotic chain diffusivity scaling approximately as . The time windows depend on segment, entanglement and terminal relaxation times, so these powers do not all describe the same time interval.
After sufficiently long times both dilute and concentrated chains undergo ordinary translational diffusion. The main long-time difference is the coefficient and its chain-length dependence, with entanglement strongly reducing mobility when the chains are long enough; concentration alone does not guarantee an entangled regime.
In a sufficiently uniform flow the Stokes force is toward trap 2. Its work over separation lowers the second effective well by . Equal effective free energies require . Hence
at 298 K. The temperature was not specified, so the numerical result scales as . This approximation treats a weak uniform drag as a constant tilt and neglects changes in barrier shape and hydrodynamic interactions with boundaries.
A moving or externally forced colloid disturbs the surrounding fluid, changing the velocity of other colloids even without direct contact. These hydrodynamic interactions are long-ranged at small Reynolds number. In an unbounded incompressible Newtonian liquid, the Stokes flow response to a point force is the Stokeslet:
Its Oseen tensor is
This is a mobility tensor, with units velocity per force, rather than a stress tensor. At separation the response parallel to the line of centres is ; the perpendicular response is . Flow is not isotropic even though the solvent is. The disturbance decays as , so treating all drag forces as independent can miss correlated motion.
For well-separated equal spheres of radius , the leading hydrodynamic mobility matrix gives
The first term is the isolated-sphere self mobility; the remaining terms couple distinct particles. At thermal equilibrium the fluctuation-dissipation theorem also predicts correlated short-time increments, with the configuration locally fixed. Thus correlated Brownian displacements can reveal hydrodynamic coupling without an applied force.
The point-force approximation requires and negligible inertia. It diverges at the origin and cannot replace the finite self mobility. Near-contact lubrication, walls, many-body corrections, and a viscoelastic solvent require more complete mobility laws. A force-free particle can instead have a faster-decaying leading disturbance, so the claim specifically concerns a net-force-driven object.
Figure 1.
Oseen mobility directions
. Details.
In a dilute solution, separate coils diffuse through solvent with long-ranged hydrodynamic coupling among their segments. The Zimm model gives whole-chain diffusion
where is a hydrodynamic radius and the solvent viscosity. The coil continually changes shape while its centre of mass performs long-time Brownian motion; treating it as a rigid sphere is only a scaling analogy.
At high concentration, solvent-mediated hydrodynamic interactions are screened on sufficiently large scales and times. If chains are not entangled, additive local or blob friction gives Rouse model behaviour, at fixed concentration. For long entangled chains, neighboring chains prevent rapid transverse escape. In the reptation tube picture the chain moves mainly along its contour until it renews the surrounding tube. Contour diffusion scales as , the tube-renewal time as , and the random-coil size squared as . Therefore
Both cases ultimately give diffusive centre-of-mass motion, but entangled concentrated chains move much more slowly and retain constraints for longer. The inequality establishes strong overlap, not by itself entanglement: the latter also requires sufficiently many segments relative to an entanglement length. No numerical entanglement scale or concentrated-solution diffusion coefficient is supplied, so only the conditional regimes and scaling exponents can be fixed here.
Let a sphere of radius descend along the axis through a stationary unbounded Newtonian fluid. A second sphere with centre is held against translation but has zero applied torque. In the prescribed isolated-sphere incident flow, Faxén's rotational law and the vorticity of translating-sphere Stokes flow give
The speed cancels on integrating along the complete passage. Since , the result is the displayed angle about . The target radius is absent from the rotational Faxén response to this harmonic incident vorticity. This is a one-way incident-flow calculation; nonoverlap alone does not justify ignoring mutual hydrodynamic interactions between finite spheres. A fully coupled encounter requires the method of reflections for Stokes flow or an exact two-body solution.