==== Front R Soc Open Sci R Soc Open Sci RSOS royopensci Royal Society Open Science 2054-5703 The Royal Society 10.1098/rsos.201534 rsos201534 1009 73 Physics and Biophysics Research Article Identities for droplets with circular footprint on tilted surfaces Identities for Drops on Tilted SurfacesDunlop François 1 http://orcid.org/0000-0001-6330-4085Fatollahi Amir H. 2 Hajirahimi Maryam 3 Huillet Thierry 1 1 Laboratoire de Physique Théorique et Modélisation, CY Cergy Paris Université, CNRS UMR 8089, 95302 Cergy-Pontoise, France 2 Department of Physics, Faculty of Physics and Chemistry, Alzahra University, Tehran 1993891167, Iran 3 Physics Group, South Tehran Branch, Islamic Azad University, P. O. Box 11365, Tehran 4435, Iran Author for correspondence: Amir H. Fatollahi e-mail: fath@alzahra.ac.ir 11 2020 11 11 2020 11 11 2020 7 11 20153426 8 2020 20 10 2020 © 2020 The Authors.2020http://creativecommons.org/licenses/by/4.0/http://creativecommons.org/licenses/by/4.0/Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.Exact mathematical identities are presented between the relevant parameters of droplets displaying circular contact boundary based on flat tilted surfaces. Two of the identities are derived from the force balance, and one from the torque balance. The tilt surfaces cover the full range of inclinations for sessile or pendant drops, including the intermediate case of droplets on a wall (vertical surface). The identities are put under test both by the available solutions of a linear response approximation at small Bond numbers as well as the ones obtained from numerical solutions, making use of the Surface Evolver software. The subtleties to obtain certain angle-averages appearing in identities by the numerical solutions are discussed in detail. It is argued how the identities are useful in two respects. First is to replace some unknown values in the Young–Laplace equation by their expressions obtained from the identities. Second is to use the identities to estimate the error for approximate analytical or numerical solutions without any reference to an exact solution. dropshydrostaticssurface tensionYoung–LaplaceAlzahra Universityhttp://dx.doi.org/10.13039/501100012310D99/3/1447cover-dateNovember, 2020 ==== Body 1. Introduction Understanding the skewed shape of a sessile drop pinned on a flat incline has a long history in Physics, starting with [1,2]. It rises new asymmetrical problems compared to more studied situations where the substrate is horizontal. In principle, such problems can be handled while making use of the Young–Laplace nonlinear partial differential equation, translating a balance between surface tension forces and gravity acting on the drop. See [3] and references therein, where a perturbative approach to this problem at small Bond number was addressed when the footprint of the droplet is held fixed and circular. More recently, the case of pendant drops has also attracted some interest, see [4–7] and references therein. A perturbative approach to this problem at small Bond number has also been addressed in [8], in a similar set-up. An empirical relation between incline slope angle and contact angles at the front and rear of the droplet was given by [9], and further studied by many authors, see [3,10,11] and references therein. It relies on an approximation of the balance of forces equation along the substrate, at small Bond number. Balance of forces normal to the substrate also deserves interest, together with relations arising from the torque balance. The three-phase contact angles (for the various azimuthal angles) of a liquid condensed on a substrate are in direct relation with interfacial and body forces acting on sessile or pendant drops. Despite the settled role of the droplets based on different surfaces, the cases for which there are exact solutions are rare. During more than a century, different numerical and analytical methods have been developed for more efficient and finer approach to drop’s profile for cases with no exact solution. To evaluate or rate these numerical and analytical methods some criteria are needed, among which are the exact mathematical relations between the relevant parameters of the problem. An early example of these identities is the one by [12] between the volume, curvature at apex, height and contact radius of axi-symmetric drops on a flat horizontal surface (see also [13,14]). For axi-symmetric drops on curved surfaces the very same identities are derived in [15]. For droplets under the combined tangential and normal body forces the dynamical relations between shape parameters have been presented in a linear approximation recently [16]. This issue has also been studied semi-analytically in [17], where the problem of understanding the contact line evolution of slender unpinned droplets under arbitrary scenarios of forces is addressed, based on experimentally observed contact lines. Related to this point, in [18], sessile droplets at different tilting angles are experimentally subject to varying centrifugal forces in order to explore their spreading/sliding behaviour for different volumes and initial shapes (including non-axisymmetric). In particular, a test of the applicability of the Furmidge equation for the retention force is discussed. The mathematical identities are important in another respect, that is reducing the initial unknown values of the problem. This in particular proves helpful because some of these unknown values appear in the first place in the differential equation governing the profile of the drop. Whether one tries a perturbative solution of the drop’s profile or a numerical one, reducing the initially unknown values facilitates or boosts the procedure of reaching the final result. As an example, in [14] the identity is used to replace a combination of the unknown apex curvature and height in the Young–Laplace equation, in the procedure of developing a perturbative solution for lightweight drops. As another example, in a numerical solution of axi-symmetric drops with fixed contact angle one may use contact radius or apex height as a starting point and the other one as the final point. But the trouble is that both values are unknown in the first place, and in principle one has to search time-consumingly a two-dimensional parameter space to find both values that match the solution. However, by the identity and thanks to the mentioned combination, one can replace the height by the radius, reducing the procedure to a simple one-parameter shooting method [19]. The main purpose of the present work is to highlight examples of mathematical identities for droplets with circular contact boundary based on flat tilted surfaces. The identities are derived based on the force balance along parallel and normal directions of the tilted surface, as well as the torque balance of the droplet. The force balance normal to the tilted surface generalizes that of [12–14] for a horizontal substrate. The identity along the tilted surface is in fact an exact version of the approximate empirical Furmidge relation [9]. The identity by the torque balance is apparently the one that is introduced here, and remains to be verified numerically. All three identities are checked at the first-order approximation of the Bond number, the so-called linear response ansatz [3]. For the identities by the force balance various numerical tests are provided, generated by the Surface Evolver as a vertex-edge-facet element software [20]. The organization of the rest of the work is as follows. In §2 the two identities by the force balance along and normal to the surface are derived. In §3 the identity by the torque-balance condition on the droplet is derived. The check of all identities at linear response approximation is presented in §4. Numerical checks of force-balance identities are presented in §5. This requires the contact angle as a function of azimuth angle and certain averages over it, which is the subject of §6. In §7 the possible use of identities is illustrated. Section 8 is devoted to concluding remarks. 2. Identities by force balance The set-up for the pinned droplet on a tilted plane is as follows. The z-axis is perpendicular to the substrate and inward to the liquid, with x-axis along the slope downward. For the droplet with contact-circle of radius r0, the origin is set to be the centre of circle, c. The angle of the substrate with the horizontal direction is α ∈ [0, π], with α = π representing the case of a drop pendant from the ceiling. The set-up is summarized in figure 1. We use the polar angle φ on the substrate, with φ = 0 representing the x-axis, as usual. The hydrostatic pressure inside the contact-circle depends only on x, given by 2.1 p(x)=pc+ρgxsin⁡α with ρ the density of the drop, and pc the pressure at the centre of the contact-circle. Figure 1. Set-up of a droplet on a flat tilted surface. 2.1. Force balance along surface First we consider the more familiar identity, stemming from force balance along the substrate. The capillary force along the x direction reads 2.2 Fx=γ∫−ππr0 dφ cos⁡θα(φ) cos⁡φ with θα(φ) the contact-angle of the liquid-substrate at polar angle φ on a slope with angle α. By symmetry θα(φ)=θα(−φ), so we expect the following Fourier expansion: 2.3 cos⁡θα(φ)=C0+∑n=1∞Cncos⁡(nφ), where 2.4 C0=12π∫−ππ dφcos⁡θα(φ)=⟨cos⁡θα(φ)⟩ and 2.5 Cn=1π∫−ππ dφcos⁡θα(φ)cos⁡(nφ),n≥1 leading to 2.6 Fx=γ πr0 C1. The force balance along the surface then gives: 2.7 mgsin⁡α+γ πr0 C1=0 with m = ρV as the mass of the droplet with volume V. The dimensionless form of (2.7) is 2.8 2 Bo sin⁡α+πC1=0, where Bo is the modified Bond (or Eötvös) number taken as Bo = mg/(2r0γ). Consistently with [21], our use of Bo contains only the input parameters of the drop-solid system, namely volume of the drop, gravitational acceleration, fluid density, surface tension, and radius of the footprint and not data involving the dimensions of the drop unknown prior to the experiment (such as the height of the drop or radius of curvature at the drop apex). The coefficient C1, equal to twice the average over φ of cos⁡θα(φ)cos⁡φ, is remarkably linear in Bo as will be seen later. In §5, in order to confront theoretical and simulation values, we will evaluate and compare to 1 the ratio 2.9 ratio∥=−2 Bo sin⁡α/(πC1). 2.2. Force balance normal to surface The capillary force normal to the surface (downward z direction) reads 2.10 Fz=γ∫−ππr0 dφ sin⁡θα(φ). Again by symmetry we expect θα(φ)=θα(−φ), hence the Fourier expansion 2.11 sin⁡θα(φ)=A0+∑n=1∞Ancos⁡(nφ), 2.12 A0=12π∫−ππ dφsin⁡θα(φ)=⟨sin⁡θα(φ)⟩ 2.13 andAn=1π∫−ππ dφsin⁡θα(φ)cos⁡(nφ),n≥1 leading to 2.14 Fz=2γ πr0A0. The pressure force from the tilted surface is given by 2.15 N=∫−r0r0p(x) dA,dA=2r02−x2 dx and 2.16 N=πpcr02. The zero of pressure is chosen as the atmospheric pressure. Balance of forces in the normal direction leads to 2.17 N=Fz+mgcos⁡α or 2.18 πpcr02=γ 2πr0A0+ρVgcos⁡α. The pressure pc can be written in terms of the pressure at point o (intersection point of the z-axis and the surface of the droplet), as follows: 2.19 pc=po+ρghcos⁡α with h the oc height (figure 1). This pressure, with the use of the Young–Laplace equation, can also be written in terms of the mean-curvature H0 at point o 2.20 po=−2γ H0 with γ as the surface tension of liquid. Altogether the force balance (2.18) reads 2.21 πr02(−2γ H0+ρghcos⁡α)=2γπr0A0+ρVgcos⁡α or, after dividing by 2πr02γ, 2.22 −H0 +ρgh2γcos⁡α=A0r0+ρVg2πr02γcos⁡α. The above is a direct generalization of eqn (12) of [14,15,19] for a horizontal substrate (α = 0). It is convenient to have, with the help of the zero-gravity copy (spherical cap) of the droplet, a dimensionless form of the above identity. Using the spherical-cap radius R0 and contact-angle θ0 with 2.23 r0=R0sin⁡θ0,V=πR033(1−cos⁡θ0)2(2+cos⁡θ0) one defines the dimensionless quantities [3] 2.24 h~=hR0,H0~=R0H0. The dimensionless version of (2.22) is 2.25 −H~0+3h~sin⁡θ0 Bocos⁡απ(1−cos⁡θ0)2(2+cos⁡θ0)=A0sin⁡θ0+Bo cos⁡απsin⁡θ0. In §5, in order to confront theoretical and simulation values, we will evaluate and compare to 1 the ratio 2.26 ratio⊥=−H~0 sin⁡θ0+3h~sin2⁡θ0π(1−cos⁡θ0)2(2+cos⁡θ0)Bo cos⁡αA0+1π Bo cos⁡α. 3. Identity by torque balance The contribution of the weight to the torque in the y-direction (inward figure 2) is given by 3.1 τwy=mgℓcos⁡β=ρVgℓcos⁡β in which ℓ is the distance of the centre-of-mass (c.o.m.) of the droplet from the centre c, and β is the angle between the c.o.m. position vector and the horizontal direction (figure 2). The torque applied by the substrate upon the drop can be evaluated by integration over the element by the pressure 3.2 dτ py=−x p(x) dA,dA=2r02−x2 dx leading to 3.3 τ py=−2∫−r0r0(pc+ρgxsin⁡α)x r02−x2 dx and 3.4 τ py=−π4ρgr04sin⁡α. The infinitesimal capillary force at polar angle φ with contact-angle θα(φ) reads 3.5 dFγ=γ r0 dφ(cos⁡θα(φ) cos⁡φ i+cos⁡θα(φ) sin⁡φ j−sin⁡θα(φ) k) with the position vector 3.6 r=r0(cos⁡φ i+sin⁡φ j) leading to the torque element in y-direction 3.7 dτγy=(r×dFγ)y=γr02sin⁡θα(φ) cos⁡φ dφ. Figure 2. The geometry used for torque-balance identity. The integration over the above torque element gives 3.8 τγy=γr02∫−ππsin⁡θα(φ) cos⁡φ dφ. Again using Fourier expansion (2.11) for sin⁡θα(φ) one finds 3.9 τγy=πγr02 A1. All together, the balance of torques along y-direction gives 3.10 ρVgℓcos⁡β+πγr02A1=π4ρgr04sin⁡α. Again it is convenient to have a dimensionless form of the identity. Defining 3.11 ℓ~=ℓR0 and using (2.23) and the relation between two Bond numbers [8] 3.12 Bo=πB(1−cos⁡θ0)2(2+cos⁡θ0)6sin⁡θ0 the dimensionless form of identity (3.10) reads 3.13 2 Boℓ~cos⁡β+πsin⁡θ0A1=π4B sin3⁡θ0sin⁡α. In §5, the following ratio is compared to 1 by the numerical simulations 3.14 ratio↺=2 Boℓ~cos⁡β+πsin⁡θ0A1π4B sin3⁡θ0sin⁡α. 4. Linear response 4.1. Check at linear response, along surface The linear response ansatz [3] is in terms of the Bond number 4.1 B=ρgR02γ, where R0 is the radius of the spherical cap at zero gravity, see (2.23). In the linear response approximation we have [3,8] 4.2 C1=−12(cos⁡θαmin−cos⁡θαmax)+O(B2) in which θαmax=θα(0) and θαmin=θα(π). We see that the expression (2.7) in linear approximation is the famous Furmidge relation (eqn 1 of [9]), with the constant K = π/4. The check of (2.8) by explicit expressions for cos⁡θαmin and cos⁡θαmax at the linear approximation [8] is straightforward, once the relation between two Bond numbers (3.12) being used. 4.2. Check at linear response, normal to surface The linear response ansatz [3,8] implies 4.3 cos⁡θα(φ)=cos⁡θ0+λB+μBcos⁡φ+O(B2), where θ0 is the uniform contact angle at B = 0, λ and μ are some constants whose values do not matter here. Then 4.4 sin⁡θα(φ)=1−cos2⁡θα(φ)=sin⁡θ0−λBcos⁡θ0/sin⁡θ0−μBcos⁡φ cos⁡θ0/sin⁡θ0+O(B2) which implies both the average over φ 4.5 A0=⟨sin⁡θα(φ)⟩=sin⁡θ0−λBcot⁡θ0+O(B2) and the arithmetic mean between maximum and minimum 4.6 12(sin⁡θαmax+sin⁡θαmin)=sin⁡θ0−λBcot⁡θ0+O(B2). Hence 4.7 A0=12(sin⁡θαmin+sin⁡θαmax)+O(B2). We have h~=1−cos⁡θ0+O(B) [3]. By the above, identity (2.25) comes to the form 4.8 −H0~=sin⁡θαmin+sin⁡θαmax2sin⁡θ0−B6(1−cos⁡θ0)(1+2cos⁡θ0)1+cos⁡θ0cos⁡α+O(B2). The H~0, sin⁡θαmin and sin⁡θαmax are read from [3,8], by which we have 4.9 −H~0=1−B6(1−cos⁡θ0)cos⁡α+O(B2) and 4.10 12(sin⁡θαmin+sin⁡θαmax)=sin⁡θ0−Bcos⁡θ0r01′(θ0)cos⁡α+O(B2) with [3] 4.11 r01′(θ0)=−sin⁡θ06+sin⁡θ0cos⁡θ03(1+cos⁡θ0). It is a simple matter to check that the quantities (4.7)–(4.10) inserted into (2.25) satisfy (2.25) up to O(B2). 4.3. Check at linear response: torque Proceeding as in §4.2, the linear approximation yields 4.12 A1=−12(sin⁡θαmin−sin⁡θαmax)+O(B2) 4.13 =13B(1−cos⁡θ0)(2+cos⁡θ0)1+cos⁡θ0cos⁡θ0sin⁡α+O(B2). For a check at linear order, it is enough to insert the spherical cap droplet values in the terms having B or Bo, namely the first and last terms in (3.13). The c.o.m. of a spherical cap is known. Subtracting R0cosθ0 leads to 4.14 ℓ~=3(1+cos⁡θ0)24(2+cos⁡θ0)−cos⁡θ0+O(B). Also for the spherical cap the c.o.m. lays on the z-axis, for which we have β = π/2 − α. By the relation between the two Bond numbers (3.12), it is easy to see that identity (3.13) is satisfied up to O(B2). 5. Numerical check of identities To test the identities (2.8) and (2.25), the solutions of the Young–Laplace equation are developed by the Surface Evolver software, at different slope angles αs, spherical cap contact angles θ0s and modified Bond numbers Bo = mg/(2r0γ). The identities take the form ratio∥=1,ratio⊥=1,ratio↺=1 with the ratios defined in (2.9), (2.26) and (3.14). In order to measure the mean curvature H~0 and height h~ above the origin, we export from Surface Evolver the list of vertices such that in cylindrical coordinates r~<0.06 r~0, that is 6 per cent of the contact radius. This yields between 80 and 110 vertices. We then fit a quadratic surface 5.1 z~=ax~2+by~2+cx~+h~ to obtain h~ and 5.2 −H~0=(a+b(1+c2))(1+c2)−3/2+O(B2). The measurement of A0 and C1 implies more work, described in §6. In figures 3–5 we display A0, H~0, h~, ratio⊥, −C1 and ratio ∥ as functions of Bo for α = 90° and θ0 = 60°, 90°, 120°. Tests were made with the following values: 5.3 α=45∘,90∘,135∘θ0=60∘,90∘,120∘} and many different values of the Bo number, away from the singularity, giving consistent values in all cases. As a sample of the numerical values, the data for the case with α = 90° and θ0 = 90° are given in table 1. Figure 3. Testing the identities for α = 90° and θ0 = 60°. Figure 4. Testing the identities for α = 90° and θ0 = 90°. Figure 5. Testing the identities for α = 90° and θ0 = 120°. Table 1. Testing the identities for α = 90° and θ0 = 90°. Bo ratio ∥ ratio⊥ 0.63971 1.0000 1.0006 0.95957 1.0010 1.0001 1.2794 1.0008 0.99987 1.4073 1.0023 0.99609 1.4713 1.0017 1.0020 In order to check the torque-balance identity numerically it is needed to find the c.o.m. of the droplet, for which out of the irregular distribution of vertices one has to extract regular dx × dy-mesh and dy × dz-mesh, for calculating z~cm and x~cm, respectively. The result of making meshes for one of the samples is given in figures 6 and 7. The values of l~ and angle β are then obtained by 5.4 β=tan−1zcmxcm−α and 5.5 ℓ~=x~cm2+z~cm2. To obtain the average-value A1 the azimuth-angle is needed; the computation is postponed to §6. The results of torque identity checks are summarized in table 2. Figure 6. dx × dy mesh for sample α = 30° and θ0 = 45°, with Bo = 1.015. Mesh size: 0.01 × 0.01. Figure 7. dy × dz mesh for the sample of figure 6. Mesh size: 0.01 × 0.005. Points x > xzmax in blue and x < xzmax in orange. Table 2. Testing the torque-balance identity. α θ0 Bo ratio 30° 45° 1.02 0.999 30° 90° 0.512 1.01 45° 30° 0.453 0.999 60° 60° 0.401 1.00 6. Contact angle as function of azimuth In order to check the Surface Evolver simulations against the exact identities, we have to compute C0 and C1, or the average over φ ∈ (0, 2π) of sinθ(φ) and the average over φ ∈ (0, 2π) of cosφcosθ(φ), where θ(φ) is the contact angle at azimuth φ. A fluid interface in contact with a solid surface has a contact angle θ ∈ (0, π) which may vary along the contact line. The solid surface must be smooth at the macroscopic scale, so that a unique normal vector is defined at every point. We have performed the Surface Evolver simulations with Dirichlet boundary conditions: the displacement is zero on the contact line, a fixed circle of radius r0. Physically, other than a droplet on or below an incline, it may represent a pocket of liquid made with an elastic membrane fixed onto a circular metallic wire, opening to a reservoir of liquid. Initially, at zero gravity, the pocket has the shape of a spherical cap of contact angle θ0, giving the desired volume. The pocket lies entirely on one side, say {z ≥ 0}, of the plane containing the contact line circle. Upon switching on the gravity, the pocket deforms as a solution of the Young–Laplace equation under constant volume constraint and the given Dirichlet boundary conditions. Eventually the interface or membrane may go partly into {z < 0} region. Of course, so far as the interface remains in {z ≥ 0}, it represents as well a liquid drop on or below a solid plane substrate, with contact line pinned on the circle and contact angle θ ∈ (0, π) varying along the circle. However, where the interface goes into {z < 0}, representing a pocket, a contact angle θ > π or θ < 0 will be found. This occurs e.g. for α = π/2, θ0 = π/3, Bo = 1.46 around φ = π (figure 8) and for α = π/2, θ0 = 2π/3, Bo = 1.32 around φ = 0 (figure 9). At such parameters, a droplet on the incline would have unpinned, dewetting from the top around φ = π or overflowing at the bottom around φ = 0. The Young–Laplace equation then has to be solved with moving boundary, which we have not done. Figure 8. Cosine of contact angle θ as function of azimuth φ for α = 90° and θ0 = 60° at Bond numbers Bo = 0.50, 1.00, 1.45. Linear response approximation is plotted at Bo = 0.50. Figure 9. Cosine of contact angle θ as function of azimuth φ for α = 90° and θ0 = 120° at Bond numbers Bo = 0.44, 0.88, 1.32. Linear response approximation is plotted at Bo = 0.44. The figures also show the linear response approximation, namely cosθ as a linear function of cosφ, at the smallest Bond number for which data are displayed. As another noticeable feature, it is observed that to some approximation, for each choice of α and θ0, there is an azimuth φ at which θ≃θ0 ∀Bo. The number of mesh vertices in the Surface Evolver simulation is up to 700 000, with finer mesh near the contact line and finer mesh also in the more delicate cases where the Bond number Bo approaches the instability threshold. In order to measure the contact angle θ as function of the azimuth φ, we export from Surface Evolver the list of vertices such that |z~|<0.03 (1−cos⁡θ0), that is 3 per cent of the height of the drop at zero Bond number. This yields between 1300 and 20 000 vertices. We divide the range of azimuth into a hundred equal sectors. Each sector contains between four and six hundred vertices. In each sector we fit with gnuplot a surface containing the contact line in the form, in cylindrical coordinates, r(z)=r0−zcot⁡θ+az2orz(r)=−(r−r0)tan⁡θ+b(r−r0)2. The fit is performed, independently in each sector, in terms of cot⁡θ and a or tanθ and b. 7. Uses of identities 7.1. In first place: an illustrative example The identity by force balance in normal direction can be used in the first place in the Young–Laplace equation. It helps to replace some unknown values in the equation from the beginning. This has been used in [14] to develop a perturbation solution and in [19] to reduce the numerical procedure to a simple shooting method, for sessile drops on horizontal surface with fixed contact-angle. Here we show that the identity can be used to replace the unknown quantities in the Young–Laplace equation for drops with fixed contact radius (pinned drops) on tilted surfaces as well. First one recognizes that the Young–Laplace equation can be written in the form 7.1 −2γH0+ρg(h−z)cos⁡α=−2γH in which H0 and h are unknown in the first place. Now by the identity (2.22) one has for the combination as follows: 7.2 −2γH0+ρghcos⁡α=2γA0r0+ρVgπr02cos⁡α bringing the Young–Laplace to the form 7.3 2γA0r0+ρVgπr02cos⁡α−ρgzcos⁡α=−2γH. Now in the r.h.s. only A0 is unknown. To develop a perturbative analytical solution, such as in [3], the above form has the advantage that it is involved by less to be calculated values. At linear order of B number, the above form finds its full advantage, because at the beginning even A0 is known by the linear ansatz (4.7). As mentioned earlier, the numerical solution in reduced form helps lowering the parameters space dimension [19]. 7.2. Estimation of error by identities As mentioned earlier, the exact mathematical relation may be used to make an estimation error for an analytical or numerical approach. As an example, we use identity (2.25) by force balance normal to the substrate to estimate the error of −H~0 by the linear approximation of [8], that is expected to be of second order in B. First let us rearrange the identity in the form that the apex curvature would be calculated by other values, namely 7.4 [−H~0]Identity=−3h~sin⁡θ0 Bocos⁡απ(1−cos⁡θ0)2(2+cos⁡θ0)+A0sin⁡θ0+Bo cos⁡απsin⁡θ0. The above gives the value of −H~0 required to hold the force-balance condition. It was seen earlier in §4.2 that, when only the first order is kept the above equality holds. Now 7.5 h~=1−cos⁡θ0+B6(1−cos⁡θ0+2ln1+cos⁡θ02)cos⁡α+O(B2). When expressed based on the number Bo via relation (3.12) (also (8) of [8]), the estimation of error is given by: 7.6 [−H~0]identity−[−H~0]linear solR~0−1 =−Bo23π2sin2⁡θ0(1−cos⁡θ0)4(2+cos⁡θ0)2 ×(1−cos⁡θ0+2ln1+cos⁡θ02)cos2⁡α in which we have used R~0=1. The remarkable observation by [8] is that, for relatively large values of the Bond number Bo the linear approximation solution is quite close to the numerical solution. In figure 10,we see that for Bond numbers up to Bo = 2.0 the error is less than 5%. Figure 10. The estimation of error of −H~0 versus Bo for α = 45° and θ0 = 120°. 8. Conclusion The main concern of the present work is to highlight exact mathematical identities for droplets, even in cases where an exact Young–Laplace profile solution is not available. As a special case, droplets with circular contact boundary based on flat tilted surfaces are considered. Two of the identities are derived as the requirements of the force balance, and one of the torque balance. The identities involve some relevant values of droplets. On the side of output values, some are available once an analytical or numerical solution would be given, for instance the curvature at apex or the vertical height. Some of these output values appearing in the identities are somehow indirectly available, such as certain azimuthal-angle averages, or the location of the centre-of-mass of the droplet. The identities are put under test both by the available solutions of a linear response approximation as well as the ones obtained from exact numerical solutions. We stress that the problem analysed here should not only be seen as an abstract theoretical construction. It is a real problem having an experimental counterpart, e.g. putting the droplet on a disc-shaped asperity on the substrate. Supplementary Material Reviewer comments Data accessibility This article has no additional data. Authors' contributions F.D. participated to define the problem, developed numerical analysis, drafted parts of text; A.H.F. suggested the study, contributed to maths derivations, drafted parts of text; M.H. contributed to maths derivations and fixing the results, critical reading the manuscript; T.H. participated in shaping the problem, contributed in checking and completing maths expressions, drafting parts of text and making final form of paper. All authors gave final approval for publication and agree to be held accountable for the work performed therein. Competing interests We declare we have no competing interests. Funding The work by A.H.F. is supported by the Research Council of Alzahra University. ==== Refs References 1 Macdougall G , Ockrent C 1942 Surface energy relations in liquid/solid systems I. The adhesion of liquids to solids and a new method of determining the surface tension of liquids . Proc. R. Soc. Lond. A 180 , 151 –173 . (10.1098/rspa.1942.0031 ) 2 Frenkel YI 1948 On the behavior of liquid drops on a solid surface 1. The sliding of drops on an inclined surface . Zh. Eksp. Teor. Fiz. 18 , 659 (Translated by V. Berejnov: https://arxiv.org/abs/physics/0503051) . 3 De Coninck J , Dunlop F , Huillet T 2017 Contact angles of a drop pinned on an incline . Phys. Rev. E 95 , 052805 (10.1103/PhysRevE.95.052805 )28618574 4 Milinazzo F , Shinbrot M 1988 A numerical study of a drop on a vertical wall . J. Colloid Interface Sci. 121 , 254 –264 . (10.1016/0021-9797(88)90429-8 ) 5 Bhutani G , Khandekar S , Muralidhar KV 2012 Contact angles of pendant drops on rough surfaces. In Proceedings of the 39th National Conf. on Fluid Mechanics and Fluid Power, Surat, Gujarat, India, 13–15, December . 6 de la Madrid R 2015 Comparison of the lateral retention forces on sessile and pendant water drops on a solid surface . Am. J. Phys. 83 , 531 –538 . (10.1119/1.4906413 ) 7 de la Madrid R , Garza F , Kirk J , Luong H , Snowden L , Taylor J , Vizena B 2019 Comparison of the lateral retention forces on sessile, pendant, and inverted sessile drops . Langmuir 35 , 2871 –2877 . (10.1021/acs.langmuir.8b03780 )30724570 8 De Coninck J , Fernandez-Toledano J , Dunlop F , Huillet T , Sodji A 2021 Shape of pendent droplets under a tilted surface . Physica D 415 , 132765 . 9 Furmidge C 1962 Studies at phase interfaces I. The sliding of liquid drops on solid surfaces and a theory for spray retention . J. Colloid Sci. 17 , 309 –324 . (10.1016/0095-8522(62)90011-9 ) 10 ElSherbini A , Jacobi A 2006 Retention forces and contact angles for critical liquid drops on non-horizontal surfaces . J. Colloid Interface Sci. 299 , 841 –849 . (10.1016/j.jcis.2006.02.018 )16542670 11 Dussan VEB , Chow RTP 1983 On the ability of drops or bubbles to stick to non-horizontal surfaces of solids . J. Fluid Mech. 137 , 1 –29 . (10.1017/S002211208300227X ) 12 Pitts E 1974 The stability of pendent liquid drops. Part 2. Axial symmetry . J. Fluid Mech. 63 , 487 –508 . (10.1017/S0022112074001741 ) 13 Sumesh PT , Govindarajana R 2010 The possible equilibrium shapes of static pendant drops . J. Chem. Phys. 133 , 144707 (10.1063/1.3494041 )20950030 14 Fatollahi AH 2012 On the shape of a lightweight drop on a horizontal plane . Phys. Scr. 85 , 045401 (10.1088/0031-8949/85/04/045401 ) 15 Hajirahimi M , Mokhtari F , Fatollahi AH 2015 Exact identities for sessile drops . Appl. Math. Mech. 36 , 293 –302 . (10.1007/s10483-015-1916-6 ) 16 Kostoglou M , Karapantsios T 2020 An analytical two-dimensional linearized droplet shape model for combined tangential and normal body forces . Colloids Interfaces 4 , 35 (10.3390/colloids4030035 ) 17 Kostoglou M , Karapantsios T 2019 Contact angle profiles for droplets on omniphilic surfaces in the presence of tangential forces . Colloids Interfaces 3 , 60 (10.3390/colloids3040060 ) 18 Ros-Lpez I , Evgenidis S , Kostoglou M , Zabulis X , Karapantsios T 2018 Effect of initial droplet shape on the tangential force required for spreading and sliding along a solid surface . Colloids Surf. A Physicochem. Eng. Asp. 549 , 164 –173 . (10.1016/j.colsurfa.2018.04.004 ) 19 Fatollahi AH , Hajirahimi M 2013 Making sessile drops easier. (http://arxiv.org/abs/1304.6366) . 20 Brakke K 2013 Surface evolver manual. Susquehanna University See http://facstaff.susqu.edu/brakke/evolver/evolver.html . 21 Berim G , Ruckenstein E 2020 Bond number revisited: axisymmetric macroscopic pendant drop . Langmuir 36 , 6512 –6520 . (10.1021/acs.langmuir.0c00878 )32423220