A Steady Method for the Analysis of Evaporation ... - ACS Publications

utilizing a piezoelectric dispenser to feed microscale droplets ( ≈ 9 µm) to a larger evaporating ... Using our steady technique, we studied water ...
0 downloads 0 Views 5MB Size
Subscriber access provided by PEPPERDINE UNIV

Article

A Steady Method for the Analysis of Evaporation Dynamics Ahmet Alperen Gunay, Soumyadip Sett, Junho Oh, and Nenad Miljkovic Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.7b02821 • Publication Date (Web): 26 Sep 2017 Downloaded from http://pubs.acs.org on September 28, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Langmuir is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

A Steady Method for the Analysis of Evaporation Dynamics A. Alperen Günay1, Soumyadip Sett1, Junho Oh1, Nenad Miljkovic1,2,3,*

1

Department of Mechanical Science and Engineering, University of Illinois at Urbana – Champaign, 1206 W Green St, Urbana, Illinois 61801 USA

2

Frederick Seitz Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, USA 3

International Institute for Carbon Neutral Energy Research (WPI-I2CNER), Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan

* Corresponding Author. E-mail address: [email protected]

1 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Droplet evaporation is an important phenomenon governing many man-made and natural processes. Characterizing the rate of evaporation with high accuracy has attracted the attention of numerous scientists over the past century. Traditionally, researchers have studied evaporation by observing the change in the droplet size in a given time interval. However, the transient nature coupled with the significant mass-transfer governed gas-dynamics occurring at the droplet three-phase contact line make the classical method crude. Furthermore, the intricate balance played by the internal and external flows, evaporation kinetics, thermocapillarity, binary-mixture dynamics, curvature, and moving contact lines make the decoupling of these processes impossible with classical transient methods. Here, we present a method to measure the rate of evaporation of spatially and temporally steady droplets. By utilizing a piezoelectric dispenser to feed microscale droplets ( ≈ 9 µm) to a larger evaporating droplet at a prescribed frequency, we can both create variable-sized droplets on any surface, and study their evaporation rate by modulating the piezoelectric droplet addition frequency. Using our steady technique, we studied water evaporation of droplets having base radii ranging from 20 µm to 250 µm on surfaces of different functionalities (45˚ ≤ , ≤ 162˚, where , is the apparent advancing contact angle). We benchmarked our technique with the classical un-steady method showing an improvement of 140% in evaporation rate measurement accuracy. Our work not only characterizes the evaporation dynamics on functional surfaces, it also provides an experimental platform to finally enable the decoupling of the complex physics governing the ubiquitous droplet evaporation process. Keywords: Evaporation, Droplet, Hydrophobic, Hydrophilic, Contact Angle, Steady, Steady State, Superhydrophobic, Contact Line, Non-Condensable Gases.

2 ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

1. Introduction The evaporation of droplets is an important phenomenon for both natural and man-made processes such as the hydrological cycle1-4, oil-water separation5-8, liquid hydrocarbon combustion9-12, phase change heat transfer13-18, and power generation19-22. Thus, scientists have been studying droplet evaporation over the past century23-31. The main physical phenomenon characterizing evaporation is the rate of evaporation, which determines the amount of liquid mass that is converted into vapor in a fixed amount of time. Numerous studies have been performed in an attempt to more accurately characterize the rate of evaporation on functional surfaces15, 32-39. However, the developed techniques typically rely on detecting the non-steady size changes in a droplet for a given amount of time, which include transient effects governed by vapor dynamics outside the droplet, interfacial dynamics at the liquid-vapor interface, and liquid dynamics inside the droplet. In addition, droplet evaporation dynamics are made more complex by the intricate and length-scale-dependent balance played by the bulk flow inside the droplet40-44, evaporation kinetics at the liquid-vapor interface15,

45-46

, thermocapillary stability47-48, and binary-mixture dynamics49-51, making the

decoupling of these processes impossible with classical transient methods. Furthermore, the internal, interfacial, and external flows are complicated further by the presence of moving contact lines and high curvature regions, both of which govern evaporation52. In this study, we develop a steady method capable of decoupling droplet length scale from the evaporation measurement, to examine droplet evaporation dynamics for a range of droplet sizes (20 ≤  ≤ 250 m, where  =  sin , is the base radius) and apparent advancing contact angles (45˚ ≤ , ≤ 162˚). The steady technique involves the feeding of microscale droplets ( ≈ 9 µm, where  is the radius of curvature of the droplet) to a larger evaporating droplet (≫9 µm) residing on a functional surface at a prescribed frequency. At the instant the droplet reaches the required size of interest, the microscale droplet injection rate is modulated 3 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

in order to ensure that the evaporating droplet shape and size are constant. At this instant, the rate of liquid evaporation is balanced by the rate of liquid microdroplet deposition. To ensure stability of our technique, we monitor the droplet evaporation for long times (>120 seconds). We benchmarked our technique with classical non-steady methods of evaporation characterization on functional surfaces, showing significant improvement (≈ 140%) in evaporation rate measurement accuracy with our steady method stemming from the elimination of transient effects on a droplet evaporating in the steady state. Our results show that the rate of droplet evaporation in non-condensable gas (NCG) leaden environments increases linearly with increasing droplet base radius due to the linear dependence of contact line length with base radius. Furthermore, we show that the droplet evaporative flux decreases linearly for increasing contact angles at fixed droplet sizes for both hydrophilic and hydrophobic surfaces. The outcomes of this work not only show the effects of surface functionality, droplet length scale, and contact line dynamics on the droplet evaporation rate, but also establish a powerful characterization platform that can be used to decouple the intricate and length-scale-dependent balance played by internal and external flows, evaporation kinetics at the liquid-vapor interface, thermocapillary motion, and binary-mixture dynamics, and provide avenues for the development of novel, and versatile functional surface technologies for a wide range of applications utilizing droplet evaporation.

2. Experimental Methods 2.1. Experimental Procedure. To benchmark our steady method, and to study the effects of droplet size scale, base area, contact line length, and contact angle, we fabricated six different samples with distinct and wide ranging wetting characteristics, including a gold coated silicon wafer (Si-Au, Figure 1a, , = 45.1˚ ± 5.2˚), polished copper (Cu, Figure 1b, , = 85.0˚ ± 1.9˚), smooth silicon wafer coated with a fluorinated silane (TFTS) (Si-TFTS, Figure 1c, , = 110.6˚ ± 2.8˚), carbon nanotube mats coated with a 90 nm layer of P2i 4 ACS Paragon Plus Environment

Page 4 of 36

Page 5 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

fluoropolymer (CNT1, Figure 1d, , = 130.6˚ ± 2.9˚), carbon nanotube mats coated with a 30 nm layer of P2i fluoropolymer (CNT2, Figure 1e, , = 141.3˚ ± 3.1˚), and HDTS coated copper oxide microstructures (CuO-HDTS, Figure 1f, , = 162.0˚ ± 2.9˚). Note, lower contact angles approaching superhydrophilicity were not studied here due to requirement for wide angle imaging techniques to capture the large base radius. To study droplet evaporation and to eliminate transient effects of moving interfaces, we interfaced a high speed camera (Fastcam SA2, Photron) with a frequency-controlled piezoelectric micro-goniometer (MCA-3, Kyowa Interface Science). All experiments were performed in ambient conditions ( = 25 ± 0.5 ℃,  = 50% ± 5%, where  is the relative humidity). To grow droplets, the piezoelectric dispenser was placed above our functional surfaces with a spacing ranging from 5 to 10 mm (ℎ), and turned on to dispense monodisperse microscale droplets (Figure 2). 2.2. Surface Fabrication. For CNT growth by chemical vapor deposition (CVD), a 20 nm thick Al2O3 diffusion barrier and a 5 nm thick film of Fe catalyst layer were deposited on silicon growth substrates using electron beam deposition. After purging in He/H2 atmosphere in a 2.54 cm quartz furnace for 15 min, the silicon substrate was heated to 750°C in a 2.54 cm quartz furnace and annealed for 3 min under a flow of H2 and He at 400 and 100 sccm, respectively. Randomly aligned CNTs were then grown by flowing C2H4 at 200 sccm for 1 min. These thermally grown CNT had characteristic diameter of  ≈ 7 nm53. CuO nanostructured surfaces were created on commercially available copper (Cu) tabs (99.9% purity, 25 mm x 25 mm x 0.8 mm). Similar to silicon substrates, each Cu tab was thoroughly rinsed with acetone, ethanol, isopropyl alcohol (IPA), and deionized (DI) water. The cleaned tabs were then dipped in 2.0 M Hydrochloric acid solution for 2 min to remove native oxide film on the surface, rinsed with DI water and dried with nitrogen gas. 5 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nanostructured CuO films were formed by immersing the cleaned tabs in hot (96 ± 3 ˚C) alkaline solution composed of NaClO2, NaOH, Na3PO4·12H2O, and DI water (3.75:5:10:100 wt%). This oxidation process leads to initial formation of a thin (≈300 nm) Cu2O layer which reoxidizes to form sharp, knife-like CuO oxide structures ( ≈ 1µm). 2.3. Surface Functionalization. Trichloroperfluorooctyl-silane (TFTS, CAS no. 78560-45-9, Sigma Aldrich) was deposited on silicon growth substrate by liquid phase deposition (LPD). Initially, the obtained silicon substrate was thoroughly cleaned by rinsing it with acetone, methanol, ethanol, isopropyl alcohol (IPA) and deionized (DI) water. After drying the silicon substrate with nitrogen gas, it was plasma cleaned for 10 minutes. For hydrophobic functionalization, a drop of TFTS was added to 500 mL hexane and the resulting solution was vigorously mixed before placing the clean silicon substrate inside the solution. The coated silicon substrate (referred to as Si-TFTS in the text) was removed from the solution after 24 hours and was rinsed with ethanol and DI water and dried in nitrogen gas. Goniometric measurements (MCA-3, Kyowa Interface Science Ltd.) of ≈ 100 nL droplets on TFTS coated silicon substrates showed advancing contact angle of  = 110.6˚ ± 2.8˚. For hydrophilic functionalization, the plasma cleaned silicon substrates were sputter coated with gold having 10 nm thickness (referred to as Si-Au in the text). Goniometric measurements (MCA-3, Kyowa Interface Science Ltd.) of ≈ 100 nL droplets on Si-Au substrates showed advancing contact angle of , = 45.1˚ ± 5.2˚. Hexadecyltriethoxysilane (HDTS, CAS no. 16415-12-6, Sigma Aldrich) was deposited by liquid phase deposition on CuO nanostructured surfaces. Before deposition, the nanostructured CuO tab was dried in clean nitrogen gas. A solution containing a drop of HDTS in 500 mL of hexane was prepared and vigorously mixed, before placing the CuO tab in it. The coated Cu tab (referred to as CuO-HDTS in the text) was removed from the solution after 24 hours, rinsed with ethanol, DI water and dried with clean nitrogen gas before placing 6 ACS Paragon Plus Environment

Page 6 of 36

Page 7 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

it in an atmospheric oven at 80 °C for 90 mins. Goniometric measurements (MCA-3, Kyowa Interface Science Ltd.) of ≈ 100 nL droplets on HDTS coated CuO substrates showed advancing contact angle of , = 162˚ ± 2.9˚. Plasma-enhanced vapor deposition was used for obtaining P2i hydrophobic coatings on the CNT surfaces. Within a vacuum chamber at room temperature and low pressure, the P2i fluoropolymer coating is introduced as vapor and ionized. By varying the time of deposition, two different thicknesses of P2i coating, 90 nm and 30 nm were obtained. Goniometric measurements (MCA-3, Kyowa Interface Science Ltd.) of ≈ 100 nL droplets on the thicker (90 nm, referred as CNT1 in text) and thinner (30 nm, referred as CNT2 in text) P2i coated sample showed advancing contact angles of , = 130.6˚ ± 2.9˚ and , = 141.3˚ ± 3.1˚, respectively. Note, the functionalized CNT samples were used extensively in condensation experiments prior to this study54, resulting in degradation and lower than expected advancing contact angle behavior. 2.4. Surface Characterization. Contact angle measurements of ≈100 nL droplets on all samples were performed using a microgoniometer (MCA-3, Kyowa Interface Science). Field emission scanning electron micrographs (QUANTA FEG 450 ESEM, FEI) were performed on all samples at an imaging voltage of 5.0 kV with 2.0 nm of spot size. 2.5. Experimental Setup. Rate of evaporation measurements were performed on a microgoniometer (MCA-3, Kyowa Interface Science) combined with a high speed camera (Fastcam SA2, Photron) (Figure 2). The samples were horizontally placed on a movable stage with 3D spatial control ability (X-Y-Z). Illumination was supplied by a LED source attached to the dispenser (AITECSYSTEM, TSPA22x8). Droplets were added on the sample through a 2D piezoelectric dispenser (X-Z) with voltage and frequency control. The experiments were carried out at 6 V, 5-120 Hz (), and a sample-dispenser spacing of 5-10 mm (ℎ). Imaging 7 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

was performed at 7×, 12×, 18× and 24× magnification depending on the droplet size. Image acquisition was performed at 200 fps. All the experiments were performed under ambient conditions ( = 25 ± 0.5 ℃,  = 50% ± 5%).

3. Results and Discussion 3.1. Steady Measurement Technique. To study droplet evaporation in ambient conditions, a high speed camera was interfaced with a frequency-controlled piezoelectric microgoniometer. The piezoelectric dispenser was placed 5 to 10 mm above the functional surfaces to grow droplets (Figure 2). The surface-dispenser spacing was chosen to have close control of the dispensed droplet trajectory, which is not possible at larger spacings due to the presence of random air currents which act to deflect microscale droplets. For the spacing range studied (5 to 10 mm), the evaporation dynamics of steady droplets were independent of dispenser to surface spacing. Due to the relatively small size of departing droplets (≈9 µm), the shape of droplets remained spherical during flight. This assumption is justified given that the Bond, Weber, and Capillary numbers are all much less than one (Bo = #$ % & /( ≪ 1, We = #- .& /( ≪ 1, Ca = - ./( ≪ 1, where % is the gravitational constant, . is the dispensing speed of droplets (≈ 1.5 cm/s), - is the vapor dynamic viscosity, #$ and #- are the liquid and vapor densities of water, respectively, and ( is the water liquid-vapor surface tension). For a detailed list of physical property values used in this study, please see Table S1 of the Supplementary Information. As the microscale droplets landed on the functional surface, they began to accumulate and grow to form a single larger droplet (Figure 3a). The camera focal plane was located at the growing droplet mid-plane for the desired optical zoom level. Subsequently, while the droplet was growing, contact angle measurements on the surface were performed in order to characterize functionality via the apparent advancing contact angle. Careful stroboscopic imaging of the dispensing droplets 5 mm beneath the piezoelectric dispenser confirmed the production of spherical and monodisperse droplets ( = 8.77 ± 0.56 8 ACS Paragon Plus Environment

Page 8 of 36

Page 9 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

µm, Figure 2b) coming from the piezoelectric dispenser that were independent of both dispensing frequency (2 Hz <  < 10 kHz) and driving voltage (3-9 V). After surface characterization and droplet formation, the frequency of deposition was modulated in order to find the dispensing frequency at which the stationary droplet size becomes spatially and temporally steady (Figure 3b). Once the droplet was steady, the frequency was noted, as well as the standard deviation from measurements made on multiple locations on the same sample. The frequency of deposition () that resulted in a steady droplet for a certain droplet diameter was then related to the corresponding rate of evaporation 12 = 1, where 1 is the mass of dispensed microscale droplets, 1 = 4/35 6 #$ . For each surface, multiple frequencies were probed in order to evaluate the droplet size and dependent evaporation dynamics. The piezoelectric head is capable of dispensing picoliter-scale droplets at frequencies of 2 – 10,000 Hz during measurement. The potential to introduce vibrations in the evaporating droplet by dispensing droplets for the piezoelectric head, which may cause the contact line to oscillate, was considered. The Ohnesorge number for a water droplet of the size observed during evaporation measurements (diameter 7 ≈ 50 µm) was Oh = $ /:#$ (7 ≈ 0.01. Therefore, viscous effects are neglected and a balance between inertial and surface tension forces is considered. Determining the droplet natural frequency for the first resonant mode in this case55-56:  = ( [6(ℎA, B]/[#D A1 − cos , BA2 + cos , B] ≈ 6 kHz, where , is the apparent advancing contact angle of the droplet residing on the surface, and ℎA, B is a geometry factor57. The experimental data points were obtained with droplet dispensing frequencies of 5 – 120 Hz 20 µm, indicating that the steady approach is suitable for high accuracy characterization of droplet evaporation. It is important to note, due to the finite size of dispensing droplets ( ≈ 9 µm), the minimum recommended steady droplet size is ≈ 20 µm due to the added possibility of incoming dispensed droplets (at any frequency) affecting the local vapor diffusion dynamics surrounding the droplet residing on the surface. For the steady analysis of smaller droplets (< 20 µm), a piezoelectric dispenser having an ejection droplet size of  ≪ 9 µm should be utilized to avoid coupling between the deposition and steady droplets. The steady technique of droplet evaporation characterization presents a few significant advantages when compared to classical transient based methods. The removal of moving liquid-vapor interfaces and three-phase contact lines during evaporation allows for the elimination of liquid-gas momentum exchange (i.e. motion of the liquid-vapor interface) during evaporation and the characterization of a true evaporation rate at a fixed geometry. The elimination of boundary motion is particularly important in evaporation processes where vapor diffusion is a limiting factor, such as that observed in NCGs58-59. Furthermore, the ubiquitous presence of contact angle hysteresis and pinning on surfaces of interest results in difficulty in obtaining consistent evaporation results at a fixed contact angle or droplet base radius. Instead, an amalgamation of data is obtained with varying droplet parameters which 10 ACS Paragon Plus Environment

Page 10 of 36

Page 11 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

must be post processed in order to extract meaningful information. The steady technique on the other hand can probe the effects of droplet size, droplet base area, contact line length, and contact angle, with high precision, high accuracy, and elimination for the need to post-process data in order to obtain limited results. Lastly, the consistent addition of liquid at a constant enthalpy (temperature) from the piezoelectric dispenser ensures that the temperature profile inside the droplet remains steady and that self-cooling due to evaporation does not introduce secondary transient effects during the measurement27. 3.2. Droplet Radius. Figure 4a shows the droplet evaporation rate (12) as a function of the droplet radius of curvature (R) as determined by the steady technique. A linear trend was observed for all surfaces studied. Interestingly, the results show that the droplet contact angle did not play a large role in the rate of evaporation when the considered size parameter is droplet radius of curvature. Due to the inability to control for contact line length and interfacial area for the varying contact angles, we re-plotted our evaporation data as a function of droplet base radius as determined using the droplet radius and the contact angle information. Figure 4b shows the droplet evaporation rate as a function of droplet base radius, showing a linear dependence for each tested sample. However, the contact angle now plays a significant role. As the hydrophobicity of the sample increases (i.e. , increases), the rate of evaporation increases for the same droplet base radius. This was due to the increase in the exposed liquid-vapor surface area as the advancing contact angle increases, resulting in higher heat transfer with the environment60. It is important to note, however, that the rate of evaporation of the hydrophilic materials (Cu and Si-Au) are higher than Si-TFTS, which was hydrophobic, indicating that a balance between droplet conduction, contact line length, and liquid-vapor interfacial area govern the evaporation process. 3.3. Evaporative Flux. In order to have a better understanding of the droplet evaporation dynamics, the evaporative mass flux (1") was studied. The evaporative flux was evaluated by 11 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 36

normalizing the rate of evaporation (12, Figure 4a, b) with the liquid-vapor surface area in order to include the contact angle effects in the analysis (Figure 5a). Examining Figure 5a, the effect of the advancing contact angle on the evaporation process can be understood as a balance between liquid-vapor interfacial area for vapor diffusion and efficient heat transfer through the droplet required to facilitate evaporation. The results show that the evaporative flux increases with decreasing contact angle for the same droplet base radius, indicating that efficient heat transfer between the surface and liquid-vapor interface is a key enabler for efficient vaporization. The thermal resistance due to heat conduction through the droplet can be approximated as61-63: IJKLM ≅

1 , 4O$  sin(, )

(1)

where O$ is the liquid water thermal conductivity. Assuming that vapor diffusion and evaporative cooling limitations are negligible compared to droplet conduction (large radius limit), the evaporative flux can be determined from 1" ~ ∆/IJKLM / & ~ 1/ ~ 1/ , in excellent agreement with the results of Figure 5a at larger length scales ( > 100 µm, for a log-log plot of Figure 5a, please see Figure S1 of the Supplementary Information). In addition, Eq. (1) shows that as the contact angle increases from hydrophilic to hydrophobic, the droplet conduction resistance increases and 1" decreases, as shown in Figure 5a. Note, Eq. (1) is a valid approximation for droplets residing on hydrophobic surfaces, and does not accurately predict thermal resistances of droplets on hydrophilic substrates, which continue to decrease as contact angle decreases, unlike that suggested by Eq. (1). 3.4. Vapor Diffusion. Droplet conduction may not be the only limit governing droplet evaporation. Indeed, vapor diffusion and evaporative cooling at the liquid-vapor interface have been shown to be a critical limit to evaporation in the presence of NCGs58-59 and at high heat fluxes64-65. In order to gain further insight, we compared the experimental results with the

12 ACS Paragon Plus Environment

Page 13 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

power law exponent model which can be utilized to explain vapor-diffusion controlled condensation and evaporation.66 In these experiments, due to the presences of NCGs, we expect vapor diffusion limitations to play a key role in droplet evaporation dynamics at low length scales where droplet conduction is not important ( < 100 µm) as well as at much larger length scales where droplet conduction affects vapor diffusion via the droplet temperature distribution, and the two couple to create non-uniform diffusive fields surrounding the evaporating droplets. At small droplet length scales ( < 100 µm) as the droplet evaporates, direct evaporation of water molecules from the liquid-vapor interface induces a 3D concentration gradient of water molecules inside the vapor surrounding the droplet. For a constant flux of water molecules leaving the water droplet surface, the concentration gradient should vary as 1/ &, resulting in a scaled droplet volumetric evaporation rate of: R  1 = 45 & ~ 45 & T & U 7V& S S 

(2)

where, R is the volume of water droplet, 7V& is mutual diffusion coefficient of water in air. Rearranging Equation 2 reveals: 1" ~ (R/S)/ & ~ 1/ & ~ 1/ & , in excellent agreement with the results of Figure 5a at small length scales ( < 100 µm) (see Figure S.1). Note, the scaling analysis used to obtain Equation 2 is only an order-of-magnitude analysis and therefore is devoid of constants related to the droplet geometry. For larger evaporating droplets ( > 100 µm) which are conduction limited (Eq. 1), the majority of the evaporation occurs at the contact line as observed by the convergence of the evaporative fluxes at large base radii (Figure 5a). The resulting 2D concentration gradient varies as 1/. Thus, scaling analysis yields: R  1 = 45 & ~ 45 & T U 7V& S S  13 ACS Paragon Plus Environment

(3)

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 36

Rearranging Eq. (3) reveals: 1" ~ (R/S)/ & ~ 1/ ~ 1/ , in excellent agreement with the results of Figure 5a at very large length scales ( > 200 µm, see Figure S.1). The power law results reinforce the critical importance that NCGs play in the dynamics of droplet evaporation. 3.5. Natural Convection. To investigate the effects of natural convection with the surroundings of the evaporating droplet, the Rayleigh number (Ra) of the droplet for the experimental conditions was evaluated. As the droplet size has an important effect on the Rayleigh number (~ 6 ), the Rayleigh number had a broad range for the droplets studied here, calculated as67: %X7 6 ∆ Ra = , YZ

(4)

where % is the gravitational constant, X is the thermal expansion coefficient, 7 is the diameter of the droplet, ∆ is the temperature difference between the ambient air and the evaporating droplet, Y is the thermal diffusivity and Z is the kinematic viscosity of air. Utilizing Eq. (4) yields Ra ≈ 0.08 for the largest droplet of interest ( ≈ 500 µm), revealing that the heat transport with the surroundings is governed by conduction since the maximum Rayleigh number was much smaller than the critical Rayleigh number for air (Ra ≪ RaJ = 1708)68-69. 3.6. Evaporative Cooling. Evaporative cooling of droplets has been shown to play an important role in the dynamics of droplet evaporation16-17. During evaporation, the absorption of latent heat results in the cooling of the droplet interface, which in turn reduces the local liquid temperature, vapor pressure, and hence evaporation rate. Evaporative cooling effects are especially more pronounced in high contact angle (> 90°) droplets due to larger liquidvapor interfacial area, leading to increased resistance to evaporation.

14 ACS Paragon Plus Environment

Page 15 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

For droplets undergoing evaporation, the relative effect of evaporative cooling compared to vapor diffusion can be characterized from the dimensionless Evaporative cooling number (Ec), defined as70: Ec =

ℎ\] 7- ^ , O$

(5)

where ℎ\] is the latent heat of evaporation of water, 7- is the diffusion coefficient of water vapor (subscript v) in air (subscript a), ^ = _`a ()/, and _`a () is the saturated vapor concentration of liquid water at temperature . For liquid water at room temperature (295 K) and pressure (101 kPa), Ec ≈ 0.11, below the critical value, EcJa = 1, implying that evaporative cooling had negligible influence on the droplet evaporation experiments conducted here. 3.7. Effect of Contact Angle. To further emphasize the scaling results and the delicate balance between vapor diffusion and droplet conduction, and to help visualize the effect of the contact angle, the evaporative flux data at constant base radii ( = 60 µm and  = 100 µm) were compared (Figure 5b). The results show that droplets having lower contact angle are more sensitive to droplet size effects due to the low conduction resistance (Eq. 1), dominating 3D vapor diffusion resistance (Eq. 2) and mitigation of evaporative cooling due to the flat shape of the droplets16. For droplets having larger contact angles (> 90°), droplet size effects are less pronounced due to the conduction-limit enabling 2D vapor diffusion resistance (Eq. 3) which acts to govern evaporation by contact line length. The results indicate that the two regimes of evaporation governing behavior are delineated between hydrophilic and hydrophobic surfaces, in excellent agreement with numerical results obtained by previous studies16. 3.8. Comparison with Classical Transient Methods. In order to validate our method and to establish a comparison with the traditional derivative based technique, we performed 15 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

additional experiments on the hydrophilic (Cu), hydrophobic (CNT1) and a superhydrophobic surfaces (CuO-HDTS) with the traditional method to obtain rates of evaporation. Briefly, the evaporation rate was calculated from image processing and using the rate approximation of (b/bS ≈ Δ/ΔS). Figure 6 shows the evaporation rate as a function of droplet base radius for the classical and steady approach. The transient method overestimates the rate of evaporation for each case due to the uncertainties associated with it such as variable contact angle during the process, pinning, and interfacial motion. In fact, pinning to the surface was observed during the experiments which led to a constant droplet base radius for most of the cases with varying contact angle; hence, an accurate calculation of the rate of evaporation cannot be performed due to this effect. Furthermore, the transient method experiments are conducted by turning the droplet flow off, and allowing evaporation to take place as the droplet recedes. Thus, the rate of evaporation is obtained for receding contact angle values. On the contrary, the steady method developed here utilizes continuous droplet flow and is hence characteristic of the droplet evaporation rate in the advancing state with little variability due to pinning. A secondary drawback of the classical transient method is need to measure time as an additional parameter in the calculations since transient effects are observed due to the changes in the size and shape, especially for a small droplets. The need for temporal and spatial information can therefore increase error in the calculation. However, for the steady technique, we eliminate the need for including time as a parameter and only match the frequency of deposition to the frequency of evaporation to get a higher fidelity estimation of the rate of evaporation. Finally, due to the changes in the droplet size, vapor diffusion from the droplet also induces errors in the classical transient based method. However, for droplet evaporating at steady state; vapor effects are steady as well since no driving forces exist to induce changes on the vapor dynamics in the vicinity of the droplet. Hence, for applications

16 ACS Paragon Plus Environment

Page 16 of 36

Page 17 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

requiring steadiness of evaporating droplets, the newly developed steady method corrects the shortcomings of the traditional unsteady methods. The steady method developed here presents a unique platform to study the evaporation dynamics of droplets. In addition to higher accuracy, the approach allows for the study of a multitude of working fluids at length scales not previously possible using transient based methods. For example, low surface tension fluids such as ethanol, methanol, and toluene are notoriously difficult to dispense in small quantities on omniphobic surfaces. The use of a piezoelectric dispenser eliminates this problem by quantizing the dispensed fluid to a small volume and ‘building’ it on the surface independently of the dispenser. In the future, it would be interesting to use the steady approach to build an evaporation rate database for a multitude of fluids and conditions. Furthermore, the inexorable presence of contact angle hysteresis on any real surface introduces a range of contact angles a droplet can exist in for a single droplet base radius in the pinned state. The steady approach overcomes this limitation by enabling the study of droplet evaporation on a surface for the entire range of contact angles between the receding and advancing state without the need for crude approximations. Future studies should use the data presented here as a benchmark for 3D numerical simulations of droplet evaporation,

which classically use non-moving boundaries and the

quasi-steady

approximation71-72. Furthermore, the steady approach developed here has a plethora of other applications where the volume flux of a working fluid needs to be measured, whether it be permeation through a membrane73, or the Wickability of a superhydrophilic surface74-75. Indeed, the use of the steady approach to characterize Wickability of surfaces would eliminate the need for contact based methods which may introduce errors in the measurement. 4. Conclusions In this study, a steady evaporation method is developed for surfaces of different functionalizations (i.e. hydrophilic or hydrophobic). By utilizing a piezoelectric dispenser to 17 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

feed microscale droplets ( ≈ 9 µm) to a larger evaporating droplet at a prescribed frequency, we studied water evaporation of droplets having base radii ranging from 20 µm to 250 µm on surfaces of different functionalities (45˚ ≤ , ≤ 162˚, where , is the apparent advancing contact angle). Utilizing our steady method we showed that the rate of evaporation has a linear dependence on the droplet base radius; and the evaporative flux has a linear dependence on contact angle for both hydrophobic and hydrophilic samples. We benchmarked our technique with the classical unsteady method showing significant improvement (≈ 140%) in evaporation rate measurement accuracy for droplets evaporating in steady state where the contact line and interfacial effects are negligible. Our work not only sheds light into the evaporation dynamics for various droplet sizes and contact angles on functional surfaces, it provides an experimental platform that enables the decoupling of the intricate and length-scale-dependent balance played by internal and external flows, evaporation kinetics, thermocapillary motion, and binary-mixture dynamics. ASSOCIATED CONTENT Supporting Information The supporting information is available free of charge on the ACS Publications website at DOI: Logarithmic plots of the evaporative flux, and a table of values of the parameters used in the calculation of non-dimensional numbers. Author Contributions A.A.G. and N. M. conceived the initial idea for this study. N. M. guided the work. A.A.G and S.S. fabricated, functionalized and characterized the samples and carried out the experiments. J.O. carried out SEM imaging. A.A.G, S.S. and J.O analyzed the data. All authors contributed to the writing of the manuscript and have given approval to the final version of this manuscript. 18 ACS Paragon Plus Environment

Page 18 of 36

Page 19 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Notes The authors declare no competing financial interest.

Acknowledgments The authors would like to thank Xiao Yan for his support in image analysis. The authors gratefully acknowledge funding support from the National Science Foundation under Award No. 1554249. N.M. gratefully acknowledges funding support from the International Institute for Carbon Neutral Energy Research (WPI-I2CNER), sponsored by the Japanese Ministry of Education, Culture, Sports, Science and Technology. SEM imaging was carried out at the Beckman Institute for Advanced Science and Technology, University of Illinois.

19 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References 1. Ramanathan, V.; Crutzen, P.; Kiehl, J.; Rosenfeld, D., Aerosols, climate, and the hydrological cycle. science 2001, 294 (5549), 2119-2124. 2. Worden, J.; Noone, D.; Bowman, K.; Beer, R.; Eldering, A.; Fisher, B.; Gunson, M.; Goldman, A.; Herman, R.; Kulawik, S. S., Importance of rain evaporation and continental convection in the tropical water cycle. Nature 2007, 445 (7127), 528-532. 3. Brutsaert, W.; Parlange, M., Hydrologic cycle explains the evaporation paradox. Nature 1998, 396 (6706), 30-30. 4. Black, T. A.; Gardner, W. R.; Thurtell, G. W., The Prediction of Evaporation, Drainage, and Soil Water Storage for a Bare Soil1. Soil Science Society of America Journal 1969, 33 (5), 655-660. 5. Fakhru’l-Razi, A.; Pendashteh, A.; Abdullah, L. C.; Biak, D. R. A.; Madaeni, S. S.; Abidin, Z. Z., Review of technologies for oil and gas produced water treatment. Journal of hazardous materials 2009, 170 (2), 530-551. 6. Shang, Y.; Si, Y.; Raza, A.; Yang, L.; Mao, X.; Ding, B.; Yu, J., An in situ polymerization approach for the synthesis of superhydrophobic and superoleophilic nanofibrous membranes for oil–water separation. Nanoscale 2012, 4 (24), 7847-7854. 7. Alex, R.; Bodmeier, R., Encapsulation of water-soluble drugs by a modified solvent evaporation method. I. Effect of process and formulation variables on drug entrapment. Journal of Microencapsulation 1990, 7 (3), 347-355. 8. Nihant, N.; Schugens, C.; Grandfils, C.; Jérôme, R.; Teyssié, P., Polylactide microparticles prepared by double emulsion/evaporation technique. I. Effect of primary emulsion stability. Pharmaceutical Research 1994, 11 (10), 1479-1484. 9. Abramzon, B.; Sirignano, W. A., Droplet vaporization model for spray combustion calculations. International Journal of Heat and Mass Transfer 1989, 32 (9), 1605-1618. 10. Crespo, A.; LiÑAn, A., Unsteady Effects in Droplet Evaporation and Combustion. Combustion Science and Technology 1975, 11 (1-2), 9-18. 11. Sazhin, S.; Abdelghaffar, W.; Sazhina, E.; Heikal, M., Models for droplet transient heating: effects on droplet evaporation, ignition, and break-up. International journal of thermal sciences 2005, 44 (7), 610-622. 12. Kitano, T.; Nishio, J.; Kurose, R.; Komori, S., Evaporation and combustion of multicomponent fuel droplets. Fuel 2014, 136, 219-225. 13. Boreyko, J. B.; Collier, C. P., Delayed Frost Growth on Jumping-Drop Superhydrophobic Surfaces. ACS Nano 2013, 7 (2), 1618-1627. 14. Duan, F.; Ward, C. A., Investigation of Local Evaporation Flux and Vapor-Phase Pressure at an Evaporative Droplet Interface. Langmuir 2009, 25 (13), 7424-7431. 15. Hu, H.; Larson, R. G., Evaporation of a Sessile Droplet on a Substrate. The Journal of Physical Chemistry B 2002, 106 (6), 1334-1344. 16. Pan, Z.; Weibel, J. A.; Garimella, S. V., Influence of Surface Wettability on Transport Mechanisms Governing Water Droplet Evaporation. Langmuir 2014, 30 (32), 9726-9730. 17. Dash, S.; Garimella, S. V., Droplet Evaporation Dynamics on a Superhydrophobic Surface with Negligible Hysteresis. Langmuir 2013, 29 (34), 10785-10795. 18. Shahidzadeh-Bonn, N.; RafaÏ, S.; Azouni, A. Z. A.; Bonn, D., Evaporating droplets. Journal of Fluid Mechanics 2006, 549, 307-313. 19. He, M.; Zhang, X.; Zeng, K.; Gao, K., A combined thermodynamic cycle used for waste heat recovery of internal combustion engine. Energy 2011, 36 (12), 6821-6829. 20. Bombarda, P.; Invernizzi, C. M.; Pietra, C., Heat recovery from Diesel engines: A thermodynamic comparison between Kalina and ORC cycles. Applied Thermal Engineering 2010, 30 (2), 212-219. 20 ACS Paragon Plus Environment

Page 20 of 36

Page 21 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

21. He, C.; Liu, C.; Gao, H.; Xie, H.; Li, Y.; Wu, S.; Xu, J., The optimal evaporation temperature and working fluids for subcritical organic Rankine cycle. Energy 2012, 38 (1), 136-143. 22. Churi, N.; Achenie, L. E., The optimal design of refrigerant mixtures for a twoevaporator refrigeration system. Computers & chemical engineering 1997, 21, S349-S354. 23. Penman, H. L. In Natural evaporation from open water, bare soil and grass, Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, The Royal Society: 1948; pp 120-145. 24. Priestley, C.; Taylor, R., On the assessment of surface heat flux and evaporation using large-scale parameters. Monthly weather review 1972, 100 (2), 81-92. 25. Monteith, J. L. In Evaporation and environment, Symp. Soc. Exp. Biol, 1965; p 4. 26. Brutsaert, W., Evaporation into the atmosphere: theory, history and applications. Springer Science & Business Media: 2013; Vol. 1. 27. Hubbard, G.; Denny, V.; Mills, A., Droplet evaporation: effects of transients and variable properties. International Journal of Heat and Mass Transfer 1975, 18 (9), 10031008. 28. Plawsky, J.; Fedorov, A.; Garimella, S.; Ma, H.; Maroo, S.; Chen, L.; Nam, Y., Nanoand microstructures for thin-film evaporation—a review. Nanosc Microsc Therm 2014, 18 (3), 251-269. 29. Birdi, K. S.; Vu, D. T.; Winter, A., A study of the evaporation rates of small water drops placed on a solid surface. The Journal of Physical Chemistry 1989, 93 (9), 3702-3703. 30. McHale, G.; Aqil, S.; Shirtcliffe, N. J.; Newton, M. I.; Erbil, H. Y., Analysis of Droplet Evaporation on a Superhydrophobic Surface. Langmuir 2005, 21 (24), 11053-11060. 31. Ranjan, R.; Murthy, J. Y.; Garimella, S. V., Analysis of the Wicking and Thin-Film Evaporation Characteristics of Microstructures. Journal of Heat Transfer 2009, 131 (10), 101001-101001-11. 32. Nešić, S.; Vodnik, J., Kinetics of droplet evaporation. Chemical Engineering Science 1991, 46 (2), 527-537. 33. Picknett, R. G.; Bexon, R., The evaporation of sessile or pendant drops in still air. Journal of Colloid and Interface Science 1977, 61 (2), 336-350. 34. Bradley, R. S.; Evans, M. G.; Whytlaw-Gray, R. W., The Rate of Evaporation of Droplets. Evaporation and Diffusion Coefficients, and Vapour Pressures of Dibutyl Phthalate and Butyl Stearate. Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 1946, 186 (1006), 368. 35. Crafton, E. F.; Black, W. Z., Heat transfer and evaporation rates of small liquid droplets on heated horizontal surfaces. International Journal of Heat and Mass Transfer 2004, 47 (6–7), 1187-1200. 36. Hubbard, G. L.; Denny, V. E.; Mills, A. F., Droplet evaporation: Effects of transients and variable properties. International Journal of Heat and Mass Transfer 1975, 18 (9), 10031008. 37. Kincaid, D.; Longley, T., A water droplet evaporation and temperature model. Trans. ASAE 1989, 32 (2), 457-463. 38. Gelderblom, H.; Marín, Á. G.; Nair, H.; van Houselt, A.; Lefferts, L.; Snoeijer, J. H.; Lohse, D., How water droplets evaporate on a superhydrophobic substrate. Physical Review E 2011, 83 (2), 026306. 39. Deegan, R. D.; Bakajin, O.; Dupont, T. F.; Huber, G.; Nagel, S. R.; Witten, T. A., Contact line deposits in an evaporating drop. Physical Review E 2000, 62 (1), 756-765. 40. Deegan, R. D.; Bakajin, O.; Dupont, T. F.; Huber, G., Capillary flow as the cause of ring stains from dried liquid drops. Nature 1997, 389 (6653), 827.

21 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

41. Deegan, R. D.; Bakajin, O.; Dupont, T. F.; Huber, G.; Nagel, S. R.; Witten, T. A., Contact line deposits in an evaporating drop. Physical review E 2000, 62 (1), 756. 42. Hu, H.; Larson, R. G., Analysis of the microfluid flow in an evaporating sessile droplet. Langmuir 2005, 21 (9), 3963-3971. 43. Hu, H.; Larson, R. G., Analysis of the effects of Marangoni stresses on the microflow in an evaporating sessile droplet. Langmuir 2005, 21 (9), 3972-3980. 44. Hu, H.; Larson, R. G., Marangoni effect reverses coffee-ring depositions. The Journal of Physical Chemistry B 2006, 110 (14), 7090-7094. 45. Picknett, R.; Bexon, R., The evaporation of sessile or pendant drops in still air. Journal of Colloid and Interface Science 1977, 61 (2), 336-350. 46. Popov, Y. O., Evaporative deposition patterns: spatial dimensions of the deposit. Physical Review E 2005, 71 (3), 036313. 47. Sefiane, K.; Moffat, J.; Matar, O.; Craster, R., Self-excited hydrothermal waves in evaporating sessile drops. Applied Physics Letters 2008, 93 (7), 074103. 48. Sobac, B.; Brutin, D., Thermocapillary instabilities in an evaporating drop deposited onto a heated substrate. Phys Fluids 2012, 24 (3), 032103. 49. Christy, J. R.; Hamamoto, Y.; Sefiane, K., Flow transition within an evaporating binary mixture sessile drop. Physical review letters 2011, 106 (20), 205701. 50. Kim, H.; Boulogne, F.; Um, E.; Jacobi, I.; Button, E.; Stone, H. A., Controlled uniform coating from the interplay of Marangoni flows and surface-adsorbed macromolecules. Physical review letters 2016, 116 (12), 124501. 51. Tan, H.; Diddens, C.; Lv, P.; Kuerten, J. G.; Zhang, X.; Lohse, D., Evaporationtriggered microdroplet nucleation and the four life phases of an evaporating Ouzo drop. Proceedings of the National Academy of Sciences 2016, 113 (31), 8642-8647. 52. Sáenz, P.; Wray, A.; Che, Z.; Matar, O.; Valluri, P.; Kim, J.; Sefiane, K., Dynamics and universal scaling law in geometrically-controlled sessile drop evaporation. Nat Commun 2017, 8. 53. Enright, R.; Miljkovic, N.; Sprittles, J.; Nolan, K.; Mitchell, R.; Wang, E. N., How coalescing droplets jump. ACS nano 2014, 8 (10), 10352-10362. 54. Cha, H.; Xu, C.; Sotelo, J.; Chun, J. M.; Yokoyama, Y.; Enright, R.; Miljkovic, N., Coalescence-induced nanodroplet jumping. Physical Review Fluids 2016, 1 (6), 064102. 55. Bormashenko, E.; Pogreb, R.; Whyman, G.; Erlich, M., Resonance Cassie−Wenzel Wetting Transition for Horizontally Vibrated Drops Deposited on a Rough Surface. Langmuir 2007, 23 (24), 12217-12221. 56. Noblin, X.; Buguin, A.; Brochard-Wyart, F., Vibrated sessile drops: Transition between pinned and mobile contact line oscillations. The European Physical Journal E: Soft Matter and Biological Physics 2004, 14 (4), 395-404. 57. Celestini, F.; Kofman, R., Vibration of submillimeter-size supported droplets. Physical Review E 2006, 73 (4), 041602. 58. Aoki, K.; Takata, S.; Kosuge, S., Vapor flows caused by evaporation and condensation on two parallel plane surfaces: Effect of the presence of a noncondensable gas. Phys Fluids 1998, 10 (6), 1519-1533. 59. Murase, M.; Kataoka, Y.; Fujii, T., Evaporation and condensation heat transfer with a noncondensable gas present. Nuclear engineering and design 1993, 141 (1-2), 135-143. 60. Stauber, J. M.; Wilson, S. K.; Duffy, B. R.; Sefiane, K., Evaporation of droplets on strongly hydrophobic substrates. Langmuir 2015, 31 (12), 3653-3660. 61. Weisensee, P. B.; Wang, Y.; Hongliang, Q.; Schultz, D.; King, W. P.; Miljkovic, N., Condensate droplet size distribution on lubricant-infused surfaces. International Journal of Heat and Mass Transfer 2017, 109, 187-199.

22 ACS Paragon Plus Environment

Page 22 of 36

Page 23 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

62. Kim, S.; Kim, K. J., Dropwise condensation modeling suitable for superhydrophobic surfaces. Journal of heat transfer 2011, 133 (8), 081502. 63. Miljkovic, N.; Enright, R.; Wang, E. N., Modeling and optimization of superhydrophobic condensation. Journal of Heat Transfer 2013, 135 (11), 111004. 64. Kandlikar, S. G.; Steinke, M. E., Contact angles and interface behavior during rapid evaporation of liquid on a heated surface. International Journal of Heat and Mass Transfer 2002, 45 (18), 3771-3780. 65. Armstrong, R.; O’Rourke, P.; Zardecki, A., Vaporization of irradiated droplets. The Physics of fluids 1986, 29 (11), 3573-3581. 66. Beysens, D.; Steyer, A.; Guenoun, P.; Fritter, D.; Knobler, C. M., How Does Dew Form? Phase Transitions 1991, 31 (1-4), 219-246. 67. Bergman, T. L.; Incropera, F. P., Fundamentals of heat and mass transfer. John Wiley & Sons: 2011. 68. Koschmieder, E. L., Bénard cells and Taylor vortices. Cambridge University Press: 1993. 69. Kays, W. M., Convective heat and mass transfer. Tata McGraw-Hill Education: 2012. 70. Xu, X.; Ma, L., Analysis of the effects of evaporative cooling on the evaporation of liquid droplets using a combined field approach. Sci Rep-Uk 2015, 5. 71. Miller, R.; Harstad, K.; Bellan, J., Evaluation of equilibrium and non-equilibrium evaporation models for many-droplet gas-liquid flow simulations. International Journal of Multiphase Flow 1998, 24 (6), 1025-1055. 72. Aggarwal, S.; Peng, F., A review of droplet dynamics and vaporization modeling for engineering calculations. TRANSACTIONS-AMERICAN SOCIETY OF MECHANICAL ENGINEERS JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER 1995, 117, 453-453. 73. Bouwmeester, H.; Kruidhof, H.; Burggraaf, A., Importance of the surface exchange kinetics as rate limiting step in oxygen permeation through mixed-conducting oxides. Solid State Ionics 1994, 72, 185-194. 74. Rahman, M. M.; Olceroglu, E.; McCarthy, M., Role of wickability on the critical heat flux of structured superhydrophilic surfaces. Langmuir 2014, 30 (37), 11225-11234. 75. Son, H. H.; Seo, G. H.; Jeong, U.; Shin, D. Y.; Kim, S. J., Capillary wicking effect of a Cr-sputtered superhydrophilic surface on enhancement of pool boiling critical heat flux. International Journal of Heat and Mass Transfer 2017, 113, 115-128.

23 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figures

Figure 1. Scanning electron microscopy images of (a) gold coated silicon wafer (Si-Au) (de = 45.1° ± 5.2°), (b) polished copper (Cu) (de = 85˚ ± 1.9˚), (c) TFTS coated silicon wafer (SiTFTS) (de = 110.6˚ ± 2.8˚), (d) 90 nm P2i Carbon Nanotubes (CNT1) (de,eff = 130.6˚ ± 2.9˚), (e) 30 nm P2i Carbon Nanotubes (CNT2) (de,eff = 141.3° ± 3.1°), and (f) HDTS Coated Copper Oxide (CuO-HDTS) (de,eff = 162° ± 2.9°).

24 ACS Paragon Plus Environment

Page 24 of 36

Page 25 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 2. (a) Schematic of the experimental setup. The substrate is placed on top of the stage at ambient temperature, and droplets grown on the substrate. The piezoelectric dispenser is placed a distance g above the substrate, and forms small microdroplets (≈ 9 µm diameter) to feed the bigger droplet growing on the substrate. The droplet growth process is monitored by a high speed camera from the side. Images from the high speed camera of the (b) strobed dispensing microdroplets, and the growing droplet on the substrate (c) right before and (d) right after incoming droplet coalescence. The small droplet on top of the large droplet residing on the surface in figure (c) is the droplet dispensed by the piezoelectric dispenser.

25 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Side view time lapse images of droplet growth for frequency of microdroplet deposition of (a) h = 5 Hz, h = 10 Hz, and h = 15 Hz. (b) Droplet radius as a function of time for different frequency of deposition. It is seen that the droplet asymptotes a constant size as the deposition time increases due to the matching of rate of evaporation and rate of deposition at a certain droplet radius.

26 ACS Paragon Plus Environment

Page 26 of 36

Page 27 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 4. Rate of evaporation as a function of the (a) droplet radius and (b) droplet base radius for all tested samples. The rate of evaporation varies linearly with increasing droplet size due to the linear dependence of contact line length with base radius.

27 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. Droplet evaporative flux as a function of the (a) droplet base radius (ij ) and (b) droplet advancing contact angle (de). A cut-off exists at de = kl° due to the change in the cause of evaporation (for a hydrophilic sample, the heat transfer from the base governs evaporation; whereas evaporation is governed by vapor diffusion for a hydrophobic sample). For log-log plots of (a) with associated slopes, please see supplementary information, Figure S1.

28 ACS Paragon Plus Environment

Page 28 of 36

Page 29 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 6. Droplet rate of evaporation characterized by the traditional transient method (solid symbols) and the steady method (hollow symbols). For all surfaces tested, the traditional method overestimates the rate of evaporation by approximately 140% due to the inner transient effects and the liquid-gas momentum exchange at the interface. The main sources of error included measurement error of droplet size for the traditional method, and error due to spread in the dispensed droplet size (≈ ±6%) for the steady method.

29 ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ToC Image

30 ACS Paragon Plus Environment

Page 30 of 36

Page 31 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 1. Scanning electron microscopy images of (a) gold coated silicon wafer (Si-Au) (θ = 45.1° ± 5.2°), (b) polished copper (Cu) (θ = 85˚ ± 1.9˚), (c) TFTS coated silicon wafer (Si-TFTS) (θ = 110.6˚ ± 2.8˚), (d) 90 nm P2i Carbon Nanotubes (CNT1) (θ = 130.6˚ ± 2.9˚), (e) 30 nm P2i Carbon Nanotubes (CNT2) (θ = 141.3° ± 3.1°), and (f) HDTS Coated Copper Oxide (CuO-HDTS) (θ = 162° ± 2.9°). 437x192mm (96 x 96 DPI)

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. (a) Schematic of the experimental setup. The substrate is placed on top of the stage at ambient temperature, and droplets grown on the substrate. The piezoelectric dispenser is placed a distance h above the substrate, and forms small microdroplets (≈ 9 µm diameter) to feed the bigger droplet growing on the substrate. The droplet growth process is monitored by a high speed camera from the side. Images from the high speed camera of the (b) strobed dispensing microdroplets, and the growing droplet on the substrate (c) right before and (d) right after incoming droplet coalescence. The small droplet on top of the large droplet residing on the surface in figure (c) is the droplet dispensed by the piezoelectric dispenser. 155x115mm (120 x 120 DPI)

ACS Paragon Plus Environment

Page 32 of 36

Page 33 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 3. Side view time lapse images of droplet growth for frequency of microdroplet deposition of (a) f = 5 Hz, f = 10 Hz, and f = 15 Hz. (b) Droplet radius as a function of time for different frequency of deposition. It is seen that the droplet asymptotes a constant size as the deposition time increases due to the matching of rate of evaporation and rate of deposition at a certain droplet radius. 180x227mm (120 x 120 DPI)

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Rate of evaporation as a function of the (a) droplet radius and (b) droplet base radius for all tested samples. The rate of evaporation varies linearly with increasing droplet size due to the linear dependence of contact line length with base radius. 274x106mm (120 x 120 DPI)

ACS Paragon Plus Environment

Page 34 of 36

Page 35 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 5. Droplet evaporative flux as a function of the (a) droplet base radius (R_b) and (b) droplet advancing contact angle (θ_a). A cut-off exists at θ_a=90° due to the change in the cause of evaporation (for a hydrophilic sample, the heat transfer from the base governs evaporation; whereas evaporation is governed by vapor diffusion for a hydrophobic sample). For log-log plots of (a) with associated slopes, please see supplementary information, Figure S1. 260x107mm (120 x 120 DPI)

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Droplet rate of evaporation characterized by the traditional transient method (solid symbols) and the steady method (hollow symbols). For all surfaces tested, the traditional method overestimates the rate of evaporation by approximately 140% due to the inner transient effects and the liquid-gas momentum exchange at the interface. The main sources of error included measurement error of droplet size for the traditional method, and error due to spread in the dispensed droplet size (≈ ±6%) for the steady method. 101x70mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 36 of 36