Gibbsian Thermodynamic Study of Capillary Meniscus Depth

In the presence of gravity or other external fields, liquid surface curvature deviates from a spherical shape and the surface configuration can be found by numerical integration of the Young–Laplace equation and the typical initial point for integration is the apex of the interface. The meniscus shape in large Bond number systems, which have the central portion of the interface flattened, cannot be determined with the apex as the initial point for integration. Here we find the depth of capillary menisci by considering an initial point for integration to be at the three-phase-contact-line (TPCL) and evaluate the curvature at the TPCL by free energy analysis and inspect the effect of different parameters on the interface shape. A new parameter—which is the deviation of equilibrium curvature at the TPCL from the spherical shape (SR)—is introduced and inspected and it was found that at a Bond number of 13 the maximum deviation, approximately 0.8 of spherical curvature, takes place while for large enough Bond numbers the curvature at the three-phase contact line is near the spherical shape (0.95 < SR < 1). A potential application of this approach is to measure the capillary rise at the TPCL to find the surface tension in high Bond number systems such as those with low surface/interfacial tensions.

Finding the equilibrium shape of fluid interfaces has been frequently studied due to its importance in surface science [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17] . The fluid surface tension and the contact angle that the fluid makes with a solid are two important surface parameters that directly relate to the fluid interface shapes. Sessile drops, pendent drops, bubbles, foams and capillary menisci are some of the important fluid interfaces that we are involved with ubiquitously. The Young-Laplace equation is the most famous relation that governs the equilibrium shape of fluid interfaces. It connects a capillary pressure (the inner and outer phase pressure difference across the curved interface) to the interfacial tension between the two fluid phases and the geometric property of the interface-its mean curvature. When two fluid phases meet a solid phase, the competition between their surface energies leads to the formation of a capillary meniscus.
The calculation of equilibrium shapes of capillary menisci has been done for two main reasons: (a) determination of surface/interfacial tensions and contact angles, and (b) determination of capillary forces.
The first calculation of fluid interface shape is attributed to Bashforth and Adams 18 . They performed numerical integration of the Young-Laplace equation for sessile drops. They non-dimensionalized the Young-Laplace equation with b-the mean curvature of the surface at the apex-and used the apex as the starting point for numerical integration and provided the results of their calculation in several tables for different values of ρ γ Δ gb LV 2 where g, Δρ and γ LV are gravitational acceleration, liquid-vapor density difference, and liquid-vapor surface tension, respectively. Later, Padday and Pitt 14 completed their work by changing the starting point to be at turning angle ϕ = 270° and non-dimensionalized the equation with surface curvature at this starting point. Padday and Pitt mentioned that for menisci that never intersect the axis of symmetry, such as liquid bridges between two solids, the mean curvature at the apex has no real significance and that this difficulty can be overcome by moving the integration starting point. Bucher 1 published an excellent review for interface configurations and classified these systems into sessile and pendent drops, emergent and captive bubbles, fluid bridges and holms (the bridging fluid between a solid and a bulk fluid). He summarized the integration starting point for these interfaces, and for sessile and pendent drops and emergent and captive bubbles the apex of the interface is chosen to be the initial point for integration. Rothenberg et al. 5 minimized an objective function-the difference between the experimental and theoretical shapes of the sessile and pendent drops-and as a result of this optimization obtained the surface tensions and contact angles of fluid interfaces. Grzybowski et al. 10 considered several floating objects of approximately millimeter size and different hydrophobicities at the interface of water and perfluorodecalin. They calculated the interface shape of menisci by two methods: analytic solution of the Young-Laplace equation and use of SCIEnTIFIC REPORTS | (2019) 9:657 | DOI: 10.1038/s41598-018-36514-w the finite element method in the case of no analytical solution with the former method. They found the energy profile and capillary forces of the system from the calculated interface shape. Asekomhe and Elliott 19 computed the shapes of fluid menisci used to infer line tensions from capillary rise in a conical tube and showed that considering the effect of gravity on the liquid-vapor interface shape changed the inferred solid-liquid-vapor line tension values by up to 50%. Danov et al. 2 investigated the capillary meniscus around an electrically charged particle deforming a fluid interface and found that the meniscus deformation for small enough particles is the sum of deformations due to gravity and electric forces individually. Kuchin et al. 11 carefully examined the shape of the capillary meniscus in the transition zone-where the meniscus meets a thin liquid film at the wall so that both the disjoining/conjoining pressure and the capillary pressure play important roles in deforming the meniscus shape. They mentioned that disjoining/conjoining pressure is the dominant surface force in the case of thin liquid films and capillary pressure is the prevailing force in the spherical part of the meniscus. Rayleigh 20 investigated the capillary rise method in order to measure surface tension and proposed two approximate relations: one for systems with Bond numbers smaller than 0.04 and one for systems with Bond numbers larger than 0.04.
Chatterjee 21 investigated two coupled menisci in a capillary and by fitting the upper meniscus curvature to the experimental data predicted the shape of the bottom meniscus. Ratcliffe et al. 17 explored the shape of a drop trapped in a constriction. They solved the Young-Laplace equation for the fluid-fluid interface with the constriction imposing its shape on the solid-liquid interface. Lubarda and Talke 12 assumed ellipsoidal shapes for sessile drops in the presence of gravity and did the minimization of free energy to get the equilibrium height and droplet spreading. They compared their ellipsoidal results with the numerical solution and found that the ellipsoidal shape assumption is accurate for droplets with contact angles smaller than 120° and for droplets with sizes on the order of the capillary length.
Hartland and Hartley calculated the external meniscus shape of fluid which has a horizontal part far from the deforming body 22 . They considered the vertical axis of the deforming body and the horizontal axis along the horizontal part of the meniscus to be the z and x axes of their system. In order to be able to solve the problem they simplified the problem by considering dz/dx ≪1 and found the solution to be a Bessel function. They further did the integration by starting the integration from arbitrarily small ϕ and tabulated their results. In order to use their calculation approach a parameter must be fitted to experimental measurement of the system of interest.
Calculation of meniscus shape and specifically the shape of a liquid droplet is also of great interest in optical science 3,23,24 . Due to the symmetry of liquid droplets they act like lenses, and today they attract attention as adaptive lenses especially in miniaturized systems where surface forces are dominant. Electrowetting and dielectrophoretic operations are two manipulations that can change meniscus curvature and consequently the focal length of these lenses. Ren et al. 3 did several experiments for small droplets and found that even for vertical lenses, the gravity will affect the symmetry of droplets negligibly. As a result, any two immiscible liquids regardless of their density mismatch can be used as an optical lens.
Contact angle measurement in biochemical systems-where only a small volume of matter is available-is of great importance. In order to prevent contact line pinning due to evaporation and also to avoid contamination of the biochemical material, instead of using the sessile drop method, microscopic imaging of the capillary meniscus in a retainer such as a test tube, is suggested 24 . However, in case of a cylindrical retainer there are some difficulties due to image distortion. Cheong et al. showed that knowing the shape of the meniscus (radius and height) is enough to measure the contact angle. However, they assumed a spherical shape for the meniscus and verified their calculated contact angles with measured ones.
Generally speaking, meniscus shape calculation has been done by various methods and for numerous applications. In small capillaries (small Bond numbers) a capillary meniscus forms part of a sphere and the depth of the meniscus scales with the capillary radius. As the Bond number increases, part of the meniscus becomes flattened due to gravitational effects. In order to find the capillary depth in menisci that are flattened in the center, the initial point for numerical integration cannot be the meniscus apex. As a result, we choose the initial point to be at the three-phase contact line and we inspect the curvature at this point by introducing a parameter which shows the deviation of equilibrium curvature at the three-phase contact line from the spherical shape (SR).

Theory
The system of interest is shown in Fig. 1 which is a capillary partially filled with a single component i as a liquid phase which coexists with its vapor phase. The system is symmetric so the axis of symmetry is the z axis and the three-phase contact line in the two-dimensional system is the vertical location of the x axis. We are going to find the meniscus depth of fluid (z h = h). At equilibrium the system has constant temperature (T L = T V ) and the Young-Laplace equation governs the pressure difference across the fluid interface: where P V and P L are the vapor and liquid pressures, z is the elevation with respect to the three-phase contact line, R 1 and R 2 are the meridional and azimuthal principal radii of curvature of the interface with an axis of rotational symmetry, R m is the mean radius of surface curvature and subscript I indicates that a property is referred to the interface. The following equation is the third equilibrium condition which relates the chemical potential to the elevation: In order to use these equilibrium conditions for our system we have to use the equations of state. If we assume the liquid phase to be an incompressible liquid and the gas phase to be an ideal gas, then we will have: where ∞ v L is the liquid molar volume at the saturation condition and R u is the universal gas constant. Substituting Equation (2) in (3) results in: where ρ L is the liquid phase density. And substituting Equation (2) in (4) leads to: Since we are dealing with small (z − z ref ) and at small x, exp(x) ≈ 1 + x, we can simplify Equation (6) to: where ρ V is the vapor phase density. By subtracting Equation (5) from (7) we will have: If we combine Equation (8) with Equation (1) then it can be used for the liquid-vapor interface and the meniscus shape can be calculated from where R m0 is the mean radius of curvature at the reference point. Following the technique initiated by Bashforth 19 and Adams and used later by Ward and Sasges 8 , the following relations relate the meridional and azimuthal principal radii of curvature: where ϕ is the turning angle. Substituting Equations (10)(11)(12) in Equation (9), the following differential equations govern the interfacial shape: where R m0 is the curvature of the meniscus at the three-phase contact line: We assume the reference point to be at the three-phase contact line (z ref = 0) and that | = , the saturation pressure, and then by making use of Equation (7) ρ . In addition, by non-dimensionalizing the terms x, z, R m0 with respect to the capillary radius i.e. ′ = x x r Equations (13) and (14) will simplify to: where B is the Bond number. When in a system both the capillary force and the gravitational force are appreciable, there is a competition between them. This competition can be evaluated by a non-dimensional number-the Bond number-and it is defined as follows: where L is a characteristic length scale of the system which in capillaries is usually assumed to be the capillary radius. Most of the previous calculations 5,7,8,12,16,19,21,24 use the curvature at the meniscus apex as an iterated boundary condition to solve the differential equations. However, this approach cannot be used to predict the shape of a meniscus with a flat portion at its center which is the case for systems with large Bond numbers. The challenge in numerically solving this problem upon moving the initial point for integration is to devise an integration procedure using an iterated boundary condition and stopping conditions that yields a unique solution. Now, in order to fully overcome this shortcoming, we consider the starting point for integration to be at the three-phase contact line (z = 0). As a result, the equilibrium meniscus shape will be found by numerically integrating the differential Equations (16) and (17)  . However, in order to do the integration we still have one more unknown which is the curvature at the three-phase contact line, ′ R mo . We try to find it by assuming an initial value for ′ R mo and doing the integration. If the correct shape is found should be equal to zero in a smooth way. However, this approach will not lead to a unique ′ R mo hence it is not adequate. Here we found an interval for ′ R mo which satisfies all the equilibrium conditions while in SCIEnTIFIC REPORTS | (2019) 9:657 | DOI:10.1038/s41598-018-36514-w order to find a unique mathematical and physical solution, the solution should both satisfy all the equilibrium conditions and obtain the lowest free energy subject to the system constraints. In order to choose an appropriate ′ R mo from this interval we use the free energy analysis [25][26][27][28][29] . It is worth mentioning that considering curvature length (s) as an integrating variable a similar result will be obtained. In this case the following relations exist between the x, z, ϕ and s: Free energy calculation method. In order to choose an appropriate ′ R m0 from the acceptable interval (for which all conditions for equilibrium are satisfied) we analyze the free energy of the system to find which ′ R m0 results in the lowest system free energy among them and hence is the stable ′ R m0 . Figure 2 shows the system of interest which consists of liquid and vapor phases and LV, SL and SV interfaces. There are two meniscus profiles in this figure which we will name states A and B. We are going to analyze the free energy difference of F B − F A near equilibrium where both cases of A and B have the same volume of water with two different values of ′ R m0 . By changing the value of ′ R m0 the corresponding F B is changed and F A acts as a reference free energy. x h and z h are the values of x and z at the point closest to the capillary wall at which the turning angle ϕ = 0. In the case of a meniscus without a flat portion this point would be at the capillary center, but for a meniscus with a flat portion this point is not at the center but rather at the furthest extent of the flat portion. The reference case A is chosen in a way that it has the lowest z h . Since both states A and B have the same volume of water with two different ′ R m0 , Y is the displacement of the origin of the coordinate system of meniscus B with respect to the origin of the coordinate system of meniscus A due to the liquid volume being the same for A and B: In Fig. 2 The system has constant volume therefore Helmholtz free energy (F) plus potential energy (E) 30 governs the system: where for the bulk volume (liquid and vapor phase) the intensive Helmholtz energy (f) is given by the following and n is the number of moles per unit volume of bulk phases. Using Equations (3-7) we will have: For the interfaces (j = SL, SV and LV) the intensive Helmholtz energy is given by: where n is the adsorbed number of moles per unit area and since we neglect any adsorption effects in the system, for all the three interfaces The intensive potential energy (e) for the bulk phases is: At equilibrium the Young equation governs the contact angle at the three-phase contact line:

SV SL LV
By having a fixed coordinate system of the reference case for all menisci and considering menisci displacement with respect to that fixed coordinate system, based on Equations (26) and (27) we will have for bulk phases j = L,V: For more clarity we separate the total energy difference into three parts: with a portion due to differences in interfacial free energies, ΔF surf a portion due to differences in bulk liquid and vapor phase free energies, ΔF bulk , and a portion due to differences in potential energies, ΔF pot . The surface energy difference is the main part of the total energy difference and is based on the combination of interface energies. By using Equations (24), (25), (28) and (30) and non-dimensionalizing the length scales the surface free energy difference is found to be equal to: The bulk energy difference relates to Helmholtz free energy of the bulk phases. Based on Equations (31) and (32) the f j do not depend on z so we can take them out of the integral. Since V L and total volumes are constant for all the menisci we can simplify the bulk energy as follows: By making use of Equations (3) and (4), the bulk free energy difference term will be equal to zero: SCIEnTIFIC REPORTS | (2019) 9:657 | DOI:10.1038/s41598-018-36514-w The potential energy difference is obtained by: By calculating the total free energy and plotting it for different values of ′ R m0 , the most stable fluid curvature at the three-phase contact line will be evaluated.

Results and Discussion
Liquid water -water vapor system. Here we found the meniscus shape for liquid water in a capillary with capillary radius r and a contact angle of θ. The physical properties of the system corresponding to T = 25 °C are 31  . There is another unknown -′ R m0 -in addition to z′ and x′ which has to be fixed before doing the integration. We use free energy analysis to find the appropriate ′ R m0 , and consequently the meniscus shape and depth. As is shown in Fig. 3 the free energy of the system is plotted for two different capillary sizes and contact angles. It shows that the potential energy is almost negligible especially for smaller capillary sizes and the surface energy almost governs the total energy difference. The trend of free energies shows that the largest acceptable ′ R m0 has the lowest free energy and consequently is more stable. As a result, we can generalized that for any size of capillary, the largest acceptable ′ R m0 is the correct meniscus curvature at the three-phase contact line. In order to verify this generalization, we redid the calculation by considering "s" as the integrating variable. Figure 4 shows the capillary meniscus profile for a capillary with r = 0.1 m and zero contact angle for the liquid water -water vapour system. In Fig. 4(a), ϕ is the integrating variable and the meniscus profile is plotted for a series of ′ R m0 . As is clear in Fig. 4(a) all ′ ≤ . R 0 038 m0 satisfy the conditions for equilibrium while ′ = . R 0 0382 m0 is not acceptable. As a result, the correct reduced meniscus curvature at the three-phase contact line which is the largest acceptable ′ R m0 would be between 0.038 < ′ R m0 < 0.0382. Figure 4(b,c) shows the meniscus profile for the same system while s is the integrating variable. In the former case ′ s c = 0.1 m and in the latter case ′ = . s 0 2 c . It is shown that when 0.038 < ′ R m0 < 0.0382 the ′ x I will tend to zero if ′ s c increases properly. In Fig. 5(a) the appropriate reduced curvature at the three-phase contact line is plotted for various capillary sizes and contact angles. The figure shows that ′ R m0 is decreasing with increasing capillary size and decreasing contact angle. In Fig. 5(b) the curvature at the three-phase contact line (R m0 ) is reduced by θ r cos . When the entire meniscus has the shape of part of a sphere then the spherical curvature is = θ R mo sph r cos , . It can be seen from Fig. 5 going to one as capillary size becomes smaller and smaller. After finding the appropriate ′ R m0 , the calculated meniscus profile can be plotted. Figure 6 shows z′ versus x′ for a capillary with radius of 1 cm for four different contact angles of 0°, 25°, 50° and 75°. It illustrates that at constant capillary size, the higher the contact angle, the smaller the reduced meniscus depth. Figure 7 shows z′ versus x′ for a zero contact angle meniscus at five different capillary radii of 0.1 mm, 1 mm, 5 mm, 10 mm and 100 mm. It can be seen that for small capillaries the meniscus forms part of a sphere and has the largest reduced meniscus depth while for large capillaries the inner majority of the meniscus is flat and the reduced meniscus depth is insignificant compared to capillary radius. Figure 8 shows the highest reduced depth of capillary menisci at various contact angles and different capillary sizes. The highest depth of meniscus for each capillary happens at ϕ = 0, = | ϕ In order to find the extent of deviation of the meniscus curvature from the spherical shape, we define a parameter SR. If the whole meniscus were part of sphere, as was explained in Fig. 5( . This is the case for small enough capillaries. In this case, the maximum depth of the capillary is: However, if only a part of the meniscus around the three-phase contact line forms part of a sphere then the maximum depth of the capillary is: . These cases are shown in Fig. 9. As a result, the param- is an illustration of the deviation of the three phase contact line curvature from the spherical shape. Figure 10 shows the SR versus the Bond number on a logarithmical scale (Fig. 10b) and versus capillary size (Fig. 10a) at various contact angles. It shows that for all contact angles the highest deviation from spherical curvature at the three-phase contact line (the higher SR) happens for a capillary size of around 1 cm or more generally at Bond numbers of around 13. As was shown previously in Fig. 5b, for small capillaries the whole menisci form part of a sphere. Therefore, the SR is equal to 1 for small capillary sizes. By increasing the size of the capillary the three phase contact line curvature deviates more and more from the spherical shape and the SR becomes smaller than one. At some specific capillary size the three-phase contact line curvature has the maximum deviation from spherical and after that the deviation decreases more and more until the SR again approaches approximately to one. For large capillaries, although there is a large flat portion of the interface around the center of the meniscus, there exists approximately spherical curvature at the three phase contact line. Figure 11 shows the capillary depth versus the capillary size at various contact angles.
In order to validate the results, for large enough capillaries we can estimate our result by the meniscus depth facing a vertical plate which is: ver plate LV and was verified by experiment 32 . For small capillaries, as we expect the depth of the meniscus is equal to the capillary radius for zero contact angle and in general it is equivalent to Equation (38). The numerical values for these validations are shown in Tables 1 and 2 and Fig. 12 shows this validation for the liquid water -water vapor system with contact angle of zero.
Water -soybean oil system. Here we make calculations for a meniscus with very different physical properties, a liquid-liquid meniscus rather than a liquid-vapor meniscus. Recently, the meniscus between oil and       water has new important applications such as in slippery lubricant impregnated surfaces 33,34 or microfluidic applications 35 where usually vegetable oils such as soybean oil, silicon oil, or peanut oil are used as a lubricant. We found the meniscus shape for water and soybean oil in a capillary with radius r and contact angle of the watersoybean oil meniscus with the capillary wall of θ measured through the water (Fig. 13). The physical properties of this system 36 are: where ρ W is the density of the   liquid water phase and ρ SO is the soybean oil density. It is found that for both systems of liquid water -water vapor and water -soybean oil the plots of R′ mo cosθ versus the Bond number for different contact angles follow the general trend which is shown in Fig. 14. Figure 15 for the water -soybean oil meniscus is similar to Fig. 10 for the liquid water -water vapor meniscus. The figure illustrates that for either of these systems the maximum deviation from spherical curvature at the three-phase contact line happens around a Bond number of 13, while because of the different system physical properties the maximum deviation for the water -soybean oil system happens at a capillary radius of 2 cm. Since the Bond number includes all the physical properties involved in determining meniscus shape-including the capillary size, surface tensions and density differences between the phases-it is the most comprehensive variable to be investigated in these types of problems. Figure 16 shows the depth of water menisci with an oil phase for various contact angles and different capillary sizes and the trends are similar to those for the water -water vapor system. Figure 17 indicates the reduced capillary depth vs. the capillary size for several contact angles for water -soybean oil menisci. By increasing the capillary size the reduced capillary depth increases while the capillary depth decreases.
Potential application: Measurement of surface tension by the capillary rise method (low surface tension or large capillaries). One of the oldest methods to measure surface tension is the capillary rise method. In cases for which Bo < 0.01 this method is easily accessible. However, usually this criteria is hard to achieve and for larger Bond numbers, Sugden 37 and later Lane 38 -based on the Bashforth and Adams work-prepared a table/graph of r/a versus r/b where r, a and b are the capillary radius, capillary length 2 LV and curvature of the meniscus at the apex. However, their work is restricted to a contact angle of zero and did not cover the ranges of system parameters for which a portion of the interface is flattened. Importantly this restriction means that the capillary rise method could not be used to measure low surface tension values-more expensive specialized equipment such as the spinning drop tensiometer is typically used, which introduces its own difficulties 39 .      As shown in Fig. 18, the capillary rise can be determined from the three-phase contact line (h rise,t ) or from the meniscus apex (h rise,b ). The latter is usually measured due to the use of the Bashforth and Adams table which tabulated the meniscus curvature at the apex. Although working with the former capillary rise is not usual, it can be used in measuring the surface tension in small capillaries and even in larger capillaries in which the meniscus has a flattened center where h rise,b = 0, which will be the case for low surface tension values. These two capillary rises are related by following equation: In fact, in the case of a flat interface we do not have any capillary rise at the apex (h rise,b = 0) but there is still a capillary rise at the three-phase contact line (h rise,t ≠0). Our approach provides a means to measure the surface tension even in cases that there is a big container, or a low surface tension. Figure 19 shows the inverse of meniscus curvature versus r/a for different contact angles. Figure 19(b) is reproduced from the work of Sugden and Lane's and the meniscus curvature is at the apex (b) while in Fig. 19(a) the meniscus curvature is at the three-phase contact line. It should be noted that since all the menisci with various contact angles have a spherical shape at small sizes-to have them all start at the same position we have multiplied R mo with cosθ. Figure 20 shows the capillary rise versus the Bond number. In Fig. 20(a), the capillary rise at the triple contact line is plotted for various contact angles and in Fig. 20(b) both h rise,t and h rise,b are plotted for a contact angle of zero. As we expect, for smaller capillaries h rise,b and h rise,t coincide since the gravity effect is negligible and the pressure difference throughout the whole interface is the same. However, by increasing the size of the capillary, or decreasing the surface tension, both of which increase the Bond number, the deviation between h rise,b and h rise,t becomes important and in the case of a flatted central part of the interface, although h rise,b is equal to zero and cannot be used for surface tension measurement, h rise,t has a determined value that can be used to measure interfacial tension. Figure 21(a) shows the reduced curvature at the apex (based on the work of Sugden and Lane) and the reduced curvature at the three-phase contact line. By increasing the Bond number the apex curvature increases more and more to reach the flat interface-with radius of curvature going to infinity- (Fig. 21-b-1) while the curvature at the three-phase contact line increases at such a slower rate (Fig. 21-b-2) that it leads to the decrease in its reduced curvature.

Conclusion
In this study the capillary meniscus depth is evaluated for various capillary radii and different contact angles for two systems with different properties: the water liquid -water vapor system and the water -soybean oil system. In order to integrate the Young-Laplace equation numerically, we have changed the typical starting point for the integration from the apex of the interface as is used in typical interface shape calculations to be instead at the three-phase contact line, since the apex point cannot be used as the starting point for integration to predict the meniscus shape for systems with larger Bond numbers where the meniscus has some central portion flattened. The curvature at this new integration starting point-at the three-phase contact line-is evaluated by free energy  calculation. It was found that the curvature at the three-phase contact line would be the highest acceptable curvature and this value correspond to the smallest system free energy. The same result is obtained if the curvature length (s) is the integrating variable instead of turning angle (ϕ) as the integrating variable. The deviation of curvature at the three-phase contact line from that of a spherical shape is examined and it was found that the highest deviation which is about SR = 0.83 happens at a Bond number of approximately 13 for both systems studied. As a result, if the curvature at the three phase contact line is unknown-which is the pre-requirement to solve the Young-Laplace equation-it can be estimated with spherical shape with the deviation reported here. Finally, using the numerical methods presented here to calculate the shapes of capillary menisci in the two systems for various values of capillary radius and contact angle, insight is gained into the dependence of meniscus depths on the various variables. One of the potential applications of this approach is to measure the capillary rise at the three-phase contact line to find the surface tension in high Bond number system such as systems with large capillaries or low surface/interfacial tensions.