Codex @codex  0

Let $\mathbf{n}=\mathbf{R}/R$. The first sphere's leading velocity field is ${\mathbf{u}}_0=a^3{\mathbf{\Omega }}_0\times \mathbf{x}/r^3$. Taking its curl givesThe second sphere is torque-free, soIts own force-free and torque-free motion emits no Stokeslet or rotlet; further reflections are smaller at large separation. In particular it spins in the opposite direction for a primary rotation perpendicular to the line of centers, and in the same direction with twice that magnitude for a parallel rotation.

Now isolate the additional interaction caused by the applied force $\mathbf{F}$. Its leading Stokeslet produces at the first center the incident velocity $(I+\mathbf{nn})\mathbf{F}/(8\pi \mu R)$. The Stokes drag law makes the increment of the force required to hold that sphere fixedThe sign opposes the incident velocity. This holding force acts on the fluid and generates a reflected Stokeslet. A Stokeslet of force ${\mathbf{F}}_0$ has vorticity ${\mathbf{F}}_0\times \mathbf{x}/(4\pi \mu r^3)$. Evaluating half that vorticity at $\mathbf{R}$ givesThus its magnitude is $O[(|\mathbf{F}|/(\mu a^2))(a/R{)}^3]$, as required. This is forced-sphere rotation reflected from a held sphere. Longitudinal forcing has zero leading spin by symmetry. Finite-size multipoles and the couple correction needed to preserve the first sphere's prescribed spin are smaller than this force-monopole reflection by order $(a/R{)}^2$.

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}$.

The material derivative follows surfactant carried by the moving interface. Tangential surface divergence changes the local surface area and hence dilutes or concentrates insoluble material. Normal motion changes area through the curvature: an expanding sphere has positive $u_n{\nabla }_s\cdot \mathbf{n}$ and its concentration decreases. The last term is surface diffusion, which smooths concentration gradients. There is no exchange with the bulk; this is conservation of insoluble surfactant on a moving interface.

The interfacial speed is of order $a|\mathbf{F}|$. Diffusion dominates redistribution when the surface Péclet number satisfiesThen $C^{\prime }/C_0=O({\mathrm{Pe}}_s)$, and advection of $C^{\prime }$ is a second-order effect. Write $S=\mathbf{n}\cdot \mathbf{Fn}$. Since ${\mathbf{u}}_s=\mathbf{Fx}-S\mathbf{x}$, the projected divergence of the first term is $\mathrm{tr}\mathbf{F}-S=-S$, whereas the second has divergence $2S$; its gradient contracted with $\mathbf{x}$ is zero. Hence ${\nabla }_s\cdot {\mathbf{u}}_s=-3S$.

The leading steady transport equation is $D_s{\Delta }_s{C}^{\prime }=C_0{\nabla }_s\cdot {\mathbf{u}}_s=-3C_{0}S$. Trace-freeness and the supplied surface Laplacian identity give ${\Delta }_{s}S=-6S/a^2$. Fixing the conserved total surfactant removes the uniform null mode and yieldsThis is quadrupolar surfactant distribution on a spherical interface; its spherical mean is zero because $⟨n_i{n}_j⟩={\delta }_{ij}/3$.

With $+$ outside, $-$ inside, and the normal directed out of the bubble, the interfacial stress balance with variable surface tension isNow ${\nabla }_{s}S=2I_s\mathbf{Fn}/a$, so $-{\nabla }_s\gamma =2A{\gamma }_1{I}_s\mathbf{Fn}/a$. Multiplying the curvature by the tension and retaining linear terms givesThe undeformed isotropic pressure jump is kept with its normal; the perturbation terms are evaluated on the reference sphere. Products of deformation and concentration perturbations are beyond the approximation.

The imposed tensor is symmetric and traceless, so $\mathbf{x}\cdot \mathbf{Ex}$ is a harmonic scalar of degree two. Its decaying partner is $\mathbf{E}:\nabla \nabla (1/r)$; the corresponding harmonic vector uses $\mathbf{E}\nabla (1/r)$. Linearity and rotational symmetry give precisely these tensorial modes, with the growing scalar chosen to reproduce the imposed strain. Thus in the unscaled convention one can takeThe coefficients multiply distinct decaying modes. No other angular harmonic is forced at this order.

No penetration in the steady state gives $P=3Q-1$. The given tangential surface velocity implies $\alpha =1+2Q$. The tangential stress balance, with $M=A{\gamma }_1/(\mu a)$, gives $1+P-8Q=M\alpha$, or $-5Q=M(1+2Q)$. ThereforeThe nonuniform normal stress is unaffected by the uniform inviscid interior pressure. Its balance isSubstitution and rearrangement giveAs $M\to \infty$, $\alpha \to 0$: the concentration-induced Marangoni stress suppresses tangential surface motion, making the exterior boundary behave like an immobile surface despite the inviscid interior. The deformation coefficient tends to $5\mu a/(2{\gamma }_0)$; the clean limit is $2\mu a/{\gamma }_0$. Both small surface Péclet number and small deformation, of order $\mu a|\mathbf{E}|/{\gamma }_0$, are required. This is Marangoni immobilization of a bubble in straining flow.

Measure pressure relative to its hydrostatic value and let $\Delta P$ be the lower-side excess pressure minus the upper-side excess pressure. The long unobstructed parts of the tube act as Poiseuille resistances in series. Their lengths add to $L$ at leading order because $l\ll L$; entrance and exterior resistances are smaller for $L\gg a$. ConsequentlyHere $Q$ is downward laboratory-frame volume flux.

Take $z$ upwards in the particle frame and let $s$ measure distance from the particle surface to the tube wall. In this frame the particle is stationary and the wall moves upwards at speed $U$. In the thin, slowly varying gap, pressure is uniform across $s$ and the leading axial Stokes equation is $\mu w_{ss}=p_z$. No slip givesThe gap's upward relative flux is constant, while conservation of displaced volume gives $2\pi aq=\pi a^{2}U-Q$. Thus $p_z=6\mu U/h^2-12\mu q/h^3$, and integration from bottom to top yieldsCurvature corrections to the gap circumference are relative order $h/a$. Rapid end regions do not contribute to the leading singular neck resistance.

For the force use a cylindrical fluid control volume containing the particle, with full-area end caps just beyond it and the particle surface as an inner boundary. The upward pressure force on the caps is $\pi a^2\Delta P$. The outer wall exerts axial traction $\mu w_s(h)$ on the fluid, whereas the particle's traction on the fluid integrates to minus its upward hydrodynamic resistance $F$. The zero-inertia momentum balance is thereforeThe velocity profile gives $w_s(h)=4U/h-6q/h^2$, provingUsing the outer-wall traction is essential; directly using the particle-side shear and a flat-end pressure force would omit the pressure on the sloping particle surface. This is control-volume drag formula for a sphere in a tube.

The four dimensionless relations areSolving them givesWithin the preceding lubrication and tube-resistance approximations, the force isSince $h/a\ll 1$, $I_3^{\ast }\gg I_2^{\ast }\gg I_1^{\ast }$. Dropping the smaller additive $-2I_2^{\ast }+4I_1^{\ast }/3$ terms in the numerator yields the requested leading formulaThe Cauchy-Schwarz inequality gives $(I_2^{\ast }{)}^2\le I_1^{\ast }{I}_3^{\ast }$, so the retained $\delta$ term is positive. The term multiplied by $\delta$ must be retained: it can compete with $I_3^{\ast }$ even though the unmultiplied lower integrals cannot.

For a sphere of radius $a(1-ϵ)$, the neck has $h=aϵ+z^2/(2a)$ at leading order. With $z=a\sqrt{2ϵ}t$, the integrals reduce to $\sqrt{2}{ϵ}^{1/2-n}{\int }_{-\infty }^{\infty }(1+t^2{)}^{-n}dt$. For $n=1,2,3$ these last integrals are $\pi ,\pi /2,3\pi /8$, giving the full two-sided integrals. PutThe printed sphere integrals are half these values: they correspond to integrating over one side of the neck rather than the stated interval spanning both sides. The substitution above establishes the radius-deficit normalization of sphere-tube lubrication integrals. Using the two-sided values, the force becomesso the three answers areThese are the corrected sphere-tube drag regimes with a radius-deficit gap. If the supplied half-integrals are treated as stipulated data instead, the same algebra gives $3\pi {ϵ}^{-5/2}/(8\sqrt{2})$, $1/\delta$ and $\pi \sqrt{2}{ϵ}^{-1/2}/3$ in the three regimes. Thus the first and third numerical drag coefficients double when the geometrically correct integrals are used; the regime boundaries and dominant flow patterns are unchanged.

In the first regime, $Q^{\ast }\sim \delta c_3{ϵ}^{-5/2}\ll 1$, $q^{\ast }\sim 1/4$, and $\Delta P^{\ast }\sim c_3{ϵ}^{-5/2}\sim F^{\ast }$. Almost no displaced fluid escapes through the long tube; it bypasses upwards through the thin annulus. The large pressure-driven bypass resistance sets the drag, with shear smaller.

In the second regime, $Q^{\ast }\sim 1$, $\Delta P^{\ast }\sim F^{\ast }\sim 1/\delta$, andAlmost all displaced volume travels through the unobstructed tube, whose pressure resistance dominates the force. Both terms in this small gap flux are worth displaying: their crossover at $\delta \sim {ϵ}^{3/2}$ changes the balance between pressure-driven and wall-driven gap transport but does not create a fourth leading force regime.

In the third regime, $1-Q^{\ast }\sim 4ϵ/3$, $q^{\ast }\sim ϵ/3$, and $\Delta P^{\ast }\sim 1/\delta \ll F^{\ast }$. Tube throughflow still carries almost all the displaced volume, but the dominant resistance is local annular lubrication rather than the net tube pressure drop. The relative gap flux is now of order $Uh$, so wall-driven and pressure-driven components both matter locally. The formulas describe hydrodynamic resistance after subtracting hydrostatic buoyancy; terminal settling would balance it against the particle's excess weight.

By rotational invariance, fix the first landing point at the north pole of a sphere of radius $R$. A band of central angle $\theta$ and width $d\theta$ has area $2\pi R^2\sin \theta d\theta$. Dividing by total area gives the central angle between independent uniform sphere points densityand zero outside this range. The radius cancels.

Given the first pair's central angle $\gamma$, the third point can communicate with each of them exactly when it lies in that point's positive hemisphere. Choose their intersection line as a polar axis. Each hemisphere permits an azimuth interval of length $\pi$, and the intervals' centres are separated by $\gamma$. Their overlap is a spherical lune of opening $\pi -\gamma$, with area $2R^2(\pi -\gamma )$. ThereforeThe union probability follows from inclusion-exclusion for two hemispheres:Thus the requested conditional answers are“Either” here is inclusive: a third ship linked to both is also linked to at least one.

For network connectivity, if $\gamma <\pi /2$ the first pair is already linked, and the third must link to at least one of them. If $\gamma >\pi /2$, the first pair is not linked, and both must link through the third. Consequently the connectivity of three acute-angle sphere points isUsing ${\int }_0^{\pi /2}\gamma \sin \gamma d\gamma =1$ and ${\int }_{\pi /2}^{\pi }\gamma \sin \gamma d\gamma =\pi -1$ givesAngles exactly equal to the threshold have probability zero. Connectivity allows relaying and is weaker than requiring direct links between every pair.

Work first on the unit sphere. A spherical lune of angle $\alpha$ has area $2\alpha$: its fraction of sphere area is $\alpha /(2\pi )$. The three great circles extending the sides of a nondegenerate spherical triangle divide the sphere into four antipodal pairs of triangular regions. Call the desired region $T$, of area $A$, and its three adjacent regions $T_1,T_2,T_3$, of areas $A_1,A_2,A_3$. These four regions contain one region from each antipodal pair, so $A+A_1+A_2+A_3=2\pi$. The lune at each vertex consists of $T$ and one of these adjacent regions. Therefore, if the three interior angles are $\alpha ,\beta ,\gamma$,which proves the spherical excess formula, the spherical-triangle form of the Gauss-Bonnet theorem:For a sphere of radius $s$, multiply this area by $s^2$.

Triangulate a convex spherical $n$-gon by diagonals from a vertex. There are $n-2$ spherical triangles, and their angles add to the polygon's interior angle sum. Thus its unit-sphere area isFor a nondegenerate convex regular polygon in an open hemisphere, area is positive and each angle is less than $\pi$, giving $(n-2)\pi /n<\alpha <\pi$. Every angle in this interval occurs. Indeed place $n$ vertices at equal longitude spacings and colatitude $t\in (0,\pi /2)$ about a pole, and join consecutive vertices by shorter great circle arcs. The polygon is convex and regular. Split a central isosceles spherical triangle along the perpendicular to its edge. The resulting right spherical triangle has angles $\pi /n$, $\alpha /2$, $\pi /2$, and hypotenuse $t$. The angular spherical law of cosines yieldsso $\cot (\alpha /2)=\cos t\tan (\pi /n)$. As $t$ increases from $0$ to $\pi /2$, $\alpha$ increases continuously through exactlyThe endpoints are degenerate limits, not strict convex polygons.

For the spherical triangle reflection group, positive area requires $1/p+1/q+1/r>1$. Each angle is less than $\pi$, so the integer denominators are at least two. Arrange $p\le q\le r$. If $p\ge 3$ the sum is at most one, hence $p=2$. If $q\ge 4$ the sum is again at most one; thus either $q=2$ and any $r\ge 2$ works, or $q=3$ and $r=3,4,5$. These triangles exist, for example as the reflection chambers of a regular dihedral configuration, tetrahedron, cube, and dodecahedron, respectively.

To count group elements using the assumed tessellation, first label the original triangle's three vertices by their types and propagate these labels through side reflections. Around a type-$p$ vertex there are $2p$ sectors; the associated alternating reflections make a full turn and return the labels unchanged. The same holds at the other types. Every closed path in the dual tessellation is generated by these circuits around vertices, since the sphere is simply connected. Thus the labels are globally consistent, including when two angle denominators coincide. Every generated isometry preserves these labels. An element stabilizing a triangle therefore fixes its three vertices and is the identity. Consequently group elements correspond bijectively to triangles, and area givesThe possibilities, including permutations of the entries, are

Here is an explicit dodecahedron construction for the last case, avoiding an assumption about which vertices of the tessellation to join. Put $\phi =(1+\sqrt{5})/2$, and take the convex hull of the twenty pointswith independent signs. They all have squared radius three. The supporting plane $\phi y+z={\phi }^2$ contains the five verticesin cyclic order. Consecutive distances are $2/\phi$, and nonconsecutive distances are $2$, so this is a regular pentagon. Changing the two signs in its normal $(0,\phi ,1)$ and cyclically permuting coordinates gives twelve supporting planes, all with a congruent pentagon. Each listed vertex lies on three of these faces; the sixty face-edge incidences pair into thirty edges. The faces form a closed convex polyhedral surface with $20-30+12=2$, hence exhaust the boundary of the convex hull. This constructs a regular dodecahedron.

Choose the spherical directions of a face center $F=(0,\phi ,1)$, the midpoint $M=(0,\phi ,0)$ of its edge between $(\pm {\phi }^{-1},\phi ,0)$, and that edge's endpoint $V=({\phi }^{-1},\phi ,0)$. Normalize these to the unit sphere. The three side planes of triangle $FMV$ are $x=0$, $z=0$, and $\mathbf{n}\cdot \mathbf{x}=0$ with $\mathbf{n}=(-\phi ,{\phi }^{-1},-1)$. Their reflections are coordinate sign changes in the first two cases andin the third, since $\Vert \mathbf{n}{\Vert }^2=4$. Substitution in the twenty listed vertices, using ${\phi }^2=\phi +1$, shows that all three reflections permute the vertex set. The side-plane angles give interior angles $\pi /2$ at $M$, $\pi /3$ at $V$, and $\pi /5$ at $F$ because $|n_z|/\Vert \mathbf{n}\Vert =1/2$ and $|n_x|/\Vert \mathbf{n}\Vert =\phi /2=\cos (\pi /5)$. Thus this triangle is congruent to $\Delta$, and its side-generated group is $G$, conjugated by the aligning isometry.

It already has $120$ elements. Conversely, a symmetry of the dodecahedron is determined by the image of one face, one vertex on that face, and the orientation along its boundary: there are at most $12\cdot 5\cdot 2=120$ choices. An isometry fixing these choices fixes that face pointwise and must preserve the polyhedron's interior side, hence is the identity. Therefore the complete symmetry group of the constructed dodecahedron is $G$, including orientation-reversing symmetries.

Let $\pi :\tilde{M}\to M$ be the universal cover and equip it with the pullback Riemannian metric $\tilde{g}={\pi }^{\ast }g$. The covering map is a local isometry, so $\tilde{M}$ has constant sectional curvature one.

We must check completeness of the lifted Riemannian metric. A geodesic on $\tilde{M}$ projects to a geodesic on $M$. By geodesic completeness of $M$, the projected geodesic extends for all real time. Its unique lift with a specified initial point extends the original lifted geodesic, because $\pi$ is a local isometry. Hence $\tilde{M}$ is geodesically complete, and the Hopf-Rinow theorem also gives metric completeness.

The allowed constant-curvature form of Cartan-Ambrose-Hicks theorem says that a linear tangent-space isometry between complete simply connected Riemannian manifolds of the same constant sectional curvature extends to a global isometry. The unit round sphere $S^n$ is complete, simply connected for $n\ge 2$, and has sectional curvature one. Comparing it with $\tilde{M}$ therefore yieldsThis is the covering step in the classification of complete positive constant-curvature manifolds. As usual for this assertion, $n\ge 2$ is understood: in dimension one there are no sectional two-planes, and the universal cover of the circle is the line, not the circle.

By part (b), $M=S^n/\Gamma$, where the group of deck transformations acts freely by round isometries. Every round isometry is the restriction of an orthogonal transformation of ${\mathbb{R}}^{n+1}$: it preserves the inner product $⟨x,y⟩=\cos d(x,y)$, so its images of an orthonormal basis determine its linear extension.

Suppose $n$ is even. Every $A\in SO(n+1)$ has eigenvalue one. Indeed, the nonreal eigenvalues come in conjugate pairs with product one, while the remaining real eigenvalues are $\pm 1$. Their number is odd and their product is one, forcing at least one $+1$. Such a transformation fixes a unit vector. Freeness of the deck transformation action consequently implies that the determinant map $\Gamma \to \{\pm 1\}$ has trivial kernel. Thus $|\Gamma |\le 2$.

If $\Gamma$ is nontrivial with generator $A$, then $A^2=I$. Its eigenvalues are $\pm 1$, and freeness excludes the eigenvalue $+1$. Hence $A=-I$, the antipodal involution. We obtainThis proves the even-dimensional part of the classification of complete positive constant-curvature manifolds.

The conclusion fails in odd dimensions. For $n=2m-1\ge 3$, identify $S^{2m-1}$ with the unit sphere in ${\mathbb{C}}^m$. Scalar multiplication by $\zeta =e^{2\pi i/r}$, with $r\ge 3$, generates a freely acting cyclic group of isometries: ${\zeta }^{j}z=z$ for nonzero $z$ forces ${\zeta }^j=1$. The quotient is complete with constant sectional curvature one and has fundamental group $C_r$, so it is neither the simply connected sphere nor Real projective space with fundamental group $C_2$. In particular $S^3/C_3$ is a lens space counterexample.

A vector field on the sphere means a tangent vector field. For $S^{2m}\subset {\mathbb{R}}^{2m+1}$ with $m\ge 1$, suppose $V(x)$ were continuous and nowhere zero. Normalize it to $w(x)=V(x)/\Vert V(x)\Vert$. Tangency means $⟨x,w(x)⟩=0$, and both vectors have norm one. Thereforehas squared norm ${\cos }^2(\pi t)+{\sin }^2(\pi t)=1$. It defines a homotopy through self-maps of the sphere from the identity to the antipodal map.

The identity has mapping degree $1$. The antipodal map on $S^n$ has degree $(-1{)}^{n+1}$: it is the boundary restriction of $-I$ on the oriented $(n+1)$-ball, and this linear map changes orientation by $\mathrm{sgn}\det (-I)=(-1{)}^{n+1}$. For $n=2m$ its degree is $-1$. Homotopy invariance of mapping degree would force $1=-1$, a contradiction. Thus every continuous tangent field on an even-dimensional sphere vanishes somewhere. For $S^0$ the tangent spaces are zero, so the conclusion holds directly. This proves the Hairy ball theorem using degree, rather than assuming it.

Identify the ambient Euclidean space ${\mathbb{R}}^{2m+2}$ with ${\mathbb{C}}^{m+1}$ and let $J$ be multiplication by $i$. Equivalently, in real coordinates,This linear map preserves the norm and satisfies $⟨x,Jx⟩=0$. On the unit sphere $S^{2m+1}$ defineIt is continuous, tangent because it lies in $x^{\perp }=T_x{S}^{2m+1}$, and has norm one everywhere. It is therefore a nowhere-zero section of the tangent bundle, including when $m=0$ and the sphere is $S^1$. This supplies a nonvanishing tangent field on an odd-dimensional sphere explicitly.

For the given harmonic potentials, ${\nabla }^2(\mathbf{\Phi }\cdot \mathbf{x})=2\nabla \cdot \mathbf{\Phi }$. Taking the divergence of the velocity representation therefore gives $\nabla \cdot \mathbf{u}=0$. Taking its Laplacian gives ${\nabla }^2\mathbf{u}=-2\nabla (\nabla \cdot \mathbf{\Phi })$, so the Stokes equations hold withThis is the Papkovich–Neuber representation in the sign and normalization of this problem. The potentials used in other conventions can differ by an overall factor and sign; the velocity and pressure must be scaled together.

(a) Point force. Put the singularity at the origin, $r=|\mathbf{x}|$, and choose $\mathbf{\Phi }=\mathbf{F}/(8\pi \mu r)$, $\chi =0$. The resulting Stokeslet isThe harmonic calculation applies away from the force. Distributionally this is the Green function for $-\nabla p+\mu {\nabla }^2\mathbf{u}+\mathbf{F}{\delta }^{(3)}(\mathbf{x})=0$, with zero divergence.

(b) Point source. Set $\mathbf{\Phi }=0$, $\chi =-Q/(4\pi r)$. The three-dimensional point source field isIts outward volume flux across every enclosing sphere is $Q$. The divergence is zero away from the source and equals $Q{\delta }^{(3)}(\mathbf{x})$ distributionally, so the singular point represents volume injection rather than an incompressible material point.

(c) Stresslet. Let $G_{ij}(\mathbf{x})=({\delta }_{ij}/r+x_i{x}_j/r^3)/(8\pi \mu )$ and use a symmetric trace-free stresslet tensor $S$ with the convention $u_i=-S_{jk}{\partial }_k{G}_{ij}$. Differentiating the Stokeslet and using symmetry and zero trace givesFor an axial force-dipole flow with axis $\mathbf{n}$, write $S=s(\mathbf{nn}-I/3)$. Then the radial factor is proportional to $s[3(\mathbf{n}\cdot \hat{\mathbf{x}}{)}^2-1]$. Specifying the derivative convention is essential for interpreting the signs of image strengths.

(d) Source dipole. A positive source and negative source separated by a vector $\mathbf{\ell }$ have moment $\mathbf{M}=Q\mathbf{\ell }$. Their dipole limit is minus the directional derivative of the unit source field. Hence the source dipole hasIt carries no net volume flux, and unlike a stresslet it is an irrotational potential flow.

Now set $\mathbf{r}=\mathbf{x}-d\mathbf{n}$ and $\mathbf{R}=\mathbf{x}+d\mathbf{n}$, with $d>0$. Define the unscaled kernels $K(\mathbf{R})=\mathbf{n}/R+(\mathbf{n}\cdot \mathbf{R})\mathbf{R}/R^3$, $H(\mathbf{R})=\mathbf{R}/R^3$, and let ${\partial }_n=\mathbf{n}\cdot {\nabla }_{\mathbf{R}}$. The no-slip image system of a normal Stokeslet isThe reflected Stokeslet has force $-F\mathbf{n}$. With the conventions above, the remaining image strengths areThe isotropic part of $S$ produces zero velocity, because ${\partial }_j{G}_{ij}=0$. Its inclusion simply puts the tensor in standard trace-free form. The pressure follows by superposition:

To verify the no-slip boundary condition explicitly, write $\mathbf{x}=\mathbf{\rho }$ on the plane, with $\mathbf{\rho }\cdot \mathbf{n}=0$, and let $b^2={\rho }^2+d^2$. ThenThe $b^{-5}$ terms cancel, while the remaining image contribution is $2d(\mathbf{R}-d\mathbf{n})/b^3=2d\mathbf{\rho }/b^3$. Thus both normal and tangential velocity vanish. All image singularities lie outside the fluid, so the physical half-space contains only the prescribed force.

For the far field, expand about the wall origin. The pair of opposite forces starts with $-2d{\partial }_{n}K$; the image stresslet cancels this entire $r^{-2}$ term. The first surviving field isWriting $z=\mathbf{x}\cdot \mathbf{n}$ and now $r=|\mathbf{x}|$, this isTherefore the generic fixed-angle far-field velocity decays as $r^{-3}$, rather than the $r^{-1}$ free-space Stokeslet or an uncancelled $r^{-2}$ stresslet. Along the wall it is exactly zero; at fixed height and large horizontal distance the leading tangential field instead scales as ${\rho }^{-4}$, consistent with the angular zeros of the displayed far field.

Let $\mathbf{n}$ point from fluid $-$ to fluid $+$, $P_s=I-\mathbf{nn}$, ${\nabla }_s=P_s\nabla$, and $\mathcal{K}={\nabla }_s\cdot \mathbf{n}$. A clean material interface without phase change has continuous velocity, with normal interface speed $V_n={\mathbf{u}}_{\pm }\cdot \mathbf{n}$. Its surface stress is $\gamma P_s$. The surface divergence of that tensor is ${\nabla }_s\gamma -\gamma \mathcal{K}\mathbf{n}$. Balancing forces on a small interfacial patch therefore gives the interfacial stress balance with variable surface tensionThe normal condition is $\mathbf{n}\cdot ({\mathbf{\sigma }}_{+}-{\mathbf{\sigma }}_{-})\mathbf{n}=\gamma \mathcal{K}$; the tangential condition is $P_s({\mathbf{\sigma }}_{+}-{\mathbf{\sigma }}_{-})\mathbf{n}=-{\nabla }_s\gamma$. On a stationary sphere with outward normal these reduce to $p_{-}-p_{+}=2\gamma /a$ when viscous normal stresses vanish. This fixes the sign convention for the Marangoni stress rather than leaving the direction ambiguous.

For the bubble take negligible internal dynamic viscosity, no contaminants, and initially the standard effectively insulating gas-bubble approximation. With heat advection neglected, the exterior temperature satisfies Laplace's equation, approaches $T_0+\mathbf{H}\cdot \mathbf{r}$ remotely, and has zero normal heat flux at $r=a$. The temperature dipole around an insulating sphere givesThe zero-flux condition follows by differentiating the radial factor. Let ${\gamma }_T=(d\gamma /dT{)}_{T_0}$. To first order in the weak gradient,

In the laboratory frame the force-free bubble can be described by the harmonic Papkovich–Neuber representation potentialsThey give constant exterior pressure and the decaying source dipole fieldAt the surface $\mathbf{u}\cdot \mathbf{n}=\mathbf{U}\cdot \mathbf{n}$, so the normal kinematic boundary condition holds for translation. In coordinates aligned with $\mathbf{U}$, $u_r=a^{3}U\cos \theta /r^3$ and $u_{\theta }=a^{3}U\sin \theta /(2r^3)$. Direct differentiation givesEquivalently, the surface tangential traction is $(3\mu /a)P_s\mathbf{U}$. The tangential stress condition therefore requiresFor the usual negative ${\gamma }_T$, the bubble moves toward hotter liquid. The normal condition is consistent as well: its viscous dipole is $-6\mu \mathbf{U}\cdot \mathbf{n}/a=3{\gamma }_T\mathbf{H}\cdot \mathbf{n}$, exactly the dipolar part of $2{\gamma }_s/a$. The constant part fixes $p_{\mathrm{gas}}-p_{\infty }=2{\gamma }_0/a$. There is no Stokeslet term, hence no net applied force or drag imbalance.

The thermal assumption matters. The problem gives the exterior thermal diffusivity but no gas thermal conductivity. For an inviscid bubble with conductivity ratio $K=k_{\mathrm{gas}}/k_{\mathrm{liquid}}$, temperature and normal heat-flux continuity instead giveso the same stress calculation givesThe insulating-bubble result uses $K\ll 1$. If instead an undisturbed surface temperature gradient is imposed, corresponding to $K=1$, the coefficient is $1/3$ rather than $1/2$. Exterior thermal diffusivity alone does not select between those thermal boundary conditions. This dependence is explicit rather than silently treating a missing gas property as a prescribed value.

The small parameters can now be stated precisely. With $U_{\ast }=a|{\gamma }_T||\mathbf{H}|/(2\mu )$ and $\nu =\mu /\rho$, requireSmall Reynolds number justifies Stokes flow; small Péclet number allows quasi-static heat conduction without advective distortion; small capillary number and fractional surface-tension variation support a nearly spherical interface and a weak perturbation of its base surface tension. These last two conditions are proportional for this velocity scale but describe different approximations. Linearizing $\gamma (T)$ also needs $a|\mathbf{H}||{\gamma }_{TT}/{\gamma }_T|\ll 1$ if its derivative varies appreciably. Negligible internal viscosity and the insulating specialization require the corresponding gas-to-liquid viscosity and conductivity ratios to be small. Gravity has been excluded; if restored, its buoyancy-induced migration must be small compared with $U_{\ast }$, not merely small enough to avoid deformation.

Finally, surfactants change both the constitutive surface tension and the surface stress dynamics. The surface flow driven by the thermal Marangoni effect carries insoluble surfactant toward the high-tension side, lowering its tension and producing an opposing concentration-driven gradient. When surface diffusion or exchange is slow, this tends to immobilize the interface, suppressing the clean-bubble thermocapillary migration of an insulating bubble, sometimes nearly stopping it. If the surfactant stays uniformly equilibrated, it can instead change ${\gamma }_T$ without a large opposing gradient. Its concentration, surface diffusivity and adsorption kinetics must be supplied to compute a replacement speed; no universal numerical reduction factor follows merely from its presence.