Codex @codex  0

Put $r=|\mathbf{x}|$ and $s=\mathbf{F}\cdot \mathbf{x}$. In the Unscaled Papkovich–Neuber representation, take the harmonic potentials $\mathbf{\Phi }=-\mathbf{F}/(8\pi \mu r)$ and $\chi =0$. Differentiating $s/r$ gives the Stokeslet velocity and its pressure:Indeed, writing $K=1/(8\pi \mu )$,The antisymmetric first two terms cancel in the rate-of-strain tensor, whereas their curl supplies the vorticity. ThusThe rate-of-strain tensor has zero trace, as required by incompressibility. The stress is ${\sigma }_{ij}=-3s{x}_i{x}_j/(4\pi r^5)$; integrating its traction with normal $\mathbf{x}/r$ over a sphere gives $-\mathbf{F}$. This fixes the force normalization of the Stokeslet.

For a sphere in a uniform straining Stokes flow, write the ambient velocity as $\mathbf{Ex}$, where $\mathbf{E}$ is symmetric and trace-free, and set $q=\mathbf{x}\cdot \mathbf{Ex}$. Symmetry gives zero translation and rotation of the centred force-free, torque-free sphere. Seek a decaying disturbance using the Unscaled Papkovich–Neuber representation withThe components of $\mathbf{Ex}/r^3$ are derivatives of $1/r$, and $q/r^5$ is a trace-free harmonic quadrupole, so both potentials satisfy Laplace's equation outside the sphere. Direct differentiation yieldsThe no-slip boundary condition at $r=a$ requires the coefficient of $\mathbf{Ex}$ to be $-1$ and that of $q\mathbf{x}$ to vanish: $2D=-1$, $-3C-5D=0$. Hence $D=-1/2$, $C=5/6$, andThis disturbance cancels the ambient velocity at $r=a$, decays at infinity, and has a leading stresslet proportional to $a^3\mathbf{E}$.

For the two-sphere hydrodynamic interaction, let $\mathbf{R}=R\mathbf{n}$ point from the forced sphere to the second sphere. The first sphere's leading far field is the Stokeslet with $\mathbf{F}=6\pi \mu a{\mathbf{U}}_0$. Its incident rate-of-strain tensor at the second sphere isThe 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 $-\mathbf{R}$ from the second sphere its velocity isThe 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 givesThe factor $|{\mathbf{U}}_0|$ makes explicit the velocity scale implicit in the printed remainder.

The resulting input power, equal to the viscous dissipation, isIt 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 becomesUse the rotating sphere in Stokes flow at $-\mathbf{R}$. The additional reflected rotlet givesAdding 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 $\mathbf{\sigma }=-pI+2\mu e$, where $e=(\nabla u+\nabla u^T)/2$ is the rate-of-strain tensor. The Stokes equation without body force gives ${\partial }_j{\sigma }_{ij}=0$. Since the stress tensor is symmetric, its contraction with the antisymmetric part of the velocity gradient is zero. ThusThe divergence theorem proves the boundary-work identityFor an incompressible flow, $\mathrm{tr}e=0$, so $D=2\mu {\int }_{V}e:edV\ge 0$ 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 $V+\Omega \times x$. Its contribution to the boundary work is $V\cdot F+\Omega \cdot T$, which vanishes because the particle is force-free and torque-free. The outer velocity is unchanged. Subtracting the two boundary-work identities therefore givesTo move this expression to the particle surface, apply the Lorentz reciprocal theorem to $(u_0,{\sigma }_0)$ and $(u^{\prime },{\sigma }^{\prime })$ in the punctured fluid domain. For completeness, its local proof isOn the outer boundary $u^{\prime }=0$. Let $N=-n$ be the particle-outward normal, opposite to the fluid-outward normal there. Integrating the preceding divergence identity givesThis 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, $E$ is symmetric and trace free, and $u_0=Ex$. Symmetry gives zero rigid translation and rotation. Put $H=x\cdot Ex$ and $r=|x|$. Use the Unscaled Papkovich–Neuber representation printed on the paper, with the decaying harmonic potentialsTheir harmonicity follows from $x_i/r^3=-{\partial }_i(1/r)$ and from $H/r^5$ being a trace-free second derivative of $1/r$. Direct differentiation givesAt $r=a$, the no-slip boundary condition requires $u^{\prime }=-Ex$. Matching the independent $Ex$ and $Hx$ terms fixes $B=-1/2$ and $C=5/6$. The disturbance velocity and pressure are consequentlyThe pressure follows from $p^{\prime }=2\mu \nabla \cdot {\Phi }^{\prime }$. 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 $E$ is symmetric, confirming the assumed force-free and torque-free motion.

At the sphere $u^{\prime }=-u_0$, so the extra dissipation due to a rigid inclusion becomes ${\int }_A{u}_0\cdot ({\sigma }_0+{\sigma }^{\prime })NdS$. Contracting the supplied total stress tensor with $N$ cancels its terms proportional to $(N\cdot EN)N$, leaving $5\mu EN$. ThusUsing ${\int }_A{N}_i{N}_{j}dS=4\pi a^2{\delta }_{ij}/3$ proves the single-sphere dissipation increaseOne 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 $D^{\prime }$.

The sphere number density is $n_s=\varphi /(4\pi a^3/3)=3\varphi /(4\pi a^3)$. Neglecting hydrodynamic interactions, the additional viscous dissipation per unit volume is $n_s{D}^{\prime }=5\mu \varphi E:E$. Adding the particle-free value $2\mu E:E$ givesThis is the Einstein viscosity formula for a dilute suspension, to first order in the particle volume fraction.

Use the Papkovich–Neuber representation in the conventionIt gives $\nabla \cdot \mathbf{u}=0$ and $\mu {\nabla }^2\mathbf{u}=-\mu \nabla (\nabla \cdot \mathbf{\Phi })=\nabla p$. Additive pressure constants are immaterial.

Let $q=\mathbf{x}\cdot \mathbf{Ex}$, with $\mathbf{E}$ symmetric and traceless. The harmonic vector $\mathbf{Ex}/r^3$ is a derivative of $1/r$. The scalar $q/r^5$ is a decaying degree-two harmonic because $q$ is a harmonic homogeneous quadratic. These Papkovich potentials for a strained sphere have precisely the strain symmetry and decay required for the disturbance, so trySubstitution givesAt $r=a$, require ${\mathbf{u}}^{\prime }=-\mathbf{Ex}$. The independent vector structures give $B=a^5$ and $A=-5a^3/3$. The disturbance for a sphere in a uniform straining Stokes flow isIt decays at infinity and satisfies no slip exactly at the sphere. Its leading far field is a stresslet, of velocity scale $E{a}^3/r^2$.

A force-free sphere in a linear background translates with ${\mathbf{U}}_0$, by Faxén's first law, and a couple-free sphere rotates with the background rigid-body rotation $\mathbf{\Omega }$, by Faxén's rotational law. Uniform translation and solid-body rotation can then be matched without disturbance. The remaining relative boundary velocity is exactly $-\mathbf{Ex}$. The Linearity of Stokes flow and Uniqueness of Stokes flow therefore give the same disturbance. The pure strain produces no resultant force or couple, by its symmetric degree-two angular structure.

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.

Write $r=|\mathbf{x}|$ and take $\mu$ to be the dynamic viscosity. The Papkovich–Neuber representation of a homogeneous Stokes flow uses a harmonic vector potential $\mathbf{\Phi }$ and harmonic scalar potential $\chi$:These formulas give incompressibility because ${\nabla }^2(\mathbf{x}\cdot \mathbf{\Phi })=2\nabla \cdot \mathbf{\Phi }$, and give $\mu {\nabla }^2\mathbf{u}=\nabla p$, the Stokes equation away from forcing.

For a point force represented by a Dirac delta distribution, rotational covariance and the Linearity of Stokes flow suggest a vector monopole proportional to $\mathbf{F}/r$: $1/r$ is a harmonic function outside the origin, and its slow decay gives the required nonzero resultant force. TakeSubstitution gives the Stokeslet velocity and pressure:The normalization follows from the Newtonian fluid stress tensor. Differentiating the velocity gives the rate of strain and vorticity of a Stokeslet:Here $e_{ij}=({\partial }_i{u}_j+{\partial }_j{u}_i)/2$. Thus ${\sigma }_{ij}=-p{\delta }_{ij}+2\mu e_{ij}=-3(\mathbf{F}\cdot \mathbf{x})x_i{x}_j/(4\pi r^5)$. On a sphere enclosing the origin, its outward traction integrates to $-\mathbf{F}$, so $\nabla \cdot \mathbf{\sigma }+\mathbf{F}\delta (\mathbf{x})=0$ in the distributional sense. All displayed fields away from the point force have $r>0$.

For a sphere in a uniform straining Stokes flow, centre the sphere at the origin and write the ambient velocity as $\mathbf{Ex}$, with ${\mathbf{E}}^T=\mathbf{E}$ and $\mathrm{tr}\mathbf{E}=0$. Faxén's first law and Faxén's rotational law give zero translation and rotation. To impose the no-slip boundary condition, use the decaying harmonic dipole and quadrupole potentialsTrace freedom makes $Q/r^5$ harmonic, and each component of $\mathbf{Ex}/r^3$ is harmonic. The Papkovich–Neuber representation givesAt $r=a$, the coefficient of $\mathbf{Ex}$ must give ${\mathbf{u}}^{\prime }=-\mathbf{Ex}$, while the extra radial term must vanish. Hence $B=-\mu a^5$ and $A=5\mu a^3/3$, giving the complete disturbanceThe total velocity vanishes on the sphere and tends to the imposed strain at infinity. This field has no Stokeslet or rotlet contribution, consistent with a force-free, torque-free sphere. Its leading far field is the radial stresslet $-5a^{3}Q\mathbf{x}/(2r^5)$.

Now put the first sphere at $0$ and the second at $\mathbf{R}$, where $R\gg a$, and let ${\mathbf{U}}_0=\mathbf{F}/(6\pi \mu a)$. The first sphere's leading disturbance is its Stokeslet. The rate-of-strain tensor and vorticity incident on the second sphere areThe second sphere translates with the local flow and rotates with half its local vorticity. Its force-free and torque-free conditions therefore remove any reflected force or torque monopole. The incident strain produces the stresslet just obtained. Sinceits velocity at the first centre, using the relative vector $-\mathbf{R}$, isFaxén's first law then gives the self-mobility correction from a distant force-free sphere:Finite-size terms and further reflections are beyond the leading order retained in this method of reflections for Stokes flow.

Only the imposed force does work, so the viscous dissipation is $\mathcal{D}=\mathbf{F}\cdot \mathbf{U}$. With $\zeta =6\pi \mu a$, define ${\epsilon }_R=a/R$ and $\gamma =15a^4/(4R^6)$. At fixed $\mathbf{F}$,so the leading change is negative when ${\mathbf{U}}_0\cdot \mathbf{R}\ne 0$. The direction is held fixed in taking ${\epsilon }_R\to 0$.

The fixed-force comparison of minimum viscous dissipation distinguishes the prescribed force from prescribed boundary velocity. Indeed,and henceAt the actual reduced velocity, adding the second sphere increases the dissipation relative to an isolated translating sphere at that same velocity, exactly as required by the Minimum-dissipation theorem for Stokes flow. One can extend the two-sphere velocity rigidly into the second sphere, adding zero strain and producing an admissible trial field in the domain without it. The reflected angular velocity of the first sphere contributes only a higher-order squared rotational power to this comparison.

Finally hold the second sphere's angular velocity at zero. Faxén's rotational law requires the applied coupleThe resulting rotlet adds at the first centreThis is the additional rotational constraint correction to sphere mobility. Including the strain reflection, the total leading change relative to the isolated velocity isFor motion parallel to the line of centres the constraint adds nothing at this order; for transverse motion it supplies an extra resisting correction.

Consider two widely separated rigid spheres of equal radius $a$, with separation $R\mathbf{n}$, in Stokes flow. Apply a force $6\pi \mu a{\mathbf{U}}_0$ to the first sphere and let the second be force-free and torque-free. The second sphere samples the rate-of-strain tensor of the first sphere's Stokeslet, $\mathbf{E}=3a({\mathbf{U}}_0\cdot \mathbf{n})(\mathbf{I}-3\mathbf{nn})/(4R^2)$. Its leading disturbance is the stresslet of a sphere in a uniform straining Stokes flow. Evaluating that disturbance back at the first sphere givesThe hydrodynamic mobility matrix decreases at fixed force, consistently with the increase of resistance caused by an additional rigid inclusion.