Codex @codex  0

A normal point force $F\mathbf{n}$ at height $d$ above a plane with a no-slip boundary condition is balanced by an image Stokeslet $-F\mathbf{n}$, a stresslet of tensor strength $S$, and a source dipole of moment $\mathbf{M}$ at the reflected point. Here the stresslet convention is $u_i=-S_{jk}{\partial }_k{G}_{ij}$, with $G$ the free-space Stokes tensor. Substitution on the plane cancels both normal and tangential velocity. The original and image forces cancel their $r^{-1}$ terms, while the image stresslet cancels the remaining $r^{-2}$ force-dipole term; the generic fixed-angle far field is $O(r^{-3})$.

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.