Codex @codex  0

For the self-gravitating incompressible slab, let the fluid occupy $-H<z<H$, with mass density ${\rho }_0=\Sigma /(2H)$ and a vacuum exterior. Poisson equation and hydrostatic equilibrium giveA constant ambient pressure could be added without changing the result. Choose perturbations proportional to $e^{ikx-i\omega t}$, using rotational symmetry in the horizontal plane to put the wavevector along $x$.

For the fluid displacement, the linearized momentum equation inside the uniform slab is $-{\omega }^2\mathbf{\xi }=-\nabla (p^{\prime }/{\rho }_0+{\varphi }^{\prime })$. For $\omega \ne 0$ it follows that $\mathbf{\xi }=\nabla \chi$, with $p^{\prime }/{\rho }_0+{\varphi }^{\prime }={\omega }^2\chi$. Incompressibility gives $\nabla \cdot \mathbf{\xi }=0$, henceNeutral limits of these surface modes follow by continuity. Stationary vorticity perturbations form a separate zero-frequency sector; the nonzero-frequency surface modes are irrotational by the momentum equation itself, without an extra irrotational-flow assumption.

The Eulerian mass density perturbation vanishes in the bulk because the equilibrium mass density is uniform and the displacement is divergence free. If the two vertical surface displacements are ${\eta }_{+}={\xi }_z(H)$ and ${\eta }_{-}={\xi }_z(-H)$, the surface density perturbation of a displaced uniform interface givesThus Poisson equation reduces to Laplace equation for ${\varphi }^{\prime }$ inside and outside. The potential is continuous, while integrating through each surface givesThe perturbing potential must decay at vertical infinity. The kinematic conditions are ${\eta }_{\pm }={\partial }_z\chi (\pm H)$. The free surface pressure condition is the vanishing Lagrangian pressure perturbation, $p^{\prime }+{\eta }_{\pm }{p}_0^{\prime }=0$; at the top this gives $p^{\prime }(H)/{\rho }_0=g_s{\eta }_{+}$.

Put $q=kH>0$. For even scalar-potential symmetry, chooseThen ${\eta }_{-}=-{\eta }_{+}$: this is a thickness-changing mode, not a vertical bending mode. Write $\eta ={\eta }_{+}$; the top derivative jump gives $-kC(\cosh q+\sinh q)=4\pi G{\rho }_0\eta$. Since $\cosh q+\sinh q=e^q$,Substitute these values into ${\omega }^2\chi (H)=g_s\eta +{\varphi }^{\prime }(H)$. Using $4\pi G{\rho }_0=g_s/H$ gives the even-mode dispersion relation for the surface modes of a self-gravitating incompressible slab:The positive term comes from the restoring hydrostatic free-surface pressure in the background gravitational field. The negative term is the attraction due to the perturbed gravitational potential of the displaced surfaces, which tends to reinforce the disturbance. Both effects ultimately involve the slab's gravity, but enter the boundary condition through different terms.

For odd scalar-potential symmetry, instead takeNow ${\eta }_{-}={\eta }_{+}$, so the surfaces bend in the same direction. The jump again gives $C=-4\pi G{\rho }_0\eta /(k{e}^q)$, but $\chi (H)=\eta \tanh q/k$ and ${\varphi }^{\prime }(H)=-4\pi G{\rho }_0\eta /[k(1+\coth q)]$. The same pressure condition yieldsThe bottom kinematic, jump and pressure conditions follow with the indicated parity and give the same dispersion relation, so no boundary condition has been discarded.

To determine the gravitational instability of an incompressible slab, factor the two dimensionless expressions asFor the even mode the bracket is strictly increasing, with derivative $1+e^{-2q}>0$, is negative at zero and tends to positive infinity. It therefore has a unique positive zeroThe even mode is unstable for $0<kH<q_c$, neutral at $q_c$, and oscillatory for $kH>q_c$. Negative ${\omega }_e^2$ gives a growing mode with $\omega =i\sqrt{-{\omega }_e^2}$ under the stated time convention. For the odd mode, $1-e^{-2q}<2q$ for every $q>0$, so the bracket is always positive: all odd modes with positive wavenumber are stable.

At long wavelengths $F_e=-q+2q^2+O(q^3)$ and $F_o=q-2q^2/3+O(q^3)$: a thickness or column-density disturbance has a gravitational instability, whereas a bending disturbance is restored. At short wavelengths both approach $q-1/2$, so both surface branches are stable. The slab equilibrium is therefore linearly unstable overall, because it admits the long-wavelength even modes despite the stability of its odd sector.

The PDF specifies decay as $y\to \infty$; the local TeX drops this limit from its displayed boundary data. Take $a=4\pi$, so the conformal map is $w=\sin (z/4)$. Its components areFor $-2\pi <x<2\pi$ and $y>0$, $\eta >0$. The bottom edge maps monotonically to $(-1,1)$, and the two vertical sides map to the real rays $(-\infty ,-1)$ and $(1,\infty )$. The derivative $\cos (z/4)/4$ never vanishes in the interior. Injectivity follows from the identities for equal sines: the alternatives $z_1/4=z_2/4+2\pi n$ and $z_1/4=\pi -z_2/4+2\pi n$ give only the identical point within this strip. Surjectivity can be seen directly: for fixed $\xi$ and $v=y/4>0$, solve $\sin (x/4)=\xi /\cosh v$; thenOn $v>\max (0,\mathrm{arcosh}|\xi |)$, with the lower endpoint read as $0$ when $|\xi |\le 1$, the right side increases strictly from zero to infinity. Thus every point of the upper half-plane has a unique preimage.

By conformal invariance of harmonicity, write the transformed harmonic function as $\Phi (\xi ,\eta )$. Its real-axis boundary data areThe inverse sine on $(-1,1)$ is its real principal branch. For bounded integrable data, take the Fourier transform convention $\hat{\Phi }(k,\eta )={\int }_{\mathbb{R}}\Phi (\xi ,\eta )e^{-ik\xi }d\xi$. The Laplace equation becomes ${\hat{\Phi }}_{\eta \eta }-k^2\hat{\Phi }=0$, with the bounded decaying solution $\hat{\Phi }=\hat{F}(k)e^{-|k|\eta }$. Inverting the Fourier transform givesConsequently the Poisson kernel for the upper half-plane yieldsThe decay is uniform across the strip: $|w{|}^2={\sin }^2(x/4)+{\sinh }^2(y/4)\ge {\sinh }^2(y/4)$ and, for $|w|>1$, $|\Phi |\le \Vert F{\Vert }_{L^1}|w|/(\pi (|w|-1{)}^2)\to 0$. This verifies the printed decay condition. Boundary limits give the prescribed values at continuity points away from the corners.

For the given sine data,Values assigned at $s=\pm 1$ do not alter the integral. There is no solution continuous on the entire closed half-strip for these data: the bottom edge has corner limits $\pm 1$, but the side edges have corner limits zero. The formula is the bounded Dirichlet problem solution with the specified limits away from those two corners. This is the usual interpretation of the printed problem; boundedness excludes additional boundary-singular harmonic functions, while decay alone without a regularity or boundedness condition does not establish uniqueness.

The general-order Cauchy-Kovalevskaya theorem applies to analytic data on a non-characteristic hypersurface. In local coordinates $(t,z)$, a real analytic system of order $m$ in an unknown vector $U$ is in normal form whenwith $F$ real analytic near the initial jet. Analytic data ${\partial }_t^{j}U(0,z)=h_j(z)$ for $0\le j<m$ determine a unique real analytic solution near each initial point. More general systems, including different orders for the components, have the same conclusion after solving for the highest normal derivatives: the corresponding coefficient matrix, or highest-derivative Jacobian for a nonlinear system, must be invertible. This is the system's non-characteristic hypersurface condition. The uniqueness asserted by the Cauchy-Kovalevskaya theorem is initially uniqueness in the real analytic class.

A real analytic coordinate change flattens a real analytic initial hypersurface. Non-characteristicity permits solving for its highest normal derivatives by the real analytic implicit function theorem. Introduce all derivatives through order $m-1$ as additional unknowns. Their normal derivatives are either another jet variable, a tangential first derivative of a jet variable, or the right-hand side of the original equation. This gives a first-order real analytic system. The derivative-compatibility identities have zero initial data and are preserved by the system; equivalently, its real analytic coefficient recurrence reproduces the derivatives of $U$. This explains both reductions, rather than assuming arbitrary real analytic hypersurfaces are already flat.

For the first example, the principal symbol of the Laplace equation at the conormal $d{x}_1$ is $1$, so the plane is a non-characteristic hypersurface. The theorem therefore applies to real analytic $u{|}_{x_1=0}$ and ${\partial }_{1}u{|}_{x_1=0}$:This does not give Hadamard well-posedness in smooth or Sobolev space norms. For $n\ge 2$, the harmonic functionshave zero value data and normal derivative data $e^{-\sqrt{N}}\cos (N{x}_2)$ tending to zero with every fixed tangential derivative. At any fixed $x_1=\delta >0$, however, $u_N(\delta ,0)$ diverges. Thus the elliptic Cauchy problem for a partial differential equation is unstable despite real analytic solvability. The exceptional case $n=1$ is the elementary affine ordinary differential equation and has no such tangential high-frequency instability.

Fix the Riemann curvature tensor conventionFor the limiting affine connection, every product of two nonzero connection coefficients is zero, because its only nonzero coefficients have spatial upper index and both lower indices $t$. Its only possibly nonzero curvature components, apart from antisymmetry in the last two indices, areThe Ricci tensor is thereforeAll other components vanish, so the limiting connection is a connection with vanishing Ricci tensor exactly when $\varphi$ satisfies the Laplace equation on spatial Euclidean space. This is Ricci-flatness, not necessarily vanishing curvature: $\varphi =x^2-y^2$ is a harmonic function but has nonzero Hessian matrix and hence a nonflat limiting connection. Reversing the curvature sign convention reverses the displayed Ricci sign but leaves the equivalence unchanged.

Let the free surface be at $r=R_1$ and take the Newtonian gravitational potential to vanish at infinity. Mass conservation fixes $R_1=(3M_1/(4\pi {\rho }_1){)}^{1/3}$. The enclosed mass is $m(r)=4\pi {\rho }_1{r}^3/3$, and hydrostatic equilibrium requires $d{p}_0/dr=-{\rho }_{1}Gm(r)/r^2$. With vacuum outside, $p_0(R_1)=0$, givingThe pressure is zero outside. A prescribed constant external pressure would simply add that constant to $p_0$.

The bulk Eulerian and Lagrangian fluid perturbations obey $\delta \rho =0$ because the mass density is uniform and the Lagrangian displacement has zero divergence. The linearized Euler momentum equation in the interior isTaking the curl gives ${\omega }^2\nabla \times \mathbf{\xi }=0$. Thus every nonzero-frequency normal mode has an irrotational flow displacement; the simply connected space formed by the interior admits $\mathbf{\xi }=\nabla U$ and incompressibility gives the Laplace equation ${\nabla }^{2}U=0$. For a regular center, one spherical harmonic component is $U=A_l{r}^l{Y}_l^m$, with $|m|\le l$.

It is essential to retain the surface gravity perturbation of a uniform-density sphere. Although the bulk mass density perturbation is zero, the displaced density discontinuity givesThe Poisson equation implies that $\delta \Phi$ is continuous and its outward radial derivative has the jumpUsing regularity at the center and decay at infinity, writeThe derivative jump is $-(2l+1)C_l{R}_1^{l-1}{Y}_l^m$, so $C_l=-4\pi G{\rho }_{1}l{A}_l/(2l+1)$.

Integrating the bulk Euler momentum equation gives $\delta p/{\rho }_1={\omega }^{2}U-\delta \Phi$ for $l\ge 1$. At the free surface the Lagrangian pressure perturbation vanishes:Combining these relations yields the Kelvin stellar mode frequencyThe $2l+1$ choices of $m$ are degenerate because of spherical symmetry. For $l=1$, $r{Y}_1^m$ is a linear combination of coordinates in a Cartesian coordinate system, so $\nabla U$ is a constant displacement: the zero frequency is rigid translation of the whole isolated body. Such a translation has no restoring force. The constant $l=0$ potential generates no displacement; a radial breathing motion is excluded by incompressibility and a regular center.

The printed irrotational assertion needs the nonzero-frequency qualification, or restriction to the potential incompressible stellar surface modes. It is false for all neutral displacements: $\mathbf{\xi }=\mathbf{\Omega }\times \mathbf{r}$ has zero divergence, zero normal displacement at the surface, and $\delta p=\delta \Phi =0$, so it is a zero-frequency displacement, but $\nabla \times \mathbf{\xi }=2\mathbf{\Omega }\ne 0$. This does not change the Kelvin stellar mode spectrum just derived.

Interpret the malformed printed domain inequality as the first quadrant $x>0,y>0$, as indicated by the two axis boundary conditions. There is also a genuine missing hypothesis: decay of the boundary data does not require the harmonic function itself to be bounded. For zero data, both $q=0$ and $q=xy$ satisfy the Laplace equation and both edge values, but the Wirtinger derivative of the latter isThus the data do not determine an arbitrary solution's derivative without a growth restriction. We derive the integral representation for the bounded solution, which is also the standard decaying solution when the continuous edge data decay. Any other admissible unbounded solution can differ by a zero-boundary harmonic function.

The conformal map $w=z^2$ sends the quadrant bijectively to the complex upper half-plane: the positive real axis maps to $s>0$, and the positive imaginary axis maps to $s<0$. Define its real-axis boundary function byThe compatibility at the corner makes $f$ continuous there. The Dirichlet Poisson integral in a quadrant follows from the upper-half-plane Poisson kernel:It is bounded by the supremum of the bounded edge data, and the kernel's approximate identity property supplies both Dirichlet boundary conditions. More explicitly, for $z^2=a+ic$, where $a=x^2-y^2$ and $c=2xy>0$, this integral isIt displays directly which part of the real-line boundary corresponds to each original axis.

To obtain its derivative, first take real data and use the Schwarz integral formula. A holomorphic function whose real part is the Poisson solution in the $w$-plane isThe regularizing subtraction is independent of $w$ and affects only the imaginary normalization; it permits bounded data without an artificial integrability assumption. Then $q=\mathrm{Re}F(z^2)$. By the chain rule and the Wirtinger derivative,The factor is $z$, not $2z$: taking the real part contributes the factor $1/2$ in the complex derivative. Differentiation under the integral sign is justified uniformly on compact interior subsets because the denominator stays away from the real boundary and decays quadratically at infinity. Splitting the real integral at zero and substituting $s=r^2$ or $s=-r^2$ gives the holomorphic derivative of a quadrant Dirichlet solution:This converges for bounded data, and in particular for the stated sufficiently decaying data. No interior pole is encountered since $\mathrm{Im}(z^2)=2xy>0$. For complex-valued data the same formula holds by applying the real-data construction separately to their real and imaginary parts and using linearity.

An equivalent form, often convenient when the data's derivatives decay, comes from integration by parts. Sincethe two endpoint terms are $-g_2(0)/z^2$ and $+g_1(0)/z^2$. They cancel by the specified compatibility, leavingThis is the harmonic quadrant solution from tangential boundary derivatives applied to the actual derivatives $g_1^{\prime },g_2^{\prime }$; the supplied $g_1,g_2$ in this question are boundary values, not derivatives. Omitting those primes would solve a different boundary problem.

For completeness, bounded Dirichlet uniqueness in a quadrant selects the solution just constructed. The difference of two bounded solutions is zero on both axes. Odd reflection across each axis gives a bounded harmonic extension to the punctured plane; the removable singularity for a bounded harmonic function fills the corner. The harmonic Liouville theorem makes the extension constant, and oddness makes it zero. Decaying data give a decaying Poisson extension: split the real-axis integral into a fixed compact interval and a tail with arbitrarily small data; the compact contribution tends to zero as $|z^2|\to \infty$, while the tail is bounded by that small supremum.

As literally printed without boundedness or a growth condition, the most one can assert isThe zero-boundary harmonic growth in a quadrant includes $v=cxy$ and more generally $v=\mathrm{Im}(z^{2m})$. This is an actual nonuniqueness of the printed problem, not a defect in the integral for the normalized bounded solution.

Separate the uniform stream and an arbitrary time-dependent gauge from each velocity potential: ${\varphi }_j=U_{j}x+f_j(t)+{\varphi }_j^{\prime }$. Only gradients of the velocity potential determine the velocity field. The remote fluid must remain in its unperturbed uniform stream, soAfter fixing the harmless gauge, the nonzero-wavenumber perturbation velocity potentials themselves vanish at the corresponding infinities. Incompressible flow and irrotational flow give the Laplace equation for each perturbation. For a normal mode with $k>0$, the vertical amplitude solves ${\Phi }_j^{\prime \prime }-k^2{\Phi }_j=0$, henceThe opposite exponentials are excluded because they make the perturbation grow at infinity. It is the perturbation, rather than the total potential $U_{j}x+f_j(t)$, that decays.

Take $z$ upwards and write $\zeta =z+Wt$, so that the unperturbed interface is $\zeta =0$, with carbon dioxide above and water below. Use the sharp-interface, incompressible flow version of Darcy's law, with constant porosity and permeability of a porous medium, no capillary pressure, and effectively unbounded layers. The distinction between Darcy flux and interface speed matters: the base flux is $q_z=-\varphi W$, since the printed $W$ is the speed of the interface. Thus the base pressure gradients satisfyThe mass density ordering is stabilizing, whereas injection of the less viscous fluid into the more viscous fluid produces the Saffman–Taylor instability.

For linear stability analysis, displace the interface by $\eta =\hat{\eta }{e}^{\sigma t+i\alpha x}$. Taking the divergence of Darcy's law shows that the perturbation pressure satisfies the Laplace equation. Decay away from the interface gives, for $a_{\alpha }=|\alpha |>0$,The kinematic boundary condition and continuity of normal Darcy flux giveMeanwhile pressure continuity, evaluated on the displaced interface, givesEliminating $A,B$ yields the buoyancy-modified Darcy fingering dispersion relation:Every nonzero wavenumber decays when $W<W_c$, grows when $W>W_c$, and is neutral at equality, whereThe zero wavenumber is a uniform displacement and is neutral. Without a short-scale regularization such as surface tension, this model has no fastest-growing finite wavenumber: on the unstable side its growth rate increases without bound as $|\alpha |$ increases. Finite layer depths would change the dispersion relation; none are specified here.

Taking $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ and a year of $365.25$ days, the data giveThe numerical paragraph supplies no porosity, so it cannot determine a unique numerical interface speed. If $W$ were instead intended to mean Darcy velocity, the numerical critical value would be $103.2\mathrm{m}\mathrm{y}{\mathrm{r}}^{-1}$, and the interface would move at $W/\varphi$. These two velocity conventions must not be conflated.

Use the linearized boundary condition: evaluate at $z=0$ and discard products of perturbations. The kinematic boundary condition and Young–Laplace equation becomeIncompressible flow and irrotational flow give the Laplace equation for each potential. For a normal mode, the wall conditions selectWith $\eta =B{e}^{i(kx-\omega t)}$, the kinematic conditions give $A_1=i\omega B/[k\sinh (k{L}_1)]$ and $A_2=-i\omega B/[k\sinh (k{L}_2)]$. Hence ${\hat{\varphi }}_1(0)=i\omega B\coth (k{L}_1)/k$ and ${\hat{\varphi }}_2(0)=-i\omega B\coth (k{L}_2)/k$. Substitution into the dynamic condition gives the finite-depth Rayleigh-Taylor dispersion relation:Its denominator is positive. Assuming $g>0$ and $\gamma >0$, negative ${\omega }^2$ therefore produces a growing branch $\omega =i\sigma$ with $\sigma >0$, as in Rayleigh-Taylor instability.

Take $k>0$; reflection in $x$ supplies the result for negative wavenumber. Since each region has uniform base velocity, the perturbation is an incompressible, irrotational flow with a velocity potential satisfying the Laplace equation. Its decaying outer forms areFor a varicose mode of a planar jet, reflection symmetry makes the inner potential proportional to $\cosh (ky)$; for a sinuous mode of a planar jet, it is proportional to $\sinh (ky)$.

At the upper interface, the linearized kinematic boundary condition isThe first equation gives ${\varphi }_{+}(h)=-\sigma \eta /k$. The second givesThe linearized Unsteady Bernoulli equation gives $p^{\prime }=-\rho ({\partial }_t+U{\partial }_x)\varphi$. Applying pressure continuity therefore produces the equal-density top-hat planar-jet dispersion relationSolving the quadratic equation givesEquivalently,The plus sign is the unstable Kelvin-Helmholtz instability branch. In particular, the two modes have exactly the same growth rate,Its limiting forms areThus neither type grows faster in this equal-density idealization, although their phase speeds differ. The unbounded large-$k$ growth is the familiar short-wave pathology of an inviscid vortex sheet.

For an axisymmetric separated solution of the Laplace equation whose angular factor is the Legendre polynomial $P_n(\cos \theta )$, the radial ordinary differential equation has the two powers $r^n$ and $r^{-n-1}$. Thusfor constants $A,B$.

Write $x=\cos \theta$. Since $3\cos (2\theta )=6x^2-3=4P_2(x)-P_0(x)$, only the spherical harmonics of degrees zero and two occur. The degree-zero radial coefficient satisfying its two Dirichlet boundary conditions is $1-2/r$. For degree two, write $C{r}^2+D{r}^{-3}$; thenso $C=-4/31$ and $D=128/31$. By linearity and uniqueness of this boundary value problem,Direct substitution at $r=1,2$ verifies both boundary values.

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.

The far-field composition $-{\varphi }_B$ has chemical potential zero. Local curved-interface coexistence imposes $\mu (R)={\mu }_s=\gamma /({\varphi }_{B}R)$. The decaying spherically symmetric solution of the exterior Laplace equation isWith constant mobility, its outward current isThe linearized exterior composition is $\varphi (r)=-{\varphi }_B+\mu (r)/\alpha$, consistent with the surface shift in the previous part.

For diffusion-controlled droplet growth, the interior is uniform to leading order and has no diffusive current. Replacing exterior material by interior material requires the composition jump $2{\varphi }_B$, so conservation gives $2{\varphi }_B\dot{R}=-J_r(R)$. HenceThe negative sign represents evaporation. The total outward current is $4\pi R^2{J}_r(R)=4\pi M\gamma /{\varphi }_B$, independent of $R$ in this leading quasi-static capillary approximation.

The flat interface is an equilibrium reservoir with $\mu =0$, while the droplet has the constant surface value ${\mu }_s=\gamma /({\varphi }_{B}R)$. The exterior chemical potential satisfies the Laplace equation, tends to zero far away, and obeys these two surface values. These are exactly the electrostatic boundary conditions at a conductor for a sphere held at potential ${\mu }_s$ above a grounded infinite plane.

Introduce an arbitrary permittivity $ϵ$. In the electrostatic problem, $Q=-ϵ{\int }_{S_R}{\partial }_n\mu dS=C(R,h){\mu }_s$. Therefore the diffusion-capacitance analogy gives the outward currentFor an isolated sphere, its capacitance is $C(R,\infty )=4\pi ϵR$, so $I(R,\infty )=4\pi MR{\mu }_s=4\pi M\gamma /{\varphi }_B$. Dividing givesThus sphere-plane capacitance determines the evaporation enhancement. The assumed spherical shape gives the corresponding law $\dot{R}=-M\gamma c(R/h)/(2{\varphi }_B^2{R}^2)$ for droplet evaporation near a planar reservoir. The local surface current is generally nonuniform, even though its total is described by a single capacitance.

Expand $g$ in its Fourier sine seriesBy separation of variables, the decaying solution of Laplace equation in each half-strip has $n$th mode proportional to $e^{-k_n|x|}\sin (k_{n}y)$. Continuity at $x=0$ and integration of the equation across the Dirac delta function require the normal-derivative jumpThe jump of $C_n{e}^{-k_n|x|}$ is $-2k_n{C}_n$, so