 Solution

Use the Unscaled Papkovich–Neuber representationwith harmonic vector and scalar potentials. For a torque-driven sphere chooseIts components are harmonic for $r>a$, and both $\mathbf{x}\cdot \mathbf{\Phi }$ and $\nabla \cdot \mathbf{\Phi }$ vanish. HenceThe last equality follows from the no-slip boundary condition at $r=a$. To check that $\mathbf{G}$ is indeed the applied couple, the traction on the sphere, with $\mathbf{n}=\mathbf{x}/a$ directed into the fluid, is ${\mathbf{\sigma }}_G\mathbf{n}=-3\mu {\mathbf{\Omega }}_G\times \mathbf{n}$. Its hydrodynamic torque isHere $\int \mathbf{nn}dS=(4\pi a^2/3)I$. The applied couple balances this resisting torque. The resulting field is the rotlet of a rotating sphere in Stokes flow.

For the Lorentz reciprocal theorem for Stokes flow, let ${\mathbf{n}}_D$ point out of the fluid domain $D$ and take the same viscosity in both flows. With $\nabla \cdot {\mathbf{\sigma }}_i+{\mathbf{f}}_i=0$, the theorem isIt follows by taking the divergence of the boundary integrand: incompressibility removes pressure-work terms, the two Newtonian strain contractions are equal, and the remaining stress divergences give the displayed body-force terms. This convention matters: on the inner boundary of an exterior domain, ${\mathbf{n}}_D=-\mathbf{n}$.