%~Mouliné par MaN_auto v.0.27.3 2023-07-24 10:21:32
\documentclass[CRMECA, Unicode, XML, screen]{cedram}
%
%\TopicEN{Section information missing}
%\TopicFR{Information section manquante}


\usepackage{amssymb}
\usepackage{graphicx}
\usepackage{subcaption}
\usepackage[labelformat=simple]{subcaption}
\renewcommand\thesubfigure{(\alph{subfigure})}

\newcommand{\dis}{\displaystyle}

\newcommand{\lb}[1]{\label{#1}}
\newcommand{\bq}{\begin{equation}}
\newcommand{\eq}{\end{equation}}

\newcommand*\R{\mathbb R}
\newcommand{\clE}{\mathcal{E}}
\newcommand{\bend}{\mathit{bending\,stress}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\graphicspath{{./figures/}}

\newcommand*{\mk}{\mkern -1mu}
\newcommand*{\Mk}{\mkern -2mu}
\newcommand*{\mK}{\mkern 1mu}
\newcommand*{\MK}{\mkern 2mu}

%\hypersetup{urlcolor=purple, linkcolor=blue, citecolor=red}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\RubricEF{Physical Sciences in Microgravity within the Thematic Group Fundamental and Applied Microgravity}{Sciences physiques en microgravité au sein du GDR Micropesanteur Fondamentales et appliquée}


\title{A direct relation between bending energy and contact angles for capillary bridges}

\alttitle{Une relation directe entre l'énergie de flexion et les angles de contact pour les ponts capillaires}

\keywords{Distortion of nonaxisymmetric capillary bridges, Mean and Gaussian curvatures impact, Euler characteristic, Generalized Young--Laplace equation, Bending effects, Fenchel's theorem in differential geometry, Gauss--Bonnet Theorem, Geodesic curvature, Bending stress, Influence of the contact angles}

\altkeywords{Distorsion des ponts capillaires non axisymétriques, impact des courbures moyennes et gaussiennes, caractéristique d'Euler, équation de Young--Laplace généralisée, effets de flexion, théorème de Fenchel en géométrie différentielle, théorème de Gauss--Bonnet, courbure géodésique, contrainte de flexion, influence des angles de contact}
\subjclass{49N45, 53A10, 58E12, 74F10, 74G05, 74G15, 53Z05}


\author{\firstname{Olivier} \lastname{Millet}}
\address{LaSIE, UMR-CNRS 7356, Université de La Rochelle, avenue Michel Crépeau, 17042 La Rochelle cedex 1, France.}
\email{olivier.millet@univ-lr.fr}

\author{\firstname{Gérard} \lastname{Gagneux }\IsCorresp}
\address[1]{LaSIE, UMR-CNRS 7356, Université de La Rochelle, avenue Michel Crépeau, 17042 La Rochelle cedex 1, France.}

\email{gerard.maryse.gagneux@gmail.com}

\begin{abstract} 
The didactic object of these developments on differential geometry of curves and surfaces is to present fine and convenient mathematical strategies, adapted to the study of capillary bridges. The common thread is to be able to calculate accurately in any situation the bending stress over the free surface, represented mathematically by the integral of the Gaussian curvature over the surface (called the total curvature) involved in the generalized Young--Laplace equation. We prove in particular that the resultant of the bending energy is directly linked to the wetting angles at the contact line.
\end{abstract}

\begin{altabstract}
L'objet didactique de ces développements basés sur la géométrie différentielle des courbes et des surfaces est de présenter des stratégies mathématiques adaptées à l'étude des ponts capillaires. Le fil conducteur est de pouvoir calculer avec précision, dans n'importe quelle situation, la contrainte de flexion de la surface libre d’un pont capillaire, représentée mathématiquement par l'intégrale de courbure de Gauss (courbure totale) de la surface libre intervenant dans l'équation de Young--Laplace généralisée. Nous établissons en particulier un résultat très général suivant lequel la résultante de l'énergie de flexion est directement liée aux angles de mouillage au niveau de la ligne de contact.
\end{altabstract}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\dateposted{2023-09-29}
\begin{document}
\maketitle
\section{Introduction}


Under boundaryless manifold condition, the most common formulation of the~Gauss--Bonnet integration theorem indicates that for a spherical drop or a soap bubble without contact, with or without bump, the integral of the Gaussian curvature\footnote{The product of the main curvatures. } over the surface, proportional to the bending energy, is invariant if one bends and deforms the surface (this value is~a~topological invariant).

Indeed, according to the Gauss--Bonnet integration theorem, for a closed free surface, the energy contribution of the Gaussian curvature during deformation is directly proportional to the Euler constant (see Eq.~\eqref{intK}), and therefore is constant as long as the topology of the surface, described by the~Euler characteristic of the surface, does not change~\cite{w}, and can be ignored when determining the shape of such a membrane.


This has probably favored the in-depth study of constant mean curvature~surfaces, excluding gravity effects~\cite{RFINN,GMMY,DEG,GIWI,Orr,SAUV,SAUVI}. In the case of capillary bridges, the presence of contact surfaces does not allow this simplification (the total geodesic curvature of the boundary is to be taken into account~to apply the Gauss--Bonnet integration theorem) and makes it a priori necessary to take into account the Gaussian curvature, to establish a hierarchy of the various configurations with regard to the bending effects and to introduce the generalized Young--Laplace equation. Some other works proposed a general law for continuum media with interface or a derivation of thermo balance equation for systems with interface~\cite{Isola87}. A second order thermodynamical approach has been used to model surface tension of bubbles leading to a generalization of Young--Laplace theory~\cite{Isola95}, but without involving the Gaussian curvature of the interface.


In this work, we present various results and complementary strategies of mathematical ana\-lysis that can be applied to concrete capillary bridges problems, concerning in a new way, the Gauss--Bonnet and Fenchel's theorems to establish various analytical formulas easy to use for capillary bridges. In a first step, we will focus on surfaces of revolution (circular boundaries, which makes these cases much easier, especially for the explicit calculation of the total geodesic curvatures of the boundaries, rarely possible in practice by the integral calculus). We prove that, in the general case including particular axisymmetric capillary bridges, the resultant of the bending energy is directly linked to the wetting angles at the contact line. We also highlight the determining parameters and their respective influence in the bending energy and its variation.



Then the approach is extended to the rather delicate modelling of nonaxisymmetric capillary bridges distortions. The key to achieving generalization is a direct consequence of the Fenchel's theorem in differential geometry which avoids a lot of dead-end integration calculations. These developments relate to surfaces of revolution on the basis of an unit speed reparameterization (or by arc length) for a regular curve, in this case, the semi-meridian. For detailed presentations of the subject, the reader may refer to~\cite[p.~161-164]{carmo}, \cite[p.~161-162]{MONT} and also to~\cite{car,GRAY}.



\section{Generalized Young--Laplace equation and associated}

The generalized Young--Laplace equation concerns the strong distortions for which the bending effects are modeled by an additional curvature-related term involving the Gaussian curvature $K$ though a multiplier coefficient $C_{K}$ which has the dimension of a force and stands for the bending stress~\cite{BORU,GAYD,SCHO}. According to generalized Young--Laplace equation, the downward vertical measurement $x$ in relation with the value $\Delta p_{0}$ at $x=0$ (a spontaneous unknown value), may be linked to the mean an Gaussian curvature according to~\cite{DEL,GAYD,mazz,SAUV,SAUVI}:
\begin{equation}\label{YLgen}
\gamma\left(\frac{1}{\rho_{c}}+\frac{1}{N}\right) +C_{K}\frac{1}{\rho _{c}\ N}=\Delta p_{0}-\Delta\rho\ gx\text{,}
\end{equation}
where $C_{K}$ has a dimension of a force\footnote{ $C_{K}$ results from the physics at the interface molecular scale. Its sign a priori depends on the dynamics of the wetting (advancing or recessing wetting angle) and on the value of the wetting angle (see Eq. \eqref{Ebend2}).}. Therefore the term $C_{K}/ \rho_{c}N$ stands \emph{via }a pressure for the local bending stress, where $\rho_{c}$ and $N$ denote the principal radii of curvature (evaluated algebraically, positively when the curvature is turned into the interior of the capillary bridge). Finally, $\Delta p_{0}$ is the pressure deficiency at $x=0$ and $\Delta\rho$ the difference of the densities between the fluid and the gas.
%\medskip

It is assumed that the different cefficients, implicit unknown a priori, as $\Delta p_{0}$, resulting from the final equilibrium, have been previously identified \emph{\ in situ} from exprerimental data, by solving a linear system, well posed and numerically stable (for example, thanks to a first integral and a principle of conservation) \cite{GAG,GG,MMGY}, \cite{Ng,Ng2019a,Ng2019b,Ng2019c,Ng2020a,Ng2021}.
%\medskip

It is extremely noteworthily~\cite{rodval}, that this strongly nonlinear differential equation is mathematically isomorphic (the same structure) but with different variables and physical units, to the {Gullstrand equation} of geometrical optics, which relates the optic power $P_{op}^{\prime}$ of a thick lens\linebreak (in diopters, the reciprocal of the equivalent focal length) to its geometry and the properties of the media. For example, the superficial tension $\gamma$ is equivalent to the refractivity $\frac{n_{1}}{n_{2}}-1$, where $n_{i}$ is a refractive index, $C_{K}$ is analogous to the expression $-(\frac{n_{1}}{n_{2}}-1)^{2}\frac{n_{2}}{n_{1}}d$, $d$ the lens thickness and $\Delta p_0$ corresponds to $P_{op}^{\prime}$.
%\medskip

Shear or free energy problems and the longitudinal bending stress of ship hulls have an analogy with the subject~\cite{carr,millet1997,millet2003,OK,QUI}. The mathematical modeling and simulations of the petroleum engineering are also concerned by this theoretical topic, in order to obtain for media with periodic microstructure an ``equivalent'' macroscopic representation, by some statistical or homogenization methods~\cite[Chapter~1]{GMT}.


The bending stress over the free surface $\Sigma$ may be represented in the following integral form, at the dimension of a force:
\[
\clE_{\bend}=C_{K} \int_{\Sigma} K\ d\Sigma,
\]
where $K$ is the Gaussian curvature of the free surface $\Sigma$, intrinsic value, in particular independent of the choice of the unit normal vector, and the nondimensional integral is the total curvature\footnote{For example, the total curvature of the catenoid whose axis is of infinite length is $-4\pi$, the total curvature of the sphere of radius $r$ is $4\pi$ and the torus $0$.}.


Concerning the capillary tension forces, by term by term integrating over the free surface $\Sigma$ the generalized Young--Laplace equation, we have for example the relationship between various forces:
\[
\gamma {\displaystyle\int_{\Sigma}}
\dis{ \left(\frac{1}{\rho_{c}}+\frac{1}{N}\right) } 
 d\Sigma=-C_{K} {\displaystyle\int_{\Sigma}} K\ d\Sigma \ + {\displaystyle\int_{\Sigma}} \left(\Delta p_{0}-\Delta\rho\ gx\right)d\Sigma,
\]
with the particular situation:
\[
\gamma {\displaystyle\int_{\Sigma}}
\dis{ \left(\frac{1}{\rho_{c}}+\frac{1}{N}\right) } d\Sigma=-C_{K} {\displaystyle\int_{\Sigma}} K\ d\Sigma+ \Delta p_{0} \text{ area }(\Sigma)
\]
when neglecting gravity effects.
%\medskip

This would allow to have a reasoned opinion on the relative importance of the bending forces; according to an objective criterion, either by relative value or by intrinsic value.


\section{Homotopic surfaces, Euler characteristic and Gauss--Bonnet theorem}\label{secGB}

\subsection{General theory}

Recall that the Euler characteristic (or Euler--Poincaré characteristic) is a topological invariant, an integer that describes, according to precise axiomatic principles, the shape or a structure of a topological space {regardless of how it is bent} according to the formula: number of vertices$-$number of edges$+$ number of faces with the property of invariance by homeomorphy. It is commonly denoted by $\chi$ or $\chi(M) $. As examples for surfaces in homological algebra, we have $\chi(M)=2$ for a sphere, $\chi(M) =4 $ for two spheres (not connected), $\chi(M) = 0 $ for a torus and $\chi(M) = -2$ for a two-holed torus.
%\medskip

To speak very figuratively, quite approximately, the Euler--Poincaré characteristic is {an integer}, invariant when the size and the shape of a geometrical object change by an effect of a ``plastic'' deformation.
%\medskip

This invariance property makes it a providential tool in the context of this study on the bending effects, associated to the Gauss--Bonnet theorem, a deep relationship between surfaces in differential geometry, connecting the Gaussian curvature of a surface to its Euler characteristic.
%\medskip

The Euler characteristic of the right cylinder is zero, thus so is that of the cylinder with one or two boundaries. These following free surfaces with two circular boundaries and whose meridian is an arc of Delaunay roulette are considered topologically equivalent (same common topological genus), because it is possible to continuously move one to obtain the other: portion of concave or convex, catenoid or unduloid (the right cylinder being the transition case). Accordingly, these axisymmetric surfaces have in common the same Euler characteristic, in this case, the value zero. It is the same for their continuous axisymmetric smooth deformations by distorting effect of bending or gravity~\cite{mazz,MYSH}.
 

The Gauss--Bonnet theorem is reputed to be one of the most profound and elegant results of the study of surfaces~\cite{car,carmo,CIA,w}. It has no surprisingly many applications in Physics. It is used in sectors of activity where the problems of bending beams surely arises (civil engineering, naval architecture, shell theory to predict the stress and the displacement arising in an elastic shell, \cite{CIA, EqComp, Hamdouni2003b, Hamdouni2006}, etc\dots).
%\smallskip
%\smallskip

In fact, it unexpectedly links two completely different ways of studying a surface: one geometric, the other topological. Indeed, for any compact, boundaryless two-dimensional Riemannian manifold $\Sigma$, the integral of the Gaussian curvature $K$ over the entire manifold with respect to area measure is 2$\pi$ times the Euler characteristic of $\Sigma$, also called the Euler number of the manifold,~i.e.
\begin{equation}\label{intK}
\int_{\Sigma} K\ d\Sigma=2\pi\chi\left(\Sigma\right).
\end{equation}

For example, for a sphere $\Sigma$ of radius $R$ in $\R^3$, it comes:
\begin{equation}\label{int-R}
\int_{\Sigma} K\ d\Sigma=\dfrac{1}{R^{2}} 4\pi R^{2}=4\pi \text{ and here} \ \chi\left(\Sigma\right) =2
\end{equation}

Suppose now that $M$ is a compact two-dimensional Riemannian manifold with a boundary $\delta M$ and let $k_{g}$ the signed geodesic curvature of $\delta M.$ Then, in nondimensional writing,
\begin{equation}\label{gaussB1}
\int_{M} K\ dM+ \int_{\delta M}
\ k_{g} ds= 2\pi\chi\left(M\right).
\end{equation}



%\medskip


We recall that the {geodesic curvature} $\ k_{g}$, of an arbitrary curve at a point $P$ on a smooth surface, is defined as the curvature at $P$ of the orthogonal projection of the curve onto the plane tangent to the surface at $P$ and we have:
\begin{equation}\label{defKg}
k_g=k \cos \theta_g,
\end{equation}
where $\theta_g$ is the angle between the osculating plane of $C$ and the tangent plane $Q$ at point $P$, which corresponds to the local contact angle for capillary bridges wetting a plane surface. In particular, when the curve $C$ representing $\delta M$ is a circle of radius $R$, we have $k=\frac{1}{R}$ and therefore:
\begin{equation}\label{intKg}
\int_{\delta M}\ k_{g} ds= 2\pi \cos \theta_g
\end{equation}


\begin{figure}[!h]
\centering
\includegraphics[width=0.65\linewidth]{Figures/Geodesic.png}
\caption{Representation of the geodesic curvature. }\label{geodesic}
\end{figure}

%\medskip



\subsection{Application to droplets and capillary bridges} 
As explained in the introduction, the case of a spherical droplet (closed surface) of radius $R$ has no interest, as the contribution of the Gaussian curvature is directly proportional to the Euler constant $\chi(M) $ whose value $\chi(M) =2$ for a sphere\footnote{ $\chi(M) =1$ for a half-sphere or a cap.}. In that case, relation~\eqref{gaussB1} gives no supplementary information, as no wetting area exists.

%\medskip



Let us consider now the case of a spherical droplet lying on a plane surface as represented on Fig.~\ref{fig:goutte1}. The application of Gauss--Bonnet theorem leads to calculate both contributions in the right hand side of equation~\eqref{gaussB1}. The first one is directly link the external surface area of the spherical cap representing the droplet and is equal to $ 2 \pi (1- \cos \alpha)$, where $\alpha$ is the opening angle (Fig.~\ref{fig:goutte1}). The contribution of the geodesic curvature is equal to $2\pi\cos \theta $, according to~\eqref{intKg}, since the contact line is circular and the droplet lies on a plane surface~\eqref{fig:goutte1}. In that case, \eqref{gaussB1} leads to the obvious and known geometric relation $ \theta = \alpha$.
For a droplet lying on the sphere, as the wetted surface is not plane, the geodesic curvature is $k_g= \frac{1}{r_c} \cos \theta_g $ with $\theta_g \neq \theta$. In that case, the Gauss--Bonnet theorem leads to the relation $\theta_g = \frac{\pi}{2} - \delta$.

\begin{figure}[!h]
\centering
\subcaptionbox{Representation of a droplet on a plane surface.\label{fig:goutte1}}
{\includegraphics[scale=0.7]{Figures/petite_goutte_rc.pdf}}
\hfill
\subcaptionbox{Representation of a droplet on a sphere.\label{fig:goutte2}}{\includegraphics[scale=0.7]{Figures/calotte_sur_sphere.pdf}}
\caption{}
\end{figure}


For axisymmetric capillary bridges whose contact lines are circles, typically capillary bridges between two parallel planes (Fig.~\ref{fig:plane}), we obtain a general expression of the bending stress for a surface of revolution\footnote{We have denoted $\theta_{1}$ and $\theta_{2}$ the upper and lower contact angles assumed to be constant. In the case of a plane wetted surface, $\theta_{g}$ corresponds to the wetting angle $\theta$.}
\begin{equation}\label{Ebend2}
\clE_{bending\ stress}=C_{k} \int_{M} K\ dM=-2\pi C_{k}\left(\cos\theta_{1}+\cos\theta_{2}\right),
\end{equation}
with here the Euler characteristic $\chi(M) =0 $ and when the total geodesic curvature at the boundaries is $2\pi(\cos\theta_{1}+\cos\theta_{2}) $ according to~\eqref{intKg}.

\begin{figure}[!h]
\centering
\includegraphics[width=0.65\linewidth]{Figures/planes.png}
\caption{Capillary bridge between two parallel planes. }\label{fig:plane}
\end{figure}

For symmetric profiles with $\theta_1=\theta_2=\theta$, we obtain a relationship between the contact angle $\theta $ and the bending stress:
\[
\theta =\arccos \frac{\displaystyle \int_{M} - K\ dM}{4\pi }.
\]

Coming back to~\eqref{Ebend2} with the same contact angle $\theta$, the relative finite variation of the bending stress, as function of the contact angle $\theta $, is then given by the formula:
\begin{equation}\label{delataEb}
\frac{\delta \left(\clE_{\bend}\right)}{\clE _{\bend}}=-\tan \theta \ \delta \theta \ , \ \theta \neq \frac{\pi}{2}.
\end{equation}

%\medskip
In summary, it should be kept in mind that the value of the bending stress depends, besides physical constants, only on the observed values of the contact angles, whereas these angles result in part implicitly from the final equilibrium of the device.



%\medskip

In the general case of capillary bridges, using the classification and the associated parameterization of~\cite{GAG}, the Gauss--Bonnet theorem leads to a supplementary relation linking the geometric properties of the capillary bridge. In the particular case of a catenoid\footnote{Transition between unduloid or nodoid shape.}, where the parameterization of the meridian is given by
\[
y(x)= y^* \cosh \left(\frac{x}{y^*}\right),
\]
$y^*$ denoting the gorge radius (Fig.~\ref{fig:plane}), we may obtain an explicit and useful relation between between $\theta$, $y^*$ and the length $D$ of the capillary bridge.

%\medskip



In the next, we want to generalize these expressions {in any situation}, for capillary bridges with non necessarily circular contact lines in non-axisymmetric cases (this is the case for instance when coalescence of capillary bridges occurs~\cite{GG}).


%\medskip



\subsection{The prevaling roles of the contact angles values and of the convexity/nonconvexity of the boundaries.}


For didactic purposes, we mainly limited ourselves to the detailed case of axisymmetric capillary bridges to show the interest of the simultaneous use of the Gauss--Bonnet formula and of the topological notion of Euler characteristic to evaluate the importance of the bending stress. The obtained result clearly shows in an explicit way, the major role of the {contact angles values} after distortion effects, eventually distinct (Eqs.~\eqref{Ebend2} and~\eqref{delataEb}). Therefore, all the factors determining the contact angle have consecutively an influence on the bending stress (surface roughness and heterogeneity, influence of gravity, contact angle hysteresis~\cite{carr,DEG,H,VANH,wang}).
 
 
It is well known that the contact angle value is determined by the balance between adhesive and cohesive forces on the rigid supports. As the tendency of a drop to spread out over a flat, solid surface increases, the contact angle decreases. Thus, the contact angle provides an inverse measure of wettability. In this context, the case of the right cylinders is still a borderline case.
%\medskip

A contact angle less than $\frac{\pi}{2}$ (low contact angle) usually indicates that wetting of the surface is very favorable, and the fluid will spread over a large area of the surface. Contact angles greater than $\frac{\pi}{2}$ (high contact angle) generally mean that wetting of the surface is unfavorable. It should be quoted that a certain number of terms of the generalized Young--Laplace equation are spontaneous values, resulting from instantaneous equilibrium, and are therefore implicit unknowns. This is a difficulty for the mathematical resolution of this nonlinear differential boundary problem.

%\medskip
In addition, the Fenchel’s theorem sheds light on the importance of the convexity or the nonconvexity of the outer edges in calculating exactly the value of the total curvatures.

Consequently, even a limited displacement of the surface boundary can modify the bending stress, by local modifications of the contact angles or affecting the local curvature of the outline curve and then, the total geodesic curvature. The contact angle hysteresis can also be significant.


\section{The aim and effective convenience of the arc length reparameterization strategy}

The generalization of expression~\eqref{Ebend2} for non necessarily circular contact lines and possibly non closed will be performed using Fenchel's theorem based on an unit speed reparameterization (or by arc length) of the contact line to calculate the total geodesic curvature.



\subsection{Surface of revolution}

The rather paradoxical aspect of this method is that it is only very rarely easy to get in practice an explicit calculation formula. However, it leads to general quantitative results in the form of analytical formulas very easy to use from the experimental data, via a very convenient expression of the Gaussian curvature of a surface of revolution (the speed and acceleration vectors are then orthogonal).


To illustrate what we are talking about, let us consider a smooth curve of the half-plan $\{y>$ $0,$ $z=0\}$ parametrized by the arc length. The surface of revolution resulting in ${\R}^{3}$ from the rotation of the curve around the $x$-axis, $\psi$ being the angle of rotation, is parametrized by:

\[
M\left(\left(s,\psi\right) \right) =\left(x\left(s\right),\ y\left(s\right) \cos\psi,\ y\left(s\right) \sin\left(
\psi\right)\right) \ 0\leq\psi\leq2\pi, 0\leq s\leq L.
\]
As the meridian portion is parameterized by arc length, we have\emph{ ipso facto} the following remarkable and convenient relations and convenient expressions for angular (in radians) and trigonometric values as well as for the Gaussian curvature of the surfaces of revolution:
\begin{equation}\label{xprime}
x^{\prime2}+y^{\prime2}=1 \ \text{ at any point},
\end{equation}
and therefore, by differentiating, the orthogonality relationship
\[
x^{\prime}x^{\prime\prime}+y^{\prime}y^{^{\prime\prime}}=0,
\]
that is to say that $T.\frac{dT}{ds}=0$ where the dot denotes the scalar product of $\R^3$ and where $T(s)=\frac{dM(s)}{ds}$ is the unit tangent vector to the curve $M(s)=M(s,0)$.



The Gaussian curvature $K$ \ of the surface of revolution has then the very convenient expression (see~\cite[p.~162, eq.~(9)]{car}):
\[
K\left(s,\theta\right) =-\dfrac{y^{^{\prime\prime}}\left(s\right)}{y\left(s\right)}.
\]

In this context, the expression of the mean curvature $H$ of a surface of revolution is less attractive (see~\cite[p.~162, eq.~(11)]{car}).


A remarkable illustrative example of the arc length reparameterization strategy is the determination of the axisymmetric surfaces of constant Gaussian curvature\footnote{Problem studied more extensively by Gaston Darboux 1890. Among the solutions, surfaces are found that look like a hyperboloid.}. We consider then the classical differential equation, linear, of the second order, homogeneous:
\[
y^{^{\prime\prime}}\left(s\right) +K y\left(s\right) -0, \; 0 \leq s\leq L,
\] 
with the three following cases: $K<0,$ $K=0,$ $K>0$. Then, introducing the general form of the corresponding solutions in $y$, we consider the resulting differential equations
\begin{equation}
x^{\prime2}=1-y^{\prime2}, 0 \leq s\leq L
\end{equation}
resulting from~\eqref{xprime}.



Moreover, $\varphi$ being the angle that the tangent to the profile curve makes with the $x$-axis, we have the following relationships:
\[
\mathit{sin}\,\varphi=\dfrac{y^{\prime}\left(s\right) }{\sqrt{x^{\prime2}\left(s\right) +y^{\prime2}\left(s\right) }},\text{i.e.} \ \ \sin
\varphi=\ y^{\prime}\left(s\right) \ \ \text{and} \ \cos \varphi=x^{\prime }\left(s\right).
\]

To compute the global bending stress in this context, we have to consider succesively:
\[
\clE_{\bend}=C_{k} \int_{M} K\ dM=\dis{ C_{k}\int\int-\frac{y^{^{\prime\prime}}\left(s\right) }{y\left(s\right) }\ y\left(s\right) d\psi ds }
\]
and therefore
\[
\dis{ C_{k}\int\int-y^{^{\prime\prime}}\left(s\right) d\psi ds=-2\pi
C_{k}\ \left(y^{\prime}\left(L\right) -y^{\prime}\left(0\right)
\right) }
\]
so that
\[
\clE_{\bend} =
\dis{ - 2\pi(C_{k}\left(\sin\left(\varphi\left(L\right) \right) -\sin\left(
\varphi(0\right) \right),}
\]
$\varphi$ being the angle that the tangent to the profile curve makes with the $x-$axis, the axis of rotation.
%\medskip


Let us quote that this case of a surface of revolution around $x$-axis may correspond to a capillary bridge between two parallel planes at $x=0$ and $x=L$. With the notation of Fig.~\ref{fig:plane}, we have $ \varphi(L) = \pi /2 -\theta_1 $ and $ \varphi(0) = \theta_2 - \pi /2 $ and we recover the general expression~\eqref{Ebend2} for axisymmetric capillary bridges.


%\medskip

\subsection{The conclusive Fenchel’s theorem for the general case}

In the case of non circular contact lines, for instance portions of an ellipse, the parameterization involves elliptic integrals, rarely possible to explain in practice, so that the parameters would have to be sought numerically (spline interpolation) \cite{HA,MESMOU}. This computational difficulty is overcame by knowing the {Fenchel’s theorem}~\cite{car,HOR}, which shows the complementarity of the three methods leading to a generalization of expression~\eqref{Ebend2} for possibly non convex or even non closed contact lines.

%\bigskip

To well illustrate the interest of Fenchel’s theorem associated with the theorems of {Gauss and Bonnet}, let us consider a reparameterization by arc length of the curve $\Gamma$. 


\subsubsection{Close plane curve}
According to Fenchel's theorem\footnote{The Fary--Milnor theorem concerning the total curvature of the knotted closed curves does not seem appropriate for the subject of this study.} (1929), the value of the total curvature
\[
\int_{\Gamma} k(s)ds 
\]
of any smooth closed space curve ${\Gamma} $ is at least $\ 2\pi$\emph{, i.e} $ \int_{\Gamma} k(s)ds\geq2\pi$. The equality holds if and only if the curve is a {convex plane} curve. In other words, the average curvature of a closed convex plane curve equals $2\pi/L$, where $L$ is the length (the perimeter) of the curve\footnote{For a given arc of a plane curve, the { local} average curvature quantifies the ratio of the change in inclination of the tangent to the curve over the arc length.}.

By the Fenchel’s theorem, without calculation of primitive functions, often tedious or ineffective, we deduce directly, for any closed convex plane curve $\Gamma$ (i.e. the curve is the boundary of a convex set in the Euclidean plane), that
\[
\dis{ \int_{\Gamma} k(s)ds=2\pi.}
\]
This case is certainly the most encountered in practice when the boundaries of the capillary bridge (the contact lines) are two closed plane convex curves, not necessarily circular. We then recover expression~\eqref{Ebend2} which is still valid in this more general case\footnote{ When the wet solid surface is not plane, relation~\eqref{Ebending2} involves $\theta_{g1}$ and $\theta_{g2} $ which corresponds to $\theta_{1}$ and $\theta_{2}$ to within a constant linked to the geometry of the wet solid interface.}:
\begin{equation}\label{Ebending2}
\clE_{\bend}=C_{k} \int_{M} K\ dM=-2\pi C_{k\ }\left(\cos\theta_{1}+\cos\theta_{2}\right).
\end{equation}


In the case of a closed nonconvex plane curve, we are led to conclude by defining of the notion of the {winding index}, a topological argument, in what follows.



\subsubsection{Open plane curve}

Let us give some classical preliminary elements of differential geometry related to smooth boundaries of surfaces, parametrized by arc length. The curvature of a plane curve parametrized by arc length is the rate of turning of the tangent line with respect to an \emph{ad hoc} frame along the curve.



Let $\varphi (s)$ be the angle of inclination of the unit tangent vector $T=T(s)$ with respect to a fixed frame of reference, for instance $x-$axis. Considered then as a rate of turning for the tangent line when one moves along the curve at unit speed, the curvature $k(s)$ becomes
\[
k(s) =\dfrac{d\varphi}{ds}(s)= \varphi^{\prime}(s).
\]

It follows that the total curvature of a smooth curve $C$ is then given by the formula depending only of the initial and final states:
\begin{equation}\label{ks}
 \int_{C} k(s)ds=\varphi(\textit{ending}) -\varphi(\textit{starting}),  (\text{in radians}).
\end{equation}

For a piecewise smooth curve parametrized by arc length, then we need to deal with the exterior angles at the corners according to the orientation of the curve in the turning motion. However, up to now, to our knowledge such capillary bridges with non convex or open contact line have are not considered in literature.

\section{The general case and its implementation}

In the general cases of nonaxisymmetric capillary bridges between two supports, possibly of distinct natures, the method remains applicable in principle. The difficulty is not conceptual in dealing with the general case but rather calculative. We must then, in any given case, engage in a delicate exercise in differential and analytical geometries to explicitly calculate the total signed geodesic curvature of the boundaries by the classical methods of analytical geometry.
%\medskip

\subsection{The calculation procedure is as follows}

At any point $P$ of the border liquid-solid, one considers the tangent plane in $P$ to the free surface (that supposes an adequate local regularity). One then considers the orthogonal projection of each edge into this tangent plane. The curvature in $P$ of the projected curve is then calculated, what introduces the important role of the cosine of the local contact angle and leads to expression~\eqref{defKg} of the the curvature $k$ which related related to the geodesic curvature $k_g$ at $P$ by the relationship:
\[
k_g=k\cos \theta_g
\]
where $\theta_g$ is linked to the local contact angle $\theta$ to within a constant depending on the geometry of the wet solid surface. When it is plane, $\theta_g$ corresponds to the wetting angle $\theta$, that will be considered in the next to simplify the developments.

%\medskip

When the contact angle is constant on the considered contact surface, we have the particularly simple relationship:
\[
\int_{\Gamma_{i}} \ k_{g}(s) ds=cos\theta \int_{\Gamma_{i}} k(s)\ ds.
\]
%
%\medskip
\subsection{The special situation of heterogeneous contact surfaces}

When the contact angles are separately variable on each of the contact surfaces, i.e. $\theta=\theta_{1}(M)$ and $\theta=\theta_{2}(M)$ according to the physical conditions of the two surfaces (non ideal smooth surfaces), the integral along each boundary, corresponding to the total geodesic curvature, in fact, of the kind
\begin{equation}\label{nonconstT}
\int_{\Gamma_{i}} cos\theta_{i}(M(s))\ k(s)\ ds,
\end{equation}
is more complicated to calculate with computational prediction of wetting (at our knowledge, an open problem for the probably most realistic case). The use of a mean theorem would likely be imprecise (effects of surface roughness). 
%\medskip


\subsection{The general case of homogeneous contact surfaces}

When multiplied by the coefficient $(-C_{K})$ at the dimension of a force, the dimensionless integral of these curvature values along the reunion of the two contact edges gives finally the value of the resulting bending stress by the Fenchel’s theorem (the cornerstone of the method).

The three possible scenarios then arise according to the geometry of the boundaries (closed plane convex or nonconvex curves) are the following, the surfaces having in common, without loss of generality, the same Euler characteristic, in this illustrative case, the value zero.
%\medskip

By introducing the contact angles $\theta_{1}$ and $\theta_{2\text{ }}$(in radians) on each outline of contact surfaces, we proved that, at least theoretically, the wettability being evaluated, here, by constant contact angles, separately on each contact support.
%\bigskip

\paragraph{Case 1:}
The boundaries are two closed plane convex curves. Then,
\[
\clE_{\bend}=C_{k} \int_{M} K\ dM=-2\pi C_{k\ }\left(\cos\theta_{1}+\cos\theta_{2}\right).
\]

\paragraph{Case 2:}
The boundaries are two closed plane curves, one convex and the other nonconvex. Then,
\[
\clE_{\bend}=C_{k} \int_{M} K\ dM=-C_{k\ }\left(2k_{1}\pi\cos\theta_{1}+2\pi\cos\theta_{2}\right),
\]
the observed integer $k_{1}$, $k_{1}\geq2$, being the winding number of the nonconvex curve (the winding index in algebraic topology).


\paragraph{Case 3:}
The boundaries are two closed disjoint plane curves, nonconvex. Then,
\[
\clE_{\bend}=C_{k} \int_{M} K\ dM=-C_{k\ }\left(2k_{1}\pi\cos\theta_{1}+2k_{2}\pi\cos\theta_{2}\right),
\]
$k_{1}$ and $k_{2}$ being the integers, $\geq2$, winding numbers of the curves, observed \;and known\emph{ in situ}.

In the rather theoretical case, where the value of the Euler characteristic is non-zero, it should be necessary to write:
\[
\clE_{\bend}=C_{k} \int_{M} K\ dM=2\pi\chi\left(M\right) -C_{k\ }\left(2k_{1}\pi\cos\theta_{1} +2k_{2}\pi\cos\theta_{2}\right).
\]

It must be emphasized that, when the contact angles are separately variable on each of the contact surfaces according to the physical conditions of the two surfaces (non ideal smooth surfaces), the integral along each boundary, $\ $corresponding to the total geodesic curvature of the plane and closed boundaries, seems a serious difficulty to explain. The question might interest specialists in differential geometry.


\section{Conclusion}

The developments obtained here for surfaces of revolution and their generalization to more general surfaces representing the shape of capillary bridges, result from concepts in differential geometry and geometric analysis with applications to Lagrangian Mechanics, without resorting to differential calculus and integral calculus. The methods of {Euler's} characteristic, associated to the Gauss--Bonnet--Binet theorem and the strongly complementary Fenchel’s theorem apply immediately to the cases of the nonaxisymmetric surfaces, with explicit, easy-to-use, results formulations.
%\medskip


We proved that in the general way, the value of the bending stress depends, besides physical constants, only on the observed values of the contact angles, whereas these angles result in part implicitly from the final equilibrium of the device. Therefore, all the factors determining the contact angle have an influence on the bending stress (surface roughness and heterogeneity, influence of gravity, contact angle hysteresis. 

It would be interesting to reconsider, in taking into account these new results concerning the bending effects, the important role of the contact curves geometry and the Gauss--Bonnet and Fenchel theorems, an analytical framework for reassessing the cohesion effects of coalescence between saddle shaped capillary bridges~\cite{GG}. Finally, by creating a support material having a nonconvex region with high wettability and a complementary region with very low wettability, the experimenter could illustrate the theory by experimentation.

\section*{Conflicts of interest}
The authors have no conflict of interest to declare.

\section*{Acknowledgments}
The authors would like to express their sincere gratitude to the French National Centre for Space Studies (CNES) for having supported this work. They also thank Antoine Moreau (LaSIE) and Marc Medale (IUSTI) for valuable discussions.


\bibliographystyle{crunsrt}
\bibliography{crmeca20220694}
\end{document}




