Unexpected Suppression of Leidenfrost Phenomenon on Superhydrophobic Surfaces
1 Unexpected Suppression of Leidenfrost Phenomenon on Superhydrophobic Surfaces
Meng Shi, Ratul Das, Sankara Arunachalam, Himanshu Mishra* Interfacial Lab (iLab), Water Desalination and Reuse Center (WDRC), Biological and Environmental Science and Engineering (BESE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia *Corresponding author: [email protected] Abstract
The Leidenfrost phenomenon entails the levitation of a liquid droplet over a superheated surface, cushioned by its vapor layer. For water, superhydrophobic surfaces are believed to suppress the Leidenfrost point ( T L )—the temperature at which this phenomenon occurs. The vapor film obstructs boiling heat transfer in heat exchangers, thereby compromising energy efficiency and safety. Thus, it is desirable to realize superhydrophobicity without suppressing T L . Here we demonstrate that the T L of water on microtextured superhydrophobic surfaces comprising doubly reentrant pillars (DRPs) can exceed those on hydrophilic and even superhydrophilic surfaces. We disentangle the contributions of microtexture, heat transfer, and surface chemistry on T L and reveal how superhydrophobicity can be realized without suppressing T L . For instance, silica surfaces with DRPs facilitate ~300% greater heat transfer to water droplets at 200 ° C in comparison with silica surfaces coated with perfluorinated-nanoparticles. Thus, superhydrophobic surfaces could be harnessed for energy efficient thermal machinery.
Keywords:
Leidenfrost phenomenon; superhydrophobic surfaces; doubly reentrant pillars; water interfaces; adhesion; heat transfer; microtextured surfaces Introduction
Known for over two centuries, the Leidenfrost phenomenon entails the levitation of an evaporating liquid droplet on a hot surface such that the droplet’s weight is counterbalanced by the pressure of the vapor film formed underneath . This phenomenon has attracted interest across many disciplines
2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 , and it has been applied to engineering applications, such as directional droplet propulsion
7, 9 , electricity generation , and green chemistry , among others . For liquids with strong intermolecular forces and a high boiling point such as water, the Leidenfrost phenomenon occurs—the Leidenfrost point ( T L )—at temperatures significantly higher than at normal temperature and pressure (NTP: 293 K and 1 atm). For example, for a smooth and flat silica surface, the T L of water is in the range 270 ° C–300 ° C
16, 17 . As a corollary, the Leidenfrost phenomenon can be realized at NTP with liquids that have ultralow boiling points, such as liquid nitrogen and 1,1,1,2-tetrafluoroethane with boiling points of −196 ° C and −26.1 ° C , respectively. In either scenario, a vapor film with a thickness of ~100 µ m levitates the droplet above the surface, after which the droplet slowly evaporates, akin to film boiling
3, 9, 18 . Roughening of surfaces can have profound consequences on T L : for example, hierarchical micro/nano-texturing furthers the wettability of a silica surface with water, increasing the temperature needed for the water droplet to achieve the Leidenfrost state (e.g., T L > 400 ° C)
16, 17, 19, 20 . Conversely, when hydrophobic coatings are applied to micro/nano-textured surfaces, they can give rise to superhydrophobicity, characterized by water contact angles, q r ≥ Dq ≤
21, 22, 23 . Superhydrophobicity lowers the T L of water by facilitating (i) low adhesion energy per unit area at the solid–liquid–vapor (SLV) interface given by the Young-Dupré equation, 𝑊 !" =𝛾 × (1 + cos𝜃 % ) , where g LV is the surface tension, and 𝜃 % > 90° is the actual (or intrinsic) contact angle of the liquid (water) , and (ii) reduced liquid–solid contact area (hence adhesion) due to the interfacial entrapment of air
4, 25, 26, 27, 28 . For instance, after a flat metallic surface was TM coating (hereafter referred to as Glaco), the T L of water reduced from 210°C to 130°C; similarly, Glaco-coated superhydrophobic steel balls and rotors exhibited T L approaching 100°C and also afforded significant reduction in frictional drag . Taken together, these findings give credence to conventional wisdom that superhydrophobicity lowers the T L of water
25, 27, 28, 29 . Superhydrophobic surfaces have been exploited to enhance pool boiling ; however, their applications in thermal machinery operating at high temperatures remain limited as superhydrophobicity generally lowers the T L of water. Specifically, the vapor film prevents the interfacial heat transfer because of its lower thermal conductivity than that of the liquid , for example, in boiler pipes, heat exchangers, and turbines. This situation, known as boiling crisis or burnout , can severely impact energy efficiency, damage machinery, and cause accidents . Hence, it is of practical interest to investigate whether the lowering of the T L of water is a universal attribute of superhydrophobicity, similarly to the water contact angle, q r > 150°, the contact angle hysteresis, Dq < 20°, and the robust entrapment of air under water on superhydrophobic surfaces, because superhydrophobicity without suppressing the T L of water could reduce frictional drag while facilitating interfacial heat transfer. Here, contrary to conventional wisdom, we demonstrate that superhydrophobic surfaces do not always lower the T L of water. We utilize silica surfaces adorned with a specialized biomimetic microtexture comprising doubly reentrant pillars (DRPs) that exhibit superhydrophobicity despite their hydrophilic surface chemistry.
22, 23, 35, 36
Counterintuitively, water’s T L on these superhydrophobic silica surfaces is significantly higher than that on flat hydrophilic silica surfaces and even higher than that on superhydrophilic microtextured silica surfaces (see Table 1 for the apparent contact angles of water droplets). We combine experiment Results
We utilized SiO /Si wafers and Glaco-coated SiO /Si wafers as hydrophilic and superhydrophobic surfaces (hereafter referred to as flat silica and Glaco-coated silica, respectively) (Figs. 1a and 1b) . Additionally, we microfabricated arrays of cylindrical pillars ( Fig. 1c ) and DRPs with diameter, D = 10 µ m, height, H = µ m, and pitch, L = 50 µ m ( Fig. 1d-f) . The details of the microfabrication processes have been reported recently; briefly, silicon wafers with a 2.4- µ m-thick silica layer were prepared by photolithography, dry etching, and thermal oxide growth
22, 23, 37, 38 ( Methods ). Fig. 1| Representative scanning electron micrographs of surfaces investigated in this study and their wettability with water at NTP. a) Flat silica; b) Glaco-coated silica that exhibits superhydrophobicity; c) tilted view (45 ° ) of cylindrical pillars with a diameter ( D ) of 10 µ m, height ( H ) of 50 µ m, and pitch ( L ) of 50 µ m; d) tilted view (45 ° ) of doubly reentrant pillars with D = 10 µ m, H = 50 µ m, and L = 50 µ m; e, f) magnified views of DRPs with the vertical length and width of the double reentrant edge, respectively, S v = 2 µ m and S h = 0.5 µ m. g–j) Representative images of water droplets (15 µ L) on the surfaces outlined in panels (a–d) and their apparent contact angles (CAs). First, we characterized the wettability of these surfaces with water by measuring the advancing ( 𝜃 & ) and receding ( 𝜃 ’ ) contact angles (CAs) of droplets (6 µ L) advanced/retracted at 0.2 µ L/s (Methods) . On flat silica, 𝜃 & = 45 ° ± 2 ° , 𝜃 ’ = 0°, and 𝜃 ( = 40 ° ± 2 ° ( Fig. 1g ); therefore we considered 40 ° to be the actual (or intrinsic) contact angle ( 𝜃 % ) for silica material in our theoretical analysis . Glaco-coated silica exhibited superhydrophobicity with q r = 156° ±
1° and contact angle hysteresis Dq = 𝜃 & − 𝜃 ’ =19° ±
2° (
Fig. 1h ). Conversely, silica surfaces with cylindrical pillars were superhydrophilic—they imbibed water, thus reaching the fully filled (or Wenzel) state ( Fig. 1i ); apparent CAs of water droplets were ultralow, q r < 10°. Silica surfaces with DRPs exhibited superhydrophobicity with q r = 150° ±
1° and Dq =19° ±
2° (
Fig. 1j ), similarly to Glaco-coated silica. The sliding angles of water droplets (10 µ L) on Glaco-coated silica and silica with DRPs were 7 ° ± 3 ° and 13 ° ± 2 ° , respectively ( Table 1 and Methods ). Table 1. Characterization of various surfaces investigated in this work
Substrates Apparent contact angles Contact angle hysteresis Nature of the substrate Sliding angle Measured Leidenfrost point, T L ( ° C ) 𝜽 𝐑 𝜽 𝐫 𝜽 𝐀 ∆𝜽 T L Flat silica ° ° ± 2 ° ° ± 2 ° ° ± 2 ° Hydrophilic / ~286
Glaco-coated silica ° ° ± 1 ° ° ± 1 ° ° ± 2 ° Superhydrophobic ° ± 3 ° ~130 Cylindrical pillars ° <10 ° <10 ° <10 ° Superhydrophilic / ~314
DRPs ° ± 2 ° ° ± 1 ° ° ± 2 ° ° ± 2 ° Superhydrophobic ° ± 2 ° ~363 ± ° C): (i) flat silica, (ii) Glaco-coated silica, and silica surfaces with arrays of (iii) cylindrical pillars and (iv) DRPs. We arbitrarily chose an initial temperature of 347°C, which is beyond the T L of water on common surfaces . Water droplets (15 µ L) were gently placed onto the surfaces, and the subsequent behaviors were recorded via high-speed imaging. The placement of droplets ensured that the Weber numbers, defined as the ratio of the droplet’s inertia to capillarity,
𝑊𝑒 = 𝜌𝑉 , 𝑅 " 𝛾 ⁄ , where r , R d , g lv , and V is the density, radius, water surface tension, and impact velocity, respectively, were as small as possible (~1.7). Contrary to our expectation, all samples except superhydrophobic DRPs exhibited the Leidenfrost phenomenon at 347°C (
Supplementary Materials, Movies S1–S4,
Figs. 2a–d ). Specifically, after the initial impact and followed by 3-6 bounces due to inertia, water droplets re-established contact with the DRPs with slight vibrations (
Supplementary Materials, Movie S4 ). To confirm this result, we ascertained T L on these surfaces by tracking the lifetime of sessile water droplets on them as a function of surface temperature to determine how long it took for the water droplets to evaporate on these surfaces. As the surface temperature ( T s ( ° C)) increases, the droplet lifetime decreases due to faster evaporation. In the vicinity of T L , the droplet lifetime starts increasing monotonically due to the formation of an insulating vapor layer underneath the droplet that diminishes the interfacial heat transfer. By gradually increasing the temperature in the range 70°C–440 ° C, the T L values for hydrophilic silica ° C and ~314 ° C, respectively, (
Fig. 2e ). Fig. 2| High-speed micrographs of water droplets exhibiting Leidenfrost phenomena after gentle placement on the following surfaces maintained at 347°C : a) flat silica (hydrophilic); b) Glaco-coated silica (superhydrophobic); c) silica surface with cylindrical pillars with D = 10 µ m, H = 50 µ m, and L = 50 µ m (superhydrophilic); and d) silica surfaces with doubly reentrant pillars (DRPs) with D = 10 µ m, H = 50 µ m, and L = 50 µ m (superhydrophobic). These snapshots are derived from Supplementary Movies S1–S4 . Lifetimes of water droplets placed on the surfaces described in (a–d) as a function of surface temperature: e) hydrophilic flat silica surface and superhydrophilic silica surfaces with cylindrical pillars and f) superhydrophobic surfaces (Glaco-coated silica and silica with DRPs).
For superhydrophobic surfaces—DRPs and Glaco-coated silica—there were no distinct peaks in the droplet lifetime curves (Fig. 2f) . Therefore, we utilized high-speed imaging to monitor vapor film formation at the droplet-solid interface to pinpoint the onset of Leidenfrost phenomena (
Supplementary Materials, Figs. S3 and S4 ). For water droplets, superhydrophobic Glaco-coated silica exhibited T L ~130 ° C. In stark contrast, however, superhydrophobic silica surfaces with DRPs did not exhibit the Leidenfrost phenomenon until ° C. This dramatic variation in the T L values of superhydrophobic surfaces challenges our current understanding of the factors and mechanisms underlying the Leidenfrost phenomenon. Fig. 3| Pinning of water droplets on silica surfaces adorned with DRPs maintained at 347°C . a) The water droplet is unable to detach completely from the DRPs, and it exhibits an oscillatory motion during which it contacts a varying number of pillars; b, c) droplet at its furthest point from the surface, yet pinned on a few DRPs. Schematics of the forces at the liquid–solid–vapor interface due to adhesion ( F Ad ), weight ( G ), and vapor force ( F v ) as the drop moves d) upward and e) downward. f–i) High-speed snapshots of the curvature of the droplet–DRPs interface moving f–h) upward and g–i) downward. Next, we observed the behavior of water droplets (15 µ L) on silica with DRPs at T = 347°C via high-speed imaging complemented with high-magnification optics ( Fig. 3 ). We found that a typical droplet exhibited an up-and-down vibratory motion under the influence of its vapor force built underneath, weight, and adhesion at the SLV interface (
Figs. 3a–3c , and
Supplementary Materials, Video S5 ). As the droplet moved upward, the number of DRPs in contact with the droplet decreased and vice versa, and the air–water interface assumed concave and convex curvatures, respectively (
Figs. 3d and 3e ). After the initial contact, when the center
10 of mass of the droplet was at the highest point, it was still pinned to ~3–10 DRPs, which preempted the onset of the Leidenfrost phenomenon (
Figs. 3b and 3c ). We considered that as a liquid droplet is placed on a hot substrate, its vapor flow is obstructed by the surface microtexture, which builds up an upward vapor force, F v (N). If F v is greater than the droplet’s weight, G (N), it tends to lift it, whereas the liquid–solid adhesive force, F Ad (N), impedes this motion ( Fig. 3d ). Conversely, if F v < G , the droplet tends to sink into the microtexture, which the adhesive force, F Ad , also counters ( Fig. 3e ). In the concave configuration, the vapor has adequate space to escape; thus, F v drops ( Figs. 3f and ). Subsequently, the droplet advances downward into the microtexture, assuming a near-flat convex shape ( Figs. 3g and ) with F Ad acting upward, which indicates F v < G . To test the crucial role of F Ad , we reduced the liquid–solid adhesion on silica with DRPs by functionalizing it with perfluorodecyltrichlorosilane (FDTS) (Methods) , which yielded apparent CAs for water droplets in the range: 𝜃 & = 164 ± 2 ° and 𝜃 ’ = 143 ± 2 ° . When we placed water droplets onto FDTS-coated DRPs, T L dropped to ~172 ° C (
Supplementary Materials, Fig. S5a, and Video S6 ). The adhesion force between a DRP and water in our experimental configuration is given by 𝐹 &" = 𝛾 × 𝜋𝐷 × (1 + cos𝜃 % ) . As the surface chemistry changed from silica to FDTS, the intrinsic CA ( 𝜃 % ) of water (on smooth and flat surfaces) increased from 40 ° for silica to 112 ° for FDTS, it caused a 64% reduction in F Ad and lowered T L significantly. This point demonstrates the sensitivity of the Leidenfrost phenomenon on adhesion rather than empirical cut-offs set for the apparent CAs of water droplets on them . Particularly, not all superhydrophobic surfaces suppress the T L of water. But why did the onset of the Leidenfrost phenomenon on superhydrophobic silica with DRPs get delayed in comparison to that on superhydrophilic cylindrical pillars despite their identical
11 surface chemistry and similar dimensions? The lower contact angle hysteresis in the former due to the robust entrapment of air should have led to a lower T L than that of the latter. To address this question, first, we pinpoint some crucial differences in the ways DRPs interact with liquids in comparison to cylindrical pillars. The bioinspired DRP geometry is peculiar in the sense that despite the hydrophilicity of the substrate, it exhibits superhydrophobicity when liquids are placed on top because (i) the mushroom shape of the pillars prevents liquid imbibition so the droplet is in touch with the pillars’ caps and air, and (ii) as the droplet recedes, it has to detach from individual pillars, which requires small energy barriers to surpass, affording low contact angle hysteresis
22, 35, 42, 43 (e.g., please, see the confocal images in
Fig. S5b-c before and after FDTS coating on DRPs and their effect of L-S contact). In silica with DRPs, the pillar cap and base are comprised of silica (2 µ m thickness) and silicon (50 µ m), respectively. When the water droplet lands onto DRPs, it only touches the pillars for a depth of ~2 µ m ( Fig. 4a ), after which the air–water interface assumes a convex shape and any further liquid movement does not increase the liquid–solid contact (
Fig. 4b ). Particularly, the liquid–solid contact on the DRPs is sensitive only to the dimensions of the pillars’ caps and not to the penetration depth ( L p ); this peculiarity would cease to exist at ultrahigh Weber numbers, where the drops would penetrate the microtexture, and those cases are beyond the scope of this study. In stark contrast to DRPs, when a water droplet rests on superhydrophilic cylindrical pillars ( Supplementary Materials, Fig. S6 , and
Video S7) , it tends to penetrate them (
Fig. 4c ), thus increasing the liquid–solid contact area and deepening the liquid penetration depth (
Fig. 4d ). The enhanced penetration depth remarkably impacts heat transfer, which we next investigate computationally by simulating the temperature distributions
12 in DRPs and cylindrical pillars, followed by presenting an analytical model to derive expressions for F v to gain mechanistic insights into the onset of the Leidenfrost phenomenon. Computational and analytical modeling
Temperature distributions in silica surfaces adorned with DRPs and cylindrical pillars as a function of water penetration depth were simulated using COMSOL ® (COMSOL Multiphysics, version 5.4 ). In this approach, we assumed that (i) a quasi-steady state, (ii) the water temperature was 99°C , (iii) the convective heat transfer coefficient for the vapor inside microtextures was 20 W/(m × K) ( Fig. S7 ), and (iv) the evaporation at the air–water interface was considered as an outflow of heat flux. The experimentally observed penetration depths of water droplets into the microtextures comprising DRPs and cylindrical pillars were 2 and 20 µ m, respectively ( Figs. 4a and ). Consequently, the top region of the DRPs ( Fig. 4f ) were significantly hotter than that of the cylindrical pillars (
Fig. 4h ), which explains why the heat transfer in the case of the DRPs was significantly lower than that in the case of cylindrical pillars (
Fig. 4e and ). The higher heat transfer in the case of cylindrical pillars generated higher vapor fluxes and hence the upward vapor force, thereby lifting the droplet at lower temperatures than that in the case of DRPs. Fig. 4| Heat transfer analysis at the liquid–solid interface as a function of penetration depth.
Representative experimental micrograph of a water droplet penetrating silica surfaces adorned with (a) DRPs and (c) cylindrical pillars of identical diameter, height, and pitch. Schematics of water droplet penetration into DRPs (b) and cylindrical pillars (d), pinpointing heat transfer rates. Here, Q eva , is the overall evaporation heat; Q is the heat transfer rate in Region 1; and Q is the heat transfer rate in Region 2 (i.e., from the hot substrate). Temperature distributions due to heat transfer at the liquid–solid interface of (e, f) DRPs and (g, h) cylindrical pillars. Next, we developed an analytical model to exploit these temperature profiles to estimate the vapor force, F v . The assumptions underlying this quasi-steady state analytical model include (i) the temperature of the pillars in Region 1, wherein water penetrates, is the same as that of the water (based on simulation results in Fig. 4e–h ), (ii) heat transfer is dominated by conduction and the contribution of convection is negligible , (iii) the flow of vapor across the microtexture is laminar; i.e., the characteristic Reynolds number is low, which is typical for the Leidenfrost phenomenon (at 100°C, the density of the water vapor is r = , the viscosity of the water vapor is µ = −5 Pa × s, and the speed of the vapor in the microtexture is U ≈
1 m/s; thus, Re = r v UH / µ is ~1)
3, 45 . The energy for water evaporation ( Q eva (W)) is received through the heat influx at the water interface, Q (W) (Region 1 in Figs. 4b and ). This heat transfer rate, Q , is also equal to the heat transfer rate absorbed from the substrate underneath Q (W) (Region 2 in Figs. 4b and ). Therefore, 𝑄 -.! = − /0 ! /1 ℎ = 𝑄 , (1) where 𝑄 , = 𝑞 𝐴 % , where 𝑞 (W/m ) is the heat flux from the hot surface. A = L (m ) is the projected area of the periodic unit ( Fig. 5a ), m e (kg) is the mass of evaporated water in a periodic cell, t (s) is time, and h fg (J/kg) is the latent heat of the evaporation of water. Accounting for heat conduction across the vertical direction, Eq. (1) yields 𝜌 𝑢 𝐴 ℎ = 𝑘
722 /8/6 𝐴 % (2)
14 where r v (kg/m ) is the density of the water vapor, u z (m/s) is the vapor velocity along the z-axis, A v = 𝐿 , − π(𝐷 2⁄ ) , (m ) is the vapor escape area inside the periodic cell, and T (°C) is the temperature. The thermal gradient, dT/dz (°C/m), depends on the penetration depth of water into the microtexture. Considering the droplet keeps moving up and down during heating ( Fig. 5a ), the average thermal gradient, dT/dz was taken as ( T s − T sat ) / ( H − L p /2) in the calculation. Here, T s (°C) and T sat (°C) are the surface temperature of substrates and the boiling point of water, respectively; H (m) is the height of DRPs, and L p is the penetration depth when the droplet moves down, ranging within, 0 < L p < H . Fig. 5| Water T L as a function of penetration depth, pillar height, and contact angle. a) Schematic representing the flow of water vapor escaping through the microtexture. b) Water T L as a function of water penetration depth (rendered dimensionless by dividing it with the pillars’ height, H = 50 µ m). c) Water T L as a function of pillar height (rendered dimensionless by dividing it with H ). d) Water T L as a function of the actual (or intrinsic) contact angle.
15 Note that k eff (W/(m × K)) is the effective heat transfer coefficient determined by the area average of air and the solid fractions : 𝑘 = 𝜀𝑘 + (1 − 𝜀)𝑘 (3) where k v (W/(m × K)) and k s (W/(m × K)) are the thermal conductivities of the vapor and the solid, respectively, and e is the ratio of the vapor area to the projected area, ε = (𝐿 , − π(𝐷/2) , )/𝐿 , . Since the radius of the droplet’s base ( R b (m), ~0.7 mm from experimental observation) is dramatically larger than the height of the pillars, i.e., R b >> H , the lubrication analysis can be applied . As the vapor is produced from the bottom of the droplet and released from the spaces between pillars, the continuity equation can be written as ∫ 𝑢 : = 𝜋𝑟 , 𝑢 , (4) where u r (m/s) is the radial vapor velocity. Considering an approximately linear temperature distribution along the pillars , we can obtain the vapor velocity along the z-axis: 𝑢 ≈ < !"" (8 >8 )@A & B "’ (4> ( /,) (5) Thus, the average vapor velocity at the r-direction ( 𝑢 : NNN ) can be obtained: 𝑢 : NNN = D4 ∫ 𝑢 : 𝑑𝑧 = :E ) ,4 = 𝑟 < !"" (8 >8 ),@A & B "’ ( /,) (6) After applying lubrication analysis, the momentum equation for flow in the porous media can be reduced to classical Brinkman Equation : 𝜇 * E + F6 * = 𝜀 /G + /: + H & @I 𝑢 : (7) where P r (Pa) is the pressure distribution at the r-direction, µ v (Pa × s) is the dynamic viscosity of vapor, and K (m ) is the permeability of the porous media . We apply the non-slip boundary conditions and estimate the average radial velocity as 𝑢 : NNN = D4 ∫ 𝑢 : 𝑑𝑧 = IH & P /G + /: Q P
J!KL(4M/,)4M/, − 1Q (8)
16 where 𝛽 = S𝜀/𝐾 . Next, we substitute Eq. (6) into Eq. (8), and applying the boundary conditions, we obtain the radial pressure profile as 𝑃 : = (𝑅 N, − 𝑟 , ) H & < !"" (8 >8 )O@A & B "’ ( /,)I P1 −
J!KL(4M/,)4M/, Q >D . (9) Thus, the vapor force, F v (N), exerted by this vapor film onto the droplet can be estimated by 𝐹 = ∫ 𝑃 : P , ; = 𝜋𝑅 NO H & < !"" (8 >8 )Q@A & B "’ ( /,)I P1 −
J!KL(4M/,)4M/, Q >D . (10) Furthermore, the meniscus adhesive force on each DRP, F adi (N), is estimated as 𝐹 R/ - = 𝜋(𝐷 + 2𝑆 B )𝛾 ST (1 + cos 𝜃 U ), (11) where D (m) is the diameter of the top surface of DRPs, S h (m) is the horizontal length of edge DRPs, g LV (N/m) is the surface tension of water, and 𝜃 % is the actual (or intrinsic) contact angle of water on the silica surface. The Leidenfrost phenomenon occurs when the vapor force ( F v ) overcomes the adhesive force ( F ad ) and gravitation force ( G ), which allows the droplet to be levitated by a continuous vapor layer: 𝐹 = 𝐺 + ∑ 𝐹 R/ - , VWXD (12) where G = mg is the gravitation force, m (kg) is the mass of the droplet, g (m/s ) is the gravitation acceleration, and n is the number of DRPs in contact with the water droplet and can be estimated as 𝑛 = YP ,* * . At the onset of the Leidenfrost phenomenon, T s is equal to T L . Substituting Eqs. (10) and (11) into Eq. (12), the temperature corresponding to T L is 𝑇 − 𝑇 = Q@A & B "’ ( /,)IYP ,. H & < !"" × \P YP ,* * Q 𝜋(𝐷 + 2𝑆 B )𝛾 ST cos 𝜃 : + 𝐺] × P1 − J!KL(4M/,)4M/, Q (13) Equation (13) captures the dependence of the T L on the thermal, structural, chemical, and inertial characteristics of the system, especially the heat conduction at the water–substrate interface that scales with the penetration length ( L p ), microtexture dimensions, liquid–solid
17 adhesion, and the vapor flow through the microtexture. As the penetration depth, L p , increases, T L drops for both microtextures (Fig. 5b). To disentangle the effects of geometry and surface chemistry, we utilized the model to predict the variation in the T L of water with pillar height and surface chemistry. We found that T L in both cases increased with increased pillar height ( Fig. 5c ) because the taller the pillars are for a given pitch, the larger is the space to escape for the vapor, lowering F v . The model also pinpointed that as the surface hydrophobicity increases, T L drops ( Fig. 5d ). Therefore, the onset of the Leidenfrost phenomenon depends on a delicate balance among F v , G , and F Ad , which cannot be satisfied by all superhydrophobic surfaces. Thus, these forces could be balanced to engineer superhydrophobic surfaces that can prevent the onset of the Leidenfrost phenomenon even better than hydrophilic and superhydrophilic surfaces. Discussion
These experimental and theoretical findings might aid the rational design of superhydrophobic surfaces in thermal power engineering, for instance, for reducing frictional drag without compromising interfacial heat transfer. We illustrate this by comparing the time it takes to evaporate a drop of water placed on silica with DRPs with that on Glaco-coated silica, both maintained at 200 ° C. Although both surfaces were superhydrophobic, the water evaporated much faster on silica with DRPs than that on Glaco-coated silica (
Fig. 6 ). Silica with DRPs facilitated a 300% higher heat transfer than Glaco-coated silica because the former suppressed the Leidenfrost phenomenon to well above 200 ° C, whereas the latter experienced film boiling at T L ≈ 130 ° C. Fig. 6| Evaporation of water droplets on superhydrophobic surfaces at a surface temperature of 200 ° C. a) Silica with DRPs. b) Glaco-coated silica. Penetration depth was recognized as a crucial factor influencing the water–substrate heat transfer, and for a given microtexture, it depends on the Weber number . At low Weber numbers, such as in the experiments reported here, DRPs pinned the liquid meniscus at the doubly reentrant edge such that any further penetration did not increase the liquid–solid contact area. The DRP architecture helped us disentangle the effects of the penetration depth and the droplet–substrate contact area on the onset of the Leidenfrost phenomenon. However, at higher Weber numbers, the liquid meniscus might de-pin from doubly reentrant edges and land onto the DRPs’ stems . Subsequently, akin to the case of cylindrical pillars, the curvature of the liquid–vapor interface will reverse, the liquid-solid contact area will increase, and the conductive heat transfer will yield the Leidenfrost phenomenon. However, as the droplet’s momentum dissipates and once it touches the solid bottom surface, DRPs’ original behavior, suppressing the Leidenfrost phenomenon, might get restored. These predictions necessitate a systematic investigation. Notably, the coating-free aspect of the DRP architecture would offer resilience against harsh operational and cleaning protocols and organic fouling as perfluorinated coatings/materials are vulnerable to these conditions . In fact, DRP-based coating-free design principles are being exploited for mitigating cavitation and water desalination via scalable manufacturing platforms .
19 In conclusion, we have demonstrated that it is possible to realize superhydrophobic surfaces that do not suppress the T L of water via a judicious choice of surface microtexture and chemistry. This design can yield superhydrophobic surfaces or coatings for simultaneously reducing frictional drag and preempting boiling crisis in thermal machinery, such as heat exchangers, boiler pipes, and turbine blades. The coating-free aspect of the bioinspired DRP architecture is also enticing for its robustness against harsh conditions and organic fouling. Methods
Sample preparation
Silicon wafers (p-doped, 4-inch diameter, and 500 µ m thickness with a 2.4- µ m-thick thermally grown silica layer) were used as substrates ( Fig. 1a ). To fabricate DRPs (
Fig. 1d ), we applied photolithography, dry etching, and thermal oxide growth, followed by the protocols reported in the previous studies.
23, 37
Besides DRPs, we microfabricated silica surface adorned with arrays of cylindrical pillars, which were rendered superhydrophilic (
Fig. 1c ), manufactured by the same materials above, and the pitch, height, and top diameter of cylindrical pillars were the same as those of the DRPs. After surface fabrications and before droplet experiments, all samples were cleaned with Piranha solution (H SO : H O = 3:1) for 15 min at 110 ° C and stored in the vacuum oven for 24 h at 50 ° C and ~0.03 bar. Consequently, the actual (or intrinsic) contact angle ( 𝜃 % ) of water on the silica surface was stabilized at ~40 ° . For comparison, Glaco coating was used to coat flat silica (same materials in the DRPs) as the common superhydrophobic surface ( Fig. 1b ). The effect of surface chemistry was studied by coating perfluorodecyltrichlorosilane (FDTS) on DRPs using a molecular vapor deposition system (ASMT 100E), as described previously . All the surface structures were examined by a scanning electron microscope (Helios Nanolab SEM). Contact angle measurements
We characterized the wettability of these surfaces with water by measuring apparent CAs of water droplets (6 µ L) and advancing and receding angles at a rate of 0.2 µ L/s. The water CAs on the surfaces of different samples were measured using a drop shape analyzer (Kruss: DSA100) interfaced with the Advance software. The sliding angles of water droplets (10 µ L) on these surfaces were measured by an automatic tilting stage, where a camera was used to record the movement of droplets. All experiments were repeated at least three times.
Droplet generation for Leidenfrost experiments
To generate water droplets, we used a Harvard Apparatus syringe pump (PHD ULTRA: 703007INT) and a 60-mL plastic syringe connected by a plastic syringe extension tube. A stainless-steel capillary (diameter: ~0.51 mm) was used to dispense water perpendicular to heated test surfaces. The flow speed was set as 100 µ L/s, and the syringe system continuously generated 15- µ L droplets (radius: ~1.5 mm, which is smaller than the capillary length of water ~2.7 mm at NTP).
Confocal experiments
A Zeiss upright laser confocal microscope (Model number: LSM710) was used to detect the meniscus at the interface between a droplet and pillars. A ~100- µ L water droplet with Rhodamine B (Acros) was gently placed on samples using a manual syringe, and a 20X micro-objective lens was immersed into the droplet before acquisition. The 520-nm laser at the Z-stack mode was used to scan the droplet–surface interface. Subsequently, using Imaris v.8.1 by Bitplane, we obtained cross-sectional images to visualize the meniscus shape (
Supplementary Materials, Fig. S5b-c ). High-speed imaging
A high-speed camera (model number: Phantom V1212) with a zoomed lens was applied to record images at a high frame rate. To ensure that all the images have the same configuration, the imaging parameters of the camera such as frame per second (FPS), image pixels, and exposure time in all experiments were set to 1000 FPS, 512 ´
512 pixels, and 200 µ s,
21 respectively. A triaxial stage was used to adjust the position of the samples with respect to the high-speed camera. The hot plate was leveled using a bubble leveler, and the images were calibrated by the stainless-steel capillary. For image processing, Phantom Camera Control 2.8 software was used.
Leidenfrost experiments
A hot plate (temperature: 20°C–540°C) was placed on a triaxial stage and balanced by a bubble level. Considering the nonequilibrium nature of the hot plate, the actual temperatures of the hot plate surface were measured by a K-type thermocouple and recorded by a temperature data logger (Omega: RDXL6SD). We moved the capillary tip with a droplet as close to the surfaces as possible to gently release droplets from the capillary onto the surfaces to minimize the droplet impact. Thus, all the droplets keep the same low Weber number ( We » µ L) were generated by a syringe pump and released from the needle tip to the heated sample surface. A high-speed camera was used to record the impact, bounce, and vibration of droplets and observe liquid–vapor–solid interactions at the interface between droplets and surfaces. All the experiments were repeated at least three times. To measure Leidenfrost points on different surfaces, the classical lifetime method was used to measure T L on these surfaces. The first peak in the curve between the lifetime and surface temperatures yields T L10, 25 . To measure the lifetimes of droplets, the temperature of the hot plate increased by 10°C at each step, and the droplet evaporation times were recorded by a timer. All experiments were repeated at least five times. Since there is usually no clear lifetime peak when the Leidenfrost transition occurs on superhydrophobic surfaces, we used the high-speed camera to detect T L on Glaco-coated flat silica and silica with DRPs by observing the moments of vapor film formation and droplet bouncing while gradually increasing the temperature. Calculations
Simulation results were obtained by utilizing the heat transfer module in COMSOL (5.4) (the detailed boundary conditions are shown in
Supplementary Materials ). Equations 1–13 were solved using the commercial software Wolfram Mathematica (11.0). The thermal properties of water were taken at 99 ° C, and other materials’ properties were taken at 225 ° C, which is an average temperature between the water and the hot substrate.
Acknowledgments
This work was supported by King Abdullah University of Science and Technology (KAUST). M.S. thanks Prof. Sigurdur Thoroddsen from KAUST and Prof. Shangsheng Feng from Xi’an Jiaotong University for fruitful discussions. References
1. Leidenfrost JG.
De aquae communis nonnullis qualitatibus tractatus . Ovenius, 1756. 2. Schutzius TM, Jung S, Maitra T, Graeber G, Köhme M, Poulikakos D. Spontaneous droplet trampolining on rigid superhydrophobic surfaces.
Nature (7576) : Annual Review of Fluid Mechanics
Nature (7415) : Nature Physics (12) : Nature Physics (11) : Nature Physics (5) : , et al. Directional transport of high-temperature Janus droplets mediated by structural topography.
Nature Physics , et al.
Self-propelled Leidenfrost droplets.
Physical review letters (15) : International Journal of heat and mass transfer (11) : Sci Rep Physics of Fluids (4) : Nature communications , et al. Green chemistry and nanofabrication in a levitated Leidenfrost drop.
Nature communications
23 15. Shrestha BR, Pillai S, Santana A, Donaldson Jr SH, Pascal TA, Mishra H. Nuclear Quantum Effects in Hydrophobic Nanoconfinement.
The Journal of Physical Chemistry Letters (18) : Appl Phys Lett (8) : Applied Physics Letters (20) : , et al. Vitrification and levitation of a liquid droplet on liquid nitrogen.
Proc Natl Acad Sci U S A (10) : , et al. Asymmetric wettability of nanostructures directs Leidenfrost droplets.
ACS nano (1) : , et al. The Leidenfrost temperature increase for impacting droplets on carbon-nanofiber surfaces.
Soft matter (13) : Surface Innovations (1) : Science (6213) : Journal of Colloid and Interface Science
Capillarity and Wetting Phenomena: Drops, Bubbles, Pearls, Waves . Springer; 2004 edition, 2003. 25. Bourrianne P, Lv C, Quéré D. The cold Leidenfrost regime.
Science advances (6) : eaaw0304. 26. Saranadhi D, Chen D, Kleingartner JA, Srinivasan S, Cohen RE, McKinley GH. Sustained drag reduction in a turbulent flow using a low-temperature Leidenfrost surface. Science advances (10) : e1600686. 27. Dupeux G, Bourrianne P, Magdelaine Q, Clanet C, Quere D. Propulsion on a superhydrophobic ratchet. Sci Rep Langmuir (42) :
24 29. Quéré D. Superhydrophobic surfaces: Leidenfrost becomes a fakir.
Nature materials (11) : Energy
Liquid-vapor phase-change phenomena : an introduction to the thermophysics of vaporization and condensation processes in heat transfer equipment , 2nd edn. Taylor and Francis: New York, 2008. 32. Theofanous T, Dinh T-N, Tu J, Dinh A. The boiling crisis phenomenon: Part II: dryout dynamics and burnout.
Experimental Thermal and Fluid Science (6-7) : International Journal of Multiphase Flow (1-2) : Experiments in Fluids (176). 35. Das R, Ahmad Z, Nauruzbayeva J, Mishra H. Biomimetic Coating-free Superomniphobicity. Scientific Reports (1) : , et al. Rendering SiO2/Si Surfaces Omniphobic by Carving Gas-Entrapping Microtextures Comprising Reentrant and Doubly Reentrant Cavities or Pillars.
JoVE : e60403. 37. Domingues EM, Arunachalam S, Mishra H. Doubly Reentrant Cavities Prevent Catastrophic Wetting Transitions on Intrinsically Wetting Surfaces. ACS Applied Materials & Interfaces (25) : Nature Communications (1) : Industrial and Engineering Chemistry (8) :
7. 40. Kim SH, Lee G, Kim H, Kim MH. Leidenfrost point and droplet dynamics on heated micropillar array surface.
International Journal of Heat and Mass Transfer , et al.
Solid-Liquid Work of Adhesion.
Langmuir : the ACS journal of surfaces and colloids (15) : , et al. Rendering SiO2/Si Surfaces Omniphobic by Carving Gas-Entrapping Microtextures
25 Comprising Reentrant and Doubly Reentrant Cavities or Pillars.
JoVE (Journal of Visualized Experiments) : e60403. 43. Mahadik GA, Hernandez-Sanchez JF, Arunachalam S, Gallo Jr A, Cheng L, Farinha AS , et al. Superhydrophobicity and size reduction enabled Halobates (Insecta: Heteroptera, Gerridae) to colonize the open ocean.
Scientific Reports (1) : COMSOL AB, Stockholm, Sweden . 45. Lide DR.
CRC Handbook of Chemistry and Physics , vol. 79. CRC Press, 1999. 46. Pinkus O, Sternlicht B, Saibel E. Theory of hydrodynamic lubrication. 1962. 47. Adera S, Raj R, Enright R, Wang EN. Non-wetting droplets on hot superhydrophilic surfaces.
Nature communications . 48. Sangani AS, Acrivos A. Slow flow past periodic arrays of cylinders with application to heat transfer. International journal of Multiphase flow (3) : , et al. Droplet impact on superheated micro-structured surfaces.
Soft Matter (12) : Advanced Materials Interfaces (22) : Journal of Membrane Science
Water research
Science Advances (13) : eaax6192. 54. Das R, Arunachalam S, Ahmad Z, Manalastas E, Mishra H. Bio-inspired gas-entrapping membranes (GEMs) derived from common water-wet materials for green desalination. Journal of Membrane Science Supplemental Information Unexpected Suppression of Leidenfrost Phenomenon on Superhydrophobic Surface
Meng Shi, Ratul Das, Sankara Arunachalam, Himanshu Mishra*
Interfacial Lab (iLab), Water Desalination and Reuse Center (WDRC), Biological and Environmental Science and Engineering (BESE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia *Corresponding author: [email protected] Figure S1 Droplets’ behaviors on flat silica surfaces maintained at different temperatures: the Leidenfrost phenomenon appeared at T L ≈ 286 ° C. Representative snapshots of droplet behavior when the surface temperature was below T L (76 ° C, 158 ° C, and 245 ° C) and above T L (347 ° C and 441 ° C). Although the droplets bounced after initially contacting the hot surface at 245 ° C, we observed transition boiling at the bottom of the droplet after several bounces. Therefore, the droplet is in transition boiling at 245 ° C, which can also be seen from the lifetime curve (
Figure 2e ). Figure S2 Droplets’ behaviors on silica with cylindrical pillars maintained at different temperatures: the Leidenfrost phenomenon appeared at T L ≈ 314 ° C. Representative snapshots of droplet behavior when the surface temperature was below T L (76 ° C, 158 ° C, and 245 ° C) and above T L (347 ° C and 441 ° C). Although droplets bounced when initially contacting the hot surface at 245 ° C, we observed transition boiling at the bottom of the droplet after several bounces. Hence, the droplet is in transition boiling at 245 ° C, which can also be deduced from lifetime curve (
Figure 2e ). Figure S3 Droplets’ behaviors on Glaco-coated flat silica maintained at different temperatures: the Leidenfrost phenomenon appeared at T L ≈ 130 ° C. Representative snapshots of droplet behavior when the surface temperature was below T L (76 ° C) and above T L (158 ° C, 245 ° C, 347 ° C, and 441 ° C). Figure S4 Droplets’ behaviors on silica with DRPs maintained at different temperatures: the Leidenfrost phenomenon appeared at T L ≈ 363 ° C. Representative snapshots of droplet behavior when the surface temperature was below T L (76 ° C, 158 ° C, 245 ° C, and 347 ° C) and above T L (441 ° C). ° C) and found that the water droplet was levitated by a vapor film and started bouncing on the FDTS-coated DRPs (
Figure S5a ). Through confocal experiments, we obtained water droplet meniscus curvatures at the interfaces of FDTS-coated and untreated DRPs. On the untreated DRPs, a water droplet meniscus goes down to the doubly reentrant edges (
Figure S5b ). After FDTS coating, a water droplet can stay only on the top of the DRPs and cannot wet the edges (
Figure S5c ). Thus, the adhesion is significantly reduced by not only increasing the actual (intrinsic) contact angle q o to 112° but also changing the contact length slightly. Thus, the local pinning effect is eliminated from the surface of the FDTS-coated DRPs. Figure S5 Effect of hydrophobic coating on the Leidenfrost phenomenon on silica with doubly reentrant pillars. a) Droplet bouncing on FDTS-coated DRPs at 347°C; b) Droplet meniscus at the interface between the droplet bottom and DRPs; c) Droplet meniscus at the interface between the droplet bottom and FDTS-coated DRPs. Figure S6 Droplet detachment from cylindrical pillars at 347°C.