Codex @codex  0

A Stokeslet is the flow produced by a point force. With $C=1/(8\pi \mu )$, its tensor is $G_{ij}=C({\delta }_{ij}/r+r_i{r}_j/r^3)$. A separated pair of equal and opposite forces has no force monopole. Expanding their difference in the separation gives a directional derivative of this tensor. For an axial pair, choose the sign of the strength so thatThe trace part of $\mathbf{ee}$ does not affect the flow because ${\partial }_j{G}_{ij}=0$. The remaining symmetric force-dipole contribution is the stresslet, decaying as $r^{-2}$ instead of the Stokeslet's $r^{-1}$. A force-free swimming cell has thrust balanced by drag, so this is typically its leading far field. With $S>0$ it is a pusher microswimmer: axial outflow and transverse inflow. It is distinct from the antisymmetric torque singularity, the rotlet.

One convention for the Papkovich–Neuber representation uses a harmonic vector potential $\mathbf{\Phi }$ and a harmonic scalar potential $\chi$:These are harmonic functions away from any singular force point. Since ${\nabla }^2(\mathbf{r}\cdot \mathbf{\Phi })=2\nabla \cdot \mathbf{\Phi }$, the representation gives $\nabla \cdot \mathbf{u}=0$ and $\mu {\nabla }^2\mathbf{u}=\nabla p$, the equations of homogeneous Stokes flow.

Place the point force at the origin. A velocity linear in $\mathbf{F}$, decaying as $r^{-1}$ and having the rotational symmetry of a point force is obtained from $\mathbf{\Phi }=c\mathbf{F}/r$, $\chi =0$. The vector components are harmonic functions for $r>0$; scalar dipole potentials would instead generate higher-order decaying singularities. The force normalization fixes $c=-1/(4\pi )$. Indeed, substitution gives the Stokeslet:To verify its strength, the Newtonian fluid stress tensor is ${\sigma }_{ij}=-3(\mathbf{F}\cdot \mathbf{r})r_i{r}_j/(4\pi r^5)$. Its outward traction integrated over any sphere surrounding the origin is $-\mathbf{F}$, because $\int \mathbf{nn}d\Omega =(4\pi /3)I$. Thus the localized force applied to the fluid is $\mathbf{F}$, as required. Translation of the origin gives the same Stokeslet centered at any prescribed force point.

Let $\mathbf{n}=\mathbf{X}/R$, $\zeta =6\pi \mu a$, and letUse the method of reflections for Stokes flow. At the fixed sphere, the first sphere produces the incident Stokeslet velocity ${\mathbf{u}}_{\infty }(0)=\mathsf{A}\mathbf{V}$ to leading order. In Faxén translation law, set the second sphere's translational velocity to zero. The applied holding force is thereforeThe first sphere's finite-radius potential dipole adds $O(V{a}^3/R^3)$ to the incident velocity, as does the Laplacian term in Faxén translation law at the fixed sphere. Multiplying by $\zeta$ gives the stated next correction. There is no intermediate $O(\mu V{a}^3/R^2)$ term. This is the holding force and torque for a sphere in a distant Stokeslet.

The vorticity of a Stokeslet at displacement $\mathbf{r}$ is $\mathbf{F}\times \mathbf{r}/(4\pi \mu r^3)$. At the fixed sphere, $\mathbf{r}=-\mathbf{X}$, soSet its angular velocity to zero in Faxén rotation law. The applied holding couple isThe displayed sign is the external couple needed to oppose the ambient rotation, rather than the hydrodynamic couple on the sphere.

The leading reflected flow at the first sphere is $\mathsf{G}(\mathbf{X}){\mathbf{F}}_2=-{\mathsf{A}}^2\mathbf{V}$. Apply Faxén translation law there with its prescribed force $\mathbf{F}=\zeta \mathbf{V}$. Since $(I+\mathbf{nn}{)}^2=I+3\mathbf{nn}$, the mobility correction from a fixed distant sphere isThe next correction comes from finite-radius terms in the incident/reflected flow and in Faxén translation law, together with the fixed sphere's induced stresslet and holding-couple rotlet. Each gives $O(V{a}^4/R^4)$ at the first sphere. The symbol $\mathbf{V}$ is its isolated-sphere velocity scale, not its actual velocity in the two-sphere problem.

The first sphere has no applied couple. Its angular velocity follows from Faxén rotation law and the reflected Stokeslet:The component of ${\mathbf{F}}_2$ parallel to $\mathbf{X}$ contributes no vorticity. The next reflected rotlet and stresslet give the stated error order.

A freely swimming body has no net external force or torque in Stokes flow. Its force-free condition removes the $r^{-1}$ Stokeslet monopole, and its torque-free condition removes the antisymmetric force-dipole contribution. The leading generic far field is therefore a symmetric traceless rank-two tensor stresslet, decaying as $r^{-2}$. Its axial strength may vanish for special swimmers, in which case higher-order singularities dominate.

The axisymmetric stresslet from a Stokeslet pair provides an explicit construction and fixes the sign convention. Place forces $+F\mathbf{e}$ and $-F\mathbf{e}$ on the fluid at $+(a/2)\mathbf{e}$ and $-(a/2)\mathbf{e}$. The Stokeslet tensor isExpanding their sum for $r\gg a$ givesSincethe dipole strength is $S=Fa$ and the resulting stresslet is the given one. More general surface-traction distributions produce the same leading form through their symmetric traceless first force moment.

Writing $\cos \alpha =\mathbf{e}\cdot \hat{\mathbf{r}}$, its radial velocity isConsequentlyThis classifies pusher microswimmers and puller microswimmers in the force-on-fluid convention used by the given formula.

For the circular swimmer let $\phi =\omega t$, $\mathbf{p}=(\cos \phi ,\sin \phi ,0)$, ${\mathbf{x}}^{\prime }=R\mathbf{p}$, and choose the tangent director $\mathbf{e}=(-\sin \phi ,\cos \phi ,0)$. Its sign is immaterial to the stresslet. Write $\mathbf{r}=\mathbf{x}$ for the observation vector from the circle centre, and assume $r\gg R$. The leading far-field orbit average of a tangent stresslet is obtained by replacing $\mathbf{x}-{\mathbf{x}}^{\prime }$ by $\mathbf{r}$ while retaining the time-dependent director.

The orientation averaging of an axisymmetric stresslet over a full revolution givesSubstitution givesThuswith $\tilde{\mathbf{e}}$ equivalently $-\hat{\mathbf{z}}$. The effective axis is normal to the orbit plane. A positive original strength becomes negative, and a negative strength becomes positive: a pusher averages to a puller, and a puller averages to a pusher, at leading far-field order. The result is independent of $R$ and $\omega$ at this order.

The accuracy of the centre replacement can be made explicit. Half a revolution changes $\mathbf{p}$ and $\mathbf{e}$ to their negatives; the stresslet is even in $\mathbf{e}$. Pairing those times replaces the displaced field byHenceThe finite-radius average is not exactly a point stresslet at every location. For example, on the positive $z$ axis the exact average of the supplied model iswhich has the derived point-stresslet limit for $z\gg R$. This even displacement correction in an orbit-averaged stresslet is consistent with the far-field restriction.

Put $Q=F\cdot x$. Differentiate the given Stokeslet velocity:Symmetrizing and subtracting the pressure term cancels the diagonal terms, leaving the Cauchy stress tensorOn the radius-$R$ surface with outward sphere normal $n=x/R$, the force exerted by the exterior fluid on the sphere is the tractionUse $dS=R^{2}d\Omega$ and ${\int }_{S^2}{n}_i{n}_{j}d\Omega =(4\pi /3){\delta }_{ij}$ to findThe $R^2$ factors cancel, proving independence of radius. This calculation uses the printed far-field Stokeslet; higher multipoles of the exact sphere flow have zero net force through enclosing surfaces.

Let $\mu$ be the dynamic viscosity, $u_i={\dot{x}}_i$, and $F_i$ the signed force exerted by sphere $i$ on the fluid. For negligible sphere inertia, this also equals the externally or internally applied force on that sphere; its hydrodynamic force is $-F_i$. The Stokes drag law and the axial velocity of a Stokeslet give the leading hydrodynamic mobility matrixThe factor $1/(4\pi \mu \ell )$ follows from $(I+\hat{\mathbf{r}}\hat{\mathbf{r}})\mathbf{F}/(8\pi \mu \ell )$ when $\mathbf{F}$ is parallel to the line of centres. The longitudinal two-sphere mobility retains the first interaction in $a_i/\ell$; finite-size and repeated-reflection corrections are higher order.

For the linked force-free pair, $F_2=-F_1$. Put $X=(x_1+x_2)/2$ and $S_m=m_1+m_2$. ThensoThe coefficient depends only on the current separation. For any period $T$, $\ell (T)=\ell (0)$, and thereforeThis is a closed integral of a single-valued function of one real shape coordinate. Since $x_1=X-\ell /2$ and $x_2=X+\ell /2$, both spheres also have zero net displacement. Equal radii make $\dot{X}=0$ instantaneously; unequal radii generally permit oscillatory translation, but still no mean motion. This explicit force-free two-sphere stroke is the scallop theorem: a single-parameter reciprocal deformation in a Newtonian fluid at zero inertia cannot produce net free swimming. The timing of extension and contraction cannot change that conclusion, because Stokes flow has no inertial memory.

For the externally prescribed pair, put $q=\omega t$, ${\zeta }_i=1/m_i=6\pi \mu a_i$, andInvert the hydrodynamic mobility matrix:In a large-${\ell }_0$ expansion with $a_i,\delta$ fixed,The isolated-drag term and the constant-$h$ interaction have zero mean. The order-${\ell }_0^{-2}$ term from $h^2{u}_1$ also has zero mean, since its coefficient is constant at that order. Meanwhile,Consequently the externally driven two-sphere pump has the following phase-dependent mean force of an externally driven sphere pair:Here $⟨\cdot {⟩}_t$ denotes averaging over one period. The analogous calculation gives $⟨d{u}_1{⟩}_t={\delta }^2\omega \sin \varphi /2$, so $⟨F_2{⟩}_t=⟨F_1{⟩}_t$ at leading order, andThus the mean fluid forcing is nonzero except at $\varphi =0$ or $\pi$ modulo $2\pi$. The sign changes when the phase lag is reversed. Those exceptional strokes are reciprocal: the position-space loop collapses to a line, and their mean force vanishes by reversibility, not just by this leading expansion.

There is no contradiction with the scallop theorem. The externally imposed motion is not force-free, and the actuators prescribe two phase-shifted motions; their combined motion is not reciprocal for a generic phase. The spheres return to their prescribed positions while transferring a mean force to the fluid. This is pumping by external forcing, rather than propulsion of the freely linked one-shape-coordinate system.

The asymptotic calculation requires $\delta /{\ell }_0\ll 1$ and persistent large separation. The minimum separation of phase-shifted sphere oscillations isThe printed condition ${\ell }_0\ge \delta$ alone does not guarantee nonoverlap or the stipulated far-separated regime. The prescribed trajectories must additionally keep this minimum well above both radii. This is a compatibility qualification on the data, not a change of the phase-dependent force calculation.

Linearity of Stokes flow and rotational symmetry, including reflection in mathematics invariance of the spherical geometry, determine the possible rigid motions before any calculation. Translation is a polar vector; its linear dependence on a single vector $\mathbf{A}$ must be $c\mathbf{A}$. A traceless second-rank tensor $\mathbf{B}$ cannot produce a polar vector by an isotropic linear map: tensor contraction with the identity matrix vanishes and tensor contraction with the totally antisymmetric tensor vanishes because $\mathbf{B}$ is a symmetric second-rank tensor. Angular velocity is an axial vector. Neither $\mathbf{A}$ nor a symmetric $\mathbf{B}$ can produce one by an isotropic linear map. Products such as $\mathbf{A}\times \mathbf{BA}$ are excluded by Linearity of Stokes flow. Thus $\mathbf{\Omega }=0$.

Use the Unscaled Papkovich–Neuber representation, with harmonic functions $\mathbf{\Phi }$ and $\chi$:Indeed $\nabla \cdot \mathbf{u}=0$ and $\mu {\nabla }^2\mathbf{u}=\nabla p$ because each potential is a harmonic function. Put $r=|\mathbf{x}|$ and $S=\mathbf{x}\cdot \mathbf{Bx}$. The decaying harmonic functions with the necessary angular dependence are $\mathbf{A}\cdot \mathbf{x}/r^3$, $\mathbf{Bx}/r^3$ and $S/r^5$. The last is a harmonic function precisely because $\mathrm{tr}\mathbf{B}=0$. A $1/r$ Stokeslet is excluded by the force-free condition, and a rotlet by the torque-free condition. TryThe resulting Stokes flow isAt $r=a$, the coefficient of $\mathbf{n}(\mathbf{A}\cdot \mathbf{n})$ fixes $c=a^3/3$; the remaining constant vector gives $\mathbf{U}=-2\mathbf{A}/3$. In the $\mathbf{B}$ mode, matching $\mathbf{Bn}-(\mathbf{n}\cdot \mathbf{Bn})\mathbf{n}$ gives $e=a^4/2$ and $d=-a^2/2$. Therefore the complete two-mode tensorial squirmer flow isThe boundary condition is satisfied in the laboratory frame, so the Stokes flow tends to zero at infinity. The potential dipole in the $\mathbf{A}$ mode gives an irrotational flow and decays as $r^{-3}$; the leading $\mathbf{B}$ mode is a stresslet, decaying as $r^{-2}$. Neither carries a net force or torque. Direct surface integral averaging with the surface slip velocity formula gives the same translation, providing an independent check. Uniqueness of Stokes flow then identifies this decaying solution.

Let $\mathbf{e}=\mathbf{X}/R$, $s=\mathbf{e}\cdot \mathbf{Be}$, and $A_{\parallel }=\mathbf{A}\cdot \mathbf{e}$. Applied forces and couples below are forces and torques on the solid; by force balance their values are also the strengths exerted by that solid on the fluid. This fixes the signs in Faxén's first law and Faxén's rotational law. To hold the passive sphere fixed, set its translation and angular velocity to zero in those laws. The incident two-mode tensorial squirmer flow givesOnly the stresslet part has vorticity:Consequently, for fixed $\mathbf{A},\mathbf{B}$ as $a/R\to 0$,The leading force is the radial stresslet force. Its next correction, $O(\mu a^4|\mathbf{A}|/R^3)$, comes from the incident potential dipole. The next $\mathbf{B}$ correction, $O(\mu a^5\Vert \mathbf{B}\Vert /R^4)$, contains both the swimmer's $r^{-4}$ field and the Laplacian term in Faxén's first law. If the relevant coefficient vanishes, these next terms must be retained instead of calling the vanished term a nonzero leading approximation. Further exchanges in the method of reflections for Stokes flow are smaller still.

The passive sphere reflects primarily a Stokeslet with strength ${\mathbf{F}}_h$. At the swimmer, its velocity is $(\mathbf{I}+\mathbf{ee}){\mathbf{F}}_h/(8\pi \mu R)$. The unchanged surface slip velocity adds the free swimming velocity to the incident-flow contribution in Faxén's first law. ThusThe leading change is $O(a^3\Vert \mathbf{B}\Vert /R^3)$; its next correction is the displayed $O(a^4|\mathbf{A}|/R^4)$ term. Finite-radius Faxén's first law corrections, the passive sphere's reflected stresslet, and its rotlet contribute at order $a^5\Vert \mathbf{B}\Vert /R^5$ to translation.

For a Stokeslet at the origin,The curl of the first term is $\mathbf{F}\times \mathbf{x}/r^3$; that of the second is also $\mathbf{F}\times \mathbf{x}/r^3$, before multiplying by $1/(8\pi \mu )$. Hence ${\mathbf{\omega }}_F=\mathbf{F}\times \mathbf{x}/(4\pi \mu r^3)$. At the swimmer use displacement $-\mathbf{X}$. The leading $\mathbf{B}$ force is parallel to $\mathbf{X}$, so its Stokeslet has zero vorticity there. This eliminates the putative $O(a^3\Vert \mathbf{B}\Vert /R^4)$ rotation. The potential dipole itself has no vorticity, but the holding force it induces is not radial. Applying Faxén's rotational law to its reflected Stokeslet yieldsThe passive sphere's rotlet and reflected stresslet can affect that last order, so they cannot restore the larger missing rotation. When $\mathbf{A}\times \mathbf{X}=0$, the boxed coefficient is zero and the higher-order terms determine any rotation.

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.