Research Article | Open Access

Yue Hu, Kaili Jiang, Kim Meow Liew, Lu-Wen Zhang, "Nanoarray-Embedded Hierarchical Surfaces for Highly Durable Dropwise Condensation", *Research*, vol. 2022, Article ID 9789657, 13 pages, 2022. https://doi.org/10.34133/2022/9789657

# Nanoarray-Embedded Hierarchical Surfaces for Highly Durable Dropwise Condensation

#### Abstract

Durable dropwise condensation of saturated vapor is of significance for heat transfer and energy saving in extensive industrial applications. While numerous superhydrophobic surfaces can promote steam condensation, maintaining discrete microdroplets on surfaces without the formation of a flooded filmwise condensation at high subcooling remains challenging. Here, we report the development of carbon nanotube array-embedded hierarchical composite surfaces that enable ultra-durable dropwise condensation under a wide range of subcooling (), which outperforms existing nanowire surfaces. This performance stems from the combined strategies of the hydrophobic nanostructures that allow efficient surface renewal and the patterned hydrophilic micro frames that protect the nanostructures and also accelerate droplet nucleation. The synergistic effects of the composite design ensure sustained Cassie wetting mode and capillarity-governed droplet mobility () as well as the large specific volume of condensed droplets, which contributes to the enhanced condensation heat transfer. Our design provides a feasible alternative for efficiently transferring heat in a vapor environment with relatively high temperatures through the tunable multiscale morphology.

#### 1. Introduction

Condensation is a ubiquitous phenomenon in nature and plays a key role across a wide range of industrial applications, such as anti-icing, power generation, sewage treatment, water desalination, and energy harvesting [1–3]. To enhance the condensation performances for higher heat transfer efficiency and further substantial savings in energy [1, 4], efficient nucleation of droplets [5] and rapid shedding [6, 7] are essential in surface designs. Considerable efforts have demonstrated that dropwise condensation, normally on hydrophobic surfaces, exhibits efficient heat transfer over filmwise condensation [3, 8]. Facilitated by lithography, etching, and the other straightforward approaches for surface modification, various anti-wetting designs are emerging, including pillars [9, 10], pit-array surfaces [11], multiscale surfaces [12–14], and hybrid surfaces [4, 15]. The customized and sophisticated designs enable control of the surface roughness, and optimization of the microtopographic parameters, thereby allowing regulation of the liquid-solid contact mode, further affecting the droplets dynamics, and manipulating the condensation performance [16].

However, sustaining droplet condensation at a high degree of subcooling (difference between the temperature of saturated vapor and the condensation surface [17]) remains challenging. Increasing subcooling will result in elevated nucleation density and condensation rate, which may exceed the droplet departure rate and form a liquid film [8, 18]. Therefore, the film-like droplets that pin on the surface will undergo the irreversible Cassie–Wenzel wetting transition, leading to the phenomenon of flooding that declines the heat transfer efficiency of the surface [19, 20]. In efforts to maintain the Cassie state of the droplet, previous studies have focused on decreasing the solid-liquid contact area [21, 22] or the number of droplet nucleation sites [23, 24] by designing increasingly smaller arrays or approximate arrays, from micro to nanoscale size, such as pine needle shapes [6], microconical architectures [25], nanoscaffolds [26], microratchet arrays [27], and vertical arrays of carbon nanotube (CNTs) [28, 29]. These surfaces indeed exhibit an enhancement in dropwise condensation by collecting a volume of water three times higher than that in the case of plate surfaces [30], achieving a 100% higher heat flux than that on the plane hydrophobic surface [31], and reaching a directional transport efficiency of approximately 80% for tiny droplets [27].

Nevertheless, two fundamental contradictions prevent the applications of micro/nanoarrays at large subcooling. The first one is the paradox of diminishing array size and mechanical robustness. Exposed array structures with high aspect ratios are considerably fragile and suffer from disruptions due to buckling, clustering, and fracture, even under a small external disturbance [32]; these disruptions eventually affect the wetting state and weaken the hydrophobic effect [33]. A recently armor-protected superhydrophobic surface presents the remarkably robustness and inspires the development of hierarchical surface [34]. The second one is the contradiction between high droplet mobility and accelerated nucleation. Hydrophobic array surfaces promise rapid droplet removal that is the key to dropwise condensation, but has a negative effect on droplet nucleation, especially at the initial stage of condensation [35, 36]. Conversely, hydrophilic surfaces promote the nucleation rate of droplets, but lead to undesired flooding.

In this study, we report an integrated design of a self-protected hierarchical hybrid surface that is composed of micropatterned surface structures and vertically arrayed carbon nanotubes (VACNTs) grown internally and directly. We propose that this approach is a complementary and enhanced strategy for achieving stable dropwise condensation efficiency at high subcooling. The robustness of our hierarchical patterned surfaces is demonstrated through a cyclic loading test. The CNT array, which is grown under control in restricted frames, is expected to act as the major hydrophobic component because of the nanoscale contact sites and the intrinsic water-repellency of CNTs. The hydrophilic micro-substrate patterned as a frame is expected to play a critical role in protecting the slender nanoarrays and providing essential hydrophilic sites that accelerate nucleation without compromising the global hydrophobicity of the surface. We perform a comprehensive analysis of the coalescence, jumping, and removal of droplets on the hierarchical surfaces. A heterogeneous nucleation model is developed for the nanoarray surface in the Cassie wetting state to reveal the underlying mechanism for dropwise condensation. Our experiments and molecular simulation are performed to show that the droplets and the density distribution of the water molecules are affected by the microtextured silicon (Si) substrate. This work enables a new microstructured surface with durable dropwise condensation and enhanced mechanical robustness at high subcooling under cyclic mechanical loading.

#### 2. Results

##### 2.1. Hierarchical CNT Array-Planted Patterned Surface

We designed and fabricated a hierarchical surface in which VACNTs are planted in a carefully patterned Si substrate (Figure 1). This bi-scale surface design involves two major steps. To construct the microscale structure, the Si wafer is patterned in triangular, circular, square, and hexagonal tessellation through photolithography. Unlike pillars, the interconnected indentations act as a stiff frame that plays the vital roles of protecting and limiting the deformation of inner nanostructures. Figure 1(a) shows the hexagonal pattern, with photolithographic cavities having a height of 25 *μ*m, an inter-cavity distance of 5 *μ*m, and cavity sizes of 50 *μ*m, 100 *μ*m, 150 *μ*m, and 200 *μ*m. The geometry and sizes of the triangular, circular, and square pits are shown in Figure S1. With regard to the nanostructure, we intend to plant vertical CNTs in the interconnected cavities by using chemical vapor-phase deposition (CVD), thus generating a hydrophobic surface with nanoscale contact sites. The efficient CVD approach for VACNTs has been explored in our previous studies [29, 37], and a detailed description is given in the methods. Undesired CNTs that grow on the walls of the microcavities can be eliminated by removing the photoresist mask, thus leaving only the VACNTs planted in the notches. This planting strategy ensures that the dense, slender nanostructures occupy most of the surface area and thus determine the global wettability of the surface. The contact sites between the intrinsically hydrophobic CNT arrays and the vapor offers reduced liquid-solid contact area, thereby promoting the mobility of droplets and increasing the condensation efficiency [38].

##### 2.2. Dropwise Condensation on the Hierarchical Surface

To understand the condensation dynamics on our hierarchical surfaces, we experimentally investigated the nucleation, growth, coalescence, and removal of the droplets under saturated water vapor. Figure 2(a) confirms stable dropwise condensation on the hierarchical surface without film formation for a wide range of subcooling, up to a maximum of 38 K; the main reason for this observed behavior is that CNT arrays contribute to the surface superhydrophobicity. The wetting state measurements using optical microscopy show that the apparent CA, , reaches a maximum value of 154 and has an average value of 151, with a roll-off angle of approximately 3.4 across a wide scenario (Figure 2(c)). As the Cassie state is maintained during the condensation process, the surface enables coalescence-induced jumping of neighboring droplets (yellow dashed lines in Figure 2(a)) on the condensed surface, subsequently exposing a completely renewed surface (blue dashed lines in Figure 2(a)) in the saturated vapor environment for a new cycle of condensation and heat exchange. Figure 2(b) depicts an entire cycle, illustrating the nucleation, coalescence, jumping, and shedding of droplets on the hierarchical surface. The coalescence of nucleated droplets through jumping contributes to droplet mobility, resulting in a fast renewal of the surface and an increase in heat transfer efficiency [39]. From the thermodynamic perspective, after the droplets in the steam environment grow and coalesce, the surface energies of the individual droplets () on the coalescence path are converted to the surface energy of a merged droplet (); a certain amount of this energy is released (), which provides the kinetic energy for the jumping of droplets and drives coalescence.

Noted that in comparing to Si nanowire, copper nanowire, and other nanoarray surfaces that undergo a transition from dropwise to filmwise condensation at an averaged surface subcooling 20 K in subcooling tolerance, our hierarchical surfaces achieve an improvement of around 18 K in subcooling tolerance (Figure 2(d)). The dropwise condensation process at =8 K−38 K is shown in Videos S1, S2, and S3. In all three cases, droplets move and depart by the coalescence-induced jumping rather than sliding, and the probability of droplet departure decreases with the increasing . This is due to the high heat flux caused by the large form uncontrolled nucleation in some microgaps between CNTs, which provides large adhesive force on droplets [31]. To explain the maintenance of dropwise condensation at high , we introduce the dimensionless Bond number () to determine the dominant driving force in the movement of a droplet. Defined in terms of the ratio between the gravitational and capillary forces, is given by (where is the condensate density, is the gravitational constant, is the droplet radius just before the jumping movement, and is the surface tension), indicating capillarity-controlled () or gravity-controlled droplet mobility (), corresponding to the desired dropwise and undesired filmwise condensation, respectively [36, 40]. In our experimental measurements, the maximum radius of the jumping droplets, , ranges from approximately 164 *μ*m () to approximately 613 *μ*m (), and falls within the range of 72.8 mN/m–67.1 mN/m as the temperature elevates from 293 K to 328 K [41], yielding a maximum value of in the range 0.004–0.055. Note that the predicted result for is within the well-known domain of capillary-governed droplet departure and dropwise condensation [36, 40]. The removal of the capillarity-dominated droplets relies on not only the gravity, but also the coalescence-induced jumping [42]. Therefore, the small number reveals the fundamental reasons for the maintenance of the spherical droplet shape and the efficiency of surface renewal at high subcooling.

To further identify the advantages of the hierarchical surface during dropwise condensation, we studied the formation of discrete droplets and compared it with the corresponding scenario for a micropatterned substrate without implanted nanoarrays. It is worth noting that the Si surface is intrinsically hydrophilic (), and even though reaches 121 at the room temperature of 24 °C (indicating hydrophobicity) after lithography (cavities with a feature size of 150 *μ*m) (Figure S2(a)), the formation of irregular films or puddles on the condensing patterned Si substrate is in sharp contrast to the dropwise condensation on CNT-implanted surfaces (Figures S2(b) and S2(c)). This discrepancy is manifested in the distribution of the diameters of discrete droplets. As evidenced in Figure 2(e), in the case of the hierarchical surfaces, 95% of the droplets have a diameter less than 60 *μ*m, whereas in the case of the pure patterned Si surface, a major portion of the surface is occupied by large irregular droplets (17% of the droplets have an equivalent diameter*μ*m) in the absence of nanoarrays. As the hierarchical surfaces have condensed droplets of smaller size with much higher , it plays a substantial role in declining the thermal resistance between vapor and the substrate during condensation, especially at high subcooling [14]. As another essential indicator of surface condensation, the thermal diffusivity coefficient of the solid substrate is quantified in Figure S3. It increases from 38 mm^{2}/s for a patterned Si surface to a maximum value of 87 mm^{2}/s and an average value of 75 mm^{2}/s for a hierarchical surface. This substantial increase of 97% in the thermal diffusivity indicates that the CNT-implanted surfaces exhibit better heat conduction performance, which endows the hierarchical surfaces with an enhanced condensation property.

##### 2.3. Heat Transfer Performance Dominated by the Nanoscale Arrays

To gain insights into the physical mechanisms underlying stable dropwise condensation, we laid emphasis on the significant contribution of nanoscale vertical arrays in enabling efficient heat transfer. Figure 3(a) shows air pockets at the solid-liquid interface, demonstrating that the CNT array maintains the Cassie wetting state, as remarked in previous studies [30, 43, 44]. Furthermore, as seen in Figure 2(a), the Si microframes do not have a pinning effect on the three-phase contact line of droplets, which remains circular rather than transforming into the shape of the pattern. These measurements evidence the role of nanoarrays in determining the wetting state and condensation efficiency at high subcooling.

Previous studies have proven that the surface properties have an obvious influence on the nucleation rate [8, 16, 45]. However, the description for vapor nucleation on the array-like surface is inadequate in terms of the surface heterogeneity. Here, we developed a heterogeneous nucleation model for a nanoarray surface in the Cassie wetting state and specifically investigated the contribution of the nanoarray on the heat transfer performance. As the Gibbs free energy for heterogeneous nucleation from classical nucleation theory is restricted to ideal smooth surfaces [46], we modified it by splitting the term for surface free energy into a term related to the spherical cap and a term related to the bottom contact area, while keeping the volume term unchanged: where is the outer surface area of the spherical cap of the droplet with radius, is the bottom area of the droplet cap, and represents the volume of the droplet with an apparent CA . , , and represent the liquid-gas, solid-liquid, and solid-gas surface tension, respectively. is the solid-liquid contact area fraction, and is the liquid-gas contact area fraction. is the difference in free energy per unit volume for water and supersaturated vapor at the same pressure. Subsequently, the heterogeneous nucleation energy barrier or the work of critical nucleation can be derived by using , which gives where satisfies based on Cassie theory, and according to Young’s equation, where is the intrinsic CA. Drawing from the expression for homogeneous critical nucleation energy from classical nucleation theory, can be reformulated as

is a critical scaling factor residing in the range [0, 1], termed as Fletcher factor. Although another form of was derived in a previous study for describing fractal surface [47], here it quantitatively characterizes the dependence of the surface geometry and the intrinsic wetting property on heterogeneous nucleation. A larger value of , i.e., closer to 1, implies that a higher energy barrier must be overcome to form condensed droplets and also suggests that critical heterogeneous nucleation tends to be extremely close to ideal homogeneous nucleation.

By leveraging the proposed model, Figure 3(b) elaborates the relationship between the Fletcher factor and the two critical characteristic parameters of a surface: the intrinsic wetting angle and the solid-liquid contact area fraction , which can be related to the density of nucleation sites in the Cassie state. The predictions indicate that surfaces with a larger or a smaller tend to suppress droplet nucleation. As CNT array surfaces have a of approximately 107 [48] and a small , which is identified to be 0.15 when measured using the image contrast approach (Figure S4), *F* is extremely close to 1.0, as shown in Figure 3(b) (represented by a red star). The large work-of-critical-nucleation value implied by the high indicates that it is difficult to pin droplets on the nanoscale arrayed surface, leading to favorable conditions for droplet mobility, surface renewal, and maintenance of dropwise condensation.

To quantify the effect of surface roughness on the condensation performance at altered subcooling, we studied the heat transfer efficiency and heat flux, which are predominantly influenced by the CNT array surfaces (Figures 3(c)–3(f)). The heat transfer models for an individual droplet and a group of droplets are used to characterize the heat transfer rate and heat flux (Supplementary Section S5). As illustrated in Figure 3(c), the heat transfer rate of a single droplet on a rough surface with (relatively small solid-liquid contact area) shows a decrease of 87% when compared with the corresponding value for a smooth surface with (relatively large solid-liquid contact area) at a high subcooling of 40 K. This decrease in the case of a rough surface is attributable to the air layer trapped in the cavities, which has a lower thermal conductivity than a smooth surface. With regard to the influence of the apparent CA and the droplet radius *r*, Figure 3(d) suggests that the condensation of a small droplet decreases the heat transfer rate to a greater extent than does the condensation of a large droplet and this decrease can be up to 1 order of magnitude for =90. For condensed droplets with the same radius, a large value of minimizes the heat transfer rate for individual droplets; therefore, the heat-sensitive surface tension of the droplet is maintained, preventing flooding of the surface.

The global heat transfer efficiency of the surface was further investigated with respect to variable . In Figure 3(e), the heat flux shows a maximum increase of 38%, when decreases from 0.8 to at =40 K. Furthermore, the surface with a smaller value of facilitates the condensation of a larger number of droplets with radius less than 100 *μ*m (Figure S5(b)), and the increased nucleation density , which denotes the number of nucleation sites per unit area on the condensing surface, increases the heat flux across the surface, as shown in Figure 3(f). Therefore, the nanoscale array structures raise the energy barrier for droplets nucleation and optimize surface heat transfer performance.

##### 2.4. Condensation Activation Governed by the Hydrophilic Microstructures

Although the uncovered narrow wall of the microframes occupies a maximum area of only 10% of the surface in our experiments (), the hydrophilic nature of Si is likely to have an advantageous effect on condensation efficiency. To determine the potential influence of this patterned Si surface occupying a small percentage of the area, we investigated the distribution of condensed droplets of different sizes as the tessellation length shortens from 200 to 50 , and the corresponding hydrophilic contact area enlarges from 3% to 10% (Figure 4).

At the very initial stage of condensation (10 s), microdroplets nucleate on the hydrophilic silicon frame, in the width of only 5 *μ*m (Figure S6). The surface with displays the largest coverage ratio of condensed droplets than others (Figure S7). After 30 s of condensation, as displayed in Figure 4(a), the size of the droplets increases considerably, and the nucleation density decreases, along with an increase in and a growing area exhibiting local hydrophilicity (original images shown in Figure S8). From the corresponding Voronoi analysis of droplet distribution, the local areal proportions around large droplets are apparently higher than those around small droplets (Figure 4(b)). As the local areal proportion is governed by the distance between adjacent droplets, the results suggest that a larger area of the heat transfer surface is exposed to the vapor environment as a greater number of large droplets condense on the hybrid hierarchical surface. Figure 4(c) shows additional statistics on the maximum diameter of condensed droplets, which decreases from 172 ± 12 *μ*m at *L* = 50 *μ*m to 97 ± 11 *μ*m at L = 200 *μ*m, indicating that the densely patterned surface favors the collection of droplets in large volumes (the entire size distribution is shown in Figure S9). However, this outcome gives rise to another issue—large condensed droplets result in a large coverage area that hinders heat transfer between the surface and vapor. To extract the dominant factor in nucleation efficiency at the initial stage of condensation, the trend for the specific volume, defined by , against pattern density is obtained, as illustrated in Figure 5(d), where denotes the spherical crown volume of the droplets, is the maximum projection area, and the average value of is 151, from the previous characterization. The results indicate that for the hierarchical surface, the specific volume when is 32% higher than the specific volume when . Although a large drop volume and a large coverage area constitute a paradox for condensation efficiency, the large-sized droplets favored by hydrophilic sites exhibit a balance between the volume of condensed droplets and the corresponding coverage area (Figure S10); thus, the hierarchical surface with a highly dense pattern favors condensation of highly efficient condensation without affecting the heat transfer performance.

To provide insights into the influence of hydrophilic microframes on condensation behaviors at the molecular level, we carried out a comparative study of a hybrid surface and a hydrophobic VACNT surface through molecular dynamics simulation (models are shown in Figure S11). The density map images show that fast nucleation of water molecules occurs in the hydrophilic domain, whereas numerous droplets tend to scatter on the VACNT (Figure 5(a)). From the kinetic energy distribution of H_{2}O near the hybrid surface (Figure 5(c)), it can be noted that the hydrophilic Si region promotes strong interaction between vapor and the substrate and exhibits apparently low kinetic energy, whereas the CNTs provide a large amount of kinetic energy for motion of H_{2}O, leading to the rebounding of water molecules and the prevention of nucleation [49]. Furthermore, the hydrophilic nucleation site affects the density distribution of water molecules away from the surface, which is observed in the statistics related to the number of water molecules per vapor bin (Figure 5(b)). This behavior is attributable to the equilibrium process of evaporation and condensation for the nucleated droplets, which promotes the collection of a great number of adjacent water molecules near the surface and leads to apparent aggregation in the subsurface (15–45 ). The increased cluster size and the rate of condensation nucleation on the hierarchical surface show that the hydrophilic microstructure provides localized activation sites that accelerate the nucleation and growth of droplets (Figures 5(d) and 5(e)). The difference between the vapor and the wall temperature increases the number of nucleation cluster on the VACNT surface with smaller cluster size, but is not apparent on the hierarchical surface. Thus, the hierarchical hybrid design that integrates hydrophilic microframes with hydrophobic CNT arrays promotes the rate of droplet nucleation and improves condensation efficiency while maintaining both surface hydrophobicity and droplet mobility.

##### 2.5. Mechanical Robustness of the Hierarchical Surface

Unlike fragile array structures featuring a large ratio between length and diameter, our hierarchical material enhances structural robustness when VACNTs are planted inside the interconnected cavities. We conducted cyclic scratching and compression tests to examine the stability of the designed hierarchical surfaces in the wetting state. Under cyclic scratching, the nanoarrays remain in the cavities with only slight damage to the hierarchical surfaces (Figure 6(c)); in contrast, on the substrate that is not protected by stiff frames, no trace of the CNTs remains (Figure 6(b)). Under cyclic compressive loading, cracks appear on the array surface (Figure 6(e)); however, the patterned hierarchical surface remains intact, without noticeable damage (Figure 6(f)). To test the robustness of the wetting property, the average CAs were measured after each scratch and compression load during cyclic testing. As shown in Figure 6(g), of the hierarchical surface remains greater than 130 after 13 rounds of cyclic scratching, whereas the bare VACNT surface loses hydrophobicity. After the compression cycle, of the hierarchical surface remains greater than 145, with a maximum decrease in value of 2.8%, which is superior to the unprotected array, whose hydrophobicity weakens significantly (to approximately 125) (Figure 6(h)). Thus, the interconnected micro-indentations offer a stiff frame for the fragile CNTs, ensuring mechanical robustness of the hierarchical structure under interruptions and preventing a wetting state transition.

#### 3. Discussion

In this study, we fully exploit the advantages of distinct scales of array structures with different wetting characteristics, creating a superhydrophobic hierarchical patterned surface by integrating hydrophobic nanoscale vertical CNT arrays and hydrophilic microscale frames. The results show that the surface exhibits high endurance in a wide range of subcooling, 8 K–38 K, and exhibits highly stable dropwise condensation with good condensation efficiency; the surface outperforms conventional hydrophobic nanoarray/microarray surfaces that switch to flooding at an average of 20 K [4, 9, 13, 14, 21–23, 26, 50]. Our experimental and theoretical study reveals that the high-endurance dropwise condensation is mainly attributable to the CNT arrays, which improve the surface properties in terms of high water repellency, better thermal performance, and reduced solid-liquid contact area owing to the nanoscale contact area fraction. Previous studies have shown the benefits of using CNT arrays in steam generation, evaporation-assisted detection, and dew harvesting [30, 51, 52], while in this study, we demonstrate the superior potential of CNT arrays for efficient condensation in extreme environments.

Through further exploring the effects of CNT arrays on critical nucleation and heat transfer, our surface with implanted CNT arrays shows a substantial enhancement of 97% in the thermal diffusivity coefficient when compared with a micropatterned Si surface without implanted CNT arrays. Furthermore, we find that the hydrophilic substrates patterned as microframes accelerate droplet nucleation by providing active nucleation sites for saturated vapor without compromising on the hydrophobicity of the entire surface. The specific volume of droplets increases by 35%, indicating a larger droplet size for the same coverage area and resulting in higher condensation efficiency, which is determined through kinetic energy analysis from molecular simulation. The interconnected microframes promise regionalization and confinement of the planted CNT array under mechanical perturbation, thereby conferring structural stability to of the surface. Thus, the integrated design results in increased droplet mobility and efficient droplet nucleation, which speeds up surface renewal; consequently, the condensation efficiency at high subcooling improves. Based on the current design, the VACNT-embedded hierarchical surfaces can be extended by different etching processes on various substrate materials to fulfill requirements for not only advanced condensation but also considerable use in anti-fogging, anti-icing, and highly efficient heat transfer applications.

#### 4. Materials and Methods

##### 4.1. Sample Preparation

Commercial 4-inch Si wafers (Suzhou Research Materials Microtech Co., Ltd., China) with a thickness of 0.5 mm were first patterned with microstructures by using photolithography; then, patterned cavities with a depth of 25 *μ*m were fabricated using dry etching. Next, the wafers were deposited with 5 nm Fe as catalyst by using an e-beam evaporator, after which a lift-off process was carried out to remove the photoresist. Finally, the microfabricated Si wafer was placed in a low-pressure CVD system to grow VACNTs by using acetylene as precursor. The reaction took place at a temperature of 973 K. Multiwalled CNTs with a diameter of were grown for 4 min. The surface samples were cut into pieces of for subsequent testing. A field-emission scanning electron microscope (SEM, FEI, Sirion 200) was used to characterize the surface morphology and the lateral profile of the VACNTs in the cavities.

##### 4.2. Condensation Experiments

Condensation tests were conducted under standard atmospheric pressure. Two types of steam conditions are prepared: ambient atmospheric conditions (~297 K) and hot steam at ~327 K. The hierarchical surfaces of different sizes were vertically placed on a cold source platform to attain temperatures of ~299 K and~289 K in different cases. Results from tests at three values of were analyzed: the cooling surface under ambient temperature conditions (= 8 K), the room-temperature surface in water-vapor steam (= 28 K), and the cooling surface in water-vapor steam (= 38 K). The condensation processes were visualized and recorded under an optical microscope. Statistics related to droplet distribution and size were obtained by frame analysis methods using the Nano Measurer 1.2 software package.

##### 4.3. Contact Angle (CA) Test

The apparent CAs were measured using a CA meter (KRUSS, DAS 100). The CA values were tested using a droplet volume of 10 *μ*L, and the average of the values at three positions on each sample was computed.

##### 4.4. Heat Conduction Measurement

The heat transfer diffusivity coefficients of the surfaces were measured at room temperature by using a laser thermal conductivity meter (NETZSCH, LFA 467 HyperFlash). The surface heating process was recorded, and the thermal diffusivity of the sample at room temperature was given by where represents the thickness of the surface and is the time required to reach half of the peak value of temperature.

##### 4.5. Mechanical Tests

Scratching tests were conducted using a steel scraper that provided tangential force, and compression tests were carried out using a 500 g weight loaded vertically on the surfaces. After each robustness test, the CA values were measured three times, and the average value was calculated. Each testing cycle was conducted 15 times.

##### 4.6. Molecular Dynamics (MD) Simulations

In the condensation simulations, a uniform VACNT surface (modeling the CNT area in cavities) and a heterogeneous surface (modeling the hydrophilic edges near hydrophobic zones) were studied and compared. A periodic boundary was applied to the simulation box with the size of . In the two surface models, the CNTs arranged in an orderly manner had a height of 15 Å, and the heterogeneous surface model had a hydrophilic zone in the middle. The hydrophilic Si region was modeled as a diamond lattice structure with a lattice constant of 5.43 Å.

For the solid surface, the carbon-carbon atom interaction was described using the AIREBO potential [53], whereas the Si-Si and Si-carbon interactions were described using Tersoff potential energy [54]. For the water droplets, the SPC/E rigid water model [55] was used; in this model, the water molecules had an internal bond angle of 109.47° and a bond length of 1 Å. Coulomb long-range interactions were applied to both hydrogen and oxygen atoms for producing, and a Lennard-Jones (LJ) potential was applied only to oxygen atoms for generating van der Waals interactions. The LJ potential was also used to describe the interaction between water molecules and solid surfaces. The LJ potential parameters were set as follows: and for Si–O interaction, and and for C–O interaction. The LJ potential and electrostatic terms were both cut off at 15 Å using the long-range solver of the particle-particle particle-mesh solver (PPPM) summation method. All simulations were carried out using the LAMMPS software [56].

A time step of 2 fs was used for all simulation processes. The solid surface atoms at the bottom of 3 were fixed, with all forces on atoms set to zero during the entire simulation. At the start, the vapor region consisting of 2093 water molecules was heated to 500 K in a canonical ensemble (NVT), and the unfixed surface atoms were relaxed at 300 K for 60 ps. Then, the vapor molecules were gradually cooled by heat transfer from the cold source surface in a microcanonical ensemble (NVE), while the surface remained under NVT thermodynamic control at 300 K. The duration of the process was 1 ns, and nuclei were formed on the surface. A cluster of H_{2}O nuclei is defined as a set of H_{2}O molecules, each of which is located at a distance of less than 3.36 Å from at least one other molecule in the cluster.

#### Data Availability

All relevant data that support the findings are available within this article and supporting information and are also available from authors upon reasonable request.

#### Conflicts of Interest

The authors declare no competing financial interest.

#### Acknowledgments

L.-W. Zhang acknowledges the support from the Marine Equipment Foresight Innovation Union Project (A-03).

#### Supplementary Materials

*Supplementary 1. *Section S1. Patterned surfaces with different shapes. Section S2. Transition from filmwise to dropwise condensation after planting of VACNTs. Section S3. Thermal diffusivity of the hierarchical surfaces. Section S4. Contact area fraction measurement on the VACNTs. Section S5. Condensation model for heat transfer rate. Section S6. Initial stage of condensation on the hierarchical surfaces. Section S7. Size distribution of condensation droplets on different surfaces. Section S8. Size-dependent analyses of condensation efficiency. Section S9. Snapshots of condensation process on the heterogeneous surface and on the as-grown VACNT surface.

*Supplementary 2. *Video S1. Condensation on the hierarchical surface at a subcooling of 8 K.

*Supplementary 3. *Video S2. Condensation on the hierarchical surface at a subcooling of 28 K.

*Supplementary 4. *Video S3. Condensation on the hierarchical surface at a subcooling of 38 K.

#### References

- K. Park, P. Kim, A. Grinthal et al., “Condensation on slippery asymmetric bumps,”
*Nature*, vol. 531, no. 7592, pp. 78–82, 2016. View at: Publisher Site | Google Scholar - X. J. Feng and L. Jiang, “Design and creation of superwetting/antiwetting surfaces,”
*Advanced Materials*, vol. 18, no. 23, pp. 3063–3078, 2006. View at: Publisher Site | Google Scholar - R. F. Wen, X. H. Ma, Y. C. Lee, and R. G. Yang, “Liquid-vapor phase-change heat transfer on functionalized nanowired surfaces and beyond,”
*Joule*, vol. 2, no. 11, pp. 2307–2347, 2018. View at: Publisher Site | Google Scholar - C. W. Lo, Y. C. Chu, M. H. Yen, and M. C. Lu, “Enhancing condensation heat transfer on three-dimensional hybrid surfaces,”
*Joule*, vol. 3, no. 11, pp. 2806–2823, 2019. View at: Publisher Site | Google Scholar - S. Gao, W. Liu, and Z. C. Liu, “Tuning nanostructured surfaces with hybrid wettability areas to enhance condensation,”
*Nanoscale*, vol. 11, no. 2, pp. 459–466, 2019. View at: Publisher Site | Google Scholar - S. L. Feng, J. Delannoy, A. Malod, H. Zheng, D. Quéré, and Z. Wang, “Tip-induced flipping of droplets on Janus pillars: from local reconfiguration to global transport,”
*Science Advances*, vol. 6, no. 28, article eabb4540, 2020. View at: Publisher Site | Google Scholar - S. Feng, P. Zhu, H. Zheng et al., “Three-dimensional capillary ratchet-induced liquid directional steering,”
*Science*, vol. 373, no. 6561, pp. 1344–1348, 2021. View at: Publisher Site | Google Scholar - X. Hu, Q. Yi, X. Kong, and J. Wang, “A review of research on dropwise condensation heat transfer,”
*Applied Sciences*, vol. 11, no. 4, p. 1553, 2021. View at: Publisher Site | Google Scholar - M. C. Lu, C. C. Lin, C. W. Lo, C. W. Huang, and C. C. Wang, “Superhydrophobic Si nanowires for enhanced condensation heat transfer,”
*International Journal of Heat and Mass Transfer*, vol. 111, pp. 614–623, 2017. View at: Publisher Site | Google Scholar - Q. Y. Zhang, L. L. He, X. F. Zhang, D. L. Tian, and L. Jiang, “Switchable direction of liquid Transportviaan anisotropic microarray surface and thermal stimuli,”
*ACS Nano*, vol. 14, no. 2, pp. 1436–1444, 2020. View at: Publisher Site | Google Scholar - Y. Lin, J. Han, M. Cai et al., “Durable and robust transparent superhydrophobic glass surfaces fabricated by a femtosecond laser with exceptional water repellency and thermostability,”
*Journal of Materials Chemistry A*, vol. 6, no. 19, pp. 9049–9056, 2018. View at: Publisher Site | Google Scholar - A. Starostin, V. Strelnikov, V. Valtsifer, I. Lebedeva, I. Legchenkova, and E. Bormashenko, “Robust icephobic coating based on the spiky fluorinated Al
_{2}O_{3}particles,”*Scientific Reports*, vol. 11, no. 1, article 5394, 2021. View at: Publisher Site | Google Scholar - R. Wen, S. Xu, D. Zhao et al., “Sustaining enhanced condensation on hierarchical mesh-covered surfaces,”
*National Science Review*, vol. 5, no. 6, pp. 878–887, 2018. View at: Publisher Site | Google Scholar - R. F. Wen, S. S. Xu, D. L. Zhao, Y. C. Lee, X. Ma, and R. Yang, “Hierarchical superhydrophobic surfaces with micropatterned nanowire arrays for high-efficiency jumping droplet condensation,”
*ACS Applied Materials & Interfaces*, vol. 9, no. 51, pp. 44911–44921, 2017. View at: Publisher Site | Google Scholar - Y. Hou, M. Yu, X. Chen, Z. Wang, and S. Yao, “Recurrent filmwise and dropwise condensation on a beetle mimetic surface,”
*ACS Nano*, vol. 9, no. 1, pp. 71–81, 2015. View at: Publisher Site | Google Scholar - X. Yan, F. Chen, S. Sett et al., “Hierarchical condensation,”
*ACS Nano*, vol. 13, no. 7, pp. 8169–8184, 2019. View at: Publisher Site | Google Scholar - B. El Fil, G. Kini, and S. Garimella, “A review of dropwise condensation: theory, modeling, experiments, and applications,”
*International Journal of Heat and Mass Transfer*, vol. 160, article 120172, 2020. View at: Publisher Site | Google Scholar - A. Goswami, S. C. Pillai, and G. McGranaghan, “Surface modifications to enhance dropwise condensation,”
*Surfaces and Interfaces*, vol. 25, article 101143, 2021. View at: Publisher Site | Google Scholar - H. J. Cho, D. J. Preston, Y. Zhu, and E. N. Wang, “Nanoengineered materials for liquid–vapour phase-change heat transfer,”
*Nature Reviews Materials*, vol. 2, no. 2, article 16092, 2016. View at: Google Scholar - X. M. Dai, N. Sun, S. O. Nielsen et al., “Hydrophilic directional slippery rough surfaces for water harvesting,”
*Science Advances*, vol. 4, no. 3, 2018. View at: Publisher Site | Google Scholar - T. Zhang, L. Mou, J. Zhang, L. Fan, and J. Li, “A visualized study of enhanced steam condensation heat transfer on a honeycomb-like microporous superhydrophobic surface in the presence of a non- condensable gas,”
*International Journal of Heat and Mass Transfer*, vol. 150, article 119352, 2020. View at: Publisher Site | Google Scholar - R. Wang, F. Wu, F. Yu, J. Zhu, X. Gao, and L. Jiang, “Anti-vapor-penetration and condensate microdrop self-transport of superhydrophobic oblique nanowire surface under high subcooling,”
*Nano Research*, vol. 14, no. 5, pp. 1429–1434, 2021. View at: Publisher Site | Google Scholar - Y. Tang, X. Yang, Y. Li, Y. Lu, and D. Zhu, “Robust micro-nanostructured superhydrophobic surfaces for long-term dropwise condensation,”
*Nano Letters*, vol. 21, no. 22, pp. 9824–9833, 2021. View at: Publisher Site | Google Scholar - J. Yong, Q. Yang, X. Hou, and F. Chen, “Nature-inspired superwettability achieved by femtosecond lasers,”
*Research*, vol. 2022, article 9895418, 51 pages, 2022. View at: Publisher Site | Google Scholar - J. Y. Long, P. Y. Zhou, Y. J. Huang, and X. Xie, “Enhancing the long-term robustness of dropwise condensation on nanostructured superhydrophobic surfaces by introducing 3D conical microtextures prepared by femtosecond laser,”
*Advanced Materials Interfaces*, vol. 7, no. 21, article 2000997, 2020. View at: Publisher Site | Google Scholar - A. Tripathy, C. W. E. Lam, D. Davila et al., “Ultrathin lubricant-infused vertical graphene nanoscaffolds for high-performance dropwise condensation,”
*ACS Nano*, vol. 15, no. 9, pp. 14305–14315, 2021. View at: Publisher Site | Google Scholar - L. Wang, J. Li, B. Zhang et al., “Counterintuitive ballistic and directional liquid transport on a flexible droplet rectifier,”
*Research*, vol. 2020, article 6472313, 11 pages, 2020. View at: Publisher Site | Google Scholar - A. I. Aria and M. Gharib, “Physicochemical characteristics and droplet impact dynamics of superhydrophobic carbon nanotube arrays,”
*Langmuir*, vol. 30, no. 23, pp. 6780–6790, 2014. View at: Publisher Site | Google Scholar - J. Wang, T. Li, B. Xia et al., “Vapor-condensation-assisted optical microscopy for ultralong carbon nanotubes and other nanostructures,”
*Nano Letters*, vol. 14, no. 6, pp. 3527–3533, 2014. View at: Publisher Site | Google Scholar - R. A. Pinheiro, F. M. Rosa, R. M. Volú, G. de Vasconcelos, V. J. Trava-Airoldi, and E. J. Corat, “Vertically aligned carbon nanotubes (VACNT) surfaces coated with polyethylene for enhanced dew harvesting,”
*Diamond and Related Materials*, vol. 107, article 107837, 2020. View at: Publisher Site | Google Scholar - R. F. Wen, S. S. Xu, X. H. Ma, Y. C. Lee, and R. G. Yang, “Three-dimensional superhydrophobic nanowire networks for enhancing condensation heat transfer,”
*Joule*, vol. 2, no. 2, pp. 269–279, 2018. View at: Publisher Site | Google Scholar - V. Kondrashov and J. Ruhe, “Microcones and nanograss: toward mechanically robust superhydrophobic surfaces,”
*Langmuir*, vol. 30, no. 15, pp. 4342–4350, 2014. View at: Publisher Site | Google Scholar - J. J. Schneider, “Vertically aligned carbon nanotubes as platform for biomimetically inspired mechanical sensing, bioactive surfaces, and electrical cell interfacing,”
*Advanced Biosystems*, vol. 1, no. 11, article 1700101, 2017. View at: Publisher Site | Google Scholar - D. Wang, Q. Sun, M. J. Hokkanen et al., “Design of robust superhydrophobic surfaces,”
*Nature*, vol. 582, no. 7810, pp. 55–59, 2020. View at: Publisher Site | Google Scholar - Y. Cheng, M. Wang, J. Sun et al., “Rapid and persistent suction condensation on hydrophilic surfaces for high-efficiency water collection,”
*Nano Letters*, vol. 21, no. 17, pp. 7411–7418, 2021. View at: Publisher Site | Google Scholar - H. Cha, H. Vahabi, A. Wu et al., “Dropwise condensation on solid hydrophilic surfaces,”
*Science Advances*, vol. 6, no. 2, article eaax0746, 2020. View at: Publisher Site | Google Scholar - K. Liu, K. Jiang, C. Feng, Z. Chen, and S. Fan, “A growth mark method for studying growth mechanism of carbon nanotube arrays,”
*Carbon*, vol. 43, no. 14, pp. 2850–2856, 2005. View at: Publisher Site | Google Scholar - R. Enright, N. Miljkovic, A. Al-Obeidi, C. V. Thompson, and E. N. Wang, “Condensation on superhydrophobic surfaces: the role of local energy barriers and structure length scale,”
*Langmuir*, vol. 28, no. 40, pp. 14424–14432, 2012. View at: Publisher Site | Google Scholar - K. X. Meng, W. L. Fan, and H. Wang, “Dynamic scenario simulation of dropwise condensation on a superhydrophobic surface with droplet jumping,”
*Applied Thermal Engineering*, vol. 148, pp. 316–323, 2019. View at: Publisher Site | Google Scholar - C. W. Lo, C. C. Wang, and M. C. Lu, “Scale effect on dropwise condensation on superhydrophobic surfaces,”
*ACS Applied Materials & Interfaces*, vol. 6, no. 16, pp. 14353–14359, 2014. View at: Publisher Site | Google Scholar - N. Vargaftik, B. Volkov, and L. Voljak, “International tables of the surface tension of water,”
*Journal of Physical and Chemical Reference Data*, vol. 12, no. 3, pp. 817–820, 1983. View at: Publisher Site | Google Scholar - Y. Suh, J. Lee, P. Simadiris et al., “A deep learning perspective on dropwise condensation,”
*Advanced Science*, vol. 8, no. 22, article 2101794, 2021. View at: Publisher Site | Google Scholar - D. W. Jeong, U. H. Shin, J. H. Kim, S. H. Kim, H. W. Lee, and J. M. Kim, “Stable hierarchical superhydrophobic surfaces based on vertically aligned carbon nanotube forests modified with conformal silicone coating,”
*Carbon*, vol. 79, pp. 442–449, 2014. View at: Publisher Site | Google Scholar - Z. Y. Zhou, T. C. Gao, S. McCarthy et al., “Parahydrophobicity and stick-slip wetting dynamics of vertically aligned carbon nanotube forests,”
*Carbon*, vol. 152, pp. 474–481, 2019. View at: Publisher Site | Google Scholar - A. Starostin, V. Valtsifer, Z. Barkay, I. Legchenkova, V. Danchuk, and E. Bormashenko, “Drop-wise and film-wise water condensation processes occurring on metallic micro-scaled surfaces,”
*Applied Surface Science*, vol. 444, pp. 604–609, 2018. View at: Publisher Site | Google Scholar - S. Karthika, T. K. Radhakrishnan, and P. Kalaichelvi, “A review of classical and nonclassical nucleation theories,”
*Crystal Growth & Design*, vol. 16, no. 11, pp. 6663–6681, 2016. View at: Publisher Site | Google Scholar - D. Yan, Q. Zeng, S. Xu, Q. Zhang, and J. Wang, “Heterogeneous nucleation on concave rough surfaces: thermodynamic analysis and implications for nucleation design,”
*Journal of Physical Chemistry C*, vol. 120, no. 19, pp. 10368–10380, 2016. View at: Publisher Site | Google Scholar - L. He, A. Karumuri, and S. M. Mukhopadhyay, “Wettability tailoring of nanotube carpets: morphology-chemistry synergy for hydrophobic-hydrophilic cycling,”
*RSC Advances*, vol. 7, no. 41, pp. 25265–25275, 2017. View at: Publisher Site | Google Scholar - D. Niu and G. Tang, “Molecular dynamics simulation of droplet nucleation and growth on a rough surface: revealing the microscopic mechanism of the flooding mode,”
*RSC Advances*, vol. 8, no. 43, pp. 24517–24524, 2018. View at: Publisher Site | Google Scholar - Q. Peng, L. Jia, J. Guo et al., “Forced jumping and coalescence-induced sweeping enhanced the dropwise condensation on hierarchically microgrooved superhydrophobic surface,”
*Applied Physics Letters*, vol. 114, no. 13, article 133106, 2019. View at: Publisher Site | Google Scholar - Z. Yin, H. Wang, M. Jian et al., “Extremely black vertically aligned carbon nanotube arrays for solar steam generation,”
*ACS Applied Materials & Interfaces*, vol. 9, no. 34, pp. 28596–28603, 2017. View at: Publisher Site | Google Scholar - J. Song, W. F. Cheng, M. T. Nie et al., “Partial leidenfrost evaporation-assisted ultrasensitive surface-enhanced Raman spectroscopy in a Janus water droplet on hierarchical plasmonic micro-/nanostructures,”
*ACS Nano*, vol. 14, no. 8, pp. 9521–9531, 2020. View at: Publisher Site | Google Scholar - S. J. Stuart, A. B. Tutein, and J. A. Harrison, “A reactive potential for hydrocarbons with intermolecular interactions,”
*Journal of Chemical Physics*, vol. 112, no. 14, pp. 6472–6486, 2000. View at: Publisher Site | Google Scholar - S. Munetoh, T. Motooka, K. Moriguchi, and A. Shintani, “Interatomic potential for Si-O systems using Tersoff parameterization,”
*Computational Materials Science*, vol. 39, no. 2, pp. 334–339, 2007. View at: Publisher Site | Google Scholar - H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, “The missing term in effective pair potentials,”
*Journal of Physical Chemistry*, vol. 91, no. 24, pp. 6269–6271, 1987. View at: Publisher Site | Google Scholar - S. Plimpton, “Fast parallel algorithms for short-range molecular dynamics,”
*Journal of Computational Physics*, vol. 117, no. 1, pp. 1–19, 1995. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2022 Yue Hu et al. Exclusive Licensee Science and Technology Review Publishing House. Distributed under a Creative Commons Attribution License (CC BY 4.0).