 Solution

For an incompressible flow, write the Newtonian fluid stress tensor as $\sigma =-pI+2\mu e$, where $e=(\nabla u+\nabla u^T)/2$ is the rate-of-strain tensor. The Stokes equation gives ${\partial }_j{\sigma }_{ij}=0$. Symmetry of the Cauchy stress tensor therefore givesIntegrating and applying the divergence theorem expresses the viscous dissipation as boundary power:At a moving rigid body, the surface velocity is $U+\Omega \times (x-x_c)$. Its boundary-power contribution is $U\cdot F+\Omega \cdot G$, 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 givesApply the Lorentz reciprocal theorem to $(u_0,{\sigma }_0)$ and $(u^{\prime },{\sigma }^{\prime })$ in the fluid outside the inclusion. On the outer boundary $u^{\prime }=0$, soWriting $n_p=-n$ 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 $E$ 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 $E$ eliminates its torque. To derive the disturbance, put $r=|x|$, $Q=x\cdot Ex$ and seekThe incompressible flow condition and Stokes equation reduce toA decaying family satisfying these equations is $f=A{r}^{-5}$, $g=B{r}^{-5}-5A{r}^{-7}/2$, $j=2B{r}^{-5}$. The no-slip boundary condition requires $f(a)=-1$ and $g(a)=0$, so $A=-a^5$ and $B=-5a^3/2$. ConsequentlyThe total velocity is zero at $r=a$, and the disturbance decays as $r^{-2}$. 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 $u^{\prime }=-u_0=-aE{n}_p$. Contracting the supplied Newtonian fluid stress tensor with $n_p$ makes its two terms proportional to $(n_p\cdot E{n}_p)n_p$ cancel, leaving $\sigma n_p=5\mu E{n}_p$. Hence the extra dissipation due to a rigid inclusion isHere ${\int }_{S^2}{n}_j{n}_{k}d\Omega =4\pi {\delta }_{jk}/3$. 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 $N_v$ denotes the number of spheres per unit volume, $N_v(4\pi a^3/3)=\varphi$. ThusAdd this to the particle-free viscous dissipation density $2\mu E:E$. Comparing with $2{\mu }_{\mathrm{eff}}E:E$ gives the Einstein viscosity formula for a dilute suspension:to first order in $\varphi$. The hydrodynamic interactions neglected here contribute beyond that dilute order.