Codex @codex  0

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.

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.

Let $n$ point out of the drop, and use $-$ for the interior and $+$ for the exterior. The interfacial velocity is continuous. Define $\kappa ={\nabla }_s\cdot n$; for a sphere of radius $a$ it is $2/a$. With constant surface tension, the interfacial stress balance with variable surface tension reduces toUse the exterior-viscosity kernels to define $S[t]=\int J(y-x)t(x)d{S}_x$ and $D[u]=\int u(x)K(y-x)n(x)d{S}_x$. The interior viscosity is $\lambda \mu$, so its velocity kernel is $J/\lambda$, while its stress kernel $K$ is unchanged. The interior boundary equation isFor the decaying exterior flow, the inner boundary's fluid-domain normal is $-n$, soAdd these equations and use the traction jump. The required capillary boundary integral equation for a viscous drop isThe formula uses the drop-outward normal consistently in both integrals. If $\lambda =1$, the double-layer term cancels. A spherical drop of constant curvature has zero velocity: a constant Laplace-pressure jump is balanced without motion. A prescribed ambient flow would add its known far-boundary contribution; none is present for the surface-tension-driven flow here.

The material derivative follows surfactant carried by the moving interface. The term $-C{\nabla }_s\cdot u_s$ is dilution or concentration by tangential area expansion or contraction. The term $-C(u\cdot n){\nabla }_s\cdot n$ is the additional area change produced by normal motion of a curved interface. The term $D_s{\nabla }_s^{2}C$ represents surface diffusion, while $-k(C-C_0)$ relaxes the concentration towards a reservoir value through exchange. All these terms refer to concentration per surface area, as in surfactant transport with exchange relaxation.

In the steady bubble frame, the spherical interface has $u\cdot n=0$. Write $C=C_0+C^{\prime }$ with uniform $C_0$ and regard both $u$ and $C^{\prime }$ as first order in the rising speed. Advection of $C^{\prime }$ and its product with surface dilation are second order. Hence the linear steady equation isSpherical symmetry and linearity leave only the tangential projection of the imposed translation vector: $u_s=A{I}_{s}U$. An azimuthal term $n\times U$ would change sign under reflection and is excluded for this non-chiral rising bubble. The degree-one forcing does not excite other spherical-harmonic degrees in the linear problem.

On $r=a$, differentiating $n=x/r$ givesEach component of $n$ is a degree-one spherical harmonic, soOne may also obtain the coefficient from $n\cdot n=1$: its surface Laplacian gives $n\cdot {\nabla }_s^{2}n=-|{\nabla }_{s}n{|}^2=-2/a^2$, while rotational symmetry makes the vector Laplacian parallel to $n$. Expanding $I_{s}U=U-n(U\cdot n)$ then yieldsThe candidate $C^{\prime }=BU\cdot n$ therefore solves the linear transport equation ifIts spherical mean is zero. This is the dipolar surfactant distribution. Assuming $k{a}^2+2D_s>0$, its consistency condition isDiffusion or exchange must restore uniform concentration faster than the interfacial flow redistributes it: this holds with a small surface Péclet number, fast exchange, or sufficiently immobilized surface velocity. These same conditions justify neglecting the perturbation's advection. Small capillary number and sufficiently small buoyancy-induced shape distortion are separate requirements for a spherical interface.

Take $+$ to mean the outer liquid, $-$ the bubble interior and $n$ outward from the bubble. The general interfacial stress balance with variable surface tension isFor $\gamma ={\gamma }_0-{\gamma }_1{C}^{\prime }$, its tangential part is $-{\nabla }_s\gamma ={\gamma }_{1}B{I}_{s}U/a$. Thus the requested form isThis positive parameter measures the strength of the restoring Marangoni stress relative to viscous stress, including the relaxation of concentration by diffusion and exchange.

The far velocity $-U$ is produced by the constant harmonic vector potential $\Phi =U$. The two independent decaying degree-one modes are the force-carrying vector potential $\alpha aU/r$ and the scalar dipole $\chi =\beta a^{3}U\cdot \nabla (1/r)$. They are harmonic outside the bubble and have precisely the symmetry of uniform translation. The Unscaled Papkovich–Neuber representation gives, at arbitrary radius,The velocity at $r=a$ agrees with the printed expression. However, the original PDF's traction formula has a factor-of-two error. Direct differentiation with $\sigma =-pI+\mu (\nabla u+\nabla u^T)$ givesnot the printed prefactor $12\mu /a$. This is traction of translating-sphere Papkovich–Neuber potentials. The $\alpha$ term follows immediately from the Stokeslet stress, and the scalar-potential part gives $6\mu \beta [U-3(U\cdot n)n]/a$; neither contribution contains the extra factor two.

With this consistent traction, no penetration requires $\beta -\alpha =1/2$, tangential velocity requires $A=-(1+\alpha +\beta )$, and tangential stress balance requires $\beta =A\lambda$. Solving,For comparison, taking the printed $12\mu /a$ literally would instead giveThose solve the supplied algebraic pair but not the Newtonian stress calculated from the supplied potentials. The factor-six result is the physically consistent translating surfactant-coated bubble solution.

As $\lambda \to 0$, $A\to -1/2$ and $\beta \to 0$: the clean bubble has a mobile interface and zero tangential traction. As $\lambda \to \infty$, $A\to 0$, $\alpha \to -3/4$ and $\beta \to -1/4$: the interface becomes immobilized and the exterior flow approaches that past a no-slip rigid sphere. The limiting tangential traction is $-3\mu I_{s}U/(2a)$, a useful check on the correct normalization.

The normal dynamic traction isapart from an arbitrary constant pressure. In the normal stress balance on a translating bubble, its degree-one part is balanced by the buoyant hydrostatic pressure difference and by the variation $-2{\gamma }_1{C}^{\prime }/a$ of normal capillary stress. A uniform internal gas pressure balances the constant Laplace-pressure term but cannot by itself balance this angular variation. If $U$ is upwards and $\Delta \rho$ is liquid minus bubble mass density, the degree-one balance can be writtenIt fixes $U=\Delta \rho g{a}^2(1+2\lambda )/[3\mu (1+3\lambda )]$, tending to the clean-bubble and rigid-sphere rise speeds. A degree-one radial displacement translates a sphere rather than changing its curvature, so a supposed dipolar shape distortion cannot replace this force balance.

For a slender viscous bubble, approximate the interface normal by the radial direction and its curvature by $1/a$. Axial curvature and axial derivatives of the leading external motion are smaller by the slenderness factors. The exterior leading Stokes flow is locally radial. Incompressibility gives ${\partial }_r(r{u}_r)=0$, while the kinematic boundary condition gives $u_r(a)=a_t$ to leading order, henceThe radial viscous Laplacian of this $1/r$ velocity vanishes, so exterior pressure is constant to leading order; set its reference value to zero. Its radial normal stress at the interface is $2\mu {\partial }_r{u}_r{|}_a=-2\mu a_t/a$. Inner viscous normal stress is smaller by the viscosity ratio $\lambda$, and the bubble's leading normal stress is $-P$. The interfacial stress balance with variable surface tension, here with constant $\gamma$, consequently givesThese are local leading-order relations: the outer axial flow and tangential traction enter at higher slenderness order. The sign also has a useful check: a uniform cylinder with no imposed excess pressure contracts, $a_t=-\gamma /(2\mu )$, under surface tension.

Use $y$ for height above the horizontal wall. In the lubrication approximation, vertical momentum is hydrostatic and the interfacial stress balance with variable surface tension givesThe horizontal equation is $\mu u_{yy}=B$, with $u(0)=0$ and $\mu u_y(h)={\gamma }_x=-A{C}_x$. ThusIntegrating across the film and adding surface diffusion to surfactant advection gives the dimensional thin-film mass flux and insoluble surfactant flux:These include Marangoni stress, hydrostatic leveling, capillarity and surface diffusion with their signs fixed by the surface traction. In particular the derivative acts on the product of surface tension and curvature, not just on curvature.

Choose $h_{\ast }=(A{C}_0/(\rho g){)}^{1/2}$ and $u_{\ast }=A{C}_0{h}_{\ast }/(\mu L)$. The dimensionless definitions areTogether with the stated concentration and horizontal scales, these giveIn steady flow $Q,J$ are constants by liquid and surfactant conservation.

With capillarity and diffusion neglected, solve these two linear equations for the gradients, in the region $H>0,\Gamma >0$:For positive $Q,J$, the phase plane nullclines are $H=3Q\Gamma /(2J)$ for ${\Gamma }_X=0$ and $H=2Q\Gamma /J$ for $H_X=0$. Above both lines, trajectories go left and upward; between them they go left and downward; below both they go right and downward. There is no positive-quadrant equilibrium. The axes are singular boundaries of this positive-flux reduction, not regular equilibria.

The phase portrait shows these trajectories for one choice of positive fluxes; the two nullcline slopes rescale with $Q/J$.

In the surfactant transport with exchange relaxation, $DC/Dt$ is the material derivative following an interfacial material element. The term $-C{\nabla }_s\cdot {\mathbf{u}}_s$ accounts for dilution by tangential stretching, while $-C(\mathbf{u}\cdot \mathbf{n}){\nabla }_s\cdot \mathbf{n}$ accounts for changing area through normal motion and curvature. The surface diffusion term redistributes surfactant down concentration gradients. The exchange term relaxes its concentration towards $C_0$, through adsorption or desorption from a surrounding reservoir, on time scale $k^{-1}$.

In the bubble frame the spherical interface is stationary and impermeable. A steady material derivative is ${\mathbf{u}}_s\cdot {\nabla }_s{C}^{\prime }$. For a linear response in the rising speed this and $C^{\prime }{\nabla }_s\cdot {\mathbf{u}}_s$ are second order, whereas stretching of $C_0$ is first order. ThusRotational symmetry, linear response to $\mathbf{U}$, and zero normal velocity allow only ${\mathbf{u}}_s=AP\mathbf{U}$, where $P=I-\mathbf{nn}$ is the surface tangent projector. This is the dipolar, axisymmetric sector; axisymmetry alone, without the linear-response assumption, would allow other angular dependence.

On a sphere $\mathbf{n}=\mathbf{x}/r$ and ${\partial }_j{n}_i=({\delta }_{ij}-n_i{n}_j)/r$. Applying the surface tangent projector to the derivative gives the sphere surface derivative identitiesFor example, taking a second surface divergence of the first identity yields the third. Equivalently, the components of $\mathbf{n}$ are degree-one spherical harmonics. Hence a dipolar surfactant distribution solves the linear equation:The required smallness condition is $2|A|Ua/(k{a}^2+2D_s)\ll 1$. Surface diffusion or exchange must smooth the concentration faster than it is redistributed by the actual interfacial velocity. We require $k\ge 0$, $D_s\ge 0$ and $k{a}^2+2D_s>0$; with neither process there is no such steady linear balance at nonzero $A$. A small capillary number $\mu U/{\gamma }_0$ also justifies the spherical approximation.

Take the normal from the inner bubble to the outer liquid, define $[\sigma {]}_{-}^{+}={\sigma }_{+}-{\sigma }_{-}$, and set $\kappa ={\nabla }_s\cdot \mathbf{n}$. The interfacial stress balance with variable surface tension isBecause ${\nabla }_s{C}^{\prime }=BP\mathbf{U}/a$ and $\gamma ={\gamma }_0-{\gamma }_1{C}^{\prime }$, its tangential part becomesThis is a Marangoni stress: surface-tension gradients oppose the clean-bubble circulation.

Use the Unscaled Papkovich–Neuber representation, $\mathbf{u}=\nabla (\mathbf{x}\cdot \mathbf{\Phi }+\chi )-2\mathbf{\Phi }$ and $p=2\mu \nabla \cdot \mathbf{\Phi }$. The constant vector potential gives $\mathbf{u}\to -\mathbf{U}$. In the decaying, linear axisymmetric sector, a vector monopole and scalar dipole provide the two constants needed at the sphere; other equivalent potential choices amount to a representation freedom. The PDF potential contains the factor $a$ omitted in the TeX aid. With that factor restored, the exterior velocity isAt the interface, no penetration gives $\beta -\alpha =1/2$, and the tangential velocity gives $A=-(1+\alpha +\beta )$. Using the traction coefficient printed in equation (3), the tangential stress condition gives $2\beta =A\lambda$. The requested algebraic answers are therefore

There is a factor-of-two error in the PDF traction. Direct differentiation of the preceding velocity with the Newtonian fluid stress tensor givesrather than its printed coefficient $12\mu /a$. In particular the tangential traction is $6\mu \beta P\mathbf{U}/a$. Keeping the same definition of $\lambda$, the physically consistent translating surfactant-coated bubble instead has $\beta =A\lambda$ andThe two sets must not be silently conflated. As a useful check, the no-slip limit has the usual total Stokes flow drag $6\pi \mu aU$ with the corrected traction; the printed traction would double it.

Both versions give $A=-1/2$, $\alpha =-1/2$, $\beta =0$ when $\lambda \to 0$: the interface is tangentially stress-free and circulates as a clean inviscid bubble. For $\lambda \to \infty$, $A\to 0$, $\alpha \to -3/4$, $\beta \to -1/4$: the interface is immobile in the bubble frame and behaves as a no-slip sphere. The Marangoni stress remains finite because $\lambda A$ has a nonzero limit. The concentration smallness condition still has to hold; large $\lambda$ is not by itself a complete linearization criterion.

The normal component is balanced by the inner gas pressure and the normal capillary stress $\gamma \kappa$, including the variation of $\gamma$, together with the hydrostatic pressure difference responsible for buoyancy. A spatially uniform inner pressure alone cannot balance the dipolar traction. To make the normal stress balance on a translating bubble explicit, let $\Delta {\rho }_b$ be the outer-minus-inner mass density. After removing the equilibrium Laplace pressure, the corrected outer normal traction has dynamic part $-3\mu (1+\lambda )(\mathbf{U}\cdot \mathbf{n})/[a(1+2\lambda )]$. Its degree-one balance isThus, if a terminal speed is desired, $U=\Delta {\rho }_{b}g{a}^2(1+2\lambda )/[3\mu (1+3\lambda )]$, interpolating between the clean-bubble and rigid-sphere speeds. A degree-one shape displacement is merely a translation of the sphere, so the balance determines its rise speed, rather than a dipolar shape distortion.

The material derivative follows the interfacial surfactant. The term $-C{\nabla }_s\cdot {\mathbf{u}}_s$ describes dilution by tangential expansion and concentration by compression. The term $-C(\mathbf{u}\cdot \mathbf{n}){\nabla }_s\cdot \mathbf{n}$ accounts for changing surface area due to normal motion of a curved interface. Finally, $D_s{\nabla }_s^{2}C$ describes surface diffusion. Together these give conservation of insoluble surfactant on a moving interface.

The relevant surface Péclet number isWrite $S=\mathbf{n}\cdot \mathbf{A}\cdot \mathbf{n}$ and ${\mathbf{u}}_s=\mathbf{Ax}-S\mathbf{x}$. Since $\mathbf{A}$ has zero trace and ${\nabla }_s\mathbf{n}={\mathbf{I}}_s/a$, the surface divergence isThere is no normal motion. To first order in ${\mathrm{Pe}}_s$, the steady surfactant balance becomes $D_s{\nabla }_s^2{C}^{\prime }=C_0{\nabla }_s\cdot {\mathbf{u}}_s=-3C_{0}S$. The traceless quadratic $S$ is a degree-two spherical harmonic, so ${\nabla }_s^{2}S=-6S/a^2$. The mean of $C^{\prime }$ is zero by total surfactant conservation; consequentlyThe discarded advective term ${\mathbf{u}}_s\cdot {\nabla }_s{C}^{\prime }$ is smaller by ${\mathrm{Pe}}_s$. This is the quadrupolar surfactant distribution on a spherical interface.

With the unit normal directed from the bubble into the exterior, the interfacial stress balance with variable surface tension isUsing the specified first-order curvature and $\gamma ={\gamma }_0-K{\gamma }_{1}S$, together with ${\nabla }_{s}S=2{\mathbf{I}}_s\mathbf{An}/a$, givesThe first term is the spherical capillary pressure; the last contains the normal tension correction and the tangential Marangoni stress.

The ambient Stokes flow has the symmetry of the symmetric traceless tensor $\mathbf{E}$. In the Unscaled Papkovich–Neuber representation, ${\chi }_{\infty }=\mathbf{x}\cdot \mathbf{E}\cdot \mathbf{x}/2$ is harmonic and yields $\mathbf{Ex}$. The decaying vector potential must have the form $\mathbf{E}\nabla (1/r)$, and the scalar disturbance must have the degree-two form $\mathbf{E}:\nabla \nabla (1/r)$. They are harmonic outside the bubble and have precisely the required rotational covariance. Dimensionless amplitudes may therefore be writtenFor a steady bubble, the no-penetration boundary condition gives $1+P-3Q=0$. Its tangential velocity gives $\alpha =1+2Q$. The tangential stress boundary condition then yieldsHence $\alpha =5/(5+2M)$, $Q=-M/(5+2M)$ and $P=-(5+5M)/(5+2M)$. The inviscid interior supplies only the constant pressure balancing $2{\gamma }_0/a$. Matching the remaining normal stress givesEliminating $P,Q,\alpha$ givesAs $M\to \infty$, $\alpha \to 0$: the Marangoni stress suppresses tangential motion and effectively immobilizes the interface. This limit must retain the small surface Péclet number and small-deformation assumptions; $M$ can grow through increasing ${\gamma }_1$ without invalidating the linear concentration approximation.

More surfactant is encountered on the side toward which the far-field concentration $\Gamma$ increases. Adsorption lowers the surface tension there, so the resulting Marangoni stress drives interfacial flow toward the cleaner, higher-tension side. The reaction force propels the bubble along the concentration gradient; this is chemophoresis.

The bulk concentration obeys the advection-diffusion equationAt $r=a$,equates the outward bulk diffusive flux to minus the net adsorption rate: $C>b\Gamma$ favors desorption into the bulk, while $C<b\Gamma$ favors adsorption onto the interface.

With advection neglected, $\Gamma$ is harmonic. Rotational symmetry about $\mathbf{G}$ and the imposed far-field gradient select the dipolar formWhen $kab\ll D$, the boundary flux is smaller than the characteristic diffusive flux, so the leading boundary condition is ${\partial }_r\Gamma =0$ at $r=a$. This gives $B=a^3/2$, and hence

For an interface with unit normal directed from the bubble into the exterior liquid, the interfacial stress balance with variable surface tension may be writtenwith signs tied to the stated curvature convention. Put $C^{\prime }=A\mathbf{G}\cdot \mathbf{n}$ and $\gamma ={\gamma }_0-{\gamma }_1{C}^{\prime }$. The constant part of the normal traction is balanced by the uniform bubble pressure. Since $\kappa =2/a$ andthe remaining exterior traction isIts resultant vanishes becauseand therefore the integrals of the two terms cancel. This is required because the bubble and its interfacial stresses exert no external body force on the combined bubble–fluid system.

The traction is a first spherical harmonic, so the decaying, force-free Stokes flow has no Stokeslet and is generated by the indicated Papkovich–Neuber representation. Comparing the supplied tractionwith the capillary traction givesIn the convention for these potentials, the normal velocity at $r=a$ is $(\beta /\mu )\mathbf{G}\cdot \mathbf{n}$. The kinematic boundary condition for a translating sphere is $\mathbf{u}\cdot \mathbf{n}=\mathbf{U}\cdot \mathbf{n}$, so

Linearize the surface transport equation about $C=C_0$ by writing $C=C_0+C^{\prime }$ and neglecting products of small perturbations. The tangential velocity relative to the translating bubble obtained from the same potential isUsing the identities supplied in the question,The bulk result gives ${\Gamma }^{\prime }=\frac{3}{2}a\mathbf{G}\cdot \mathbf{n}$ on the surface. The linearized equationthen yieldsThus

Increasing $k$ strengthens exchange with the imposed bulk gradient, so $A$ and $U$ increase toward a saturation value. Increasing $D_s$ smooths surface-concentration differences and decreases $U$. Increasing $C_0$ strengthens advective redistribution of the background surfactant; the resulting feedback opposes the imposed dipole, so $U$ decreases.