Codex @codex  0

Write $K(u,f)=\Vert u{\Vert }_{C^2(B_1)}+\Vert f{\Vert }_{C^{0,\alpha }(B_1)}$; the paper's $|u{|}_{2;B_1}$ denotes a $C^2$ norm, not an $L^2$ norm. We use a blow-up compactness proof of an interior Schauder estimate. If the estimate failed for fixed $\delta >0$, normalization would produceChoose $x_j,y_j\in B_{1/2}$ with $d_j=|y_j-x_j|$ and $|D^2{u}_j(y_j)-D^2{u}_j(x_j)|\ge \frac{1}{2}{d}_j^{\alpha }$. Since $\Vert D^2{u}_j{\Vert }_{\infty }\to 0$, necessarily $d_j\to 0$. Subtract the quadratic Taylor polynomial $P_j$ at $x_j$ and setThe domains expand to ${\mathbb{R}}^n$. At zero, $w_j,D{w}_j,D^2{w}_j$ vanish; the Hölder seminorm of $D^2{w}_j$ is at most ${\delta }^{-1}$. Integration along segments gives uniform $C^{2,\alpha }$ bounds on each fixed ball. Alsolocally uniformly, since $[f_j{]}_{\alpha }\to 0$. By compact embedding of Hölder spaces and a diagonal subsequence, $w_j$ converges locally in $C^{2,\beta }$ to an entire harmonic function $w$, with $[D^{2}w{]}_{\alpha }\le {\delta }^{-1}$ and $D^{2}w(0)=0$. The unit vectors $(y_j-x_j)/d_j$ have a convergent subsequence; its limit $e$ satisfies $|D^{2}w(e)|\ge 1/2$.

Each second derivative of $w$ is an entire harmonic function with finite global Hölder seminorm. The supplied polynomial-growth Liouville theorem for harmonic functions, at growth exponent $\alpha <1$, makes it constant, and its value at zero makes it zero. This contradicts the nonzero Hessian at $e$, proving the estimate with the small $\delta$ term.

To remove that term, rescale the estimate to all contained balls and apply the Simon absorption lemma. Explicitly, with $d(x)=1-|x|$ and $d_{xy}=\min (d(x),d(y))$, defineFor $|x-y|<d_{xy}/8$, apply the rescaled estimate on $B_{d(x)/2}(x)$; its larger-ball seminorm is controlled by $M$ because all points there have boundary distance at least $d(x)/2$. For more separated pairs use $2\Vert D^{2}u{\Vert }_{\infty }$. The result is $M\le C_0\delta M+C_{\delta }K(u,f)$. Choose $C_0\delta <1/2$ and use $d_{xy}\ge 1/2$ on $B_{1/2}$ to obtainThe forcing's Hölder seminorm cannot be discarded. For integers $k\to \infty$, takeTheir $C^2$ and forcing supremum norms stay bounded, but the $11$ entry of the Hessian matrix differs by two at $x_1=\pm \pi /(2k)$, giving $[D^2{u}_k{]}_{\alpha ;B_{1/2}}\ge 2(k/\pi {)}^{\alpha }\to \infty$.

Use the usual connected-manifold convention for this assertion. Suppose, for contradiction, that $H^2(M,{\mathcal{O}}_M)=0$. The Dolbeault theorem gives $H^{0,2}=0$, and the Hodge decomposition theorem for compact Kähler manifolds and its conjugation symmetry give $H^{2,0}=0$. Hence every real class in $H^2(M,\mathbb{R})$ has a real harmonic differential form of type $(1,1)$ as its representative.

The intersection form is $Q([u],[v])={\int }_{M}u\wedge v$. Its value on the Kähler class is positive: $Q([\omega ],[\omega ])={\int }_M{\omega }^2>0$. Every real harmonic $(1,1)$ form splits uniquely as $c\omega +u_0$, with $u_0$ primitive and $c$ constant. Indeed, $\Lambda u$ is a harmonic function, hence constant by connectedness, and $\Lambda \omega =2$; subtracting $c\omega$ kills $\Lambda u$, and part (b) makes $u_0$ primitive.

Part (c) gives $\ast u_0=-u_0$, while $\omega \wedge u_0=0$. ThusThis proves the Hodge index theorem for compact Kähler surfaces in this case: $Q$ has signature $(1,b_2-1)$.

On the alleged two-dimensional totally isotropic subspace $V$, the linear functional $v\mapsto Q(v,[\omega ])$ has a nonzero vector in its kernel. Its harmonic differential form representative is primitive, so the preceding negative-definiteness statement gives $Q(v,v)<0$, contradicting the vanishing of $Q$ on $V$. ThereforeConnectedness is essential if one's definition of a manifold allows disconnected spaces. For $M=Y\bigsqcup Y$, with $Y={\mathbb{CP}}^1\times {\mathbb{CP}}^1$, one has $H^2(M,{\mathcal{O}}_M)=0$. Take on each component the pullback of the degree-two class from its first factor. Each class has square zero, and their mutual pairing is zero because their supports lie on different components. Their span is then a two-dimensional totally isotropic subspace, disproving the unqualified disconnected version.

Let $U_t=(1,0,0)+(B_t^1,B_t^2,B_t^3)$, with three independent standard Brownian motions, and set $R_t=|U_t|$ and $M_t=R_t^{-1}$. The radial stochastic process $R$ is a three-dimensional Bessel process, which never hits zero by the Hitting-zero classification for a Bessel process. In three dimensions $u\mapsto |u{|}^{-1}$ is a harmonic function away from zero: direct differentiation gives $\Delta (|u{|}^{-1})=0$.

For the stopping times $T_n=\inf \{t:R_t\notin (1/n,n)\}$, $n\ge 2$, the Itô formula givesThe coefficients are bounded by $n^2$, so each stopped Itô integral is a true martingale on every finite horizon. The stopping times form a localizing sequence; as permitted, its divergence need not be proved here. Thus $M$ is a nonnegative continuous local martingale.

To see that $M$ is not a true martingale, integrate the three-dimensional multivariate normal density of $U_t$ in spherical coordinates. The angular integral yieldswhere $\Phi$ is the standard normal distribution function. Therefore the reciprocal of a three-dimensional Bessel process is a strict local martingale.

Set $F(w)=(\varphi (w)-\varphi (0))/w$ for $w\ne 0$. The singularity is removable, with $F(0)={\varphi }^{\prime }(0)$. Injectivity of the conformal map makes $F(w)\ne 0$ away from zero, and conformality makes $F(0)\ne 0$ as well. Since the disc is simply connected, $F$ has a holomorphic logarithm. Therefore $\psi =\log |F|$ is a harmonic function, including at zero, and $\psi (0)=\log |{\varphi }^{\prime }(0)|$.

The mean value property for harmonic functions gives, for every $0<r<1$,For general $D$, boundary values mean radial limits, rather than a continuous extension of $\varphi$ to every point of the circle. The boundary logarithmic mean of a univalent function justifies taking $r\uparrow 1$: the Koebe distortion theorem bounds $|F|$ below by $|{\varphi }^{\prime }(0)|/4$, while the standard integral-mean bound for a univalent function, ${\sup }_{r<1}\int |\varphi (r{e}^{i\theta }){|}^{p}d\theta <\infty$ for $0<p<1/2$, gives uniform integrability of its positive logarithm. Radial limits exist almost everywhere, and passage to the integral follows. Thus, writing $dw=d\theta$ for arc length on the unit circle,The expression on the left is the logarithm of the conformal radius of $D$ at $z=\varphi (0)$.

The canonical Dirichlet Green function on the unit disc with its pole at zero is $-\log |w|$. Conformal invariance of the planar Green kernel therefore givesComparing with $G_D(z,y)=-\log |z-y|-G_z(y)$ yieldsThis also explains the boundary data: when $|w|=1$, the correction is $-\log |\varphi (w)-z|$. The composition with ${\varphi }^{-1}$ is harmonic by conformal invariance of harmonicity. Its apparent singularity at $z$ is removable because the quotient defining $\psi$ extends to ${\varphi }^{\prime }(0)$.

Evaluating there proves the regular part of the planar Dirichlet Green function formulaOn an unbounded domain, boundary values alone need not specify a unique harmonic function. Here $G_z$ is the correction belonging to the canonical Dirichlet Green kernel, which fixes that ambiguity.

A Brownian excursion in the upper half-plane from $0$ to infinity is killed planar Brownian motion conditioned to stay in $\mathbb{H}$, with its entrance law at $0$. Precisely, it is the Doob h-transform for $h(z)=\mathrm{Im}z$, obtained by starting at $i\epsilon$ and letting $\epsilon \downarrow 0$. Equivalently,where $W$ is standard Brownian motion and the independent process $R$ is a dimension-three Bessel process started at zero. The transformed generator is $\frac{1}{2}\Delta +y^{-1}{\partial }_y$, which explains this representation.

First compute avoidance from an interior point $z$. Let $\tau$ be the exit time of ordinary planar Brownian motion from $\mathbb{H}\setminus A$. The harmonic correctionhas boundary values $\mathrm{Im}z$ on $A$ and zero on the real axis. The hydrodynamic normalization at infinity gives $g_A(z)=z+\mathrm{hcap}(A)/z+O(|z{|}^{-2})$. Together with the maximum principle for harmonic functions this gives the useful bound $0\le v\le {\sup }_{a\in A}\mathrm{Im}a$. Since ordinary Brownian motion hits the real axis almost surely, the optional stopping theorem applied to this bounded harmonic function givesThe Doob h-transform weights a stopped path by its final height divided by its initial height. Applying this at the first hit of $A$, with bounded-time localization and then a limit, therefore yieldsFor the entrance law the same result follows by conditioning at a small positive time and letting that time decrease to zero: the avoidance function has the boundary limit computed next, and the excursion is initially in a neighbourhood disjoint from $A$.

Because the closure of the hull avoids $0$, the Schwarz reflection principle makes $g_A$ analytic in a neighbourhood of zero, with real $g_A(0)$ and positive real $g_A^{\prime }(0)$. ThusIn fact, the reflected Taylor expansion gives $\mathrm{Im}{g}_A(z)/\mathrm{Im}z=g_A^{\prime }(0)+O(|z|)$ uniformly as $z\to 0$ inside the half-plane, so the same limit applies to the entrance process. The map in the question is ${\psi }_A=g_A-g_A(0)$, so its derivative agrees with $g_A^{\prime }$. We obtain the restriction probability of a Brownian half-plane excursion:

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.

Take $x$ in the injection direction and $y$ along the flat interface. The specified $U$ must be a Darcy velocity, so the pore-fluid velocity and unperturbed interface speed are $U/\varphi$. Write the perturbed interface as $x=X(t)+\eta (y,t)$, with $\dot{X}=U/\varphi$, and choose the normal from fluid 1 into fluid 2. The background pressure gradients from Darcy law are ${\partial }_x{p}_i^{(0)}=-{\mu }_{i}U/k$.

By incompressibility and uniform permeability of a porous medium, pressure perturbations are harmonic functions. A transverse Fourier mode with positive wavenumber $\alpha$ has the decaying formsThe kinematic boundary condition and equal normal Darcy flux implyWithout capillarity the pressure is continuous at the displaced interface. Its linear perturbation is thereforeConsequently the planar viscous-fingering dispersion relation reduces toThis is the Saffman–Taylor instability: less viscous fluid advances preferentially through protrusions. Without a short-wave regularizing mechanism the idealized rate has no finite maximum.

For the prescribed apparent capillary pressure, the chosen normal is $\mathbf{n}=(1,-{\eta }_y)/\sqrt{1+{\eta }_y^2}$, so $\nabla \cdot \mathbf{n}=-{\eta }_{yy}+O({\eta }^2)={\alpha }^2\eta +O({\eta }^2)$. Subtract the flat-front pressure jump $\gamma (1+\beta X){\kappa }_0$. Linearization of the jump at $x=X+\eta$ givesThere is no first-order product of the surface-tension perturbation and the perturbed macroscopic curvature. ThusUsing the capillary number convention $\mathrm{Ca}={\mu }_{2}U/\gamma$ and ${\mu }_2/({\mu }_1+{\mu }_2)=(1+M)/2$, this isIf $x$ is measured from the unperturbed front and $\gamma$ denotes its local apparent tension, set $X=0$ to obtain the usual fixed-coefficient expression. For a physically fixed spatial gradient, $X=Ut/\varphi$ and the displayed rate is an instantaneous, generally time-dependent growth rate; the amplitude solves $\dot{a}=\sigma (\alpha ;X(t))a$, rather than one constant exponential law for the entire displacement. The parameter $\beta$ has dimensions of inverse length; the absolute tension gradient is $\gamma \beta$.

Assume $1+\beta X>0$, as required for positive apparent surface tension, and define $A=({\mu }_2-{\mu }_1)U-k\gamma \beta {\kappa }_0$. If $A>0$, unstable modes satisfy $0<{\alpha }^2<A/[k\gamma (1+\beta X)]$. Differentiating the cubic dispersion relation gives the fastest-growing modeWhen $A\le 0$, every nonzero mode decays, with a neutral zero-wavenumber translation in the infinite-domain limit. For $\beta {\kappa }_0>0$, surface-tension-gradient stabilization of viscous fingering suppresses the instability entirely at positive injection speeds satisfyingIf $\beta {\kappa }_0\le 0$, no positive injection speed eliminates every long-wave unstable mode in an unbounded interface. The prescribed wetting variation is treated as an effective normal pressure law for porous-media flow.

For a spherically symmetric function $\Phi (r)$, the Laplacian in spherical coordinates isIntegrating from zero to $r$ and using regularity at the origin gives

For the stated density, integration and the boundary condition $\Phi (a)=0$ giveThe two pieces and their first derivatives agree at $r=b$. If another solution existed, its difference from $\Phi$ would be harmonic with zero boundary data, so part (a) proves Uniqueness of the Dirichlet problem.

On the shell $U(b,a)$, normalize the outer radial part to obtain the harmonic functionBy the Dirichlet principle, this function minimizes the energy among all $w$ with the same boundary values. Sinceits energy isTherefore

Fix $B=B(x_0,r)⋐\Omega$. Solvability of the Dirichlet problem on a ball gives a unique harmonic function $h\in C^2(B)\cap C^0(\overline{B})$ with $h=u$ on $\partial B$. Both $u$ and $h$ have the mean value property for harmonic functions, so $v=u-h$ has it as well and vanishes on $\partial B$.

If $v$ were not zero, compactness of $\overline{B}$ would give either a positive maximum or a negative minimum in the interior. Part 1(i) would make $v$ constant, contradicting its zero boundary values. Hence $u=h$ on $B$. Every point lies in such a ball, so

Let the exterior ball be $B(y,R)$, tangent at $x_0$. Then $|x-y|\ge R$ on $\overline{\Omega }$, with equality only at $x_0$. Choose $\delta =R^{-1}$ and defineIn two dimensions, $\log |x-y|$ is a harmonic function away from $y$, so $\Delta w=0\le 0$ in $\Omega$. Moreover $w(x_0)=0$ and $w>0$ on $\overline{\Omega }\setminus \{x_0\}$. Thus $w$ is a barrier for the Dirichlet problem, and $\boxed{x_0\text{ is a regular boundary point}}$.

Fix a closed ball $\overline{B}(x,r)\subset D$ and let $\tau$ be its first exit time. By the Strong Markov property, conditioning at $\tau$ givesThe orthogonal invariance of Brownian motion implies that $X_{\tau }$ is uniformly distributed on the sphere $\partial B(x,r)$. HenceThus $\varphi$ has the mean value property on every ball compactly contained in $D$. Since $0\le \varphi \le 1$, the mean-value characterization of harmonic functions yieldsThis function is the harmonic measure of $A$ viewed from $x$.

The Wirtinger derivative $u_z=(u_x-i{u}_y)/2$ is holomorphic, since ${\partial }_{\bar{z}}{u}_z=\Delta u/4=0$ by the Laplace equation. The conformal map $w=z^2$ takes the quadrant to the complex upper half-plane. Integrating the prescribed tangential boundary derivatives, choose a common corner value $c$ and define the real boundary functionThe decaying data are bounded under the usual regularity at the corner, so $f(s)=O(1+\sqrt{|s|})$. Thus its Poisson integral converges. A holomorphic function whose real part is that Poisson integral is supplied by the regularized Schwarz integral formula:For $u_{\ast }(z)=\mathrm{Re}W(z^2)$, the chain rule gives $(u_{\ast }{)}_z=z{W}^{\prime }(z^2)$. Apply integration by parts to obtainSubstituting $s=r^2$ and $s=-r^2$ gives the requested particular integral representation:Both integrals converge absolutely for an interior point of the quadrant. The Sokhotski–Plemelj formula gives $\mathrm{Re}(u_{\ast }{)}_z(x)=g_2(x)/2$ on the horizontal edge and $\mathrm{Im}(u_{\ast }{)}_z(iy)=-g_1(y)/2$ on the vertical edge, verifying the prescribed derivatives.

The printed conditions alone do not determine $u_z$ uniquely. The general answer adds a holomorphic function $H$ satisfyingEquivalently, write $H(z)=ih(z^2)/z$, where $h$ is holomorphic in the complex upper half-plane with real boundary values on the nonzero real axis. An antiderivative of $2H$ has a harmonic real part with zero tangential boundary derivatives. For example, $u=xy$ contributes $H=-iz/2$ without changing either datum. If the two edge constants differ by $c_2-c_1$, the angular harmonic function contributes $i(c_2-c_1)/(\pi z)$. The boxed expression selects the Poisson integral representative with equal edge constants; regularity and suitable growth conditions can be used to select that representative.

Use the Poisson integral to choose a smooth harmonic function $G$ on the unit disk whose boundary value is $g$. Seek $v$ in the form $v=w+G$, which is the standard homogenization of Dirichlet boundary data. Then $w\in H_0^1(\Omega )$ and, since $\Delta G=0$,Part (d) gives a unique such $w$, and henceis the unique weak solution with boundary value $g$.

SetThen $g$ is entire and never takes a real value. The continuous function $\mathrm{Im}g$ is nowhere zero; because $\mathbb{C}$ is connected, it has one sign everywhere. Thus either $\mathrm{Im}g$ or $-\mathrm{Im}g$ is a positive harmonic function. Part (i) makes it constant. The Cauchy-Riemann equations then force the real part of $g$ to be constant as well. Hence $g$ and therefore $f$ are constant.

On $\mathbb{H}\setminus (A\cup C)$, the nonnegative harmonic function $h_A+h_C$ dominates the boundary data $\mathrm{Im}z$ on $A\cup C$: at a point of $A$, for example, $h_A(z)=\mathrm{Im}z$ and $h_C(z)\ge 0$. The maximum principle for harmonic functions, or equivalently Brownian motion stopped on the union, givesComparing the coefficients at infinity proves

Set $h_A(z)=\mathrm{Im}(z-g_A(z))$. This is a nonnegative harmonic function on $\mathbb{H}\setminus A$, has boundary value $\mathrm{Im}z$ on the hull boundary, and tends to zero on the real boundary. If $\tau$ is the first exit time of planar Brownian motion from $\mathbb{H}\setminus A$, the optional stopping theorem givesThe hydrodynamic expansion yields $h_A(iy)=\mathrm{hcap}(A)/y+O(y^{-2})$, and therefore the Brownian representation of half-plane capacity

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