## Novel method for determination of contact angle of highly volatile liquids

K Nilavarasi 1, a) and V Madhurima 1, b) Department of Physics, School of Basic and Applied Sciences, Central University of Tamil Nadu, Thiruvarur, Tamil Nadu, India-610005

(Dated: 2 August 2019)

In this paper, a novel method for the measurement of equilibrium contact angle of highly volatile binary liquids is proposed. The proposed method, which combines finite element method and energy equilibration, is able to calculate the solid-liquid contact area. The calculated solid-liquid contact area can then be used to estimate the equilibrium contact angle. Using the proposed approach, the contact angles of binary liquid droplet on a microgrooved and smooth polycarbonate substrate were calculated. The proposed method can be an efficient tool for finding the contact angle of all liquids (both volatile and non-volatile).

## I. INTRODUCTION

Wettability becomes crucial for many industrial and scientific applications such as painting/coating 1-4 , surface chemistry 5 , oil recovery 1-4 and so on. Wettability describes the balance of three inter-facial interactions namely, solidliquid ( g sl ), liquid-vapor ( g lv )and solid-vapor ( g sv ) inter-facial interactions 6-8 . The balance between such three interactions are expressed by Young's equation 9-11 ,

$$\gamma _ { l v } c o s \theta = \gamma _ { s v } - \gamma _ { s l } \quad \quad \quad ( 1 ) \quad \text{liq}$$

where q is the contact angle. Although the surface tension of the liquid can be measured experimentally with satisfactory accuracy, the solid-liquid inter-facial tension cannot be measured directly, and therefore the wettability is usually described by the contact angle, which is the angle formed between the solid-liquid interface and the liquid-vapor interface. The contact angle is measured experimentally using a standard approach called sessile drop method, in which the camera focuses the liquid droplet placed over a solid substrate and the geometry of the droplet is used to obtain the contact angle 12 .

Apart from sessile droplet method, there are also other methods of measuring contact angle directly which includes "tilting plate" method 13,14 , captive bubble method, etc., 15-17 . Although the measurement of contact angle from these methods is relatively straightforward, there are several issues that require attention while using volatile liquids: 1. The inherent inaccuracy of the direct measurement techniques and 2. Simultaneous variation in the contact area and contact angle of the liquid over the solid surfaces 18 . There have been numerous studies reporting contact angles for a variety of liquids, binary systems, etc., In most such studies, the liquids used are less volatile and the binary system used contains water as one of its moiety.

a) Electronic mail: nilavarasikv@gmail.com

b) Electronic mail: madhurima@cutn.ac.in

In the present paper, an attempt is made to address the problem of measuring contact angle of highly volatile liquids. Here, equilibrium droplet shape of highly volatile binary liquids on the horizontal smooth and constrained surfaces are simulated to obtain the contact angle of binary liquids on solid surfaces. This 3D-drop shape model is used to numerically analyze the shape, contact area and contact angle of the liquid droplets over the solid surface. The effect of variation of surface tension and surface roughness on the drop shape and apparent contact angle is examined. The liquids used in the present study includes water, methanol, ethanol and different concentrations of the ethanol-methanol binary system. The simulated results are validated with the experimentally obtained data. The present study shows that the proposed method can be an efficient tool for finding the contact angle of all liquids (both volatile and non-volatile).

The contact angle measurements are also useful in determining the size of the pores formed as a result of self-assembly of condensed liquid droplets. Since, pore size is a measure of diameter of a triple phase contact line. In the present study, an attempt is also made to compare the experimental pore size of the self-assembled droplet patterns with the simulated droplet size.

## II. SIMULATION DETAILS

The equilibrium shape of a liquid drop is achieved through energy minimization. In the present study, the finite element method and gradient descent method is used to evolve the surface toward minimal energy 19,20 . The liquid droplet equilibration is achieved by minimizing various energies namely, surface tension, gravitational energy, etc., involved in the defined system. The initial geometrical parameters, energies and constraints involved in the system are given as inputs and the program minimizes the total energy of the system by modifying the surface geometry according to the defined parameters and constraints. The droplet size obtained from simulated stable equilibrium shape of various liquid droplets over specific smooth and constrained surfaces are compared with those from experiments.

FIG. 1. System evolving from initial geometry.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-3b16d645-7ae5-4496-bfcf-dcec3f9a2fee.jpg)

In the present work, the Surface Evolver is used to simulate the drop of different liquids onto a smooth and constrained surface. The gravitational effect is negligible in the present case and hence do not influence the results significantly. Therefore, gravitational energy was not taken into account in our model. The surface tension of the droplet and the inter-facial energy are specified in the data file. The initial geometry and the shape after successive evolutions are shown in Figure 1. The free energy of the system is expressed as 21 ,

$$G / \gamma _ { l v } = A _ { l v } - \int \int \cos \theta d A \quad \quad ( 2 ) \quad \ 1$$

and the contact angle ( q ) is defined by Young's equation 9-11 . Here in equation 4, g lv refers to surface tension of the liquid and Alv refers to the liquid-vapor contact area. The drop volume specified in this work is 9 m l . The bottom face of the droplet is constrained to move and this boundary condition is considered to be responsible for obtaining the shape of the droplets. For the sake of convenience, the constraints are specified to the edges which defines the three phase contact line. The successive refinements and steps concerning energy minimization computes the equilibrium shape of the system. Figure 1 gives the illustration of the steps involved in the Surface Evolver simulation.

## III. RESULTS AND DISCUSSION

## A. Smooth surfaces

Measuring the contact angle of the highly volatile liquid is experimentally difficult 5,15 . Hence, in the present study, Surface Evolver is used to model the wetting behavior of the ethanol-methanol binary liquid drop on the smooth surface. For this solid-liquid inter-facial tension calculated from experiments is used an input parameter instead of conventionally used contact angle. The variation of solid-liquid inter-facial tension for various concentration of methanol is shown in Figure 2.

The shape of the surfaces after the minimization of surface

FIG. 2. Variation of solid-liquid inter-facial tension of smooth surface for various concentration of methanol in ethanol.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-7da4acb7-6f74-44aa-8a45-6c45687d2eca.jpg)

free energy are shown in Figure 3 and 4. It is observed from the figure that the variation of concentration of binary liquid affects the width and height of the liquid drops on the smooth surface. Hence the shapes calculated from Surface Evolver helps in understanding the variation of pores shapes with varying concentration of methanol 22 .

The solid-liquid inter-facial tension for all concentration of ethanol-methanol binary system is calculated following ideal Raoult's law 23,24 and using the obtained values, the drop shapes are simulated. The results of Surface Evolver simulation are compared with the experimentally observed profile of pore shapes. The simulation results show a similar trend as observed in experimental results 22 . The complex variation at low concentration of methanol as seen in experiments is also observed in simulation results. To further understand the correlation between the concentration of methanol and the pore size, diameter of the droplets are calculated. From the obtained diameter of drops, normalized diameter is calculated. The normalized diameter versus concentration of methanol is also plotted (shown in Figure 5). Comparison is made between the normalized diameter calculated from simulation and the normalized pore size obtained from experiments and

FIG. 3. Variation of droplet diameter on smooth surfaces simulated using Surface Evolver .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-06e7e59f-ba8a-4308-9d55-587c2c26134d.jpg)

## is shown in Figure 5.

Further, the deviation in diameter obtained from experiment and simulation from the ideal case (diameter calculated from simulations using inter-facial tension calculated from Raoult's law) is calculated and is shown in Figure 6. It is observed that except at the concentration of 21% and 40% of methanol, all other concentrations showed a positive value. This concentration 21% of methanol is the concentration where the complexity in intermolecular interactions is observed 25 and 40% of methanol is the concentration, where the presence of strong hydrogen bonds were reported 25,26 .

The difference between the simulation and experimental diameter is also calculated and is shown in Figure 7. Positive deviation is observed for all concentration of methanol except at 21% of methanol. This shows that the simulated results are in good agreement with the experimental data i.e., the relative error between experimental and simulated values are small. The negative value at 21% occurs in the region where the complex behavior in inter-molecular interactions is observed. For the sake of comparison, the liquid-vapor inter-facial tension obtained from experiments and calculated from Raoult's law is plotted and shown in Figure 8. It is clear from the figure that the ethanol-methanol binary system is far from ideal.

FIG. 4. Variation of droplet diameter on smooth surfaces (From the top a)12% b) 14% c) 16% d)18% e) 19% f) 20% g) 21% h)22% i) 24% j) 26% k) 28% of methanol in ethanol-methanol binary system).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-4f5d6d27-dd6e-45bf-ab97-4372c2437334.jpg)

FIG. 5. Comparison of normalized diameter on smooth surface calculated from simulation and experiments.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-6c0789e2-f98a-4675-8ae0-8795deb87be0.jpg)

The total energy and total area of the system for various concentration of methanol are calculated from the simulations and are shown in Figure 9. The total energy and total area show a similar trend as seen from experimental result at the

FIG. 6. Deviation in experimental pore diameter and simulated droplet diameter from the ideal case (For smooth surface).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-b0a4b2cf-45b7-4c55-8cbc-b442c83ae179.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-6d569d9f-f40f-4b0b-9630-aaa4a58584e4.jpg)

FIG. 8. Surface tension of liquids obtained from experiments and Raoult's law.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-3789ceec-6124-4fd4-840a-929a74802564.jpg)

FIG. 7. Difference between diameter (on smooth surface) obtained from simulation and experiment.

lower concentration of methanol in ethanol-methanol binary system.

Hence in conclusion, the Surface Evolver results is found to corroborate the experimental findings. Surface Evolver is proved to be a tool for calculating the contact angle of highly volatile liquids for which experimental measurement of contact angle is difficult.

## B. Constrained surfaces

The solid-liquid inter-facial tension for constrained surfaces are calculated from experiments and used as an input parameter instead of conventionally used contact angle. The variation

FIG. 9. Variation of total energy and total area of the system (liquid droplet on smooth surface) for various concentration of methanol in ethanol.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-6367c56f-5f1f-4dc0-9415-f37362e6448a.jpg)

of solid-liquid inter-facial tension for various concentration of methanol is shown in Figure 10.

The shape of the droplet on constrained surfaces after the energy minimization are shown in Figure 11 and 12. It is observed from the figure that the variation of concentration of binary liquid affects the width and height of the liquid drops on the constrained surface. Hence the shapes calculated from Surface Evolver helps in understanding the variation of pores shapes with varying concentration of methanol. It is also observed that the constraints on the surface has significant influence in the wetting behavior of the surface.

Similar to smooth surfaces, the solid-liquid inter-facial tension for all concentration of ethanol-methanol binary system is calculated following ideal Raoult's law and using

FIG. 10. Variation of solid-liquid inter-facial tension of constrained surfaces with various concentration of methanol in ethanol.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-9e943d8a-8923-4b3d-939e-ea7cec38319c.jpg)

FIG. 11. Variation of droplet diameter on constrained surfaces (From the top a) 0% b) 10% c) 20% d) 30% e) 40% f) 50% g) 60% h) 70% i) 80% j) 90% k) 100% of methanol in ethanol-methanol binary system).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-3d48474d-b8a6-49a2-ba94-5e12407931f8.jpg)

the obtained values the drop shapes are simulated. The results of Surface Evolver simulation are compared with the experimentally observed profile of pore shapes on constrained surfaces. The simulation results show a similar trend as observed in experimental results. The complex variation at low concentration of methanol is also observed in the simulation results. To further understand the correlation between the concentration of methanol and the pore size, diameter of the

FIG. 12. Variation of droplet diameter on constrained surfaces (From the top a) 12% b) 14% c) 16% d) 18% e) 19% f) 20% g) 21% of methanol in ethanol-methanol binary system).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-863db705-1e32-40ba-b9be-510fbff0d53f.jpg)

droplets are calculated. From the obtained diameter of drops, normalized diameter is calculated. The normalized diameter versus concentration of methanol is also plotted (shown in Figure 13). Comparison is made between the normalized diameter calculated from simulation and the normalized pore size obtained from experiments and is shown in Figure 13.

FIG. 13. Normalized diameter of droplet on constrained surfaces obtained from simulation and experiment.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-ab8e0513-d0fc-42f5-8534-e3d94ef33fac.jpg)

Further, the deviation in diameter obtained from experiment and simulation from the ideal case (diameter calculated from simulations using inter-facial tension calculated from Raoult's law) is calculated and is shown in Figure 14. Positive deviation is observed for all concentration of methanol except at 21% of methanol and 100%, where the complexity in

FIG. 14. Deviation in experimental pore diameter and simulated droplet diameter from ideal case (constrained surfaces).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-945f0389-f25b-4120-a460-79cfe20afaca.jpg)

properties is observed. At certain concentrations, the datas obtained from experiments and simulations doesnot correlate. This is due to the presence of chloroform in the experiments which is not taken into account in simulations.

The difference between the simulation and experimental diameter for constrained surfaces is also calculated and is shown in Figure 15. It is observed that the deviations is observed to be positive for low concentration of methanol. The higher concentration shows a negative deviation, which may be due to the possibility of complex interactions as observed in 22,25 . For the sake of comparison, the liquid-vapor inter-facial tension obtained from experiments and calculated from Raoult's law is plotted and shown in Figure 8. It is clear from the figure that the ethanol-methanol binary system is far from ideal.

The total energy and total area of the constrained system

FIG. 15. Difference between the simulated and experimental diameter for constrained surfaces.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-6593ed51-f2cf-40cb-a110-242503cfbe30.jpg)

for various concentration of methanol are calculated from

FIG. 16. Variation of total area and total energy of the system (liquid droplet on constrained surfaces) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00246v1-figure-24741a07-9d25-4652-95fd-137bdc28c716.jpg)

the simulations and are shown in Figure 16. The total energy and total area show a similar trend as seen from experimental result at the lower concentration of methanol in ethanol-methanol binary system.

Similar to the simulations obtained for smooth surfaces, the Surface Evolver results of constrained surfaces is found to corroborate the experimental findings. Surface Evolver is proved to be a tool for calculating the contact angle of highly volatile liquids for which experimental measurement of contact angle is difficult.

## IV. CONCLUSION

The contact angle of highly volatile liquid droplets on horizontal smooth and constrained surfaces are calculated. This is performed by simulating the equilibrium liquid droplet shape on the specified surfaces and numerically analyzing the contact area and contact angle of the liquid droplets on the surfaces. The effect of variation of surface tension and surface roughness on the drop shape and apparent contact angle are examined through the simulating the droplets of water, methanol, ethanol and various concentration of ethanol-methanol binary system. The simulated results are validated with the experimentally obtained data. The present study shows that the Surface Evolver can be an efficient tool for finding the contact angle of all liquids (both volatile and non-volatile).

1 K.N. Prabhu, P. Fernandes, G. Kumar, Mater. Des. 2, 297 (2009)

- 2 X. Zhao, M.J. Blunta, J.J. Yao, Pet. Sci. Technol. Eng. 71, 169 (2010)

3 Y.Q.Wang, H.F. Yang, Q.G. Hang, L. Fang, S.R. Ge, Adv.Mater. Res. 154155, 1019 (2010)

4 M. Sakai, T. Yanagisawa, A. Nakajima, Y. Kameshima, K. Okada, Langmuir 25, 13 (2009)

5 H. B. Eral, D. J. C. M. 't Mannetje, and J. M. Oh., Colloid Polym. Sci., 291, 2:247-260, (2013).

- 6 Robert J. Good, Journal of Adhesion Science and Technology, 6:12, 12691302, (1992).
- 7 Gregory J. Merchant, Physics of Fluids A: Fluid Dynamics, 4, 477 (1992).
- 8
- L. M. Hawking, Physics of Fluids, 7, 2950 (1995).
- 9 T. Young, Philos. Trans. R. Soc. Lond. 95, 65 (1805)
- 10 A. M. Dupre, Theorie Mechanique de la Chaleur, Gauthier-Villars, Paris, 369, (1869).
- 11 M. E. Schrader, Langmuir, 11, 3585-3589, (1995).
- 12 W.C. Bigelow, D.L. Pickett, W.A.J. Zisman, Coll. Sci. 1, 513 (1946)
- 13 G. MacDougall, C. Ockrent, Proc. R. Soc. 180A, 151 (1942)
- 14 C.W. Extrand, Y. Kumagai, J. Colloid Interface Sci. 184, 191 (1996)
- 15 K. L. Mittal and R. Jaiswal, Particle Adhesion and Removal, Scrievener publishers, (2014).
- 16 A.F. Taggart, T.C. Taylor, C.R. Ince, Trans. Am. Inst. Min. Metall. Pet. Eng. 87, 285 (1930).
- 17 W. Zhang, M. Wahlgren, B. Sivik, Desalination, 72, 263 (1989)
- 18 Fowkes, Contact Angle, Wettability, and Adhesion Advances in Chemistry; American Chemical Society: Washington, DC, (1964)
- 19 K. A. Brakke, The Surface Evolver, Exp. Math. 1, 141, (1992)
- 20 K. A. Brakke, 2008, Surface Evolver manual: version 2.30, http://www.susqu.edu/brakke/aux/downloads/manual230.pdf
- 21 Y Chen, B He, J Lee, N A Patankar, Journal of Colloid and Interface Science, 281, 458-464, (2005)
- 22 K. Nilavarasi and V. Madhurima, Eur. Phys. J. E, 41: 82, (2018).
- 23 D. A. McQuarrie, J. D. Simon, University Science Books, (1997).
- 24 E. B. Smith, Clarendon Press, Oxford, (1993).
- 25 K. Nilavarasi, Thejus R Kartha and V. Madhurima, Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 188, 301-310, 2018.
- 26 V. H. Alvarez, S. Mattedi, M. Iglesias, R. Gonalez-Olmos and J. M. Resa, Physics and Chemistry of Liquids, 49:1, 52-71, (2011).
K Nilavarasi V. Madhurima
Novel method for determination of contact angle of highly volatile liquids
In this paper, a novel method for the measurement of equilibrium contact
angle of highly volatile binary liquids is proposed. The proposed method, which
combines finite element method and energy equilibration, is able to calculate
the solid-liquid contact area. The calculated solid-liquid contact area can
then be used to estimate the equilibrium contact angle. Using the proposed
approach, the contact angles of binary liquid droplet on a microgrooved and
smooth polycarbonate substrate were calculated. The proposed method can be an
efficient tool for finding the contact angle of all liquids (both volatile and
non-volatile).