To derive the magnitude of this net unbalanced force, let us first consider an infinitesimal patch of the interface, approximately rectangular in shape with sides $\Delta x\times\Delta y$ (see Figure 2.1a). Here we take two orthogonal directions in the tangent plane of the interface to be the $x$ and the $y$ axes of a Cartesian coordinate system, with its origin at the centre of the rectangular patch. The interface normal is then along the $z$ axis. The sum of the forces on the opposite sides, say those aligned along the $y$ direction, can be evaluated as follows. Using the sketch in Figure 2.2, the component of the interfacial force in the direction normal to the interface is

$\displaystyle\sigma\Delta y\times 2\sin\theta\approx\sigma\Delta y\times 2\theta\approx\sigma\Delta y\times\kappa_{xx}\Delta x,$

(2.5)

where $\kappa_{xx}$ is the curvature of the interface patch along the $x$ direction, such that $\theta=\kappa_{xx}\Delta x/2$. Here, the force along the tangent place to the interface vanishes. (If surface tension is not constant, then the tangential component evaluates to the Marangoni stress derived in §2.1.) Similarly, the unbalanced force on the sides aligned with the $x$ direction is

$\displaystyle\sigma\Delta x\times\kappa_{yy}\Delta y,$

(2.6)

where $\kappa_{yy}$ is the interface curvature along the $y$ direction. The net force may be written as

$\displaystyle\int_{C}\sigma\hat{\bm{n}}~{}\textrm{d}l~{}\approx~{}\sigma\Delta x\Delta y\times(\kappa_{xx}+\kappa_{yy}).$

(2.7)

Noting that $\kappa_{xx}+\kappa_{yy}$ is twice the mean curvature $\kappa_{\text{mean}}$, and in the limit of an infinitesimal patch size, we arrive at

$\displaystyle\int_{C}\sigma\hat{\bm{n}}~{}\textrm{d}l~{}\approx~{}\sigma\Delta x\Delta y\times 2\kappa_{\text{mean}}=\sigma\Delta A\times 2\kappa_{\text{mean}},$

(2.8)

where $\Delta A$ is the area of the patch. The result is independent of the choice of $x$ and $y$ axes direction, so long as they are along the tangent plane to the interface, because the mean curvature is an invariant of the interface geometry.

To see that the result is independent of the shape of the patch, and that the Marangoni component of the force is additive, let us consider the net force on an infinitesimal patch of the interface of arbitrary shape, such as shown in Figure 2.1(b). The tangent plane is again used to construct a local Cartesian coordinate system with its origin inside the patch, with the $x$ and $y$ coordinates parameterizing the patch. The curvature at the origin is, say, given by the tensor

$\displaystyle\bm{\kappa}=\begin{bmatrix}\kappa_{xx}&\kappa_{xy}\\ \kappa_{xy}&\kappa_{yy}\end{bmatrix}.$

(2.9)

The normal to the interface at any location parameterized by the coordinates $(x,y)$ is given to $O(x,y)$ by

$\displaystyle\hat{\bm{N}}=\bm{\hat{z}}+(\kappa_{xx}x+\kappa_{xy}y)\bm{\hat{x}}+(\kappa_{xy}x+\kappa_{yy}y)\bm{\hat{y}}.$

(2.10)

Noting that $\hat{\bm{n}}=\hat{\bm{t}}\times\hat{\bm{N}}$ and that $\hat{\bm{t}}=a\bm{\hat{x}}+b\bm{\hat{y}}+O(x^{2},xy,y^{2})$, i.e. it lies in the $xy-$plane, yields

$\displaystyle\hat{\bm{n}}=\hat{\bm{n}}^{\prime}+\bm{\hat{z}}\left\{\hat{\bm{n}}^{\prime}\cdot\left[(\kappa_{xx}x+\kappa_{xy}y)\bm{\hat{x}}+(\kappa_{xy}x+\kappa_{yy}y)\bm{\hat{y}}\right]\right\}+O(x^{2},xy,y^{2}),$

(2.11)

where $\hat{\bm{n}}^{\prime}=\hat{\bm{t}}\times\bm{\hat{z}}=a\bm{\hat{y}}-b\bm{\hat{x}}$ is the unit vector in the $xy$-plane that is perpendicular to $\hat{\bm{t}}$. The net force evaluates to

$\displaystyle\begin{split}\int_{C}\sigma\hat{\bm{n}}~{}\textrm{d}l&=\int_{C}\sigma\hat{\bm{n}}^{\prime}~{}\textrm{d}l+\bm{\hat{z}}\int_{C}\sigma~{}\hat{\bm{n}}\cdot\left[(\kappa_{xx}x+\kappa_{xy}y)\bm{\hat{x}}+(\kappa_{xy}x+\kappa_{yy}y)\bm{\hat{y}}\right]\textrm{d}l+O(x^{3},x^{2}y,xy^{2},y^{3})\\ &=\int_{A}\bm{\nabla}\sigma~{}\textrm{d}A+\bm{\hat{z}}\int_{A}\sigma\left[{\dfrac{\partial{(\kappa_{xx}x+\kappa_{xy}y)}}{\partial{x}}}+{\dfrac{\partial{(\kappa_{xy}x+\kappa_{yy}y)}}{\partial{y}}}\right]\textrm{d}A+O(x^{3},x^{2}y,xy^{2},y^{3})\\ &=\int_{A}\bm{\nabla}\sigma~{}\textrm{d}A+\bm{\hat{z}}\int_{A}\sigma(\kappa_{xx}+\kappa_{yy})\textrm{d}A+O(x^{3},x^{2}y,xy^{2},y^{3}),\end{split}$

(2.12)

where the order of the error term arises from the product of $O(x^{2},xy,y^{2})$ error in the integrand and the $O(x,y)$ extent of the integration domain along the length of the curve. In the limit of an infinitesimal interface patch, the integrals scale as the area of the patch , i.e. $O(x^{2},xy,y^{2})$, while the error term is one order smaller, and thus vanishes. This leads to

$\displaystyle\int_{C}\sigma\hat{\bm{n}}~{}\textrm{d}l=\underbrace{\int_{A}\bm{\nabla}\sigma~{}\textrm{d}A}_{\text{Marangoni stress}}+\bm{\hat{z}}\underbrace{\int_{A}\sigma(\kappa_{xx}+\kappa_{yy})\textrm{d}A}_{\text{Laplace pressure}}.$

(2.13)

As in §2.1, for flat interfaces, (2.13) reduces to the Marangoni stress and for constant $\sigma$, it derives its contribution only from Laplace pressure.

The results so far can be used to determine the shape of a liquid blob in equilibrium. For this determination, we will ignore the influence of gravity. Equilibrium dictates that the surface tension be uniform, i.e. agencies such as temperature and surfactants that cause surface tension to vary are uniform. In the absence of gravity, hydrostatics dictate that the pressure inside $p_{\text{in}}$ and outside $p_{\text{out}}$ the drop must also be uniform. The force balance on a small patch of the interface then only consitutes the forces of pressure and those acting on the interface as

$\displaystyle\textrm{d}\bm{F}=\int_{A}\left(p_{\text{in}}-p_{\text{out}}-2\sigma\kappa_{\text{mean}}\right)~{}\textrm{d}A=0.$

(2.14)

Consequently,

$\displaystyle p_{\text{in}}=p_{\text{out}}+2\sigma\kappa_{\text{mean}},$

(2.15)

or, in other words, the inside of the blob is pressurized relative to the outside by an amount equal to the Laplace pressure. Since the inside pressure must be uniform, the shape of the blob must be one with constant mean curvature, i.e. a sphere.

Energetic considerations may alternatively be used to arrive at the same conclusion. The thermodynamic energy of the system is derived from the surface energy of the blob. Minimizing this energy is then equivalent to determining the shape the minimizes the surface area, while maintaining a fixed volume equal to the volume of the blob. This shape is a sphere.